Search NASA⌕ Search

SEARCH · Search NASA

Results for “WRF”

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 181 records · Page 10

Evaluating Implementations of the Immersed Boundary Method in the Weather Research and Forecasting Model

The terrain-following coordinate system used by many atmospheric models can cause numerical instabilities due to discretization errors as resolved terrain slopes increase and the grid becomes highly skewed. The immersed boundary (IB) method, which does not require the grid to conform to the terrain, has been shown to alleviate these errors, and has been used successfully for high-resolution atmospheric simulations over steep terrain, including vertical building surfaces. Since many previous applications of IB methods to atmospheric models have used very fine grid resolution (5 m or less), the present study seeks to evaluate IB method performance over a range of grid resolutions and aspect ratios. Two classes of IB algorithms, velocity reconstruction and shear stress reconstruction, are tested within the common framework of the Weather Research and Forecasting (WRF) Model. Performance is evaluated in two test cases, one with flat terrain and the other with the topography of Askervein Hill, both under neutrally stratified conditions. WRF-IB results are compared to similarity theory, observations, and native WRF results. Despite sensitivity to the location at which the IB intersects the model grid, the velocity reconstruction IB method shows consistent performance when used with a hybrid RANS/LES surface scheme. The shear stress reconstruction IB method is not sensitive to the grid intersection, but is less consistent and near-surface velocity errors can occur at coarse resolutions. This study represents an initial investigation of IB method variability across grid resolutions in WRF. Future work will focus on improving IB method performance at intermediate to coarse resolutions.

54 ENVIRONMENTAL SCIENCES↗

High-Resolution Regional Atmosphere–Ocean–Wave Coupled Simulations of Hurricane Henri (2021)

To explore the integrated effects of ocean and ocean surface wave related physical processes on TC simulations, a set of three model simulations is performed. * In experiment 'A', the event is modeled using the Weather Research Forecasting (WRF) model alone with prescribed Sea Surface Temperature (SST) at 6-hour intervals. * In experiment ‘AO,’ WRF is coupled with the Finite Volume Community Ocean Model (FVCOM), enabling variable exchange between atmosphere and ocean, but without considering ocean surface wave-related physical processes. * In experiment ‘AOW’, WRF, FVCOM, and Simulating WAves Nearshore (SWAN) exchange variables with each other every hour through the OASIS3-MCT Coupler to allow direct and indirect atmosphere-ocean-wave interactions. * Observational data are also included in this dataset (Dropsonde, HRD-Radar, NDBC_wave). All simulations are initialized at 18:00 UTC on August 19, 2021, within a domain encompassing the western North Atlantic Ocean. The atmospheric domain features a horizontal resolution of 3 km. The ocean domain, which covers a substantial portion of the WRF ocean domain, employs an unstructured triangular grid with resolutions ranging from 3 km near the coast to 9 km in the open ocean, effectively resolving the complex coastline of the U.S. Northeast Coast. Initial and boundary conditions for the atmosphere model are obtained from the 6-hourly 0.25° NCEP (National Centers for Environmental Prediction) Global Forecast System (GFS; NCEP, 2015) data.

17 WIND ENERGY↗

High-Resolution Regional Atmosphere–Ocean–Wave Coupled Simulations of Hurricane Henri (2021)

To explore the integrated effects of ocean and ocean surface wave related physical processes on TC simulations, a set of three model simulations is performed. * In experiment 'A', the event is modeled using the Weather Research Forecasting (WRF) model alone with prescribed Sea Surface Temperature (SST) at 6-hour intervals. * In experiment ‘AO,’ WRF is coupled with the Finite Volume Community Ocean Model (FVCOM), enabling variable exchange between atmosphere and ocean, but without considering ocean surface wave-related physical processes. * In experiment ‘AOW’, WRF, FVCOM, and Simulating WAves Nearshore (SWAN) exchange variables with each other every hour through the OASIS3-MCT Coupler to allow direct and indirect atmosphere-ocean-wave interactions. * Observational data are also included in this dataset (Dropsonde, HRD-Radar, NDBC_wave). All simulations are initialized at 18:00 UTC on August 19, 2021, within a domain encompassing the western North Atlantic Ocean. The atmospheric domain features a horizontal resolution of 3 km. The ocean domain, which covers a substantial portion of the WRF ocean domain, employs an unstructured triangular grid with resolutions ranging from 3 km near the coast to 9 km in the open ocean, effectively resolving the complex coastline of the U.S. Northeast Coast. Initial and boundary conditions for the atmosphere model are obtained from the 6-hourly 0.25° NCEP (National Centers for Environmental Prediction) Global Forecast System (GFS; NCEP, 2015) data.

