Search NASA⌕ Search

SEARCH · Search NASA

Results for “global radiation”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Aerosol-Cloud Interactions From Aviation Soot Emissions

Current models estimate global aviation contributes approximately 5% to the total anthropogenic climate forcing, with aerosol-cloud interactions having the greatest effect. However, radiative forcing estimates from aviation aerosol-cloud interactions remain undetermined. There is an expected significant increase in aircraft emissions with aviation demand expected to rise by over 4% per year. Soot may play an important role in the ice nucleation of aircraft-induced cirrus formation due to a high emission rate, but the ice nucleating properties are poorly constrained. Understanding the microphysical processes leading to atmospheric ice crystal formation is crucial for the reliable parameterization of aerosol-cloud interactions in climate models due to their impact on precipitation and cloud radiative properties. Ice nucleation of aircraft-emitted soot is potentially affected by particle morphology with condensation of supercooled water occurring in pores followed by ice nucleation. However, soot has heterogeneous properties and undergoes atmospheric aging and oxidation that could change surface properties and contribute to complex ice nucleation processes. Further, this review synthesizes current knowledge of ice nucleation catalyzed by aviation in the cirrus regime and its effects on global radiative forcing. Further research is required to determine the ice nucleation and microphysical processes of cirrus cloud formation from aviation emissions in both controlled laboratory and field investigations to inform models for more accurate climate predictions and to provide efficient mitigation strategies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Representing Anthropogenic Dust From Agricultural Sources in E3SMv1: Implementation, Evaluation, and Assessment of the Radiative Forcing

Dust emissions related to anthropogenic activities (i.e., anthropogenic dust) is not represented in most global climate models and its radiative impact remains unassessed. In this study, we develop a new and physically based method to parameterize anthropogenic dust emission from agricultural sources (i.e., agricultural dust, or AD) based on the DOE's Energy Exascale Earth System Model version 1 (E3SMv1). This method relates AD emission to the crop land use fraction in the E3SMv1 land component. Major AD sources simulated by our parameterization include those over Central America, the Sahel, North India, and North China. The annual averaged AD emission with diameter less than 10 μm is 567 Tg yr −1 in present-day (year 2000), which contributes to 13.3% of total dust emission. Model evaluation against satellite and ground-based observations shows that the new parameterization can represent AD emissions and global dust cycle reasonably well. We find that the total dust emission increases by 13.1% (495 Tg yr −1 ) from 1850 to 2000 due to the cropland land use fraction changes without considering the impact from climate changes. This induces a net dust direct effective radiative forcing of −0.041 W m −2 at top of the atmosphere. This dust-induced cooling exceeds 10% of the total anthropogenic aerosol direct effective radiative forcing from 1750 to 2014 estimated by the Intergovernmental Panel on Climate Change Sixth Assessment Report. Our findings indicate an important role of anthropogenic dust in global climate change, which should be included in future climate change assessments.

agricultural dust↗
DOE OSTI2025

Global Simulations of Gravitational Instability in Protostellar Disks with Full Radiation Transport. II. Locality of Gravitoturbulence, Clumpy Spirals, and Implications for Observable Substructure

Spiral perturbations in a gravitationally unstable accretion disk regulate disk evolution through angular momentum transport and heating and provide an observational signature of gravitational instability (GI). We use global 3D simulations to systematically characterize and understand these spiral perturbations. The spiral perturbations and the resulting transport are overall insensitive to the cooling type, with the exception that radiative cooling, especially in the optically thick regime, reduces the amplitude of temperature perturbations. Spiral perturbations are localized around corotation, allowing transport to be approximated by a local α viscosity to zeroth order in aspect ratio (H/R), but only after averaging over multiple orbits in time and/or multiple scale heights in space. Meanwhile, large-amplitude perturbations from strong gravitoturbulence can cause $\mathcal{O}$(α 1/2 ) deviation in the cooling rate of the disk. We develop empirical prescriptions for the angular momentum transport, heating, and cooling in a gravitoturbulent disk that capture the deviation from a viscous, unperturbed disk to first order in H/R and α 1/2 . The spiral perturbations in saturated gravitoturbulence are clumpy, with dense clumps forming through the nonlinear coupling between multiple modes at different m. Observationally, the clumpy gravitoturbulence produced by saturated GI can be mistaken with observational noise or embedded companions, especially under finite resolution. Meanwhile, grand-design spirals with m-fold symmetry may be uncommon among disks in saturated gravitoturbulence, and we speculate that they may instead be a signature of recently triggered or decaying GI.

