Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variance enhancement”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Long-range Doppler lidar measurements of wind turbine wakes and their interaction with turbulent atmospheric boundary-layer flow at Perdigao 2017

As part of the Perdigão 2017 campaign, vertical RHI (range-height indicator) scans with long-range pulsed Doppler wind lidars were performed aligned with the main wind direction and a wind turbine (WT) located on a mountain ridge. The measurements are used to not only retrieve flow velocities, but also their variance and - by using the turbulent broadening of the Doppler spectrum - also turbulent kinetic energy (TKE) dissipation rate. The study shows that turbulence in the WT wake is dependent on the turbulence of the inflow, but also on atmospheric stability. In stable atmospheric conditions, wakes could be analyzed up to five rotor diameters downstream (D) and showed the maximum turbulence in the wake at 2-3 D, whereas in unstable conditions, the maximum was found at 2 D and the wake could not be detected further than 3 D. A clear dependency of wake turbulence enhancement on inflow turbulence intensity is found, which levels out to no further enhancement at turbulence intensities of 30%.

17 WIND ENERGY↗
DOE OSTI2024

Implications of stop-and-go traffic on training learning-based car-following control

Learning-based car-following control (LCC) of connected and autonomous vehicles (CAVs) is gaining significant attention with the advancement of computing power and data accessibility. While the flexibility and large model capacity of model-free architecture enable LCC to potentially outperform the model-based car-following (CF) model in improving traffic efficiency and mitigating congestion, the generalizability of LCC for traffic conditions different from the training environment/dataset is not well-understood. Herein, this study seeks to explore the impact of stop-and-go traffic in the training dataset on the generalizability of LCC. It uses the characteristics of lead vehicle trajectories to describe stop-and-go traffic, and links the theory of identifiability (i.e., obtaining a unique parameter estimation result using sensor measurements) to the generalizability of behavior cloning (BC) and policy-based deep reinforcement learning (DRL). Correspondingly, the study shows theoretically that: (i) stop-and-go traffic can enable the property of identifiability and enhance the control performance of BC-based LCC in different traffic conditions; (ii) stop-and-go traffic is not necessary for DRL-based LCC to generalize to different traffic conditions; (iii) DRL-based LCC trained with only constant-speed lead vehicle trajectories (not sufficient to ensure identifiability) can be generalized to different traffic conditions; and (iv) stop-and-go traffic increases variance in the training dataset, which improves the convergence of parameter estimation while negatively impacting the convergence of DRL to the optimal control policy. Numerical experiments validate the above findings, illustrating that BC-based LCC entails comprehensive training datasets for generalizing to different traffic conditions, while DRL-based LCC can achieve generalization with simple free-flow traffic training environments. This further suggests DRL as a more promising and cost-effective LCC approach to reduce operational costs, mitigate traffic congestion, and enhance safety and mobility, which can accelerate the deployment and acceptance of CAVs.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Synergistic epistasis enhances the co-operativity of mutualistic interspecies interactions

Abstract Early evolution of mutualism is characterized by big and predictable adaptive changes, including the specialization of interacting partners, such as through deleterious mutations in genes not required for metabolic cross-feeding. We sought to investigate whether these early mutations improve cooperativity by manifesting in synergistic epistasis between genomes of the mutually interacting species. Specifically, we have characterized evolutionary trajectories of syntrophic interactions of Desulfovibrio vulgaris (Dv) with Methanococcus maripaludis (Mm) by longitudinally monitoring mutations accumulated over 1000 generations of nine independently evolved communities with analysis of the genotypic structure of one community down to the single-cell level. We discovered extensive parallelism across communities despite considerable variance in their evolutionary trajectories and the perseverance within many evolution lines of a rare lineage of Dv that retained sulfate-respiration (SR+) capability, which is not required for metabolic cross-feeding. An in-depth investigation revealed that synergistic epistasis across pairings of Dv and Mm genotypes had enhanced cooperativity within SR− and SR+ assemblages, enabling their coexistence within the same community. Thus, our findings demonstrate that cooperativity of a mutualism can improve through synergistic epistasis between genomes of the interacting species, enabling the coexistence of mutualistic assemblages of generalists and their specialized variants.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Weak-form latent space dynamics identification