17 WIND ENERGY↗

High-Resolution Regional Atmosphere–Ocean–Wave Coupled Simulations of Hurricane Henri (2021)

To explore the integrated effects of ocean and ocean surface wave related physical processes on tropical cyclone simulations, a set of three model simulations is performed. * In experiment 'A', the event is modeled using the Weather Research Forecasting (WRF) model alone with prescribed Sea Surface Temperature (SST) at 6-hour intervals. * In experiment ‘AO,’ WRF is coupled with the Finite Volume Community Ocean Model (FVCOM), enabling variable exchange between atmosphere and ocean, but without considering ocean surface wave-related physical processes. * In experiment ‘AOW’, WRF, FVCOM, and Simulating WAves Nearshore (SWAN) exchange variables with each other every hour through the OASIS3-MCT Coupler to allow direct and indirect atmosphere-ocean-wave interactions. * Observational data are also included in this dataset (Dropsonde, HRD-Radar, NDBC_wave). All simulations are initialized at 18:00 UTC on August 19, 2021, within a domain encompassing the western North Atlantic Ocean. The atmospheric domain features a horizontal resolution of 3 km. The ocean domain, which covers a substantial portion of the WRF ocean domain, employs an unstructured triangular grid with resolutions ranging from 3 km near the coast to 9 km in the open ocean, effectively resolving the complex coastline of the U.S. Northeast Coast. Initial and boundary conditions for the atmosphere model are obtained from the 6-hourly 0.25° NCEP (National Centers for Environmental Prediction) Global Forecast System (GFS; NCEP, 2015) data. These CSV files are derived from the NetCDF files in the c0 dataset. Unlike the original format, where geographic coordinates were stored in a separate file, each CSV now embeds the corresponding latitude and longitude values alongside the measured or simulated variables.

17 WIND ENERGY↗

Evaluation and development of surface layer scheme representation of temperature inversions over boreal forests in Arctic wintertime conditions

Abstract. In this study, the Noah land surface model used in conjunction with the Mellor–Yamada–Janjić surface layer scheme (hereafter, Noah-MYJ) and the Noah multiphysics scheme (Noah-MP) from the Weather Research and Forecasting (WRF) 4.5.1 mesoscale model are evaluated with regard to their performance in reproducing positive temperature gradients over forested areas in the Arctic winter. First, simplified versions of the WRF schemes, recoded in Python, are compared with conceptual models of the surface layer in order to gain insight into the dependence of the temperature gradient on the wind speed at the top of the surface layer. It is shown that the WRF schemes place strong limits on the turbulent collapse, leading to lower surface temperature gradient at low wind speeds than in the conceptual models. We implemented modifications to the WRF schemes to correct this effect. The original and modified versions of Noah-MYJ and Noah-MP are then evaluated compared to long-term measurements at the Ameriflux Poker Flat Research Range, a forest site in interior Alaska. Noah-MP is found to perform better than Noah-MYJ because the former is a two-layer model which explicitly takes into account the effect of the forest canopy. Indeed, a non-negligible temperature gradient is maintained below the canopy at high wind speeds, leading to overall larger gradients than in the absence of vegetation. Furthermore, the modified versions are found to perform better than the original versions of each scheme because they better reproduce strong temperature gradients at low wind speeds.

Geology↗

IM3/HyperFACETS Thermodynamic Global Warming (TGW) Simulation Datasets

