Search NASA⌕ Search

SEARCH · Search NASA

Results for “Atmospheric Emissions”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Aquaplanets as a Framework for Examination of Aerosol Effects

Although fundamental to the planetary radiative balance, aerosol impacts are highly uncertain in climate simulations because of the uneven distribution of aerosol sources and the complex interactions with radiation and clouds that are difficult to represent in climate models. This study proposes that aquaplanet configurations represent an idealized framework to investigate aerosol effects. As a simple demonstration, a series of aquaplanet simulations with the Community Atmosphere Model version 6 shows that the spatial distribution of aerosol emissions changes the aerosol effective radiative forcing even with unchanged total emissions. Some statistical properties of the simulations are presented to show that relatively short model integrations yield robust results. Much of the aerosol effect is shown to arise from aerosol–cloud interactions, especially through rapid adjustments associated with the aerosol lifetime effect that alter the cloud optical thickness.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

On the nature of droplet production in DC glows with a liquid anode: mechanisms and potential applications

The interactions between plasma and liquid solutions give rise to the formation of chemically reactive species useful for many applications, but the mass transport in the interfacial region is usually limited and not fully understood. In this work, we report on the observation and explanation of droplet ejection at the plasma–liquid interface of a one-atmosphere glow discharge with the liquid anode. The impact of droplets emission on plasma properties is also analyzed by spectroscopy. The process, which is an efficient mass and charge transport mechanism, apparently occurs during discharge operation and thus constitutes a feedback vehicle between the discharge and the liquid. Distinctive from the well-known Talyor cone droplets associated with liquid cathodes, the observed droplets originate from the bubbles due to electrolysis and solvated air which does not require strong electric field at liquid surface. Instead, the droplets are ejected by bubble cavity rupture at the plasma–liquid interface and their size, initial speed are strongly dependent on the gravity, inertia and capillarity. Further, the droplets emerge near the plasma attachment and are subsequently vaporized, emitting intense UV and visible light, which originated from excited OH radicals and sodium derived from the liquid electrolyte. Spectroscopy analysis confirmed that the bursting droplets generally reduce the gas temperature while their effects on electron density depend on the composition of the liquid anode. Results also show that droplets from NaCl solution increase the plasma electron density due to the lower ionization potential of sodium. These findings reveal a new mechanism for discharge maintenance and mass transport as well as suggest a simple approach to dispersing plasma-activated liquid into the gas phase and thus enhancing plasma–liquid interaction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The Science, Engineering, and Validation of Marine Carbon Dioxide Removal and Storage

Scenarios to stabilize global climate and meet international climate agreements require rapid reductions in human carbon dioxide (CO 2 ) emissions, often augmented by substantial carbon dioxide removal (CDR) from the atmosphere. While some ocean-based removal techniques show potential promise as part of a broader CDR and decarbonization portfolio, no marine approach is ready yet for deployment at scale because of gaps in both scientific and engineering knowledge. Marine CDR spans a wide range of biotic and abiotic methods, with both common and technique-specific limitations. Further targeted research is needed on CDR efficacy, permanence, and additionality as well as on robust validation methods—measurement, monitoring, reporting, and verification—that are essential to demonstrate the safe removal and long-term storage of CO 2 . Engineering studies are needed on constraints including scalability, costs, resource inputs, energy demands, and technical readiness. Research on possible co-benefits, ocean acidification effects, environmental and social impacts, and governance is also required.

climate mitigation↗
DOE OSTI2024

Strategic Siting of Direct Air Capture Facilities in the United States

