SEARCH · Search NASA

Results for “process simulation 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 595 records · Page 33

All sources NASA NTRS DOE OSTI

DOE OSTI2019

Negative extreme events in gross primary productivity and their drivers in China during the past three decades

Climate extremes have remarkable impacts on ecosystems and are expected to increase with future global warming. However, only few studies have focused on the ecological extreme events and their drivers in China. In this study, we carried out an analysis of negative extreme events in gross primary productivity (GPP) in China and the sub-regions during 1982–2015, using monthly GPP simulated by 12 process-based models (TRENDYv6) and an observation-based model (Yao-GPP). Extremes were defined as the negative 5th percentile of GPP anomalies, which were further merged into individual extreme events using a three-dimensional contiguous algorithm. Spatio-temporal patterns of negative GPP anomalies were analyzed by taking the 1000 largest extreme events into consideration. Results showed that the effects of extreme events decreased annual GPP by 2.8% (i.e. 208 TgC year –1 ) in TRENDY models and 2.3% (i.e. 151 TgC year –1 ) in Yao-GPP. Hotspots of extreme GPP deficits were mainly observed in North China (–53 gC m –2 year –1 ) in TRENDY models and Northeast China (–42 gC m –2 year –1 ) in Yao-GPP. For China as a whole, attribution analyses suggested that extreme low precipitation was associated with 40%–50% of extreme negative GPP events. Most events in northern and western China could be explained by meteorological droughts (i.e. low precipitation) while GPP extreme events in southern China were more associated with temperature extremes, in particular with cold spells. GPP was revealed to be much more sensitive to heat/drought than to cold/wet extreme events. Here, combined with projected changes in climate extremes in China, GPP negative anomalies caused by drought events in northern China and by temperature extremes in southern China might be more prominent in the future.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

A Lake Biogeochemistry Model for Global Methane Emissions: Model Development, Site‐Level Validation, and Global Applicability

Abstract Lakes are important sentinels of climate change and may contribute over 30% of natural methane (CH 4 ) emissions; however, no earth system model (ESM) has represented lake CH 4 dynamics. To fill this gap, we refined a process‐based lake biogeochemical model to simulate global lake CH 4 emissions, including representation of lake bathymetry, oxic methane production (OMP), the effect of water level on ebullition, new non‐linear CH 4 oxidation kinetics, and the coupling of sediment carbon pools with in‐lake primary production and terrigenous carbon loadings. We compiled a lake CH 4 data set for model validation. The model shows promising performance in capturing the seasonal and inter‐annual variabilities of CH 4 emissions at 10 representative lakes for different lake types and the variations in mean annual CH 4 emissions among 106 lakes across the globe. The model reproduces the variations of the observed surface CH 4 diffusion and ebullition along the gradients of lake latitude, depth, and surface area. The results suggest that OMP could play an important role in surface CH 4 diffusion, and its relative importance is higher in less productive and/or deeper lakes. The model performance is improved for capturing CH 4 outgassing events in non‐floodplain lakes and the seasonal variability of CH 4 ebullition in floodplain lakes by representing the effect of water level on ebullition. The model can be integrated into ESMs to constrain global lake CH 4 emissions and climate‐CH 4 feedback.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2015

Ultrafast defect dynamics: A new approach to all optical broadband switching employing amorphous selenium thin films

Optical switches offer higher switching speeds than electronics, however, in most cases utilizing the interband transitions of the active medium for switching. As a result, the signal suffers heavy losses. In this article, we demonstrate a simple and yet efficient ultrafast broadband all-optical switching on ps timescale in the sub-bandgap region of the a-Se thin film, where the intrinsic absorption is very weak. The optical switching is attributed to short-lived transient defects that form localized states in the bandgap and possess a large electron-phonon coupling. We model these processes through first principles simulation that are in agreement with the experiments.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Linking Transient Voltage to Spatially-Resolved Luminescence Imaging to Understand Reliability of Perovskite Photovoltaics