Recent work in data-driven modeling has demonstrated that a weak formulation of model equations enhances the noise robustness of a wide range of computational methods. In this paper, we demonstrate the power of the weak form to enhance the LaSDI (Latent Space Dynamics Identification) algorithm, a recently developed data-driven reduced order modeling technique. We introduce a weak form-based version WLaSDI (Weak-form Latent Space Dynamics Identification). WLaSDI first compresses data, then projects onto the test functions and learns the local latent space models. Notably, WLaSDI demonstrates significantly enhanced robustness to noise. With WLaSDI, the local latent space is obtained using weak-form equation learning techniques. Compared to the standard sparse identification of nonlinear dynamics (SINDy) used in LaSDI, the variance reduction of the weak form guarantees a robust and precise latent space recovery, hence allowing for a fast, robust, and accurate simulation. We demonstrate the efficacy of WLaSDI vs. LaSDI on several common benchmark examples including viscid and inviscid Burgers', radial advection, and heat conduction. For instance, in the case of 1D inviscid Burgers' simulations with the addition of up to 100% Gaussian white noise, the relative error remains consistently below 6% for WLaSDI, while it can exceed 10,000% for LaSDI. Similarly, for radial advection simulations, the relative errors stay below 15% for WLaSDI, in stark contrast to the potential errors of up to 10,000% with LaSDI. Moreover, speedups of several orders of magnitude can be obtained with WLaSDI. For example applying WLaSDI to 1D Burgers' yields a 140X speedup compared to the corresponding full order model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Introducing a multiscale feature integration network for inpainting with applications to enhanced CMB map reconstruction

We introduce a novel neural network, SkyReconNet, which combines the expanded receptive fields of dilated convolutional layers along with standard convolutions, to capture both the global and local features for reconstructing the missing information in an image. We implement our network to inpaint the masked regions in a full-sky cosmic microwave background (CMB) map. Inpainting CMB maps is a particularly formidable challenge when dealing with extensive and irregular masks, such as galactic masks which can obscure substantial fractions of the sky. The hybrid design of SkyReconNet leverages the strengths of standard and dilated convolutions to accurately predict CMB fluctuations in the masked regions by effectively utilizing the information from surrounding unmasked areas. During training, the network optimizes its weights by minimizing a composite loss function that combines the structural similarity index measure (SSIM) and mean squared error (MSE). SSIM preserves the essential structural features of the CMB, ensuring an accurate and coherent reconstruction of the missing CMB fluctuations, while MSE minimizes the pixelwise deviations, thus enhancing the overall accuracy of the predictions. The predicted CMB maps and their corresponding angular power spectra align closely with the targets, achieving the performance limited only by the fundamental uncertainty of cosmic variance. The network’s generic architecture enables application to other physics-based challenges involving data with missing or defective pixels, systematic artifacts, etc. In conclusion, our results demonstrate its effectiveness in addressing the challenges posed by large irregular masks, offering a significant inpainting tool not only for CMB analyses but also for image-based experiments across disciplines where such data imperfections are prevalent.

Cosmic microwave background↗
DOE OSTI2022

High-Temperature 2D Optical Relaxation Visualizes Enhanced Oxygen Exchange Kinetics at Metal-Mixed Conducting Oxide Interfaces

