Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sub scale”

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,603 records · Page 89

A Diagnostic Test for Determining the Location of the GeV Emission in Powerful Blazars

An issue currently under debate in the literature is how far from the black hole is the Fermi-observed GeV emission of powerful blazars emitted. Here we present a clear diagnostic tool for testing whether the Ge V emission site is located within the sub-pc broad emission line (BLR) region or further out in the few pc scale molecular torus (MT) environment. Within the BLR the scatteri~ takes place at the onset of the Klein-Nishina regime, causing the electron cooling time to become almost energy independent and as a result, the variation of high-energy emission is expected to be achromatic. Contrarily, if the emission site is located outside the BLR, the expected GeY variability is energy-dependent and with amplitude increasing with energy. We demonstrate this using time-dependent numerical simulations of blazar variability.

Dotson, Amanda↗

Measuring the Redshift Dependence of The Cosmic Microwave Background Monopole Temperature With Planck Data

We study the capability of Planck data to constrain deviations of the cosmic microwave background (CMB) blackbody temperature from adiabatic evolution using the thermal Sunyaev-Zeldovich anisotropy induced by clusters of galaxies. We consider two types of data sets depending on how the cosmological signal is removed: using a CMB template or using the 217 GHz map. We apply two different statistical estimators, based on the ratio of temperature anisotropies at two different frequencies and on a fit to the spectral variation of the cluster signal with frequency. The ratio method is biased if CMB residuals with amplitude approximately 1 microK or larger are present in the data, while residuals are not so critical for the fit method. To test for systematics, we construct a template from clusters drawn from a hydro-simulation included in the pre-launch Planck Sky Model. We demonstrate that, using a proprietary catalog of X-ray-selected clusters with measured redshifts, electron densities, and X-ray temperatures, we can constrain deviations of adiabatic evolution, measured by the parameter α in the redshift scaling T (z) = T0(1 + z)(sup 1−alpha), with an accuracy of sigma(sub alpha) = 0.011 in the most optimal case and with sigma alpha = 0.018 for a less optimal case. These results represent a factor of 2-3 improvement over similar measurements carried out using quasar spectral lines and a factor 6-20 with respect to earlier results using smaller cluster samples.

cosmic background radiation↗

MMS Examination of FTEs at the Earth's Subsolar Magnetopause

Determining the magnetic field structure, electric currents, and plasma distributions within flux transfer event (FTE)-type flux ropes is critical to the understanding of their origin, evolution, and dynamics. Here the Magnetospheric Multiscale mission's high-resolution magnetic field and plasma measurements are used to identify FTEs in the vicinity of the subsolar magnetopause. The constant-alpha flux rope model is used to identify quasi-force free flux ropes and to infer the size, the core magnetic field strength, the magnetic flux content, and the spacecraft trajectories through these structures. Our statistical analysis determines a mean diameter of 1,700 ± 400 km (~30 ± 9 d(sub i)) and an average magnetic flux content of 100 ± 30 kWb for the quasi-force free FTEs at the Earth's subsolar magnetopause which are smaller than values reported by Cluster at high latitudes. These observed nonlinear size and magnetic flux content distributions of FTEs appear consistent with the plasmoid instability theory, which relies on the merging of neighboring, small-scale FTEs to generate larger structures. The ratio of the perpendicular to parallel components of current density, R(sub J), indicates that our FTEs are magnetically force-free, defined as R(sub J) < 1, in their core regions (<0.6 R(sub flux rope)). Plasma density is shown to be larger in smaller, newly formed FTEs and dropping with increasing FTE size. It is also shown that parallel ion velocity dominates inside FTEs with largest plasma density. Field-aligned flow facilitates the evacuation of plasma inside newly formed FTEs, while their core magnetic field strengthens with increasing FTE size.

Akhavan-Tafti, M.↗

Integrated transmission electron microscope

An integrated transmission electron microscope comprising multiple electron sources for tuned beams of ultrafast, scanning probe, and parallel illumination in varied beam energies can be alternated within sub-microseconds onto a sample with dynamic ‘transient state’ processes to acquire atomic-scale structural/chemical data with site specificity. The various electron sources and condenser optics enable high-resolution imaging, high-temporal resolution imaging, and chemical imaging, using fast-switching magnets to direct the different electron beams onto a single maneuverable objective pole piece where the sample resides. Such multimodal in situ characterization tools housed in a single microscope have the potential to revolutionize materials science.