In this work, we present a methodology to separate effects of perovskite device metastability from irreversible degradation, using stress/rest cycling under constant current bias while collecting a series of electroluminescence images and continuously monitoring voltage. We develop a simulation model and procedures for image processing to better understand the effects of ion parameters on the transient nature of voltage and evolving electroluminescence images.

electroluminescence↗

DOE OSTI2022

Development of a 250 Degree C 15kV Supercascode Switch Using SiC JFET Technology

Tetra Corporation is developing a pulsed power drilling system for the geothermal industry. The system requires a high temperature switching technology to replace the existing technology based on Si thyristors with a temperature limit of 150 degrees C. This work focuses on the development of such switch based on a Supercascode configuration using SiC JFET technology. The development presented contains thermal models, packaging design and processes, circuit simulations and prototypes of both the Supercascode switch and the high temperature packaged SiC JFETs.

cascaded switch↗

DOE OSTI2025

Model Scripts for "Old-Aged Groundwater Contributes to Mountain Hillslope Hydrologic Dynamics"

The partitioning of water inputs between deep and shallow groundwater flow paths is a fundamental processes, yet is challenging to observe. Numerical models provide a valuable tool to further develop insights on these groundwater mixing processes. This repository contains scripts to run the ParFlow-CLM and EcoSLIM integrated hydrologic models along the Pumphouse Hillslope in the East River Watershed and associated python scripts to process model outputs. Files includes input decks for the associated models, python scripts, and text and css files to support the model runs and processing. The 2-D hillslope model simulates from 2000 to 2021 using transient forcing conditions from the publicly available NLDAS-2 dataset. We develop an ensemble of models with variable hydrogeologic and soil parameters. Model outputs from the ensemble of runs are compared to to water level observations at the PLM1 and PLM6 wells published at (https://data.ess-dive.lbl.gov/view/doi:10.15485/1866836) and groundwater mean ages published at (https://data.ess-dive.lbl.gov/datasets/doi:10.15485/1960042. The model results are used to evaluate hillslope water mass-balance transience and mixing dynamics between groundwater flow paths with young (<10 year) and old (>10 year) ages.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Diagnosis of convective organization and cold pools using ARM datasets and evaluation of a unified convection parameterization (UNICON)

Tropical thunderstorms often cluster together. Studies have suggested that the degree to which the tropical thunderstorms are clustered impacts Earth's energy balance and water cycle, as well as extreme precipitation events. However, the processes controlling the spatial distribution of the thunderstorms are poorly understood and are not properly represented in most computer models for weather and climate prediction. Under the goals of better understanding how convection organizes at the mesoscale and advancing the representation of mesoscale convective organization in global climate models, we i) objectively quantified the degrees of convective organization and diagnosed cold pool processes using ARM field campaign observations, ii) examined the organization processes in storm-resolving model simulations, and iii) evaluated the impacts of mesoscale convective organization in global model simulations. The project yielded a firm reference against which the global model representation of mesoscale convective organization and cold pools can be evaluated against and shed new light into the role of parameterized convective organization in the global model simulation of the basic state and variability. Our results revealed two distinct phases of convective clustering during the two-day rain episodes (Cheng et al. 2018) and a new mechanism through which vertical wind shear in the low-troposphere can aid convective organization over tropical oceans (Cheng et al. 2020). It was demonstrated that the interactive representation of cold pools and mesoscale convective organization is key for global models to successfully simulate both the mean state and intraseasonal variability in the tropics (Ahn et al. 2019; 2020).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Oxy-Combustion System Process Optimization

The overall objective of this work is to develop a new chemical absorbent-based, high pressure, CO 2 purification system to remove the residual oxygen that currently contaminates the recovered CO 2 , and to optimize the Pressurized Oxy-Combustion (POxC) process to minimize the Cost of Electricity (COE) generated in this advanced combustion process. TDA developed and validated the performance of the oxygen removal system for CO 2 purification. In collaboration with the Advanced Power & Energy Program (APEP) of University of California, Irvine (UCI), we optimized the POxC process, including thermal management, heat integration, and power cycle optimization using process design and modeling supported with Aspen Plus® process simulations. The techno-economic analysis results indicate that the pressurized oxycombustion coal power plant with Ion-transport membrane (ITM) air separation unit (ASU) (Case 2 – 30.55%) does not show an advantage over a cryogenic ASU (Case 1 – 31.24%) while TDA’s sorption-based ASU (Case 4 - 32.61%) shows a significant advantage over the cryogenic ASU (Case 1 – 331.24%). The specific plant costs show a wide range with a low of $2544/kW for Case 11C (TDA ASU, co-sequestering the SO x , and ultra-supercritical steam cycle) to a high of $2975/kW for Case 12A (cryogenic ASU and sCO 2 cycle). In general, the ITM ASU based cases have lower specific plant costs than the corresponding cryogenic ASU based cases while the TDA ASU based cases show the lowest specific plant costs. The main reason for lowering these costs is the higher overall plant thermal efficiency which decreases the plant cost on a $ per kW basis. Next comparing the cases with different power cycle working fluid conditions in terms of temperature and pressure while all utilizing steam, similar trends as the plant costs may be observed. However, with the supercritical CO 2 (sCO 2 ) cycle, the increase in thermal efficiency of the sCO 2 cycle was not able to offset its increase in plant cost making the plant costs higher than those of the corresponding steam cycle cases. The Cost of Electricity (COE) again shows similar trends as the specific plant costs. The COE for Case 11C at $110.1/MWh is also the lowest, but among all cases that do not co-sequester the SO x , Case 12C (TDA ASU and sCO 2 cycle) has the lowest COE at $\$$117.5/MWh while the highest is for Case 8A (cryogenic ASU and supercritical steam cycle) at $130.4/MWh.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2022

LDRD-DR: Accelerating Nuclear Fuel Qualification through Integrated Multiscale and Multiphysics Models [Slides]

Nuclear fuel qualification implies ensuring that a new fuel can safely operate in a reactor under both normal and accident conditions. Historically, extensive irradiation campaigns probe fuel performance relationships and interconnections in an integral sense, which is very time-consuming and costly. The presentation focus is on introducing advanced modeling and simulation into the qualification process in order to reduce the qualification timeline by targeting fewer integral tests

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Advanced Structured Adsorbent Architectures for Transformative Carbon Dioxide Capture Performance (Final Report)

Svante is a world leader at using solid sorbents for low-cost Carbon Dioxide (CO 2 ) capture, a technology which is recognized as critical in meeting the dual mandates of energy security/reliability and the mitigation of man-made CO 2 emissions. Svante has been developing proprietary adsorbent material compositions, forming them into structured laminates, developing and optimizing process cycles, and system design for efficient capture of CO 2 from post-combustion flue gases of thermal power plants and industrial facilities. The deployment of first-generation CO 2 capture technology has been significantly hampered by high costs and energy penalties, among other barriers. Second generation CO 2 capture technologies (including the Mark I variant of Svante’s Veloxotherm™ adsorption-based technology), utilizing single adsorbent architecture, show promise for reducing the barriers to deploying CO 2 capture plants in commercially meaningful numbers. The objective of this project was to evaluate the Recipient’s transformational (Mark-II) VeloxoTherm™ Technology via the development and bench-scale testing of an advanced structured adsorbent, including novel Bi-layer, laminated adsorbent structures and segmented beds. Svante selected, synthesized, and characterized tailored solid adsorbents for computational modeling, advanced structured adsorbent development, process simulations, and dynamic bench scale (~1-10 kg/day CO 2 captured) testing using an existing single-bed VeloxoTherm™ Station (VTS) coupled with a natural gas-fired boiler. Segmented beds used the in-house, multi-bed Process Demonstration Unit (PDU) to demonstrate key performance indicators (KPIs), such as recovery, product purity, regeneration energy, and the integrated system's productivity in lifetime analysis. Segmented beds were used at a 1 tonne per day (TPD) unit at an industrial site to provide bench-scale validation of performance in an industrial setting. Svante was developing and optimizing the post-combustion CO 2 adsorption technology architectures, including the Bi-layer and segmented laminated adsorbent structure design, integrated rapid cycle temperature swing adsorption (RC-TSA) cycle, flow path architecture, and adsorbent bed construction and packaging (including gas porting) to progress towards achievement of DOE’s Transformational CO 2 Capture goals of 95% CO 2 purity and a cost of electricity at least 30% lower than a supercritical Pulverized Coal (PC) power plant with CO 2 capture, or approximately $30 per tonne of CO 2 captured ready for demonstration by 2030. The main requirements to reach the DOE target cost of carbon capture below $30/MT using Rapid-Cycling Temperature Swing Adsorption (RC-TSA) are as follows: (1) Increased capacity at different CO 2 concentrations, (2) Increased sorbents cycle life, (3) Increased O 2 resistance, and (4) Decreased steam requirement to extract 1 kg of CO 2 .

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2024

Dynamic Binary Complexes (DBC) as Super-Adjustable Viscosity Modifiers for Hydraulic Fracturing Fluids

In the preceding project year two, we refined three DBC formulations from a selection of over 50 different chemistries. The optimization study primarily encompassed testing for (i) reversibility extent, (ii) performance in the presence of chemical additives, (iii) adhesion and friction behavior during displacement in wellbores and pipelines, (iv) corrosion protection performance, and (v) injection performance with model fracture systems at the laboratory scale. Highly promising results obtained from all these tests signify the significant potential of DBCs in enhancing hydrocarbon recovery from unconventional reservoirs. The primary activities in the third project year included publishing experimental findings across multiple articles and conducting outreach initiatives. Throughout the year, we undertook tasks such as replicating experimental results, further optimizing various formulations and their associated experimental sets, and conducting additional tests to address missing components based on reviewer feedback and suggestions. We also explored the surfactant and friction-reduction aspects of selected formulations through drag reduction tests. In addition, we constructed an improved fracturing performance setup and performed flow injection tests. The specific DBC formulations focused on during this project period were A8/B1, A12/B5, and A10/B12. We also obtained results for additional DBC formulations and a select few commercial fracturing fluids for the purpose of comparison. Within the project's scope, we aim to enhance the experimental findings with the development of various models. The first two years focused on two key aspects: (i) the creation of a high-fidelity hydraulic fracturing model for non-Newtonian fluids to gain insights into the implementation of DBC fluids in fracking environments, and (ii) the development of a multiphase flow simulator for estimating total production, fluid saturation in the reservoir, and the creation of a fracture propagation model and kinetic Monte Carlo (kMC) models for diverse applications. In the third year, we delved into the fundamental nanostructural properties of DBCs, exploring aspects such as material chemistry, pH tunability, and control of DBC formation and stability. Subsequently, in the extension year, we conducted a systematic investigation of various building blocks containing primary, secondary, and tertiary amine functional groups to understand their impact on rheological and viscoelastic properties. Furthermore, we explored a Dissipative Particle Dynamics (DPD) model to simulate self-assembly processes with precision, creating a high-fidelity representation of relevant nanostructures. The tasks performed this year with the significant results obtained have been discussed in Section 2.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Linking Transient Voltage to Spatially-Resolved Luminescence Imaging to Understand Reliability of Perovskite Photovoltaics: Preprint

41 EE - Solar Energy Technologies Office (EE-4S)↗

DOE OSTI2024

Recovery of valuable metals from electronic waste using a novel ammonia-based hydrometallurgical process

The growing quantity of waste electrical and electronic equipment (WEEE), also known as electronic waste (E-waste), has been an area of growing public concern. The abundance of valuable metals contained in waste printed circuit boards (WPCBs) has made it a promising secondary resource, especially for Cu and Au. Although the recovery of metals from WPCBs via hydrometallurgical routes has been studied extensively over the past 20 years, most of the research has been limited in the laboratory. In current hydrometallurgical processes, strong acids and expensive oxidizers are often used to ensure a high recovery of metals without considering the sustainability aspects of the environment and economics. To improve upon current hydrometallurgical offerings, the current study seeks to develop an energy-saving, environmentally friendly, economic and sustainable process to efficiently recover the valuable metals from real-world end-of-life WPCBs. The new contributions presented in this study are 1) design and evaluation of a comprehensive hydrometallurgical flowsheet; 2) employment of real end-of-life PCBs as feed materials in an investigation on Cu-NH3 leaching kinetics; 3) further application of kinetic model on a counter flow process simulation; and 4) evaluation of the influences by co-existing metals in Au-S2O3 leaching and recommendation for favorable leaching conditions.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2024

