Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mixing”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Mixing across stable density interfaces in forced stratified turbulence

Understanding how turbulence enhances irreversible scalar mixing in density-stratified fluids is a central problem in geophysical fluid dynamics. While isotropic overturning regions are commonly the focus of mixing analyses, we here investigate whether significant mixing may arise in anisotropic statically stable regions of the flow. Focusing on a single forced direct numerical simulation of stratified turbulence, we analyse spatial correlations between the vertical density gradient $\partial \rho /\partial z$ and the dissipation rates of kinetic energy $\epsilon$ and scalar variance $\chi$ , the latter quantifying scalar mixing. The domain is characterized by relatively well-mixed density layers separated by sharp stable interfaces that are correlated with high vertical shear. While static instability is most prevalent within the mixed layers, much of the scalar mixing is localized to the intervening interfaces, a phenomenon not apparent if considering local static instability or $\epsilon$ alone. While the majority of the domain is characterized by the canonical flux coefficient $\varGamma \equiv \chi /\epsilon =0.2$ , often assumed in ocean mixing parametrizations, extreme values of $\chi$ within the statically stable interfaces, associated with elevated $\varGamma$ , strongly skew the bulk statistics. Our findings suggest that current parametrizations of turbulent mixing may be biased by undersampling, such that the most common, but not necessarily the most significant, mixing events are overweighted. Having focused here on a single simulation of stratified turbulence, it is hoped that our results motivate a broader investigation into the role played by stable density interfaces in mixing, across a wider range of parameters and forcing schemes representative of ocean turbulence.

Couchman, Miles M. P. (ORCID:0000000246676829)↗
DOE OSTI2022

Influences of an entrainment–mixing parameterization on numerical simulations of cumulus and stratocumulus clouds

Abstract. Different entrainment–mixing processes can occur in clouds; however, a homogeneous mixing mechanism is often implicitly assumed in most commonly used microphysics schemes. Here, we first present a new entrainment–mixing parameterization that uses the grid mean relative humidity without requiring the relative humidity of the entrained air. Then, the parameterization is implemented in a microphysics scheme in a large eddy simulation model, and sensitivity experiments are conducted to compare the new parameterization with the default homogeneous entrainment–mixing parameterization. The results indicate that the new entrainment–mixing parameterization has a larger impact on the number concentration, volume mean radius, and cloud optical depth in the stratocumulus case than in the cumulus case. This is because inhomogeneous and homogeneous mixing mechanisms dominate in the stratocumulus and cumulus cases, respectively, which is mainly due to the larger turbulence dissipation rate in the cumulus case. Because stratocumulus clouds break up during the dissipation stage to form cumulus clouds, the effects of this new entrainment–mixing parameterization during the stratocumulus dissipation stage are between those during the stratocumulus mature stage and the cumulus case. A large aerosol concentration can enhance the effects of this new entrainment–mixing parameterization by decreasing the cloud droplet size and evaporation timescale. The results of this new entrainment–mixing parameterization with grid mean relative humidity are validated by the use of a different entrainment–mixing parameterization that uses parameterized entrained air properties. This study sheds new light on the improvement of entrainment–mixing parameterizations in models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Why is height-dependent mixing observed in stratocumulus clouds?

Recent aircraft measurements in stratocumulus clouds suggest that entrainment mixing is inhomogeneous (IM) near cloud top and homogeneous (HM) within the cloud. However, this proposed height-dependence of mixing transition is uncertain because of artifacts involved in the aircraft measurements. In this study, we use the Explicit Mixing Parcel Model to simulate mixing scenarios in stratocumulus clouds and reconstruct the virtual aircraft measurements to investigate the mixing signature. Results show that, from the aircraft-measurement perspective, the mixing signature always exhibits IM characteristic near cloud top and HM characteristic within cloud, independent of the types of the local entrainment-mixing process. The appearance of the vertical IM-to-HM transition is essentially a collective behavior of multiple parcels sampled at the same height, experiencing distinct entrainment-mixing-evaporation histories. This bulk view of mixing process, which is widely used for aircraft measurements, could lead to misinterpretations of the true mixing mechanism occurring in clouds. Our result underscores the limitations of using aircraft measurements to identify the entrainment-mixing mechanism at the process level.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

