Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Environmental Models”

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 1,135 records · Page 63

Plateau to River Model Predictive Simulations for All Ensemble Realizations to Support Modeling Work in Fiscal Year 2025

The purpose of this environmental calculation file (ECF) is to document predictions of flow and hydraulic head on the Central Plateau of the Hanford Site using the Plateau-to-River (P2R) Model (CP-57037, Model Package Report for the Plateau-to-River Model: Version 9.1). This calculation documents the simulation of the groundwater for the parent model domain of the P2R Model as a basis for use in other applications of the P2R Model. This application is unique from the standpoint that it will simulate all ensemble member models of the P2R Model whereas other applications may only utilize specific ensemble members. These simulations provide results that can be used in the process of selecting an appropriate subset of ensemble members for other applications.

54 ENVIRONMENTAL SCIENCES↗

CFD-trained ANN Model for Approximating Near-occupant Condition in Real-time Simulations

The main drawback of Computational Fluid Dynamics (CFD) simulations has been the time and resource consuming nature which is not suitable for real-time applications. In this work, we first generated numerous CFD models of a given indoor space to obtain airspeed, temperature, and mean radiant temperature near an occupant as training data. Several artificial neural networks (ANN) models were trained using this CFD simulated data to approximate near real-time environmental conditions for a given occupant. This trained ANN model approach is a part of a real-time simulation of building operations using a combination of software and real hardware (HVAC equipment) approaches. The preliminary results suggest that the CFD- generated training data and the trained ANN model can accurately approximate such conditions in a real-time application, a method that has great potential in building simulation and building digital twin areas of research.

Zhang, Yun↗

Computational Fluid Dynamics Uncertainty Analysis for Payload Fairing Spacecraft Environmental Control Systems

Spacecraft thermal protection systems are at risk of being damaged due to airflow produced from Environmental Control Systems. There are inherent uncertainties and errors associated with using Computational Fluid Dynamics to predict the airflow field around a spacecraft from the Environmental Control System. This proposal describes an approach to validate the uncertainty in using Computational Fluid Dynamics to predict airflow speeds around an encapsulated spacecraft. The research described here is absolutely cutting edge. Quantifying the uncertainty in analytical predictions is imperative to the success of any simulation-based product. The method could provide an alternative to traditional"validation by test only'' mentality. This method could be extended to other disciplines and has potential to provide uncertainty for any numerical simulation, thus lowering the cost of performing these verifications while increasing the confidence in those predictions. Spacecraft requirements can include a maximum airflow speed to protect delicate instruments during ground processing. Computationaf Fluid Dynamics can be used to veritY these requirements; however, the model must be validated by test data. The proposed research project includes the following three objectives and methods. Objective one is develop, model, and perform a Computational Fluid Dynamics analysis of three (3) generic, non-proprietary, environmental control systems and spacecraft configurations. Several commercially available solvers have the capability to model the turbulent, highly three-dimensional, incompressible flow regime. The proposed method uses FLUENT and OPEN FOAM. Objective two is to perform an uncertainty analysis of the Computational Fluid . . . Dynamics model using the methodology found in "Comprehensive Approach to Verification and Validation of Computational Fluid Dynamics Simulations". This method requires three separate grids and solutions, which quantify the error bars around Computational Fluid Dynamics predictions. The method accounts for all uncertainty terms from both numerical and input variables. Objective three is to compile a table of uncertainty parameters that could be used to estimate the error in a Computational Fluid Dynamics model of the Environmental Control System /spacecraft system. Previous studies have looked at the uncertainty in a Computational Fluid Dynamics model for a single output variable at a single point, for example the re-attachment length of a backward facing step. To date, the author is the only person to look at the uncertainty in the entire computational domain. For the flow regime being analyzed (turbulent, threedimensional, incompressible), the error at a single point can propagate into the solution both via flow physics and numerical methods. Calculating the uncertainty in using Computational Fluid Dynamics to accurately predict airflow speeds around encapsulated spacecraft in is imperative to the success of future missions.

Groves, Curtis E.↗

A Predictive Model for Survival of Escherichia coli O157:H7 and Generic E. coli in Soil Amended with Untreated Animal Manure

