Parameterization of Aerosol Sinks in Chemical Transport Models

The modelers point of view is that the aerosol problem is one of sources, evolution, and sinks. Relative to evolution and sink processes, enormous attention is given to the problem of aerosols sources, whether inventory based (e.g., fossil fuel emissions) or dynamic (e.g., dust, sea salt, biomass burning). On the other hand, aerosol losses in models are a major factor in controlling the aerosol distribution and lifetime. Here we shine some light on how aerosol sinks are treated in modern chemical transport models. We discuss the mechanisms of dry and wet loss processes and the parameterizations for those processes in a single model (GEOS-5). We survey the literature of other modeling studies. We additionally compare the budgets of aerosol losses in several of the ICAP models

Near Real-Time Sub/Seasonal Prediction of Aerosol and Air Quality at the NASA Global Modeling and Assimilation Office

Version 2 of the coupled modeling and analysis system used to produce near real time subseasonal to seasonal forecasts was released almost two years ago by the NASA/Goddard Global Modeling and Assimilation Office. The model runs at approximately 1/2 degree globally in the atmosphere and ocean, contains a realistic description of the cryosphere, and includes an interactive aerosol model. The data assimilation used to produce initial conditions is weakly coupled, in which the atmosphere-only assimilated state is coupled to an ocean data assimilation system using a Local Ensemble Transform Kalman Filter. Results of aerosol-derived air quality (Particulate Matter) from an extensive series of retrospective forecasts will be shown, with particular focus on the continental United States and eastern Asia. In addition, under some circumstances, the interactive aerosol is shown to improve seasonal time scale prediction skill. Plans for a future version of the system with predicted biomass burning from fires will also be discussed

Improvement of Stratospheric Aerosol Extinction Retrieval from OMPS/LP Using a New Aerosol Model

The Ozone Mapping and Profiler Suite Limb Profiler (OMPS/LP) has been flying on the Suomi National Polar-orbiting Partnership (S-NPP) satellite since October 2011. It is designed to produce ozone and aerosol vertical profiles at 2km vertical resolution over the entire sunlit globe. Aerosol extinction profiles are computed with Mie theory using radiances measured at 675nm. The operational Version 1.0 (V1.0) aerosol extinction retrieval algorithm assumes a bimodal lognormal aerosol size distribution (ASD) whose parameters were derived by combining an in situ measurement of aerosol microphysics with the Stratospheric Aerosol and Gas Experiment (SAGE II) aerosol extinction climatology. Internal analysis indicates that this bimodal lognormal ASD does not sufficiently explain the spectral dependence of LP-measured radiances. In this paper we describe the derivation of an improved aerosol size distribution, designated Version 1.5 (V1.5), for the LP retrieval algorithm. The new ASD uses a gamma function distribution that is derived from Community Aerosol and Radiation Model for Atmospheres (CARMA)-calculated results. A cumulative distribution fit derived from the gamma function ASD gives better agreement with CARMA results at small particle radii than bimodal or unimodal functions. The new ASD also explains the spectral dependence of LP-measured radiances better than the V1.0 ASD. We find that the impact of our choice of ASD on the retrieved extinctions varies strongly with the underlying reflectivity of the scene. Initial comparisons with collocated extinction profiles retrieved at 676nm from the SAGE III instrument on the International Space Station (ISS) show a significant improvement in agreement for the LP V1.5 retrievals. Zonal mean extinction profiles agree to within 10% between 19 and 29km, and regression fits of collocated samples show improved correlation and reduced scatter compared to the V1.0 product. This improved agreement will motivate development of more sophisticated ASDs from CARMA results that incorporate latitude, altitude and seasonal variations in aerosol properties

Principle Component Analysis of the Evolution of the Saharan Air Layer and Dust Transport: Comparisons between a Model Simulation and MODIS Retrievals

The onset and evolution of Saharan Air Layer (SAL) episodes during June-September 2002 are diagnosed by applying principal component analysis to the NCEP reanalysis temperature anomalies at 850 hPa, where the largest SAL-induced temperature anomalies are located. The first principal component (PC) represents the onset of SAL episodes, which are associated with large warm anomalies located at the west coast of Africa. The second PC represents two opposite phases of the evolution of the SAL. The positive phase of the second PC corresponds to the southwestward extension of the warm anomalies into the tropical-subtropical North Atlantic Ocean, and the negative phase corresponds to the northwestward extension into the subtropical to mid-latitude North Atlantic Ocean and the southwest Europe. A dust transport model (CARMA) and the MODIS retrievals are used to study the associated effects on dust distribution and deposition. The positive (negative) phase of the second PC corresponds to a strengthening (weakening) of the offshore flows in the lower troposphere around 10deg - 20degN, causing more (less) dust being transported along the tropical to subtropical North Atlantic Ocean. The variation of the offshore flow indicates that the subseasonal variation of African Easterly Jet is associated with the evolution of the SAL. Significant correlation is found between the second PC time series and the daily West African monsoon index, implying a dynamical linkage between West African monsoon and the evolution of the SAL and Saharan dust transport