An index to characterize gas‐solid and solid‐solid mixing from average volume fraction fields

Abstract A mixing index based on solid volume fraction fields is developed for gas‐solid flows. Conventional mixing indices are based on particle realizations of granular mixing and are applicable to experimental data or discrete element method simulations. However, these indices cannot be used as‐is for multifluid models, and an index for characterizing mixing in gas‐solid flows from continuous fields is needed. The performance of the new mixing index is tested in two applications. The first is a 3D simulation of the mixing of biomass and sand in a fluidized bed reactor, and the second is a 2D simulation of binary particle segregation in a fluidized bed. The simulations are performed using OpenFOAM®. The mixing index is used to quantify gas‐solid mixing using solid volume fractions and solid‐solid mixing using solid fractions. The formulation of conventional mixing indices is extended to be used with solid volume fractions fields, and methods for performance improvement are presented.

09 BIOMASS FUELS↗
DOE OSTI2023

Estimating the Meridional Extent of Adiabatic Mixing in the Stratosphere Using Age-Of-Air

Wave-induced adiabatic mixing in the winter midlatitudes is one of the key processes impacting stratospheric transport. Understanding its strength and structure is vital to understanding the distribution of trace gases and their modulation under a changing climate. Age-of-air is often used to understand stratospheric transport, and this study proposes refinements to the vertical age gradient theory of Linz et al. (2021), https://doi.org/10.1029/2021JD035199. The theory assumes exchange of air between a well-mixed tropics and a well-mixed extratropics, separated by a transport barrier, quantifying the adiabatic mixing flux across the interface using age-based measures. These assumptions are re-evaluated and a refined framework that includes the effects of meridional tracer gradients is established to quantify the mixing flux. This is achieved, in part, by computing a circulation streamfunction in age-potential temperature coordinates to generate a complete distribution of parcel ages being mixed in the midlatitudes. The streamfunction quantifies the “true” age of parcels mixed between the tropics and the extratropics. Applying the revised theory to an idealized and a comprehensive climate model reveals that ignoring the meridional gradients in age leads to an underestimation of the wave-driven mixing flux. Stronger, and qualitatively similar fluxes are obtained in both models, especially in the lower-to-middle stratosphere. While the meridional span of adiabatic mixing in the two models exhibits some differences, they show that the deep tropical pipe, that is, latitudes equatorward of 15° barely mix with older midlatitude air. The novel age-potential temperature circulation can be used to quantify additional aspects of stratospheric transport.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Using Machine Learning to Predict Cloud Turbulent Entrainment–Mixing Processes

