Search NASA⌕ Search

SEARCH · Search NASA

Results for “digital elevation model”

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 109 records · Page 6

Processed Surface Imagery & Raster Imagery taken onboard TigerShark (U3) UAS

Land-surface orthomosaic imagery derived from the images collected using the Micasense Altum Imager on board the UAS TigerShark (U3) on March 11th, 2021 at an altitude of 2000 ft. The Altum takes images of the land surface across five visible bands and one long wave infrared thermal band. The photogrametry software Agisoft PhotoScan v 1.4 is used to align and stitch the images into a larger composet image using the technique of structure from motion image capture to construct a dense cloud and 3-D model of the surface, which is used to produce a digital elevation model of the terraine surveyed and orthomosaic imagery. This dataset contains the raster calculated orthomosaics of various vegetative indicies, which are used to indicate plant health and land surface characteristics, including Normalized Difference Vegetation Index (NDVI), Green Normalized Difference Vegetation Index (GNDVI), Normalized Difference Water Index (NDWI), surface temperature, and Enhanced_vegetation_index (EVI). These orthomosaics were collected west of the Starkville, MS, airport, and include natural, agricultural, and industrial land coverage. The images have a resolution of 0.2 m/pix.

54 ENVIRONMENTAL SCIENCES↗

New Geologic Map of the Scandia Region of Mars

We have begun work on a sophisti-cated digital geologic map of the Scandia region (Fig. 1) at 1:3,000,000 scale based on post-Viking image and to-pographic datasets. Through application of GIS tools, we will produce a map product that will consist of (1) a printed photogeologic map displaying geologic units and relevant modificational landforms produced by tectonism, erosion, and collapse/mass wasting; (2) a landform geoda-tabase including sublayers of key landform types, attributed with direct measurements of their planform and to-pography using Mars Orbiter Laser Altimeter (MOLA) altimetry data and High-Resolution Stereo Camera (HRSC) digital elevation models (DEMs) and various image datasets; and (3) a series of digital, reconstructed paleostratigraphic and paleotopographic maps showing the inferred distribution and topographic form of materi-als and features during past ages

Tanaka, K. L.↗

Bulk Insolation Models as Predictors for Locations for High Lunar Hydrogen Concentrations

In this study we consider the bulk effects of surface illumination on topography (insolation) and the possible thermodynamic effects on the Moon's hydrogen budget. Insolation is important as one of the dominant loss processes governing distributions of hydrogen volatiles on the Earth, Mars and most recently Mercury. We evaluated three types of high latitude > 65 deg., illumination models that were derived from the Lunar Observing Laser Altimetry (LOLA) digital elevation models (DEM)'s. These models reflect varying accounts of solar flux interactions with the Moon's near-surface. We correlate these models with orbital collimated epithermal neutron measurements made by the Lunar Exploration Neutron Detector (LEND). LEND's measurements derive the Moon's spatial distributions of hydrogen concentration. To perform this analysis we transformed the topographic model into an insolation model described by two variables as each pixels 1) slope and 2) slope angular orientation with respect to the pole. We then decomposed the illumination models and epithermal maps as a function of the insolation model and correlate the datasets.

Hydrogen Concentrations↗

