Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Robust stability”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Safe Physics-Informed Machine Learning for Dynamics and Control

This tutorial paper focuses on safe physics-informed machine learning in the context of dynamics and control, providing a comprehensive overview of how to integrate physical models and safety guarantees. As machine learning techniques enhance the modeling and control of complex dynamical systems, ensuring safety and stability remains a critical challenge, especially in safety-critical applications like autonomous vehicles, robotics, medical decision-making, and energy systems. We explore various approaches for embedding and ensuring safety constraints, including structural priors, Lyapunov and Control Barrier Functions, predictive control, projections, and robust optimization techniques. Additionally, we delve into methods for uncertainty quantification and safety verification, including reachability analysis and neural network verification tools, which help validate that control policies remain within safe operating bounds even in uncertain environments. The paper includes illustrative examples demonstrating the implementation aspects of safe learning frameworks that combine the strengths of data-driven approaches with the rigor of physical principles, offering a path toward the safe control of complex dynamical systems.

Drgona, Jan↗
DOE OSTI2021

Stress- and interface-compatible red phosphorus anode for high-energy and durable sodium-ion batteries

Sodium-ion battery is a promising candidate for energy storage application, but the absence of high-capacity and low-cost anode materials significantly limit their practical specific energy and cost. Red phosphorus (RP) possesses a high theoretical specific capacity but suffers from large volume change, low electronic conductivity and unstable solid-electrolyte interphase (SEI). Herein, a hierarchical micro/nanostructured antimony-doped RP/carbon anode was developed, which demonstrate extraordinary electrochemical performance with high initial Coulombic efficiency of ~90%, high areal (~1.7 mAh cm-2) capacity together with good cycle stability and rate capability. Combined experimental and computational studies consistently revealed that such a unique structural design can dramatically accommodate the mechanical stress and moreover effectively restrain the undesired decomposition of electrolyte solvents regardless of electrolytes formulation, resulting in superior structural integrity and thin/robust SEI formation during cycling. The present finding has offered an alternative strategy for stress management and interface engineering on high-capacity alloying-based anode materials.

Liu, Xiang↗
DOE OSTI2021

Bilayer Wigner crystals in a transition metal dichalcogenide heterostructure

One of the first theoretically predicted manifestations of strong interactions in many-electron systems was the Wigner crystal, in which electrons crystallize into a regular lattice. The crystal can melt via either thermal or quantum fluctuations. Quantum melting of the Wigner crystal is predicted to produce exotic intermediate phases and quantum magnetism because of the intricate interplay of Coulomb interactions and kinetic energy. However, studying two-dimensional Wigner crystals in the quantum regime has often required a strong magnetic field or a moiré superlattice potential, thus limiting access to the full phase diagram of the interacting electron liquid. Here we report the observation of bilayer Wigner crystals without magnetic fields or moiré potentials in an atomically thin transition metal dichalcogenide heterostructure, which consists of two MoSe 2 monolayers separated by hexagonal boron nitride. We observe optical signatures of robust correlated insulating states at symmetric (1:1) and asymmetric (3:1, 4:1 and 7:1) electron doping of the two MoSe 2 layers at cryogenic temperatures. Here, we attribute these features to bilayer Wigner crystals composed of two interlocked commensurate triangular electron lattices, stabilized by inter-layer interaction. The Wigner crystal phases are remarkably stable, and undergo quantum and thermal melting transitions at electron densities of up to 6 × 10 12 per square centimetre and at temperatures of up to about 40 kelvin. Our results demonstrate that an atomically thin heterostructure is a highly tunable platform for realizing many-body electronic states and probing their liquid–solid and magnetic quantum phase transitions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modeling of thermal-ion-driven internal kink in DIII-D high- T i plasmas