Transitioning from Simulation to Reality: Applying Chatter Detection Models to Real-World Machining Data

Chatter, a self-excited vibration phenomenon, is a critical challenge in high-speed machining operations, affecting tool life, product surface quality, and overall process efficiency. While machine learning models trained on simulated data have shown promise in detecting chatter, their real-world applicability remains uncertain due to discrepancies between simulated and actual machining environments. The primary goal of this study is to bridge the gap between simulation-based machine learning models and real-world applications by developing and validating a Random Forest-based chatter detection system. This research focuses on improving manufacturing efficiency through reliable chatter detection by integrating Operational Modal Analysis (OMA), Receptance Coupling Substructure Analysis (RCSA), and Transfer Learning (TL). The study applies a Random Forest classification model trained on over 140,000 simulated machining datasets, incorporating techniques like Operational Modal Analysis (OMA), Receptance Coupling Substructure Analysis (RCSA), and Transfer Learning (TL) to adapt the model for real-world operational data. The model is validated against 1600 real-world machining datasets, achieving an accuracy of 86.1%, with strong precision and recall scores. The results demonstrate the model’s robustness and potential for practical implementation in industrial settings, highlighting challenges such as sensor noise and variability in machining conditions. This work advances the use of predictive analytics in machining processes, offering a data-driven solution to improve manufacturing efficiency through more reliable chatter detection.

