Search NASA⌕ Search

SEARCH · Search NASA

Results for “Range limit”

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

Annihilation-limited long-range exciton transport in high-mobility conjugated copolymer films

A combination of ultrafast, long-range, and low-loss excitation energy transfer from the photoreceptor location to a functionally active site is essential for cost-effective polymeric semiconductors. Delocalized electronic wavefunctions along π-conjugated polymer (CP) backbone can enable efficient intrachain transport, while interchain transport is generally thought slow and lossy due to weak chain–chain interactions. In contrast to the conventional strategy of mitigating structural disorder, amorphous layers of rigid CPs, exemplified by highly planar poly(indacenodithiophene-co-benzothiadiazole) (IDT-BT) donor-accepter copolymer, exhibit trap-free transistor performance and charge-carrier mobilities similar to amorphous silicon. Here, we report long-range exciton transport in HJ -aggregated IDTBT thin-film, in which the competing exciton transport and exciton–exciton annihilation (EEA) dynamics are spectroscopically separated using a phase-cycling-based scheme and shown to depart from the classical diffusion-limited and strong-coupling regime. In the thin film, we find an annihilation-limited mechanism with ≪100% per-encounter annihilation probability, facilitating the minimization of EEA-induced excitation losses. In contrast, excitons on isolated IDTBT chains diffuse over 350 nm with 0.56 cm 2 s −1 diffusivity, before eventually annihilating with unit probability on first contact. We complement the pump–probe studies with temperature-dependent photocurrent and EEA measurements from 295 K to 77 K and find a remarkable correspondence of annihilation rate and photocurrent activation energies in the 140 K to 295 K temperature range.

Shi, Yuping (ORCID:0000000318289415)↗

Limits on the Range and Rate of Change in Power Take-Off Load in Ocean Wave Energy Conversion: A Study Using Model Predictive Control

Previous work comparing power take-off (PTO) architectures for ocean wave-powered reverse osmosis suggests that variable displacement in the wave energy converter (WEC)-driven pump does not offer a significant performance advantage. A limitation of that study is that the WEC was subject to a constant load within a given sea state (“Coulomb damping”) and did not account for controlled, moment-to-moment variation of the PTO load enabled by a variable displacement pump. This study explores the potential performance advantage of a variable PTO load over Coulomb damping. Model predictive control is used to provide optimal load control with constraints on the PTO load. The constraints include minimum and maximum loads and a limit on the rate of load adjustment. Parameter studies on these constraints enable conclusions about PTO design requirements in addition to providing an estimated performance advantage over Coulomb damping. Numerical simulation of the Oyster 1 WEC is carried out with performance weighted by historical sea state data from Humboldt Bay, CA. The results show a performance advantage of up to 20% higher yearly-average power absorption over Coulomb damping. Additionally, the parameter studies suggest that the PTO load should be adjustable down to at least 25% of the maximum load and should be adjustable between the minimum and maximum loads within a few seconds.

16 TIDAL AND WAVE POWER↗

Hybrid method for monitoring large Fabry-Pérot cavity displacements with nanometer precision

The change in length of an optical fiber-based Fabry-Pérot cavity (FPC) can be precisely measured using phase tracking, but the displacement range is limited by phase ambiguity. Period tracking techniques determine the absolute FPC length, but with larger uncertainties from tracking the spacing between multiple peaks. A hybrid method is demonstrated that identifies appropriate peaks for phase tracking using a coarse estimate obtained from the free spectral range to effectively maintain the high precision (∼1 nm) of phase tracking techniques to measure ∼24 µ m displacements, well beyond the range limitations (typically <1 µ m) of phase tracking methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Aerosol‐Jet‐Printable Covalent Organic Framework Colloidal Inks and Temperature‐Sensitive Nanocomposite Films

