Search NASA⌕ Search

SEARCH · Search NASA

Results for “Temporal modeling”

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 433 records · Page 24

APRES Prototype Mission Planner System Demonstration

Activity Planning with Resources for the Exploration of Space (APRES) is a mixed-initiative mission planning system for ground operations. APRES has been designed to support multi-spacecraft missions. The APRES Interface is browser-based and includes a plan editor, a timeline plan display, a temporal constraint editor, display of the state and numeric chronicles, and a violation resolution manager. Automation support is supplied by the APRES Service, which includes components that provide the following capabilities:(1) plan simulation, which determines the state and numeric chronicles (values of the model variables over time) and determines when "processes" are triggered and terminated based on world states in the execution trace, (2) violation detection of constraints and flight rules encoded in the domain model, and of the temporal constraints created by the user, (3) violation resolution suggestions as to how to fix the plan's violations via rescheduling. The user controls when and how to utilize the automation support. Demo video included with paper, runtime 8:54 in color with sound.

John L. Bresina↗

Robust Rapid Cellular Metabolite Sensing Using Benchtop NMR and SABRE-Hyperpolarized [1- 13 C]Pyruvate

Hyperpolarized NMR has emerged as a powerful analytical technique to significantly enhance targeted NMR signals, improving the sensitivity for investigations of unique chemical and biological dynamics. Here, we demonstrate the use of a hyperpolarization strategy based on Signal Amplification By Reversible Exchange (SABRE) to generate highly reproducible doses of a hyperpolarized [1- 13 C]pyruvate probe for benchtop characterization of yeast metabolism. This method allows rapid, scalable, and benchtop preparation of biocompatible hyperpolarized solutions suitable for live-cell experiments. We show that this production can be dove-tailed into a modular, compact workflow to characterize real-time metabolism in cell cultures, using Saccharomyces cerevisiae (Baker’s yeast) as a model organism. With high temporal resolution, we show that this method can resolve the conversion of hyperpolarized [1- 13 C]pyruvate into oxidative decarboxylation products CO 2 and bicarbonate. This conversion exhibits sustained and detectable metabolic activity for over 300 s after introduction of the agent to the cells. We model the metabolite kinetics to show decarboxylation activity and derive estimates of the pH over time from the CO 2 and bicarbonate (carbonic acid buffer system) equilibrium to probe changes in the cellular environment during active metabolism. These results highlight the utility of benchtop SABRE-hyperpolarized [1- 13 C]pyruvate as a scalable, specific probe for metabolic phenotyping of living cells using compact, low-cost instrumentation well-suited for future high-throughput applications across microbial engineering, drug response profiling, and dynamic metabolic screening.

fungi↗

Deep Learning Scene Classification Experiments in Automatic Detection of Slums on Planetscope Imagery

Population growth is increasingly happening in slum settlements of the large urban centers in the Global South. The term "slum" encompasses a wide range of communities, located mostly in underserved areas, and often exhibiting distinct structural and functional informalities with a relatively high concentration of marginalized populations. To address the issues confronting slums for effective planning and development, including the realistic estimation of the resident population, identifying them accurately is fundamental. Given the disagreements over a universal definition, diverse characteristic features, and socio-political limitations, global detection of slums is a veritable challenge. In this paper, we present experiments in slum detection using a scene classification algorithm and 3-meter spatial resolution satellite imagery. We train and evaluate the model for slum detection in Mumbai, India for the year 2023 and test the temporal generalization of the trained model on Mumbai in 2020 and 2018. In addition, we explore the pathways toward geographic generalization to Kolkata and Delhi (India). We discuss several limitations in the workflow and model, situate our findings in the existing literature, and suggest improvements and alternatives. With this, we establish baseline methods and experiments as a first step towards developing an image-based global slum detection framework and algorithm. This work adds to the community discussion on methods, data challenges, and open questions related to the detection of slums globally. With this research, we hope to improve our understanding of human settlements, especially in critical areas, improve population estimates, and help measure progress towards the sustainable development goals.

Arndt, Jacob↗

Quantitative representation of nonrepetitive temporal behavior

