Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation modelling”

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,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Glacial geology of the Hudson Mountains, Amundsen Sea sector, West Antarctica

The Hudson Mountains are situated in the eastern Amundsen Sea sector of the West Antarctic Ice Sheet, adjacent to Pine Island Glacier. They form a volcanic field of 17 stratovolcanoes and parasitic vents, preserved as nunataks. Two former tributaries of Pine Island Glacier (Larter and Lucchitta glaciers) flow through the mountains. Here we present a detailed study of the glacial geology of the area. We describe field observations and measurements of geomorphological features from 15 of the nunataks, meltwater ponds found on the surface of three nunataks and supraglacial features (ice dolines) from two sites near the present grounding line. Together these provide constraints on the past ice sheet extent, flow pathways and thermal regime, and enhance our understanding of the present hydrological regime – all of which are important as context for the observed modern ice sheet behaviour. We find evidence suggesting that all nunataks in the Hudson Mountains were covered by ice during the Last Glacial Maximum (defined here as 26.5-19 ka) and have since deglaciated. Faceted and polished erratic cobbles and boulders of exotic lithologies (syenites, alkali granites, granites, granodiorites, tonalites and gabbros) are numerous and perched on nunatak surfaces. A marked difference between the dominant erratic lithologies on nunataks adjacent to Pine Island Glacier (granite) and Lucchitta Glacier (granodiorite-tonalite) indicates that the ice sheet was transporting clasts from at least two distinct upstream source regions. The similarity in degree of weathering suggests, however, that all the erratics were transported by one phase of (warm-based) glaciation; their presence on or close to the summits of all except one nunatak indicates that the ice sheet during that time was at least 700 m thicker than present. These results are consistent with ice sheet model simulations which suggest that all nunataks in the Hudson Mountains were completely submerged by the Last Glacial Maximum ice sheet.

58 GEOSCIENCES↗
DOE OSTI2021

Role of emission controls in reducing the 2050 climate change penalty for PM2.5 in China

Previous studies demonstrated that global warming can lead to deteriorated air quality even when anthropogenic emissions were kept constant, which has been called a climate change penalty on air quality. It is expected that anthropogenic emissions will decrease significantly in the future considering the aggressive emission control actions in China. However, the dependence of climate change penalty on the choice of emission scenario is still uncertain. To fill this gap, we conducted multiple independent model simulations to investigate the response of PM2.5 to future (2050) climate warming (RCP8.5) in China but with different emission scenarios, including the constant 2015 emissions, the 2050 CLE emissions (based on Current Legislation), and the 2050 MTFR emissions (based on Maximum Technically Feasible Reduction). For each set of emissions, we estimate climate change penalty as the difference in PM2.5 between a pair of simulations with either 2015 or 2050 meteorology. Under 2015 emissions, we find a PM2.5 climate change penalty of 1.43 µg m-3 in Eastern China, leading to an additional 35,000 PM2.5-related premature deaths [95% confidence interval (CI), 21,000-40,000] by 2050. However, the PM2.5 climate change penalty weakens to 0.24 µg m-3 with strict anthropogenic emission controls under the 2050 MTFR emissions, which decreases the associated PM2.5-related deaths to 17,000. The smaller MTFR climate change penalty contributes 14% of the total PM2.5 decrease when both emissions and meteorology are changed from 2015 to 2050, and 24% of total health benefits associated with this PM2.5 decrease in Eastern China. This finding suggests that controlling anthropogenic emissions can effectively reduce the climate change penalty on PM2.5 and its associated premature deaths, even though a climate change penalty still occurs even under MTFR. Strengthened controls on anthropogenic emissions is key to attaining air quality targets and protecting human health in the context of future global climate change.

Liu, Song↗
DOE OSTI2022

Using probabilistic solar power forecasts to inform flexible ramp product procurement for the California ISO