Different turbulent entrainment–mixing mechanisms between clouds and environment are essential to cloud–related processes; however, accurate representation of entrainment–mixing in weather/climate models still poses a challenge. This study exploits the use of machine learning (ML) to address this challenge. Four ML (Light Gradient Boosting Machine [LGB], eXtreme Gradient Boosting, Random Forest, and Support Vector Regression) are examined and compared. It is found that LGB performs best, and thus is selected to understand the impact of entrainment–mixing on microphysics using simulation data from Explicit Mixing Parcel Model. Compared with traditional parameterizations, the trained LGB provides more accurate microphysical properties (number concentration and cloud droplet spectral dispersion). The partial dependences of predicted microphysics on features exhibit a strong alignment with physical mechanisms and expectations, as determined by the interpreting method, thus overcoming the limitations of the “black box” scheme. The underlying mechanisms are that the smaller number concentration and larger spectral dispersion correspond to more inhomogeneous entrainment–mixing. Specifically, number concentration after entrainment–mixing is positively correlated with adiabatic number concentration and liquid water content affected by entrainment–mixing, and inversely correlated with adiabatic volume mean radius. Spectral dispersion after entrainment–mixing is negatively correlated with liquid water content affected by entrainment–mixing, turbulent dissipation rate and relative humidity of entrained air. Sensitivity analysis further suggests that number concentration is mainly determined by cloud microphysical properties whereas spectral dispersion is influenced by both cloud microphysical properties and environmental variables. The results indicate that the LGB scheme has the potential to enhance the representation of entrainment–mixing in weather/climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Self-similar Reynolds-averaged mechanical–scalar turbulence models for reshocked Richtmyer–Meshkov instability-induced mixing in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing incompressible turbulent Richtmyer–Meshkov instability-induced mixing in planar geometry derived in the small Atwood number limit are extended to construct models for reshocked Richtmyer–Meshkov mixing. In this study, the models are based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar variance S and its dissipation rate χ modeled either differentially or algebraically. The three- and four-equation models allow for a simultaneous description of mechanical and scalar mixing, i.e., mixing layer growth and molecular mixing. Mixing layer growth parameters and other physical observables were obtained explicitly as functions of the model coefficients and were used to calibrate the model coefficients. Here, the solutions for the singly shocked Richtmyer–Meshkov case for the mixing layer width and the turbulent fields are used to construct piecewise-continuous generalizations of these quantities for times after reshock. For generality, the post-reshock mixing layer width is not assumed to grow with the same power-law as the pre-reshock width, and an impulsive approximation applied to Rayleigh–Taylor instability growth is used to establish the expression for the post-reshock width. A four-equation model is then used to illustrate the spatiotemporal behavior of the mean and turbulent fields and late-time turbulent equation budgets across the mixing layer. The reference solutions derived here can provide systematic calibrations and better understanding of mechanical–scalar turbulence models and their predictions for reshocked Richtmyer–Meshkov instability-induced turbulent mixing in the very large Reynolds number limit.

42 ENGINEERING↗
DOE OSTI2024

A study on the scale dependence of mixing indices for Eulerian multiphase models

Abstract Mixing can vary based on the scale at which the system is observed, and a mixing index that can capture the features at different length scales is desirable. In this article, we analyze the scale dependence of the mixing indices developed for Eulerian multiphase models. Relevant length scales are distinguished by filtering solid fraction fields. The scale‐dependence study is first done on manufactured fields of solid fraction to assess the performance of the mixing indices. The study is extended to a two‐dimensional CFD simulation of the segregation of a bidisperse gas–solid mixture. The local mixing index performs well in capturing the spatial variation of mixing at different scales. The scale dependence of two global mixing indices is considered in the study, where the state of mixing is defined based on statistical measures. We demonstrate that the choice of measures influences the sensitivity of mixing indices to mixing at different scales.

Nagawkar, Barlev R.↗
DOE OSTI2024

Engineering advancements in microfluidic systems for enhanced mixing at low Reynolds numbers

Mixing within micro- and millichannels is a pivotal element across various applications, ranging from chemical synthesis to biomedical diagnostics and environmental monitoring. The inherent low Reynolds number flow in these channels often results in a parabolic velocity profile, leading to a broad residence time distribution. Achieving efficient mixing at such small scales presents unique challenges and opportunities. This review encompasses various techniques and strategies to evaluate and enhance mixing efficiency in these confined environments. It explores the significance of mixing in micro- and millichannels, highlighting its relevance for enhanced reaction kinetics, homogeneity in mixed fluids, and analytical accuracy. We discuss various mixing methodologies that have been employed to get a narrower residence time distribution. The role of channel geometry, flow conditions, and mixing mechanisms in influencing the mixing performance are also discussed. Various emerging technologies and advancements in microfluidic devices and tools specifically designed to enhance mixing efficiency are highlighted. We emphasize the potential applications of micro- and millichannels in fields of nanoparticle synthesis, which can be utilized for biological applications. Additionally, the prospects of machine learning and artificial intelligence are offered toward incorporating better mixing to achieve precise control over nanoparticle synthesis, ultimately enhancing the potential for applications in these miniature fluidic systems.

Biochemistry & Molecular Biology↗
DOE OSTI2026

Exploring composition mixing in kilonova ejecta with ray-by-ray simulations