Gravitational instability↗
DOE OSTI2023

Global organic and inorganic aerosol hygroscopicity and its effect on radiative forcing

The climate effects of atmospheric aerosol particles serving as cloud condensation nuclei (CCN) depend on chemical composition and hygroscopicity, which are highly variable on spatial and temporal scales. Here we present global CCN measurements, covering diverse environments from pristine to highly polluted conditions. We show that the effective aerosol hygroscopicity, κ, can be derived accurately from the fine aerosol mass fractions of organic particulate matter (ϵ org ) and inorganic ions (ϵ inorg ) through a linear combination, κ = ϵ org · κ org + ϵ inorg · κ inorg . In spite of the chemical complexity of organic matter, its hygroscopicity is well captured and represented by a global average value of κ org = 0.12 ± 0.02 with κinorg = 0.63 ± 0.01 as the corresponding value for inorganic ions. By showing that the sensitivity of global climate forcing to changes in κ org and κ inorg is small, we constrain a critically important aspect of global climate modelling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Terrestrial radiative cooling: Using the cold universe as a renewable and sustainable energy source

Photonic materials designed at wavelength scales have enabled a range of emerging energy technologies, from solid-state lighting to efficient photovoltaics that have transformed global energy landscapes. Daytime passive radiative cooling materials shed heat from the ground to the cold universe by taking advantage of the terrestrial thermal radiation that is as large as the renewable solar energy. Newly developed photonic materials permit subambient cooling under direct sunshine, and their applications are expanding rapidly enabled by scalable manufacturing. We review here the recent advancement of daytime subambient radiative cooling materials, which allow energy-efficient cooling and are paving the way toward technologies that harvest the coldness from the universe as a new renewable energy source.

Science & Technology - Other Topics↗
DOE OSTI2021

Scattering by whistler-mode waves during a quiet period perturbed by substorm activity

We study the dynamics of radiation belt electrons during a 10-day quiet period perturbed by substorm activity and preceding a high-speed stream (HSS), aiming at a global description of the radiation belts in L-shell, L in [2, 6], and energy [0.1, 10] MeV. We combine Van Allen Probes observations and Fokker-Planck numerical simulations of pitch-angle diffusion. The Fokker-Planck model uses event-driven pitch angle diffusion coefficients from whistler-mode waves, built from the wave properties and the ambient plasma density measurements from the Van Allen Probes. We first find this event has some similar characteristics to regular quiet times previously studied; a widely extended plasmasphere within which we observe strong and varying whistler-mode waves. These ambient conditions lead to strong pitch-angle scattering, which contributes to the creation of a wide slot region as well as a significant decay of the outer radiation belts, which are observed and qualitatively well simulated. In addition, we find the substorm activity causes short duration (within ± 4h) decay of the plasma density and a lowering amplitude of the whistler-mode waves within the plasmasphere, both causing opposite effects in terms of pitch angle diffusion. This leads to a diminution of pitch-angle diffusion at the time of the main substorm activity. Conversely, whistler-mode waves become enhanced in the time periods between the substorm injections. All effects cumulated, we find an enhancement of pitch angle diffusion by whistler-mode waves above L~4.7 during the 10-day period. This directly relates to the combination of quietness and substorm activity which allows pitch angle diffusing of up to 1 MeV electrons in the outer belt. Relativistic electrons of 1–2 MeV remain trapped in the outer belt, from L~4.7 to L~5.2, forming, in both the observations and the simulations, a distinct pocket of remnant electrons.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A Multivariate Space‐Time Dynamic Model for Characterizing the Atmospheric Impacts Following the Mt. Pinatubo Eruption