Abstract This study aimed at developing a predictive model that captures the influences of a variety of agricultural and environmental variables and is able to predict the concentrations of enteric bacteria in soil amended with untreated Biological Soil Amendments of Animal Origin (BSAAO) under dynamic conditions. We developed and validated a Random Forest model using data from a longitudinal field study conducted in mid‐Atlantic United States investigating the survival of Escherichia coli O157:H7 and generic E. coli in soils amended with untreated dairy manure, horse manure, or poultry litter. Amendment type, days of rain since the previous sampling day, and soil moisture content were identified as the most influential agricultural and environmental variables impacting concentrations of viable E. coli O157:H7 and generic E. coli recovered from amended soils. Our model results also indicated that E. coli O157:H7 and generic E. coli declined at similar rates in amended soils under dynamic field conditions.The Random Forest model accurately predicted changes in viable E. coli concentrations over time under different agricultural and environmental conditions. Our model also accurately characterized the variability of E. coli concentration in amended soil over time by providing upper and lower prediction bound estimates. Cross‐validation results indicated that our model can be potentially generalized to other geographic regions and incorporated into a risk assessment for evaluating the risks associated with application of untreated BSAAO. Our model can be validated for other regions and predictive performance also can be enhanced when data sets from additional geographic regions become available.

Pang, Hao↗

Completion of initial reduced-order model for flammable refrigerant dispersal in residential spaces

Environmental regulations aimed at reducing global warming impacts of HVAC&R refrigerants have resulted in the phase-out of chlorofluorocarbons (CFCs) in 2010 and hydrochlorofluorocarbons (HCFCs) by 2030 in developed countries with additional restrictions on the use of hydrofluorocarbons (HFCs) set to take effect in 2036. Many of the remaining alternative refrigerants that have lower global warming potential (GWP) and that are suitable for use in HVAC&R systems (e.g., propane, difluoromethane) are flammable to some degree. Flammable refrigerants introduce new challenges and hazards to property and personal safety related to potential deflagration during system maintenance or due to leakage of the refrigerant accumulating in the conditioned space. Standards have been developed to set maximum charge limits for flammable refrigerants in HVAC&R systems; however, existing safety standards and codes still restrict their use. The bodies that maintain and update these codes need publicly available, science-based information to enable credible guidelines for setting safe charge limits for different flammable refrigerants in different HVAC&R applications. In 2016, the Alliance for Responsible Atmospheric Policy, the Air-Conditioning Heating and Refrigeration Institute (AHRI), ASHRAE, the U.S. Department of Energy (DOE), and the State of California began efforts to develop such information. As part of this effort, Oak Ridge National Laboratory (ORNL) began the current, ongoing project to examine imposed charge limits for flammable refrigerants and identify reasonable adjustments to these limits when found appropriate. Past tasks under this project have included development of experimentally vetted, computational fluid dynamics (CFD) simulation approaches to study the results of leakage of flammable refrigerants from various HVAC&R systems into different types of commercial and residential spaces. In this report, we discuss recent efforts to develop a predictive model of the flammable volume fraction and accumulated refrigerant mass in a single-room residential space resulting from the leak of a flammable refrigerant from a small room air conditioning (RAC) unit. The eventual goal is development of a model suitable for public release which could simulate a range of scenarios. In discussions with the AHRTI Flammable Refrigerant Subcommittee (FRS) at the beginning of this effort, a total of 9 input parameters were chosen for consideration including room area, room/door opening area, ventilation fan flow rate, unit/leak height, leak area, leak rate, total refrigerant charge, refrigerant molecular weight, and state of the unit fan.

42 ENGINEERING↗

Unified Entrainment and Detrainment Closures for Extended Eddy-Diffusivity Mass-Flux Schemes