Binary neutron star merger (BNSM) ejecta are considered a primary repository of r-process nucleosynthesis and a source of the observed heavy-element abundances. We implement composition mixing into ray-by-ray radiation-hydrodynamic simulations of BNSM ejecta, coupled with an online nuclear network (NN). We model mixing via a gradient-based mixing approximation that evolves simultaneously with the hydrodynamics. Here, we find that mixing occurs in regions where the electron fraction changes rapidly. While mixing smooths composition gradients in transition regions, it has a negligible impact on the heavy-element yields. This is because the primary r-process site (the equatorial ejecta) is initially homogeneous in free neutrons, leaving no strong gradients for mixing to act upon. In each angular ray, the abundances of the most produced elements are robust under mixing, while the less abundant ones are more affected. The total global abundances change only slightly from mixing, since each angular ray contributes its most abundant elements. Furthermore, the predicted kilonova light curves show only minor reddening, with differences below the detectability of state-of-the-art telescopes. In general, we do not observe significant effects from mixing in the time span of the r-process. Consequently, mixing only leads to minor variations in abundances and light curves in ray-by-ray simulations.

Explosive burning↗
DOE OSTI2026

Upper ocean mixing, surface heat fluxes, and heat content variability in the upper 150 m during Hurricane Laura (2020)

Observations (from thermistor sensors and a hurricane glider) and model-based temperature were used to answer questions including the following: 1) What physical and environmental conditions explain the rapid intensification of Hurricane Laura? 2) How did the pre-existing warm mixed layer at Stone mooring (StM) influence the degree of cooling in the mixed layer? 3) How did vertical mixing, surface heat fluxes, and advective processes collectively shape the mixed layer heat evolution at StM? Observations measured by thermistor sensors showed that Hurricane Laura induced a cooling of -1.2oC in the mixed layer when Hurricane Laura’s wind speed increased to 47 m/s on 26 August. This contrasts with the -1.04oC mixed layer temperature change estimated from model-based temperature. At StM, the presence of a 31oC warm mixed layer and elevated heat content (60-80kJ/cm2) effectively preconditioned the upper ocean ahead of Hurricane Laura’s passage. After Hurricane Laura passed, turbulent mixing associated with the storm transported the heat anomalies downward into subsurface layers. In addition, this study shows that a surface heat flux of 5.04 kJ/cm2 supported the intensity of Hurricane Laura as it traversed over the StM. Using a mixed layer heat balance model, this study shows that entrainment flux and surface flux mostly contributed to the observed change in the mixed layer temperature.

Tsai, Senam↗
DOE OSTI2021

A Deep Learning Modeling Framework to Capture Mixing Patterns in Reactive-Transport Systems

Prediction and control of chemical mixing are vital for many scientific areas such as subsurface reactive transport, climate modeling, combustion, epidemiology, and pharmacology. Due to the complex nature of mixing in heterogeneous and anisotropic media, the mathematical models related to this phenomenon are not analytically tractable. Numerical simulations often provide a viable route to predict chemical mixing accurately. However, contemporary modeling approaches for mixing cannot utilize available spatial-temporal data to improve the accuracy of the future prediction and can be compute-intensive, especially when the spatial domain is large and for long-term temporal predictions. To address this knowledge gap, in this work we will present in this paper a deep learning (DL) modeling framework applied to predict the progress of chemical mixing under fast bimolecular reactions. This framework uses convolutional neural networks (CNN) for capturing spatial patterns and long short-term memory (LSTM) networks for forecasting temporal variations in mixing. By careful design of the framework—placement of non-negative constraint on the weights of the CNN and the selection of activation function, the framework ensures non-negativity of the chemical species at all spatial points and for all times. Our DL-based framework is fast, accurate, and requires minimal data for training. The time needed to obtain a forecast using the model is a fraction (≈ O(-6)) of the time needed to obtain the result using a high-fidelity simulation. To achieve an error of 10% (measured using the infinity norm) for capturing local-scale mixing features such as interfacial mixing, only 24% to 32% of the sequence data for model training is required. To achieve the same level of accuracy for capturing global-scale mixing features, the sequence data required for model training is 64% to 70% of the total spatial-temporal data. Hence, the proposed approach—a fast and accurate way to forecast long-time spatial-temporal mixing patterns in heterogeneous and anisotropic media—will be a valuable tool for modeling reactive-transport in a wide range of applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Single stage clarifier and mixing assembly

