Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “range prediction”

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 451 records · Page 25

Data from: "Lab and Field Warming Similarly Advance Germination Date and Limit Germination Rate for High and Low Elevation Provenances of Two Widespread Subalpine Conifers"

This archive contains data used to draw conclusions in “Lab and Field Warming Similarly Advance Germination Date and Limit Germination Rate for High and Low Elevation Provenances of Two Widespread Subalpine Conifers”, by Kueppers et al. 2017. Cone collection and field experiments were conducted on Niwot Ridge, Colorado, with the latter completed within the Alpine Treeline Warming Experiment (ATWE). Laboratory germination experiments were conducted at the University of California Merced. Geospatial files are also provided for additional geographical context.There are four data file formats in this package: comma-separated values (.csv), Microsoft Excel (.xlsx), keyhole markup language (.kml), and ESRI shapefile (.shp). The latter two are geospatial file formats. .csv files can be opened using any simple text editor program such as TextEdit (Mac) and Notepad (Windows), and .xlsx files can be opened using Microsoft Excel or Google Sheets. The .kml file contains coordinate data, and is compatible with Google Earth and Google Maps. ESRI shapefiles (.shp) are compatible with any geospatial software able to read the file type, such as QGIS and ESRI’s ArcGIS Suite. This data user’s guide is available in as a .pdf, and can be opened using software such as Adobe Acrobat Reader.-----------------------------------------------------------------------------------------------------------------------------------Accurately predicting upslope shifts in subalpine tree ranges with warming requires understanding how future forest populations will be affected by climate change, as these are the seed sources for new tree line and alpine populations. Early life history stages are particularly sensitive to climate and are also influenced by genetic variation among populations. We tested the climate sensitivity of germination and initial development for two widely distributed subalpine conifers, using controlled-environment growth chambers with one temperature regime from subalpine forest in the Colorado Rocky Mountains and one 5 °C warmer, and two soil moisture levels. We tracked germination rate and timing, rate of seedling development, and seedling morphology for two seed provenances separated by ~300 m elevation. Warming advanced germination timing and initial seedling development by a total of ~2 weeks, advances comparable to mean differences between provenances. Advances were similar for both provenances and species; however, warming reduced the overall germination rate, as did low soil moisture, only for Picea engelmannii. A three-year field warming and watering experiment planted with the same species and provenances yielded responses qualitatively consistent with the lab trials. Together these experiments indicate that in a warmer, drier climate, P. engelmannii germination, and thus regeneration, could decline, which could lead to declining subalpine forest populations, while Pinus flexilis forest populations could remain robust as a seed source for upslope range shifts.

54 ENVIRONMENTAL SCIENCES↗

Machine Learning Software for Cylindrical Battery Design and Performance Prediction

A software that delivers optimal design parameters and performance predictions for cylindrical cells, which can range in size from micro batteries to EV batteries, is developed. The Cylindrical battery design V1.0 is comprised of three types of cylindrical batteries, Micro battery (Primary), Micro battery (Secondary) and 18650/21700/xxxxx cylindrical battery. The software was developed in MATLAB. The software has the capability to output the cell design with the capacity ranges from several mAh to several million Ah. The software utilizes machine learning and includes a graphical user interface to enable rapid prototyping to accelerate energy storage research, development, and manufacturing.

25 ENERGY STORAGE↗

Predicting High Energy Arcing Fault Zones of Influence for Aluminum Using an Arc Flash Modeling Approach: Evaluation of a model bias, uncertainty, parameter sensitivity and zone of influence estimation

This report documents the development of an arc flash hazard model to calculate the incident energy and zone of influence from high energy arcing faults involving aluminum. The NRC has identified the potential for (HEAFs) involving aluminum to increase the damage zone beyond what is currently postulated in fire probabilistic risk assessment (PRA) methodologies. To estimate the hazard from HEAFs involving aluminum an arc flash model was developed. Differences between the initial model and nuclear power plant (NPP) fire PRA scenarios were identified. Modification of the initial model established from existing literature and test data was used to minimize these differences. The developed model was evaluated against NRC datasets to understand the model prediction and relative uncertainties. Finally, a range of fire PRA zone of influences (ZOI) were developed based on the developed model, target fragility estimates and update HEAF PRA methodology. The results were developed to support an NRC LIC-504 evaluation in tandem with other modeling efforts. The report documents the effort and provides a reference for any future advancements in arc flash modeling.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Accelerating computational fluid dynamics simulation of post-combustion carbon capture modeling with MeshGraphNets