Jungjohann, Katherine L.↗

Evaluation of Oak Ridge National Laboratory Health Physics Research Reactor Operation Data for Critical Benchmark Creation [Abstract]

The Oak Ridge National Laboratory (ORNL) Health Physics Research Reactor (HPRR) was a research reactor designed and built at ORNL in 1961. The critical assembly was using a highly enriched uranium and molybdenum alloy as the fuel, and it could be operated in steady-state or burst modes. The reactor was used for about 25 years to produce a lot of publications related to dosimetry, radiobiology and radiation detectors testing before its decommissioning in 1987. In recent years, the idea of using legacy operation data from the to create a valuable critical accident alarm system shielding benchmark arose. Such a benchmark has been submitted to the International Criticality Safety Benchmark Experiment Project (ICSBEP) Technical Review Group for a potential inclusion in the 2022 version of the handbook. Another way to use the valuable data from the operation of the HPRR is to evaluate the feasibility of the creation of a subcritical or prompt supercritical benchmark for inclusion in the ICSBEP or the International Reactor Physics Experiments Evaluation Project (IRPhEP) handbooks. To initiate a burst, the HPRR had to be operated in a slightly subcritical state for a few minutes. Then, the insertion of the burst control rod would greatly increase the reactivity of the system and start the burst. No critical configuration of the HPRR critical assembly could be located. The only information available concerns stable subcritical and prompt supercritical states, found in a burst experiments’ logbook. In the recovered logbook pages, information about 8 different bursts is available. The information includes the rods positions before and during a burst, the recorded subcritical reactor period and reactivity, and the burst fission yield derived from the temperature elevation sulfur pellet irradiation analysis. By using the HPRR logbook information and the as-built drawings of the critical assembly, a highly detailed model of the HPRR was created with SCALE 6.2.4/KENO-VI. Eight KENO-VI models were created to replicate the sub-critical assembly configurations described in the eight bursts from the recovered logbook pages. KENO-VI calculates k eff and it can be linked to a reactivity value in cents by using the delayed neutron fraction B eff , also calculated by KENO-VI. KENO-VI can also be used to model the prompt super-critical configurations of the HPRR and to assess the similarity with the burst measurements by comparing the calculated k eff and the measured fission yields between each burst. Unfortunately, high uncertainty exist and the obtained discrepancies between experiments and calculation results are high, compromising the creation of a valuable critical benchmark from HPRR operation data. The reasons of the discrepancies and potential ways to solve them are explored.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Quantitative comparison between sub-millisecond time resolution single-molecule FRET measurements and 10-second molecular simulations of a biosensor protein

Molecular Dynamics (MD) simulations seek to provide atomic-level insights into conformationally dynamic biological systems at experimentally relevant time resolutions, such as those afforded by single-molecule fluorescence measurements. However, limitations in the time scales of MD simulations and the time resolution of single-molecule measurements have challenged efforts to obtain overlapping temporal regimes required for close quantitative comparisons. Achieving such overlap has the potential to provide novel theories, hypotheses, and interpretations that can inform idealized experimental designs that maximize the detection of the desired reaction coordinate. Here, we report MD simulations at time scales overlapping with in vitro single-molecule Förster (fluorescence) resonance energy transfer (smFRET) measurements of the amino acid binding protein LIV-BP SS at sub-millisecond resolution. Computationally efficient all-atom structure-based simulations, calibrated against explicit solvent simulations, were employed for sampling multiple cycles of LIV-BP SS clamshell-like conformational changes on the time scale of seconds, examining the relationship between these events and those observed by smFRET. The MD simulations agree with the smFRET measurements and provide valuable information on local dynamics of fluorophores at their sites of attachment on LIV-BP SS and the correlations between fluorophore motions and large-scale conformational changes between LIV-BP SS domains. We further utilize the MD simulations to inform the interpretation of smFRET data, including Förster radius (R 0 ) and fluorophore orientation factor (κ 2 ) determinations. The approach we describe can be readily extended to distinct biochemical systems, allowing for the interpretation of any FRET system conjugated to protein or ribonucleoprotein complexes, including those with more conformational processes, as well as those implementing multi-color smFRET.