Abstract With molecularly well‐defined and tailorable 2D structures, covalent organic frameworks (COFs) have emerged as leading material candidates for chemical sensing, storage, separation, and catalysis. In these contexts, the ability to directly and deterministically print COFs into arbitrary geometries will enable rapid optimization and deployment. However, previous attempts to print COFs have been restricted by low spatial resolution and/or post‐deposition polymerization that limits the range of compatible COFs. Here, these limitations are overcome with a pre‐synthesized, solution‐processable colloidal ink that enables aerosol jet printing of COFs with micron‐scale resolution. The ink formulation utilizes the low‐volatility solvent benzonitrile, which is critical to obtaining homogeneous printed COF film morphologies. This ink formulation is also compatible with other colloidal nanomaterials, thus facilitating the integration of COFs into printable nanocomposite films. As a proof‐of‐concept, boronate‐ester COFs are integrated with carbon nanotubes (CNTs) to form printable COF‐CNT nanocomposite films, in which the CNTs enhance charge transport and temperature sensing performance, ultimately resulting in high‐sensitivity temperature sensors that show electrical conductivity variation by 4 orders of magnitude between room temperature and 300 °C. Overall, this work establishes a flexible platform for COF additive manufacturing that will accelerate the incorporation of COFs into technologically significant applications.

2D materials↗

Advancing neutron diffraction for accurate structural measurement of light elements at megabar pressures

Over the last 60 years, the diamond anvil cell (DAC) has emerged as the tool of choice in high pressure science because materials can be studied at megabar pressures using X-ray and spectroscopic probes. In contrast, the pressure range for neutron diffraction has been limited due to low neutron flux even at the strongest sources and the resulting large sample sizes. Here, we introduce a neutron DAC that enables break-out of the previously limited pressure range. Key elements are ball-bearing guides for improved mechanical stability, gem-quality synthetic diamonds with novel anvil support and improved in-seat collimation. We demonstrate a pressure record of 1.15 Mbar and crystallographic analysis at 1 Mbar on the example of nickel. Additionally, insights into the phase behavior of graphite to 0.5 Mbar are described. These technical and analytical developments will further allow structural studies on low-Z materials that are difficult to characterize by X-rays.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Probabilistic flux limiters

The stable numerical integration of shocks in compressible flow simulations relies on the reduction or elimination of Gibbs phenomena (unstable, spurious oscillations). A popular method to virtually eliminate Gibbs oscillations caused by numerical discretization in under-resolved simulations is to use a flux limiter. A wide range of flux limiters have been studied in the literature, with recent interest in their optimization via machine learning methods trained on high-resolution datasets. The common use of flux limiters in numerical codes as plug-and-play blackbox components makes them key targets for design improvement. Even for deterministic dynamical models, numerical uncertainty is introduced via coarse-graining required by insufficient computational power to solve all scales of motion. Conventional flux limiters are deterministic and lack the capacity to address uncertainties, both aleatoric (inherent randomness) and epistemic (modeling uncertainty due to limited knowledge), which arise in coarse-grained numerical simulations. Here, we introduce a conceptually distinct type of flux limiter that is designed to handle the effects of randomness in the model and uncertainty in model parameters. Unlike traditional single-function flux limiters, these new probabilistic flux limiters incorporate multiple flux limiting functions, each applied with a learned probability drawn from high-resolution data to mitigate the effects of uncertainty in numerical simulations. This approach departs from traditional single-function limiters by explicitly modeling and incorporating uncertainty into the shock capturing process. Using the example of Burgers' equation as a testbed, we show that a machine learned, probabilistic flux limiter may be used in a shock capturing code to more accurately capture shock profiles. In particular, we show that our probabilistic flux limiter outperforms standard limiters and can be successively improved upon (up to a point) by expanding the set of probabilistically chosen flux limiting functions.

97 MATHEMATICS AND COMPUTING↗

Evaluation of Electron Tomography Capabilities for Shale Imaging

Despite the advantageous resolution of electron tomography (ET), reconstruction of three-dimensional (3D) images from multiple two-dimensional (2D) projections presents several challenges, including small signal-to-noise ratios, and a limited projection range. This study evaluates the capabilities of ET for thin sections of shale, a complex nanoporous medium. A numerical phantom with 1.24 nm pixel size is constructed based on the tomographic reconstruction of a Barnett shale. A dataset of 2D projection images is numerically generated from the 3D phantom and studied over a range of conditions. First, common reconstruction techniques are used to reconstruct the shale structure. The reconstruction uncertainty is quantified by comparing overall values of storage and transport metrics, as well as the misclassification of pore voxels compared to the phantom. We then select the most robust reconstruction technique and we vary the acquisition conditions to quantify the effect of artifacts. We find a strong agreement for large pores over the different acquisition workflows, while a wider variability exists for nanometer-scale features. Furthermore, the limited projection range and reconstruction are identified as the main experimental bottlenecks, thereby suggesting that sample thinning, advanced holders, and advanced reconstruction algorithms offer opportunities for improvement.

