Search NASA⌕ Search

SEARCH · Search NASA

Results for “MIXING RATIO”

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 91 records · Page 5

Printing noble metal alloy films with compositional gradient

Additive manufacturing has offered great promise in fabricating materials and devices with complex structures and unique properties. In addition to complex structures, additive manufacturing also allows for local changes to material compositions and/or structure, ideal for combinatorial material investigations. Besides, advantages over freedom of design, mass customization, waste minimization and the capability to fast prototyping, the on-demand control of compositional distributions in low dimensional materials remains challenging. In this work, we demonstrate a combinatorial aerosol jet printing technique to create noble metal films with gradient compositions. Nanoparticle sizes and ink chemistry are critical for aerosol jet printing and sintering process in order to achieve defect free thin films. Ruthenium (Ru) and platinum (Pt) nanoparticles with size <5 nm were synthesized and formulated into printable inks. Further, after thermal sintering at temperature as low as 400 °C, printed nanoparticles formed highly reflective films free from microstructural defect. Alloy films of gradient compositions were realized via in-situ mixing Ru and Pt ink streams and varying the mixing ratio during the printing process. This combinatorial printing method provides great potential to rapidly transform nanoparticle inks into new materials compositions and structures unobtainable via conventional methods for high-throughput materials studies.

36 MATERIALS SCIENCE↗

Model-based optimization strategies for direct hydrogenation of carbon dioxide to dimethyl ether

Here, this study discusses model-based optimization strategies for CO 2 hydrogenation to dimethyl ether (DME) over aCuZnZr (CZZ) and ferrite (FER) mixed catalyst system in a packed-bed reactor configuration. A two-dimensional axisymmetric, nonisothermal packed-bed reactor model was developed using COMSOL Multiphysics 6.2 software. The model solves two-dimensional (radial and axial) heat and mass transport equations in the packed-bed and integrates intraparticle diffusion and heat transfer in a 1D approach. This powerful feature differs from a traditional porous media approach and takes into account any heat and mass transfer limitations that may exist. Analysis shows that the heat transfer limitations are negligible, but strong internal mass transfer limitations were observed at 10 ≤ WHSV ≤ 90 h -1 on the FER catalyst and at 240 °C. The optimum catalyst composition (i.e., mixing ratio) varies depending on the operating regime. The FER catalyst weight in the mixture can be as low as 5 wt.%, but the ideal composition de-pends on the internal mass transfer limitation and its relationship with the operating regime (i.e., weight hourly space velocity, temperature). A catalyst composition of 80 wt.% CZZ and 20 wt.% FER was suggested; this composition can provide high CO 2 conversion and DME production rates at a wide range of temperatures and flow rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A general predictive methodology for fuel-mixture properties up to supercritical conditions

A predictive thermodynamic model is utilized for the calculation of fuel properties of oxymethylene dimethyl ethers (OME 3–4 ), surrogates for gasoline, diesel and aviation fuel, as well as alcohol blends with gasoline and diesel. The alcohols used for these blends are methanol, ethanol, propanol, butanol and pentanol; their mixing ratio ranges from 10 to 50% by volume. The model is based on the Perturbed-Chain Statistical Association Fluid Theory (PC-SAFT) equation of state (EoS) and Vapor Liquid Equilibrium (VLE) calculations at constant temperature, density and composition. The model includes the association term, with the assumption of two association sites (2B scheme), to enable the modeling of alcohols. The pure-component parameters are estimated based on the Group Contribution (GC) method of various sources, as well as a parametrization model specifically designed for the case of OME 3–4 . The results of the computational model for the density, vapor pressure and distillation curves at various conditions, including high-pressure, high-temperature (HPHT), are compared to experimental and computational data available in the literature. In the cases where no measurements are available for the surrogates, experimental data for the corresponding target fuel are used, taking into consideration the inherent deviation in properties between real and surrogate fuel. Overall, the results are in good agreement with the data from the literature, with the average deviation not exceeding 12% for temperature (Kelvin) on the distillation curves, 10% for density and 46% for vapor pressure and the general trend being captured successfully. The use of different pure component parameter estimation techniques can further improve the prediction quality in the cases of OME3–4 and the aviation fuel surrogate, especially for the vapor pressure, leading to an average deviation lower than 18%. These results demonstrate the predictive capabilities of the model, which extend to a wide range of fuel types and pressure/temperature conditions. Through this investigation, the present work aims to establish the limits of applicability of this thermodynamic property prediction methodology.