Toroidal modeling utilizing the non-perturbative magnetohydrodynamic-kinetic hybrid code MARS-K (Liu Y.Q. et al 2008 Phys. Plasmas 15 112503), to a DIII-D high-performance discharge with the core thermal ion temperature of T i ~ 13 keV, finds multiple branches of unstable internal kink (IK) modes triggered by the bounce resonance of trapped thermal ions. The range of the computed mode frequency aligns well with the bounce frequency of trapped thermal ions, confirming the nature of resonance drive for the instability. The computed mode frequency is close to that observed in experiments. Plasma toroidal flow, in particular the fast flow ($ > $10% of Alfvén speed) as in the considered DIII-D discharge, strongly stabilizes the fluid IK branch but not all of the thermal-ion-driven internal kink (TIK) branches, leaving one dominant unstable TIK branch under the experimental conditions. Plasma resistivity also significantly affects stability of the TIKs (but not much on the IK which is ideally strongly unstable) in combination with the drift kinetic drive by thermal ions. High resistivity can drive additional unstable TIK branches. Scanning the on-axis safety factor over a large range confirms robustness of the TIK triggering by thermal ions. The computed TIKs contain both the IK and the tearing components as the dominant contributions to the eigenmode structure. Furthermore, these findings, together with the experimental results reported in reference (Du X.D. et al 2021 Phys. Rev. Lett. 127 025001), strongly suggest the possibility of thermal-ion driven instabilities in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Efficient and Robust Dynamic Simulation of Power Systems With Holomorphic Embedding

Dynamic simulation is vitally important in power system analysis, but traditional approaches based on numerical integration over small time steps are time-consuming. Also, the Newton-Raphson method suffers from difficulty in convergence when solving nonlinear algebraic equations. In this paper, we propose a novel dynamic simulation approach based on holomorphic embedding. By obtaining a high-order approximation of system dynamics, it achieves a much larger time step and thus enhances the computational efficiency significantly. In addition, the new approach avoids non-convergence issues in solving algebraic equations, which improves robustness. The approach includes flexible modeling of synchronous generators and controllers, and we propose a method for modeling generator coordinate transformations. The approach is tested on the IEEE 39-bus, 10-generator system and a Polish 2383-bus, 327-generator system. The results demonstrate promising computational efficiency and satisfactory numerical robustness for the analysis of large-scale power systems.

42 ENGINEERING↗
DOE OSTI2023

A geospatial risk analysis graphical user interface for identifying hazardous chemical emission sources