3D scanning transmission electron microscopy↗

Room temperature colossal superparamagnetic order in aminoferrocene–graphene molecular magnets

Intensive studies are published for graphene-based molecular magnets due to their remarkable electric, thermal, and mechanical properties. However, to date, most of all produced molecular magnets are ligand based and subject to challenges regarding the stability of the ligand(s). The lack of long-range coupling limits high operating temperature and leads to a short-range magnetic order. Herein, we introduce an aminoferrocene-based graphene system with room temperature superparamagnetic behavior in the long-range magnetic order that exhibits colossal magnetocrystalline anisotropy of 8 × 10 5 and 3 × 10 7 J/m 3 in aminoferrocene and graphene-based aminoferrocene, respectively. These values are comparable to and even two orders of magnitude larger than pure iron metal. Aminoferrocene [C 10 H 11 FeN] + is synthesized by an electrophilic substitution reaction. It was then reacted with graphene oxide that was prepared by the modified Hammers method. The phase structure and functionalization of surface groups were characterized and confirmed by XRD, FT-IR, and Raman spectroscopy. To model the behavior of the aminoferrocene between two sheets of hydroxylated graphene, we have used density functional theory by placing the aminoferrocene molecule between two highly ordered hydroxylated sheets and allowing the structure to relax. Here, the strong bowing of the isolated graphene sheets suggests that the charge transfer and resulting magnetization could be strongly influenced by pressure effects. In contrast to strategies based on ligands surface attachment, our present work that uses interlayer intercalated aminoferrocene opens routes for future molecular magnets as well as the design of qubit arrays and quantum systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Comparing the influence of cation order and composition in simulated Zn(Sn, Ge)N 2 on structure, elastic moduli, and polarization for solid state lighting

Alloying and site ordering play complementary roles in dictating a material’s properties. However, deconvolving the impacts of these separate phenomena can be challenging. In this work, we simulate structures of Zn(Sn,Ge)N 2 with varied Sn content and site ordering to determine the impacts of order and composition on structural and electronic properties. We assess the formation enthalpy, lattice parameters, elastic constants, spontaneous polarization, and piezoelectric coefficients. In mostly disordered structures (order parameters ranging from 0.2 to 0.4), the formation enthalpy exhibits local extrema as a function of the order parameter, deviating from the more linear trends seen in both fully disordered and fully ordered systems. This anomalous deviation from the otherwise linear trend in formation enthalpy with order manifests in each of the other properties calculated. This range of order parameters of interest may be caused by a transition in the ordering of the quaternary material similar to phase changes seen in ternary compounds but stretched over a region inclduing 20% of the order parameter range. Most parameters calculated are more sensitive to order than to composition in the limited composition range tested; however, the lattice parameter c, piezoelectric coefficient e 33 , and elastic moduli C 12 , C 13 , and C 23 are more sensitive to composition. Of the properties compared, the piezoelectric coefficients are influenced most significantly by changes in both the composition and order parameter. Lattice parameters undergo the smallest changes with order and composition, but these small differences appear to impart large trends in the other properties. Better understanding the effects of disorder and group IV alloying in Zn(Sn,Ge)N 2 allows for more accurate modeling of characteristics of this material system for solid state lighting and other applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Imaging Wavepackets in Real-Space & Disentangling Ultrafast Solvation Dynamics with High-Energy Ultrafast X-ray Scattering