Analytical representations suitable to analyze the nonrepetitive pulse-like temporal behavior of physical quantities are derived. The representation utilizes solutions of a linear model equation in which the temporal variation is subject to time-dependent driving and dissipative forces. The property of solutions is described, and it is shown that such representations can provide a basis for quantitative comparisons of behaviors and a basis for physically meaningful interpretations of the results. Observations of solar flares in the soft X-ray flux have been analyzed with this method.

Nakagawa, Y.↗

Performance analysis and comparison of data-driven models for predicting indoor temperature in multi-zone commercial buildings

Building thermal models, which characterize the properties of a building’s envelope and thermal mass, are essential for accurate indoor temperature and cooling/heating demand prediction. Because of their flexibility and ease of use, data-driven models are increasingly used. Here, this study compared and analyzed the performance of gray-box (resistance-capacitance) and black-box (recurrent neural network) models for predicting indoor air temperature in a real multi-zone commercial building. The developed resistance-capacitance model served as a benchmark model for which full sets of temporal data and building information were used as inputs. The recurrent neural network models were trained and tested assuming various available types and amounts of temporal data and known building physical information to investigate the effects of data and information availability. Feature importance analysis was conducted to select the key variables for different prediction targets under different scenarios. This research provides guidance in selecting an appropriate building thermal response modeling method based on the measured data availability, building physical information, and application.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Impacts of Light Use Efficiency and fPAR Parameterization on Gross Primary Production Modeling

This study examines the impact of parameterization of two variables, light use efficiency (LUE) and the fraction of absorbed photosynthetically active radiation (fPAR or fAPAR), on gross primary production(GPP) modeling. Carbon sequestration by terrestrial plants is a key factor to a comprehensive under-standing of the carbon budget at global scale. In this context, accurate measurements and estimates of GPP will allow us to achieve improved carbon monitoring and to quantitatively assess impacts from cli-mate changes and human activities. Spaceborne remote sensing observations can provide a variety of land surface parameterizations for modeling photosynthetic activities at various spatial and temporal scales. This study utilizes a simple GPP model based on LUE concept and different land surface parameterizations to evaluate the model and monitor GPP. Two maize-soybean rotation fields in Nebraska, USA and the Bartlett Experimental Forest in New Hampshire, USA were selected for study. Tower-based eddy-covariance carbon exchange and PAR measurements were collected from the FLUXNET Synthesis Dataset. For the model parameterization, we utilized different values of LUE and the fPAR derived from various algorithms. We adapted the approach and parameters from the MODIS MOD17 Biome Properties Look-Up Table (BPLUT) to derive LUE. We also used a site-specific analytic approach with tower-based Net Ecosystem Exchange (NEE) and PAR to estimate maximum potential LUE (LUEmax) to derive LUE. For the fPAR parameter, the MODIS MOD15A2 fPAR product was used. We also utilized fAPAR chl, a parameter accounting for the fAPAR linked to the chlorophyll-containing canopy fraction. fAPAR chl was obtained by inversion of a radiative transfer model, which used the MODIS-based reflectances in bands 1-7 produced by Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm. fAPAR chl exhibited seasonal dynamics more similar with the flux tower based GPP than MOD15A2 fPAR, especially in the spring and fall at the agricultural sites. When using the MODIS MOD17-based parameters to estimate LUE, fAPAR chl generated better agreements with GPP (r2= 0.79-0.91) than MOD15A2 fPAR (r2= 0.57-0.84).However, underestimations of GPP were also observed, especially for the crop fields. When applying the site-specific LUE max value to estimate in situ LUE, the magnitude of estimated GPP was closer to in situ GPP; this method produced a slight overestimation for the MOD15A2 fPAR at the Bartlett forest. This study highlights the importance of accurate land surface parameterizations to achieve reliable carbon monitoring capabilities from remote sensing information.

grossprimary↗

History of Thermally Processed Solids in the Protoplanetary Disk: Reconciling Theoretical Models and Meteoritical Evidence

