Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “WRF simulation”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2020

Comparison of Surface Fluxes Derived from CYGNSS and Simulated by WRF Model: An MJO Case Study

This study focuses on ocean surface fluxes, mainly the latent heat flux, and their impact on MJO propagation and associated precipitation structures over the Indian Ocean and Maritime Continent. The Coupled-Ocean-Atmosphere-Wave-Sediment Transport (COAWST) model is used to simulate two MJO events during the 2017-2018 season: the December 10 - January 20, 2017 case, which maintained its strong precipitation signal over the Maritime Continent, and the March 1 - 20, 2018 case, which was weaker and did not propagate through the Maritime Continent. Both simulated MJO events show positive biases in surface rainfall compared with GPM IMERG data. During the MJO suppressed phase, the simulations rain more often than the observations. During the active phase, the westward propagating precipitation structures are more organized and much stronger compared with the observations, sometimes forming westward propagating cyclones that weakened the eastward precipitation signals. Two aspects of the surface flux interactions are investigated: the impact of the domain mean surface fluxes, and the impact of storm scale circulations and their interactions with local surface fluxes. Both aspects affect water vapor budget, atmosphere instability and mean flow, through which convection initiation, organization, and propagation are influenced. Model sensitivity tests with different radiation, microphysics, PBL schemes and nudging schemes indicate that in the control simulations, higher SST and surface fluxes, especially during the suppressed period, are the main reason of rainfall overestimation compared with IMERG data. The strong westward propagating signals are caused by both increased atmosphere instability and reduced mean wind shear. Unfortunately, the small differences in mean SST and surface fluxes between different model sensitivity tests are all within the satellite observation error margin, and cannot be directly corroborated by observations. One of the advantages of CYGNSS satellites is that they observe ocean surface wind and heat fluxes underneath strong rainfall events such as the convective systems associated with MJO active phases. Currently we are comparing CYGNSS level 2 surface fluxes retrievals and the model simulations in order to better understand the second aspect of the MJO and surface fluxes interactions, and how this affects MJO strengths and propagations. The interactive atmosphere-ocean-wave model also provides cases that directly comparing satellite observables (the bistatic radar cross section) and the model simulations (through CYGNSS satellite simulator). These discrepancies are more prominent in coupled ocean simulations, mainly due to higher SST and enhanced surface fluxes.

Li, Xiaowen↗
NASA NTRS2012

Two Distinct Modes in One-Day Rainfall Event during MC3E Field Campaign: Analyses of Disdrometer Observations and WRF-SBM Simulation

A unique microphysical structure of rainfall is observed by the surface laser optical Particle Size and Velocity (Parsivel) disdrometers on 25 April 2011 during Midlatitude Continental Convective Clouds Experiment (MC3E). According to the systematic differences in rainfall rate and bulk effective droplet radius, the sampling data can be divided into two groups; the rainfall mostly from the deep convective clouds has relatively high rainfall rate and large bulk effective droplet radius, whereas the reverse is true for the rainfall from the shallow wrm clouds. The Weather Research and Forecasting model coupled with spectral bin microphysics (WRF-SBM) successfully reproduces the two distinct modes in the observed rainfall microphysical structure. The results show that the up-to-date model can demonstrate how the cloud physics and the weather condition on the day are involved in forming the unique rainfall characteristic.

disdrometers↗
NASA NTRS2011

Passive Microwave Remote Sensing of Falling Snow and Associated GPM Field Campaigns

