Time-dependent entrainment of smoke presents an observational challenge for assessing aerosol–cloud interactions over the southeast Atlantic Ocean

The colocation of clouds and smoke over the southeast Atlantic Ocean during the southern African biomass burning season has numerous radiative implications, including microphysical modulation of the clouds if smoke is entrained into the marine boundary layer. NASA's ObseRvations of Aerosols above CLouds and their intEractionS (ORACLES) campaign is studying this system with aircraft in three field deployments between 2016 and 2018. Results from ORACLES-2016 show that the relationship between cloud droplet number concentration and smoke below cloud is consistent with previously reported values, whereas cloud droplet number concentration is only weakly associated with smoke immediately above cloud at the time of observation. By combining field observations, regional chemistry–climate modeling, and theoretical boundary layer aerosol budget equations, we show that the history of smoke entrainment (which has a characteristic mixing timescale on the order of days) helps explain variations in cloud properties for similar instantaneous above-cloud smoke environments. Precipitation processes can obscure the relationship between above-cloud smoke and cloud properties in parts of the southeast Atlantic, but marine boundary layer carbon monoxide concentrations for two case study flights suggest that smoke entrainment history drove the observed differences in cloud properties for those days. A Lagrangian framework following the clouds and accounting for the history of smoke entrainment and precipitation is likely necessary for quantitatively studying this system; an Eulerian framework (e.g., instantaneous correlation of A-train satellite observations) is unlikely to capture the true extent of smoke–cloud interaction in the southeast Atlantic.</p

Use of lidar aerosol extinction and backscatter coefficients to estimate cloud condensation nuclei (CCN) concentrations in the southeast Atlantic

Accurately capturing cloud condensation nuclei (CCN) concentrations is key to understanding the aerosol–cloud interactions that continue to feature the highest uncertainty amongst numerous climate forcings. In situ CCN observations are sparse, and most non-polarimetric passive remote sensing techniques are limited to providing column-effective CCN proxies such as total aerosol optical depth (AOD). Lidar measurements, on the other hand, resolve profiles of aerosol extinction and/or backscatter coefficients that are better suited for constraining vertically resolved aerosol optical and microphysical properties. Here we present relationships between aerosol backscatter and extinction coefficients measured by the airborne High Spectral Resolution Lidar 2 (HSRL-2) and in situ measurements of CCN concentrations. The data were obtained during three deployments in the NASA ObseRvations of Aerosols above CLouds and their intEractionS (ORACLES) project, which took place over the southeast Atlantic (SEA) during September 2016, August 2017, and September–October 2018. Our analysis of spatiotemporally collocated in situ CCN concentrations and HSRL-2 measurements indicates strong linear relationships between both data sets. The correlation is strongest for supersaturations (S) greater than 0.25 % and dry ambient conditions above the stratocumulus deck, where relative humidity (RH) is less than 50 %. We find CCN–HSRL-2 Pearson correlation coefficients between 0.95–0.97 for different parts of the seasonal burning cycle that suggest fundamental similarities in biomass burning aerosol (BBA) microphysical properties. We find that ORACLES campaign-average values of in situ CCN and in situ extinction coefficients are qualitatively similar to those from other regions and aerosol types, demonstrating overall representativeness of our data set. We compute CCN–backscatter and CCN–extinction regressions that can be used to resolve vertical CCN concentrations across entire above-cloud lidar curtains. These lidar-derived CCN concentrations can be used to evaluate model performance, which we illustrate using an example CCN concentration curtain from the Weather Research and Forecasting Model coupled with physics packages from the Community Atmosphere Model version 5 (WRF-CAM5). These results demonstrate the utility of deriving vertically resolved CCN concentrations from lidar observations to expand the spatiotemporal coverage of limited or unavailable in situ observations.</p

Data assimilation in atmospheric chemistry models: current status and future prospects for coupled chemistry meteorology models