How can independent system operators (ISOs) take advantage of probabilistic solar forecasts to lower generation costs and improve reliability of power systems? We discuss one three-step approach for doing so, focusing on how such forecasts might help the California Independent System Operator (CAISO) prepare unexpected net load ramps, where net load equals gross demand minus wind and solar production. First, we enhance an existing solar forecasting system to provide well-calibrated hours-ahead probabilistic forecasts. We then relate the degree of uncertainty reflected in the forecasted prediction intervals (independent variables) to error distributions for net load ramp forecasts for the CAISO real-time market (dependent variable) using machine learning and quantile regression. Projected ramp forecast errors conditioned on solar uncertainty are translated into flexible ramp requirements that therefore reflect real-time meteorological and solar conditions, improving on typical ISO procedures. Detailed descriptions are provided on the quantile regression and kth-nearest neighbor categorization methods for accomplishing that translation. Finally, a multiple time-scale look-ahead market simulation model is applied to a 118-bus IEEE Reliability Test System, modified to represent the CAISO generation mix and demand distributions. The model runs quantify how solar-conditioned ramp requirements can, first, decrease operating costs by reducing requirements compared to often conservative unconditional methods and, second, decrease generation scarcity events and consequently improve reliability by increasing flexibility requirements at times when unconditional forecast-based requirements understate actual ramp uncertainty. Solar-conditioned ramp requirements are found to reduce generation operating costs by about 2% for the test system (which would be equivalent to over $\$100$ million per year for a CAISO-size system).

14 SOLAR ENERGY↗
DOE OSTI2023

Enabling scientific machine learning in MOOSE using Libtorch

A neural-network-based machine learning interface has been developed for the Multiphysics Object-Oriented Simulation Environment (MOOSE). The interface relies on Libtorch, the C++ front-end of PyTorch, and enables an online interaction between modern machine learning algorithms and all the existing simulation, modeling, and analysis processes available in MOOSE. New capabilities in MOOSE include the native generation and training of artificial neural networks together with options to load pretrained neural networks in TorchScript format. Furthermore, the MOOSE stochastic tools module (MOOSE-STM) has been enhanced with neural network-based surrogate and reduced-order model generation options for efficient stochastic analyses. Lastly, a reinforcement learning capability has been added to MOOSE-STM for the interactive control and optimization of complex multiphysics problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

An efficient method to identify uncertainties of WRF-Solar variables in forecasting solar irradiance using a tangent linear sensitivity analysis

Uncertainty in predicting solar energy resources introduces major challenges in power system management and necessitates the development of reliable probabilistic solar forecasts. As the first part of the development of probabilistic forecasts based on the Weather Research and Forecasting model with solar extensions (WRF-Solar), this study presents a tangent linear approach to identify input variables responsible for the largest uncertainties in predicting surface solar irradiance and clouds. A tangent linear analysis is capable of efficiently investigating sensitivities of output variables with respect to various input variables of WRF-Solar because this approach avoids the computational burden of perturbing the initial conditions of individual input variables. We develop tangent linear models (TLMs) for six WRF-Solar physics packages that control the formation and dissipation of clouds and solar radiation, and we evaluate the validity of TLMs using a linearity test. The tangent linear sensitivity analysis is conducted under various scenarios based on satellite observations and model simulations to consider realistic input conditions. A simple method is used to quantify the impact of the uncertainty of input variables on the output variables from the TLMs. The results demonstrate that uncertainties in the output variables that are the focus of this study—including global horizontal irradiance, direct normal irradiance, cloud mixing ratio, cloud tendency, cloud fraction, and sensible and latent heat fluxes—are highly sensitive to uncertainties in 14 input variables. This study indicates that the tangent linear method can identify key variables of physics modules in WRF-Solar that can be stochastically perturbed to generate ensemble-based probabilistic forecasts.

14 SOLAR ENERGY↗
DOE OSTI2022

Additive manufacturing and testing of a ceramic heat exchanger for high-temperature and high-pressure applications for concentrating solar power