The extended Eddy Diffusivity Mass Flux (EDMF) scheme is used as a unied parameterization that simulates subgrid-scale dynamics across a range of regimes.This work presents unied entrainment and detrainment closures that allow this scheme to simulate the dry convective boundary layer, shallow convection and deep convection , while adjacent works present other necessary closures for this unied scheme and test its performance for stable, stratocumulus topped and convective boundary layers. Entrainment and detrainment are modeled as a combination of turbulent and dynamical processes. The turbulent component is taken as downgradient diffusion between subdomains. The dynamical component is based on a ratio of buoyancy and vertical velocity that ensures consistent treatment of the ratio of mass-flux at vanishing area fraction. Dynamical detrainment is controlled by the difference in relative humidity between subdomains when liquid water is present. These closures are rst evaluated with entrainment and detrainment diagnosed from Large Eddy Simulations where decaying tracers are used to identify the different subdomains and mass exchanges between them. When implemented in a single column model (SCM), these closures are tested in simulations of several observationally motivated cases of dry convective boundary layer, shallow convection and deep convection. Results show good agreement between the SCM and Large Eddy Simulations in key properties for climate models and profiles of first, second and third moments. The development of diurnal cycles and the onset of and precipitation (with simplified microphysics) are well predicted by the SCM.

54 ENVIRONMENTAL SCIENCES↗

Combined Effects of Stream Hydrology and Land Use on Basin‐Scale Hyporheic Zone Denitrification in the Columbia River Basin

Abstract Denitrification in the hyporheic zone (HZ) of river corridors is crucial to removing excess nitrogen in rivers from anthropogenic activities. However, previous modeling studies of the effectiveness of river corridors in removing excess nitrogen via denitrification were often limited to the reach‐scale and low‐order stream watersheds. We developed a basin‐scale river corridor model for the Columbia River Basin with random forest models to identify the dominant factors associated with the spatial variation of HZ denitrification. Our modeling results suggest that the combined effects of hydrologic variability in reaches and substrate availability influenced by land use are associated with the spatial variability of modeled HZ denitrification at the basin scale. Hyporheic exchange flux can explain most of spatial variation of denitrification amounts in reaches of different sizes, while among the reaches affected by different land uses, the combination of hyporheic exchange flux and stream dissolved organic carbon (DOC) concentration can explain the denitrification differences. Also, we can generalize that the most influential watershed and channel variables controlling denitrification variation are channel morphology parameters (median grain size (D50), stream slope), climate (annual precipitation and evapotranspiration), and stream DOC‐related parameters (percent of shrub area). The modeling framework in our study can serve as a valuable tool to identify the limiting factors in removing excess nitrogen pollution in large river basins where direct measurement is often infeasible.

54 ENVIRONMENTAL SCIENCES↗

SIMPLE-G: A multiscale framework for integration of economic and biophysical determinants of sustainability

We introduce SIMPLE-G, a Simplified International Model of agricultural Prices, Land use, and the Environment- Gridded version, which is a novel tool for evaluating sustainability policies in a global context while factoring in local heterogeneity in land and water resources and natural ecosystem services. This multi-scale model can provide boundary conditions for local decision makers, as well as capturing feedback from local policies to national and global scales. Additionally, to illustrate its value in environmental analysis, we provide two applications of the model. First, we quantify the local stresses on land and water resources due to global changes in population, income, and productivity. Second, we quantify the global impacts of local policy responses and adaptations to water scarcity.

42 ENGINEERING↗

Renewable Energy Potential Model: Geothermal Supply Curves

The Renewable Energy Potential (reV) model is a geospatial platform for estimating technical potential and developing renewable energy supply curves, initially developed for wind and solar technologies. The model evaluates deployment constraints, considering land use, environmental, and cultural factors, and estimates the distance to existing grid features to connect future plants (Maclaurin et al., 2021). A pressing deficiency in the reV model, however, is representation of geothermal electricity generation technologies. To address this gap, we developed a novel geothermal generation module for reV that allows for representation and analysis at the same level of detail as other renewable technologies. The included paper describes our process for evaluating data sources for the modeling, and presents five preliminary reV geothermal results. More specifically, we present two sets of resource data that represent upper and lower bounds for geothermal potential. We then present several sensitivity runs using the upper bound resource data; the results are encouraging that levelized cost of electricity (LCOE) can be reduced by optimizing the location and estimated capacity of the spatially diverse geothermal resource while considering the distance to existing grid infrastructure. Our preliminary supply curves and levelized cost of electricity (LCOE) results provided here should be considered with care due to the high uncertainty in geothermal resource potential data. We present median LCOE values for the conterminous U.S. for three scenarios: two hydrothermal (3.5km depth, USGS heat flow & SMU temperatures respectively) and one EGS (4.5km depth, SMU temperatures). The capital and operating costs for each respective technology are modeled. We also compare results using two different resource data sources.