DEEPEN Leapfrog Geodata Model Cleaned and Reformatted Exploration Datasets from Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. As part of the DEEPEN 3D play fairway analysis (PFA) conducted at Newberry Volcano for multiple play types (conventional hydrothermal, superhot EGS, and supercritical), existing geoscientific exploration datasets needed to be acquired, cleaned, reformatted, and assembled in Leapfrog Geothermal. This GDR submission includes all of the cleaned and reformatted (X (m), Y (m), elevation (m), processed data values) datasets used to build the Leapfrog Geodata model. Existing datasets were acquired from the GDR, from AltaRock, and from other sources. This yielded the following datasets: - Digital elevation model produced from LiDAR data by Ramsey and Bard, 2016 - MT surveys from 2006, 2011, 2014, and 2017 (including single inversions) - Gravity surveys from 2006, 2007, and 2011 (including single) - Earthquake catalogs from PNSN, LLNL, and the Newberry EGS Demonstration project - Seismic velocity model from Templeton et al., 2014 - The Frone, 2015 temperature model and a new one produced through extrapolating downhole temperature measurements and the SMU temperature at depth maps. Two versions of the new model are provided: 250 m spacing and 500 m spacing - EarthVision geologic model with alteration from Moser et al., 2016 - Well data from EGS well 55-29, deep geothermal wells, coreholes (GEO N-2 through 5) and several thermal gradient holes - "Newberry Well Data:" Location, simple lithology, directional survey data, and temperature data for the 34 wells and coreholes used in the Newberry PFA Although there are additional 2D datasets available in the area, such as aeromagnetic surveys, these were not included in the analysis. While it may be possible to project these datasets into three dimensions by assuming the surface measurements do not vary with depth, this method is associated with high uncertainty. Preexisting inversions of these data were unavailable, and inverting additional geophysical datasets is outside the scope of this project.

15 GEOTHERMAL ENERGY↗

Mars rover concept development

A structured study effort to develop an extensive, innovative set of mobility and navigation concepts for a planetary exploration vehicle along with the concomitant value system and evaluation tools is presented. A further objective is to submit these concepts to a rigorous, structured evaluation process to derive the most promising candidate systems. To support the evaluation process, a three-layer computer model of the Martian surface was developed, based on the 1/64 deg Digital Elevation Model (DEM) of Mars. Local surface roughness based on measured Martian slope distribution and power spectral density was superimposed on the DEM, and rocks based on Moore's distribution model were added. To assess performance, selected concepts were modeled using DADS, and simulations were run with the vehicle traversing the Martian surface model, including one-meter-high vertical steps and one-meter-wide crevasses. The design details of three promising candidate systems are presented, along with the discussion of their evolution with some recommendations.

Mctamaney, Louis S.↗

Topological Relationship-Based Flow Direction Modeling: Stream Burning and Depression Filling

Flow direction modeling consists of (a) an accurate representation of the river network and (b) digital elevation model (DEM) processing to preserve characteristics with hydrological significance. In part 1 of our study, we presented a mesh-independent approach to representing river networks on different types of meshes. This follow-up part 2 study presents a novel DEM processing approach for flow direction modeling. This approach consists of (a) a topological relationship-based hybrid breaching-filling method to conduct stream burning for the river network and (b) a modified depression removal method for rivers and hillslopes. Our methods reduce modifications to surface elevations and provide a robust two-step procedure to remove local depressions in DEM. They are mesh-independent and can be applied to both structured and unstructured meshes. We applied our new methods with different model configurations to the Susquehanna River Basin. The results show that topological relationship-based stream burning, and depression-filling methods can reproduce the correct river networks, providing high-quality flow direction and other characteristics for hydrologic and Earth system models.

54 ENVIRONMENTAL SCIENCES↗

Image Segmentation, Registration, Compression, and Matching

A novel computational framework was developed of a 2D affine invariant matching exploiting a parameter space. Named as affine invariant parameter space (AIPS), the technique can be applied to many image-processing and computer-vision problems, including image registration, template matching, and object tracking from image sequence. The AIPS is formed by the parameters in an affine combination of a set of feature points in the image plane. In cases where the entire image can be assumed to have undergone a single affine transformation, the new AIPS match metric and matching framework becomes very effective (compared with the state-of-the-art methods at the time of this reporting). No knowledge about scaling or any other transformation parameters need to be known a priori to apply the AIPS framework. An automated suite of software tools has been created to provide accurate image segmentation (for data cleaning) and high-quality 2D image and 3D surface registration (for fusing multi-resolution terrain, image, and map data). These tools are capable of supporting existing GIS toolkits already in the marketplace, and will also be usable in a stand-alone fashion. The toolkit applies novel algorithmic approaches for image segmentation, feature extraction, and registration of 2D imagery and 3D surface data, which supports first-pass, batched, fully automatic feature extraction (for segmentation), and registration. A hierarchical and adaptive approach is taken for achieving automatic feature extraction, segmentation, and registration. Surface registration is the process of aligning two (or more) data sets to a common coordinate system, during which the transformation between their different coordinate systems is determined. Also developed here are a novel, volumetric surface modeling and compression technique that provide both quality-guaranteed mesh surface approximations and compaction of the model sizes by efficiently coding the geometry and connectivity/topology components of the generated models. The highly efficient triangular mesh compression compacts the connectivity information at the rate of 1.5-4 bits per vertex (on average for triangle meshes), while reducing the 3D geometry by 40-50 percent. Finally, taking into consideration the characteristics of 3D terrain data, and using the innovative, regularized binary decomposition mesh modeling, a multistage, pattern-drive modeling, and compression technique has been developed to provide an effective framework for compressing digital elevation model (DEM) surfaces, high-resolution aerial imagery, and other types of NASA data.