Heat exchangers with excellent corrosion and oxidation resistances are essential for the next-generation concentrating solar power (CSP) plants using a molten salt heat transfer fluid and a supercritical CO 2 Brayton power cycle that can operate at temperatures > 700 °C for higher efficiencies. Techniques were developed for additively manufacturing ceramic materials for applications involving high temperature, high pressure, and high corrosion resistance as needed for the CSP application. Based on a previous ceramic heat exchanger design with its cross-section geometry of the flow channels optimized for heat transfer and mechanical strength, lab-scale prototype heat exchangers, with integrated headers and incorporated flow channel dimensional compensations, were fabricated by using the binder jetting process. The printed prototype heat exchangers were successfully densified through the processes of liquid polymer infiltration and pyrolysis. In addition, thermophysical properties of the densified silicon carbide parts were measured to provide necessary design information. Experimental heat transfer testing of the lab-scale prototype heat exchanger was conducted, and the experimental data agreed reasonably well with the simulated results. Furthermore, this agreement provides validation of the simulation models for their applicability to the development of a full-scale ceramic heat exchanger.

14 SOLAR ENERGY↗
DOE OSTI2026

Comparing Delay-, Distance-, and Cordon-Based Congestion Pricing Strategies Via Large-Scale Simulation

This study compares the impacts of delay-, distance-, and cordon-based congestion pricing strategies for Austin, Texas, using the POLARIS agent-based activity-based travel demand simulation model. This approach enables agent-level heterogeneity and realistic choice options (including destination, mode, and activity scheduling) for dynamic traffic assignment and congestion feedbacks across a major metro region, which are features lacking in past work. To ensure comparability, distance-based tolls were set to generate the same revenue as delay-based tolling of $3.5 M/day, averaging $1.17/resident/day or $0.42/vehicle-trip. Delay-based pricing delivers 44% lower network delay and 13% lower VHT compared to the no-toll baseline, levels unmatched by other pricing strategies. At the height of the AM peak, drivers pay up to $0.13/mile on average, though most links in the network remain untolled. Distance-based pricing is the most effective at reducing VMT (by 4%), but VHT reductions (of 6%) primarily stem from drivers selecting closer destinations, achieving only one-fourth the delay reduction of delay-based pricing. Across various implementations of delay- and distance-based pricing, the results suggest that spatial variations of tolls are far more important than temporal variations. Cordon tolls produce minimal impacts at the network-wide level, but offer substantial delay reductions inside the cordon. Other major findings include: 1) delay-based pricing increases trip-making during the PM peak period due to backward shifts in discretionary-activity start times by higher-income residents; and 2) tolls’ spatial impacts, including changes in network flows and tolls paid by residents, vary substantially between delay- and distance-based pricing strategies.

Agent-based modeling↗
DOE OSTI2023

Freeze range of a condensing gas propagating in a liquid helium-cooled tube

Understanding air propagation and condensation following a catastrophic vacuum break in particle accelerator beamlines cooled by liquid helium is crucial for maintaining the operational safety of these facilities. Previous experimental investigations on nitrogen gas propagation in both normal liquid helium (He I) and superfluid helium (He II) cooled copper tubes unveiled a nearly exponential deceleration of the gas propagation. A comprehensive theoretical model incorporating gas dynamics, heat transfer, and condensation mechanisms has been developed, which effectively reproduces various key experimental observations. An intriguing phenomenon uncovered in our model simulation is that the gas propagation appears to nearly stop beyond a certain distance from the location where condensation starts. We refer to this distance as the freeze range. In this paper, we present our systematic study of the freeze range at various inlet mass fluxes and tube diameters. We show that the results can be well described by a simple correlation. The underlying physical mechanism that supports this useful correlation is explained. Knowing the freeze range may allow accelerator engineers to develop protocols for controlling frost-layer contamination in the beamline tubes, which is of great practical importance.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Oil in Alaska North Slope Gas Hydrate Reservoir: Micro-CT and Flow Simulation Insights into Permeability