33 ADVANCED PROPULSION SYSTEMS↗

A method for generating quantitative vapor-phase infrared spectra of solids: results for phenol, camphor, menthol, syringol, dicyclopentadiene and naphthalene

Here, a method is presented to generate quantitative vapor-phase infrared spectra from substances that naturally occur as solids with moderate volatility. The solid is gravimetrically dissolved into a solvent that has few infrared spectral features, typically CS 2 and CCl 4 separately. The solution is flowed at a constant rate from a linearly pumped syringe into a metered stream of nitrogen carrier gas regulated by a mass flow controller. The analyte/solvent mix is flash vaporized by volatilizing the solution across a heated stainless-steel surface as it emanates from the syringe tip. The N 2 gas-solution mixture is flowed into a long-path White cell thermostatted at a desired temperature, the long optical path compensating for the modest analyte mixing ratio. A composite spectrum is generated from typically ten or more 760-Torr pressure-broadened spectra over the 600 to 6500 cm -1 spectral range at 0.1 cm -1 spectral resolution. The solid analytes reported here using this novel technique include dicyclopentadiene, menthol, syringol, phenol, camphor, and naphthalene.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Evidence against the wobbling nature of low-spin bands in 135 Pr

The electromagnetic character of the ΔI=1 transitions connecting the 1- to 0-phonon and the 2- to 1-phonon wobbling bands should be dominated by an E2 component, due to the collective motion of the entire nuclear charge. In the present work it is shown, based on combined angular correlation and linear polarization measurements, that the mixing ratios of all analyzed connecting transitions between low-lying bands in 135 Pr interpreted as 0-, 1-, and 2-phonon wobbling bands, have absolute values smaller than one. This indicates predominant M1 magnetic character, which is incompatible with the proposed wobbling nature. All experimental observables are instead in good agreement with quasiparticle-plus-triaxial-rotor model calculations, which describe the bands as resulting from a rapid re-alignment of the total angular momentum from the short to the intermediate nuclear axis

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Development of deformed intruder structures along the N=50 magic number: Insights from 83 Se

The 83 Se nucleus was studied in the context of shape coexistence employing a 82 Se(d,p) 83 Se reaction performed in inverse kinematics with the GALILEO-II + SPIDER setup to reach the desired selectivity on the reaction channel of interest via γ -particle coincidences. The lifetimes of the 1/2 + 540-keV and 3/2$^{+}_{3}$ 1100-keV low-lying intruder states of 83 Se were measured for the first time by employing the Recoil Distance Doppler Shift and the Doppler Shift Attenuation Method. The deduced transition probabilities, together with the mixing ratios obtained from a separate neutron-capture experiment performed at ILL with the FIPPS array, are compared and discussed in the framework of shell-model calculations performed with the PFSDG-U Hamiltonian. Large transition probabilities of hundreds of e 2 fm 4 were found for the E 2 decays of the states built on top of the intruders, pointing to the development of quadrupole collectivity and large deformation in the proximity of N = 50 shell closure.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

An efficient method to identify uncertainties of WRF-Solar variables in forecasting solar irradiance using a tangent linear sensitivity analysis

Uncertainty in predicting solar energy resources introduces major challenges in power system management and necessitates the development of reliable probabilistic solar forecasts. As the first part of the development of probabilistic forecasts based on the Weather Research and Forecasting model with solar extensions (WRF-Solar), this study presents a tangent linear approach to identify input variables responsible for the largest uncertainties in predicting surface solar irradiance and clouds. A tangent linear analysis is capable of efficiently investigating sensitivities of output variables with respect to various input variables of WRF-Solar because this approach avoids the computational burden of perturbing the initial conditions of individual input variables. We develop tangent linear models (TLMs) for six WRF-Solar physics packages that control the formation and dissipation of clouds and solar radiation, and we evaluate the validity of TLMs using a linearity test. The tangent linear sensitivity analysis is conducted under various scenarios based on satellite observations and model simulations to consider realistic input conditions. A simple method is used to quantify the impact of the uncertainty of input variables on the output variables from the TLMs. The results demonstrate that uncertainties in the output variables that are the focus of this study—including global horizontal irradiance, direct normal irradiance, cloud mixing ratio, cloud tendency, cloud fraction, and sensible and latent heat fluxes—are highly sensitive to uncertainties in 14 input variables. This study indicates that the tangent linear method can identify key variables of physics modules in WRF-Solar that can be stochastically perturbed to generate ensemble-based probabilistic forecasts.

