Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “retrievals”

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

Developing and Implementing New Waste Chemistry Controls for Hanford's Double-Shell Tanks - 20041

Waste chemistry controls for the Hanford double-shell tanks (DSTs) were established in the 1980's in response to tank failures caused by stress corrosion cracking (SCC) at the Savannah River Site (SRS). The controls were established based on a combination of SRS chemistry limits and corrosion testing results. Conditions in the Hanford DSTs have changed over time, with the temperature in the majority of the DSTs dropping below 50 deg. C. Testing indicated that the original chemistry control specification and SCC testing performed in 2010 did not protect against pitting corrosion mechanisms - what is believed to have caused the failure of tank 241-AY-102. New waste chemistry controls are necessary to protect the DSTs from both SCC and pitting corrosion. Savannah River National Laboratory (SRNL), with input from the Tank Integrity Expert Panel (TIEP) Corrosion Subgroup (CSG), performed a statistically based investigation of nitrate and halide ion induced pitting corrosion. The investigation was intended to develop a comprehensive waste chemistry envelope that minimizes the risk of both SCC and pitting caused by halide and nitrate ions. The waste chemistry envelope needed to be robust enough to address future waste retrieval and process changes which could significantly change the waste composition of the DSTs. The experimental design allowed the statistically significant ions and their effect on pitting corrosion to be determined. SRNL developed the 'pitting factor' which is an empirical relationship between the statistically significant inhibitive species (hydroxide, and nitrite), and pitting inducing species (nitrate, chloride, and fluoride). The coefficients are weighted factors that show the degree of influence the species have on the propensity for pitting. New waste chemistry control limits were recommended by SRNL, and concurred with by the TIEP CSG, that incorporated the pitting factor along with other limits. The controls developed to reduce the risk of pitting corrosion were evaluated to determine if they also protected against SCC. Hundreds of test results were reviewed and compared against the proposed limits. With the exception of one test result, instances where cracking was observed corresponded to chemistry that was outside of the proposed limits - indicating that the chemistry controls developed for pitting also adequately addressed SCC concerns. Existing tank composition estimates were evaluated against the new chemistry controls - three current tank chemistries are suspect, and may not meet the new chemistry control requirements. These included the interstitial liquid (liquid associated with solids) of DSTs 241-AN-102, 241-AN-106, 241-AN-107, and 241-AY-101. Measures are ongoing to evaluate the tanks with the suspect chemistry. These include varying combinations of: performing additional evaluations (history of use, modeling); DST specific laboratory testing; and, pursuing core samples of the DSTs. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Reducing AI RAG Hallucination by Optimizing Routing Techniques

Large Language Models (LLMs), such as ChatGPT, tend to “hallucinate”, meaning they confidently generate false information. Retrieval Augmented Generation (RAG) attempts to diminish hallucination by providing context to the LLM from data stores (indexes) containing relevant information. The LLM uses this context to formulate its response. RAG systems can still suffer from hallucination because of bad embeddings or ineffective routing. For example, a router will often return context from an irrelevant index, resulting in a hallucinated answer. In this study, we aim to minimize the frequency of routing hallucinations by optimizing Index Summary Routing.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A Processing and Analytics System for Microscopy Data Workflows: The Pycroscopy Ecosystem of Packages

Major advancements in fields as diverse as biology and quantum computing have relied on a multitude of microscopy techniques. Despite the considerable proliferation of these instruments, significant bottlenecks remain in terms of processing, analysis, storage, and retrieval of the acquired datasets. Aside from lack of file standards, individual domain-specific analysis packages are often disjoint from the underlying datasets, and thus keeping track of analysis and processing steps remains tedious for the end-user, hampering reproducibility. Here, in this study, the pycroscopy ecosystem of packages is introduced, an open-source python-based ecosystem underpinned by a common data model. The data model, termed the N-dimensional spectral imaging data format, is realized in pycroscopy's sidpy package. This package is built on top of dask arrays, thus leveraging dask array attributes, but expanding them to accelerate microscopy relevant analysis and visualization. Several examples of the use of the pycroscopy ecosystem to create workflows for data ingestion and analysis of scanning transmission electron microscopy (STEM) and scanning probe microscopy data are shown. Adoption of such standardized routines will be critical to usher in the next generation of autonomous instruments where processing, computation, and meta-data storage will be critical to overall experimental operations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Leaf traits and canopy structure together explain canopy functional diversity: an airborne remote sensing approach

