Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variability analysis”

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 1,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

DER Digital Supply Chain Gap Analysis

Solar photovoltaic (PV) cybersecurity is a growing field of research. As deployments of solar PV have increased, cyber risk has also increased. Utility solar PV installations, however, are not required to comply with the North American Electric Reliability Corporation (NERC) Critical Infrastructure Protection (CIP) plan unless they meet a minimum generation threshold of 75 MW. Individual residential-scale solar PV deployments will not meet that generation threshold and are therefore excluded from the NERC CIP requirements. With most solar installations less than 75 MW, solar PV has been deployed with minimal oversight and highly variable cybersecurity maturity. The resources that comprise the digital supply chain can include software, code, data, and other digital components. But as clean energy technologies advance, cybersecurity threats and vulnerabilities continue to evolve and grow in sophistication. Supply chain cybersecurity represents a critical area for ensuring safe operations as the U.S. moves toward a clean energy future.

cybersecurity↗
DOE OSTI2025

Phase-field modeling and experiments of dynamic fracture in single crystal quartz

Predicting the onset and characteristics of brittle fracture is important for a wide range of engineering and geological material applications. In this paper, we study important aspects of brittle fracture in α-quartz by phase-field modeling and experiments using a top-down approach. In the modeling framework, the work term in the Griffith energy balance is replaced with internal energy contributions that represent surface energy, thermal energy, and elastic strain energy stored in defects. This allows parametrization of individual energy contributions in terms of internal state variables and keeps track of energy partitioning after the onset of fracture. The path and history dependence of fracture is included in evolution laws for internal state variables, e.g., entropy evolution, while the energy remains a true potential. In the experimental part, dynamic compression experiments coupled with X-ray phase contrast imaging are performed on cube-like samples with a hole. In the top-down analysis, dynamic compression and three point bending experiments from the literature are simulated with the developed phase-field damage model. In conclusion, the fitted model highlights the strain rate, size, and stress state dependence of damage nucleation and evolution in single crystal α-quartz.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Determination of a suitable molar absorption coefficient (ε) for lignin analysis of fibrous plants using the CASA method

In this study, the benchmarking and applicability of the CASA method for the analysis of fibrous plants, such as flax, hemp, and jute, is evaluated. Lignin is a phenolic biopolymer present in plant cell walls and is composed of three primary monomeric units, designated G, S, and H. Various factors, such as the genetic variability, influence the relative proportions of these units in plant samples. Recently, the cysteine-assisted sulfuric acid (CASA) method has been introduced as a rapid method for the quantification of lignin in wood samples. The aim of this study is to establish a suitable molar absorption coefficient (ε, L·g –1 ·cm –1 ) to adapt the CASA method for use with annual plant fibers. This investigation was motivated by the technical advantages of the CASA method, including higher throughput, lower reaction temperatures, and ecological benefits due to the absence of carcinogenic, mutagenic, or reprotoxic (CMR) substances and the need for minimal sample quantities. In this method, lignin solubilization is facilitated by cysteine, which is an amino acid that enhances the reaction kinetics, using a one-hour incubation period. However, as with any spectrophotometric technique, the CASA method depends on a molar absorption coefficient (ε) that varies according to the ratio of the aromatic units within the polymer. To evaluate the suitability of CASA for quantifying lignin in economically significant plant fibers, we investigate the impact of the unit ratios on the accuracy of ε. The results are compared with those of two widely recognized lignin analytical methods: the Klason method, which is a gravimetric reference method, and the acetyl bromide soluble lignin method, which is the most commonly used spectrophotometric approach. The final ε obtained in this study reveals a relative difference in lignin content ranging from –8 % to +9 % based on a comparison between the CASA and Klason methods across different industrial hemp varieties. As a result, using known G:S ratios in annual fibrous plants, ε can be estimated from our results.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

A generalizable machine learning-assisted fast Fourier transform algorithm to simulate the large strain phenomena in polycrystalline materials