The microscopic information of solute and solute-solvent motions can be measured using ultrafast diffuse x-ray scattering, and compared to molecular dynamics simulations, as demonstrated in studies of solvated molecular systems done in X-ray free-electron lasers. However, for the typical photon energies used in these experiments, the photon momentum transfer range was limited to lower Q values where the signals arising from different parts of the studied system overlap significantly. Recently, this limitation was greatly alleviated by extending the photon energy at the Linac Coherent Light-Source (LCLS), significantly increasing the accessible momentum transfer range. This improvement is transformative for ultrafast diffuse x-ray scattering measurements, enabling for the first time to disentangle the solute scattering difference signal, that persists to higher Q values, from the bulk solvent and solute-solvent cross-terms difference signals. In addition, it enhances the prediction ability of modeling and simulation methods such as the hybrid QM/MM approach. The extended Q range also opens the way to resolve in real-space details regarding coherent wavepackets motions beyond their center-of-mass positions. Here, we present the first results on high-energy (18keV) diffuse x-ray scattering in solution, demonstrating high fidelity time-resolved scattering and analysis of the photoexcited model photocatalysts PtPOP ([Pt2(POP)4]4-), and IrDimen ([Ir2(dimen)4]2+) in several solvents. These complexes provide ideal systems for demonstrating the ability of high-Q ultrafast scattering and QM/MM simulations to decompose ultrafast solvation dynamics into specific changes in the solute-solvent pair distribution function with a particular focus on how electronically excited states change the interaction between the solvent and photo-catalytically active metal sites. We introduce an approach to further utilize the high-energy capability and demonstrate a single-shot ultra-wide-angle X-ray Scattering modality using two perpendicular detectors, spanning a scattering angle range of more than 100 degrees, allowing to extend the accessible momentum transfer range up to Q~14 Å-1. We develop a model-free method to invert the x-ray scattering signals and enable the recovery of multiple pair-density motions that happen simultaneously, allowing us to directly measure in real-space nuclear wavepacket motions. . [1] Natan, Adi. "Real-Space Inversion and Super-Resolution of Ultrafast X-ray Scattering using Natural Scattering Kernels." arXiv preprint arXiv:2107.05576 (2021)

Natan, Adi↗

An experimental investigation of flame and autoignition behavior of propane

We report autoignition delay time data are one important means to develop, quantify, and validate fundamental understanding of combustion chemistry at low temperatures (T<1200 K). However, low-temperature chemistry often has higher uncertainties and scatter in the experimental data compared with high-temperature ignition data (T>1200 K). In this study, autoignition properties of propane and oxygen mixtures were investigated using the University of Michigan rapid compression facility in order to understand the effects of ignition regimes on low-temperature ignition data. For the first time for propane, autoignition delay times were determined from pressure histories, and autoignition characteristics were simultaneously recorded using highspeed imaging of the test section through a transparent end-wall. Propane mixtures with fuel-to-O 2 equivalence ratios of φ = 0.25 and φ = 0.5 and O 2 -to-inert gas molar ratios of 1:3.76 were studied over the pressure range of 8.9 to 11.3 atm and the temperature range of 930 – 1070 K. The results showed homogeneous or strong autoignition occurred for all φ = 0.25 experiments, and inhomogeneous or mixed autoignition occurred for all φ = 0.5 experiments. While a limited temperature range is covered in the study, importantly the data span predicted transitions in autoignition behavior, allowing validation of autoignition regime hypotheses. Specifically, the results agree well with strong-autoignition limits proposed based on the Sankaran Criterion. The autoignition delay time data at the strong-ignition conditions are in excellent agreement with predictions using a well-validated detailed reaction mechanism from the literature and a zerodimensional modeling assumption. However, the experimental data at the mixed autoignition conditions were systematically faster than the model predictions, particularly at lower temperatures (T< ~970 K). The results are an important addition to the growing body of data in the literature that show mixed autoignition phenomena are important sources of the higher scatter observed in the low-temperature autoignition data for propane and other fuels. The results are discussed in terms of different methods to capture the effects of pre-autoignition heat release associated with mixed autoignition conditions and thereby address some of the discrepancies between kinetic modeling and experimental measurements.

42 ENGINEERING↗

Analysis of biofilm assembly by large area automated AFM