Direct air capture (DAC) systems that capture carbon dioxide (CO 2 ) directly from the atmosphere are garnering considerable attention for their potential role as negative emission technologies in achieving net-zero CO 2 emission goals. Common DAC technologies are based either on liquid–solvent (L-DAC) or solid–sorbent (S-DAC) to capture CO 2 . A comprehensive multi-factor comparative economic analysis of the deployment of L-DAC and S-DAC facilities across various United States (U.S.) cities is presented in this paper. The analysis considers the influence of various factors on the favorability of DAC deployment, including local climatic conditions such as temperature, humidity, and CO 2 concentrations; the availability of energy sources to power the DAC system; and costs for the transport and storage of the captured CO 2 along with the consideration of the regional market and policy drivers. The deployment analysis in over 70 continental U.S. cities shows that L-DAC and S-DAC complement each other spatially, as their performance and operational costs vary in different climates. L-DAC is more suited to the hot, humid Southeast, while S-DAC is preferrable in the colder, drier Rocky Mountain region. Strategic deployment based on regional conditions and economics is essential for promoting the commercial adoptability of DAC, which is a critical technology to meet the CO 2 reduction targets.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Benchmark estimate of the effect of anthropogenic emissions on the ocean surface

Abstract Investigations into the role of anthropogenic emissions in the occurrence of extreme weather often use a method that compares simulations of atmospheric climate models run under a factual scenario of historical boundary conditions observed during the period of the event against simulations run under a counterfactual scenario of what those boundary conditions might naturally have been over that same period in the absence of anthropogenic emissions. A particular requirement for this experiment design is an accurate estimation of ocean surface boundary conditions for use by the counterfactual natural simulations. Here we use output from the CMIP5 multi‐climate‐model archive to develop a robust estimate of sea surface temperatures and sea ice conditions for use in counterfactual natural simulations, intended as a benchmark estimate to facilitate comparison across climate models and across studies. This development includes tests to ensure that the final estimate is stable from year‐to‐year and stable against other perturbations to the methodology, as well as consideration of the strengths and weaknesses in comparison to other available attributable warming estimates. While this estimate is tailored specifically for the International CLIVAR C20C+ Detection and Attribution Project, it can be used by related projects as well.

Stone, Dáithí A.↗
DOE OSTI2022

Particle emissions from mobile sources: Discussion of ultrafine particle emissions and definition

There is no universally agreed upon definition for ultrafine particles (UFP). Commonly used definitions for UFP are either particle number below 100 nm or total particle number, but without an agreed upon lower cut point. For example, a lower cut point of 3 nm compared to 10 nm could result in a substantially higher count. Another definition for UFP is total particle mass but without a commonly agreed upon aerodynamic diameter upper cut point, e.g., below 100 nm, 200 nm, 300 nm, etc. Yet another definition is lung deposited surface area weighted by lung deposition fraction, found mainly in the particle mobility diameter range from 20 to 400 nm. It is clear from these definitions that there are inconsistencies in the way UFP is used and defined in the literature. Sometimes these metrics are well correlated, sometimes not. In this paper we suggest three exposure metrics: UFP-N, UFP-M, and UFP-S, that we believe will add clarity. These metrics represent total number, mass, and surface area below 500 nm, respectively. For surface area and mass, the 500 nm cut point can be either aerodynamic or mobility diameter depending upon measurement methodology. For all metrics, this cut point captures nearly all of the primary particle emissions from mobile sources. Furthermore, UFP-N would include a lower cut point of 3–6 nm and would not require an upper size cut point because there is very little particle number above 500 nm or even above 100 nm. Thus, our definition of UFP-N is consistent with the current definition of ultrafine number except for, importantly, the specification of a lower cut point. These exposure metrics can help facilitate consistency in the characterization of both short- and long-term UFP ambient exposures and associated health effects in epidemiological studies.

42 ENGINEERING↗
DOE OSTI2021

Simulation of radon-222 with the GEOS-Chem global model: Emissions, seasonality, and convective transport