In this talk we assess theoretical models of the radial, temporal, and thermal evolution of nebula solids, and their ultimate accretion into planetesimals such as we see today, using meteorite evidence as a guide. Each class of chondrites contains a characteristic suite of chondrules and CAIs that may have formed over a period of several Myr during which planetesimals were accreting in the disk. Details of the various models for transient melting of chondrules and igneous CAIs will be left to others. However, high-temperature processes of different kinds evaporation, alteration, etc did affect these constituents and their environment over this time span. Here we describe evolutionary scenarios consistent with a large time gap between CAI and chondrule formation and the presence of distinctive suites of chondrules and CAIs in each chondrite class. Particle-gas dynamical processes transport particles of all relevant sizes (microns to many meters) within the nebula and affect their evolution in a variety of important ways. Turbulent radial diffusion spreads particles radially down their concentration gradients - as one example, it can prevent CAIs from being lost into the sun on several Myr timescales [1]. Vertical diffusion spreads the dense midplane particle layer, determining its volume density, which in turn affects the particle growth rate and even the dominant growth process [2-4]. Turbulent concentration selects aerodynamically sorted particles for orders-of-magnitude density enhancement, and is applicable to porous, fluffy particles of appropriate size as well as to solid chondrules [5]. Inward radial drift under gas drag brings a surprisingly large amount of material to regions where it evaporates; these evaporation fronts cause significant chemical modification of the nebula gas over a wide range of radii [6]. Radial transport by stellar winds can be important for small particles [7].

Cuzzi, J. N.↗

Systems and methods for modeling water quality

A system, method, device and computer-readable medium for creating an ensemble model of water quality. The ensemble model is generated by determining a set of optimal component models for spectral regions of a body of water, and combining the optimal models. The optimal models can be based on remote sensing data, including satellite imagery. A K-fold partition approach or a global approach can be used to determine the optimal component models, and the optimal component models can be combined through spectral space partition rules to generate an ensemble model of water quality. The ensemble model not only has improved water quality prediction ability, but also has strong spatial and temporal extensibility. The spatial and temporal extensibility of the ensemble model is fundamentally important and desirable for long-term and large-scale remote sensing monitoring and assessment of water quality.

Liu, Hongxing↗

Integrating data types to estimate spatial patterns of avian migration across the Western Hemisphere

For many avian species, spatial migration patterns remain largely undescribed, especially across hemispheric extents. Recent advancements in tracking technologies and high-resolution species distribution models (i.e., eBird Status and Trends products) provide new insights into migratory bird movements and offer a promising opportunity for integrating independent data sources to describe avian migration. Here, we present a three-stage modeling framework for estimating spatial patterns of avian migration. First, we integrate tracking and band re-encounter data to quantify migratory connectivity, defined as the relative proportions of individuals migrating between breeding and nonbreeding regions. Next, we use estimated connectivity proportions along with eBird occurrence probabilities to produce probabilistic least-cost path (LCP) indices. In a final step, we use generalized additive mixed models (GAMMs) both to evaluate the ability of LCP indices to accurately predict (i.e., as a covariate) observed locations derived from tracking and band re-encounter data sets versus pseudo-absence locations during migratory periods and to create a fully integrated (i.e., eBird occurrence, LCP, and tracking/band re-encounter data) spatial prediction index for mapping species-specific seasonal migrations. To illustrate this approach, we apply this framework to describe seasonal migrations of 12 bird species across the Western Hemisphere during pre- and postbreeding migratory periods (i.e., spring and fall, respectively). We found that including LCP indices with eBird occurrence in GAMMs generally improved the ability to accurately predict observed migratory locations compared to models with eBird occurrence alone. Using three performance metrics, the eBird + LCP model demonstrated equivalent or superior fit relative to the eBird-only model for 22 of 24 species–season GAMMs. In particular, the integrated index filled in spatial gaps for species with over-water movements and those that migrated over land where there were few eBird sightings and, thus, low predictive ability of eBird occurrence probabilities (e.g., Amazonian rainforest in South America). This methodology of combining individual-based seasonal movement data with temporally dynamic species distribution models provides a comprehensive approach to integrating multiple data types to describe broad-scale spatial patterns of animal movement. Further development and customization of this approach will continue to advance knowledge about the full annual cycle and conservation of migratory birds.

59 BASIC BIOLOGICAL SCIENCES↗

Standing waves observed at the dayside magnetopause