Retrievals of falling snow from space represent one of the next important challenges for the atmospheric, hydrological, and energy budget scientific communities. Historically, retrievals of falling snow have been difficult due to the relative insensitivity of satellite rain-based channels as used in the past. We emphasize the use of high frequency passive microwave channels (85-200 GHz) since these are more sensitive to the ice in clouds and have been used to estimate falling snow from space. While satellite-based remote sensing provides global coverage of falling snow events and the science is relatively new, retrievals are still undergoing development with challenges remaining. There are several current satellite sensors, though not specifically designed for estimating falling snow, are capable of measuring snow from space. These include NOAA's AMSU-B, the MHS sensors, and CloudSat radar. They use high frequency (greater than 85 GHz) passive and active microwave and millimeter-wave channels that are sensitive to the scattering from ice and snow particles in the atmosphere. Sensors with water vapor channels near 183 GHz center line provide opaqueness to the Earth's surface features that can contaminate the falling snow signatures, especially over snow covered surface. In addition, the Global Precipitation Measurement (GPM) mission scheduled for launch in 2013 is specifically designed to measure both liquid rain and frozen snow precipitation. Since falling snow from space is the next precipitation measurement challenge from space, information must be determined in order to guide retrieval algorithm development for these current and future missions. This information includes thresholds of detection for various sensor channel configurations, snow event system characteristics, and surface types. For example, can a lake effect snow system with low cloud tops having an ice water content (IWC) at the surface of 1.0 gram per cubic meter be detected? If this information is known, we can focus retrieval efforts on detectable storms and concentrate advances on achievable results. In this work, the focus is to determine thresholds of detection for falling snow for various snow conditions over land and lake surfaces. The results rely on simulated Weather Research Forecasting (WRF) simulations of falling snow cases [9] since simulations provide all the information to determine the measurements from space and the ground truth. This analysis relies on data from the Canadian CloudSat/CALIPSO Validation Program (C3VP) field campaign held from October 31, 2006 through March 1, 2007. In January 2012 GPM will return to the C3VP area for the GPM Cold Precipitation Experiment (GCPEx). This presentation will describe the thresholds-of-detection procedure and results, as well as the field campaign details.

Skofronick-Jackson, Gail↗
NASA NTRS2012

Detection Thresholds of Falling Snow from Satellite-Borne Active and Passive Sensors

Precipitation, including rain and snow, is a critical part of the Earth's energy and hydrology cycles. Precipitation impacts latent heating profiles locally while global circulation patterns distribute precipitation and energy from the equator to the poles. For the hydrological cycle, falling snow is a primary contributor in northern latitudes during the winter seasons. Falling snow is the source of snow pack accumulations that provide fresh water resources for many communities in the world. Furthermore, falling snow impacts society by causing transportation disruptions during severe snow events. In order to collect information on the complete global precipitation cycle, both liquid and frozen precipitation must be collected. The challenges of estimating falling snow from space still exist though progress is being made. These challenges include weak falling snow signatures with respect to background (surface, water vapor) signatures for passive sensors over land surfaces, unknowns about the spherical and non-spherical shapes of the snowflakes, their particle size distributions (PSDs) and how the assumptions about the unknowns impact observed brightness temperatures or radar reflectivities, differences in near surface snowfall and total column snow amounts, and limited ground truth to validate against. While these challenges remain, knowledge of their impact on expected retrieval results is an important key for understanding falling snow retrieval estimations. Since falling snow from space is the next precipitation measurement challenge from space, information must be determined in order to guide retrieval algorithm development for these current and future missions. This information includes thresholds of detection for various sensor channel configurations, snow event system characteristics, snowflake particle assumptions, and surface types. For example, can a lake effect snow system with low (approx 2.5 km) cloud tops having an ice water content (IWC) at the surface of 0.25 g / cubic m and dendrite snowflakes be detected? If this information is known, we can focus retrieval efforts on detectable storms and concentrate advances on achievable results. Here, the focus is to determine thresholds of detection for falling snow for various snow conditions over land and lake surfaces. The results rely on simulated Weather Research Forecasting (WRF) simulations of falling snow cases since simulations provide all the information to determine the measurements from space and the ground truth. Sensitivity analyses were performed to better ascertain the relationships between multifrequency microwave and millimeter-wave sensor observations and the falling snow/underlying field of view. In addition, thresholds of detection for various sensor channel configurations, snow event system characteristics, snowflake particle assumptions, and surface types were studied. Results will be presented for active radar at Ku, Ka, and W-band and for passive radiometer channels from 10 to 183 GHz.

Gail Skofronick-Jackson↗
NASA NTRS2012

