Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spatial statistics”

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,387 records · Page 77

An Initial Assessment of a SMAP Soil Moisture Disaggregation Scheme Using TIR Surface Evaporation Data over the Continental United States

The Soil Moisture Active Passive (SMAP) mission is dedicated toward global soil moisture mapping. Typically, an L-band microwave radiometer has spatial resolution on the order of 36-40 km, which is too coarse for many specific hydro-meteorological and agricultural applications. With the failure of the SMAP active radar within three months of becoming operational, an intermediate (9-km) and finer (3-km) scale soil moisture product solely from the SMAP mission is no longer possible. Therefore, the focus of this study is a disaggregation of the 36-km resolution SMAP passive-only surface soil moisture (SSM) using the Soil Evaporative Efficiency (SEE) approach to spatial scales of 3-km and 9-km. The SEE was computed using thermal-infrared (TIR) estimation of surface evaporation over Continental U.S. (CONUS). The disaggregation results were compared with the 3 months of SMAP-Active (SMAP-A) and Active/Passive (AP) products, while comparisons with SMAP-Enhanced (SMAP-E), SMAP-Passive (SMAP-P), as well as with more than 180 Soil Climate Analysis Network (SCAN) stations across CONUS were performed for a 19 month period. At the 9-km spatial scale, the TIR-Downscaled data correlated strongly with the SMAP-E SSM both spatially (r = 0.90) and temporally (r = 0.87). In comparison with SCAN observations, overall correlations of 0.49 and 0.47; bias of −0.022 and −0.019 and unbiased RMSD of 0.105 and 0.100 were found for SMAP-E and TIR-Downscaled SSM across the Continental U.S., respectively. At 3-km scale, TIR-Downscaled and SMAP-A had a mean temporal correlation of only 0.27. In terms of gain statistics, the highest percentage of SCAN sites with positive gains (>55%) was observed with the TIR-Downscaled SSM at 9-km. Overall, the TIR-based downscaled SSM showed strong correspondence with SMAP-E; compared to SCAN, and overall both SMAP-E and TIR-Downscaled performed similarly, however, gain statistics show that TIR-Downscaled SSM slightly outperformed SMAP-E.

Disaggregation↗

Impact of solar geoengineering on temperatures over the Indonesian Maritime Continent

Climate change has been projected to increase the intensity and magnitude of extreme temperature in Indonesia. Solar radiation management (SRM) has been proposed as a strategy to temporarily combat global warming, buying time for negative emissions. Though the global impacts of SRM have been extensively studied in recent years, regional impacts, especially in the tropics, have received much less attention. This paper investigates the potential stratospheric sulfate aerosol geoengineering to modify mean and extreme temperature, as well as the relative humidity and wet bulb temperature (WBT) change over Indonesian Maritime Continent (IMC) based on simulations from three different earth system models. We applied statistical downscaling and corrected the bias of model output to reproduce historical temperatures and relative humidity over IMC. We compared SRM as specified by the Geoengineering Model Intercomparison Project (GeoMIP) G4 experiment – 5 Tg of SO 2 into the lower stratosphere between 2020 and 2069 – concurrent with the RCP4.5 emissions scenario. SRM is able to significantly reduce the temperature means and extremes, and although differences in magnitude of response and spatial pattern occur, there is a generally consistent response. Additionally, the spatial response of changes forced by RCP4.5 scenario and SRM impact is notably heterogeneous in the archipelago, highlighting uncertainties that would be critical in assessing socio-economic consequences of both doing, and not doing SRM. In general, SRM has bigger impacts in reducing temperatures over land than oceans, and the southern monsoon region shows more variability. SRM is also effective at reducing the likelihood of WBT > 27OC events compared to RCP4.5 after some years of SRM deployment as well as during the post-termination period of SRM. Regional downscaling may be an effective tool in obtaining policy-relevant information about local effects of different future scenarios involving SRM.

54 ENVIRONMENTAL SCIENCES↗

Comparing multi-source urban flood indicators: satellite, simulation, and citizen-reported data