42 ENGINEERING↗

DOE OSTI2026

Assessing the Accuracy of Property Model Predictions for Cost Optimization of Desalination Technologies

Accurate modeling of seawater thermophysical and thermodynamic properties is critical for optimizing desalination processes. This study compares three seawater property models, a Reaktoro multicomponent model, the thermophysical seawater properties library from the Massachusetts Institute of Technology, and a simplified sodium chloride model, in the context of levelized cost of water (LCOW) minimization for reverse osmosis (RO) and mechanical vapor compression systems. Process simulations and cost optimizations reveal that although all three models yield comparable LCOW and specific energy consumption (SEC) estimates under baseline conditions, deviations among their predictions increase with salinity. Relative differences in LCOW and SEC reach up to 6% and 8%, respectively. RO results show greater variability due to differences in osmotic pressure predictions, which affect pressure constraints at high recoveries. Computational performance varies substantially; specifically, Reaktoro simulations are up to 28 times slower than empirical models due to their detailed equilibrium calculations. These results suggest that empirical models offer acceptable accuracy for routine desalination process design, while Reaktoro provides advantages in scenarios requiring detailed speciation, such as scaling or pH adjustment studies. These findings underscore the importance of selecting appropriate property models based on the modeling objective of desalination applications and motivate future work integrating thermodynamic rigor with empirical efficiency.