Plant functional diversity is strongly connected to photosynthetic carbon assimilation in terrestrial ecosystems. However, many of the plant functional traits that regulate photosynthetic capacity, including foliar nitrogen concentration and leaf mass per area, vary significantly between and within plant functional types and vertically through forest canopies, resulting in considerable landscape-scale heterogeneity in three dimensions. Hyperspectral imagery has been used extensively to quantify functional traits across a range of ecosystems but is generally limited to providing information for top-of-canopy leaves only. On the other hand, lidar data can be used to retrieve the vertical structure of forest canopies. Because these data are rarely collected at the same time, there are unanswered questions about the effect of forest structure on the 3-D spatial patterns of functional traits across ecosystems. In the United States, the National Ecological Observatory Network's Airborne Observation Platform (NEON AOP) provides an opportunity to address this structure-function relationship by collecting lidar and hyperspectral data together across a variety of ecoregions. With a fusion of hyperspectral and lidar data from the NEON AOP and field-collected foliar trait data, we assessed the impacts of forest structure on spatial patterns of N. In addition, we examine the influence of abiotic gradients and management regimes on top of canopy %N and total canopy N (i.e. the total amount of N (g/m 2 ) within a forest canopy) at a NEON site consisting of a mosaic of open longleaf pine and dense broadleaf deciduous forests. Here, our resulting maps suggest that in contrast with top-of-canopy values, total canopy N variation is dampened across this landscape resulting in relatively homogeneous spatial patterns. At the same time, we found that leaf functional diversity and canopy structural diversity showed distinct dendritic patterns related to the spatial distribution of plant functional types.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

On the numerical accuracy in finite-volume methods to accurately capture turbulence in compressible flows

The goal of the present article is to understand the impact of numerical schemes for the reconstruction of data at cell faces in finite-volume methods, and to assess their interaction with the quadrature rule used to compute the average over the cell volume. Here, third-, fifth- and seventh-order WENO-Z schemes are investigated. On a problem with a smooth solution, the theoretical order of convergence rate for each method is retrieved, and changing the order of the reconstruction at cell faces does not impact the results, whereas for a shock-driven problem all the methods collapse to first-order. Here, study of the decay of compressible homogeneous isotropic turbulence reveals that using a high-order quadrature rule to compute the average over a finite-volume cell does not improve the spectral accuracy and that all methods present a second-order convergence rate. However the choice of the numerical method to reconstruct data at cell faces is found to be critical to correctly capture turbulent spectra. In the context of simulations with finite-volume methods of practical flows encountered in engineering applications, it becomes apparent that an efficient strategy is to perform the average integration with a low-order quadrature rule on a fine mesh resolution, whereas high-order schemes should be used to reconstruct data at cell faces.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A Study on Recycled Polymers Recovered from Multilayer Plastic Packaging Films by Solvent–Targeted Recovery and Precipitation (STRAP)

Multilayer plastic film use increased in multiple packaging applications due to its versatility and overall increased performance over monolayer structures. However, the performance gains from multiple layers also make recycling difficult because they contain multiple polymers that can be immiscible and burdensome to traditional mechanical recycling operations. A possible solution is the solvent-targeted recovery and precipitation (STRAP) process, but the effect on the retrieved polymer is still unknown. Here, the STRAP process is applied to two different multilayer films and samples of recovered polymers are evaluated for physical, molecular, and thermal properties. Changes in the molecular weight are not significant, but differences in thermal properties are reported along with the coprecipitation of different polymers. Solvent retention in the polymer matrix from STRAP reduced the glass transition temperature of samples, but enhanced drying recovered it. Heavy metals, such as Cd, Cr, and Pb are not detected, indicating regulatory compliance for different applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Drivers for the cracking of multilayer polyamide‐based backsheets in field photovoltaic modules: In‐depth degradation mapping analysis