Publication For a thorough description of the methods, see the peer-reviewed paper: Jones, A.D., Rastogi, D., Vahmani, P. et al. Continental United States climate projections based on thermodynamic modification of historical weather. Sci Data 10, 664 (2023). https://doi.org/10.1038/s41597-023-02485-5 Overview The IM3 / HyperFACETS climate simulations provide 40-year historical (1980-2019) as well as four 80-year future simulations (2020-2099) over the U.S. The future simulations are split into near (2020-2059) and far future (2060-2099) segments. The future scenarios span a range of plausible changes in future climate (both Global Circulation Model (GCM) and Representative Concentration Pathways/Shared Socioeconomic Pathway (RCP/SSP) dimensions). The simulations provide climate variables with high spatiotemporal resolution (25 hourly variables and 207 3-hourly variables at 12 km2). The datasets are generated using dynamical downscaling with the WRF (Weather Research and Forecasting) model (version 4.2.1) and therefore preserve physical consistency across variables. WRF is a state-of-the-art, fully compressible, non-hydrostatic, mesoscale numerical weather prediction model. WRF is coupled with an urban canopy model (UCM), which resolves urban surfaces. The future scenarios were developed using a thermodynamic global warming approach where past events are replayed under a range of future warming conditions. These scenarios therefore provide a perspective on potential increases in extreme event intensity, geographic scope, and duration, with previously non-extreme conditions potentially crossing new thresholds to be considered extreme by today's standards. This approach is not intended to estimate future changes in extreme event frequency that might result from changes in large-scale atmospheric dynamics. This dataset has NOT been bias corrected. A bias corrected version of selected variables is under development and will be released here when available. Scenarios Files Data for each scenario is provided in weekly NetCDF files. 25 variables are available at hourly resolution, and 207 variables are available at three-hourly resolution. Spatial resolution is 12km and spans the conterminous United States (CONUS), including some areas of Canada and Mexico, resulting in a grid of 424 by 299 cells. The spatial projection is a Lambert Conformal Conic with the following proj-string: "+proj=lcc +lat_0=40.0000076293945 +lon_0=-97 +lat_1=30 +lat_2=45 +x_0=0 +y_0=0 +R=6370000 +units=m +no_defs". The available scenarios and simulation periods are listed below: historical | 1980 - 2019 rcp45cooler | 2020 - 2059 rcp45cooler | 2060 - 2099 rcp45hotter | 2020 - 2059 rcp45hotter | 2060 - 2099 rcp85cooler | 2020 - 2059 rcp85cooler | 2060 - 2099 rcp85hotter | 2020 - 2059 rcp85hotter | 2060 - 2099 * The first year (1979, 2019, and 2059) of data within each scenario represents a model warmup period and should not be used. These are located in the `spinup_files` directory. Historical year 2020 is considered an extra year of data beyond the simulation period and can be found in the `additional_files` directory. For information on specific variables and a more in-depth discussion of methodology, please refer to the data landing page at https://tgw-data.msdlive.org. Delta Warming Files The global and CONUS warming deltas for each scenario are provided in degrees Celsius annually and monthly. Restart Files Yearly restart files are provided for each scenario which can be used to restart the WRF model at a particular point in time. Spinup Files The first year of data within each simulation period represents a model warmup period and should not be used. The files are provided here for the sake of reproducibility. Additional Files Additional years of data are provided as an extension of the historic simulation.

Jones, Andrew D.↗

Urban Parameters Arizona Urban Corridor 100m

132 Urban parameters based on building physical dimensions and location were generated for the cities in six Arizona Counties at 100m resolution using the NATURF model. To use the binary file with WRF, the binary file and the index file must be placed in their own directory in WRF_GEOG and accessed in the same way NUDAPT44 would be accessed.

Dumas, Melissa [ORNL] (ORCID:0000000233190846)↗

Urban Parameters Washington, DC 100m

132 Urban parameters based on building physical dimensions and location were generated for Washington, District of Columbia at 100m resolution using the NATURF model. To use the binary file with WRF, the binary file and the index file must be placed in their own directory in WRF_GEOG and accessed in the same way NUDAPT44 would be accessed.

Dumas, Melissa [ORNL] (ORCID:0000000233190846)↗

Developing a 20-year high-resolution wind data set for Puerto Rico

The purpose of this study is to develop a high-resolution wind resource data set for Puerto Rico as part of the wind resource assessment in the Puerto Rico Grid Resilience and Transition to 100 % Renewable Energy Study (PR100) project. The Weather and Research and Forecasting (WRF) is used to model 20-years (2001-2020) of wind resource data on a 3-km grid for the Puerto Rico region. Because accurately representing the planetary boundary layer (PBL)-physics in the WRF is key to accurately model low-level wind speed, 11 different PBL schemes are examined and evaluated to find a WRF configuration that can provide the most accurate wind data. The modeled wind speeds are validated with buoy observations, and the Shin-Hong PBL scheme is selected to generate the final wind data set. Here, the overall results demonstrate that the high-resolution data correctly represents the climate of Puerto Rico (e.g., a dominant easterly trade wind). At 160 m above ground, the southern regions of Puerto Rico show strong offshore wind speed (>8 m/s on average) throughout the daytime and nighttime. For land-based wind, average wind speed ranges from 5 m/s-7 m/s during daytime, whereas notable high wind speeds (>9 m/s) appear in the mountain regions during nighttime.