Detection Thresholds of Falling Snow from Satellite-Borne Active and Passive Sensors

Precipitation, including rain and snow, is a critical part of the Earth's energy and hydrology cycles. In order to collect information on the complete global precipitation cycle and to understand the energy budget in terms of precipitation, uniform global estimates of both liquid and frozen precipitation must be collected. Active observations of falling snow are somewhat easier to estimate since the radar will detect the precipitation particles and one only needs to know surface temperature to determine if it is liquid rain or snow. The challenges of estimating falling snow from passive spaceborne observations still exist though progress is being made. While these challenges are still being addressed, knowledge of their impact on expected retrieval results is an important key for understanding falling snow retrieval estimations. Important information to assess falling snow retrievals includes knowing thresholds of detection for active and passive sensors, various sensor channel configurations, snow event system characteristics, snowflake particle assumptions, and surface types. For example, can a lake effect snow system with low (~2.5 km) cloud tops having an ice water content (Iwe) at the surface of 0.25 g m-3 and dendrite snowflakes be detected? If this information is known, we can focus retrieval efforts on detectable storms and concentrate advances on achievable results. Here, the focus is to determine thresholds of detection for falling snow for various snow conditions over land and lake surfaces. The analysis relies on simulated Weather Research Forecasting (WRF) simulations of falling snow cases since simulations provide all the information to determine the measurements from space and the ground truth. Results are presented for active radar at Ku, Ka, and W-band and for passive radiometer channels from 10 to 183 GHz (Skofronick-Jackson, et al. submitted to IEEE TGRS, April 2012). The notable results show: (1) the W-Band radar has detection thresholds more than an order of magnitude lower than the future GPM sensors, (2) the cloud structure macrophysics influences the thresholds of detection for passive channels, (3) the snowflake microphysics plays a large role in the detection threshold for active and passive instruments, (4) with reasonable assumptions, the passive 166 GHz channel has detection threshold values comparable to the GPM DPR Ku and Ka band radars with ~0.05 g m-3 detected at the surface, or an ~0.5-1 mm hr-l melted snow rate (equivalent to 0.5-2 cm hr-l solid fluffy snowflake rate).

Jackson, Gail↗
NASA NTRS2014

The Sensitivity of WRF Daily Summertime Simulations over West Africa to Alternative Parameterizations. Part 1: African Wave Circulation

The performance of the NCAR Weather Research and Forecasting Model (WRF) as a West African regional-atmospheric model is evaluated. The study tests the sensitivity of WRF-simulated vorticity maxima associated with African easterly waves to 64 combinations of alternative parameterizations in a series of simulations in September. In all, 104 simulations of 12-day duration during 11 consecutive years are examined. The 64 combinations combine WRF parameterizations of cumulus convection, radiation transfer, surface hydrology, and PBL physics. Simulated daily and mean circulation results are validated against NASA's Modern-Era Retrospective Analysis for Research and Applications (MERRA) and NCEP/Department of Energy Global Reanalysis 2. Precipitation is considered in a second part of this two-part paper. A wide range of 700-hPa vorticity validation scores demonstrates the influence of alternative parameterizations. The best WRF performers achieve correlations against reanalysis of 0.40-0.60 and realistic amplitudes of spatiotemporal variability for the 2006 focus year while a parallel-benchmark simulation by the NASA Regional Model-3 (RM3) achieves higher correlations, but less realistic spatiotemporal variability. The largest favorable impact on WRF-vorticity validation is achieved by selecting the Grell-Devenyi cumulus convection scheme, resulting in higher correlations against reanalysis than simulations using the Kain-Fritch convection. Other parameterizations have less-obvious impact, although WRF configurations incorporating one surface model and PBL scheme consistently performed poorly. A comparison of reanalysis circulation against two NASA radiosonde stations confirms that both reanalyses represent observations well enough to validate the WRF results. Validation statistics for optimized WRF configurations simulating the parallel period during 10 additional years are less favorable than for 2006.

variability↗
NASA NTRS2008

Application of the CloudSat and NEXRAD Radars Toward Improvements in High Resolution Operational Forecasts