Urban flooding arises from complex mechanisms, making it challenging to capture accurately with a single detection method. This study evaluates three complementary approaches to detect flooding across three Chicago neighborhoods: (i) Sentinel-1 synthetic aperture radar (SAR), offering weather-independent, high-resolution (10 m) imagery of surface inundation; (ii) the storm water management model (SWMM), simulating combined sewer overflow and drainage performance; and (iii) citizen-generated 311 service requests, capturing observed flooding impacts. By analyzing six storms ranging from severe to mild, we examine how each source uniquely contributes to identifying urban flood events. SAR imagery effectively identifies standing water but can miss brief flooding due to satellite revisit constraints. SWMM provides detailed insights into system-wide drainage behavior yet may underestimate localized street-level flooding. Meanwhile, 311 calls reflect real-world flooding impacts but are vulnerable to underreporting. Statistical overlap analysis highlights chronic flood hotspots repeatedly identified across multiple detection methods, indicating persistent infrastructure and topographic vulnerabilities. Temporal analysis further reveals that while SWMM flooding aligns closely with rainfall peaks, 311 calls typically precede or persist beyond these peaks. Our findings emphasize the value of using satellite observations, hydrological modeling, and resident-reported data in a complementary manner to better interpret patterns in flood timing, severity, and spatial distribution—providing insights that can inform targeted infrastructure improvements and contribute to urban flood resilience planning.

311↗

Sustainability Trait Modeling of Field-Grown Switchgrass ( Panicum virgatum ) Using UAV-Based Imagery

Unmanned aerial vehicles (UAVs) provide an intermediate scale of spatial and spectral data collection that yields increased accuracy and consistency in data collection for morphological and physiological traits than satellites and expanded flexibility and high-throughput compared to ground-based data collection. In this study, we used UAV-based remote sensing for automated phenotyping of field-grown switchgrass (Panicum virgatum), a leading bioenergy feedstock. Using vegetation indices calculated from a UAV-based multispectral camera, statistical models were developed for rust disease caused by Puccinia novopanici, leaf chlorophyll, nitrogen, and lignin contents. For the first time, UAV remote sensing technology was used to explore the potentials for multiple traits associated with sustainable production of switchgrass, and one statistical model was developed for each individual trait based on the statistical correlation between vegetation indices and the corresponding trait. Also, for the first time, lignin content was estimated in switchgrass shoots via UAV-based multispectral image analysis and statistical analysis. The UAV-based models were verified by ground-truthing via correlation analysis between the traits measured manually on the ground-based with UAV-based data. The normalized difference red edge (NDRE) vegetation index outperformed the normalized difference vegetation index (NDVI) for rust disease and nitrogen content, while NDVI performed better than NDRE for chlorophyll and lignin content. Overall, linear models were sufficient for rust disease and chlorophyll analysis, but for nitrogen and lignin contents, nonlinear models achieved better results. As the first comprehensive study to model switchgrass sustainability traits from UAV-based remote sensing, these results suggest that this methodology can be utilized for switchgrass high-throughput phenotyping in the field.

60 APPLIED LIFE SCIENCES↗

Kaon Femtoscopy with Lévy-Stable Sources from $\sqrt{s_{NN}}$ = 200 GeV Au+Au Collisions at RHIC

Femtoscopy has the capacity to probe the space-time geometry of the particle-emitting source in heavy-ion collisions. In particular, femtoscopy of like-sign kaon pairs may shed light on the origin of non-Gaussianity of the spatial emission probability density. The momentum correlations between like-sign kaon pairs are measured in data recorded by the STAR experiment, from $\sqrt{s_{NN}}$ = 200 GeV Au+Au collisions at RHIC, BNL. Preliminary results hint at the possible existence of non-Gaussian, Lévy-stable sources, with the likely presence of an anomalous diffusion process in the signal for the identically charged kaon pairs so produced. More statistically significant studies at lower centre-of-mass energies may contribute to the search for the critical end point of QCD.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Simulation of ocean SAR images via phase history generation