17 WIND ENERGY↗

A Machine‐Learning‐Assisted Stochastic Cloud Population Model as a Parameterization of Cumulus Convection

Abstract A machine‐learning‐assisted stochastic cloud population model is coupled with the Advanced Research Weather Research and Forecasting (WRF) model to represent fluctuations in the cloud‐base mass flux associated with the life cycles and interactions among cumulus convection cells. In this cloud population model, the size distribution and the associated cloud‐base mass flux of the convective cells are related to their previous state and to the change in the total convective area via a transition function. The convective area tendency in turn is assumed to depend on the cloud‐base mass flux that is resolved by the host WRF model. The transition function is represented by a single hidden‐layer neural network trained by the evolution of convective cell size distributions in a 1‐km grid‐spacing WRF simulation run over the Australian Monsoon region. At every grid point of the host model, the cloud population model predicts the cell size and cloud‐base mass flux distributions from which a random sample of cells is fed to an entraining parcel model that calculates precipitation as well as the associated liquid water potential temperature and total moisture tendencies. These tendencies are averaged over the cells and provided to the host model. Several regional simulations are performed over tropical and midlatitude domains to test this as a potential approach to scale‐aware parameterization. It is shown that such an approach could be a new promising path to simulating realistic precipitation statistics and propagation of precipitation associated with the Madden‐Julian Oscillation while maintaining realistic depictions of the diurnal cycle over both land and ocean.

54 ENVIRONMENTAL SCIENCES↗

Data and scripts associated with “Sequential Precipitation Input Tagging (SPIT) to Estimate Water Transit Times and Hydrologic Tracer Dynamics within Water-Tagging Enabled Hydrologic Models” (v3)

This data package is associated with the publication “Sequential Precipitation Input Tagging (SPIT) to Estimate Water Transit Times and Hydrologic Tracer Dynamics within Water-Tagging Enabled Hydrologic Models” submitted to Journal of Advances in Modeling Earth Systems (Butler et al. 2025). This study developed the Sequential Precipitation Input Tagging (SPIT) framework to tag input precipitation and estimate water transit times and hydrologic tracers. SPIT tags all precipitation events at regular intervals over an extended period (monthly tags over seven years) in a hydrologic model from 2016-2022. SPIT is applied at six National Ecological Observatory Network (NEON) sites across the continental United States to calculate transit time distributions (TTD) and derive from these mean transit times (MTT), fractions of young water (Fyw), and hydrologic tracer concentrations in stream water (δ18O) within a water-tagging enabled version of the Weather Research and Forecast (WT-WRF-Hydro) model with national water model (NWM) configurations. We go on to validate WT-WRF-Hydro estimates against Butler et al. (2023), who analyzed the same NEON sites using stable water isotope data to estimate water transit times. This new tracking method provides a detailed picture of water movement and helps improve predictions about water availability in the future. This data package was originally published in January 2025. It was updated May 2025 (v2; new and modified files) and October 2025 (v3; new and modified files). File and folder names were not revised to indicate changes. See the change history section in the readme for more details. This data package contains the data and scripts used to develop the SPIT framework WT-WRF-Hydro (Water Tagging Weather Research and Forecasting Hydrologic) model and is associated with the following GitHub repository: https://github.com/zbutler33/SPIT-Framework. This data package contains five parent folders: (1) “Manipulated_outputs”, (2) “Metadata”, (3) “Observed”, (4) “Outputs”, and (5) “Scripts”. Each of these parent folders contains additional subfolders and files. Please see the FLMD (“v*_Butler_2024_WT_WRF_Hydro_flmd.csv”) for a list of all the files contained in this data package and descriptions for each. See the data dictionary (“v*_Butler_2024_WT_WRF_Hydro_dd.csv”) for definitions and units of all of the tabular (files ending in “.csv” and ".tsv") column headers.

54 ENVIRONMENTAL SCIENCES↗

wrfhydro

This dataset contains a continuous WRF-Hydro simulation for the CACTI field campaign. The simulated period is 1-Aug-2018 through 21-Mar-2019 with hourly model output. The simulation is made with two grids that match the two outer LASSO-CACTI grids, with 7.5 km and 2.5 km grid spacings. Forcing data for WRF-Hydro was obtained from the ERA5-Land product. Values from this WRF-Hydro simulation were used to initialize the soil state for each LASSO-CACTI simulation.

54 ENVIRONMENTAL SCIENCES↗