59 BASIC BIOLOGICAL SCIENCES↗

The influence of burn severity on dissolved organic carbon concentrations across a stream network differs based on seasonal wetness conditions

Large, high-severity wildfires in many regions across the globe have increased concerns about their impacts on carbon cycling in watersheds. Altered sources of carbon and changes in catchment hydrology after wildfire can lead to shifts in dissolved organic carbon (DOC) concentrations in streams, which can have negative impacts on aquatic ecosystem health and downstream drinking-water treatment. Despite its importance, post-fire DOC responses remain relatively unconstrained in the literature, and we lack critical knowledge of how burn severity, landscape elements, and climate interact to affect DOC concentrations. To improve our understanding of the impact of burn severity on DOC concentrations, we measured DOC at 129 sites across a stream network extending upstream, within, and downstream of a large, high-severity wildfire in Oregon, USA. We collected samples across the study sub-basin during four distinct seasonal wetness conditions. We used our high-spatial-resolution data to develop spatial stream network (SSN) models to predict DOC across the stream network and to improve our understanding of the controls on DOC concentrations. Spatially, we found no obvious wildfire signal – instead, we observed a pattern of increasing DOC concentrations from the high-elevation headwaters to the sub-basin outlet, while the mainstem maintained consistently low DOC concentrations. This suggests that effects from large wildfires may be “averaged” out at higher stream orders and larger spatial scales. When we grouped DOC concentrations by burn severity group, we observed a significant decrease in the variability of DOC concentrations in the moderate and high burn severity sub-catchments. However, our SSN models were able to predict decreases in DOC concentrations with increases in burn severity across the stream network. Decreases in DOC concentrations were also highly variable across seasonal wetness conditions, with the greatest (–1.40 to –1.64 mg L –1 ) decrease occurring in the high-severity group during the wetting season. Additionally, our models indicated that in all seasons, baseflow index was more influential in predicting DOC concentrations than burn severity was, indicating that groundwater discharge can obscure the impacts of wildfire in a stream network. Overall, our results suggested that landscape characteristics can regulate the DOC response to wildfire. Moreover, our results also indicated that the seasonal timing of sampling can influence the observed response of DOC concentrations to wildfire.

54 ENVIRONMENTAL SCIENCES↗

A fully Lagrangian, non-parametric bias model for dark matter halos

Here, we present a non-parametric Lagrangian biasing model and fit the ratio of the halo and mass densities at the field level using the mass-weighted halo field in the AbacusSummit simulations at z=0.5. Unlike the perturbative halo bias model that has been widely used in interpreting the observed large-scale structure traced by galaxies, we find a non-negative halo-to-mass ratio that increases monotonically with the linear overdensity δ 1 in the initial Lagrangian space. The bias expansion, however, does not guarantee non-negativity of the halo counts, and may lead to rising halo number counts at negative overdensities. The shape of the halo-to-mass ratio is unlikely to be described by a polynomial function of δ 1 and other quantities. Especially for massive halos with 6×10 12 h -1 M ⊙ , the halo-to-mass ratio starts soaring up at δ 1 >0, substantially different from the predictions of the bias expansion. We show that for the halo masses we consider (M>3×10 11 h -1 M ⊙ ) a non-parametric halo-to-mass ratio as a function of δ 1 and its local derivative ∇^2δ 1 can recover the halo power spectra to sub-percent level accuracy for wavenumbers k=0.01-0.1 h Mpc -1 given a proper smoothing scale to filter the initial density field, even though we do not fit the power spectrum directly. However, there is mild dependence of the recovery of the halo power spectrum on the smoothing scale and other input parameters. At k<0.01 h Mpc -1 and for massive halos with M>6x10 12 h -1 M ⊙ , our non-parametric model leads to a few percent overestimation of the halo power spectrum, indicating the need for larger or multiple smoothing scales. The halo-to-mass ratios obtained qualitatively agree with intuitions from extended Press-Schechter theory. We compare our framework to the bias expansion and discuss possible extensions.

79 ASTRONOMY AND ASTROPHYSICS↗

Groups and the Entropy Floor: XMM-Newton Observations of Two Groups