Packed columns are commonly used in post-combustion processes to capture CO 2 emissions by providing enhanced contact area between a CO 2 -laden gas and CO 2 -absorbing solvent. To study and optimize solvent-based post-combustion carbon capture systems (CCSs), computational fluid dynamics (CFD) can be used to model the liquid–gas countercurrent flow hydrodynamics in these columns and derive key determinants of CO 2 -capture efficiency. However, the large design space of these systems hinders the application of CFD for design optimization due to its high computational cost. In contrast, data-driven modeling approaches can produce fast surrogates to study large-scale physics problems. We build our surrogates using MeshGraphNets (MGN), a graph neural network framework that efficiently learns and produces mesh-based simulations. We apply MGN to a random packed column modeled with over 160K graph nodes and a design space consisting of three key input parameters: solvent surface tension, inlet velocity, and contact angle. Our models can adapt to a wide range of these parameters and accurately predict the complex interactions within the system at rates over 1700 times faster than CFD, affirming its practicality in downstream design optimization tasks. This underscores the robustness and versatility of MGN in modeling complex fluid dynamics for large-scale CCS analyses.

97 MATHEMATICS AND COMPUTING↗

A Universal Power-law Prescription for Variability from Synthetic Images of Black Hole Accretion Flows

We present a framework for characterizing the spatiotemporal power spectrum of the variability expected from the horizon-scale emission structure around supermassive black holes, and we apply this framework to a library of general relativistic magnetohydrodynamic (GRMHD) simulations and associated general relativistic ray-traced images relevant for Event Horizon Telescope (EHT) observations of Sgr A*. We find that the variability power spectrum is generically a red-noise process in both the temporal and spatial dimensions, with the peak in power occurring on the longest timescales and largest spatial scales. When both the time-averaged source structure and the spatially integrated light-curve variability are removed, the residual power spectrum exhibits a universal broken power-law behavior. On small spatial frequencies, the residual power spectrum rises as the square of the spatial frequency and is proportional to the variance in the centroid of emission. Beyond some peak in variability power, the residual power spectrum falls as that of the time-averaged source structure, which is similar across simulations; this behavior can be naturally explained if the variability arises from a multiplicative random field that has a steeper high-frequency power-law index than that of the time-averaged source structure. We briefly explore the ability of power spectral variability studies to constrain physical parameters relevant for the GRMHD simulations, which can be scaled to provide predictions for black holes in a range of systems in the optically thin regime. We present specific expectations for the behavior of the M87* and Sgr A* accretion flows as observed by the EHT.

79 ASTRONOMY AND ASTROPHYSICS↗

A Hybrid Heavy-Duty Diesel Power System for Off-Road Applications - Concept Definition

A multi-year Power System R&D project was initiated with the objective of developing an off-road hybrid heavy-duty concept diesel engine with front end accessory drive-integrated energy storage. This off-road hybrid engine system is expected to deliver 15-20% reduction in fuel consumption over current Tier 4 Final-based diesel engines and consists of a downsized heavy-duty diesel engine containing advanced combustion technologies, capable of elevated peak cylinder pressures and thermal efficiencies, exhaust waste heat recovery via SuperTurbo™ turbocompounding, and hybrid energy recovery through both mechanical (high speed flywheel) and electrical systems. The first year of this project focused on the definition of the hybrid elements using extensive dynamic system simulation over transient work cycles, with hybrid supervisory controls development focusing on energy recovery and transient load assist, in Caterpillar’s DYNASTY™ software environment. Three key off-road applications were the focus of the hybrid concept definition with an aim of understanding the system’s modular capability for the diverse off-road heavy-duty market. Core engine performance 1D and 3D simulations isolated the efficiency contributions from the downsized engine, turbocompounding, and in-cylinder thermal barrier coatings. A fuel consumption improvement range of 14 to 24% was predicted, resulting in successful project progression to the design and experimental validation phase. Furthermore, an overview of the experimental engine and hybrid system status concludes the discussion along with the multi-year project’s next steps.