Yadegar, Jacob↗

Utilizing Satellite Based Observations and Physical Hydrological Modeling for Freshwater Ecosystem Health in the Lower Mekong River Basin

Freshwater availability is necessary to promote economic growth through agriculture, fisheries, transport, environmental health, and social equity.The National Aeronautics and Space Administration (NASA) and the Conservation International (CI) are partnering to use remote sensing Earth observations to improve regional efforts that assess natural resources for conservation and sustainable management. (Vollmer et al.,2018) have presented the social-ecological framework named the Freshwater Health Index (FHI), which takes account of the interplay between governance, stakeholders, freshwater ecosystems and the ecosystem services they provide.In this work, we develop decision support and making tools for natural resources conservation in the Lower Mekong by leveraging the FHI framework, multiple data products, and hydrological modeling capabilities (Mohammed et al., 2018). Modeling capabilities enable the integration of satellite-based daily gridded precipitation, air temperature, digital elevation model, soil characteristics, and land cover and land use information to simulate water flux framework.

Mohammed, Ibrahim N.↗

Orthographic Stereo Correlator on the Terrain Model for Apollo Metric Images

A stereo correlation method on the object domain is proposed to generate the accurate and dense Digital Elevation Models (DEMs) from lunar orbital imagery. The NASA Ames Intelligent Robotics Group (IRG) aims to produce high-quality terrain reconstructions of the Moon from Apollo Metric Camera (AMC) data. In particular, IRG makes use of a stereo vision process, the Ames Stereo Pipeline (ASP), to automatically generate DEMs from consecutive AMC image pairs. Given camera parameters of an image pair from bundle adjustment in ASP, a correlation window is defined on the terrain with the predefined surface normal of a post rather than image domain. The squared error of back-projected images on the local terrain is minimized with respect to the post elevation. This single dimensional optimization is solved efficiently and improves the accuracy of the elevation estimate.

Terrain Model↗

Hierarchical approach to reconstruct surfaces by using iteratively rectified imagery

A new approach to reconstruct the 3D surface of the object space from digital images is described. All the object points obtained by an automatic orientation procedure lead to a first approximation of the surface. Edges are computed for one image and matched to the other image by grey level correlation or least-squares matching through the scale space. To every discrete step in the scale space there exists the digital stereopair (image pyramid), the corresponding surface digital elevation model (DEM), and the warped images. The warped images in this discrete scale space representation correspond to digital orthophotos obtained from the DEMs that result from matching the image pyramid. The warped images will be used on every successive level in the image pyramid in order to reduce the foreshortening problems associated with any area-based matching method.

Schenk, Toni↗

Assessing Heterogeneity of Surface Water Temperature Following Stream Restoration and a High-Intensity Fire from Thermal Imagery