A method for simulating an ocean synthetic-aperture radar (SAR) image is illustrated for a simple internal wave current pattern. The method calculates both the amplitude image and the radar signal history. The simulation model used consists of six stages. In the first stage, the full wave spectrum in two spatial coordinates is calculated from the wind speed and direction using the action spectral density equation. In the second stage, the pixel size is selected and used as the basis for dividing the spectrum into large and small scale motions. Realizations for the large-scale ocean surface height and velocity are then calculated. In the third stage, the sensor wavelength and geometry are used to calculate the small scale statistics (radar cross section, coherence time, root-mean square (RMS), radial velocity and RMS slope). In the fourth stage, an autoregressive method is used to generate a realization of the surface reflectivity history that is consistent with the radial velocity and the radial velocity variance. In the fifth stage, the signal history is generated by summing the reflectivities at the proper times with the appropriate antenna weighting function for the SAR. As a consequence of this process, speckle is automatically included in the signal. In the final stage, the SAR image can be created by using any of the traditional ways to process the signal history, including variable focusing and multilook processing.

Bennett, John R.↗

Supernova remnant rich fields in the Carina spiral arm

The analysis of the ROSAT PSPC data on five fields containing supernova remnants several months ago was completed. Dr. Una Hwang, prepared a paper describing our results which was published last August ('An X-ray Study of Five Supernova Remnants in the Carina Spiral Arm', by Hwang and Markert, 1994, Ap. J., 431, p. 819). Hwang's earlier analysis of this data became part of her PhD thesis ( 'X-ray Studies of Supernova Remnants', February 1994, MIT). A copy of the Hwang and Markert paper is appended. The results of the study are well-summarized in the Hwang and Markert paper: the best spatial-spectral X-ray study yet made of the intriguing SNR G296.1-0.7 was obtained. This study showed interesting spectral variations over the surface of the object.It was also determined the gross physical properties of G296 based on its X-ray emission. Four other fields were also examined. For three of these the upper limits to the radio object were determined, and in one case found a weak, but statistically significant X-ray object coincident with the peak of the radio flux.

Markert, Thomas↗

Spatial Distribution of Large Cloud Drops

By analyzing aircraft measurements of individual drop sizes in clouds, we have shown in a companion paper (Knyazikhin et al., 2004) that the probability of finding a drop of radius r at a linear scale l decreases as l(sup D(r)) where 0 less than or equal to D(r) less than or equal to 1. This paper shows striking examples of the spatial distribution of large cloud drops using models that simulate the observed power laws. In contrast to currently used models that assume homogeneity and therefore a Poisson distribution of cloud drops, these models show strong drop clustering, the more so the larger the drops. The degree of clustering is determined by the observed exponents D(r). The strong clustering of large drops arises naturally from the observed power-law statistics. This clustering has vital consequences for rain physics explaining how rain can form so fast. It also helps explain why remotely sensed cloud drop size is generally biased and why clouds absorb more sunlight than conventional radiative transfer models predict.

Marshak, A.↗

Passive/active Microwave Soil Moisture Change Disaggregation Using Smapvex12 Data

The SMAPVEX12 (Soil Moisture Active Passive (SMAP) Validation Experiment 2012) experiment was conducted during June-July 2012 in Manitoba, Canada with the goal of collecting remote sensing data and ground measurements for the development and testing of soil moisture retrieval algorithms under varying vegetation and soil conditions for the SMAP satellite. The aircraft based soil moisture data provided by the passive/active microwave sensor PALS (Passive and Active L-band System) has a nominal spatial resolution of 1600 m. However, this resolution is not compatible with agricultural, meteorological and hydrological studies that require high spatial resolutions and this issue can be solved by soil moisture disaggregation. The soil moisture disaggregation algorithm integrates radiometer soil moisture retrievals and high-resolution radar observations and it can provide soil moisture estimates at a finer scale than the radiometer data alone. In this study, a change detection algorithm was used for disaggregation of coarse resolution passive microwave soil moisture retrievals with radar backscatter coefficients obtained from the higher spatial resolution UAVSAR (Unmanned Air Vehicle Synthetic Aperture Radar) at crop field scale. The accuracy of the disaggregated change in soil moisture was evaluated using ground based soil moisture measurements collected during SMAPVEX12 campaign. The results showed that soil moisture spatial variabilities were better characterized by the disaggregated change in soil moisture estimates at 5 m / 800 m resolution as well as good agreement with in situ measurements. It also showed that VWC (Vegetation Water Content) did not have a big impact on disaggregation algorithm performance, with R2 of the disaggregated results ranging 0.628-0.794. The 5 m and 800m resolution disaggregated soil moisture did no show significant difference in statistical performance variables.