Abstract. Data assimilation is used in atmospheric chemistry models to improve air quality forecasts, construct re-analyses of three-dimensional chemical (including aerosol) concentrations and perform inverse modeling of input variables or model parameters (e.g., emissions). Coupled chemistry meteorology models (CCMM) are atmospheric chemistry models that simulate meteorological processes and chemical transformations jointly. They offer the possibility to assimilate both meteorological and chemical data; however, because CCMM are fairly recent, data assimilation in CCMM has been limited to date. We review here the current status of data assimilation in atmospheric chemistry models with a particular focus on future prospects for data assimilation in CCMM. We first review the methods available for data assimilation in atmospheric models, including variational methods, ensemble Kalman filters, and hybrid methods. Next, we review past applications that have included chemical data assimilation in chemical transport models (CTM) and in CCMM. Observational data sets available for chemical data assimilation are described, including surface data, surface-based remote sensing, airborne data, and satellite data. Several case studies of chemical data assimilation in CCMM are presented to highlight the benefits obtained by assimilating chemical data in CCMM. A case study of data assimilation to constrain emissions is also presented. There are few examples to date of joint meteorological and chemical data assimilation in CCMM and potential difficulties associated with data assimilation in CCMM are discussed. As the number of variables being assimilated increases, it is essential to characterize correctly the errors; in particular, the specification of error cross-correlations may be problematic. In some cases, offline diagnostics are necessary to ensure that data assimilation can truly improve model performance. However, the main challenge is likely to be the paucity of chemical data available for assimilation in CCMM

Diurnal variation of aerosol optical depth and PM_2.5 in South Korea: a synthesis from AERONET, satellite (GOCI), KORUS-AQ observation, and the WRF-Chem model

Spatial distribution of diurnal variations of aerosol properties in South Korea, both long term and short term, is studied by using 9 AERONET (AErosol RObotic NETwork) sites from 1999 to 2017 and an additional 10 sites during the KORUS-AQ (Korea–United States Air Quality) field campaign in May and June of 2016. The extent to which the WRF-Chem (Weather Research and Forecasting coupled with Chemistry) model and the GOCI (Geostationary Ocean Color Imager) satellite retrieval can describe these variations is also analyzed. On a daily average, aerosol optical depth (AOD) at 550&thinsp;nm is 0.386 and shows a diurnal variation of 20 to −30&thinsp;% in inland sites, which is larger than the AOD of 0.308 and diurnal variation of ±20&thinsp;% seen in coastal sites. For all the inland and coastal sites, AERONET, GOCI, and WRF-Chem, and observed PM2.5 (particulate matter with aerodynamic diameter less than 2.5&thinsp;µm) data generally show dual peaks for both AOD and PM2.5, one in the morning (often at  ∼ 08:00–10:00&thinsp;KST, Korea Standard Time, especially for PM2.5) and another in the early afternoon ( ∼ 14:00&thinsp;KST, albeit for PM2.5 this peak is smaller and sometimes insignificant). In contrast, Ångström exponent values in all sites are between 1.2 and 1.4 with the exception of the inland rural sites having smaller values near 1.0 during the early morning hours. All inland sites experience a pronounced increase in the Ångström exponent from morning to evening, reflecting an overall decrease in particle size in daytime. To statistically obtain the climatology of diurnal variation of AOD, a minimum requirement of  ∼ 2 years of observation is needed in coastal rural sites, twice as long as that required for the urban sites, which suggests that the diurnal variation of AOD in an urban setting is more distinct and persistent. While Korean GOCI satellite retrievals are able to consistently capture the diurnal variation of AOD (although it has a systematically low bias of 0.04 on average and up to 0.09 in later afternoon hours), WRF-Chem clearly has a deficiency in describing the relative change of peaks and variations between the morning and afternoon, suggesting further studies for the diurnal profile of emissions. Furthermore, the ratio between PM2.5 and AOD in WRF-Chem is persistently larger than the observed counterparts by 30&thinsp;%–50&thinsp;% in different sites, but spatially no consistent diurnal variation pattern of this ratio can be found. Overall, the relatively small diurnal variation of PM2.5 is in high contrast with large AOD diurnal variation, which suggests the large diurnal variation of AOD–PM2.5 relationships (with the PM2.5&thinsp;∕&thinsp;AOD ratio being largest in the early morning, decreasing around noon, and increasing in late afternoon) and, therefore, the need to use AOD from geostationary satellites to constrain either modeling or estimate of surface PM2.5 for air quality application.</p