As computational power increases, operational forecast models are performing simulations with higher spatial resolution allowing for the transition from sub-grid scale cloud parameterizations to an explicit forecast of cloud characteristics and precipitation through the use of single- or multi-moment bulk water microphysics schemes. investments in space-borne and terrestrial remote sensing have developed the NASA CloudSat Cloud Profiling Radar and the NOAA National Weather Service NEXRAD system, each providing observations related to the bulk properties of clouds and precipitation through measurements of reflectivity. CloudSat and NEXRAD system radars observed light to moderate snowfall in association with a cold-season, midlatitude cyclone traversing the Central United States in February 2007. These systems are responsible for widespread cloud cover and various types of precipitation, are of economic consequence, and pose a challenge to operational forecasters. This event is simulated with the Weather Research and Forecast (WRF) Model, utilizing the NASA Goddard Cumulus Ensemble microphysics scheme. Comparisons are made between WRF-simulated and observed reflectivity available from the CloudSat and NEXRAD systems. The application of CloudSat reflectivity is made possible through the QuickBeam radiative transfer model, with cautious application applied in light of single scattering characteristics and spherical target assumptions. Significant differences are noted within modeled and observed cloud profiles, based upon simulated reflectivity, and modifications to the single-moment scheme are tested through a supplemental WRF forecast that incorporates a temperature dependent snow crystal size distribution.

Molthan, A. L.↗
NASA NTRS2015

A Module for Assimilating Hyperspectral Infrared Retrieved Profiles into the Gridpoint Statistical Interpolation System for Unique Forecasting Applications

Hyperspectral infrared sounder radiance data are assimilated into operational modeling systems however the process is computationally expensive and only approximately 1% of available data are assimilated due to data thinning as well as the fact that radiances are restricted to cloud-free fields of view. In contrast, the number of hyperspectral infrared profiles assimilated is much higher since the retrieved profiles can be assimilated in some partly cloudy scenes due to profile coupling other data, such as microwave or neural networks, as first guesses to the retrieval process. As the operational data assimilation community attempts to assimilate cloud-affected radiances, it is possible that the use of retrieved profiles might offer an alternative methodology that is less complex and more computationally efficient to solve this problem. The NASA Short-term Prediction Research and Transition (SPoRT) Center has assimilated hyperspectral infrared retrieved profiles into Weather Research and Forecasting Model (WRF) simulations using the Gridpoint Statistical Interpolation (GSI) System. Early research at SPoRT demonstrated improved initial conditions when assimilating Atmospheric Infrared Sounder (AIRS) thermodynamic profiles into WRF (using WRF-Var and assigning more appropriate error weighting to the profiles) to improve regional analysis and heavy precipitation forecasts. Successful early work has led to more recent research utilizing WRF and GSI for applications including the assimilation of AIRS profiles to improve WRF forecasts of atmospheric rivers and assimilation of AIRS, Cross-track Infrared and Microwave Sounding Suite (CrIMSS), and Infrared Atmospheric Sounding Interferometer (IASI) profiles to improve model representation of tropopause folds and associated non-convective wind events. Although more hyperspectral infrared retrieved profiles can be assimilated into model forecasts, one disadvantage is the retrieved profiles have traditionally been assigned the same error values as the rawinsonde observations when assimilated with GSI. Typically, satellitederived profile errors are larger and more difficult to quantify than traditional rawinsonde observations (especially in the boundary layer), so it is important to appropriately assign observation errors within GSI to eliminate potential spurious innovations and analysis increments that can sometimes arise when using retrieved profiles. The goal of this study is to describe modifications to the GSI source code to more appropriately assimilate hyperspectral infrared retrieved profiles and outline preliminary results that show the differences between a model simulation that assimilated the profiles as rawinsonde observations and one that assimilated the profiles in a module with the appropriate error values.

Berndt, Emily↗
NASA NTRS2013

Atmospheric Turbulence Estimates from a Pulsed Lidar