The June 1991 Mt. Pinatubo eruption resulted in a massive increase of sulfate aerosols in the atmosphere, absorbing radiation and leading to global changes in surface and stratospheric temperatures. A volcanic eruption of this magnitude serves as a natural analog for stratospheric aerosol injection, a proposed solar radiation modification method to combat a warming climate. The impacts of such an event are multifaceted and region-specific. Our goal is to characterize the multivariate and dynamic nature of the atmospheric impacts following the Mt. Pinatubo eruption. We developed a multivariate space-time dynamic linear model to understand the full extent of the spatially- and temporally-varying impacts. Specifically, spatial variation is modeled using a flexible set of basis functions for which the basis coefficients are allowed to vary in time through a vector autoregressive (VAR) structure. This novel model is cast in a Dynamic Linear Model (DLM) framework and estimated via a customized MCMC approach. We demonstrate how the model quantifies the relationships between key atmospheric parameters prior to and following the Mt. Pinatubo eruption with reanalysis data from MERRA-2 and highlight when such a model is advantageous over univariate models.

Dynamic Linear Model↗
DOE OSTI2025

Assessing the performance of solar radiation management geoengineering simulations

Offsetting the global warming caused by anthropogenic increases in atmospheric greenhouse gases by deliberate injection of aerosols into the stratosphere is the most studied of solar radiation management geoengineering schemes. The long-term success or failure of such schemes in achieving their stated goals is assessed by comparing simulated geoengineered temperature, precipitation and tropical cyclones metrics to equivalent fields in the simulated targeted climate simulations. Results using available data sets from three single model stabilized climate target experiments and three multimodel climate change reduction experiments are presented and compared against a measure of internal variability. While all but one experimental scheme is successful in achieving their targeted global mean annual surface temperature, their success at regional scales varies significantly and is often larger than the internal variability metric used here.

climate model evaluation↗
DOE OSTI2024

Derivations of the Total Radiation Belt Electron Content

Abstract We present multiple derivations of the Total Radiation Belt Electron Content (TRBEC), an indicator of the global number of electrons that instantaneously occupy the radiation belts. Derived from electron flux measurements, the TRBEC reduces the spatial information into a scalar quantity that concisely describes global aspects of the system. This index provides a simple, global, and long‐term assessment of the radiation belts that enables systematic analysis. In this work, we examine the TRBEC using the adiabatic invariants of which has been used in previous articles as this coordinate system removes reversible adiabatic effects. We then introduce a new expression to compute the TRBEC using the non‐adiabatic coordinates of , relevant in the contexts of energetic electron precipitation, chorus, and hiss scattering where adiabatic invariant quantities are no longer conserved. From both expressions of the TRBEC we demonstrate that an erroneous factor of that appeared in previous works using the adiabatic derivation led to an overestimate of the reported electron populations. In addition, we quantify electron loss in the outer radiation belt via a case study using the Van Allen Probes data over a 20‐day period from March 2013 specifying particle populations both in terms of the aforementioned adiabatic and non‐adiabatic variables. The total number of electrons in the outer radiation belt reached upwards of electrons at the peak of the storm, a rest mass of roughly 10 g.

Astronomy & Astrophysics↗
DOE OSTI2024

An analytic theory for the degree of Arctic Amplification

Arctic Amplification (AA), the amplified surface warming in the Arctic relative to the global mean, is a robust and impactful feature of climate change. While the basic physical picture of AA has been depicted, a clear understanding of how the degree of AA is determined has not been established. Here, by deciphering the intricate role of atmospheric heat transport (AHT), we build a two-box energy-balance model of AA and derive that the degree of AA is a simple nonlinear function of the Arctic and global feedbacks, the meridional heterogeneity in radiative forcing, and the partial sensitivities of AHT to global mean warming and meridional warming gradient. The formula captures the varying degree of AA in individual climate models and attributes the variation to specific physical factors. It further conveys a concise picture of how essential physics mutually determine the degree of AA and limit the range within 1.5~3.5. Our results articulate AHT as both forcing and feedback to AA, highlighting its partial sensitivities instead of total change as the key parameters for understanding AA. We also find that the effect of lapse rate feedback, a widely-recognized major contributor to AA, is fully offset by the effect of water vapor feedback.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The Global Atmosphere‐aerosol Model ICON‐A‐HAM2.3–Initial Model Evaluation and Effects of Radiation Balance Tuning on Aerosol Optical Thickness