Effects of Spatial Resolution on the Simulated Dust Aerosol Lifecycle: Implications for Dust Event Magnitude and Timing in the NASA GEOS-5 AGCM

The NASA GEOS-5 atmospheric transport model simulates global aerosol distributions with an online aerosol module. GEOS-5 may be run at various horizontal spatial resolutions depending on the research application. For example, long integration climate simulations are typically run at 2 deg or 1 deg grid spacing, whereas aerosol reanalysis and forecasting applications may be performed at O.5 deg or 0.25 deg resolutions. In this study, we assess the implications of varying spatial resolution on the simulated aerosol fields, with a particular focus on dust. Dust emissions in GEOS-5 are calculated with one of two parameterizations, one based on the Goddard Chemistry, Aerosol, Radiation, and Transport (GO CART) model and another based on the Dust Entrainment and Deposition (DEAD) model. Emission fluxes are parameterized in terms of the surface wind speed, either the 10-m (GO CART) or friction (DEAD) wind speed. We consider how surface wind speeds and thus the dust emission rates are a function of the model spatial resolution. We find that spatial resolution has a significant effect on the magnitude of dust emissions, as higher resolution versions of the model have typically higher surface wind speeds. Utilizing space-borne observations from MISR, MODIS, and CALIOP, we find that simulated Aerosol Optical Thickness (AOT) distributions respond differently to spatial resolution over the African and Asian source regions, highlighting the need to regional dust emission tuning. When compared to ground-based observations from AERONET, we found improved timing of dust events with as spatial resolution was increased. In an attempt to improve the representation of the dust aerosol lifecycle at coarse resolutions, we found that incorporating the effects of sub-grid wind variability in a course resolution simulation led to improved agreement with observed AOT magnitudes, but did not impact the timing of simulated dust events

Differences in the Vertical and Microphysical Evolution of Volcanic and Pyrocumulonimbus Stratospheric Aerosol Plumes as Observed by CALIOP and CATS Satellite Lidar

For some time, volcanic eruptions have been thought to be the only significant direct injectors of aerosols in the stratosphere. However, recent fire seasons have featured fire events resulting in large volcanic-sized pyrocumulonimbus plumes of smoke aerosols reaching many kilometers into the lower stratosphere. To understand and model the effects of these pyrocumulonimbus events on stratospheric composition and climate, a natural analogy lies with better studied volcanic events; however, differences in plume composition may limit this comparison. Using satellite lidar from both CATS and CALIOP, we show that the stratospheric aerosol plumes from the record-setting Pacific Northwest pyrocumulonimbus event of 2017 and the Calbuco volcanic eruption of 2015 evolve differently both vertically and microphysically. Specifically, depolarization ratios indicate that this pyrocumulonimbus events aerosol particles became more irregularly shaped over time in contrast to volcanic aerosols which become more spherical over time. Accounting for these changes in aerosol properties may be significant in assessing the effects of pyrocumulonimbus events on the Earths radiative balance and aid in refining stratospheric aerosol typing algorithms to differentiate volcanic from pyrocumulonimbus plumes

The Fate of Saharan Dust Across the Atlantic and Implications for a Central American Dust Barrier

Saharan dust was observed over the Caribbean basin during the summer 2007 NASA Tropical Composition, Cloud, and Climate Coupling (TC4) field experiment. Airborne Cloud Physics Lidar (CPL) and satellite observations from MODIS suggest a barrier to dust transport across Central America into the eastern Pacific. We use the NASA GEOS-5 atmospheric transport model with online aerosol tracers to perform simulations of the TC4 time period in order to understand the nature of this barrier. Our simulations are driven by the Modem Era Retrospective-Analysis for Research and Applications (MERRA) meteorological analyses. We evaluate our baseline simulated dust distributions using MODIS and CALIOP satellite and ground-based AERONET sun photometer observations. GEOS-5 reproduces the observed location, magnitude, and timing of major dust events, but our baseline simulation does not develop as strong a barrier to dust transport across Central America as observations suggest. Analysis of the dust transport dynamics and lost processes suggest that while both mechanisms play a role in defining the dust transport barrier, loss processes by wet removal of dust are about twice as important as transport. Sensitivity analyses with our model showed that the dust barrier would not exist without convective scavenging over the Caribbean. The best agreement between our model and the observations was obtained when dust wet removal was parameterized to be more aggressive, treating the dust as we do hydrophilic aerosols