Radon-222 (222Rn) is a short-lived radioactive gas naturally emitted from land surface, and has long been used to assess convective transport in atmospheric models. In this study, we simulate 222Rn using the GEOS-Chem chemical transport model with aims to improve our understanding of 222Rn emissions and surface concentration seasonality, and to characterize convective transport associated with two Goddard Earth Observing System (GEOS) meteorological products, MERRA and GEOS-FP. We evaluate four available 222Rn emission scenarios by comparing model results with surface observations at 51 global surface sites. The default emission scenario in GEOS-Chem yields a moderate agreement with global surface observations (< 70% data within a factor of 2) and a large underestimate of wintertime surface 222Rn concentrations at Northern Hemisphere mid- and high-latitudes due to an oversimplified formulation of 222Rn emission fluxes (1 atom cm-2 s-1 over land with a reduction by a factor of 3 under freezing conditions). We compose a new global 222Rn emission scenario based on Zhang et al. (2011) and show its potential to improve simulated surface 222Rn concentrations and seasonality. The regional components of this emission scenario include spatially and temporally varying emission fluxes derived from previous measurements of soil radium content and soil exhalation models, which are key factors to 222Rn emission flux rates. However, large model underestimates of surface 222Rn concentrations still exist in Asia, suggesting unusually high regional 222Rn emissions. We propose a conservative up-scaling factor of 1.2 for 222Rn emission fluxes in China, as also constrained by the observed deposition fluxes of 210Pb (decay daughter of 222Rn). With this modification, the model shows better agreement with the observations in Europe and North America (>80% data within a factor of 2), and reasonable agreement in Asia (close to 70%). Further constraints on 222Rn emissions would require additional observations of surface 222Rn concentrations and emission fluxes in central U.S., Canada, Africa, and Asia. We also compare and assess convective transport in model simulations driven by MERRA and GEOS-FP using observed 222Rn vertical profiles during northern mid-latitude summertime and from two short-term airborne campaigns. While the simulations with both GEOS products are able to capture the observed vertical gradient of 222Rn concentrations in the lower troposphere (0-4 km), neither correctly represents the level of convective detrainment, resulting in biases in the middle and upper troposphere. Compared to GEOS-FP, MERRA leads to stronger convective transport of 222Rn, which is partially compensated by its weaker large-scale vertical advection, resulting in similar global vertical distributions of 222Rn concentrations between the two simulations.

Zhang, Bo↗
DOE OSTI2021

Estimation of Possible Primary Biological Particle Emissions and Rupture Events at the Southern Great Plains ARM Site

We use 10 years of data from the Department of Energy (DoE) Atmospheric Radiation Measurements (ARM) United States Southern Great Plains (SGP) site with nearby regional pollen and fungal spore measurements to indirectly estimate the seasonal influence of these two primary biological aerosol particles (PBAP). We estimate possible primary emissions of larger PBAP and PBAP rupture events, which form submicron organic aerosol during precipitation or high relative humidity. High pollen counts at two urban stations near SGP occur during late winter/early spring (day of year (DOY) 50–120) and late summer (DOY 240–310). Around 4–19 days per year show possible pollen events (PPE) when near-surface lidar observations of daily linear particle depolarization ratio >0.1 are coincident with high organic aerosol fraction. For PPE days with rainfall, aerosol size distribution observations show enhanced submicron particle concentrations consistent with pollen rupture events. For fungal spores, high fungal spore counts occur during late spring/early summer (DOY 110–195) and late summer/autumn (DOY 220–340). Based on size distribution observations, up to 7% of days have possible fungal spore rupture events (PFE) with higher aerosol number count specifically over the range expected for fungal spore fragment mobility diameter (20–50 nm). These short-lived PFE correlate with rainfall or occur after prolonged exposure to rainfall (e.g., >10 h). While the SGP site lacks direct measurements of bioaerosol and large particle sizes, this analysis suggests that PBAP primary emissions and rupture events could occur about 32 days per year, representing an important component of the aerosol budget during seasonal emissions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Novel gas exposure system for the controlled exposure of plants to gaseous hydrogen fluoride

Abstract Plants can serve as sensitive bioindicators of the presence of contaminant vapors in the atmosphere. This work describes a novel laboratory-based gas exposure system capable of calibrating plants as bioindicators for the detection and delineation of the atmospheric contaminant hydrogen fluoride (HF) as a preparatory step for monitoring release emissions. To evaluate changes in plant phenotype and stress-induced physiological effects attributed to HF alone, the gas exposure chamber must have additional controls to simulate otherwise optimal plant growth conditions including variables such as light intensity, photoperiod, temperature, and irrigation. The exposure system was designed to maintain constant growth conditions during a series of independent experiments that varied between optimal (control) and stressful (HF exposure) conditions. The system was also designed to ensure the safe handling and application of HF. An initial system calibration introduced HF gas into the exposure chamber and monitored HF concentrations by cavity ring-down spectroscopy for a 48-h period. Stable concentrations inside the exposure chamber were observed after approximately 15 h, and losses of HF to the system ranged from 88 to 91%. A model plant species ( Festuca arundinacea ) was then exposed to HF for 48 h. Visual phenotype stress-induced responses aligned with symptoms reported in the literature for fluoride exposure (tip dieback and discoloration along the dieback transition margin). Fluoride concentrations in exposed tissues compared to control tissues confirmed enhanced fluoride uptake due to HF exposure. The system described herein can be applied to other reactive atmospheric pollutants of interest in support of bioindicator research.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