Quantifying sensitivity in numerical weather prediction-modeled offshore wind speeds through an ensemble modeling approach

A decade of research has shown that numerical weather prediction (NWP)-modeled wind speeds can be highly sensitive to the inputs and setups within the NWP model. For wind resource characterization applications, this sensitivity is often addressed by constructing a range of setups and selecting the one that best validates against observations. However, this approach is not possible in areas that lack high-quality hub height observations, especially offshore wind areas. In such cases, techniques to quantify and disseminate confidence in NWP-modeled wind speeds in the absence of observations are needed. We address this need in the present study and propose best practices for quantifying the spread in NWP-modeled wind speeds. We implement an ensemble approach in which we consider 24 different setups to the Weather Research and Forecasting (WRF) model. We construct the ensemble by considering variations in WRF version, WRF namelist, atmospheric forcing, and sea surface temperature (SST) forcing. Our analysis finds that the standard deviation produces more consistent estimates compared to the interquartile range and tends to be the more conservative estimator for ensemble variability. We further find that model spread increases closer to the surface and on shorter time scales. In conclusion, we explore methods to attribute total ensemble variability to the different ensemble components (e.g., atmospheric forcing and SST product) and find that contributions by components also vary depending on time scale. We anticipate that the methods and results presented in this paper will provide a reasonable foundation for further research into ensemble-based wind resource data sets.

17 WIND ENERGY↗

Biochemical and structural characterization of enzymes in the 4-hydroxybenzoate catabolic pathway of lignin-degrading white-rot fungi

White-rot fungi (WRF) are the most efficient lignin-degrading organisms in nature. However, their capacity to use lignin-related aromatic compounds, such as 4-hydroxybenzoate, as carbon sources has only been described recently. Previously, the hydroxyquinol pathway was proposed for the bioconversion of these compounds in fungi, but gene- and structure-function relationships of the full enzymatic pathway remain uncharacterized in any single fungal species. Here, we characterize seven enzymes from two WRF, Trametes versicolor and Gelatoporia subvermispora, which constitute a four-enzyme cascade from 4-hydroxybenzoate to β-ketoadipate via the hydroxyquinol pathway. Furthermore, we solve the crystal structure of four of these enzymes and identify mechanistic differences with the closest bacterial and fungal structural homologs. Overall, this research expands our understanding of aromatic catabolism by WRF and establishes an alternative strategy for the conversion of lignin-related compounds to the valuable molecule β-ketoadipate, contributing to the development of biological processes for lignin valorization.

59 BASIC BIOLOGICAL SCIENCES↗

Comparative transcriptomics of CAZy enzymes in white- and brown-rot agaricomycetes: Evolutionary insights into lignocellulose degradation and the relevance of GH16 glycoside hydrolase functional divergence

White-rot and brown-rot fungi (WRF and BRF, respectively) decompose lignocellulose, the main structural component of plant biomass, through distinct mechanisms. This study examines the transcriptomic responses of three WRF species (Pleurotus ostreatus, Phanerochaete chrysosporium, and Heterobasidion irregulare) and two BRF species (Fomitopsis schrenkii and Rhodonia placenta) grown on poplar wood (W) and glucose (G) as sole carbon sources. RNA-seq analysis revealed upregulation of carbohydrate-active enzymes (CAZymes) linked to lignocellulose degradation. WRF displayed a broader enzymatic repertoire, whereas BRF employed a more selective approach. Among these responses, GH16 glycoside hydrolases were consistently upregulated across all species, including BRF. Since GH16 enzymes are involved in both plant hemicellulose modification and cell wall remodeling, their wood-induced expression may reflect multiple processes rather than a single conserved wood-decay mechanism. Structural and phylogenetic analyses revealed species-specific divergence, consistent with potential functional specialization. These findings broaden our understanding of fungal enzymatic strategies and highlight GH16 enzymes as candidates for further study. Beyond the immediate context of wood decay, this work has broader implications for fungal ecology, evolutionary biology, and biotechnological applications such as biomass conversion and sustainable bioenergy.

Fungal ecological strategies↗

Observation and Simulation of Methane Plumes During the Morning Boundary Layer Transition