Workshop Summary: International Cooperative for Aerosol Prediction Workshop On Aerosol Forecast Verification

The purpose of this workshop was to reinforce the working partnership between centers who are actively involved in global aerosol forecasting, and to discuss issues related to forecast verification. Participants included representatives from operational centers with global aerosol forecasting requirements, a panel of experts on Numerical Weather Prediction and Air Quality forecast verification, data providers, and several observers from the research community. The presentations centered on a review of current NWP and AQ practices with subsequent discussion focused on the challenges in defining appropriate verification measures for the next generation of aerosol forecast systems

Quantifying the Aerosol Semi-Direct Effect in the NASA GEOS-5 AGCM

Aerosols such as black carbon, dust, and some organic carbon species both scatter and absorb incoming solar radiation. This direct aerosol radiative forcing (DARF) redistributes solar energy both by cooling the surface and warming the atmosphere. As a result, these aerosols affect atmospheric stability and cloud cover (the semi-direct effect, or SDE). Furthermore, in regions with persistent high loadings of absorbing aerosols (e.g. Asia), regional circulation patterns may be altered, potentially resulting in changes in precipitation patterns. Here we investigate aerosol-climate coupling using the NASA Goddard Earth Observing System model version 5 (GEOS-5) atmospheric general circulation model (AGCM), in which we have implemented an online version of the Goddard Chemistry, Aerosol, Radiation and Transport (GOCART) model. GOCART includes representations of the sources, sinks, and chemical transformation of externally mixed dust, sea salt, sulfate, and carbonaceous aerosols. We examine a series of free-running ensemble climate simulations of the present-day period (2000-2009) forced by observed sea surface temperatures to determine the impact of aerosols on the model climate. The SDE and response of each simulation is determined by differencing with respect to the control simulation (no aerosol forcing). In a free-running model, any estimate of the SDE includes changes in clouds due both to atmospheric heating from aerosols and changes in circulation. To try and quantify the SDE without these circulation changes we then examine the DARF and SDE in GEOS-5 with prescribed meteorological analyses introduced by the MERRA analysis. By doing so, we are able to examine changes in model clouds that occur on shorter scales (six hours). In the GEOS-5 data assimilation system (DAS), the analysis is defined as the best estimate of the atmospheric state at any given time, and it is determined by optimally combining a first-guess short-term GCM forecast with all available observations. The Incremental Analysis Update (IAU) is added to the model forecast tendencies to align them with the analysis every six hours, thus preventing longer timescale feedbacks due to the aerosol forcing. We calculate the SDE by comparing model runs with and without aerosols, and the difference in the IAU between these runs is a useful metric with which to evaluate the impact of the SDE on the model atmosphere and clouds. Decreasing the IAU indicates that the aerosol direct and semi-direct effects act to reduce the bias between the model and observations and vice versa

Impact of Satellite Viewing-Swath Width on Global and Regional Aerosol Optical Thickness Statistics and Trends

We use the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite aerosol optical thickness (AOT) product to assess the impact of reduced swath width on global and regional AOT statistics and trends. Alongtrack and across-track sampling strategies are employed, in which the full MODIS data set is sub-sampled with various narrow-swath (approximately 400-800 km) and single pixel width (approximately 10 km) configurations. Although view-angle artifacts in the MODIS AOT retrieval confound direct comparisons between averages derived from different sub-samples, careful analysis shows that with many portions of the Earth essentially unobserved, spatial sampling introduces uncertainty in the derived seasonal-regional mean AOT. These AOT spatial sampling artifacts comprise up to 60%of the full-swath AOT value under moderate aerosol loading, and can be as large as 0.1 in some regions under high aerosol loading. Compared to full-swath observations, narrower swath and single pixel width sampling exhibits a reduced ability to detect AOT trends with statistical significance. On the other hand, estimates of the global, annual mean AOT do not vary significantly from the full-swath values as spatial sampling is reduced. Aggregation of the MODIS data at coarse grid scales (10 deg) shows consistency in the aerosol trends across sampling strategies, with increased statistical confidence, but quantitative errors in the derived trends are found even for the full-swath data when compared to high spatial resolution (0.5 deg) aggregations. Using results of a model-derived aerosol reanalysis, we find consistency in our conclusions about a seasonal-regional spatial sampling artifact in AOT Furthermore, the model shows that reduced spatial sampling can amount to uncertainty in computed shortwave top-ofatmosphere aerosol radiative forcing of 2-3 W m(sup2). These artifacts are lower bounds, as possibly other unconsidered sampling strategies would perform less well. These results suggest that future aerosol satellite missions having significantly less than full-swath viewing are unlikely to sample the true AOT distribution well enough to obtain the statistics needed to reduce uncertainty in aerosol direct forcing of climate