14 SOLAR ENERGY↗

Rapid Droplet Sampling Interface for Low-Volume, High-Throughput Mass Spectrometry Analysis

Here, we present a rapid droplet sampling interface (RDSI) electrospray ionization mass spectrometry (ESI-MS) system as a high-throughput, low-volume, noncontact, and minimal-carryover approach for characterization of liquids. Liquid characterization was achieved by combining droplet ejection with an open-face microflow capillary with a 2.5 μL/min continuous flow of carrier solvent. Through this implementation, single 0.3 nL droplets containing the analyte effectively mix with 4–8 nL of carrier solvent and create a combined electrospray plume. The carrier solvent continuously cleaned the system, eliminating carryover. A sampling rate of 5 Hz was achieved for droplets containing 1 μM propranolol or 5 μM leu-enkephalin with each droplet fully baseline-resolved (138 ± 32 ms baseline-to-baseline). Using a SCIEX API4000 mass spectrometer, a lower limit of quantification (LLOQ) of propranolol was 15 nM, corresponding to 1.16 fg of propranolol in the droplet, and was linear across 3 orders of magnitude. Further, quantitation could be achieved by adding an isotopically labeled internal standard, as done in conventional ESI. Signal transients were faster than the acquisition speed of the mass spectrometer, resulting in artificially high reproducibility of 15–30% RSD droplet-to-droplet. Analyte–solvent mixing ratios could be controlled by adjusting droplet positioning along the open-face capillary with an optimal position about 0.4 mm from the tip end. The range of analyte coverage was exemplified by measures of peptides and drugs in methanol, water, and buffer solutions. In a comparison to the Open Port Sampling Interface (OPSI) implemented on the same system, the RDSI had 78× greater sensitivity, 6× greater throughput and used significantly less carrier solvent.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Inorganic Alterations in Unconventional Shale Reservoirs: Importance of Additive and Base Fluid Chemistry

The effective development of unconventional petroleum systems requires the use of significant water resources. In an effort to reduce the consumption of freshwater resources for hydraulic fracturing, highly saline produced waters are increasingly recycled for use as a base fluid. However, there are significant knowledge gaps regarding potential water–rock interactions resulting from the introduction of produced waters and associated additives into shale reservoirs such as formation and deposition of mineral scale, which can negatively affect hydrocarbon production through wellbore restriction and damage to hydraulically generated fractures. To assess the impacts of field stimulation practices in the subsurface, a series of laboratory experiments were completed using (a) three distinct sedimentary rock formations of the Midland Basin (Texas, USA) and (b) additives with two different base fluids: municipal fresh water and clean brine. The experimental approach used relevant injection sequences and mixing ratios in specialized reactors for 3 weeks. Static pressurized experiments and nonpressurized time-resolved experiments were undertaken. The resulting solids and liquids were analyzed by using a variety of laboratory- and synchrotron-based techniques. The use of an acid spearhead (15% HCl) resulted in texturing of both clay-rich and calcareous shales, which can temporarily enhance porosity but subsequently result in mineral scale deposition. The primary matrix scale was Fe(III)-bearing phases, which occurred in all experiments regardless of base fluid chemistry. Additionally, strontium sulfate (SrSO 4 ) precipitated on shale surfaces when clean brines were used. It was concluded that clean brine was the main source of Sr 2+ species, while persulfate breaker degradation and oxidation of pyrite were the sources of SO 4 2– . Sulfate scaling was more pronounced in clay-rich shales, suggesting that Sr sorption is important for promoting celestite formation. This work demonstrates that mineral scale deposition is a complex phenomenon, whereby the type and proportions of various mineral phases are determined from reservoir alteration processes and coprecipitation of constituents from injection fluids. In conclusion, the experimental results shown here should be considered when evaluating different base fluids and additives in order to mitigate mineral precipitation in unconventional shale reservoirs, which could result in reservoir degradation.