A single stage clarifier and mixing assembly includes a housing, a mixing section within the housing and a clarifier section within the housing. The mixing section includes a mixing chamber having (a) an inlet, adapted for delivering an inlet stream to the mixing chamber, at an upper end, and (b) a mixing section outlet at a lower end. The clarifier section extends concentrically around the mixing section. The single stage clarifier and mixing assembly also includes an agitator adapted for mixing the inlet stream in the mixing chamber.

Werner, Joshua Michael↗
DOE OSTI2025

Implementation and Evaluation of Physics-Driven Dynamic Entrainment-Mixing Parameterization in a Climate Model and Its Impact on Low-Cloud Simulation

The turbulent entrainment-mixing process in the Community Earth System Model version 1.2 (CESM1.2) is assumed to follow the extremely inhomogeneous entrainment-mixing. However, different entrainment-mixing scenarios can occur in real clouds. To address this deficiency, a unifying parameterization that represents different entrainment-mixing processes is implemented and evaluated in CESM1.2. The results indicate that the homogeneous mixing degree values simulated by the new parameterization in CESM1.2 are predominantly greater than 50%, suggesting a tendency toward homogeneous mixing. Compared to the extremely inhomogeneous mixing mechanism, the new parameterization increases the cloud droplet number concentration (Nc). More importantly, the new parameterization improves low-cloud fraction (CLDLOW) simulation in Northwest Pacific (NWP) and Southeast Pacific (SEP) regions, with relative improvements of 2.95% and 4.17%, respectively. Furthermore, the improvements reach up to 44.6% and 16.2% in the NWP and SEP regions, respectively, when considering the relationship between N c and CLDLOW. Further analysis reveals that the new parameterization enhances cloud optical depth, longwave radiative cooling effect, net condensation rate, cloud water mixing ratio, lower-troposphere stability, and CLDLOW by increasing N c . Additionally, these results underscore the importance of improving entrainment-mixing parameterization in climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Little prairie under the panel: testing native pollinator habitat seed mix establishment at three utility-scale solar sites in Minnesota

Abstract As more land is being utilized for large-scale solar energy projects, there are increasing discussions from stakeholders on how to utilize land under solar panels to promote biodiversity. One path is to plant habitat beneficial to pollinators and other insects, but there have been few long-term studies that examine how different vegetation and seed mixes establish underneath solar panels. This study addresses a scientific gap to determine whether native pollinator seed mixes successfully establish over time under solar arrays using a systematic assessment of eight seed mixes planted at three utility-scale solar sites in Minnesota. We assess establishment with a percent native coverage metric, which is an assessment of native species observations compared to total observations during percent cover analyses in our vegetative test plots. The percent native coverage metric allows for a measurement of how the seed mix established and how the seed mix persists over time. The percent native coverage under and in between the solar photovoltaic (PV) arrays rose from 10% after one year of planting to 58% after three years across all sites, while the native coverage of the full sun control area rose from 9.6% to 70% under the same period, showing that native prairie and pollinator plants successfully established under the array, although to a lesser extent than in full sun conditions. Percent native coverage under the PV arrays rose 5- to 8-fold for each of the three sites from over the course of the study, while the coverage of weeds decreased for all three sites over the same period. Percent native coverage varied by seed mix over the project years, but every seed mix experienced a higher percent native coverage year after year under the PV arrays. Our results did not indicate a difference in establishment across placement within the array; the center, west, and east portions of the areas in between panels had similar establishment rates at two out of three sites, indicating that the same seed mix can be applied throughout the array. Out of 101 plant species seeded, we observed the establishment of 68 species in our vegetative test plots, and we detailed the top 20 observed species to inform future seed mix development. Based on these findings, native pollinator vegetation can establish over time at solar arrays, and it can be suitable for creating habitat at utility-scale solar sites.

14 SOLAR ENERGY↗
DOE OSTI2022