Solid-state heterointerfaces are of interest for emergent local behavior that is distinct from either bulk parent compound. One technologically relevant example is the case of mixed ionic/electronic conductor (MIEC)–metal interfaces, which play an important role in electrochemistry. Metal–MIEC composite electrodes can demonstrate improved catalytic activity vs single-phase MIECs, improving fuel cell efficiency. Similarly, MIEC surface reaction kinetics are often evaluated using techniques that place metal current collectors in contact with the surface under evaluation, potentially altering the response vs the native surface. Here, techniques enabling direct and local in situ observation of the behavior at and around such heterointerfaces are needed. Here, we develop a spatially resolved optical transmission relaxation (2D-OTR) method providing continuous evaluation of local, high-temperature, controlled atmosphere defect kinetics across a ~1 cm2 sample area simultaneously in a contact-free manner. We apply it to observe the spatial variance of oxygen incorporation and evolution rates at ~525–620 °C, in response to step changes in oxygen partial pressure, on MIEC SrTi 0.65 Fe 0.35 O 3–x films as a function of distance from porous Pt and Au layers. Using this model geometry, we find significant enhancements in kinetics adjacent to the metals that decay over a few millimeter distance. To extract kinetic parameters, we fit the short-term optical data (initial portion of relaxations) with an exponential decay function appropriate for surface-exchange-limited kinetics, yielding apparent surface exchange coefficients (k chem ) with spatial resolution, decreasing with distance from the metal. To understand the kinetic processes governing the complete (long-term) optical relaxations, we performed COMSOL simulations, which demonstrated that a combination of laterally varying k chem and in-plane diffusion controls the observed kinetics over the full time range. Further support for spatially varying kchem comes from demonstrations of changing surface and bulk chemistry vs distance from the metal–MIEC interface, by X-ray photoelectron and optical absorption spectroscopies, respectively. Although microporous Pt and Au are not excellent electrodes in isolation, both metals exert a synergistic effect on the oxygen surface exchange rate in the presence of the mixed conducting film.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Temporal covariance structure of multi-spectral phenotypes and their predictive ability for end-of-season traits in maize

Vegetation indices (VIs) derived from multi-spectral imaging (MSI) platforms can be used to study properties of crop canopy, providing non-destructive phenotypes that could be used to better understand growth curves throughout the growing season. To investigate the amount of variation present in several VIs and their relationship with important end-of-season traits, genetic and residual (co)variances for VIs, grain yield and moisture were estimated using data collected from maize hybrid trials. The VIs considered were Normalized Difference Vegetation Index (NDVI), Green NDVI, Red Edge NDVI, Soil-Adjusted Vegetation Index, Enhanced Vegetation Index and simple Ratio of Near Infrared to Red (Red) reflectance. Genetic correlations of VIs with grain yield and moisture were used to fit multi-trait models for prediction of end-of-season traits and evaluated using within site/year cross-validation. To explore alternatives to fitting multiple phenotypes from MSI, random regression models with linear splines were fit using data collected in 2016 and 2017. Heritability estimates ranging from (0.10 to 0.82) were observed, indicating that there exists considerable amount of genetic variation in these VIs. Furthermore, strong genetic and residual correlations of the VIs, NDVI and NDRE, with grain yield and moisture were found. Considerable increases in prediction accuracy were observed from the multi-trait model when using NDVI and NDRE as a secondary trait. Finally, random regression with a linear spline function shows potential to be used as an alternative to mixed models to fit VIs from multiple time points.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Inferring Plant Acclimation and Improving Model Generalizability With Differentiable Physics‐Informed Machine Learning of Photosynthesis

Net photosynthesis (A N ) is a key component of the global carbon cycle influencing climate feedback over decadal scales. Although plant acclimation to environmental changes can modify A N , traditional vegetation models in Earth system models (ESMs) often rely on plant functional type (PFT)-specific parameterizations or simplified acclimation assumptions limiting generalizability across time, space, and PFTs. In this study, we developed a differentiable photosynthesis model to learn the environmental dependencies of V c,max25 (maximum carboxylation rate at 25°C, representing photosynthetic capacity), as this genre of hybrid physics-informed machine learning can seamlessly train neural networks and process-based equations together. Compared to PFT-specific parameterization of V c,max25 , learning the environment dependencies of key photosynthetic parameters improved model spatiotemporal generalizability. Applying environmental acclimation to V c,max25 led to substantial variations in global mean A N indicating the need to address acclimation in ESMs. The model effectively captured multivariate observations (V c,max25 , A N , and stomatal conductance (g s )) simultaneously with multivariate constraints, improving generalization across space and PFTs. It also learned sensible acclimation relationships of V c,max25 to different environmental conditions. The model explained more than 54%, 57%, and 62% of the variance of A N , g s , and V c,max25 , respectively, presenting a first global-scale spatial test benchmark of A N and g s . These results highlight the potential for differentiable modeling to enhance process-based modules in ESMs and effectively leverage information from large, multivariate data sets.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