Jew, Adam D. [SLAC National Accelerator Laboratory↗

Kinetics of Direct Olefin Synthesis from Syngas over Mixed Beds of Zn–Zr Oxides and SAPO-34

A packed bed containing a physical mixture of both Zn–Zr mixed oxide catalyst and SAPO-34 converts syngas directly into a mixture of C 2 –C 5 olefins and paraffins. Specifically, the mixed oxide catalyst is responsible for intermediate oxygenate synthesis from syngas while the molecular sieve catalyzes olefin synthesis from the oxygenate intermediates. Kinetic measurements with cofed propylene over each catalyst independently confirm olefin hydrogenation activity over both components of the composite bed. The addition of either water or CO to the feed drops the activity of propylene hydrogenation over the Zn–Zr oxide. In sum, under reaction conditions of syngas feed and produced water, olefin hydrogenation predominantly occurs over the SAPO-34 catalyst, rather than over the catalyst responsible for hydrogenating CO into oxygenate intermediates. Here, a developed kinetic model consistent with this conclusion describes measurements at differing feed compositions, temperatures, space velocities, and bed catalyst mixing ratios. Technoeconomic analysis of the process indicates that the olefin-to-paraffin ratio is a key performance metric for commercial scale syngas conversion and highlights the importance of considering olefin hydrogenation rates over the molecular sieve component.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Understanding Lithium Local Environments in LiMn 0.5 Ni 0.5 O 2 Cathodes: A DFT-Supported 6 Li Solid-State NMR Study

LiMn 0.5 Ni 0.5 O 2 has been plagued by the structural failure driven by the irreversible phase change upon subsequent cycling and the structural degradation caused by the seemingly unavoidable Li/Ni exchange. This begs for a better understanding of the correlation between lithium local structures, synthesis conditions, and overall electrochemical properties of materials. In this study, we use density functional theory (DFT) to assist experimental 6 Li solid-state magic angle spinning (MAS) nuclear magnetic resonance (NMR) spectroscopy to understand the nature and change in lithium local structures with extended annealing times (168 vs 15 h) synthesized at 900 °C. To correlate the NMR peak changes after different synthesis conditions, we calculated the NMR spectra of representative low-energy model configurations of LiMn 0.5 Ni 0.5 O 2 composition–zigzag, row, and flower structures–as well as the influence of Ni/Mn ordering and Li/Ni exchange on each spectrum using a supercell that is large enough to allow us to consider the effect of different relative configurations of adjacent transition metal (TM) layers. The analysis based on a combination of the calculated NMR and the experimental spectra suggests that most configurations contributing to the experimental spectra are in fact composed of a blend of configurations higher in energy than the classical well-ordered ground-state structures. Further, the extended annealing of LiMn 0.5 Ni 0.5 O 2 slightly enhances Ni/Mn and Li/Mn orderings while reducing the Li/Ni mixing ratio. This study shows how DFT calculations are crucial in providing a better understanding of the experimental characterization data and therefore local environments and domain structures in oxide materials, such as Li-ion battery cathodes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Defect-Limited Mobility and Defect Thermochemistry in Mixed A-Cation Tin Perovskites: (CH3NH3)1-xCsxSnBr3

Hybrid organic-inorganic semiconductors crystallizing in the perovskite structure present a significant opportunity for realizing defect-tolerant semiconductors. In this work, we examine the solid solution, (CH3NH3)1-xCsxSnBr3, and identify the thermochemistry dictating the intrinsic carrier concentrations and how local structural distortions influences this electronic behavior. This family of compounds exhibits the expected systematic trend in decreasing optical gap with the cesium to methylammonium A-site mixing ratio, x, in the visible region. However, the carrier mobility, as determined from time-resolved microwave conductivity measurements, trends opposite to that expected from first-principles calculations combined with Boltzmann scattering theory calculations of the carrier mobility. We propose that this is a result of increasing carrier scattering with x in (CH3NH3)1-xCsxSnBr3 due to a significant increase in the carrier density with x. By examining the dependence of the carrier density as a function of x, we infer the compositional-dependence of the average enthalpy and nonconfigurational entropy per defect. While diffraction reveals a cubic aristotypic perovskite structure as a function of x at room temperature, the pair distribution functions obtained from synchrotron X-ray total scattering from these materials are better described by symmetry-adapted displacement modes of the Pm3m crystal structure, which we attribute to a large degree of anharmonic dynamics of the atom positions. This analysis shows that CH3NH3+-rich compositions retain mostly linear Sn-Br-Sn bonding environments through the displacements, while Cs+-rich compositions lead to more significantly bent Sn-Br-Sn environments. We propose that these bent bonding arrangements in Cs-rich compositions yield a higher propensity for defect formation. This also provides a rationale for carrier trapping that gives rise to anomalous microwave transients. Together, these results provide insight into the structure-dynamics-properties relationships in this highly anharmonic system with high amplitude atomic motions and low defect formation energies.

defect thermochemistry↗

Polariton Dynamics in Two-Dimensional Ruddlesden–Popper Perovskites Strongly Coupled with Plasmonic Lattices

Strong coupling between light and matter can produce hybrid eigenstates known as exciton-polaritons. Although polariton dynamics are important photophysical properties, the relaxation pathways of polaritons in different coupling regimes have seen limited attention. This paper reports the dynamics of hybridized states from 2D Ruddlesden-Popper perovskites coupled to plasmonic nanoparticle lattices. The open cavity architecture of Al lattices enables the coupling strength to be modulated by varying either the lead halide perovskite film thickness or the superstrate refractive index. Both experiments and finite-difference time-domain simulations of the optical dispersion diagrams showed avoided crossings that are a signature of strong coupling. Furthermore, our analytical model also elucidated the correlation between the exciton/plasmon mixing ratio and polariton coupling strength. Using fs-transient absorption spectroscopy, we found that both the upper and lower polaritons have shorter lifetimes than the excitons and that polaritons can show faster excited state dynamics when they have access to additional energy transfer channels.

2D perovskite↗

Drizzle, Turbulence, and Density Currents Below Post Cold Frontal Open Cellular Marine Stratocumulus Clouds

Nine cases of postcold frontal marine stratocumulus clouds exhibiting open cellular mesoscale organization are analyzed to characterize the drizzle, turbulence, and density currents below them. Data collected by the vertically pointing Doppler cloud radar and multiple lidars part of the Atmospheric Radiation Measurement (ARM) Eastern North Atlantic (ENA) site are used in these analyses. A total of 251 drizzle shafts passed over the site with 76 shafts sampled for more than 15 min by the vertically pointing instruments. On average the drizzle shafts were ~18 km wide with below-cloud drizzle water path of ~40 g m −2 , cloud base rain rate of 7.78 mm day −1 , and cloud base drizzle modal diameter of 393 μm. The widths of the drizzle shafts did not exhibit any relationship with the cloud base rain rate, below-cloud drizzle water path, and the total water removed from the cloud. The downdrafts in the lowest 500 m within the drizzle cells strengthened with increasing cloud base rain rates, while the updrafts within drizzle cells did not exhibit this behavior. On average, the surface air density, pressure, water vapor mixing ratio, and wind speed along the background wind during the drizzle shaft were higher by 1.88 g m −3 , 8.60 Pa, 0.24 g kg −1 , and 1.22 m s −1 , respectively, from their background values. On the other hand, the surface air temperature was lower by 0.43 K during the drizzle shaft than its background value. The coincident measurements highlight the drizzle-turbulence-surface coupling in these cloud systems.

59 BASIC BIOLOGICAL SCIENCES↗

An Overview of Atmospheric Features Over the Western North Atlantic Ocean and North American East Coast - Part 1: Analysis of Aerosols, Gases, and Wet Deposition Chemistry

The Western North Atlantic Ocean (WNAO) and adjoining East Coast of North America are of great importance for atmospheric research and have been extensively studied for several decades. This broad region exhibits complex meteorological features, diverse cloud types, and a wide range of conditions associated with gas and particulate species that originate from many sources regionally and from other continents. As Part 1 of a 2-part paper series, this work provides a characterization of quantities associated with atmospheric chemistry, including gases, aerosols, and wet deposition, based on satellite observations, ground-based data, model simulations, and reanalysis products. Key results include spatial and seasonal differences in composition along the North American East Coast and over the WNAO associated with varying sources of smoke and dust and meteorological drivers such as temperature, moisture, and precipitation. Spatial and seasonal variations of tropospheric carbon monoxide and ozone highlight different pathways towards accumulation of these species in the troposphere. Spatial distributions of speciated aerosol optical depth and vertical profiles of aerosol mass mixing ratios show a clear seasonal cycle highlighting the influence of both natural and anthropogenic sources in addition to the impact of intercontinental transport on their vertical profiles. Analysis of long-term climate model simulations of aerosol species and satellite observations of carbon monoxide confirm that there has been a significant decline in recent decades among anthropogenic constituents owing to regulatory activities. The results can inform current and future lines of research for the region.

Corral, Andrea↗

A Lagrangian perspective on tropical anvil cloud lifecycle in present and future climate

The evolution of tropical anvil clouds from their origin in deep convective cores to their slow decay determines the climatic effects of clouds in tropical convective regions. Despite the relevance of anvil clouds for climate and responses of clouds to global warming, processes dominating their evolution are not well understood. Currently available observational data reveal instantaneous snapshots of anvil cloud properties, but cannot provide a process-based perspective on anvil evolution. We therefore conduct simulations with the high resolution version of the Exascale Earth System Model in which we track mesoscale convective systems over the Tropical Western Pacific and compute trajectories that follow air parcels detrained from peaks of convective activity. With this approach we gain new insight into the anvil cloud evolution both in present day and future climate. Comparison with geostationary satellite data shows that the model is able to simulate maritime mesoscale convective systems reasonably well. Trajectory results indicate that anvil cloud lifetime is about 15 hours with no significant change in a warmer climate. The anvil ice mixing ratio is larger in a warmer climate due to a larger source of ice by detrainment and larger depositional growth leading to a more negative net cloud radiative effect along detrained trajectories. However, the increases in sources are counteracted by increases in sinks of ice, particularly snow formation and sedimentation. Furthermore, we find that the mean anvil cloud feedback along trajectories is positive and consistent with results from more traditional cloud feedback calculation methods.

climate, weather, mesocale cellular convection↗

Evaluation of the CAM6 Climate Model Using Cloud Observations at McMurdo Station, Antarctica

A comparative analysis between observational data from McMurdo Station, Antarctica and the Community Atmosphere Model version 6 (CAM6) simulation is performed focusing on cloud characteristics and their thermodynamic conditions. Ka-band Zenith Radar (KAZR) and High Spectral Resolution Lidar (HSRL) retrievals are used as the basis of cloud fraction and cloud phase identifications. Radiosondes released at 12-hour increments provide atmospheric profiles for evaluating the simulated thermodynamic conditions. Our findings show that the CAM6 simulation consistently overestimates (underestimates) cloud fraction above (below) 3 km in four seasons of a year. Normalized by total in-cloud samples, ice and mixed phase occurrence frequencies are underestimated and liquid phase frequency is overestimated by the model at cloud fractions above 0.6, while at cloud fractions below 0.6 ice phase frequency is overestimated and liquid-containing phase frequency is underestimated by the model. The cloud fraction biases are closely associated with concurrent biases in relative humidity (RH), i.e., high (low) RH biases above (below) 2 km. Frequencies of correctly simulating ice and liquid-containing phase increase when the absolute biases of RH decrease. Cloud fraction biases also show a positive correlation with RH biases. Water vapor mixing ratio biases are the primary contributor to RH biases, and hence, likely a key factor controlling the cloud biases. Finally, this diagnosis of the evident shortfalls of representations of cloud characteristics in CAM6 simulation at McMurdo Station brings new insight in improving the governing model physics therein.

54 ENVIRONMENTAL SCIENCES↗