33 ADVANCED PROPULSION SYSTEMS↗

A Combined Machine Learning and High-Energy X-ray Diffraction Approach to Understanding Liquid and Amorphous Metal Oxides

Determining the structure-property relations of liquid and amorphous metal oxides is challenging, due to their variable short-range order and polyhedral connectivity. To predict chemically realistic structures, we have developed a Machine Learned, Gaussian Approximation Potential (GAP) for HfO 2 , with a focus on enhanced sampling of the training database and accurate density functional theory calculations. By using training datasets for the GAP model at the level of Density Functional Theory-Strongly Constrained and Appropriately Normed (DFT-SCAN) level of theory, our results show that the topology of both the low viscosity liquid and the amorphous form are dominated by edge-shared chains and small corner-shared rings of polyhedra. This topology is shown to be consistent with the structure of other liquid and amorphous transition metal oxides of variable ion size, such as TiO 2 and ZrO 2 . Current limitations of the ML-GAP modeling method for obtaining glass structures and future perspectives are also discussed.

36 MATERIALS SCIENCE↗

Application of an Empirical Density Law via Python for Aqueous Plutonium Chloride Systems in MCNP6

Current aqueous plutonium processing models for criticality safety often contain significant bias due to material modeling assumptions. These solutions include plutonium chloride solutions, which are modeled as fictitious plutonium metal-water mixtures because little is known about the actual density of the solution. Furthermore, there is no current predictive capability for modeling plutonium metal-water mixtures that is approved for use at Los Alamos National Laboratory (LANL). Recent density measurements for aqueous plutonium chloride systems (PuCl 3 -HCl-H 2 O) now allow for the development and application of a more realistic density law. This work develops a Python-based density law for this ternary solution using an empirical method. This code can be used in conjunction with an MCNP6 input to determine the density and composition of a solution based upon user inputs of plutonium concentration, hydrochloric acid concentration, and temperature. The tool allows users to input plutonium and acid content of a solution in terms of molality, molarity, or concentration, and predicts density within the current data range within 1.4% of experimental data. Current preliminary MCNP6 calculations utilizing this tool have demonstrated a minimum decrease in system reactivity of 5% in comparison to the current modeling conventions. Thus, this tool enables more accurate criticality safety operational limits by better crediting chorine content while still maintaining necessary conservatism.

Density Law↗

Decayheatml

This code is designed to predict and analyze the decay heat generated in molten salt reactors (MSRs) using a hybrid approach that combines machine learning and segmented polynomial fitting. The accurate prediction of decay heat is essential for reactor safety and the optimization of spent fuel storage. The code operates through several key components: 1) Data Architecture: It incorporates a modular data architecture that handles various MSR-specific operational parameters such as power density, humidity content, and air ingress. These parameters are sampled using Sobol sequences to ensure comprehensive coverage of operational uncertainties. 2) Machine Learning Framework: The code employs a diverse set of machine learning models, including polynomial regression, decision trees, random forests, gradient boosting, support vector regression, k-nearest neighbors, multi-layer perceptrons, and symbolic regression. These models are trained to predict decay heat over a wide temporal range, from immediate shutdown up to 10,000 years. 3) Region-Optimized Training: The temporal domain is divided into multiple regions, each modeled separately to capture distinct decay heat characteristics across different time scales. This approach significantly improves the accuracy and interpretability of predictions. 4) Segmented Polynomial Interpretation (SPI): The SPI method translates machine learning predictions into piecewise polynomial equations. These equations are physically interpretable and can be directly integrated into existing engineering workflows and safety analyses. 5) Front-End Interfaces: The code includes both a Jupyter notebook interface for research development and a Streamlit web application for operational deployment. These interfaces allow users to interactively explore decay heat predictions, adjust operational parameters, and visualize results in real-time. 6) Applications: The framework supports various applications, including safety system validation and spent fuel container optimization. It enables real-time evaluation of worst-case decay heat scenarios, informing the design of passive safety systems and optimizing container designs for long-term storage. Overall, this code provides a robust, accurate, and user-friendly tool for predicting decay heat in MSRs, enhancing reactor safety, and optimizing spent fuel management.