Machine learning methods have shown initial promise in constitutive modeling for single crystals or homogenized polycrystals, delivering notable computational efficiency. However, existing machine learning-based constitutive models often lack generalizability, limiting their application across diverse boundary value problems. This study introduces a thermodynamics-informed artificial neural network model to accelerate rate-tangent crystal plasticity fast Fourier transform simulations for cross-scale deformation behaviors of polycrystals under complex loading. Our model integrates microstructural variability and local interactions effectively. To address local effects in each grain, we employ K-means clustering to group Gauss points within the microstructure into clusters assumed to be in similar mechanical states. This approach, based on self-clustering analysis, extends model scope from macroscopic stress response to the granular level, capturing mechanical responses and orientation evolution across grains. This reduces the number of nonlinear problems to solve, with cluster responses propagated throughout each group. The thermodynamics-based artificial neural network-extracted features are further processed using local material state clusters to account for history-dependent deformation and evolving microstructures. Additionally, representative volume element simulations with rate-tangent crystal plasticity fast Fourier transform provide reliable datasets for model training. The proposed model demonstrates high efficiency, accuracy, self-consistency, and enhanced generalizability in predicting strain–stress responses and orientation evolution at both individual grain and aggregate scales under complex loading conditions, such as biaxial tension and arbitrary loading scenarios.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A Numerical Modeling Framework for Flocculation and Cohesive Sediment Transport in the Wave Bottom Boundary Layer

Flocculation, a critical process in coastal and estuarine systems, plays a significant role in sediment transport, nutrient cycling, and ecological health. This study develops a cohesive sediment transport modeling framework tailored to the wave bottom boundary layer under dilute and equilibrium conditions, explicitly incorporating flocculation effects via a Population Balance Equation (PBE). Using Direct Numerical Simulation, six baseline cases, each with a distinct sediment concentration profile resulting from a constant settling velocity and critical erosion shear stress, are generated to drive the PBE flocculation model for given floc yield strength and stickiness. Results reveal that flocculation significantly influences sediment concentration profiles promoting three distinct stages, well‐mixed, transition to lutocline, and well‐developed lutocline. At low concentrations with well‐mixed profiles, cohesive floc properties are less significant, and turbulence is a main flocculation driver. In contrast, as concentration increases, cohesive floc properties become crucial, facilitating lutocline formation. Here, the analysis also highlights limitations of depth‐averaged settling velocity as a parameterization. It is suitable for well‐mixed and transitional profiles but fails in well‐developed lutoclines, where empirical formulations that explicitly incorporate turbulent shear rate and sediment concentration better capture variability. This study underscores the necessity of incorporating flocculation effects into sediment transport models to enhance predictions of sediment dynamics in wave bottom boundary layers.