Abstract There is a lack of understanding on the root cause of cracking of photovoltaic (PV) backsheets due to the challenge of multilayer characterization and the complicated failure modes at the submodule level. In this work, the in‐depth degradation mapping of field‐exposed polyamide‐based (PA‐based) PV module backsheets was studied, with the major focus on the identification of underlying drivers for the through cracking (in between the solar cells). PV modules were retrieved from five different locations, comprising a variety of climates, including humid subtropical, hot‐summer Mediterranean, tropical savanna climate, and hot arid. A suite of microscale cross‐sectional characterizations, including chemical changes, fluorescence intensity, and modulus as a function of distance from the air surface of the backsheet, was performed. Results showed more advanced signs of degradation of the inner layer than the outer layer in the cracking region. Increases in the modulus were identified as the major indicator for the cracking. Moreover, a rudimentary test by immersion in acetic acid, which forms during photodegradation of the ethylene‐vinyl acetate copolymer (EVA) encapsulant, showed the first‐time direct evidence that acetic acid can largely accelerate the chemical degradation and facilitate the cracking of PA inner layer. This study suggests that the field cracking of PA‐based backsheet can be attributed to the combined effects of chemical degradation and physical reorganization (chemi‐crystallization) under cyclic thermomechanical stresses. Therefore, it is important to consider effects of local microclimates and interlayer infection to understand the heterogeneous nature of backsheet failures.

Lyu, Yadong↗
DOE OSTI2022

Predicting Search Task Difficulty through a Discrete‐Time Action Log Representation on Spectrum Kernel

ABSTRACT Predicting perceived difficulty on a web search task is an open problem in the interactive information retrieval field. A common approach to tackle it, is through features obtained from full search sessions, which are then used to train classification models. In this poster we attempt to predict perceived task difficulty at different stages of the search process. To do so, we use the spectrum kernel for support vector machine (SVM) classification. Our preliminary results suggest that by using behavioral data from the first query segment, it is possible to provide timely classifications of whether a search task is perceived as hard or easy.

Gacitúa, Daniel↗
DOE OSTI2023

Hollow Ptychography: Toward Simultaneous 4D Scanning Transmission Electron Microscopy and Electron Energy Loss Spectroscopy

With the recent development of high-acquisition-speed pixelated detectors, 4D scanning transmission electron microscopy (4D-STEM) is becoming routinely available in high-resolution electron microscopy. 4D-STEM acts as a “universal” method that provides local information on materials that is challenging to extract from bulk techniques. It extends conventional STEM imaging to include super-resolution techniques and to provide quantitative phase-based information, such as differential phase contrast, ptychography, or Bloch wave phase retrieval. However, an important missing factor is the chemical and bonding information provided by electron energy loss spectroscopy (EELS). 4D-STEM and EELS cannot currently be acquired simultaneously due to the overlapping geometry of the detectors. Here, the feasibility of modifying the detector geometry to overcome this challenge for bulk specimens is demonstrated, and the use of a partial or defective detector for ptycholgaphic structural imaging is explored. Here, results show that structural information beyond the diffraction-limit and chemical information from the material can be extracted together, resulting in simultaneous multi-modal measurements, adding the additional dimensions of spectral information to 4D datasets.