Retamales, Mauricio Eduardo Tano [Idaho National L↗

Micrometer: Micromechanics transformer for predicting full field mechanical responses of heterogeneous materials

Predicting mechanical responses of heterogeneous materials across scales remains a significant challenge. Traditional computational methods often struggle with complex and multiscale nature of these materials, limiting their effectiveness in real-world applications. Here, in this paper, we introduce Micrometer, a vision transformer based deep learning model designed to predict full field mechanical responses of heterogeneous materials, bridging the gap between computer vision and solid mechanics problems. We show that Micrometer, trained on a large-scale high-resolution dataset of 2D fiber-reinforced composites, can achieve state-of-the-art performance in predicting microscale strain fields across a wide range of material properties and loading conditions. Our model demonstrates accuracy and computational efficiency in applications such as computational homogenization and multiscale modeling, reducing computational time by up to two orders of magnitude compared to conventional numerical solvers while maintaining less than 1 % errors in predicting macroscale stress fields. Furthermore, we showcase Micrometer’s adaptability through transfer learning experiments on new materials with limited data, highlighting its potential to tackle diverse scenarios in computational solid mechanics. These results represent a significant step towards AI-driven innovation in materials science, addressing the limitations of traditional numerical methods and paving the way for more efficient simulations of heterogeneous materials across various industrial applications.

Composite materials↗

Multitude Characterization and Prediction of DOE Advanced Biofuels Properties

Advanced multitude of experiments ranging from the liquid fuel to combustion are conducted on Co-OPTIMA fuels to aid the characterization of the fuels. The series of targeted experiments characterized Co-OPTIMA fuel spray atomization, flame topology, flame speed, autoignition, volatility, viscosity, soot/coking, and compatibility. The fuels are selected and prioritized based on input from national lab members. The research characterized and predicted biomass-based, low greenhouse gas fuels and blends combustion, autoignition, and physical properties of mixtures of identified compounds at engine-relevant conditions, in particular those properties that blend non-linearly. The main challenge of Co-OPTIMA is the evaluation of a variety of biofuels and blends in all the reaction conditions that might be encountered in new high-efficiency engines. Despite the improved high-throughput experimental techniques, it seems unlikely that all of the performance metrics could be measured for all relevant petroleum derived, bio-derived molecules and mixtures, and reaction conditions. A series of targeted experiments ranging from the liquid fuel to the combustion process is required, and to extract the maximum information from each experiment. These targeted experiments evaluated how a specific fuel will perform in an engine. The series of targeted experiments are as follows: (a) Spray Atomization, Vaporization and Droplet Formation (b) Combustion Flame and Local Fuel/Air Image-Based Measurements (c) Laminar Flame Speed Measurements (d) Autoignition and Soot Measurements (e) Synchrotron Coupled Fundamental Autoignition Experiments (f) Fuel Coking and Hot Surface Deposit g) Fuel Volatility Measurements (h) Viscosity Measurements (i) Seal Flexible Fuel Compatibility These experimental processes provide an essential pathway for the prediction of fuel behaviors in engines and systematic process for fuel down select.

09 BIOMASS FUELS↗

Accelerating Informed Local Climate Action by Improving Access to DOE Predictive Modeling and Science-Based Tools

The threat of climate change is becoming increasingly urgent, and local governments across the United States are working to develop and implement climate action plans to reduce greenhouse gas emissions and mitigate the impacts of climate change. At the same time, the U.S. Department of Energy's Office of Science, Earth and Environmental Systems Sciences Division (EESSD) has been developing a range of tools and technologies to help predict the potential multi-faceted impacts from climate change given varying environmental and socioeconomic scenarios. However, local governments have been slow to engage with these scientific assets due to a lack of understanding and awareness. The DOE researchers interviewed highlighted three primary types of resources that could benefit local governments and climate practitioners: climate models, modeled datasets, and scientific journal articles. Although many of these resources are open source and available to the public, there were several challenges identified in translating these resources and their research insights directly to decision-makers.

54 ENVIRONMENTAL SCIENCES↗

Dynamic Modeling and Simulation of a Subcritical Coal-Fired Power Plant under Load-Following Conditions