Energetic (greater than 24 keV) ion distributions are used in remotely sounding an ion trapping boundary at the dayside magnetopause for two inbound crossings of ISEE-2 on days 347 and 312, 1977, at 0830 and 1100 LT. Applying single particle trajectory models to the observed temporal relationships between trapping boundary distance and orientation parameters, it is concluded that longitudinal standing waves occur on the boundary, having typical wavelengths of 1500-3000 km, amplitudes of 500 km, and periods of 200-400 sec. The standing wave nodes move along the node axis with velocities ranging from 1.5 to 6 km/sec. The data are found to be inconsistent with models of traveling waves moving tailward along the boundary.

Fahnenstiel, S. C.↗

Towards an A-train Aerosol Assimilation System: Assimilation of MODIS Aerosol Optical Thickness Retrievals into a Global Aerosol Transport and Radiation Model

To understand the role of atmospheric aerosols in Earth's climate system we ideally require a global, three-dimensional, and temporally continuous system of measurements. The view of aerosols provided by a single polar-orbiting satellite instrument (e.g., MODIS), while providing the greatest spatial coverage of any measurement system available, can offer only a single relevant observation per day of a particular location. Furthermore, the quantitative interpretation of those observations is subject to considerable uncertainty because of we do not know the aerosol optical properties or vertical profiles well enough. Ground-based and in situ measurements can provide missing detail about aerosol vertical and optical properties, but lack the spatial (and sometimes temporal) coverage desired. Chemical transport models treating the aerosol lifecycle complement this picture by providing three-dimensional, global, and continuous aerosol fields, but suffer from poor constraints on their initial conditions (e.g., the aerosol source, sink, and transformation mechanisms). Here we present a global, three- dimensional aerosol transport model and assimilation system. The model assimilates the MODIS aerosol optical thickness product to constrain its aerosol distributions and burdens. We provide an analysis of the model aerosol distributions as constrained by ground-based sun photometer measurements and aerosol products from other satellite measurements.

Colarco, P.↗

Improved Satellite Retrievals of NO2 and SO2 over the Canadian Oil Sands and Comparisons with Surface Measurements

Satellite remote sensing is increasingly being used to monitor air quality over localized sources such as the Canadian oil sands. Following an initial study, significantly low biases have been identified in current NO2 and SO2 retrieval products from the Ozone Monitoring Instrument (OMI) satellite sensor over this location resulting from a combination of its rapid development and small spatial scale. Air mass factors (AMFs) used to convert line-of-sight "slant" columns to vertical columns were re-calculated for this region based on updated and higher resolution input information including absorber profiles from a regional-scale (15 km × 15 km resolution) air quality model, higher spatial and temporal resolution surface reflectivity, and an improved treatment of snow. The overall impact of these new Environment Canada (EC) AMFs led to substantial increases in the peak NO2 and SO2 average vertical column density (VCD), occurring over an area of intensive surface mining, by factors of 2 and 1.4, respectively, relative to estimates made with previous AMFs. Comparisons are made with long-term averages of NO2 and SO2 (2005-2011) from in situ surface monitors by using the air quality model to map the OMI VCDs to surface concentrations. This new OMI-EC product is able to capture the spatial distribution of the in situ instruments (slopes of 0.65 to 1.0, correlation coefficients of greater than 0.9). The concentration absolute values from surface network observations were in reasonable agreement, with OMI-EC NO2 and SO2 biased low by roughly 30%. Several complications were addressed including correction for the interference effect in the surface NO2 instruments and smoothing and clear-sky biases in the OMI measurements. Overall these results highlight the importance of using input information that accounts for the spatial and temporal variability of the location of interest when performing retrievals.

surface concentration↗

Spatio-Temporal Surrogates for Interaction of a Jet with High Explosives: Part II - Clustering Extremely High-Dimensional Grid-Based Data

Building an accurate surrogate model for the spatio-temporal outputs of a computer simulation is a challenging task. A simple approach to improve the accuracy of the surrogate is to cluster the outputs based on similarity and build a separate surrogate model for each cluster. This clustering is relatively straightforward when the output at each time step is of moderate size. However, when the spatial domain is represented by a large number of grid points, numbering in the millions, the clustering of the data becomes more challenging. In this report, we consider output data from simulations of a jet interacting with high explosives. These data are available on spatial domains of different sizes, at grid points that vary in their spatial coordinates, and in a format that distributes the output across multiple files at each time step of the simulation. We first describe how we bring these data into a consistent format prior to clustering. Borrowing the idea of random projections from data mining, we reduce the dimension of our data by a factor of thousand, making it possible to use the iterative k-means method for clustering. We show how we can use the randomness of both the random projections, and the choice of initial centroids in k-means clustering, to determine the number of clusters in our data set. Our approach makes clustering of extremely high dimensional data tractable, generating meaningful cluster assignments for our problem, despite the approximation introduced in the random projections.