Radar Backscatter↗

Conversions between gas-phase metallicities in MaNGA

ABSTRACT We present polynomial conversions between each of 11 different strong line gas-phase metallicity calibrations, each based on ∼1.1 million star-forming spaxels in the public Sloan Digital Sky Survey (SDSS) Data Release 15 (DR15) Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. For this sample, which is ∼20 times larger than previous works, we present fifth-order polynomial fits for each of 110 possible calibration conversions, for both Small Magellanic Cloud-type and Milky Way-type dust corrections. The typical 2σ scatter around our polynomial fits is 0.1 dex; we present the range over which the metallicities are valid. Conversions between metallicities which rely on the same set of line ratios, or a heavily shared set of emission lines, have reduced scatter in their conversions relative to those conversions with little overlap in required emission lines. Calibration conversions with less consistent sets of emission lines also have increased galaxy-to-galaxy variability, and this variability can account for up to 35 per cent of the total scatter. We also compare our conversions to previous work with the single fibre SDSS DR7 spectra along with higher spatial resolution data from the TYPHOON Integral Field Spectroscopy survey, resulting in comparison samples with spatial resolutions from several kpc down to ∼100 pc. Our metallicity conversions, obtained with the large sample of MaNGA, are robust against the influence of diffuse ionized gas, redshift, effective radius, and spatial blurring, and are therefore consistent across both integrated spectra and the high-resolution integral field spectroscopy data.

79 ASTRONOMY AND ASTROPHYSICS↗

Integrated Modeling of U.S. Agricultural Soil Emissions of Reactive Nitrogen and Associated Impacts on Air Pollution, Health, and Climate

Agricultural soils are leading sources of reactive nitrogen (Nr) species including nitrogen oxides (NO x ), ammonia (NH 3 ), and nitrous oxide (N 2 O). Here, the propensity of NO x and NH 3 to generate ozone and fine particulate matter (PM) and associated impacts on health are highly variable, whereas the climate impacts of long-lived N 2 O are independent of emission timing and location. However, these impacts have rarely been compared on a spatially-resolved, monetized basis. In this study, we update the nitrogen scheme in an agroecosystem model to simulate the agricultural emissions of Nr from across the United States. We then apply a reduced-form air pollution health effects model to assess air quality impacts from NO x and NH 3 and a social cost of N 2 O to assess the climate impacts. Assuming an $\$$8.2 million Value of a Statistical Life (VSL) and a $\$$13,100/ton social cost of N 2 O, the air quality impacts are a factor of ~7 to 15 times as large as the climate impacts in heavily populated coastal regions, whereas the ratios are closer to 2.5 in sparsely populated regions. Our results show that air pollution, health, and climate should be considered jointly in future assessments of how farming practices affect Nr emissions.

54 ENVIRONMENTAL SCIENCES↗

Dynamic self-assembly of compartmentalized DNA nanotubes

Abstract Bottom-up synthetic biology aims to engineer artificial cells capable of responsive behaviors by using a minimal set of molecular components. An important challenge toward this goal is the development of programmable biomaterials that can provide active spatial organization in cell-sized compartments. Here, we demonstrate the dynamic self-assembly of nucleic acid (NA) nanotubes inside water-in-oil droplets. We develop methods to encapsulate and assemble different types of DNA nanotubes from programmable DNA monomers, and demonstrate temporal control of assembly via designed pathways of RNA production and degradation. We examine the dynamic response of encapsulated nanotube assembly and disassembly with the support of statistical analysis of droplet images. Our study provides a toolkit of methods and components to build increasingly complex and functional NA materials to mimic life-like functions in synthetic cells.

59 BASIC BIOLOGICAL SCIENCES↗

HFTS-2 Completions Design and State-of-the-Art Diagnostics Results