Dynamic models for power plants that capture realistic general process trends and effects of manipulated variables are needed to improve load-following, while minimizing carbon footprint. In this work, a dynamic modeling approach and simulation results for subcritical coal-fired power plant components are presented. These encompass simulation of the dynamics in the fireside, including the effects of fuel, air combustion, and the dynamics of the entire waterside and power generation sections. This model development enables the simulation and analysis of the important short and long-time scale dynamics of components such as heaters, evaporative loop, and power generation units. Furthermore, additional variables in the power generation section are introduced to improve model accuracy, extending the prediction capability of subcritical power plant models and opening new opportunities for research in operator training, optimization, and advanced model-based controller design that are based on these models. The change in process gain for different ramp rates associated with disturbance signals that affect process variables is also explored and a correlation developed. This provides opportunities to study disturbance rejection control implementation and adaptation for scenarios with such variations in ramp rates. The prediction capabilities of selected components are compared to data available in literature, with the obtained root mean squared error ranges that reflect the model performance and quality of predictions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamics of dark matter misalignment through the Higgs portal

A light singlet scalar field feebly coupled through the super-renormalizable Higgs portal provides a minimal and well-motivated realization of ultra-light bosonic dark matter. We study the cosmological production of dark matter in this model by elucidating the dynamics of two sources of scalar field misalignment generated during the radiation era. For large scalar masses (above $\mathcal{O}$ (10 −3 eV)), dark matter is produced through thermal misalignment, by which the scalar field is driven towards large field values as a result of the finite-temperature effective potential. The dominance of thermal misalignment in this mass range leads to a sharp relic abundance prediction which is, to a significant extent, insensitive to the initial conditions of the scalar field. On the other hand, for low mass scalars (below $\mathcal{O}$ (10 −5 eV)), dark matter is produced via VEV misalignment, which is caused by the induced scalar field vacuum expectation value triggered by the electroweak phase transition. We show that the relic abundance in this low mass range is sensitive to the scalar field initial conditions. In the intermediate mass range, the relic abundance is a consequence of a competition between thermal misalignment and VEV misalignment, which can potentially lead to novel forced resonance effects which cause a recurring enhancement and suppression in the late time oscillation amplitude as a function of the scalar mass. We compare our relic abundance predictions with constraints and projections from equivalence principle and inverse square law tests, stellar cooling, resonant molecular absorption, and observations of extra-galactic background light and diffuse X-ray backgrounds. New experimental ideas are needed to probe most of the cosmologically motivated regions of parameter space.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Automated Operational Forecasting of Monsoon Low Pressure Systems

Monsoon low pressure systems (LPSs) are the dominant rain-bearing weather system of South Asia, often producing extreme precipitation and hydrological disasters in a region inhabited by nearly two billion people. Despite the importance of these storms, no operational system has automatically identified and tracked LPS in real time in numerical weather prediction model output; many commonly used vortex-tracking algorithms are ill suited for monsoon LPS because of the weak winds and cold cores of these systems. Here, we describe a new system that uses optimized algorithms to identify monsoon LPS in short- to medium-range forecasts from the U.S. Global Ensemble Forecast System (GEFS) and a version of the deterministic Global Forecast System (GFS) adapted and used operationally by the Indian Institute of Tropical Meteorology (IITM). We also assess the historical performance of these models in forecasting South Asian monsoon LPS, comparing this with the performance of the Integrated Forecasting System of the ECMWF. We assess the accuracy of model predictions of LPS genesis, position, intensity, and precipitation rates for forecast lead times of 1–5 days, yielding quantitative information on model biases to guide operational forecasters and disaster managers. The system we introduce here could be extended to other low-latitude regions affected by dynamically weak, heavily precipitating atmospheric vortices that are often not included in tropical cyclone inventories.

54 ENVIRONMENTAL SCIENCES↗

Advanced surrogate model for electron-scale turbulence in tokamak pedestals