15 GEOTHERMAL ENERGY↗

Biological indicators for monitoring water quality of MTF canals system

Biological models, diversity indexes, were developed to predict environmental effects of NASA's Mississippi test facility (MTF) chemical operations on canal systems in the area. To predict the effects on local streams, a physical model of unpolluted streams was established. The model is fed by artesian well water free of background levels of pollutants. The species diversity and biota composition of unpolluted MTF stream was determined; resulting information will be used to form baseline data for future comparisons. Biological modeling was accomplished by adding controlled quantities or kinds of chemical pollutants and evaluating the effects of these chemicals on the biological life of the stream.

Sethi, S. L.↗

Performance Considerations for the SIMPL Single Photon, Polarimetric, Two-Color Laser Altimeter as Applied to Measurements of Forest Canopy Structure and Composition

The Slope Imaging Multi-polarization Photon-counting Lidar (SIMPL) is a multi-beam, micropulse airborne laser altimeter that acquires active and passive polarimetric optical remote sensing measurements at visible and near-infrared wavelengths. SIMPL was developed to demonstrate advanced measurement approaches of potential benefit for improved, more efficient spaceflight laser altimeter missions. SIMPL data have been acquired for wide diversity of forest types in the summers of 2010 and 2011 in order to assess the potential of its novel capabilities for characterization of vegetation structure and composition. On each of its four beams SIMPL provides highly-resolved measurements of forest canopy structure by detecting single-photons with 15 cm ranging precision using a narrow-beam system operating at a laser repetition rate of 11 kHz. Associated with that ranging data SIMPL provides eight amplitude parameters per beam unlike the single amplitude provided by typical laser altimeters. Those eight parameters are received energy that is parallel and perpendicular to that of the plane-polarized transmit pulse at 532 nm (green) and 1064 nm (near IR), for both the active laser backscatter retro-reflectance and the passive solar bi-directional reflectance. This poster presentation will cover the instrument architecture and highlight the performance of the SIMPL instrument with examples taken from measurements for several sites with distinct canopy structures and compositions. Specific performance areas such as probability of detection, after pulsing, and dead time, will be highlighted and addressed, along with examples of their impact on the measurements and how they limit the ability to accurately model and recover the canopy properties. To assess the sensitivity of SIMPL's measurements to canopy properties an instrument model has been implemented in the FLIGHT radiative transfer code, based on Monte Carlo simulation of photon transport. SIMPL data collected in 2010 over the Smithsonian Environmental Research Center, MD are currently being modelled and compared to other remote sensing and in situ data sets. Results on the adaptation of FLIGHT to model micropulse, single'photon ranging measurements are presented elsewhere at this conference. NASA's ICESat-2 spaceflight mission, scheduled for launch in 2016, will utilize a multi-beam, micropulse, single-photon ranging measurement approach (although non-polarimetric and only at 532 nm). Insights gained from the analysis and modelling of SIMPL data will help guide preparations for that mission, including development of calibration/validation plans and algorithms for the estimation of forest biophysical parameters.

Dabney, Philip W.↗

Development of a Geothermal Module in reV: Quantifying the Geothermal Potential while Accounting for the Geospatial Intersection of the Grid Infrastructure and Land Use Characteristics

The Renewable Energy Potential (reV) model is a geospatial platform for estimating technical potential and developing renewable energy supply curves, initially developed for wind and solar technologies. The model evaluates deployment constraints, considering land use, environmental, and cultural factors, and estimates the distance to existing grid features to connect future plants (Maclaurin et al., 2021). A pressing deficiency in the reV model, however, is representation of geothermal electricity generation technologies. To address this gap, we developed a novel geothermal generation module for reV that allows for representation and analysis at the same level of detail as other renewable technologies. This paper describes our process for evaluating data sources for the modeling, and presents five preliminary reV geothermal results. More specifically, we present two sets of resource data that represent upper and lower bounds for geothermal potential. We then present several sensitivity runs using the upper bound resource data; the results are encouraging that levelized cost of electricity (LCOE) can be reduced by optimizing the location and estimated capacity of the spatially diverse geothermal resource while considering the distance to existing grid infrastructure. Our preliminary supply curves and levelized cost of electricity (LCOE) results should be considered with care due to the highly uncertainty in geothermal resource potential data. We present median LCOE values for the conterminous U.S. for five scenarios: four hydrothermal (3.5km depth) and one EGS (4.5km depth). The capital and operating costs for each respective technology are modeled. We also compare results using two different resource data sources.