Penaloza‐Giraldo, Jorge A. [Oak Ridge National Lab↗
DOE OSTI2023

Assessing residential PM 2.5 concentrations and infiltration factors with high spatiotemporal resolution using crowdsourced sensors

Building conditions, outdoor climate, and human behavior influence residential concentrations of fine particulate matter (PM 2.5 ). To study PM 2.5 spatiotemporal variability in residences, we acquired paired indoor and outdoor PM 2.5 measurements at 3,977 residences across the United States totaling >10,000 monitor-years of time-resolved data (10-min resolution) from the PurpleAir network. Time-series analysis and statistical modeling apportioned residential PM 2.5 concentrations to outdoor sources (median residential contribution = 52% of total, coefficient of variation = 69%), episodic indoor emission events such as cooking (28%, CV = 210%) and persistent indoor sources (20%, CV = 112%). Residences in the temperate marine climate zone experienced higher infiltration factors, consistent with expectations for more time with open windows in milder climates. Likewise, for all climate zones, infiltration factors were highest in summer and lowest in winter, decreasing by approximately half in most climate zones. Large outdoor–indoor temperature differences were associated with lower infiltration factors, suggesting particle losses from active filtration occurred during heating and cooling. Absolute contributions from both outdoor and indoor sources increased during wildfire events. Infiltration factors decreased during periods of high outdoor PM 2.5 , such as during wildfires, reducing potential exposures from outdoor-origin particles but increasing potential exposures to indoor-origin particles. Time-of-day analysis reveals that episodic emission events are most frequent during mealtimes as well as on holidays (Thanksgiving and Christmas), indicating that cooking-related activities are a strong episodic emission source of indoor PM 2.5 in monitored residences.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Deconvolving the roles of E × B shear and pedestal structure in the energy confinement quality of super H-mode experiements

Analysis of “super H-mode” experiments on DIII-D has put forward that high plasma toroidal rotation, not high pedestal, plays the essential role in achieving energy confinement quality H 98y2 >>1. Recently, super H-mode experiments with variable input torque have confirmed that high rotation shear discharges have very high levels of H 98y2 (>1.5), independent of the pedestal height, and that high pedestal discharges with low rotation shear have levels of H 98y2 only slightly above 1 (≤1.2). Although some increase in stored energy with higher pedestal occurs, the energy confinement quality mainly depends on the toroidal rotation shear, which varies according to different levels of injected neutral beam torque per particle. Here, quasi-linear gyrofluid modeling achieves a good match of the experiment when including the E×B shear; without including plasma rotation, the modeling predicts a confinement quality consistent with the empirical observation of H 98y2 ~1.2 at low rotation. Nonlinear gyrokinetic transport modeling shows that the effect of E×B turbulence stabilization is far larger than other mechanisms, such as the so-called hot-ion stabilization (T i /T e ) effect. Consistent with these experimental and modeling results are previous simulations of the ITER Baseline Scenario using a super H-mode pedestal solution, which showed the potential to exceed the Q=10 target if the pedestal density could be increased above the Greenwald limit. A close look at these simulations reveals that the predicted energy confinement quality is below 1 even at the highest pedestal pressure. The improvement in Q at higher pedestal density is due to the improved fusion power generation at the higher core density associated with higher pedestal density, not to an improved energy confinement quality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Detailed space–time variations of the seismic response of the shallow crust to small earthquakes from analysis of dense array data

SUMMARY We compute high-resolution space–time variations of subsurface seismic properties from autocorrelation functions (ACF’s) of noise and local earthquakes, recorded by the Sage Brush Flat dense array deployed around the Clark branch of the San Jacinto fault. The resolved temporal changes are referred to as apparent velocity changes because they reflect both nonlinear response and variations of material properties such as cracking and damage. Apparent velocity changes are estimated at four frequency bands (10–15, 10–20, 15–30 and 20–40 Hz) for two local earthquake data sets. In one analysis, ACF’s from P- and S-wave windows of 31 small events with magnitudes below 3.1 are used to compute the apparent velocity variations with respect to the mean ACF of each phase, and we also use the mean ACF of noise data as reference to estimate the changes. In a further analysis, the temporal evolution of properties is computed using moving time windows in continuous waveform over one-hour long data with noise and earthquake signals. The apparent velocity changes and recovery times are frequency dependent and present a strong spatial variability across the array. The resolved changes are larger and recovery time shorter with data associated with higher frequencies. At frequencies larger than 15 Hz, and using the mean ACF of noise data as a reference, the apparent average velocity changes across the array during the passage of the P and S waves from the small local events are 2.5 per cent and 6 per cent, respectively. The apparent velocity changes decrease by one order of magnitude when the earthquake data are used as a reference. The relatively large changes in response to very low ground motion have important implications on nonlinear processes involving degradation and healing of the subsurface material during common earthquake shaking.

Geochemistry & Geophysics↗
DOE OSTI2024

Foundations of automatic feature extraction at LHC–point clouds and graphs

Abstract Deep learning algorithms will play a key role in the upcoming runs of the Large Hadron Collider (LHC), helping bolster various fronts ranging from fast and accurate detector simulations to physics analysis probing possible deviations from the Standard Model. The game-changing feature of these new algorithms is the ability to extract relevant information from high-dimensional input spaces, often regarded as “replacing the expert” in designing physics-intuitive variables. While this may seem true at first glance, it is far from reality. Existing research shows that physics-inspired feature extractors have many advantages beyond improving the qualitative understanding of the extracted features. In this review, we systematically explore automatic feature extraction from a phenomenological viewpoint and the motivation for physics-inspired architectures. We also discuss how prior knowledge from physics results in the naturalness of the point cloud representation and discuss graph-based applications to LHC phenomenology.

Bhardwaj, Akanksha↗
DOE OSTI2021

Polarimetric Properties of Event Horizon Telescope Targets from ALMA

We present the results from a full polarization study carried out with the Atacama Large Millimeter/submillimeter Array (ALMA) during the first Very Long Baseline Interferometry (VLBI) campaign, which was conducted in 2017 April in the λ3 mm and λ1.3 mm bands, in concert with the Global mm-VLBI Array (GMVA) and the Event Horizon Telescope (EHT), respectively. We determine the polarization and Faraday properties of all VLBI targets, including Sgr A*, M87, and a dozen radio-loud active galactic nuclei (AGNs), in the two bands at several epochs in a time window of 10 days. We detect high linear polarization fractions (2%–15%) and large rotation measures (RM > 10 3.3 –10 5.5 rad m -2 ), confirming the trends of previous AGN studies at millimeter wavelengths. We find that blazars are more strongly polarized than other AGNs in the sample, while exhibiting (on average) order-of-magnitude lower RM values, consistent with the AGN viewing angle unification scheme. For Sgr A* we report a mean RM of (-4.2 ± 0.3) × 10 5 rad m -2 at 1.3 mm, consistent with measurements over the past decade and, for the first time, an RM of (–2.1 ± 0.1) × 10 5 rad m -2 at 3 mm, suggesting that about half of the Faraday rotation at 1.3 mm may occur between the 3 mm photosphere and the 1.3 mm source. We also report the first unambiguous measurement of RM toward the M87 nucleus at millimeter wavelengths, which undergoes significant changes in magnitude and sign reversals on a one year timescale, spanning the range from -1.2 to 0.3 × 10 5 rad m -2 at 3 mm and -4.1 to 1.5 × 10 5 rad m -2 at 1.3 mm. Given this time variability, we argue that, unlike the case of Sgr A*, the RM in M87 does not provide an accurate estimate of the mass accretion rate onto the black hole. We put forward a two-component model, comprised of a variable compact region and a static extended region, that can simultaneously explain the polarimetric properties observed by both the EHT (on horizon scales) and ALMA (which observes the combined emission from both components). These measurements provide critical constraints for the calibration, analysis, and interpretation of simultaneously obtained VLBI data with the EHT and GMVA.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Updates to Relevance Vector Machine: Multiclass Classification, Variable Selection, and Proof-of-Concept Application to Safeguards Fresh Fuel Verification using List-Mode Neutron Collar Data

To expand the capabilities of safeguards authorities to verify the integrity of fresh fuel assemblies, Oak Ridge National Laboratory has retrofit the existing electronics of the JCC-71 uranium neutron coincidence collar, which contains 18 3 He neutron detectors and an external 241 AmLi(α, n) neutron interrogation source arranged to surround a fresh nuclear fuel assembly. The new electronics system allows analysts to record list-mode neutron multiplicity data in addition to the singles and doubles rates that are currently measured. Based on previous proof-of-concept research, analysis of these new data will identify off-normal fuel configurations in an assembly and characterize or localize the specific partial fuel defects. The purpose of this report it to document the analysis algorithm development and then to demonstrate its capability for the safeguards verification of fresh fuel assemblies using list mode neutron collar data. To analyze the complex list-mode data collected with the upgraded uranium neutron collar, multivariate classification algorithms are being developed using a novel classification method, the relevance vector machine. This approach may be applied to multiclass problems to estimate the probability that test data belongs to one of many possible classes of data. In addition, our method identifies the most useful variables/channels for making predictions, which illuminates the basis for the model’s predictions, and this interpretability is largely unique among data analytics methods. Variable selection occurs during model training and parameter tuning and does not need any external hyperparameter tuning routines. Finally, we apply the modified relevance vector machine to a simulated dataset of list-mode neutron collar data generated with the radiation transport code MCNP. The method can correctly identify off-normal fuel configurations, categorize the data according to four fuel defect scenarios, and rank the channels in the data according to prediction utility. For nuclear safeguards applications, it is concluded that this method has the potential to increase the sensitivity and reliability to detect missing fuel rods from a standard 17 x 17 Pressurized Water Reactor (PWR) fresh fuel assembly. Within this analysis, “off-normal” (i.e., missing fuel rods) were correctly classified in 17 simulated test scenarios with one quarter (25%) of the fresh fuel rods missing using a training data set of 58 simulated measurements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Impact of Geothermal District Heating System on Flexibility of Microgrid in Tuttle, Oklahoma

Flexibility is the capability of the power grid to maintain a balance between electricity generation and variable demand. This study presents preliminary results evaluating the impact of geothermal district heating systems on the flexibility of microgrid in Tuttle, Oklahoma. Heating demand profiles were modeled using EnergyPlus for the district that includes two schools and 250 single-family houses. Then, geothermal energy production was modeled using GEOPHIRES to estimate how much heating demand in the district can be supplied by five different geothermal system scenarios. The results indicated that geothermal energy production varied depending on the resource temperature at different depths, system configurations, and flow rates. For the grid flexibility analysis, electricity consumptions in the five geothermal systems were estimated for pump operations to circulate water from the wells to radiators, while electricity consumption by air-source heat pump in the base case was estimated to supply the same heating load. Electricity consumption in the geothermal systems was significantly lower than those in base cases. The electricity saved by the geothermal system was then incorporated into the microgrid electrical load profiles where variable renewable electricity generation is significantly high. The results visually showed that geothermal district heating system can improve grid flexibility as a baseload during the winter season. The results also highlighted potential opportunities to save energy costs that will be further analyzed in future study.

Cambium↗
DOE OSTI2024

Large Ensemble Exploration of Global Energy Transitions Under National Emissions Pledges

Global climate goals require a transition to a deeply decarbonized energy system. Meeting the objectives of the Paris Agreement through countries' nationally determined contributions and long-term strategies represents a complex problem with consequences across multiple systems shrouded by deep uncertainty. Robust, large-ensemble methods and analyses mapping a wide range of possible future states of the world are needed to help policymakers design effective strategies to meet emissions reduction goals. This study contributes a scenario discovery analysis applied to a large ensemble of 5,760 model realizations generated using the Global Change Analysis Model. Eleven energy-related uncertainties are systematically varied, representing national mitigation pledges, institutional factors, and techno-economic parameters, among others. The resulting ensemble maps how uncertainties impact common energy system metrics used to characterize national and global pathways toward deep decarbonization. Results show globally consistent but regionally variable energy transitions as measured by multiple metrics, including electricity costs and stranded assets. Larger economies and developing regions experience more severe economic outcomes across a broad sampling of uncertainty. The scale of CO 2 removal globally determines how much the energy system can continue to emit, but the relative role of different CO 2 removal options in meeting decarbonization goals varies across regions. Previous studies characterizing uncertainty have typically focused on a few scenarios, and other large-ensemble work has not (to our knowledge) combined this framework with national emissions pledges or institutional factors. Our results underscore the value of large-ensemble scenario discovery for decision support as countries begin to design strategies to meet their goals.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2026

Co-Simulation Meets AI: MCP-Driven Power System Analysis

GridGPT, a fine-tuned Generative AI model is designed for on-premise use in grid control rooms. This presentation will demonstrate how eGridGPT can seamlessly integrate with control room solutions to offer operators, engineers, and corporate users enhanced guidance and decision support. It is to show how this innovative AI solution can improve state estimation, boost variable energy forecasting, and optimize grid operations. By leveraging eGridGPT's unique features, audience will learn to unlock new levels of automation, predictive analytics, and reliability within their power systems, ultimately leading to reduced downtime and improved operational efficiency.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Visible-light Phase Curves from the Second Year of the TESS Primary Mission

We carried out a systematic study of full-orbit phase curves for known transiting systems in the northern ecliptic sky that were observed during Year 2 of the TESS primary mission. We applied the same methodology for target selection, data processing, and light-curve fitting as we did in our Year 1 study. Out of the 15 transiting systems selected for analysis, seven—HAT-P-7, KELT-1, KELT-9, KELT-16, KELT-20, Kepler-13A, and WASP-12—show statistically significant secondary eclipses and day–night atmospheric brightness modulations. Small eastward dayside hot-spot offsets were measured for KELT-9b and WASP-12b. KELT-1, Kepler-13A, and WASP-12 show additional phase-curve variability attributed to the tidal distortion of the host star; the amplitudes of these signals are consistent with theoretical predictions. We combined occultation measurements from TESS and Spitzer to compute dayside brightness temperatures, TESS-band geometric albedos, Bond albedos, and phase integrals for several systems. The new albedo values solidify the previously reported trend between dayside temperature and geometric albedo for planets with 1500 K < T {sub day} < 3000 K. For Kepler-13Ab, we carried out an atmospheric retrieval of the full secondary eclipse spectrum, which revealed a noninverted temperature–pressure profile, significant H{sub 2}O and K absorption in the near-infrared, evidence for strong optical atmospheric opacity due to sodium, and a confirmation of the high geometric albedo inferred from our simpler analysis. We explore the implications of the phase integrals (ratios of Bond to geometric albedos) for understanding exoplanet clouds. We also report updated transit ephemerides for all of the systems studied in this work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Impact of Geothermal District Heating System on Flexibility of Microgrid in Tuttle, Oklahoma: Preprint

Flexibility is the capability of the power grid to maintain a balance between electricity generation and variable demand. This study presents preliminary results evaluating the impact of geothermal district heating systems on the flexibility of a conceptual microgrid in Tuttle, Oklahoma. Heating demand profiles were modeled using EnergyPlus for the district that includes two schools and 250 single-family houses. Then, geothermal energy production was modeled using GEOPHIRES to estimate how much heating demand in the district can be supplied by five different geothermal system scenarios. The results indicated that geothermal energy production varied depending on the resource temperature at different depths, system configurations, and flow rates. For the grid flexibility analysis, electricity consumptions in the five geothermal systems were estimated for pump operations to circulate water from the wells to radiators, while electricity consumption by air-source heat pump in the base case was estimated to supply the same heating load. Electricity consumption in the geothermal systems was significantly lower than those in base cases. The electricity saved by the geothermal system was then incorporated into the microgrid electrical load profiles where variable renewable electricity generation is significantly high. The results visually showed that geothermal district heating system can improve grid flexibility as a baseload during the winter season. The results also highlighted potential opportunities to save energy costs that will be further analyzed in future study.

Cambium↗
DOE OSTI2022

Stand and environmental conditions drive functional shifts associated with mesophication in eastern US forests

There is a growing body of evidence that mesic tree species are increasing in importance across much of the eastern US. This increase is often observed in tandem with a decrease in the abundance and importance of species considered to be better adapted to disturbance and drier conditions (e.g., Quercus species). Concern over this transition is related to several factors, including the potential that this transition is self-reinforcing (termed “mesophication”), will result in decreased resiliency of forests to a variety of disturbances, and may negatively impact ecosystem functioning, timber value, and wildlife habitat. Evidence for shifts in composition provide broad-scale support for mesophication, but we lack information on the fine-scale factors that drive the associated functional changes. Understanding this variability is particularly important as managers work to develop site-and condition-specific management practices to target stands or portions of the landscape where this transition is occurring or is likely to occur in the future. To address this knowledge gap and identify forests that are most susceptible to mesophication (which we evaluate as a functional shift to less drought or fire tolerant, or more shade tolerant, forests), we used data from the USDA Forest Service Forest Inventory and Analysis program to determine what fine-scale factors impact the rate (change through time) and degree (difference between the overstory and midstory) of change in eastern US forests. We found that mesophication varies along stand and environmental gradients, but this relationship depended on the functional trait examined. For example, shade and drought tolerance suggest mesophication is greatest at sites with more acidic soils, while fire tolerance suggests mesophication increases with soil pH. Mesophication was also generally more pronounced in older stands, stands with more variable diameters, and in wetter sites, but plots categorized as “hydric” were often highly variable. Our results provide evidence that stand-scale conditions impact current and potential future changes in trait conditions and composition across eastern US forests. We provide a starting point for managers looking to prioritize portions of the landscape most at risk and developing treatments to address the compositional and functional changes associated with mesophication.

Woodbridge, Margaret↗
DOE OSTI2020

Comparison of four learning-based methods for predicting groundwater redox status

Knowing the location where groundwater denitrification occurs, or by proxy the groundwater redox status (oxic, mixed, and anoxic), is valuable information for assessing and managing potential agricultural land-use impacts on freshwater quality. We compare the efficacy of supervised (Linear Discriminant Analysis LDA; Boosted Regression Trees, BRT; and Random Forest, RF) and unsupervised (Modified Self-Organizing Map, MSOM) learning-based methods to predict groundwater redox status in the agriculturally dominated Tasman, Waikato, and Wellington regions of New Zealand. Thresholds applied to regional groundwater-quality samples provide redox status variables and learn heuristics constrained by these variables and applied to spatial factors (climate, elevation, geologic, hydrology soils, and well depth) identify optimal sets of regional predictor variables. A split-sample approach is used to train and test the learning methods ability to predict redox status using the optimal predictor variables. Overall, the supervised methods demonstrate a prediction bias toward oxic conditions and inability to perform statistically well when using independent regional data; for example, consider kappa statistics for BRT (Tasman: 0.42, Waikato: 0.38, Wellington: 0.17), RF (Tasman: 0.42, Waikato: 0.47, Wellington: 0.17 and LDA (Tasman: 0.46, Waikato: 0.32, Wellington: 0.17). By contrast, the unsupervised method performs statistically well when predicting oxic, mixed, and anoxic conditions and corresponding depths when using independent regional data; for example, consider MSOM kappa statistics for Tasman: 0.78, Waikato: 0.80, Wellington: 0.76. The unsupervised learning method provides the added benefits of being (1) able to combine predictions into 3D regional anoxic probability plots for interpreting the spatial influence of paleosols and groundwater flowpaths on redox status, and (2) readily extended to map 3D redox status across New Zealand and other countries despite data bias and sparsity.

machine learning algorithms, Random Forest, booste↗