Various completion schemes were tested as part of the Hydraulic Fracture Test Site-2 (HFTS-2) project. The wells were instrumented with permanent fiber optic (FO) cable, and unique near-wellbore and far-field region surveys were acquired to evaluate spatial and temporal hydraulic fracture characteristics. The state-of-the-art diagnostics included substantial monitoring of both the well stimulation and longer-term production performance. The objective was to evaluate and improve stimulation distribution effectiveness (SDE) for different completion designs based on the diagnostics results. A hydraulic fracturing (HF) profile (i.e., proppant/fluid allocation among clusters) was obtained using near-wellbore Distributed Acoustic Sensing (DAS) data acquired during stimulation. Statistical analysis provided a quantitative assessment of SDE for the different completion designs. A novel workflow was developed to assess time-lapse fracture property changes at the cluster and stage levels based on Distributed Strain Sensing (DSS) data acquired during the production phase. These time-lapse fracture property changes were compared with the HF profile. Based on the statistical analysis, Rate Per Cluster (RPC) during proppant slurry placement appears to be a primary completion design variable, rather than the typical Rate Per Perf (RPP). If RPC was not maintained at a certain level, the likelihood of uneven distribution increased dramatically. The stage configurations that created low RPC while maintaining high RPP did not show improved SDE. The strain change results based on DSS data during the production phase indicated near-wellbore fracture property changes (e.g., aperture changes during well shut-in and flow periods) and aligned well with near-wellbore DAS signals acquired during stimulation. More importantly, fracture property changes based on DSS data acquired during two tests over a 7-month period show high consistency and good alignment with the near-wellbore HF profile. This workflow quantitatively evaluates SDE for different completion designs using diagnostic results and could help optimize completion strategy for future development projects. Here, the analysis and interpretation of results demonstrate the technical feasibility of using permanent fiber optic cable for stimulation and production monitoring and indic

58 GEOSCIENCES↗

CMIP6-based Multi-model Streamflow Projections over the Conterminous US

This dataset presents an ensemble of streamflow projections based on the hydroclimate projections dataset supporting the SECURE Water Act Section 9505 Assessment for the US Department of Energy (DOE) Water Power Technologies Office (WPTO). The six-member General Climate Model (GCM) ensemble from the Coupled Models Intercomparison Project phase 6 (CMIP6) downscaled using statistical (i.e., DBCCA) and dynamical (i.e., RegCM) and bias-corrected using two meteorological reference observations (Daymet & Livneh) are driven through two calibrated hydrologic models (VIC & PRMS) to simulate projected future hydrologic responses. This leads to the production of an ensemble of hydroclimate projections, including total runoff (surface runoff and baseflow) for 1980–2019 baseline and 2020–2059 near-term future periods under multiple emission scenarios at 1/24° (~4 km) spatial resolution across the CONUS. The total runoff projections are routed through the Routing Application for Parallel computatIon of Discharge (RAPID) routing model to produce an ensemble of streamflow projections for both periods across 2.7 million NHDPlusV2 stream reaches in the CONUS.

13 HYDRO ENERGY↗

Retrieval of Aerosol Optical Depth Under Thin Cirrus from MODIS: Application to an Ocean Algorithm

A strategy for retrieving aerosol optical depth (AOD) under conditions of thin cirrus coverage from the Moderate Resolution Imaging Spectroradiometer (MODIS) is presented. We adopt an empirical method that derives the cirrus contribution to measured reflectance in seven bands from the visible to shortwave infrared (0.47, 0.55, 0.65, 0.86, 1.24, 1.63, and 2.12 μm, commonly used for AOD retrievals) by using the correlations between the top-of-atmosphere (TOA) reflectance at 1.38 micron and these bands. The 1.38 micron band is used due to its strong absorption by water vapor and allows us to extract the contribution of cirrus clouds to TOA reflectance and create cirrus-corrected TOA reflectances in the seven bands of interest. These cirrus-corrected TOA reflectances are then used in the aerosol retrieval algorithm to determine cirrus-corrected AOD. The cirrus correction algorithm reduces the cirrus contamination in the AOD data as shown by a decrease in both magnitude and spatial variability of AOD over areas contaminated by thin cirrus. Comparisons of retrieved AOD against Aerosol Robotic Network observations at Nauru in the equatorial Pacific reveal that the cirrus correction procedure improves the data quality: the percentage of data within the expected error +/-(0.03 + 0.05 ×AOD) increases from 40% to 80% for cirrus-corrected points only and from 80% to 86% for all points (i.e., both corrected and uncorrected retrievals). Statistical comparisons with Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) retrievals are also carried out. A high correlation (R = 0.89) between the CALIOP cirrus optical depth and AOD correction magnitude suggests potential applicability of the cirrus correction procedure to other MODIS-like sensors.