exclusions↗

Development of a Geothermal Module in reV: Quantifying the Geothermal Potential While Accounting for the Geospatial Intersection of the Grid Infrastructure and Land Use Characteristics: Preprint

The Renewable Energy Potential (reV) model is a geospatial platform for estimating technical potential and developing renewable energy supply curves, initially developed for wind and solar technologies. The model evaluates deployment constraints, considering land use, environmental, and cultural factors, and estimates the distance to existing grid features to connect future plants (Maclaurin et al., 2021). A pressing deficiency in the reV model, however, is representation of geothermal electricity generation technologies. To address this gap, we developed a novel geothermal generation module for reV that allows for representation and analysis at the same level of detail as other renewable technologies. This paper describes our process for evaluating data sources for the modeling, and presents five preliminary reV geothermal results. More specifically, we present two sets of resource data that represent upper and lower bounds for geothermal potential. We then present several sensitivity runs using the upper bound resource data; the results are encouraging that levelized cost of electricity (LCOE) can be reduced by optimizing the location and estimated capacity of the spatially diverse geothermal resource while considering the distance to existing grid infrastructure. Our preliminary supply curves and levelized cost of electricity (LCOE) results should be considered with care due to the highly uncertainty in geothermal resource potential data. We present median LCOE values for the conterminous U.S. for five scenarios: four hydrothermal (3.5km depth) and one EGS (4.5km depth). The capital and operating costs for each respective technology are modeled. We also compare results using two different resource data sources.

exclusions↗

A Multiyear Gridded Data Ensemble of Surface Biogenic Carbon Fluxes for North America: Evaluation and Analysis of Results

Accurate and fine-scale estimates of biogenic carbon fluxes are critical for measuring and monitoring the biosphere's responses and feedback to the climate system. Currently available data products from flux towers and model-intercomparison projects struggle to adequately represent spatiotemporal dynamics of surface biogenic carbon fluxes, and to quantify their uncertainties, which also are crucial to atmospheric inversion systems. To address these gaps, we introduce a new perturbed-parameter model ensemble with the CASA model to estimate surface biogenic carbon fluxes at monthly and 3-hourly scales for North America at ~500-m and 5-km resolutions. We first use the Extended Fourier Amplitude Sensitivity Testing to choose the three most sensitive parameters to be perturbed, maximum light use efficiency (E max ), optimal temperature of photosynthesis (T opt ), and temperature response of respiration (Q 10 ). The initial range for each parameter is broadly sampled for the L1 ensemble, but then we pruned E max with site-level primary productivity to derive an L2 ensemble with narrower uncertainty ranges. Ensembles are strongly correlated with site-level results at both monthly and 3-hourly scales, and the spread across L1/L2 ensemble members encompasses the range of AmeriFlux observations. Monthly variability in the L2 ensemble mean is 85% of the observed variability. The L2 ensemble outperforms diverse data products with the highest Taylor skill scores at diurnal to annual scales. Finally, the ensemble's seasonality agrees well with other models for most biome types and in high and middle latitudes, but inconsistencies are found in subtropical and tropical ecoregions and for annual totals over North America.

54 ENVIRONMENTAL SCIENCES↗

Late day measurement of excised branches results in uncertainty in the estimation of two stomatal parameters derived from response curves in Populus deltoides Bartr. x Populus nigra L.