Abstract Methane (CH 4 ) contributes significantly to global warming. However, accurate identification of CH 4 sources for reducing CH 4 emissions is often hampered by inadequate accuracy and spatiotemporal coverage of CH 4 detection, and lack of accurate CH 4 forward modeling used in top‐down inversion systems. In this study, a field experiment was conducted in Pampa, Texas using two CH 4 sensors (LI‐COR and OGI camera) to detect CH 4 releases. We investigated whether high‐resolution simulations using the Weather Research and Forecasting (WRF) model with greenhouse gases (WRF‐GHG) could accurately simulate the CH 4 plumes in the presence of evolving atmospheric boundary layer from sunrise to noon. CH 4 plumes showed substantial variation in time. At a release rate of ∼17.5 kg hr −1 , the maximum enhancement of CH 4 measured by LI‐COR was 2.6 ppm at sunrise (7:36 a.m.), 250 m from the release location. Within half an hour after sunrise, this enhancement decreased to 0.3–0.4 ppm. The enhancement was 0.2 ppm by 10:00 a.m. and further dropped to less than 0.1 ppm after 11:30 a.m. Due to the low temperature at sunrise, the OGI camera failed to detect the CH 4 plume. The WRF‐GHG large‐eddy simulation (LES) with 32 m grid spacing successfully reproduced these CH 4 enhancements. In situ measurements together with numerical simulations illustrate the impact of the transition from a stable boundary layer in the early morning to a convective boundary layer at noon on the dispersion of CH 4 plumes. Additionally, CH 4 plumes from a cattle farm in Oklahoma are briefly examined using the same modeling approach.

Hu, Xiao‐Ming [Center for Analysis and Prediction ↗

Evaluating Sea Breezes and Associated Convective Cloud Evolution in the Model Gray Zone

We characterize convective clouds associated with sea‐breeze circulations (SBC) using multi‐agency observations and multi‐case ensemble model simulations. The focus is on assessing convective cloud lifecycle properties and their merging behavior, as well as the environmental conditions they are embedded in, particularly SBC features. In total, 46 SBC days over the Houston‐Galveston region are selected and simulated using the Weather Research and Forecasting (WRF) model at a gray zone scale with a forecast‐like parameterization setup. Advanced techniques, including change‐point detection, a Lagrangian cloud tracking method, and a newly developed cell merging and splitting detection algorithm, are applied and/or developed for this study. Our findings indicate that the WRF model at 1 km grid spacing well represents the thermodynamic conditions over the region, as well as SBC timing and intensity. However, for the associated convective cells, WRF overestimates the 30‐dBZ echo top height, cell area, and maximum radar reflectivity compared to radar observations. This overestimation is potentially due to under‐resolved entrainment processes, an overestimated merging frequency, and the overestimation of updraft intensity. Furthermore, the model exhibits a deficiency in simulating congestus clouds, showing a more rapid transition from shallow to deep convection compared to observed behavior. Moreover, observations indicate stronger, deeper, and wider clouds when merging happens. Conversely, in simulations, the merging process does not necessarily lead to higher or longer‐lived cells, as many cases experience rapid and frequent merging and splitting which may result in more variance in convective updraft velocity during the convection lifetime.

54 ENVIRONMENTAL SCIENCES↗

Physics informed deep neural network embedded in a chemical transport model for the Amazon rainforest

Secondary organic aerosols (SOA) are fine particles in the atmosphere, which interact with clouds, radiation and affect the Earth’s energy budget. SOA formation involves chemistry in gas phase, aqueous aerosols, and clouds. Simulating these chemical processes involve solving a stiff set of differential equations, which are computationally expensive steps for three-dimensional chemical transport models. Deep neural networks (DNNs) are universal function approximators that could be used to represent the complex nonlinear changes in aerosol physical and chemical processes; however, key challenges such as generalizability to extended time periods, preservation of mass balance, simulating sparse model outputs, and maintaining physical constraints have limited their use in atmospheric chemistry. Here, we develop an approach of using a physics-informed DNN that overcomes previous such challenges and demonstrates its applicability for the chemical formation processes of isoprene epoxydiol SOA (IEPOX-SOA) over the Amazon rainforest. The DNN is trained with data generated by simulating IEPOX-SOA over the entire atmospheric column, using the Weather Research and Forecasting Model coupled with Chemistry (WRF-Chem). The trained DNN is then embedded within WRF-Chem to replace the computationally expensive default solver of IEPOX-SOA formation. The trained DNN predictions generalizes well with the default model simulation of the IEPOX-SOA mass concentrations and its size distribution (20 size bins) over several days of simulations in both dry and wet seasons. The embedded DNN reduces the computational expense of WRF-Chem by a factor of 2. Our approach shows promise in terms of application to other computationally expensive chemistry solvers in climate models.

54 ENVIRONMENTAL SCIENCES↗