cirrus correction↗

Satellite-Based Drought Reporting on the Navajo Nation

The Navajo Nation (NN) is the largest reservation in the US, and faces challenges related to water management during long-term and widespread drought episodes. The Navajo Nation is a federally recognized tribe, which has boundaries within Arizona, New Mexico, and Utah. The Navajo Nation has a land area of over 70,000 square kilometers. The Navajo Nation Department of Water Resources (NNDWR) reports on drought and climatic conditions through the use of regional Standardized Precipitation Index (SPI) values and a network of in-situ rainfall, streamflow, and climate data. However, these data sources lack the spatial detail and consistent measurements needed to provide a coherent understanding of the drought regime within the Nation's regional boundaries. This project, as part of NASA's Western Water Applications Office (WWAO), improves upon the recently developed Drought Severity Assessment Tool (DSAT) to ingest satellite-based precipitation data to generate SPI values for specific administrative boundaries within the reservation. The tool aims to: (1) generate SPI values and summary statistics for regions of interest on various timescales, (2) to visualize SPI values within a web-map application, and (3) produce maps and comparative statistical outputs in the format required for annual drought reporting. The co-development of the DSAT with NN partners is integral to increasing the sustained use of Earth Observations for water management applications. This tool will provide data to support the NN in allocation of drought contingency dollars to the regions most adversely impacted by declines in water availability.

Satellite-based↗

Visualization of hydraulic fracture using physics-informed clustering to process ultrasonic shear waves

Ultrasonic transmission is sensitive to the spatial variation in mechanical properties of materials due to the presence of cracks/fractures. Wave propagation through fractured media introduces changes in the frequency content, travel time and transmission coefficient of the wave. A workflow based on physics-informed unsupervised learning is developed to process the transmitted ultrasonic-shear waveforms to non-invasively visualize the geomechanical alterations due to hydraulic fracturing. Novelty of the work involves the assignment of both statistically consistent and physically consistent clusters to the measurements of shear waveforms acquired across the one axial and two frontal planes. Physically consistent/relevant information is incorporated by considering the travel time of the peak of spectral energy and transmission coefficient of the transmitted waveform. The proposed workflow generates maps of geomechanical alterations across the frontal and axial planes of the sample. The outputs of the workflow are in good agreement with independent techniques viz. acoustic emission and X-ray computed tomography. Finally, the proposed workflow can be adapted for improved fracture characterization in the subsurface when processing sonic-logging, cross-wellbore seismic or surface seismic waveform data.

42 ENGINEERING↗

Stream Temperature Prediction in a Shifting Environment: Explaining the Influence of Deep Learning Architecture

Stream temperature is a fundamental control on ecosystem health. Recent efforts incorporating process guidance into deep learning models for predicting stream temperature have been shown to outperform existing statistical and physical models. This performance is in part because deep learning architectures can actively learn spatiotemporal relationships that govern how water and energy propagate through a river network. However, exploration of how spatiotemporal awareness and process guidance influence a model's generalizability under shifting environmental conditions such as climate change is limited. Here, we use Explainable Artificial Intelligence (XAI) to interrogate how differing deep learning architectures affect a model's learned spatial and temporal dependencies, and how those learned dependencies affect a model's ability to maintain high accuracy when applied to unseen environmental conditions. Using the Delaware River Basin in the northeastern United States as a test case, we compare two spatiotemporally aware process–guided deep learning models for predicting stream temperature (a recurrent graph convolution network—RGCN, and a temporal convolution graph model—Graph WaveNet). Both models achieve equally high predictive performance when testing data are well represented in the training data (test root mean squared errors of 1.64°C and 1.65°C); however, Graph WaveNet significantly outperforms RGCN in 4 out of 5 experiments where test partitions represent different types of unseen environmental conditions. XAI results show that the architecture of Graph WaveNet leads to learned spatial relationships with greater fidelity to physical processes, and that this fidelity improves the generalizability of the model when applied to shifting and/or unseen environmental conditions.

54 ENVIRONMENTAL SCIENCES↗