Abstract The Hamburg Aerosol Module version 2.3 (HAM2.3) from the ECHAM6.3‐HAM2.3 global atmosphere‐aerosol model is coupled to the recently developed icosahedral nonhydrostatic ICON‐A (icon‐aes‐1.3.00) global atmosphere model to yield the new ICON‐A‐HAM2.3 atmosphere‐aerosol model. The ICON‐A and ECHAM6.3 host models use different dynamical cores, parameterizations of vertical mixing due to sub‐grid scale turbulence, and parameter settings for radiation balance tuning. Here, we study the role of the different host models for simulated aerosol optical thickness (AOT) and evaluate impacts of using HAM2.3 and the ECHAM6‐HAM2.3 two‐moment cloud microphysics scheme on several meteorological variables. Sensitivity runs show that a positive AOT bias over the subtropical oceans is remedied in ICON‐A‐HAM2.3 because of a different default setting of a parameter in the moist convection parameterization of the host models. The global mean AOT is biased low compared to MODIS satellite instrument retrievals in ICON‐A‐HAM2.3 and ECHAM6.3‐HAM2.3, but the bias is larger in ICON‐A‐HAM2.3 because negative AOT biases over the Amazon, the African rain forest, and the northern Indian Ocean are no longer compensated by high biases over the sub‐tropical oceans. ICON‐A‐HAM2.3 shows a moderate improvement with respect to AOT observations at AERONET sites. A multivariable bias score combining biases of several meteorological variables into a single number is larger in ICON‐A‐HAM2.3 compared to standard ICON‐A and standard ECHAM6.3. In the tropics, this multivariable bias is of similar magnitude in ICON‐A‐HAM2.3 and in ECHAM6.3‐HAM2.3. In the extra‐tropics, a smaller multivariable bias is found for ICON‐A‐HAM2.3 than for ECHAM6.3‐HAM2.3.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Why Does Atmospheric Radiative Heating Weaken Midlatitude Cyclones?

Abstract Recent work has indicated that atmospheric radiative heating reduces the kinetic energy of large‐scale eddies in the midlatitudes. However, a physical mechanism that connects radiation to the midlatitude eddy kinetic energy is still uncertain. Using a high‐resolution general circulation model we perform an experiment in which the radiative cooling profile at each model time step is overwritten with the climatological mean, computed from a control simulation. This approach separates the mean and transient effects of radiative heating on the extratropical circulation. We find that, when radiative heating is fixed, the globally‐averaged eddy kinetic energy is enhanced by ∼6%. We show that thermal radiation dampens temperature anomalies near the surface and tropopause in low‐pressure systems, destroying eddy available potential energy and eddy kinetic energy. We identify this as a possible mechanism by which atmospheric radiative heating weakens midlatitude cyclones.

Mischell, Eric↗
DOE OSTI2023

A global daily soil moisture dataset derived from Chinese FengYun Microwave Radiation Imager (MWRI)(2010–2019)

Surface soil moisture (SSM) is an important variable in drought monitoring, floods predicting, weather forecasting, etc. and plays a critical role in water and heat exchanges between land and atmosphere. SSM products from L-band observations, such as the Soil Moisture Active Passive (SMAP) Mission, have proven to be optimal global estimations. Although X-band has a lower sensitivity to soil moisture than that of L-band, Chinese FengYun-3 series satellites (FY-3A/B/C/D) have provided sustainable and daily multiple SSM products from X-band since 2008. This research developed a new global SSM product (NNsm-FY) from FY-3B MWRI from 2010 to 2019, transferred high accuracy of SMAP L-band to FY-3B X-band. The NNsm-FY shows good agreement with in-situ observations and SMAP product and has a higher accuracy than that of official FY-3B product. With this new dataset, Chinese FY-3 satellites may play a larger role and provide opportunities of sustainable and longer-term soil moisture data record for hydrological study.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Quantifying the Role of EMIC Wave Scattering During the 27 February 2014 Storm by RAM–SCB Simulations