Biofilms are complex microbial communities critical in medical, industrial, and environmental contexts. Understanding their assembly, structure, genetic regulation, interspecies interactions, and environmental responses is key to developing effective control and mitigation strategies. While atomic force microscopy (AFM) offers critically important high-resolution insights on structural and functional properties at the cellular and even sub-cellular level, its limited scan range and labor-intensive nature restricts the ability to link these smaller scale features to the functional macroscale organization of the films. We begin to address this limitation by introducing an automated large area AFM approach capable of capturing high-resolution images over millimeter-scale areas, aided by machine learning for seamless image stitching, cell detection, and classification. Large area AFM is shown to provide a very detailed view of spatial heterogeneity and cellular morphology during the early stages of biofilm formation which were previously obscured. Using this approach, we examined the organization of Pantoea sp. YR343 on PFOTS-treated glass surfaces. Our findings reveal a preferred cellular orientation among surface-attached cells, forming a distinctive honeycomb pattern. Detailed mapping of flagella interactions suggests that flagellar coordination plays a role in biofilm assembly beyond initial attachment. Additionally, we use large-area AFM to characterize surface modifications on silicon substrates, observing a significant reduction in bacterial density. This highlights the potential of this method for studying surface modifications to better understand and control bacterial adhesion and biofilm formation.

59 BASIC BIOLOGICAL SCIENCES↗

Structure of Sulfuric Acid Solutions Using Pair Distribution Function Analysis

Solvation and mesoscale ordering of sulfuric acid and other strong acid solutions leads to suppressed freezing points and strong rheological changes with varying concentration. While the solid-state structures are well-understood, studies focused on the evolving solvation structure in the solution phase have probed a limited concentration range (~1–6 M). This study applies a total scattering approach in both the wide-angle X-ray scattering (WAXS) and pair distribution function (PDF) regimes to elucidate the evolving solvation structure over its full range of acid concentration (0–18 M). The emergence of a prepeak in the WAXS regime at intermediate concentrations indicates a transition from noninteracting sulfate molecules in the dilute limit to sterically limited sulfates at concentrations near its deep eutectic point. Fits to the PDF data quantify this trend, showing a transition from octahedrally hydrated sulfates up to 6–7 M concentrations, followed by gradual dehydration, and eventually reaching a solution structure similar to that of water-in-salt electrolyte systems at high acid concentrations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comparative characterization of mixed spectra and thermal neutron shielded irradiated tungsten

The effect of mixed spectra and thermal neutron shielded irradiation on tungsten was evaluated with plasma exposure in the tritium plasma experiment followed by thermal desorption spectroscopy, X-ray photoelectron spectroscopy, and transmission electron microscopy. The two different irradiation campaigns were performed at the High Flux Isotope Reactor to 0.39–0.74 displacement per atom (dpa) in the 894–1379 K temperature range. A neutron spectrum influence on the void size and void number density was not observed. However, a strong correlation was found between void size and void number density with temperature, but not with dpa in the limited dpa range of this study. Thermal neutron shielding significantly reduced the transmutation to Re+Os. Higher irradiation temperature will lead to larger voids with lower number density, which reduces deuterium retention. In conclusion, grain growth was also observed for high-temperature irradiation of over ~1300 K within the limited grains visible in the transmission electron microscopy specimens.

36 MATERIALS SCIENCE↗

Pythia8 Quark and Gluon Jets (float8 e4m3FN)

A float8 (e4m3FN) quantized version of the quark and gluon jet dataset originally published by Komiske, Metodiev, and Thaler (Zenodo record 3164691). Only the 20-file subset without charm and bottom quark jets is included here. All simulation parameters and jet selection criteria are identical to the original: Pythia 8.226, √s = 14 TeV Quarks from WeakBosonAndParton:qg2gmZq, gluons from WeakBosonAndParton:qqbar2gmZg with the Z decaying to neutrinos FastJet 3.3.0, anti-k_t jets with R = 0.4 p_T^jet ∈ [500, 550] GeV, |y^jet| < 1.7 There are 20 files, each in compressed NumPy format (QG_jets_fp8e4m3fn_0.npz through QG_jets_fp8e4m3fn_19.npz). Each file contains two arrays: X: (100000, M, 4) — 50k quark and 50k gluon jets, randomly sorted, padded to max multiplicity M, with particle features (pt, rapidity, azimuthal angle, pdgid) y: (100000,) — jet labels, gluon = 0, quark = 1 Since NumPy has no native fp8 dtype, X is stored as float32, but the values have been quantized through TensorFlow's float8_e4m3fn type and carry only fp8 precision. The quantization procedure is as follows: a global per-channel scale factor is computed from the absolute maximum value across all 20 chunks (with FP8_MAX = 448.0, the maximum representable value of e4m3FN). Each chunk is then scaled into the fp8 dynamic range, round-tripped through tf.experimental.float8_e4m3fn, and scaled back. This global scaling ensures a consistent quantization grid across the full dataset. The y labels are unchanged. Users should be aware that e4m3FN has limited dynamic range and precision. We recommend verifying this format is appropriate for your application; for a less aggressive reduction see the float16 and float32 versions linked below. If you use this dataset, please cite the original Zenodo record and its associated paper: Komiske, Metodiev, Thaler, Energy Flow Networks: Deep Sets for Particle Jets, JHEP 01 (2019) 121, arXiv:1810.05165