Using XMM-Newton spatially resolved X-ray imaging spectroscopy we obtain the temperature, density, entropy, gas mass, and total mass profiles for two groups of galaxies out to approximately 0.3 R(sub vir)(R(sub vir), the virial radius). Our density profiles agree well with those derived previously, and the temperature data are broadly consistent with previous results but are considerably more precise. Both of these groups are at the mass scale of 2x10(exp 13) M(solar mass), but have rather different properties. Both have considerably lower gas mass fractions at r < 0.3 R(sub vir), than the rich clusters. NGC2563, one of the least luminous groups for its X-ray temperature, has a very low gas mass fraction of approximately 0.004 inside 0.1 R(sub vir), which increases with radius. NGC4325, one of the most luminous groups at the same average temperature, has a higher gas mass fraction of 0.02. The entropy profiles and the absolute values of the entropy as a function of virial radius also differ, with NGC4325 having a value of approximately 100 keV cm(exp -2) and NGC2563 a value of approximately 300 keV cm(exp -2) at r approximately 0.1 R(sub vir). For both groups the profiles rise monotonically with radius and there is no sign of an entropy 'floor'. These results are inconsistent with pre-heating scenarios that have been developed to explain a possible entropy floor in groups, but are broadly consistent with models of structure formation that include the effects of heating and/or the cooling of the gas. The total entropy in these systems provides a strong constraint on all models of galaxy and group formation, and on the poorly defined feedback process that controls the transformation of gas into stars and thus the formation of structure in the universe.

Mushotzky, R. F.↗

A Multilevel Approach For SolvingLarge-Scale QUBO Problems With Noisy Hybrid Quantum Approximate Optimization

Quantum approximate optimization is one ofthe promising candidates for useful quantum computation,particularly in the context of finding approximate solutionsto Quadratic Unconstrained Binary Optimization (QUBO)problems. However, the existing quantum processing units(QPUs) are of relatively small size, and canonical mappingsof QUBO via the Ising model require one qubit per vari-able, rendering direct large-scale optimization infeasible.In classical optimization, a general strategy for addressingmany large-scale problems is via multilevel/multigrid meth-ods, where the large target problem is iteratively coarsenedand the global solution is constructed from multiple small-scale optimization runs. In this work, we experimentallytest how existing QPUs perform when used as a sub-solverwithin such a multilevel strategy. To this aim, we com-bine and extend (via additional classical processing steps)the recently proposed Noise-Directed Adaptive Remapping(NDAR) and Quantum Relax&Round (QRR) algorithms.We first demonstrate the effectiveness of our heuristicextensions on Rigetti’s superconducting transmon deviceAnkaa-2. We find approximate solutions to10instances offully connected82-qubit Sherrington-Kirkpatrick graphswith random integer-valued coefficients obtaining normal-ized approximation ratios (ARs) in the range∼0.98−1.0,and the same class with real-valued coefficients (ARs∼0.94−1.0). Then, we implement the extended NDAR andQRR algorithms as subsolvers in the multilevel algorithmfor6large-scale graphs with at most∼27,000variables.In practice, the QPU (with classical post-processing steps)is used to find approximate solutions to dozens of at most82-qubit problems, which are iteratively used to constructthe global solution. We observe that quantum optimizationresults are competitive in terms of the quality of solutionswhen compared to classical heuristics used as subsolverswithin the multilevel approach.Reproducibility: source code and data are available at[TBA upon acceptance]

quantum computing↗

Continuous fiber malleable thermoset composites with sub-1-minute dwell times; validation of impact performance and evaluation of the efficacy of the compression forming process. CRADA Final Report