We derive an advanced surrogate model for predicting turbulent transport at the edge of tokamaks driven by electron temperature gradient (ETG) modes. Our derivation is based on a recently developed sensitivity-driven sparse grid interpolation approach for uncertainty quantification and sensitivity analysis at scale, which informs the set of parameters that define the surrogate model as a scaling law. Our model reveals that ETG-driven electron heat flux is influenced by the safety factor q, electron beta β e and normalized electron Debye length λ D , in addition to well-established parameters such as the electron temperature and density gradients. To assess the trustworthiness of our model's predictions beyond training, we compute prediction intervals using bootstrapping. The surrogate model's predictive power is tested across a wide range of parameter values, including within-distribution testing parameters (to verify our model) as well as out-of-bounds and out-of-distribution testing (to validate the proposed model). Overall, validation efforts show that our model competes well with, or can even outperform, existing scaling laws in predicting ETG-driven transport.

fusion plasma↗

Group additivity values for entropy and heat capacities of C 2 –C 8 alkanes, alkyl hydroperoxides, and their radicals

Group additivity values for the thermodynamic properties of oxygenated radicals are poorly determined due to the absence of high quality reference data. Here, a set of 58 group additive values (GAV) for the standard entropy and heat capacity of relevance to alkane oxidation is derived from fits to an extensive and accurate database of standard entropies (298.15 K) and heat capacities (300–3000 K) recently calculated with the “STAR-1D” formalism. The 192 species in this database represent the alkanes (RH), alkyl radicals ($\dot{R}$), alkyl hydroperoxides (RO 2 H), alkyl-peroxy (R$\dot{O}$ 2 ) and hydroperoxy-alkyl ($\dot{Q}$OOH) radicals for all of the isomers of C 2 –C 5 alkane fuels and a select number of isomers of C 6 –C 9 species. The STAR-1D thermochemical data for this set of species was previously obtained from a coupling of scaled B2PLYPD3/cc-pVTZ vibrational analyses with scaled ωB97X-D/cc-pVTZ one-dimensional hindered rotor corrections. The 2σ uncertainties in the GAV results relative to the STAR-1D data set are 2.4 cal K –1 mol –1 for the entropies and at most 2.0 cal K –1 mol –1 for the heat capacities in the temperature range 500–800 K. The 2σ fitting uncertainties in the heat capacity gradually reduce at higher temperatures reaching a value of 0.8 cal K –1 mol –1 at 2000 K decreasing to only 0.5 cal K –1 mol –1 at 3000 K. The high degree of accuracy for the GAV representations is obtained through the introduction of various new group terms, together with the re-optimization of existing group terms. Among the full set of 58 GAV terms, 25 include non-next-nearest neighbor interactions (NNI) and β-corrections. The updated GAVs can be applied in the prediction of entropies and heat capacities for a wide range of hydrocarbons and hydroperoxide species and their radicals, which is important to the accurate prediction of fuel reactivity at low-temperatures in the range 600–1000 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Riparian vegetation shade restoration and loss effects on recent and future stream temperatures

River temperatures are expected to increase this century harming species requiring cold-water habitat unless restoration activities protect or improve habitat availability. Local shading by riparian vegetation can cool water temperatures, but uncertainty exists over the scaling of this local effect to larger spatial extents. Here, we evaluate this issue using a regional spatial stream network temperature model with covariates representing shade effects to predict mean August stream temperatures across 78,195 km of tributaries flowing into the Columbia River in the northwestern United States. We evaluate nine scenarios predicting stream temperatures for three riparian shade conditions (current, restored, and no riparian vegetation) within three different climate periods (2000s, 2040s, and 2080s). Results suggest riparian shade restoration (2000s climate) could decrease mean August stream temperatures by 0.62°C across the study network. Under the same restored shade conditions, temperature predictions for tributaries at their confluence with the Columbia River range from 0.02 to 2.08°C cooler than under current shade conditions. The climate warming effect predicted for the 2040s and 2080s, however, is greater than the cooling effect from restoring riparian shade. Streams less than 10-m bankfull width cooled more frequently with riparian shade restoration. In Oregon, the proportion of fish habitat for salmon and trout rearing and migration that meet temperature numeric water quality criteria could be increased by 20% under restored shade conditions although net habitat declines may still occur in the future. We conclude riparian vegetation restoration could partially mitigate future warming and help maintain cold-water habitats that function as thermal refuges if implemented strategically.

54 ENVIRONMENTAL SCIENCES↗