Thermal heterogeneity of rivers is essential to support freshwater biodiversity. Salmon behaviorally thermoregulate by moving from patches of warm water to cold water. When implementing river restoration projects, it is essential to monitor changes in temperature and thermal heterogeneity through time to assess the impacts to a river’s thermal regime. Lightweight sensors that record both thermal infrared (TIR) and multispectral data carried via unoccupied aircraft systems (UASs) present an opportunity to monitor temperature variations at high spatial (<0.5 m) and temporal resolution, facilitating the detection of the small patches of varying temperatures salmon require. Here, we present methods to classify and filter visible wetted area, including a novel procedure to measure canopy cover, and extract and correct radiant surface water temperature to evaluate changes in the variability of stream temperature pre- and post-restoration followed by a high-intensity fire in a section of the river corridor of the South Fork McKenzie River, Oregon. We used a simple linear model to correct the TIR data by imaging a water bath where the temperature increased from 9.5 to 33.4 °C. The resulting model reduced the mean absolute error from 1.62 to 0.35 °C. We applied this correction to TIR-measured temperatures of wetted cells classified using NDWI imagery acquired in the field. We found warmer conditions (+2.6 °C) after restoration (p < 0.001) and median absolute deviation for pre-restoration (0.30) to be less than both that of post-restoration (0.85) and post-fire (0.79) orthomosaics. In addition, there was statistically significant evidence to support the hypothesis of shifts in temperature distributions pre- and post-restoration (KS test 2009 vs. 2019, p < 0.001, D = 0.99; KS test 2019 vs. 2021, p < 0.001, D = 0.10). Moreover, we used a Generalized Additive Model (GAM) that included spatial and environmental predictors (i.e., canopy cover calculated from multispectral NDVI and photogrammetrically derived digital elevation model) to model TIR temperature from a transect along the main river channel. This model explained 89% of the deviance, and the predictor variables showed statistical significance. Collectively, our study underscored the potential of a multispectral/TIR sensor to assess thermal heterogeneity in large and complex river systems.

Barker, Matthew I. (ORCID:0000000252864930)↗

Peru Water Resources: Integrating NASA Earth Observations into Water Resource Planning and Management in Perus La Libertad Region

Developing countries often struggle with providing water security and sanitation services to their populations. An important aspect of improving security and sanitation is developing a comprehensive understanding of the country's water budget. Water For People, a non-profit organization dedicated to providing clean drinking water, is working with the Peruvian government to develop a water budget for the La Libertad region of Peru which includes the creation of an extensive watershed management plan. Currently, the data archive of the necessary variables to create the water management plan is extremely limited. Implementing NASA Earth observations has bolstered the dataset being used by Water For People, and the METRIC (Mapping EvapoTranspiration at High Resolution and Internalized Calibration) model has allowed for the estimation of the evapotranspiration values for the region. Landsat 8 imagery and the DEM (Digital Elevation Model) from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) sensor onboard Terra were used to derive the land cover information, and were used in conjunction with local weather data of Cascas from Peru's National Meteorological and Hydrological Service (SENAMHI). Python was used to combine input variables and METRIC model calculations to approximate the evapotranspiration values for the Ochape sub-basin of the Chicama River watershed. Once calculated, the evapotranspiration values and methodology were shared Water For People to help supplement their decision support tools in the La Libertad region of Peru and potentially apply the methodology in other areas of need.

Padgett-Vasquez, Steve↗

Snow distribution patterns revisited: A physics-based and machine learning hybrid approach to snow distribution mapping in the sub-Arctic Supporting Data

Snow in the Arctic and sub-Arctic is highly variable at fine scales, with deep drifts and shallow scoured areas creating a complex pattern of snow on the landscape. This fine-scale variation in snow is driven primarily by landscape and vegetation properties. Some landscape features, such as river beds, will rapidly fill in with snow during the wintertime due to high winds, while shrubs will trap blowing snow, resulting in drifts. Meanwhile, snow will blow off of exposed areas, resulting in abnormally shallow snow. These complex interactions between wind, vegetation, and terrain are difficult to represent well with physically-based models, but machine learning techniques have shown promise in the past. Here, we propose a hybrid modeling approach, where we use machine learning derived snow pattern maps to inform SnowModel, a physically-based snow process model. We develop and test this technique at the Teller 27 Seward Peninsula NGEE-Arctic study site. This dataset includes 5 *.nc files of model inputs and outputs plus one user guide (*.pdf). We present the data we use to drive SnowModel (vegetation type, digital elevation model), the snow pattern maps used to inform SnowModel (Standardized Depth Values maps, Machine Learning Snow Distribution Pattern), and our machine learning, SnowModel, and Hybrid snow depth and snow water equivalent results. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Improved LOLA Elevation Maps for South Pole Landing Sites: Error Estimates and Their Impact on Illumination Conditions