Gas hydrate-bearing sands on the Alaska North Slope (ANS) host minor volumes of crude oil whose impact on formation permeability has never been quantified relative to a gas hydrate reservoir system. Here, we combine in situ pressure-core microcomputed-tomography (μ-CT), thermogravimetric analysis (TGA), gas-chromatography–mass-spectrometry (GC–MS), scanning-electron microscopy coupled with energy-dispersive spectroscopy (SEM-EDS), and Stokes-flow simulation to (i) quantify oil saturation, (ii) infer its source, and (iii) evaluate its influence on permeability after gas-hydrate dissociation. μ-CT and SEM-EDS imaging identify the host sediment as silty, with some grain-coating clay present. Eight representative 3003-voxel subvolumes extracted from a preserved pressure core (Core 15P-3c, from the B1 sand (unit B), HYDRATE 02 Geo Data Well (GDW), 886.96–887.29 m measured depth, MD) exhibit porosities of 0.34–0.37 and an average oil saturation of 0.09 ± 0.03, in comparison to an independent TGA-based estimate of 0.16. GC-MS fingerprints obtained from an adjacent depressurized core (Core 17P-1, 891.37–891.44 m MD), together with oil pore habit revealed by μ-CT, suggest that the oil is partly native but also includes contributions from mineral oil-based drilling fluid contamination. Further analysis of oil–porewater interactions and flushing tests confirms that native oil saturation in the analyzed samples ranges from 0.04 to 0.08 with an average of approximately 0.06, and the oil is highly immobile. Flow simulations demonstrate that the native oil saturation of ≈0.06 reduces permeability by approximately 50% due to both pore blockage and increased flow-path tortuosity. These findings reveal the potential presence of native oil within the ANS gas hydrate reservoir and suggest that future simulation models may need to account for its impact on permeability to improve long-term performance predictions of gas and water production.

02 PETROLEUM↗
DOE OSTI2025

Techno-Economic Analysis and Life Cycle Assessment for the Separation of 2,3-Butanediol from Fermentation Broth Using Liquid–Liquid Extraction

It is energy-intensive to separate dilute 2,3-butanediol (2,3-BDO) (<10 wt %) from the aqueous phase of fermentation broth for sustainable aviation fuel (SAF) using conventional distillation. Liquid–liquid extraction (LLE) using oleyl alcohol as a solvent in a membrane extractor to extract BDO from water can significantly reduce the energy cost and minimize the potential emulsion. In an Aspen Plus model simulation, 95.2% BDO recovery and 97.1% BDO purity have been achieved using this LLE method with solvent recovery and heat integration. Here, the thermal energy cost was estimated to be 4.57 MJ/kg BDO, which is only about 16.8% of the lower heating value (LHV) of the BDO. This method consumes 81% less energy than the cascade distillation and reduces about $0.46/GGE (gasoline gallon equivalent) to the minimal fuel selling price. Meanwhile, the greenhouse gas (GHG) emission is 62% lower than petroleum-based jet fuel production and 34% less than using the cascade distillation.

09 BIOMASS FUELS↗
DOE OSTI2025

Combining Reactive Quantum-Mechanical Molecular-Dynamics Simulations with Mutagenesis, Crystallography, and Enzyme Kinetics to Reveal Plausible Steps of Isocyanide Hydratase Catalysis

A complete understanding of enzyme mechanisms requires atomistic details of chemical reactions. Quantum-based molecular dynamics simulations (QMD) are a potential source of this information, but trade-offs between accuracy and computational cost have limited their use. We previously developed extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) methods that leverage a negligible compromise in accuracy to substantially decrease the cost of QMD simulations. Here, we develop a reactive QMD approach using the latest XL-BOMD formulation, which enables efficient simulations of highly reactive systems, and use it to investigate mechanisms of intermediate formation in isocyanide hydratase (ICH) catalysis. In QMD simulations, molecular analogs of ICH active site residues reacted with para-nitrophenyl isocyanide, forming a thioimidate. Analysis of simulated atomic configurational and charge dynamics revealed a pathway where protonation of the isocyanide carbon occurs prior to thioimidate formation and suggested a possible role of Asp17 as a proton donor in the early phase of ICH catalysis. To test whether the pathway seen using the reactive QMD approach might be relevant to ICH catalysis, we performed X-ray crystallography and pre-steady-state enzyme kinetics studies of wild-type and D17N mutant ICH. Both the structure and kinetics are sensitive to the D17N mutation in a manner that is consistent with the order of the reaction steps seen in the simulations. Mobile protons play essential roles in many enzymes, yet they are difficult to observe experimentally, making the ordering of proton-dependent steps ambiguous in many enzyme mechanisms. The ability to directly simulate model reactions for the design of experiments that provide information about enzyme mechanisms involving mobile protons demonstrates the significance of our reactive QMD approach and motivates further biological applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Inorganic Sulfur Species Formed upon Heterogeneous OH Oxidation of Organosulfates: A Case Study of Methyl Sulfate