Surface Dimming by the 2013 Rim Fire Simulated by a Sectional Aerosol Model

The Rim Fire of 2013, the third largest area burned by fire recorded in California history, is simulated by a climate model coupled with a size-resolved aerosol model. Modeled aerosol mass, number and particle size distribution are within variability of data obtained from multiple airborne in-situ measurements. Simulations suggest Rim Fire smoke may block 4-6 of sunlight energy reaching the surface, with a dimming efficiency around 120-150 W m(exp -2) per unit aerosol optical depth in the mid-visible at 13:00-15:00 local time. Underestimation of simulated smoke single scattering albedo at mid-visible by 0.04 suggests the model overestimates either the particle size or the absorption due to black carbon. This study shows that exceptional events like the 2013 Rim Fire can be simulated by a climate model with one-degree resolution with overall good skill, though that resolution is still not sufficient to resolve the smoke peak near the source region

Understanding and improving model representation of aerosol optical properties for a Chinese haze event measured during KORUS-AQ

KORUS-AQ was an international cooperative air quality field study in South Korea that measured local and remote sources of air pollution affecting the Korean Peninsula during May–June 2016. Some of the largest aerosol mass concentrations were measured during a Chinese haze transport event (24 May). Air quality forecasts using the WRF-Chem model with aerosol optical depth (AOD) data assimilation captured AOD during this pollution episode but overpredicted surface particulate matter concentrations in South Korea, especially PM2.5, often by a factor of 2 or larger. Analysis revealed multiple sources of model deficiency related to the calculation of optical properties from aerosol mass that explain these discrepancies. Using in situ observations of aerosol size and composition as inputs to the optical properties calculations showed that using a low-resolution size bin representation (four bins) underestimates the efficiency with which aerosols scatter and absorb light (mass extinction efficiency). Besides using finer-resolution size bins (8–16 bins), it was also necessary to increase the refractive indices and hygroscopicity of select aerosol species within the range of values reported in the literature to achieve better consistency with measured values of the mass extinction efficiency (6.7 m2 g−1 observed average) and light-scattering enhancement factor (f(RH)) due to aerosol hygroscopic growth (2.2 observed average). Furthermore, an evaluation of the optical properties obtained using modeled aerosol properties revealed the inability of sectional and modal aerosol representations in WRF-Chem to properly reproduce the observed size distribution, with the models displaying a much wider accumulation mode. Other model deficiencies included an underestimate of organic aerosol density (1.0 g cm−3 in the model vs. observed average of 1.5 g cm−3) and an overprediction of the fractional contribution of submicron inorganic aerosols other than sulfate, ammonium, nitrate, chloride, and sodium corresponding to mostly dust (17 %–28 % modeled vs. 12 % estimated from observations). These results illustrate the complexity of achieving an accurate model representation of optical properties and provide potential solutions that are relevant to multiple disciplines and applications such as air quality forecasts, health impact assessments, climate projections, solar power forecasts, and aerosol data assimilation

Revealing important nocturnal and day-to-day variations in fire smoke emissions through a multiplatform inversion

We couple airborne, ground-based, and satellite observations; conduct regional simulations; and develop and apply an inversion technique to constrain hourly smoke emissions from the Rim Fire, the third largest observed in California, USA. Emissions constrainedwithmultiplatform data show notable nocturnal enhancements (sometimes over a factor of 20), correlate better with daily burned area data, and are a factor of 2–4 higher than a priori estimates, highlighting the need for improved characterization of diurnal profiles and day-to-day variability when modeling extreme fires. Constraining only with satellite data results in smaller enhancements mainly due to missing retrievals near the emissions source, suggesting that top-down emission estimates for these events could be underestimated and a multi-platform approach is required to resolve them. Predictions driven by emissions constrained with multi-platform data present significant variations in downwind air quality and in aerosol feedback on meteorology, emphasizing the need for improved emissions estimates during exceptional events