97 MATHEMATICS AND COMPUTING↗

Gas Generator Feedline Orifice Sizing Methodology: Effects of Unsteadiness and Non-Axisymmetric Flow

Engine LH2 and LO2 gas generator feed assemblies were modeled with computational fluid dynamics (CFD) methods at 100% rated power level, using on-center square- and round-edge orifices. The purpose of the orifices is to regulate the flow of fuel and oxidizer to the gas generator, enabling optimal power supply to the turbine and pump assemblies. The unsteady Reynolds-Averaged Navier-Stokes equations were solved on unstructured grids at second-order spatial and temporal accuracy. The LO2 model was validated against published experimental data and semi-empirical relationships for thin-plate orifices over a range of Reynolds numbers. Predictions for the LO2 square- and round-edge orifices precisely match experiment and semi-empirical formulas, despite complex feedline geometry whereby a portion of the flow from the engine main feedlines travels at a right-angle through a smaller-diameter pipe containing the orifice. Predictions for LH2 square- and round-edge orifice designs match experiment and semi-empirical formulas to varying degrees depending on the semi-empirical formula being evaluated. LO2 mass flow rate through the square-edge orifice is predicted to be 25 percent less than the flow rate budgeted in the original engine balance, which was subsequently modified. LH2 mass flow rate through the square-edge orifice is predicted to be 5 percent greater than the flow rate budgeted in the engine balance. Since CFD predictions for LO2 and LH2 square-edge orifice pressure loss coefficients, K, both agree with published data, the equation for K has been used to define a procedure for orifice sizing.

Rothermel, Jeffry↗

Ectomycorrhizal effects on decomposition are highly dependent on fungal traits, climate, and litter properties: A model-based assessment. Dataset.

To simulate the effects of mycorrhizal fungi on soil organic matter cycling, we incorporated mycorrhizal processes into the Carbon, Organisms, Rhizosphere, and Protection in the Soil Environment (CORPSE) model to develop a new soil model Myco-CORPSE. The new model was calibrated and evaluated against soil measurements taken at temperate forests in New Hampshire (NH) and Georgia (GA). A series of scenario analysis were also conducted to explore the conditions under which ectomycorrhizal (ECM) N acquisition processes can induce different soil C accumulation in ECM systems compared to arbuscular (AM) systems.In this data package, we included:-The Python codes of the standard Myco-CORPSE model we developed: "Standard Myco_CORPSE python codes.zip". The main program is the "gradient_sim.py" which calculates the bulk soil microbes and CN content along a user defined gradient of clay, soil temperature, soil moisture and mycorrhizal dominance, and relies on two subprograms "CORPSE_deriv.py" and "CORPSE_integrate.py". "CORPSE_deriv.py" calculated the changes in all simulated soil stock within every time step and "CORPSE_integrate.py" integrate the changes in all simulated soil stock within simulated time period. The program "Plot.py" is used to plot the major outputs produced by the main program "gradient_sim.py".-The modified Python codes of Myco-CORPSE models with site-level environmental inputs (in NH and GA) used to conduct simulations in NH and GA sites: "NH_GA model simulations.zip". -The Python codes used to evaluate the Myco-CORPSE simulation outputs in NH and GA sites against site-level measurements: "Plot NH_GA simulation against measurements.zip". It includes both the evaluation Python code, the model outputs on NH and GA sites, and the measured soil properties in both sites.-The modified Python codes of Myco-CORPSE models "Scenario analysis_model simulations.zip" that is used to conduct scenario analysis of how different litter properties, mycorrhizal fungal traits, climate, and seasonal variation in temperature and vegetation phenology impact the mycorrhizal effects on soil CN properties. The sub file folder "Scenario analysis_litter traits" contains the codes for scenario analysis of different litter properties; The sub file folder "Scenario analysis_ECM types" contains the codes for scenario analysis of different ECM fungal traits; The sub file folder "Scenario analysis_climate&seasonality" contains the codes for scenario analysis of different climate and seasonalities;-"Scenario analysis_model results and plotting codes.zip" contains all the output files from the the scenario analysis of Myco-CORPSE model as described above and the plotting codes used to the generate the heatmaps and scatterplots shown in the manuscript "Ectomycorrhizal effects on decomposition are highly dependent on fungal traits, climate, and litter properties: A model-based assessment"The majority of the model outputs did not have specific geographic information or temporal coverage because the analysis we conducted are mainly hypothetical model simulations. We only provided geographic description, coordinates and temporal coverage for those soil measurements which we used for model evaluations (included in the "Plot NH_GA simulation against measurements.zip").