Mallinda is developing polyimine malleable thermoset prepreg composite materials which have excellent mechanical properties (100 GPa tensile modulus, 2 GPA tensile strength, 2.4% elongation at break) and high operating temperatures (Tg>200°C). At scale, polyimine resins are commensurate in price with commodity epoxy resins. What distinguishes malleable thermoset prepreg from traditional thermoset prepreg materials, is that they are fully cured during Mallinda’s roll-to-roll production of prepreg laminate. This results in 5 key value-differentiating benefits. First, it simplifies manufacturing logistics by enabling ambient transportation and storage, and by significantly extending out-life and shelf-life almost indefinitely. Second, elimination of autoclave curing reduces the economic and energy costs to the customer. Third, scrap rates can be reduced as malleable thermoset prepreg materials are directly reusable. Fourth, the manufacturing consolidation step can be roughly 10x faster than traditional thermosets, because the resin is already cured. Parts can be made via compression forming by the application of heat and pressure to quickly vitrify and consolidate a multilayer part – easily leading to sub 3-minute cycle times (at lab scale we have demonstrated a 20 second dwell time, with room for further optimization). Finally, the closed-loop cradle-to-cradle solution-based recyclability of malleable thermoset composites can also contribute significantly to the future of sustainable lightweight materials. The focus of this project was the development, optimization and validation of malleable thermoset composite materials which exhibit manufacturing cycle times of 3-minute or less, high speed impact performance on par with incumbent technologies, and defect-free consolidation of 3 dimensional parts.

36 MATERIALS SCIENCE↗

Gate-All-Around Nanowire Field-Effect Transistors: A Historical Perspective

The development of transistor architectures, evolving from 2D planar metal-oxide-semiconductor field-effect transistors (MOSFETs) to FinFETs and then to gate-all-around nanowire (GAANW) FETs, plays a crucial role in downscaling technology nodes in the semiconductor industry. This perspective reviews the concept of MOSFETs and summarizes this historical development with particular emphasis on GAANW transistors due to their importance in next-generation technology for nodes below 3 nm. Specifically, the concept of GAANW transistors and their advantages over planar and FinFET devices for further scaling are presented, along with a discussion of their transition from early conceptual ideas to laboratory demonstrations and, ultimately, industrial adoption. Furthermore, potential solutions, such as complementary FETs (CFETs) and 2D semiconductor-based FETs, and their associated challenges for the future generation, known as the Angstrom Era, are discussed in a technological roadmap. In conclusion, this perspective may inspire both industry and academia in future research directions for realizing the Angstrom Era.

gate-all-around transistors↗

A petabyte size electronic library using the N-Gram memory engine

A model library containing petabytes of data is proposed by Triada, Ltd., Ann Arbor, Michigan. The library uses the newly patented N-Gram Memory Engine (Neurex), for storage, compression, and retrieval. Neurex splits data into two parts: a hierarchical network of associative memories that store 'information' from data and a permutation operator that preserves sequence. Neurex is expected to offer four advantages in mass storage systems. Neurex representations are dense, fully reversible, hence less expensive to store. Neurex becomes exponentially more stable with increasing data flow; thus its contents and the inverting algorithm may be mass produced for low cost distribution. Only a small permutation operator would be recalled from the library to recover data. Neurex may be enhanced to recall patterns using a partial pattern. Neurex nodes are measures of their pattern. Researchers might use nodes in statistical models to avoid costly sorting and counting procedures. Neurex subsumes a theory of learning and memory that the author believes extends information theory. Its first axiom is a symmetry principle: learning creates memory and memory evidences learning. The theory treats an information store that evolves from a null state to stationarity. A Neurex extracts information data without a priori knowledge; i.e., unlike neural networks, neither feedback nor training is required. The model consists of an energetically conservative field of uniformly distributed events with variable spatial and temporal scale, and an observer walking randomly through this field. A bank of band limited transducers (an 'eye'), each transducer in a bank being tuned to a sub-band, outputs signals upon registering events. Output signals are 'observed' by another transducer bank (a mid-brain), except the band limit of the second bank is narrower than the band limit of the first bank. The banks are arrayed as n 'levels' or 'time domains, td.' The banks are the hierarchical network (a cortex) and transducers are (associative) memories. A model Neurex was built and studied. Data were 50 MB to 10 GB samples of text, data base, and images: black/white, grey scale, and high resolution in several spectral bands. Memories at td, S(m(sub td)), were plotted against outputs of memories at td-1. S(m(sub td)) was Boltzman distributed, and memory frequencies exhibited self-organized criticality (SOC); i.e., 'l/f(sup beta)' after long exposures to data. Whereas output signals from level n may be encoded with B(sub output) = O(-log(2)f(sup beta)) bits, and input data encoded with B(sub input) = O((S(td)/S(td-1))(sup n)), B(sup output)/B(sub input) is much less than 1 always, the Neurex determines a canonical code for data and it is a lossless data compressor. Further tests are underway to confirm these results with more data types and larger samples.