In this report recent studies reveal that organosulfates at the particle surface can be oxidized by gas-phase OH radicals with significant rates. Inorganic sulfur species, such as the bisulfate ion (HSO 4 - ) and sulfate ion (SO 4 2- ), can be formed upon these heterogeneous oxidation processes through the formation and subsequent reactions of sulfate radical anions (SO 4 •- ) in the particle phase. However, the amount of inorganic sulfur species produced in these heterogeneous oxidation reactions is not known. We investigate the heterogeneous OH oxidation of sodium methyl sulfate (CH 3 SO 4 Na), the smallest organosulfate detected in atmospheric particles, using an oxidation flow reactor at a relative humidity of 75%. We quantify the kinetics by measuring the decay of CH 3 SO 4 Na and the amount of HSO 4 - and SO 4 2- formed upon oxidation using ion chromatography. Kinetic measurements determine the heterogeneous OH reaction rate to be (5.72 ± 0.14) x 10 -13 cm 3 molecule -1 s -1 , with an effective OH uptake coefficient, γ eff , of 0.31 ± 0.06. The molar yield of inorganic sulfur species, defined as the total number of moles of HSO 4 - and SO 4 2- formed per mole of CH 3 SO 4 Na consumed upon oxidation, is found to be significant and has an average value of 0.62 ± 0.18 upon oxidation. A kinetic model is developed to describe the kinetics and inorganic sulfur species formation upon oxidation. Model simulations suggest that CH 3 SO 4 Na tends to decompose rapidly into formaldehyde and SO 4 •- , and the reaction of SO 4 •- with CH 3 SO 4 Na plays a significant role in both governing the kinetics and the formation of inorganic sulfur species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Tight Coupling of Surface and In-Plant Biochemistry and Convection Governs Key Fine Particulate Components over the Amazon Rainforest

We report that throughout the troposphere over the Amazon rainforest, isoprene epoxydiol-secondary organic aerosols (IEPOX-SOA) are key components of submicrometer particle mass. IEPOX-SOA is thought to form by multiphase chemical pathways in the atmosphere, but using detailed model simulations, we show that these pathways are strongly inhibited by the solid particle phase state prevailing in the upper troposphere. We find that 2-methyltetrol gases formed by in-plant biochemical oxidation pathways, their partitioning to fine particles, and dissolution in cloud droplets can explain over 90% of IEPOX-SOA mass concentrations in the upper troposphere. In agreement with previous PM10 chiral measurements, our simulations indicate that plant biochemical processes are the dominant source of particulate 2-methyltetrols in the lower troposphere as well (>70%). Our results imply that in-plant biochemistry plays a central but previously unidentified role in aerosol processes and atmosphere-biosphere-climate interactions over the rainforest.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Analysis of Brackish Water Desalination for Municipal Uses: Case Studies on Challenges and Opportunities

Brackish water is a widely available, nontraditional water resource that can augment limited freshwater supplies. Although brackish water desalination has been continuously implemented in the United States and worldwide, it is necessary to reduce further its energy consumption, costs, and environmental impacts. This study conducted technoeconomic analyses to evaluate the current desalination and brine management technologies, focusing on the key factors and opportunities for sustainable brackish water desalination for municipal uses. Three case studies were selected as baseline representative of different geographic and operational conditions, including water quality, plant size, pretreatment, desalination, and concentrate management. Here, the technoeconomic analyses and model simulations identified challenges, opportunities, and research priorities to achieve specific pathways for enhanced brackish water desalination regarding levelized costs of water, electricity intensity, water recovery, zero liquid discharge, and brine valorization.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Thermal Stability of π-Conjugated n -Ethylene-Glycol-Terminated Quaterthiophene Oligomers: A Computational and Experimental Study