Electromagnetic Ion Cyclotron (EMIC) wave scattering has been proved to be responsible for the fast loss of both radiation belt (RB) electrons and ring current (RC) protons. However, its role in the concurrent dropout of these two co–located populations remains to be quantified. In this work, we study the effect of EMIC wave scattering on both populations during the 27 February 2014 storm by employing the global physics–based RAM–SCB model. Throughout this storm event, MeV RB electrons and 100s keV RC protons experienced simultaneous dropout following the occurrence of intense EMIC waves. By implementing data–driven initial and boundary conditions, we perform simulations for both populations through the interplay with EMIC waves and compare them against Van Allen Probes observations. Notably, the results indicate that by including EMIC wave scattering loss, especially by the He–band EMIC waves, the model aligns closely with data for both populations. Additionally, we investigate the simulated pitch angle distributions (PADs) for both populations. Including EMIC wave scattering in our model predicts a 90° peaked PAD for electrons with stronger losses at lower pitch angles, while protons exhibit an isotropic PAD with enhanced losses at pitch angles above 40°. Furthermore, our model predicts considerable precipitation of both particle populations, predominantly confined to the afternoon to midnight sector (12 hr < MLT < 24 hr) during the storm's main phase, corresponding closely with the presence of EMIC waves.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Harnessing stratospheric diffusion barriers for enhanced climate geoengineering

Stratospheric sulfate aerosol geoengineering is a proposed method to temporarily intervene in the climate system to increase the reflectance of shortwave radiation and reduce mean global temperature. In previous climate modeling studies, choosing injection locations for geoengineering aerosols has, thus far, only utilized the average dynamics of stratospheric wind fields instead of accounting for the essential role of time-varying material transport barriers in turbulent atmospheric flows. Here we conduct the first analysis of sulfate aerosol dispersion in the stratosphere, comparing what is now a standard fixed-injection scheme with time-varying injection locations that harness short-term stratospheric diffusion barriers. We show how diffusive transport barriers can quickly be identified, and we provide an automated injection location selection algorithm using short forecast and reanalysis data. Within the first 7 d days of transport, the dynamics-based approach is able to produce particle distributions with greater global coverage than fixed-site methods with fewer injections. Additionally, this enhanced dispersion slows aerosol microphysical growth and can reduce the effective radii of aerosols up to 200–300 d after injection. While the long-term dynamics of aerosol dispersion are accurately predicted with transport barriers calculated from short forecasts, the long-term influence on radiative forcing is more difficult to predict and warrants deeper investigation. Statistically significant changes in radiative forcing at timescales beyond the forecasting window showed mixed results, potentially increasing or decreasing forcing after 1 year when compared to fixed injections. We conclude that future feasibility studies of geoengineering should consider the cooling benefits possible by strategically injecting sulfate aerosols at optimized time-varying locations. Our method of utilizing time-varying attracting and repelling structures shows great promise for identifying optimal dispersion locations, and radiative forcing impacts can be improved by considering additional meteorological variables.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The Role of Radiative Interactions in Tropical Cyclone Development under Realistic Boundary Conditions

Abstract The impact of radiative interactions on tropical cyclone (TC) climatology is investigated using a global, TC-permitting general circulation model (GCM) with realistic boundary conditions. In this model, synoptic-scale radiative interactions are suppressed by overwriting the model-generated atmospheric radiative cooling rates with their monthly varying climatological values. When radiative interactions are suppressed, the global TC frequency is significantly reduced, indicating that radiative interactions are a critical component of TC development even in the presence of spatially varying boundary conditions. The reduced TC activity is primarily due to a decrease in the frequency of pre-TC synoptic disturbances (“seeds”), whereas the likelihood that the seeds undergo cyclogenesis is less affected. When radiative interactions are suppressed, TC genesis shifts toward coastal regions, whereas TC lysis locations stay almost unchanged; together the distance between genesis and lysis is shortened, reducing TC duration. In a warmer climate, the magnitude of TC reduction from suppressing radiative interactions is diminished due to the larger contribution from latent heat release with increased sea surface temperatures. These results highlight the importance of radiative interactions in modulating the frequency and duration of TCs.