Bugajski, Joseph M.↗

Sub-daily sea ice and deformation from RADARSAT observations

Periodic correlations in ice motion and deformation can be seen in lenght scales from 10 km and above, and suggest a 12 hour oscillation that is more likely associated with inertial rather than tidal frequencies.

oceanography ice mechanicscoriolis effects↗

A New Limit on Planck Scale Lorentz Violation from Gamma-ray Burst Polarization

Constraints on possible Lorentz invariance violation (UV) to first order in E/M(sub Plank) for photons in the framework of effective field theory (EFT) are discussed, taking cosmological factors into account. Then. using the reported detection of polarized soft gamma-ray emission from the gamma-ray burst GRB041219a that is indicative' of an absence of vacuum birefringence, together with a very recent improved method for estimating the redshift of the burst, we derive constraints on the dimension 5 Lorentz violating modification to the Lagrangian of an effective local QFT for QED. Our new constraints are more than five orders of magnitude better than recent constraints from observations of the Crab Nebula.. We obtain the upper limit on the Lorentz violating dimension 5 EFT parameter absolute value of zeta of 2.4 x 10(exp -15), corresponding to a constraint on the dimension 5 standard model extension parameter. Kappa (sup 5) (sub (v)oo) much less than 4.2 X 10(exp -3)4 / GeV.

Stecker, Floyd W.↗

Imaging local luminescence variations in CdSe 0.1 Te 0.9 thin films grown by the colossal grain growth process

Clear visualization and understanding of luminescence properties of grain interiors and grain boundaries in polycrystalline thin-film photovoltaic materials are crucial to achieving high-performance solar cells. Luminescence-based measurements, for example, reveal sub-bandgap energetic disorder or impurity-related radiative transitions in these materials, and when performed in the absolute scale, they can show the extent of voltage-limiting non-radiative recombination. Here we use a hyperspectral imaging technique to investigate inter- and intra-grain photoluminescence variations in polycrystalline CdSe 0.1 Te 0.9 films made by the colossal grain growth process. We observe that grain boundaries show higher luminescence than grain interiors for all energies below the bandgap energy. Multiple distinct sub-bandgap defect transitions are detected and determined to be related to donor-acceptor pair impurities. A comparison of the excitation-intensity-dependent photoluminescence at the grain boundaries with grain interiors suggests that carrier concentration is higher at grain boundaries. This finding is supported by photoemission electron microscopy imaging of the films, showing approximately 70 meV upward energy shift in the band structure at the grain boundaries. Finally, we investigate a dramatic thermal quenching behavior related to the donor-acceptor pairs luminescence over a temperature range and use a rate-equation model to attribute it to the sudden redirection of the radiative recombination channels to non-radiative channels with increasing temperature.

14 SOLAR ENERGY↗

Moisture sources for precipitation variability over the Arabian Peninsula

Abstract We apply the Lagrangian-based moisture back trajectory method to two reanalysis datasets to determine the moisture sources for wet season precipitation over the Arabian Peninsula, defined as land on the Asian continent to the south of the Turkish border and west of Iran. To accomplish this, we make use of the evaporative source region between 65°W–120°E and 30°S–60°N, which is divided into twelve sub-regions. Our comparison of reanalyses and multiple observations allows us to validate datasets and highlight broad-scale similarities in characteristics, notwithstanding some inconsistencies in the southwest AP. The results indicate north-to-south spatiotemporal heterogeneity in the characteristics of dominant moisture sources. In the north, moisture for precipitation is mainly sourced from midlatitude land and water bodies, such as the Mediterranean and Caspian Seas. Areas further south are dependent on moisture transport from the Western Indian Ocean and parts of the African continent. The El Niño-Southern Oscillation (ENSO) exhibits an overall positive but sub-seasonally varying influence on the precipitation variability over the region, with noticeable moisture anomalies from all major source regions. A significant drying trend exists over parts of the Peninsula, which both reanalyses partially attribute to anomalies in the moisture advection from the Congo Basin and South Atlantic Ocean. However, considerable uncertainty in evaporation trends over the terrestrial evaporative sources in observations warrants additional modeling studies to further our understanding of key processes contributing to the negative trends.

54 ENVIRONMENTAL SCIENCES↗