This work represents a joint computational and experimental study on a series of $\textit{n}$-ethylene glycol (PEO$\textit{n}$)-terminated quaterthiophene (4T) oligomers for 1 < $\textit{n}$ < 10 to elucidate their self-assembly behavior into a smectic-like lamellar phase. This study builds on an earlier study for $\textit{n}$ = 4 that showed that our model predictions were consistent with experimental data on the melting behavior and structure of the lamellar phase, with the latter consisting of crystal-like 4T domains and liquid-like PEO4 domains. The present study aims to understand how the length of the terminal PEO$\textit{n}$ chains modulates the disordering temperature of the lamellar phase and hence the relative stability of the ordered structure. In this work, a simplified bilayer model, where the 4T domains are not explicitly described, is put forward to efficiently estimate the disordering effect of the PEO domains with increasing n; this method is first validated by correctly predicting that layers of alkyl (PE)-capped 4T oligomers (for 1 < $\textit{n}$ < 10) stay ordered at room temperature. Both 4T-domain implicit and explicit model simulations reveal that the order-disorder temperature decreases with the length of the PEO capping chains, as the associated increase in conformational entropy drives a tendency toward disorder that overtakes the cohesive energy, keeping the ordered packing of the 4T domains.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Revisiting Recent U.S. Heat Waves in a Warmer and More Humid Climate

The frequency and intensity of heat waves in the United States is projected to increase in the 21st century. We investigate dry and humid heat waves in a pair of high-resolution model simulations that constrain large-scale atmospheric circulation, to isolate the thermodynamic impacts on characteristics of present and future heat waves over the United States. The two kinds of heat waves show differences in mean intensity, amplitude, duration, and frequency over the Southeast, Northeast, and Midwest, while their characteristics are largely similar in the drier central and western United States. In a warmer climate, relative humidity is projected to decrease during dry heat waves, whereas it remains unchanged during humid heat waves. However, the overall increase in daily maximum temperature intensifies the heat stress during future humid and dry heat waves across all regions. With large-scale circulation constrained, these simulations emphasize the importance of thermodynamic drivers in determining future heat wave characteristics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Modeling the Sources and Transport Processes During Extreme Ammonia Episodes in the U.S. Corn Belt

Abstract Atmospheric ammonia (NH 3 ) is the primary form of reactive nitrogen (N r ) and a precursor of ammonium (NH 4 + ) aerosols. Ammonia has been linked to adverse impacts on human health, the loss of ecosystem biodiversity, and plays a key role in aerosol radiative forcing. The midwestern United States is the major NH 3 source in North America because of dense livestock operations and the high use of synthetic nitrogen fertilizers. Here, we combine tall‐tower (100 m) observations in Minnesota and Weather Research and Forecasting model coupled with Chemistry (WRF‐Chem) modeling to investigate high and low NH 3 emission episodes within the U.S. Corn Belt to improve our understanding of the distribution of emission sources and transport processes. We examined observations and performed model simulations for cases in February through November of 2017 and 2018. The results showed the following: (1) Peak emissions in November 2017 were enhanced by above‐normal air temperatures, implying a Q 10 (i.e., the change in NH 3 emissions for a temperature increase of 10°C) of 2.5 for emissions. (2) The intensive livestock emissions rom northern Iowa, approximately 400 km away from the tall tower, accounted for 17.6% of the abundance in tall‐tower NH 3 mixing ratios. (3) Ammonia mixing ratios in the innermost domain 3 frequently (i.e., 336 hr, 48% of November 2017) exceeded 5.3 ppb, an important air quality health standard. (4) In November 2017, simulated NH 3 net ecosystem exchange (the difference between NH 3 emissions and dry deposition) accounted for 60–65% of gross NH 3 emissions for agricultural areas and was 2.8–3.1 times the emissions of forested areas. (5) We estimated a mean annual NH 3 net ecosystem exchange of 1.60 ± 0.06 nmol · m −2 · s −1 for agricultural lands and −0.07 ± 0.02 nmol · m −2 · s −1 for forested lands. These results imply that future warmer fall temperatures will enhance agricultural NH 3 emissions, increase the frequency of dangerous NH 3 episodes, and enhance dry NH 3 deposition in adjacent forested lands.

Hu, Cheng↗