We present new high-resolution topographic models of 4 high-priority lunar south pole landing sites based exclusively on the laser altimetry data acquired by the Lunar Orbiter Laser Altimeter (LOLA) onboard the Lunar Reconnaissance Orbiter. By iteratively adjusting the LOLA tracks to the LOLA-based digital elevation model (LDEM) in a self-consistent fashion, we reduce the orbital geolocation errors by over a factor of 10 such that the new ground track geolocation uncertainty is ~10–20 ​cm horizontally and ~2–4 ​cm vertically over each 16 ​× ​16 km region. These new and improved 5 ​m/pix LDEMs will be useful to constrain higher-resolution topographic models derived from imagery, which are not as well controlled geodetically and which can be hindered by shadows. We developed a method to estimate surface height uncertainty in the new LDEMs, which accounts for the reduced orbital errors and interpolation errors by assuming a fractal behavior for the short-scale topography. The LDEM surface height and slope uncertainties have typical RMS values of ~0.30–0.50 ​m and ~1.5–2.5°, respectively. Finally, we examine how height uncertainties propagate to variations in horizon elevation and thus the predicted illumination conditions at these polar latitudes, and we show how this error characterization can inform landing site studies.

Michael K Barker↗

Flight Test Results of a Synthetic Vision Elevation Database Integrity Monitor

This paper discusses the flight test results of a real-time Digital Elevation Model (DEM) integrity monitor for Civil Aviation applications. Providing pilots with Synthetic Vision (SV) displays containing terrain information has the potential to improve flight safety by improving situational awareness and thereby reducing the likelihood of Controlled Flight Into Terrain (CFIT). Utilization of DEMs, such as the digital terrain elevation data (DTED), requires a DEM integrity check and timely integrity alerts to the pilots when used for flight-critical terrain-displays, otherwise the DEM may provide hazardous misleading terrain information. The discussed integrity monitor checks the consistency between a terrain elevation profile synthesized from sensor information, and the profile given in the DEM. The synthesized profile is derived from DGPS and radar altimeter measurements. DEMs of various spatial resolutions are used to illustrate the dependency of the integrity monitor s performance on the DEMs spatial resolution. The paper will give a description of proposed integrity algorithms, the flight test setup, and the results of a flight test performed at the Ohio University airport and in the vicinity of Asheville, NC.

deHaag, Maarten Uijt↗

Application of digital terrain data to quantify and reduce the topographic effect on LANDSAT data

Integration of LANDSAT multispectral scanner (MSS) data with 30 m U.S. Geological Survey (USGS) digital terrain data was undertaken to quantify and reduce the topographic effect on imagery of a forested mountain ridge test site in central Pennsylvania. High Sun angle imagery revealed variation of as much as 21 pixel values in data for slopes of different angles and aspects with uniform surface cover. Large topographic effects were apparent in MSS 4 and 5 was due to a combination of high absorption by the forest cover and the MSS quantization. Four methods for reducing the topographic effect were compared. Band ratioing of MSS 6/5 and MSS 7/5 did not eliminate the topographic effect because of the lack of variation in MSS 4 and 5 radiances. The three radiance models examined to reduce the topographic effect required integration of the digital terrain data. Two Lambertian models increased the variation in the LANDSAT radiances. The nonLambertian model considerably reduced (86 per cent) the topographic effect in the LANDSAT data. The study demonstrates that high quality digital terrain data, as provided by the USGS digital elevation model data, can be used to enhance the utility of multispectral satellite data.

Justice, C. O.↗

Applying Transfer Learning for Street-Scale Nuisance Flood Forecasting in Coastal-Urban Cities