Estimates of the eddy dissipation rate (EDR) were obtained from measurements made by a coherent pulsed lidar and compared with estimates from mesoscale model simulations and measurements from an in situ sonic anemometer at the Denver International Airport and with EDR estimates from the last observation time of the trailing vortex pair. The estimates of EDR from the lidar were obtained using two different methodologies. The two methodologies show consistent estimates of the vertical profiles. Comparison of EDR derived from the Weather Research and Forecast (WRF) mesoscale model with the in situ lidar estimates show good agreement during the daytime convective boundary layer, but the WRF simulations tend to overestimate EDR during the nighttime. The EDR estimates from a sonic anemometer located at 7.3 meters above ground level are approximately one order of magnitude greater than both the WRF and lidar estimates - which are from greater heights - during the daytime convective boundary layer and substantially greater during the nighttime stable boundary layer. The consistency of the EDR estimates from different methods suggests a reasonable ability to predict the temporal evolution of a spatially averaged vertical profile of EDR in an airport terminal area using a mesoscale model during the daytime convective boundary layer. In the stable nighttime boundary layer, there may be added value to EDR estimates provided by in situ lidar measurements.

Pruis, Matthew J.↗
NASA NTRS2021

Cold Season Performance of the NU-WRF Regional Climate Model in the Great Lakes Region

As Earth’s largest collection of freshwater, the Laurentian Great Lakes have enormous ecological and socio-economic value. Their basin has become a regional hotspot of climatic and limnological change, potentially threatening its vital natural resources. Consequentially, there is a need to assess the current state of climate models regarding their performance across the Great Lakes region and develop the next generation of high-resolution regional climate models to address complex limnological processes and lake-atmosphere interactions. In response to this need, the current paper focuses on the generation and analysis of a 20-member ensemble of 3-km National Aeronautics and Space Administration (NASA)-Unified Weather Research and Forecasting (NU-WRF) simulations for the 2014-2015 cold season. The study aims to identify the model’s strengths and weaknesses; optimal configuration for the region; and the impacts of different physics parameterizations, coupling to a 1D lake model, time-variant lake-surface temperatures, and spectral nudging. Several key biases are identified in the cold-season simulations for the Great Lakes region, including an atmospheric cold bias that is amplified by coupling to a 1D lake model but diminished by applying the Community Atmosphere Model radiation scheme and Morrison microphysics scheme; an excess precipitation bias; anomalously early initiation of fall lake turnover and subsequent cold lake bias; excessive and overly persistent lake ice cover; and insufficient evaporation over Lakes Superior and Huron. The research team is currently addressing these key limitations by coupling NU-WRF to a 3D lake model in support of the next generation of regional climate models for the critical Great Lakes Basin.

Michael Notaro↗
NASA NTRS2008

Use of Vertically Integrated Ice in WRF-Based Forecasts of Lightning Threat

Previously reported methods of forecasting lightning threat using fields of graupel flux from WRF simulations are extended to include the simulated field of vertically integrated ice within storms. Although the ice integral shows less temporal variability than graupel flux, it provides more areal coverage, and can thus be used to create a lightning forecast that better matches the areal coverage of the lightning threat found in observations of flash extent density. A blended lightning forecast threat can be constructed that retains much of the desirable temporal sensitivity of the graupel flux method, while also incorporating the coverage benefits of the ice integral method. The graupel flux and ice integral fields contributing to the blended forecast are calibrated against observed lightning flash origin density data, based on Lightning Mapping Array observations from a series of case studies chosen to cover a wide range of flash rate conditions. Linear curve fits that pass through the origin are found to be statistically robust for the calibration procedures.

McCaul, E. W., jr.↗
NASA NTRS2012

Translation of Land Surface Model Accuracy and Uncertainty into Coupled Land-Atmosphere Prediction