DiLullo, Nicholas [Brown University] (ORCID:000000↗

Factors That Influence Variability in Stress-Drop Measurements Using Spectral Decomposition and Spectral-Ratio Methods for the 2019 Ridgecrest Earthquake Sequence

Stress drop is a fundamental parameter related to earthquake source physics, but is hard to measure accurately. To better understand how different factors influence stress-drop measurements, we compare two different methods using the Ridgecrest stress-drop validation data set: spectral decomposition (SD) and spectral ratio (SR), each with different processing options. Here, we also examine the influence of spectral complexity on source parameter measurement. Applying the SD method, we find that frequency bandwidth and time-window length could influence spectral magnitude calibration, while depth-dependent attenuation is important to correctly map stress-drop variations. For the SR method, we find that the selected source model has limited influence on the measurements; however, the Boatwright model tends to produce smaller standard deviation and larger magnitude dependence than the Brune model. Variance reduction threshold, frequency bandwidth, and time-window length, if chosen within an appropriate parameter range, have limited influence on source parameter measurement. For both methods, wave type, attenuation correction, and spectral complexity strongly influence the result. The scale factor that quantifies the magnitude dependence of stress drop show large variations with different processing options, and earthquakes with complex source spectra deviating from the Brune-type source models tend to have larger scale factor than earthquakes without complexity. Based on these detailed comparisons, we make a few specific suggestions for data processing workflows that could help future studies of source parameters and interpretations.

58 GEOSCIENCES↗

Modeling the effect of habitat availability and quality on endangered winter-run Chinook salmon ( Oncorhynchus tshawytscha ) production in the Sacramento Valley

Reconciliation between water uses and sustaining aquatic species populations requires an effort to identify and quantify essential habitat characteristics for ecosystem health and valuation methods to predict ecosystem response to restoration actions. This process is particularly challenging for anadromous fish species such as California's Sacramento River winter-run Chinook salmon, due to their limited geographic range and diverse life history habitat requirements. Tools, such as life-cycle models, are needed to manage population dynamics and quantify the composite effects of processes across space and time. Nevertheless, complex institutions can hinder result interpretation and communication, and limit model use in decision-making. This paper focuses on the federally endangered and endemic Sacramento River winter-run Chinook (Oncorhynchus tshawytscha) by developing a Winter-Run Habitat-based Population Model (WRHAP). Additionally, WRHAP is a conceptual, freshwater rearing stage model that includes alternative rearing habitats reported in the literature (e.g., floodplains, off-channel and tributaries), defines rules of habitat use based on instream conditions, and incorporates a juvenile growth module that combines bioenergetics modeling with empirical growth rates. Model outputs reasonably follow observed out-migration patterns and provide a realistic smolt size distribution arriving at the San Francisco Bay. This effort demonstrates the importance of currently “non-critical” habitats (as defined by the Endangered Species Act) for juvenile development, with floodplain habitat contributing to a quarter of out-migrating biomass (despite < 18% availability and <10% of total rearing days), and off-channel growth being one of the most sensitive parameters (explaining ~13% of average juvenile weight variance). The model shows the utility of a simple population model to explore relationships between habitat quality/quantity and juvenile development, and to assist water/environmental management and decision-making processes focused on species recovery.

54 ENVIRONMENTAL SCIENCES↗