An important challenge with Machine Learning (ML) is its transferability; that is, whether a ML model trained on one set of data can be applied to a second set of data without requiring a full re-training of the model. Transfer Learning (TL) addresses this challenge by transferring knowledge learned in the source domain (the data it was trained on) to the target domain (a second set of data that is statistically different but related, which the model was not trained on). This study investigates the use of TL for street-scale nuisance flood forecasting by exploring whether a ML model trained on data collected for one set of streets can effectively forecast flooding for another set of streets in the same city using TL. The envisioned use case is a city deploying a new flood depth monitoring sensor on a street and using TL to apply a ML model, trained on sensor data from an existing flood depth sensor network, to this new street. Eventually, the new flood depth sensor will have a sufficient dataset for training its own ML model, but TL can be used to fill the gap in time while this new dataset is being generated. This method is explored using a Long Short-Term Memory (LSTM) model trained on data for the flood-prone streets of Norfolk City, Virginia. The data used for training includes environmental time series (rainfall, tide), topographic features (Digital Elevation Model (DEM), Topographic Wetness Index (TWI), Depth To Water (DTW)), and street-scale flood depth time series obtained from a high-fidelity physics-based model, acting as a synthetic street-scale stream depth sensor dataset since actual stream depth sensor data is generally unavailable for most cities. A set of 180 flood-prone streets was used to train a base model, while another set of 180 flood-prone streets was used to re-train that model using different TL strategies. The results show that full-weight re-training proved most effective and minimal re-training of only the output layer was insufficient. The advantage of TL was most pronounced when target data was limited, meaning data collected at the new water depth sensor location included generally less than 18 flood events. As target data increased beyond 18 flood events, the benefit of TL diminished relative to training a ML model directly on the local flood events. These findings can assist cities as they implement street-scale flood sensing systems to create accurate forecasts for new sensing locations that do not yet have sufficient data records to train a local ML model.

Roy, Binata [Univ. of Virginia, Charlottesville, V↗

2D Dam-Break Analysis of L Lake and PAR Pond Dams Using HEC-RAS

In 1991 a dam-break study was conducted for the high hazard dams located at L Lake and PAR Pond on the Savannah River Site. Two scenarios were considered, over topping from a Probable Maximum Flood (PMF), and a fair weather dam-break for either or both dams. Unfortunately, no inundation map was developed from the study. The purpose of this project was to redo the original dam-break study with improved data and methodology to generate Inundation maps to assist with emergency response and evacuation plans. The Hydrologic Engineering Center's River Analysis System (HEC-RAS) is a free to download river analysis modeling program developed by the US Army Corps of Engineers capable of 1D and 2D hydraulic calculations. Version 5.0.7 (released March 2019) was used for this project. Digital elevation models for the area were retrieved from the US Geological Survey database and converted to a .hdf file within the program. Then the 2D flood area was identified from the contours. Both L Lake and PAR Pond were inputted as 1D storage areas because DEM data does not contain elevation values under water bodies. An elevation vs volume curve was available for both storage areas. Initial elevations were set for both scenarios. Both dams are earthen dams. The Steel Creek dam at L Lake has 6 ft diameter conduit with an upper and lower sluice gate. PAR pond dam consist of a weir and sluice gate connected to an 8 x 8 ft channel. Both outlets were modeled with a pool elevation vs discharge curve. The steel creek dam sluice gates were assumed to be fully open in all cases. As in the previous study, the dams were set to breach when they were overtopped by 1.5 ft during PMF conditions (Figure 1). A fair weather breach was set to be due to a piping failure (Figure 1). In the dual dam break during fair weather conditions the PAR pond dam fails 3 hours after Steel Creek to achieve maximum flooding in the down stream reaches. Simulation was run 3 days for each case and with a 1 minute computational interval. Maximum flooding occurs under PMF conditions with the failure of both dams. PAR pond dam fails first 16 hours and 32 minutes after the start of the simulation with the Steel creek dam failing 6 minutes later. In all cases, the bridges and roads spanning Steel Creek and Lower Three Runs will be inundated and potentially washed away. The Burtons Ferry Highway south of the storage areas will be partially flooded during PMF failure, dual fair weather failure, and PAR pond failure under fair weather conditions.

54 ENVIRONMENTAL SCIENCES↗