Physical and chemical properties↗

DOE OSTI2017

Metal Big Area Additive Manufacturing: Process Modeling and Validation

Metal Big Area Additive Manufacturing (mBAAM) is a new additive manufacturing (AM) technology for printing large-scale 3D objects. mBAAM is based on the gas metal arc welding process and uses a continuous feed of welding wire to manufacture an object. An electric arc forms between the wire and the substrate, which melts the wire and deposits a bead of molten metal along the predetermined path. In general, the welding process parameters and local conditions determine the shape of the deposited bead. The sequence of the bead deposition and the corresponding thermal history of the manufactured object determine the long range effects, such as thermal-induced distortions and residual stresses. Therefore, the resulting performance or final properties of the manufactured object are dependent on its geometry and the deposition path, in addition to depending on the basic welding process parameters. Physical testing is critical for gaining the necessary knowledge for quality prints, but traversing the process parameter space in order to develop an optimized build strategy for each new design is impractical by pure experimental means. Computational modeling and optimization may accelerate development of a build process strategy and saves time and resources. Because computational modeling provides these opportunities, we have developed a physics-based Finite Element Method (FEM) simulation framework and numerical models to support the mBAAM process s development and design. In this paper, we performed a sequentially coupled heat transfer and stress analysis for predicting the final deformation of a small rectangular structure printed using the mild steel welding wire. Using the new simulation technologies, material was progressively added into the FEM simulation as the arc weld traversed the build path. In the sequentially coupled heat transfer and stress analysis, the heat transfer was performed to calculate the temperature evolution, which was used in a stress analysis to evaluate the residual stresses and distortions. In this formulation, we assume that physics is directionally coupled, i.e. the effect of stress of the component on the temperatures is negligible. The experiment instrumentation (measurement types, sensor types, sensor locations, sensor placements, measurement intervals) and the measurements are presented. The temperatures and distortions from the simulations show good correlation with experimental measurements. Ongoing modeling work is also briefly discussed.