The DESI One-Percent Survey: Modelling the clustering and halo occupation of all four DESI tracers with U CHUU

We present results from a set of mock lightcones for the DESI One-Percent Survey, created from the UCHUU simulation. This 8 h −3 Gpc 3 N-body simulation comprises 2.1 trillion particles and provides high-resolution dark matter (sub)haloes in the framework of the Planck-based ΛCDM cosmology. Employing the subhalo abundance matching (SHAM) technique, we populated the UCHUU (sub)haloes with all four DESI tracers – Bright Galaxy Survey (BGS), luminous red galaxies (LRGs), emission line galaxies (ELGs), and quasars (QSOs) – to z = 2.1. Our method accounts for redshift evolution as well as the clustering dependence on luminosity and stellar mass. The two-point clustering statistics of the DESI One-Percent Survey generally agree with predictions from UCHUU across scales ranging from 0.3 h −1 Mpc to 100 h −1 Mpc for the BGS and across scales ranging from 5 h −1 Mpc to 100 h −1 Mpc for the other tracers. We observed some differences in clustering statistics that can be attributed to incompleteness of the massive end of the stellar mass function of LRGs, our use of a simplified galaxy-halo connection model for ELGs and QSOs, and cosmic variance. We find that at the high precision of UCHUU, the shape of the halo occupation distribution (HOD) of the BGS and LRG samples is smaller bias values, likely due to cosmic variance. The bias dependence on absolute magnitude, stellar mass, and redshift aligns with that of previous surveys. These results provide DESI with tools to generate high-fidelity lightcones for the remainder of the survey and enhance our understanding of the galaxy-halo connection.

cosmology↗
DOE OSTI2023

Spin Squeezing by Rydberg Dressing in an Array of Atomic Ensembles

Here, we report on the creation of an array of spin-squeezed ensembles of cesium atoms via Rydberg dressing, a technique that offers optical control over local interactions between neutral atoms. We optimize the coherence of the interactions by a stroboscopic dressing sequence that suppresses super-Poissonian loss. We thereby prepare squeezed states of N = 200 atoms with a metrological squeezing parameter ξ 2 = 0.77⁢(9) quantifying the reduction in phase variance below the standard quantum limit. We realize metrological gain across three spatially separated ensembles in parallel, with the strength of squeezing controlled by the local intensity of the dressing light. Our method can be applied to enhance the precision of tests of fundamental physics based on arrays of atomic clocks and to enable quantum-enhanced imaging of electromagnetic fields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A Machine Learning–Based Tire Life Prediction Framework for Increasing Life of Commercial Vehicle Tires

In the commercial freight industry, tire retreading decisions are often conservative due to limited knowledge of a tire’s remaining service life. This practice leads to increased costs and material waste. This paper proposes a machine learning–based approach for estimating tire casing life and retreadability, focusing on usage data rather than wear information. This approach could extend the tire’s lifespan and reduce landfill waste. Data integration from diverse tire casing measurement sources presents challenges, including imbalanced removal data. Our methodology addresses these challenges by using historical inspection, telematics, and finite element modeling (FEM) datasets. We introduce “Tire Casing Energy” as a comprehensive usage input and apply a Variance-Reduction Synthetic Minority Oversampling Technique (VR-SMOTE) for data imbalance rectification. A random forest model is used to estimate the state of the tire casing and the casing removal probability, with Bayesian optimization applied for hyperparameter tuning, enhancing model accuracy. Here, the proposed prediction framework is able to differentiate different truck fleets and tire locations based on their usage parameters. With the aid of this machine learning model, the importance and sensitivity of different tire usage parameters can be obtained, which is beneficial to maximize tire life.

Data balancing↗
DOE OSTI2020

Dust Growth by Accretion of Molecules in Supersonic Interstellar Turbulence