Meteorology & Atmospheric Sciences↗
DOE OSTI2020

Using CESM-RESFire to understand climate–fire–ecosystem interactions and the implications for decadal climate variability

Large wildfires exert strong disturbance on regional and global climate systems and ecosystems by perturbing radiative forcing as well as the carbon and water balance between the atmosphere and land surface, while short- and long-term variations in fire weather, terrestrial ecosystems, and human activity modulate fire intensity and reshape fire regimes. The complex climate–fire–ecosystem interactions were not fully integrated in previous climate model studies, and the resulting effects on the projections of future climate change are not well understood. Here we use the fully interactive REgion-Specific ecosystem feedback Fire model (RESFire) that was developed in the Community Earth System Model (CESM) to investigate these interactions and their impacts on climate systems and fire activity. We designed two sets of decadal simulations using CESM-RESFire for present-day (2001–2010) and future (2051–2060) scenarios, respectively, and conducted a series of sensitivity experiments to assess the effects of individual feedback pathways among climate, fire, and ecosystems. Our implementation of RESFire, which includes online land–atmosphere coupling of fire emissions and fire-induced land cover change (LCC), reproduces the observed aerosol optical depth (AOD) from space-based Moderate Resolution Imaging Spectroradiometer (MODIS) satellite products and ground-based AErosol RObotic NETwork (AERONET) data; it agrees well with carbon budget benchmarks from previous studies. We estimate the global averaged net radiative effect of both fire aerosols and fire-induced LCC at –0.59±0.52 W m –2 , which is dominated by fire aerosol–cloud interactions (–0.82±0.19 W m –2 ), in the present-day scenario under climatological conditions of the 2000s. The fire-related net cooling effect increases by ~170 % to –1.60±0.27 W m –2 in the 2050s under the conditions of the Representative Concentration Pathway 4.5 (RCP4.5) scenario. Such considerably enhanced radiative effect is attributed to the largely increased global burned area (+19 %) and fire carbon emissions (+100 %) from the 2000s to the 2050s driven by climate change. The net ecosystem exchange (NEE) of carbon between the land and atmosphere components in the simulations increases by 33 % accordingly, implying that biomass burning is an increasing carbon source at short-term timescales in the future. High-latitude regions with prevalent peatlands would be more vulnerable to increased fire threats due to climate change, and the increase in fire aerosols could counter the projected decrease in anthropogenic aerosols due to air pollution control policies in many regions. We also evaluate two distinct feedback mechanisms that are associated with fire aerosols and fire-induced LCC, respectively. On a global scale, the first mechanism imposes positive feedbacks to fire activity through enhanced droughts with suppressed precipitation by fire aerosol–cloud interactions, while the second one manifests as negative feedbacks due to reduced fuel loads by fire consumption and post-fire tree mortality and recovery processes. These two feedback pathways with opposite effects compete at regional to global scales and increase the complexity of climate–fire–ecosystem interactions and their climatic impacts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Coupled Climate Responses to Recent Australian Wildfire and COVID-19 Emissions Anomalies Estimated in CESM2

In this work, multiple 50-member ensemble simulations with the Community Earth System Model version 2 are performed to estimate the coupled climate responses to the 2019–2020 Australian wildfires and COVID-19 pandemic policies. The climate response to the pandemic is found to be weak generally, with global-mean net top-of-atmosphere radiative anomalies of +0.23 ± 0.14 W m –2 driving a gradual global warming of 0.05 ± 0.04 K by the end of 2022. While regional anomalies are detectable in aerosol burdens and clear-sky radiation, few significant anomalies exist in other fields due to internal variability. In contrast, the simulated response to Australian wildfires is a strong and rapid cooling, peaking globally at –0.95 ± 0.15 W m –2 in late 2019 with a global cooling of 0.06 ± 0.04 K by mid-2020. Transport of fire aerosols throughout the Southern Hemisphere increases albedo and drives a strong interhemispheric radiative contrast, with simulated responses that are consistent generally with those to a Southern Hemisphere volcanic eruption.

54 ENVIRONMENTAL SCIENCES↗