Additive Manufacturing↗

DOE OSTI2015

Evaluation of D and D alternatives for a DOE nuclear material processing facility - 15265

Ten scenarios for SRS Building 235-F were evaluated relative to groundwater protection through use of a 235-F GoldSim{sup TM} fate and transport model. The scenarios included two extreme scenarios, a no action scenario, four D and D grouting scenarios, and three D and D inventory removal scenarios. The fate and transport model simulated contaminant release from four 235-F process areas and the 294-2F Sand Filter and fate and transport through the vadose zone, the Upper Three Runs (UTR) aquifer, and the UTR creek. The model was a stochastic model for which 1,000 realizations were run over 100,000 years for each scenario. In association with the seven D and D alternatives only three groundwater standards out of 25 were exceeded, with the combined radium standard showing the greatest exceedance. The standards were only exceeded at the edge of the building after 10,000 years. It was estimated that there was a 90% likelihood that the combined radium standard would be exceeded for all four D and D grouting scenarios and that a 95% inventory removal from process cells 1-5 would result in an essentially zero percent likelihood that any standards would be exceeded. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2025

Data-model files associated with the manuscript "Modeling the Effects of Wetland Restoration on Coastal Hydrology: A Case Study of Elkhorn Slough Watershed, California"

This package contains the data, simulation setups, notebooks and figures used in “Modeling the Effects of Wetland Restoration on Coastal Hydrology: A Case Study of Elkhorn Slough Watershed, California” (Xu et al., 2025). In this study, we selected Elkhorn Slough, a tidal estuary, in California, to investigate the impact of wetland restoration and sea level rise on coastal hydrology using the process-based coastal hydrologic model, Advanced Terrestrial Simulator (ATS), informed by site-specific data. We designed a novel modeling workflow for incorporating wetland restoration features into land cover and soil properties for the model parameterization. The validation results demonstrate a strong agreement between modeled and observed data. We studied the characteristics of coastal watershed hydrology, then focused on the surface water dynamics at two wetland sites within Elkhorn Slough, a reference site and a restored site. Our simulation results indicate that the restored site successfully maintains surface elevation, resulting in reduced surface inundation. We also examined the impact of wetland restoration under expected sea level rise over the next few decades. The low-lying Yampah Marsh, the reference site, is likely to be inundated due to future sea level rise when highest tides arrive; while a higher percentage of Hester Marsh, the restored site, would retain marsh vegetation in coming decades, regardless of tidal conditions. Our study provides important information for examining the outcome of restoration practices that include surface elevation in tidal wetlands under climate changes.Several files can be found from this data package.1. README.md: This file describes the title, journal, co-authors, abstract, repository structure and model version.2. Simulation_Setups.zip: The file contains the model configuration files (XML format) for ATS. 3. Notebooks.zip: The file contains the Jupyter notebooks for generating the pre- and post-restoration meshes and the meshes of future scenarios. 4. Figures.zip: The file contains the figures used in the manuscript.5. Data.zip: The file contains the data used to drive the model simulations, including watershed and wetlands boundaries, mesh files and references to additional datasets (e.g., meteorological forcing, tidal dataset, DEMs, land cover, soil properties). Also, it contains water level observations at the restored wetland.

54 ENVIRONMENTAL SCIENCES↗