Land-atmosphere (L-A) Interactions playa critical role in determining the diurnal evolution of both planetary boundary layer (PBL) and land surface heat and moisture budgets, as well as controlling feedbacks with clouds and precipitation that lead to the persistence of dry and wet regimes. Recent efforts to quantify the strength of L-A coupling in prediction models have produced diagnostics that integrate across both the land and PBL components of the system. In this study, we examine the impact of improved specification of land surface states, anomalies, and fluxes on coupled WRF forecasts during the summers of extreme dry (2006) and wet (2007) land surface conditions in the U.S. Southern Great Plains. The improved land initialization and surface flux parameterizations are obtained through the use of a new optimization and uncertainty estimation module in NASA's Land Information System (US-OPT/UE), whereby parameter sets are calibrated in the Noah land surface model and classified according to a land cover and soil type mapping of the observation sites to the full model domain. The impact of calibrated parameters on the a) spinup of the land surface used as initial conditions, and b) heat and moisture states and fluxes of the coupled WRF Simulations are then assessed in terms of ambient weather and land-atmosphere coupling along with measures of uncertainty propagation into the forecasts. In addition, the sensitivity of this approach to the period of calibration (dry, wet, average) is investigated. Finally, tradeoffs of computational tractability and scientific validity, and the potential for combining this approach with satellite remote sensing data are also discussed.

Santanello, Joseph A.↗
NASA NTRS2016

Impact of Soil Moisture Assimilation on Land Surface Model Spin-Up and Coupled LandAtmosphere Prediction

Advances in satellite monitoring of the terrestrial water cycle have led to a concerted effort to assimilate soil moisture observations from various platforms into offline land surface models (LSMs). One principal but still open question is that of the ability of land data assimilation (LDA) to improve LSM initial conditions for coupled short-term weather prediction. In this study, the impact of assimilating Advanced Microwave Scanning Radiometer for EOS (AMSR-E) soil moisture retrievals on coupled WRF Model forecasts is examined during the summers of dry (2006) and wet (2007) surface conditions in the southern Great Plains. LDA is carried out using NASAs Land Information System (LIS) and the Noah LSM through an ensemble Kalman filter (EnKF) approach. The impacts of LDA on the 1) soil moisture and soil temperature initial conditions for WRF, 2) land-atmosphere coupling characteristics, and 3) ambient weather of the coupled LIS-WRF simulations are then assessed. Results show that impacts of soil moisture LDA during the spin-up can significantly modify LSM states and fluxes, depending on regime and season. Results also indicate that the use of seasonal cumulative distribution functions (CDFs) is more advantageous compared to the traditional annual CDF bias correction strategies. LDA performs consistently regardless of atmospheric forcing applied, with greater improvements seen when using coarser, global forcing products. Downstream impacts on coupled simulations vary according to the strength of the LDA impact at the initialization, where significant modifications to the soil moisture flux- PBL-ambient weather process chain are observed. Overall, this study demonstrates potential for future, higher-resolution soil moisture assimilation applications in weather and climate research.

soil moisture↗
NASA NTRS2020

Evaluation of Stratospheric Intrusions and Biomass Burning Plumes on the Vertical Distribution of Tropospheric Ozone Over the Midwestern U.S.

Naturally occurring ozone rich Stratosphere-to-Troposphere Transport (STT) intrusions and biomass burning plumes reaching the surface can contribute to exceedances of the U.S. National Ambient Air Quality Standards (NAAQS) for ground-level ozone (70 ppbv implemented in 2015). In addition, fires can inject significant pollution into the free troposphere where it can be transported long distances. The combined air quality impacts from these sources on ozone has only been analyzed in a few case studies for the Midwest U.S. Here we study ozone impacts in a major Midwestern city, for the first time in St. Louis, Missouri, using a series of ozonesonde profiles taken during the SEAC4RS (Studies of Emissions and Atmospheric Composition, Clouds and Climate Coupling by Regional Surveys) field campaign in August-September 2013. All ozonesondes showed some enhancements above the background profile levels (~55 ppbv) throughout each tropospheric column. Two models were used to estimate and quantify ozone origins within the columns. A chemical transport model identified enhancements from STT equivalent to 10to15 ppbv over the background with a 10 to15% contribution overall to the column. Two FLEXPART-WRF simulations, one with smoke in the boundary layer and the other with smoke above, identified enhancements from biomass burning equivalent to 10to 80 ppbv. Overall, the total biomass burning contribution is 15to 30% of the total column. Five ozonesondes showed signatures of mixed biomass burning plumes and STT intrusions. During this study period, biomass burning in the western U.S. contributed 70% to ozone enhancements in the total column compared to 3% from the central U.S and 27% from other areas.

