While stand alone satellite and model aerosol products see wide utilization, there is a significant need in numerous atmospheric and climate applications for a fused product on a regular grid. Aerosol data assimilation is an operational reality at numerous centers, and like meteorological reanalyses, aerosol reanalyses will see significant use in the near future. Here we present a standardized 2003–2013 global 1 × 1° and 6-hourly modal aerosol optical thickness (AOT) reanalysis product. This data set can be applied to basic and applied Earth system science studies of significant aerosol events, aerosol impacts on numerical weather prediction, and electro-optical propagation and sensor performance, among other uses. This paper describes the science of how to develop and score an aerosol reanalysis product. This reanalysis utilizes a modified Navy Aerosol Analysis and Prediction System (NAAPS) at its core and assimilates quality controlled retrievals of AOT from the Moderate Resolution Imaging Spectroradiometer (MODIS) on Terra and Aqua and the Multi-angle Imaging SpectroRadiometer (MISR) on Terra. The aerosol source functions, including dust and smoke, were regionally tuned to obtain the best match between the model fine- and coarse-mode AOTs and the Aerosol Robotic Network (AERONET) AOTs. Other model processes, including deposition, were tuned to minimize the AOT difference between the model and satellite AOT. Aerosol wet deposition in the tropics is driven with satellite-retrieved precipitation, rather than the model field. The final reanalyzed fine- and coarse-mode AOT at 550 nm is shown to have good agreement with AERONET observations, with global mean root mean square error around 0.1 for both fine- and coarse-mode AOTs. This paper includes a discussion of issues particular to aerosol reanalyses that make them distinct from standard meteorological reanalyses, considerations for extending such a reanalysis outside of the NASA A-Train era, and examples of how the aerosol reanalysis can be applied or fused with other model or remote sensing products. Finally, the reanalysis is evaluated in comparison with other available studies of aerosol trends, and the implications of this comparison are discussed.

This study examines major wildfires in the western United States between 2005 and 2015 to determine which species exhibit the highest percent change in mass concentration on day of peak fire influence relative to preceding nonfire days. Forty‐one fires were examined using the Environmental Protection Agency (EPA) Interagency Monitoring of Protected Visual Environments (IMPROVE) data set. Organic carbon (OC) and elemental carbon (EC) constituents exhibited the highest percent change increase. The sharpest enhancements were for the volatile (OC1) and semivolatile (OC2) OC fractions, suggestive of secondary organic aerosol formation during plume transport. Of the noncarbonaceous constituents, Cl, P, K, NO₃⁻, and Zn levels exhibited the highest percent change. Dust was significantly enhanced in wildfire plumes, based on significant enhancements in fine soil components (i.e., Si, Ca, Al, Fe, and Ti) and PM_coarse (i.e., PM₁₀–PM_2.5). A case study emphasized how transport of wildfire plumes significantly impacted downwind states, with higher levels of fine soil and PM_coarse at the downwind state (Arizona) as compared to the source of the fires (California). A global model (Navy Aerosol Analysis and Prediction System, NAAPS) did not capture the dust influence over California or Arizona during this case event because it is not designed to resolve dust dynamics in fires, which motivates improved treatment of such processes. Significant chloride depletion was observed on the peak EC day for almost a half of the fires examined. Size‐resolved measurements during two specific fires at a coastal California site revealed significant chloride reductions for particle aerodynamic diameters between 1 and 10 μm.

Since the first International Cooperative for Aerosol Prediction (ICAP) multi‐model ensemble (MME) study, the number of ICAP global operational aerosol models has increased from five to nine. An update of the current ICAP status is provided, along with an evaluation of the performance of ICAP‐MME over 2012–2017, with a focus on June 2016–May 2017. Evaluated with ground‐based Aerosol Robotic Network (AERONET) aerosol optical depth (AOD) and data assimilation quality MODerate‐resolution Imaging Spectroradiometer (MODIS) retrieval products, the ICAP‐MME AOD consensus remains the overall top‐scoring and most consistent performer among all models in terms of root‐mean‐square error (RMSE), bias and correlation for total, fine‐ and coarse‐mode AODs as well as dust AOD; this is similar to the first ICAP‐MME study. Further, over the years, the performance of ICAP‐MME is relatively stable and reliable compared to more variability in the individual models. The extent to which the AOD forecast error of ICAP‐MME can be predicted is also examined. Leading predictors are found to be the consensus mean and spread. Regression models of absolute forecast errors were built for AOD forecasts of different lengths for potential applications. ICAP‐MME performance in terms of modal AOD RMSEs of the 21 regionally representative sites over 2012–2017 suggests a general tendency for model improvements in fine‐mode AOD, especially over Asia. No significant improvement in coarse‐mode AOD is found overall for this time period.