Mixing Pathways in Simple Box Models

Abstract Tracer variance budgets can be used to estimate bulk mixing in a control volume. For example, simple, analytical, bulk formulations of salt mixing, defined here as the destruction of salinity variance, can be found for estuaries with a riverine source of freshwater and a two-layer exchange flow at the mouth using salinity as a representative tracer. For a steady case, the bulk salt mixing M can be calculated as , where S in and S out are the representative salinities in the estuarine exchange flow, and Q R and Q in are the river and landward volume fluxes, respectively. The bulk salt mixing M can be considered as the sum of mixing pathways, where each pathway has a mixing of Q (Δ S ) 2 , where Q is the volume transport and Δ S is the salinity difference across the pathway. For the estuary case, one mixing path is associated with the river inflow, and the other is associated with the inflow of salty, oceanic water. This concept of linking mixing to input–output pathways is extended, in simple box models, from estuaries to scenarios with multiple inputs/outputs, as might be found in a complex estuarine/fjord network, in a region on a continental shelf, or any other control volume with multiple exchanges. This approach allows for the estimation of the relative contributions of each input–output pathway to the total mixing within a control volume.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Destabilization of binary mixing layer in supercritical conditions

Compressible mixing layer instabilities are of importance to a wide range of environmental and industrial applications. Past studies have focused on either ideal-gas or real-fluid thermodynamic regimes of single-species mixing layers. However, mixing layers of binary mixtures at supercritical conditions, commonly encountered in fuel injection systems, introduce additional complexities due to the added compositional degree-of-freedom. Moreover, the effect of strong variations in thermodynamic response functions across the Widom line on the binary mixing layer stability remains poorly understood. Thus, the objective of this study is to examine the coupling between the hydrodynamic instability and the real-fluid thermodynamics across the Widom line and its effects on the overall binary mixing layer dynamics. To this end, we develop a linear stability analysis of the full binary-species compressible transport equations coupled with the PC-SAFT equation of state. Analysis shows the existence of a novel instability mechanism that arises from juxtapositioning of the Widom-line transition and the hydrodynamic inflection point. This novel thermodynamically induced instability mechanism has the net effect of destabilizing the binary mixing layer at lowering supercritical conditions towards the critical pressure point. This is in contrast to previous stability analyses of supercritical single-species mixing layers, where increasing pressure destabilizes the flow due to its effect on reducing the density stratification. The discovered thermodynamically induced instability mechanism of binary mixing flows highlights the need for an extension of classical instability criteria to incorporate the effect of strong variations in the thermodynamic response functions across the Widom line on mixing layer instability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Consideration of Initial Cloud Droplet Size Distribution Shapes in Quantifying Different Entrainment-Mixing Mechanisms

Entrainment-mixing process significantly affects cloud micro/macrophysics and the evaluation of aerosol indirect effects. However, it remains an open question as to how to quantify entrainment-mixing mechanism for broad cloud droplet size distributions (CDSDs). Here, CDSDs with different spectral widths are used to initialize the Explicit Mixing Parcel Model to address this problem, and microphysical properties are compared for different CDSD widths. For relatively broad CDSDs, as the number concentration and liquid water content decrease, the volume-mean radius and mean radius increase, which is opposite to the scenario for relatively narrow CDSDs and the homogeneous/inhomogeneous conceptual model. For the relatively broad CDSDs, such relationships could be mistaken as inhomogeneous mixing with subsequent ascent in analysis of the conventional microphysical mixing diagram. This then causes traditional homogenous mixing degrees to be unrealistically negative and not applicable for relatively broad CDSDs. New measures are introduced to explicitly account for the effect of CDSD spectral shapes on quantifying entrainment-mixing mechanism. The new measures yield reasonable ranges of values that conform to the dynamical measures such as the Damköhler and transitional scale numbers, and their temporal variation can be explained physically. This study extends the measures of homogeneous mixing degrees from narrow to broad CDSDs, and sheds new light on the consideration of CDSD shapes in parameterizations of entrainment-mixing mechanism.

54 ENVIRONMENTAL SCIENCES↗