COVID-19 pandemic underscores role of green space in urban carbon dynamics

For Detroit Michigan the arrival of COVID-19 led to intensive measures to prevent further spread of the virus resulting in consequent changes in traffic and energy use. Here, we take advantage of these different emission scenarios to explore CO 2 dynamics in a postindustrial city with a declining population and increasing green space. We present atmospheric CO 2 concentration and net urban ecosystem exchange of CO 2 (NUE) from a typical eddy covariance system and canopy greenness from a field camera on the Wayne State University campus in midtown Detroit. We categorized our study period (January 18, 2020–July 31, 2020) into three subperiods associated with the state-wide shelter-in-place order. Our results support that the city was a net carbon source throughout the period, particularly during the shelter-in-place period, although reduced traffic lowered CO 2 concentrations and NUE. However, during the post-order period when traffic was highest, atmospheric CO 2 concentrations and NUE were lowest, suggesting that the greening of urban vegetation may have greater carbon mitigation potential than lowering anthropogenic carbon emissions through traffic reductions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Important role of stratospheric injection height for the distribution and radiative forcing of smoke aerosol from the 2019–2020 Australian wildfires

More than 1 Tg smoke aerosol was emitted into the atmosphere by the exceptional 2019–2020 southeastern Australian wildfires. Triggered by the extreme fire heat, several deep pyroconvective events carried the smoke directly into the stratosphere. Once there, smoke aerosol remained airborne considerably longer than in lower atmospheric layers. The thick plumes traveled eastward, thereby being distributed across the high and mid-latitudes in the Southern Hemisphere, enhancing the atmospheric opacity. Due to the increased atmospheric lifetime of the smoke plume, its radiative effect increased compared to smoke that remains in lower altitudes. Global models describing aerosol-climate impacts lack adequate descriptions of the emission height of aerosols from intense wildfires. Here, we demonstrate, by a combination of aerosol-climate modeling and lidar observations, the importance of the representation of those high-altitude fire smoke layers for estimating the atmospheric energy budget. Through observation-based input into the simulations, the Australian wildfire emissions by pyroconvection are explicitly prescribed to the lower stratosphere in different scenarios. Based on our simulations, the 2019–2020 Australian fires caused a significant top-of-atmosphere (TOA) hemispheric instantaneous direct radiative forcing signal that reached a magnitude comparable to the radiative forcing induced by anthropogenic absorbing aerosol. Up to +0.50 W m -2 instantaneous direct radiative forcing was modeled at TOA, averaged for the Southern Hemisphere (+0.25 W m -2 globally) from January to March 2020 under all-sky conditions. At the surface, on the other hand, an instantaneous solar radiative forcing of up to -0.81 W m -2 was found for clear-sky conditions, with the respective estimates depending on the model configuration and subject to the model uncertainties in the smoke optical properties. Since extreme wildfires are expected to occur more frequently in the rapidly changing climate, our findings suggest that high-altitude wildfire plumes must be adequately considered in climate projections in order to obtain reasonable estimates of atmospheric energy budget changes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Life cycle greenhouse gas emissions of different CO 2 supply options for an algal biorefinery