We present an Arctic aerosol optical depth (AOD) climatology and trend analysis for 2003–2019 spring and summertime periods derived from a combination of multi-agency aerosol reanalyses, remote-sensing retrievals, and ground observations. This includes the U.S. Navy Aerosol Analysis and Prediction System ReAnalysis version 1 (NAAPS-RA v1), the NASA Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2), and the Copernicus Atmosphere Monitoring Service ReAnalysis (CAMSRA). Spaceborne remote-sensing retrievals of AOD are considered from the Moderate Resolution Imaging Spectroradiometer (MODIS), the Multi-angle Imaging SpectroRadiometer (MISR), and the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP). Ground-based data include sun photometer data from AErosol RObotic NETwork (AERONET) sites and oceanic Maritime Aerosol Network (MAN) measurements. Aerosol reanalysis AODs and spaceborne retrievals show consistent climatological spatial patterns and trends for both spring and summer seasons over the lower Arctic (60–70^∘ N). Consistent AOD trends are also found for the high Arctic (north of 70^∘ N) from reanalyses. The aerosol reanalyses yield more consistent AOD results than climate models, can be verified well with AERONET, and corroborate complementary climatological and trend analysis. Speciated AODs are more variable than total AOD among the three reanalyses and a little more so for March–May (MAM) than for June–August (JJA). Black carbon (BC) AOD in the Arctic comes predominantly from biomass burning (BB) sources in both MAM and JJA, and BB overwhelms anthropogenic sources in JJA for the study period. AOD exhibits a multi-year negative MAM trend and a positive JJA trend in the Arctic during 2003–2019, due to an overall decrease in sulfate/anthropogenic pollution and a significant JJA increase in BB smoke. Interannual Arctic AOD variability is significantly large, driven by fine-mode and, specifically, BB smoke, with both smoke contribution and interannual variation larger in JJA than in MAM. It is recommended that climate models should account for BB emissions and BB interannual variabilities and trends in Arctic climate change studies.

An ensemble-based forecast and data assimilation system has been developed for use in Navy aerosol forecasting. The system makes use of an ensemble of the Navy Aerosol Analysis Prediction System (ENAAPS) at 1 × 1°, combined with an ensemble adjustment Kalman filter from NCAR's Data Assimilation Research Testbed (DART). The base ENAAPS-DART system discussed in this work utilizes the Navy Operational Global Analysis Prediction System (NOGAPS) meteorological ensemble to drive offline NAAPS simulations coupled with the DART ensemble Kalman filter architecture to assimilate bias-corrected MODIS aerosol optical thickness (AOT) retrievals. This work outlines the optimization of the 20-member ensemble system, including consideration of meteorology and source-perturbed ensemble members as well as covariance inflation. Additional tests with 80 meteorological and source members were also performed. An important finding of this work is that an adaptive covariance inflation method, which has not been previously tested for aerosol applications, was found to perform better than a temporally and spatially constant covariance inflation. Problems were identified with the constant inflation in regions with limited observational coverage. The second major finding of this work is that combined meteorology and aerosol source ensembles are superior to either in isolation and that both are necessary to produce a robust system with sufficient spread in the ensemble members as well as realistic correlation fields for spreading observational information. The inclusion of aerosol source ensembles improves correlation fields for large aerosol source regions, such as smoke and dust in Africa, by statistically separating freshly emitted from transported aerosol species. However, the source ensembles have limited efficacy during long-range transport. Conversely, the meteorological ensemble generates sufficient spread at the synoptic scale to enable observational impact through the ensemble data assimilation. The optimized ensemble system was compared to the Navy's current operational aerosol forecasting system, which makes use of NAVDAS-AOD (NRL Atmospheric Variational Data Assimilation System for aerosol optical depth), a 2-D variational data assimilation system. Overall, the two systems had statistically insignificant differences in root-mean-squared error (RMSE), bias, and correlation relative to AERONET-observed AOT. However, the ensemble system is able to better capture sharp gradients in aerosol features compared to the 2DVar system, which has a tendency to smooth out aerosol events. Such skill is not easily observable in bulk metrics. Further, the ENAAPS-DART system will allow for new avenues of model development, such as more efficient lidar and surface station assimilation as well as adaptive source functions. At this early stage of development, the parity with the current variational system is encouraging.

Analysis of Sun photometer measured and satellite retrieved aerosol optical depth (AOD) data has shown that major aerosol pollution events with very high fine mode AOD (>1.0 in midvisible) in the China/Korea/Japan region are often observed to be associated with significant cloud cover. This makes remote sensing of these events difficult even for high temporal resolution Sun photometer measurements. Possible physical mechanisms for these events that have high AOD include a combination of aerosol humidification, cloud processing, and meteorological covariation with atmospheric stability and convergence. The new development of Aerosol Robotic Network Version 3 Level 2 AOD with improved cloud screening algorithms now allow for unprecedented ability to monitor these extreme fine mode pollution events. Further, the spectral deconvolution algorithm (SDA) applied to Level 1 data (L1; no cloud screening) provides an even more comprehensive assessment of fine mode AOD than L2 in current and previous data versions. Studying the 2012 winter‐summer period, comparisons of Aerosol Robotic Network L1 SDA daily average fine mode AOD data showed that Moderate Resolution Imaging Spectroradiometer satellite remote sensing of AOD often did not retrieve and/or identify some of the highest fine mode AOD events in this region. Also, compared to models that include data assimilation of satellite retrieved AOD, the L1 SDA fine mode AOD was significantly higher in magnitude, particularly for the highest AOD events that were often associated with significant cloudiness.