54 ENVIRONMENTAL SCIENCES↗

HCl in rocket exhaust clouds - Atmospheric dispersion, acid aerosol characteristics, and acid rain deposition

Both measurements and model calculations of the temporal dispersion of peak HCl (g + aq) concentration in Titan III exhaust clouds are found to be well characterized by one-term power-law decay expressions. The respective coefficients and decay exponents, however, are found to vary widely with meteorology. The HCl (g), HCl (g + aq), dewpoint, and temperature-pressure-altitude data for Titan III exhaust clouds are consistent with accurately calculated HCl/H2O vapor-liquid compositions for a model quasi-equilibrated flat surface aqueous aerosol. Some cloud evolution characteristics are also defined. Rapid and extensive condensation of aqueous acid clearly occurs during the first three min of cloud rise. Condensation is found to be intensified by the initial entrainment of relatively moist ambient air from lower levels, that is, from levels below eventual cloud stabilization. It is pointed out that if subsequent dilution air at stabilization altitude is significantly drier, a state of maximum condensation soon occurs, followed by an aerosol evaporation phase.

Pellett, G. L.↗

Machine Learning Surrogates of a Fuel Matrix Degradation Process Model for Performance Assessment of a Nuclear Waste Repository

Spent nuclear fuel repository simulations are currently not able to incorporate detailed fuel matrix degradation (FMD) process models due to their computational cost, especially when large numbers of waste packages breach. The current paper uses machine learning to develop artificial neural network and k-nearest neighbor regression surrogate models that approximate the detailed FMD process model while being computationally much faster to evaluate. Further, using fuel cask temperature, dose rate, and the environmental concentrations of CO 3 2- , O 2 , Fe 2+ , and H 2 as inputs, these surrogates show good agreement with the FMD process model predictions of the UO 2 degradation rate for conditions within the range of the training data. A demonstration in a full-scale shale repository reference case simulation shows that the incorporation of the surrogate models captures local and temporal environmental effects on fuel degradation rates while retaining good computational efficiency.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Deeper snow increases the net soil organic carbon accrual rate in moist acidic tussock tundra: 210Pb evidence from Arctic Alaska

The net change in the carbon inventory of arctic tundra remains uncertain as global warming leads to shifts in arctic water and carbon cycles. To better understand the response of arctic tundra carbon to changes in winter precipitation amount, we investigated soil depth profiles of carbon concentration and radionuclide activities ( 7 Be, 137 Cs, 210 Pb, and 241 Am) in the active layer of a twenty-two-year winter snow depth manipulation experiment in moist acidic tussock tundra at Toolik Lake, Alaska. Depth correlations of cumulative carbon dry mass (g cm -2 ) vs. unsupported 210 Pb activity (mBq g -1 ) were examined using a modified constant rate of supply (CRS) model. Results were best fit by two-slope CRS models indicating an apparent step temporal increase in the accumulation rate of soil organic carbon. Most of the best-fit model chronologies indicated that the increase in carbon accumulation rate apparently began and persisted after snow fence construction in 1994. The inhomogeneous nature of permafrost soils and their relatively low net carbon accumulation rates make it challenging to establish robust chronologic records. Nonetheless, the data obtained in this study support a decadal-scale increase in net soil organic carbon accumulation rate in the active layer of arctic moist acidic tussock tundra under conditions of increased winter precipitation.

54 ENVIRONMENTAL SCIENCES↗