We show that the growth rate of dust grains in cold molecular clouds is enhanced by the high degree of compressibility of a turbulent, dilute gas. By means of high-resolution (1024 3 ) numerical simulations, we confirm the theory that the spatial mean growth rate is proportional to the gas-density variance. This also results in broadening of the grain-size distribution (GSD) due to turbulence-induced variation of the grain-growth rate. We show, for the first time in a detailed numerical simulation of hydrodynamic turbulence, that the GSD evolves toward a shape that is a reflection of the gas-density distribution, regardless of the initial distribution. That is, in case of isothermal, rotationally forced turbulence, the GSD tends to be a lognormal distribution. We also show that in hypersonic turbulence, decoupling of gas and dust becomes important and that this leads to an even further accelerated grain growth.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

APSO-enhanced algebraic derivative estimation approach for real-time traffic flow prediction on critical road sections during wildfire evacuation

In rapid-onset disaster scenarios such as wildfires, evacuation traffic often significantly deviates from historical patterns, rendering conventional data-driven forecasting methods less effective. To address this challenge, we propose an improved algebraic derivative estimation (ADE) incorporating particle swarm optimization (PSO) for real-time traffic flow prediction. Our approach dynamically adjusts the ADE prediction time window at each step by minimizing a cost function based on the mean and variance of accumulated forecasting errors within the window, thereby balancing bias and variability. We evaluate the method using traffic data from the January 2025 California wildfires, focusing on key road segments critical for large-scale evacuations. The results demonstrate that our approach surpasses established machine learning and deep learning models—XGBoost, LSTM, and GRU—in predictive accuracy and maintains high computational efficiency. Notably, the proposed method eliminates the need for offline model training. Moreover, rapid PSO-based tuning enables real-time deployment, which provides a crucial advantage in scenarios where evacuation timings and road closures change dynamically. In conclusion, these findings highlight the benefits of the PSO-enhanced ADE framework for emergency traffic management, where rapid, data-sparse forecasts are essential for effective evacuation planning.

Algebraic derivative estimation↗
DOE OSTI2023

Temperature-dependent solvent reorganization entropies, free energies, and transition dipole strengths for the photoexcitation of Reichardt’s dye B30

Absorption spectra of the solvatochromic dye 2,6-diphenyl-4-2,4,6-triphenyl-1-pyridinophenolate (B30) were measured in seven solvents of varying polarity over temperature ranging from each solvent’s freezing point to 300 K. The excitation energies and their variances allowed calculations of the solvent reorganization energies, reorganization free energies and reorganization entropies as functions of temperature. The entropies of solvent packing around the chromophore are found to make major contributions to the reorganization free energies. The variances of the excitation energies depend only weakly on temperature, in disagreement with an expression that is often used for solvent reorganization free energies. Polar solvents reduce the transition dipole strength of B30’s long-wavelength absorption band, probably because interactions with the solvent enhance the charge-transfer character of the transition. The dipole strength drops further at low temperatures.

Chemistry↗
DOE OSTI2025

PQML: Enabling the Predictive Reproducibility on NISQ Machines for Quantum ML Applications

Quantum computing represents a groundbreaking approach to high-performance computing. In recent years, quantum computers have progressed from single-qubit processors to systems boasting over 400 qubits. The presence of such a large number of qubits offers significant advantages, including enhanced computational speed—a capability beyond classical computing methods. However, the current stage of quantum computing is referred to as the noisy intermediate-scale quantum (NISQ) era. The existence of noise in this era presents challenges in testing quantum computing applications, leading to considerable variance in application results. Furthermore, the diverse noise characteristics observed across different machines exacerbate this issue, complicating the selection of the appropriate machine for application execution. In response to these challenges, we introduce our Predictive Quantum Machine Learning (PQML) tool. This tool is designed to predict outcomes when executing identical quantum machine learning applications—specifically, a critical suite of variational quantum algorithms—across various quantum computers during the NISQ era. This effort relies on data collected over a 12-month period. To the best of our knowledge, this study represents the first attempt to ensure reproducibility across quantum computers for complex circuits. Additionally, we have developed a model capable of forecasting the accuracy of quantum computers for variational quantum algorithms, with a particular emphasis on quantum machine learning as a case study.

Senapati, Priyabrata [Kent State University]↗
DOE OSTI2025

Multi-qubit nanoscale sensing with entanglement as a resource