National Ambient Air Quality Standards for ground‐↗
NASA NTRS2020

Assessment of Natural and Anthropogenic Aerosol Air Pollution In the Middle East Using MERRA-2, CAMS Data Assimilation Products, and High-Resolution WRF-Chem Model Simulations

Modern-Era Retrospective analysis for Research and Applications v.2 (MERRA-2), Copernicus Atmosphere Monitoring Service Operational Analysis (CAMS-OA), and a high-resolution regional Weather Research and Forecasting model coupled with chemistry (WRF-Chem) were used to evaluate natural and anthropogenic particulate matter (PM) air pollution in the Middle East (ME) during 2015–2016. Two Moderate Resolution Imaging Spectrometer (MODIS) retrievals – combined product Deep Blue and Deep Target (MODIS-DB&DT) and Multi-Angle Implementation of Atmospheric Correction (MAIAC) – and Aerosol Robotic Network (AERONET) aerosol optical depth (AOD) observations as well as in situ PM measurements for 2016 were used for validation of the WRF-Chem output and both assimilation products. MERRA-2 and CAMS-OA assimilate AOD observations. WRF-Chem is a free-running model, but dust emission in WRF-Chem is tuned to fit AOD and aerosol volume size distributions obtained from AERONET. MERRA-2 was used to construct WRF-Chem initial and boundary conditions both for meteorology and chemical and aerosol species. SO2 emissions in WRF-Chem are based on the novel OMI-HTAP SO2 emission dataset. The correlation with the AERONET AOD is highest for MERRA-2 (0.72–0.91), MAIAC (0.63–0.96), and CAMSOA (0.65–0.87), followed by MODIS-DB&DT (0.56–0.84) and WRF-Chem (0.43–0.85). However, CAMS-OA has a relatively high positive mean bias with respect to AERONET AOD. The spatial distributions of seasonally averaged AODs from WRF-Chem, assimilation products, and MAIAC are well correlated with MODIS-DB&DT AOD product. MAIAC has the highest correlation (R = 0.8), followed by MERRA-2 (R = 0.66), CAMS-OA (R = 0.65), and WRF-Chem (R = 0.61). WRF-Chem, MERRA-2, and MAIAC underestimate and CAMS-OA overestimates MODIS-DB&DT AOD. The simulated and observed PM concentrations might differ by a factor of 2 because it is more challenging for the model and the assimilation products to reproduce PM concentration measured within the city. Although aerosol fields in WRF-Chem and assimilation products are entirely consistent, WRF-Chem is preferable for analysis of regional air quality over the ME due to its higher spatial resolution and better SO2 emissions. The WRF-Chem’s PM background concentrations exceed the World Health Organization (WHO) guidelines over the entire ME. Mineral dust is the major contributor to PM (≈ 75%–95%) compared to other aerosol types. Near and downwind from the SO2 emission sources, non-dust aerosols (primarily sulfate) contribute up to 30% to PM(sub 2.5). The contribution of sea salt to PM in coastal regions can reach 5%. The contributions of organic matter, black carbon and organic carbon to PM over the Middle East are insignificant. In the major cities over the Arabian Peninsula, the 90th percentile of PM(sub 10) and PM(sub 2.5) (particles with diameters less than 10 and 2.5 μm, respectively) daily mean surface concentrations exceed the corresponding Kingdom of Saudi Arabia air quality limits. The contribution of the non-dust component to PM(sub 2.5) is < 25%, which limits the emission control effect on air quality. The mitigation of the dust effect on air quality requires the development of environment-based approaches like growing tree belts around the cities and enhancing in-city vegetation cover. The WRF-Chem configuration presented in this study could be a prototype of a future air quality forecast system that warns the population against air pollution hazards.

Alexander Ukhov↗