4D-scanning transmitting electron microscopy (STEM↗
DOE OSTI2024

The Central Role of Oxo Clusters in Zirconium‐Based Esterification Catalysis

Oxo clusters are a unique link between oxide nanocrystals and Metal‐Organic Frameworks (MOFs), representing the limit of downscaling each of the respective crystals. Herein, the superior catalytic activity of clusters, compared to zirconium MOF UiO‐66 and nanocrystals is shown. Focus is on esterification reactions given their general importance in consumer products and the challenge of converting large substrates. Oxo clusters have a higher surface‐to‐volume ratio than nanocrystals, rendering them more active. For large substrates, for example, oleic acid, MOF UiO‐66 has negligible catalytic activity while clusters provide almost quantitative conversion, a fact we ascribe to limited diffusion of large substrates through the MOF pores. Clusters do not suffer from limited mass transfer and we also obtain high conversion in solvent‐free reactions with sterically hindered alcohols (hexanol, 2‐ethyl hexanol, benzyl alcohol, and neopentyl alcohol). The cluster catalyst can be recovered and shows identical activity when reused. The structural integrity of the cluster is confirmed using X‐ray total scattering and pair distribution function analysis. Moreover, when homogeneous zirconium alkoxides are used as catalysts, the same oxo cluster is retrieved, showing that oxo clusters are the active catalytic species, even in previously assumed homogeneously catalyzed reactions.

Pulparayil Mathew, Jikson↗
DOE OSTI2023

Effective climate sensitivity distributions from a 1D model of global ocean and land temperature trends, 1970–2021

Abstract Current theoretically based Earth system models (ESMs) produce Effective Climate Sensitivities (EffCS) that range over a factor of three, with 80% of those models producing stronger global warming trends for 1970–2021 than do observations. To make a more observationally based estimate of EffCS, a 1D time-dependent forcing-feedback model of temperature departures from energy equilibrium is used to match measured ranges of global-average surface and sub-surface land and ocean temperature trends during 1970–2021. In response to two different radiative forcing scenarios, a full range of three model free parameters are evaluated to produce fits to a range of observed surface temperature trends (± 2σ) from four different land datasets and three ocean datasets, as well as deep-ocean temperature trends and borehole-based trend retrievals over land. Land-derived EffCS are larger than over the ocean, and EffCS is lower using the newer Shared Socioeconomic Pathways (SSP245, 1.86 °C global EffCS, ± 34% range 1.48–2.15 °C) than the older Representative Concentration Pathway forcing (RCP6, 2.49 °C global average EffCS, ± 34% range 2.04–2.87 °C). The strongest dependence of the EffCS results is on the assumed radiative forcing dataset, underscoring the role of radiative forcing uncertainty in determining the sensitivity of the climate system to increasing greenhouse gas concentrations from observations alone. The results are consistent with previous observation-based studies that concluded EffCS during the observational period is on the low end of the range produced by current ESMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Observations of the Marine Atmospheric Boundary Layer’s Response to a Solar Eclipse

The atmospheric response to the solar eclipse of 8 April 2024 in North America is investigated with a specific focus on the marine atmospheric boundary layer (MABL). We leverage measurements collected during the Third Wind Forecast Improvement Project (WFIP3), including Doppler lidars, sonic anemometers, and thermodynamic profiler data to investigate the atmospheric response across sites that experienced partial eclipse conditions with nearly 90% obscuration. Using these measurements, we examine eclipse-induced changes in key meteorological parameters, such as temperature, wind speed, and turbulent fluxes. Most previous eclipse studies have been conducted over land, whereas this study provides new observations for both coastal and marine environments, offering additional insight into eclipse-driven variability in the MABL. The findings confirm a notable decrease in downwelling shortwave radiation during the eclipse, which results in rapid cooling of surface air. The temperature reduction ranges from $1.2^\circ \text {C}$ to $1.4^\circ \text {C}$ in coastal regions and from $0.3^\circ \text {C}$ to $0.5^\circ \text {C}$ over the ocean. This analysis suggests that the MABL’s higher thermal inertia compared to coastal regions moderates the temperature decrease during the eclipse. Wind speed exhibits a more complex behavior, as it is influenced by both the MABL and preexisting synoptic conditions. Although a reduction in wind speed is observable up to approximately 140 m above ground level (AGL) at more inland sites, at other locations closer to the coast, this reduction is constrained to the lowest 100 m AGL. Turbulence parameters retrieved from sonic anemometers, such as turbulence kinetic energy, turbulent heat flux, and friction velocity, decrease during the eclipse at coastal sites, accompanied by a brief transition of atmospheric stability from unstable to neutral or weakly stable conditions. For the open-ocean sites, the variability in turbulence statistics and atmospheric stability is minimal during the occurrence of the eclipse.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

A phase transition for finding needles in nonlinear haystacks with LASSO artificial neural networks

To fit sparse linear associations, a LASSO sparsity inducing penalty with a single hyperparameter provably allows to recover the important features (needles) with high probability in certain regimes even if the sample size is smaller than the dimension of the input vector (haystack). More recently learners known as artificial neural networks (ANN) have shown great successes in many machine learning tasks, in particular fitting nonlinear associations. Small learning rate, stochastic gradient descent algorithm and large training set help to cope with the explosion in the number of parameters present in deep neural networks. Yet few ANN learners have been developed and studied to find needles in nonlinear haystacks. Driven by a single hyperparameter, our ANN learner, like for sparse linear associations, exhibits a phase transition in the probability of retrieving the needles, which we do not observe with other ANN learners. To select our penalty parameter, we generalize the universal threshold of Donoho and Johnstone (Biometrika 81(3):425–455, 1994) which is a better rule than the conservative (too many false detections) and expensive cross-validation. In the spirit of simulated annealing, we propose a warm-start sparsity inducing algorithm to solve the high-dimensional, non-convex and non-differentiable optimization problem. We perform simulated and real data Monte Carlo experiments to quantify the effectiveness of our approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Mapping the kinetic evolution of metastable grain boundaries under non-equilibrium processing

The kinetic evolution of a multiplicity of metastable grain boundaries (GBs) under fast driving conditions are studied by atomistic modeling. Assisted with an enhanced statistical analysis, the energetic evolution of GBs over a broad metastability-temperature space is mapped out, wherein two distinct regimes—an ageing regime and a rejuvenating regime—are retrieved with high fidelity. By comparing the results under various conditions (e.g. random perturbations, isothermal annealing, and fast heating/cooling), it is shown that such ageing/rejuvenating mechanism map is universal, irrespective of the actual stimuli used to elicit the metastable GBs. Here, the ageing/rejuvenating phenomena are demonstrated to stem from the energy imbalance of uphill climbing and downhill dropping during sequential transitions in the system's potential energy landscape. Without the necessity of introducing free parameters, such model can reconcile experimentally measured hardness variation of nanocrystalline metals subjected to femto-second laser irradiation, and it therefore provides a novel perspective on achieving a plurality of interfacial states and facilitating previously inaccessible property regimes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

AlphaBuilding ResCommunity: A multi-agent virtual testbed for community-level load coordination

Training and validating algorithms in a simulation testbed can accelerate research and applications of optimal control of residential loads to improve energy flexibility and grid resilience. We developed an open-source simulation environment, AlphaBuilding ResCommunity, that can be used to train and validate algorithms to control a single thermostatically controlled load (TCL) or coordinate a group of TCLs. We used reduced-order models to simulate the thermodynamics of TCLs, and the parameter values were determined from the connected smart thermostat data of real households. The environment was built upon the standardized OpenAI Gym interface. Ancillary functions, such as retrieving the parameters and weather forecasts, are provided to facilitate control strategies that require predictive information. Compared with existing efforts, AlphaBuilding ResCommunity has three advantages: (1) more realistic model settings because the parameter values are identified from actual household operating data, and modelling and measurement uncertainty are considered; (2) passive thermal storage control; and (3) ease of use due to a simple software dependency and standardized interface. We demonstrated the applications of the environment by implementing a Kalman Filter and Model Predictive Control on a single TCL and a Priority-Stack-Based Control and Alternating Direction Method of Multipliers to coordinate multiple TCLs for load tracking.

Wang, Z↗
DOE OSTI2022

Registration and fusion of large-scale melt pool temperature and morphology monitoring data demonstrated for surface topography prediction in LPBF

In-situ monitoring technologies for laser powder bed fusion (LPBF) additive manufacturing often face one key challenge, extracting the ultrafast melt pool (MP) signatures for understanding the localized part properties. Further, the spatial information of each monitored MP signature is essential for correlating the MP – part property. This spatial information is often unavailable especially from commercial LPBF printers. Many MP monitoring methods have been reported and utilized. However, very few of these have the MP’s spatial information. To overcome this challenge, in this work we report a method for spatially registering the key MP signatures (MP intensity, temperature, and area) to the monitored print parts. The MP signatures are obtained from our coaxial high-speed single-camera based two-wavelength imaging pyrometry (STWIP) system and the MP spatial information is obtained from an off-axis camera system. A machine learning aided image analysis method is employed to retrieve the spatial distribution of MPs within the corresponding part’s coordinates system. Then, the MP signature maps (MPSMs) are reconstructed by mapping the STWIP measured MP signatures to the registered MP coordinates. Further, a long short-term memory (LSTM) neural network is developed for estimating the layer surface topography from the registered MPSMs. The obtained results indicate that the layer surface topography can be more accurately estimated by using MP temperature signature rather than MP intensity and/or area signatures as in common practice. Finally, our developed methods for MP monitoring, registration, and MP-surface topography prediction offer advanced capabilities for the online detection of process anomalies and part defects.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A review of coarse mineral dust in the Earth system

Mineral dust particles suspended in the atmosphere span more than three orders of magnitude in diameter, from <0.1 µm to more than 100 µm. This wide size range makes dust a unique aerosol species with the ability to interact with many aspects of the Earth system, including radiation, clouds, hydrology, atmospheric chemistry, and biogeochemistry. This review focuses on coarse and super-coarse dust aerosols, which we respectively define as dust particles with a diameter of 2.5–10 µm and 10–62.5 µm. We review several lines of observational evidence indicating that coarse and super-coarse dust particles are transported farther than previously expected and that the abundance of these particles is substantially underestimated in current global models. We synthesize previous studies that used observations, theories, and model simulations to highlight the impacts of coarse and super-coarse dust aerosols on the Earth system, including their effects on dust-radiation interactions, dust-cloud interactions, atmospheric chemistry, and biogeochemistry. Specifically, coarse and super-coarse dust aerosols produce a net positive direct radiative effect (warming) at the top of the atmosphere and can modify temperature and water vapor profiles, influencing the distribution of clouds and precipitation. In addition, coarse and super-coarse dust aerosols contribute a substantial fraction of ice-nucleating particles, especially at temperatures above –23 °C. They also contribute a substantial fraction to the available reactive surfaces for atmospheric processing and the dust deposition flux that impacts land and ocean biogeochemistry by supplying important nutrients such as iron and phosphorus. Furthermore, we examine several limitations in the representation of coarse and super-coarse dust aerosols in current model simulations and remote-sensing retrievals. Because these limitations substantially contribute to the uncertainties in simulating the abundance and impacts of coarse and super-coarse dust aerosols, we offer some recommendations to facilitate future studies. Overall, we conclude that an accurate representation of coarse and super-coarse properties is critical in understanding the impacts of dust aerosols on the Earth system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Influence of sun zenith angle on canopy clumping and the resulting impacts on photosynthesis

Addressing the impact of vegetation architecture on shortwave radiation transfer in land surface models is important for accurate weather forecasting, carbon budget estimates, and climate predictions. This paper investigates to what extent it is possible to retrieve structural parameters of two different parameterization schemes from direct transmittance derived from digital hemispherical photography and 3D radiative transfer modeling for two study sites with different vegetation canopy architectures. Neglecting the representation of 3D canopy structure in radiative transfer schemes leads to significant errors in shortwave radiation partitioning (up to 3.5 times more direct transmittance in the 3D model). Structural parameters, referred to as whole canopy ‘clumping indices’, were obtained in order to evaluate the impact of angular variation in clumping on shortwave radiation transfer. Impacts on photosynthesis were evaluated at site level with the UKESM land surface model, JULES. A comparison between flux tower derived and modeled photosynthesis indicates that considering zenith angular variations of structural parameters in the radiative transfer scheme of the UKESM land surface model significantly improves photosynthesis prediction in light limited ecosystems (from RMSE = 2.91 μmol CO 2 .m -2 .s -1 to RMSE = 1.51 μmol CO 2 .m -2 .s -1 , 48% smaller), typically with enhanced photosynthesis from bottom layers.

54 ENVIRONMENTAL SCIENCES↗