Many terrestrial biosphere models depend on an understanding of the relationship between stomatal conductance and photosynthesis. However, unlike the measurement of photosynthetic parameters, such as the maximum carboxylation capacity, where standard methods (e.g., CO 2 response or ACi curves) are widely accepted, a consensus method for empirically measuring parameters representing stomatal response has not yet emerged. Most models of stomatal response to environment represent stomatal conductance as being bounded by a lower intercept parameter (g 0 ), and linearly scaled based on a multivariate term described by the stomatal slope parameter (g 1 ). Here we employ the widely used Unified Stomatal Optimization model, to test whether g 1 and g 0 parameters are impacted by the choice of measurement method, either on an intact branch, or a cut branch segment stored in water. We measured paired stomatal response curves on intact and excised branches of a hybrid poplar clone (Populus deltoides Bartr. x Populus nigra L. OP367), measured twice over a diurnal period. We found that predawn branch excision did not significantly affect measured g 0 and g 1 when measured within 4 h of excision. Measurement in the afternoon resulted in significantly higher values of g 1 and lower values of g0, with values changing by 55% and 56% respectively. Excision combined with afternoon measurement resulted in a marked effect on parameter estimates, with g1 increasing 89% from morning to afternoon and a 25% lower g 1 for cut branches than those measured in situ. We also show that in hybrid poplar the differences in parameter estimates obtained from plants measured under different conditions can directly impact models of canopy function, reducing modeled transpiration by 18% over a simulated 12.5-hour period. While these findings are only for a single isohydric woody species, our findings suggest that stomatal optimality parameters may not remain constant throughout the day.

54 ENVIRONMENTAL SCIENCES↗

Environmental impact of national and subnational carbon policies in China based on a multi-regional dynamic CGE model

Carbon policies are important not only for mitigating global climate change, but also for controlling local and regional air pollution. The large regional disparities in economic development and air pollution across China calls for regionally differentiated policies. Previous studies have not systematically investigated the environmental impacts of regional carbon policies in China covering different spatial scales. This study constructs a multi-regional dynamic computable general quilibrium (CGE) model covering 30 provinces of China to assess the impacts of national and subnational carbon polices on CO2 emissions and co-emitted air pollutants from 2020 to 2050. We consider one national carbon policy which aims to achieve China’s 2030 national CO2 abatement target, as well as three regional policies with the same policy stringency as the national one but only applied to eastern China (EP), the Jiangsu-Shanghai-Zhejiang area (JSZP), and the Beijing-Tianjin-Hebei area (BTHP), respectively. We find that regional policies (EP, JSZP, and BTHP) are as effective in reducing CO2 emissions in their targeted regions as the national policy. However, they lead to an increase in CO2 emissions in untargeted regions (the so-called “emissions leakage”). The CO2 leakage rates, which depend on the policy spatial coverage, are 4%, 13%, and 65% for EP, JSZP, and BTHP, respectively, in 2050. Compared with CO2, changes in air pollutant emissions, including sulfur dioxide (SO2), nitrogen oxides (NOx), volatile organic compounds (VOC), and fine particulate matter (PM2.5), exhibit a similar pattern under all policy scenarios, but the magnitude of change is significantly smaller. Transportation, thermal power, and some energy intensive industries are the three largest contributors to CO2 and air pollutant emission reductions. Our results suggest that regional carbon policy is effective in reducing CO2 and air pollutant emissions in the targeted regions, and extending the spatial coverage or increasing policy stringency can largely inhibit emissions leakage.

Zhang, Wenwen↗

Synchronization of Alternative Models in a Supermodel and the Learning of Critical Behavior

“Supermodeling” climate by allowing different models to assimilate data from one another in run time has been shown to give results superior to those of any one model and superior to any weighted average of model outputs. The only free parameters, connection strengths between corresponding variables in each pair of models, are determined using some form of machine learning. It is demonstrated that supermodeling succeeds because near critical states, interscale interactions are important but unresolved processes cannot be effectively represented diagnostically in any single parameterization scheme. In two examples, a pair of toy quasigeostrophic (QG) channel models of the midlatitudes and a pair of ECHAM5 models of the tropical Pacific atmosphere with a common ocean, supermodels dynamically combine parameterization schemes so as to capture criticality, associated critical structures, and the supporting scale interactions. The QG supermodeling scheme extends a previous configuration in which two such models synchronize with intermodel connections only between medium-scale components of the flow; here the connections are trained against a third “real” model. Intermittent blocking patterns characterize the critical behavior thus obtained, even where such patterns are missing in the constituent models. In the ECHAM-based climate supermodel, the corresponding critical structure is the single ITCZ pattern, a pattern that occurs in neither of the constituent models. In conclusion, for supermodels of both types, power spectra indicate enhanced interscale interactions in frequency or energy ranges of physical interest, in agreement with observed data, and supporting a generalized form of the self-organized criticality hypothesis.

54 ENVIRONMENTAL SCIENCES↗