This paper characterizes the greenhouse gas (GHG) emissions from the provision of CO 2 for twelve hours of the day from a range of sources operating on a twenty-four-hour basis. Twelve CO 2 supply cases spanning from fossil fuel power plant stack gases to direct air capture (DAC) systems are modeled. The results are compared in a life cycle assessment (LCA) framework. The CO 2 supply from the combustion or gasification of biomass has GHG emissions of –1.62 gCO 2 eq./g CO 2 delivered and is the case with the lowest carbon footprint for supplied CO2. The emissions for all the CO 2 supply cases were found to be negative, either because of capturing more CO 2 from the atmosphere than is released from the CO 2 production process, or because of production of excess electricity that is exported to the grid (replacing emissions from grid electricity). A purpose-built Natural Gas Combined Cycle (NGCC) power plant CO 2 source has the highest GHG emissions of -0.40 gCO 2 eq./gCO 2 delivered of the cases considered because of nighttime emissions and natural gas supply chain losses. The impact of a diurnal cycle can be mitigated by using capture and refrigeration systems as part of the CO 2 management system. The results are sensitive to the GHG emissions of the grid electricity, which is imported or exported in different cases. The results provide a benchmark for the comparison of different CO 2 supply options for the establishment of a sustainable algal biorefinery.

42 ENGINEERING↗
DOE OSTI2023

Modeled production, oxidation, and transport processes of wetland methane emissions in temperate, boreal, and Arctic regions

Abstract Wetlands are the largest natural source of methane (CH 4 ) to the atmosphere. The eddy covariance method provides robust measurements of net ecosystem exchange of CH 4 , but interpreting its spatiotemporal variations is challenging due to the co‐occurrence of CH 4 production, oxidation, and transport dynamics. Here, we estimate these three processes using a data‐model fusion approach across 25 wetlands in temperate, boreal, and Arctic regions. Our data‐constrained model—iPEACE—reasonably reproduced CH 4 emissions at 19 of the 25 sites with normalized root mean square error of 0.59, correlation coefficient of 0.82, and normalized standard deviation of 0.87. Among the three processes, CH 4 production appeared to be the most important process, followed by oxidation in explaining inter‐site variations in CH 4 emissions. Based on a sensitivity analysis, CH 4 emissions were generally more sensitive to decreased water table than to increased gross primary productivity or soil temperature. For periods with leaf area index (LAI) of ≥20% of its annual peak, plant‐mediated transport appeared to be the major pathway for CH 4 transport. Contributions from ebullition and diffusion were relatively high during low LAI (<20%) periods. The lag time between CH 4 production and CH 4 emissions tended to be short in fen sites (3 ± 2 days) and long in bog sites (13 ± 10 days). Based on a principal component analysis, we found that parameters for CH 4 production, plant‐mediated transport, and diffusion through water explained 77% of the variance in the parameters across the 19 sites, highlighting the importance of these parameters for predicting wetland CH 4 emissions across biomes. These processes and associated parameters for CH 4 emissions among and within the wetlands provide useful insights for interpreting observed net CH 4 fluxes, estimating sensitivities to biophysical variables, and modeling global CH 4 fluxes.

Ueyama, Masahito↗
DOE OSTI2020

Controls on Soil Organic Matter Degradation and Subsequent Greenhouse Gas Emissions Across a Permafrost Thaw Gradient in Northern Sweden