Background: Performing back trajectory and forward trajectory using the Hybrid Single-Particle Lagrangian Integrated Trajectory Model (HYSPLIT) is a reliable approach for assessing particle transport after release among mid-field atmospheric models. HYSPLIT has an externally facing online interface that allows non-expert users to run the model trajectories without requiring extensive training or programming. However, the existing HYSPLIT interface is limited if simulations have a large amount of meteorological data and timesteps that are not coincident. The objective of this study is to design and develop a more robust tool to rapidly evaluate hazard transport conditions and to perform risk analysis, while still maintaining an intuitive and user-friendly interface. Methods: HYSPLIT calculates forward and backward trajectories of particles based on wind speed, wind direction, and the corresponding location, timestamp, and Pasquill stability classes of the regions of the atmosphere in terms of the wind speed, the amount of solar radiation, and the fractional cloud cover. The computed particle transport trajectories, combined with the online Proton Transfer Reaction-Mass Spectrometry (PTR-MS) data (https://figshare.com/articles/dataset/ARL_Data_from_PROS_station_at_Hanford_site/19993964), can be used to identify and quantify the sources and affected area of the hazardous chemicals’ emission using the potential source distribution function (PSDF). PSDF is an improved statistical function based on the well-known potential source contribution function (PSCF) in establishing the air pollutant source and receptor relationship. Performing this analysis requires a range of meteorological and pollutant concentration measurements to be statistically meaningful. The existing HYSPLIT graphical user interface (GUI) does not easily permit computations of trajectories of a dataset of meteorological data in high temporal frequency. To improve the performance of HYSPLIT computations from a large dataset and enhance risk analysis of the accidental release of material at risk, a geospatial risk analysis tool (GRAT-GUI) is created to allow large data sets to be processed instantaneously and to provide ease of visualization. Results: The GRAT-GUI is a native desktop-based application and can be run in any Windows 10 system without any internet access requirements, thus providing a secure way to process large meteorological datasets even on a standalone computer. GRAT-GUI has features to import, integrate, and convert meteorological data with various formats for hazardous chemical emission source identification and risk analysis as a self-explanatory user interface. The tool is available at https://figshare.com/articles/software/GRAT/19426742.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

An extension of the localized artificial diffusivity method for immiscible and high density ratio flows

The localized artificial diffusivity (LAD) method is widely regarded as the preferred multi-material regularization scheme for the compact finite difference method, because it is conservative, easy to implement, and generally robust for a wide range of multi-material problems. However, traditional LAD methods face significant challenges when applied to flows with large density ratios and when maintaining thermodynamic equilibrium across material interfaces. These limitations arise from the formulation of the artificial diffusivity flux and the reliance on enthalpy diffusion for interface regularization. Additionally, traditional LAD methods struggle to ensure stability under large density ratio conditions, fail to maintain a finite interface thickness, and are therefore unsuitable for modeling immiscible interfaces. Here, in this work, we discuss the origins of these issues in traditional LAD methods and propose modifications which enable the simulation of large density ratio and immiscible flows. The proposed method targets the artificial diffusion fluxes at gradients and ringing in the volume fraction, rather than the mass fraction in traditional methods, to consistently regularize large density ratio interfaces. Furthermore, the proposed method introduces an artificial bulk density diffusion term to enforce equilibrium conditions across interfaces. To address the challenge of modeling immiscible flows, a conservative diffuse interface term is incorporated into the formulation to ensure a finite interface thickness. Specific consideration is taken in the design of the method to ensure that these crucial properties are maintained for N -material flows. The effectiveness of the proposed method is demonstrated through a series of canonical test cases, and its accuracy is validated by comparison with experimental data on micro-bubble collapse in water. These results highlight the method’s robustness and its ability to overcome the limitations of traditional LAD approaches.

Artificial diffusivity↗
DOE OSTI2025

Unified Universal Control and Coordination of Inverter-Based Resources, and Validation for a PV + Battery Hybrid Plant

As renewable energy deployment grows, hybrid power plants (HPPs) combining photovoltaic (PV) and battery systems must evolve to offer both energy and grid stability services. These systems typically include a mix of grid-following (GFL) and grid-forming (GFM) inverters, presenting unique coordination and control challenges. This Department of Energy–funded project developed and validated a Unified Universal Control and Coordination (UUCC) framework for such PV + battery hybrid plants, enabling seamless and stable operation, including ultrafast black start, autonomous synchronization, and robust frequency and voltage regulation, under different grid conditions. The project significantly advanced the understanding of inverter-based resource (IBR) control by developing and validating three complementary system-level approaches for hybrid GFL/GFM operation: 1. A combined Virtual Resistance (VR)-based GFL and Virtual Oscillator Control (VOC)-based GFM method, where each inverter type is governed by a specialized control strategy. Together, these achieve stable, fast-response coordination, eliminating inrush current and enabling smooth black start and grid synchronization across a wide range of grid strengths. 2. A Deadbeat-based UUCC strategy, which uses discrete-time, switching-cycle-level control for both GFL and GFM inverters. This approach replaces traditional PI/PLL control with a control parameter-free, high-bandwidth framework that supports stable LVRT and instantaneous synchronization under all conditions. 3. A benchmark comparison with Siemens’ commercial GFM microgrid controller, which provided a fast baseline platform. The commercial approach decoupled v & f control was implemented on a commercial microgrid controller.The baseline commercial benchmark helped highlight superior transient response and black start performance offered by the deadbeat and VOC approaches. These technical contributions offer substantial improvements over conventional inverter control schemes, which often rely on slow phase-locked loop (PLL)-based synchronization, require careful control parameters tuning, and prone to unstable in weak grids with GFL inverters and in stiff grid with GFM inverters therefore challenging for hybrid GFL+GFM under all grid conditions. The deadbeat-based UUCC framework enables simpler, faster, and more robust operation of hybrid IBR systems using wide-bandgap (WBG) devices such as SiC power semiconductors. The rapid expansion of hybrid distributed energy resources (DERs), including residential and commercial PV-BESS installations such as Tesla Powerwall, PV with vehicle-to-grid (V2G) capability, and other integrated configurations, presents complex operational challenges for medium-voltage radial distribution feeders. These networks are subject to frequent disturbances such as faults, switching operations, rapid reclosing sequences, and feeder reconfigurations, all of which introduce dynamic stress on IBRs. In addition, planned feeder segmentation and deliberate islanding for resilience will require DERs that can autonomously perform blackstart, establish voltage and frequency references, and resynchronize with the main grid. The advanced deadbeat-based UUCC control and blackstart functionalities developed in this project directly address these requirements, enabling decentralized and autonomous operation of inverter-dominated DERs in distribution systems under a wide range of fault and reconfiguration scenarios. From a public benefit perspective, these innovations enable more reliable and cost-effective integration of renewable energy into distribution networks. The ability to autonomously black start and stabilize grids under varying grid conditions support accelerates recovery from outages and support decentralized resilient energy systems. By reducing system complexity and improving performance, this project lays critical groundwork for future inverter-dominated power grids that are clean, reliable, and accessible to all.

14 SOLAR ENERGY↗
DOE OSTI2023

Solvent-Mediated Formation of Quasi-2D Dion-Jacobson Phases on 3D Perovskites for Inverted Solar Cells Over 23% Efficiency

2D-on-3D (2D/3D) perovskite heterostructures present a promising strategy to realize efficient and stable photovoltaics. However, their applicability in inverted solar cells is limited due to the quantum confinement of the 2D-layer and solvent incompatibilities that disrupt the underlying 3D layer, hampering electron transport at the 2D/3D interface. Herein, solvent-dependent formation dynamics and structural evolution of 2D/3D heterostructures are investigated via in situ X-ray scattering. In this study, it is revealed that solvent interaction with the 3D surface determines the formation sequence and spatial distribution of quasi-2D phases with n = 2–4. Isopropanol (IPA) reconstructs the perovskite into a PbI 2 -rich surface, forming a strata with smaller n first, followed by a thinner substratum of larger n. In contrast, 2,2,2-Trifluoroethanol (TFE) preserves the 3D surface, promoting the formation of uniformly distributed larger n domains first, and smaller n last. Leveraging these insights, Dion–Jacobson perovskites are used with superior charge transport properties and structural robustness to fabricate 2D/3D heterostructures dominated by n ≥ 3 and engineer a favorable energy landscape for electron tunneling. Inverted solar cells based on 3-Aminomethylpyridine and TFE achieve a champion efficiency of 23.60%, with V oc and FF of 1.19 V and 84.5%, respectively, and superior stabilities with t 94 of 960 h under thermal stress.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

How non-ohmic contact-layer diodes in perovskite pinholes affect abrupt low-voltage reverse-bias breakdown and destruction of solar cells

Perovskite solar cells (PSCs) rapidly degrade under reverse bias, a condition that may occur during partial shading. Here, in this study, we use electrical measurements, electron microscopy, and optical and thermal imaging to investigate abrupt breakdown and hotspotting under low reverse potentials (<|-2| V). We show that microscopic pinholes in the perovskite layer cause rapid, destructive breakdown under reverse bias despite minimally reducing power conversion efficiencies. Measurements on miniature (200-micrometer diameter) PSCs and perovskite-free transport-layer diodes indicate that abrupt, low-voltage breakdown occurs in nanoscale to micrometer-scale defects and that metal migration and filamentation are unlikely causes. Reverse-bias stability substantially improves when pinholes in the perovskite and transport layers are eliminated. Atomic layer deposition of tin oxide prevents abrupt breakdown by ensuring physical separation between electrodes-not by blocking metal ion migration. Perovskite researchers should adopt cleaner, more uniform deposition techniques to enable robust PSCs for further research and commercial applications.

14 SOLAR ENERGY↗
DOE OSTI2022

Salton Sea Geothermal Development: Nontechnical Barriers to Entry – Analysis and Perspectives

Geothermal energy offers an opportunity to generate baseload, renewable energy that can help support the transition to an energy economy with reduced impacts on climate change and replace older, more expensive, nonrenewable, and more resource-impacting energy-generation facilities. The United States has the largest known geothermal resource in the world, with over 31 GW of conventional geothermal potential. However, due to market conditions, an inability to properly quantify both electrical grid benefits and resource stability, and the difficulty of exploring and developing the geothermal resource, few new geothermal projects have come online over the past three decades. The Salton Sea, in Imperial County, California, provides a prime location and opportunity to develop new geothermal resources. The Salton Sea contains a robust, well-mapped, geothermal resource, with opportunities for concurrent development of lithium and other mineral resources. This report describes the history of geothermal development at the Salton Sea and compares geothermal to other renewable energy sources in the area. The report then uses a techno-economic analysis (TEA) model to analyze the relative benefits and costs of various challenges and opportunities and provides recommendations for streamlining geothermal development at the Salton Sea and elsewhere. The challenges and opportunities analyzed in the TEA model were informed by stakeholder interviews and literature reviews. Based upon the identified challenges and opportunities and the results of the TEA model, primary findings are that certain nontechnical barriers such as permitting costs play only a minor role in determining the viability of development of the geothermal resource at the Salton Sea. Other barriers such as permitting timelines, government/agency coordination, and the potential co-location of lithium extraction with a geothermal plant may result in much larger impacts on project viability.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Accelerated Selectrion of Optimal Perovskite Alloys for Solar PV using a Combined Quantum and Machine Learning Hierachiral Approach

The project aims to: (i) accelerate the discovery of “Missing HP alloys” by combining quantum mechanics and artificial intelligence machine learning approaches, and (ii) analyze the stabilities of candidate alloys, including those that do not pass selection filters (and are thus expected to degrade over time) to decipher the nature of the instabilities to guide the development of durable solar cell materials. Successful candidates will be subjected to validation experiments at NREL's state-of-the-art facilities. The discoveries this effort will provide will be directly testable and implementable and will greatly impact U.S. progress in HP PV as they will provide clear direction and motivation for experimental studies including specific material synthetic targets, device optimization, and device stability protocols. A key advantage of this effort is the feedback and guidance provided by the Industry Collaborative Work Group that we established to coordinate academic and national lab research with industry needs. The proposed work will provide a basis for and direct the development of robust and reliable HP PV. It will also provide a timely, valuable and extensive roadmap to the experimental HP PV community to enable it to focus its efforts on improving and fine-tuning promising HP compositions that this effort predicts will likely be the best performers rather than wandering in the vast chemical space for decades spending enormous resources mostly evaluating unpromising candidate materials.

14 SOLAR ENERGY↗
DOE OSTI2026

MULTI-FIDELITY MODELING AND UNCERTAINTY QUANTIFICATION OF INVERTER BASED RESOURCES IN INTEGRATED T&D SYSTEMS

Uncertainty quantification plays a pivotal role in improving the accuracy and reliability of inverter operation within modern power systems that are increasingly dominated by inverter-based resources (IBRs). IBRs, especially those operating under grid forming (GFM) control, rely heavily on a complex set of control parameters and system measurements to maintain voltage, frequency, and power balance. Traditional deterministic modeling approaches often fail to capture these parameter deviations, potentially resulting in suboptimal control actions, reduced system stability, or even instability under high penetration of IBRs. In this paper, we demonstrate the application of model calibration and uncertainty quantification (UQ) principles to an integrated transmission and distribution (T&D) model involving a GFM converter and provide a framework for prioritizing control improvements, guiding robust design, and informing adaptive strategies that can accommodate real-time variability in system conditions. The proposed approach could be valuable in enhancing the robustness of current and future power systems under increased IBR penetrations.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Artificial cathode electrolyte interphase for improving high voltage cycling stability of thick electrode with Co-free 5 V spinel oxides

Spinel-type cathode LiNi 0.5 Mn 1.5 O 4 (LNMO) has intrigued the transportation industry due to its high operating voltage and total elimination of the expensive cobalt element. However, LNMO cathode with high mass loading (> 3 mAh/cm 2 in areal capacity) has suffered from excessive capacity degradation upon long cycling. Here, a robust Al 2 O 3 surface layer is introduced to the thick LNMO electrode via atomic layer deposition (ALD). The capacity retention in full cells with the graphite anode is improved from 46.3% to 75.3% after 300 cycles with cutoff voltage up to 4.85 V, while enabling average Coulombic efficiency of 99.9% during the cycling. The post-mortem analyses reveal that the Al 2 O 3 surface layer would convert to Al-O-F /Al-F species upon cycling, offering stable interphase to protect the cathode material. Furthermore, these results demonstrate the significance of surface modification enabling high voltage cathode for next-generation LIBs.

25 ENERGY STORAGE↗
DOE OSTI2021

Development and experimental qualification of novel disruption prevention techniques on DIII-D

Novel disruption prevention solutions spanning a range of control regimes are being developed and tested on DIII-D to enable ITER success. First, a new real-time control algorithm has been developed and tested for regulating nearness to stability limits and maintaining safety-margins. Its first application has been for reliable prevention of vertical displacement events (VDEs) by adjusting plasma elongation (κ) and the inner-gap between the plasma and inner-wall in response to real- time open-loop VDE growth rate (γ) estimators. VDEs were robustly prevented up to average open-loop growth rates of 800 rad/s with initial tunings, with only applying shape modification when near safety limits. Second, the disruption risk during fast, emergency shutdown after large tearing and locked modes can be significantly improved by transitioning to a limited topology during shutdown. More than 50% of emergency limited shutdowns after locked modes reach a final normalized current I_N < 0.3 before terminating, scaling to the 3 MA ITER requirement. Furthermore, this is in contrast to diverted shutdowns, the majority of which disrupt at I_N > 0.8. Despite improvements, these results highlight the critical importance of early prevention. Third, a novel emergency shut down method has been developed which excites instabilities to form a warm, helical core post-thermal quench. The current quench extends to ~100ms and avoids VDEs and runaway electron generation. Novel real-time machine learning disruption prediction has been integrated with the DIII-D proximity controller, and a real- time compatible multi-mode MHD spectroscopy technique has been developed. Results presented here were enabled by a focused effort, the Disruption Free Protocol, in DIII-D’s 2019-20 campaign to complement disruption prevention experiments with a large piggy-back program. In addition to testing novel techniques, it is estimated to have helped avoid 32 potential disruptions in piggyback operations with rapid, early shutdowns after large rotating n=1 or locked modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Monitoring leaf phenology in moist tropical forests by applying a superpixel-based deep learning method to time-series images of tree canopies

Tropical leaf phenology-particularly its variability at the tree-crown scale-dominates the seasonality of carbon and water fluxes. However, given enormous species diversity, accurate means of monitoring leaf phenology in tropical forests is still lacking. Time series of the Green Chromatic Coordinate (GCC) metric derived from tower-based red-green-blue (RGB) phenocams have been widely used to monitor leaf phenology in temperate forests, but its application in the tropics remains problematic. To improve monitoring of tropical phenology, we explored the use of a deep learning model (i.e. superpixel-based Residual Networks 50, SP-ResNet50) to automatically differentiate leaves from non-leaves in phenocam images and to derive leaf fraction at the tree-crown scale. To evaluate our model, we used a year of data from six phenocams in two contrasting forests in Panama. Here, we first built a comprehensive library of leaf and non-leaf pixels across various acquisition times, exposure conditions and specific phenocams. We then divided this library into training and testing components. We evaluated the model at three levels: 1) superpixel level with a testing set, 2) crown level by comparing the model-derived leaf fractions with those derived using image-specific supervised classification, and 3) temporally using all daily images to assess the diurnal stability of the model-derived leaf fraction. Finally, we compared the model-derived leaf fraction phenology with leaf phenology derived from GCC. Our results show that: 1) the SP-ResNet50 model accurately differentiates leaves from non-leaves (overall accuracy of 93%) and is robust across all three levels of evaluations; 2) the model accurately quantifies leaf fraction phenology across tree-crowns and forest ecosystems; and 3) the combined use of leaf fraction and GCC helps infer the timing of leaf emergence, maturation and senescence, critical information for modeling photosynthetic seasonality of tropical forests. Collectively, this study offers an improved means for automated tropical phenology monitoring using phenocams.

54 ENVIRONMENTAL SCIENCES↗