Revealing important nocturnal and day-to-day variations in fire smoke emissions through a multiplatform inversion

We couple airborne, ground-based, and satellite observations; conduct regional simulations; and develop and apply an inversion technique to constrain hourly smoke emissions from the Rim Fire, the third largest observed in California, USA. Emissions constrainedwithmultiplatform data show notable nocturnal enhancements (sometimes over a factor of 20), correlate better with daily burned area data, and are a factor of 2–4 higher than a priori estimates, highlighting the need for improved characterization of diurnal profiles and day-to-day variability when modeling extreme fires. Constraining only with satellite data results in smaller enhancements mainly due to missing retrievals near the emissions source, suggesting that top-down emission estimates for these events could be underestimated and a multi-platform approach is required to resolve them. Predictions driven by emissions constrained with multi-platform data present significant variations in downwind air quality and in aerosol feedback on meteorology, emphasizing the need for improved emissions estimates during exceptional events

On the differences in the vertical distribution of modeled aerosol optical depth over the southeastern Atlantic

The southeastern Atlantic is home to an expansive smoke aerosol plume overlying a large cloud deck for approximately a third of the year. The aerosol plume is mainly attributed to the extensive biomass burning activities that occur in southern Africa. Current Earth system models (ESMs) reveal significant differences in their estimates of regional aerosol radiative effects over this region. Such large differences partially stem from uncertainties in the vertical distribution of aerosols in the troposphere. These uncertainties translate into different aerosol optical depths (AODs) in the planetary boundary layer (PBL) and the free troposphere (FT). This study examines differences of AOD fraction in the FT and AOD differences among ESMs (WRF-CAM5, WRF-FINN, GEOS-Chem, EAM-E3SM, ALADIN, GEOS-FP, and MERRA-2) and aircraft-based measurements from the NASA ObseRvations of Aerosols above CLouds and their intEractionS (ORACLES) field campaign. Models frequently define the PBL as the well-mixed surface-based layer, but this definition misses the upper parts of decoupled PBLs, in which most low-level clouds occur. To account for the presence of decoupled boundary layers in the models, the height of maximum vertical gradient of specific humidity profiles from each model is used to define PBL heights. Results indicate that the monthly mean contribution of AOD in the FT to the total-column AOD ranges from 44 % to 74 % in September 2016 and from 54 % to 71 % in August 2017 within the region bounded by 25∘ S–0∘ N–S and 15∘ W–15∘ E (excluding land) among the ESMs. ALADIN and GEOS-Chem show similar aerosol plume patterns to a derived above-cloud aerosol product from the Moderate Resolution Imaging Spectroradiometer (MODIS) during September 2016, but none of the models show a similar above-cloud plume pattern to MODIS in August 2017. Using the second-generation High Spectral Resolution Lidar (HSRL-2) to derive an aircraft-based constraint on the AOD and the fractional AOD, we found that WRF-CAM5 produces 40 % less AOD than those from the HSRL-2 measurements, but it performs well at separating AOD fraction between the FT and the PBL. AOD fractions in the FT for GEOS-Chem and EAM-E3SM are, respectively, 10 % and 15 % lower than the AOD fractions from the HSRL-2. Their similar mean AODs reflect a cancellation of high and low AOD biases. Compared with aircraft-based observations, GEOS-FP, MERRA-2, and ALADIN produce 24 %–36 % less AOD and tend to misplace more aerosols in the PBL. The models generally underestimate AODs for measured AODs that are above 0.8, indicating their limitations at reproducing high AODs. The differences in the absolute AOD, FT AOD, and the vertical apportioning of AOD in different models highlight the need to continue improving the accuracy of modeled AOD distributions. These differences affect the sign and magnitude of the net aerosol radiative forcing, especially when aerosols are in contact with clouds.</p