Nitrogen vacancy (NV) centres in diamond are widely deployed as local magnetic sensors, using single-qubit control to measure both time-averaged fields and noise with nanoscale spatial resolution. Moving beyond single qubits to multi-qubit control enables new sensing modalities such as measuring nonlocal spatiotemporal correlators or using entangled states to enhance measurement sensitivity. Here, in this work, we describe protocols for using optically unresolved NV centre pairs and nuclear spins as multi-qubit sensors for measuring correlated noise at nanometre length scales. For noninteracting NV centres, we implement a phase-cycling protocol that disambiguates magnetic correlations from variance fluctuations, leveraging the presence of a third qubit, a 13 C nucleus, to effect coherent single-NV spin flips and enable phase cycling even for co-aligned NV centres that are spectrally unresolved. For length scales around 10 nm, we create maximally entangled Bell states through dipole–dipole coupling between two NV centres and use these entangled states to directly read out the magnetic field correlation, rather than reconstructing it from independent measurements of unentangled NV centres. Importantly, this changes the scaling of sensitivity with readout noise from quadratic to linear. For conventional off-resonant readout of the NV centre spin state (for which the readout noise is roughly 30 times the quantum projection limit), this results in more than an order of magnitude improvement in sensitivity. Finally, we demonstrate methods for detecting high spatial- and temporal-resolution correlators with pairs of strongly interacting NV centres.

condensed-matter physics↗
DOE OSTI2026

Optimizing aluminum oxide passivation layers—Laser-induced plasmas for layer formation, real-time diagnostics, and layer evaluation

Molten salts have beneficial thermophysical and electrolytic properties, but the aggressive nature of molten salts requires corrosion mitigation strategies, such as structural material surface treatment. In this study, a laser-induced breakdown spectroscopy (LIBS) system was used to simultaneously ablate and monitor de-excitation spectra to form and identify dense alumina passivation layers. Additionally, high-frequency (kHz) LIBS imaging allowed rapid elemental mapping and depth profiling for surface O/Al ratios via minimally destructive analysis. The paired t-test rejected the null hypothesis that air and Ar cover gases were equivalent (p = 0.00012), and the analysis of variance (ANOVA) showed that overlap between shots (p = 0.03112) and cover gas flow rate (p = 0.02729) were statistically significant. At 75% overlap, the O/Al ratio increased to 10% ± 2%, and 90% overlap rose further to 31% ± 2%. Regardless of overlap, tuning the gas flow from 0.5 to 2.5 L min −1 enhanced the O/Al ratio by 24% ± 4%. LIBS depth profiling showed a 3× increase in layer thickness from 75% to 90% overlap. Molecular band peaks of AlO, a precursor to Al 2 O 3 found during treatment, suggest that band intensity could be used for real-time Al 2 O 3 optimization. Ultimately, by adjusting laser spot overlap and cover gas flow rate, a 266 nm laser was shown to form Al 2 O 3 layers on an Al6061 sample, with cover gas flow rate primarily affecting the O/Al ratio and layer thickness correlating with laser spot overlap. This combination of molecular spectra collection, imaging, and depth profiling demonstrated the robust layer analysis capabilities of LIBS.

Corrosion↗
DOE OSTI2022

Examining the changes in the spatial manifestation and the rate of arrival of large tornado outbreaks

This study presents an assessment of the spatial and temporal characteristics of large tornado outbreak (LTOs) days, in which several counties were impacted by tornadoes rated F2(EF2) or greater on the Fujita (Enhanced Fujita) scale in one day. A statistical evaluation of changes in the LTO clusters for two periods, 1950–1980 and 1989–2019, has been performed. There is a geographical shift of the nucleus (central impact location) towards the southeast United States. This spatial shift is also accompanied by reduced spatial variance, suggesting LTOs have become less dispersed (or more localized) in the recent period. The overall inter-arrival rate of LTOs, and how it changed during successive 31-year climatological blocks between 1950–2019 was investigated using an exponential probability model. The arrival rate has changed from 124 days during 1950–1980 to 164 days during 1977–2007 and remained relatively constant during later periods, indicating that LTOs are becoming less frequent.

54 ENVIRONMENTAL SCIENCES↗