Warming-induced permafrost thaw could enhance microbial decomposition of previously stored soil organic matter (SOM) to carbon dioxide (CO 2 ) and methane (CH 4 ), one of the most significant potential feedbacks from terrestrial ecosystems to the atmosphere in a changing climate. The environmental parameters regulating microbe-organic matter interactions and greenhouse gas (GHG) emissions in northern permafrost peatlands are however still largely unknown. The objective of this work is to understand controls on SOM degradation and its impact on porewater GHG concentrations across the Stordalen Mire, a thawing peat plateau in Northern Sweden. Here, we applied high-resolution mass spectrometry to characterize SOM molecular composition in peat soil samples from the active layers of a Sphagnum-dominated bog and rich fen sites in the Mire. Microbe-organic matter interactions and porewater GHG concentrations across the thaw gradient were controlled by aboveground vegetation and soil pH. An increasingly high abundance of reduced organic compounds experiencing greater humification rates due to enhanced microbial activity were observed with increasing thaw, in parallel with higher CH 4 and CO 2 porewater concentrations. Bog SOM however contained more Sphagnum-derived phenolics, simple carbohydrates, and organic- acids. The low degradation of bog SOM by microbial communities, the enhanced SOM transformation by potentially abiotic mechanisms, and the accumulation of simple carbohydrates in the bog sites could be attributed in part to the low pH conditions of the system associated with Sphagnum mosses. We show that Gibbs free energy of C half reactions based on C oxidation state for OM can be used as a quantifiable measure for OM decomposability and quality to enhance current biogeochemical models to predict C decomposition rates. We found a direct association between OM chemical diversity and δ 13 C-CH 4 in peat porewater; where higher substrate diversity was positively correlated with enriched δ 13 C-CH 4 in fen sites. Oxidized sulfur-containing compounds, produced by Sphagnum, were further hypothesized to control GHG emissions by acting as electron acceptors for a sulfate-reducing electron transport chain, inhibiting methanogenesis in peat bogs. These results suggest that warming-induced permafrost thaw might increase organic matter lability, in subset of sites that become wetlands, and shift biogeochemical processes toward faster decomposition with an increasing proportion of carbon released as CH 4 .

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Glacial‐Interglacial and Millennial‐Scale Changes in Nitrous Oxide Emissions Pathways and Source Regions

During the transition from the Last Glacial Maximum (LGM) to the Holocene, the atmospheric N 2 O mole fraction increased by 80 nmol mol −1 . Using ice core measurements of N 2 O isotopomer ratios, we show that this increase was driven by increases in both nitrification and denitrification, with the relative partitioning between both production pathways depending on the assumed isotopic end-member source signatures. Similarly, we also attribute a 35 nmol mol −1 N 2 O mole fraction increase during the Heinrich Stadial 4/Dansgaard Oeschger 8 (HS4/DO8) millennial-scale event to increases in both N 2 O production pathways. In contrast, the 25 nmol mol −1 N 2 O mole fraction decrease during the Younger Dryas was driven almost exclusively by a decrease in nitrification. The deglacial and HS4/DO8 increases in N 2 O production occurred in both marine and terrestrial environments, with the terrestrial source responding faster to warming by about two centuries. Constraints on changes in nitrification and denitrification emissions are robust and consistent with previous studies showing the sensitivity of N 2 O emissions to abrupt Northern Hemisphere warming. This study demonstrates for the first time the importance of both denitrification and nitrification pathways in driving source changes. Absolute emissions are more uncertain due to uncertainty about source isotopomer signatures. For instance, the contribution of denitrification to emissions at the LGM shifts from (65 ± 10) % to (91 ± 6) % when factoring in isotope enrichment due to partial reduction of N 2 O to N 2 during denitrification. Reducing uncertainty in source signatures will increase the power of ice core N 2 O isotope records in deducing environmental change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Energy Emissions Accounting Methods Can Determine Whether Direct Air Capture with Storage Achieves Net Removal

The voluntary carbon market within the United States has expanded rapidly in recent years and enabled private companies and other organizations to provide revenue streams to carbon dioxide removal (CDR) technologies. For a CDR technology to participate in the voluntary carbon market (VCM), the emissions associated with constructing and operating the technology must be less than the CO 2 captured from the atmosphere. Assessing the extent to which this is true for direct air capture with storage (DACS), a relatively energy-intensive CDR technology, strongly depends on the accounting method used to assess the emissions intensity of purchased energy. We simulate the hourly weather-dependent operation of sorbent- and solvent-based DACS in California, Louisiana, Texas, and Wyoming, representing a wide range of local weather and electric and natural gas grid compositions. In all cases, the single most important emissions accounting decision is the method used to estimate the emissions intensity of purchased grid electricity, which varies the calculated net removal by −1049% to +108%. All other factors influencing net removal introduce a variation of at most ±14%. No electricity emissions accounting method is universally conservative across all scenarios, and none is objectively more accurate. High-spatiotemporal-resolution, high-quality, publicly available data sets and models for electricity emissions accounting do not currently exist and are urgently needed to enable standardization of emissions accounting methods to more accurately determine the true emissions impacts of DACS and other energy-intensive facilities.

54 ENVIRONMENTAL SCIENCES↗