Abstract. This paper introduces and evaluates the second version of the global aerosol-climate model ECHAM-HAM. Major changes have been brought into the model, including new parameterizations for aerosol nucleation and water uptake, an explicit treatment of secondary organic aerosols, modified emission calculations for sea salt and mineral dust, the coupling of aerosol microphysics to a two-moment stratiform cloud microphysics scheme, and alternative wet scavenging parameterizations. These revisions extend the model's capability to represent details of the aerosol lifecycle and its interaction with climate. Nudged simulations of the year 2000 are carried out to compare the aerosol properties and global distribution in HAM1 and HAM2, and to evaluate them against various observations. Sensitivity experiments are performed to help identify the impact of each individual update in model formulation. Results indicate that from HAM1 to HAM2 there is a marked weakening of aerosol water uptake in the lower troposphere, reducing the total aerosol water burden from 75 Tg to 51 Tg. The main reason is the newly introduced κ-Köhler-theory-based water uptake scheme uses a lower value for the maximum relative humidity cutoff. Particulate organic matter loading in HAM2 is considerably higher in the upper troposphere, because the explicit treatment of secondary organic aerosols allows highly volatile oxidation products of the precursors to be vertically transported to regions of very low temperature and to form aerosols there. Sulfate, black carbon, particulate organic matter and mineral dust in HAM2 have longer lifetimes than in HAM1 because of weaker in-cloud scavenging, which is in turn related to lower autoconversion efficiency in the newly introduced two-moment cloud microphysics scheme. Modification in the sea salt emission scheme causes a significant increase in the ratio (from 1.6 to 7.7) between accumulation mode and coarse mode emission fluxes of aerosol number concentration. This leads to a general increase in the number concentration of smaller particles over the oceans in HAM2, as reflected by the higher Ångström parameters. Evaluation against observation reveals that in terms of model performance, main improvements in HAM2 include a marked decrease of the systematic negative bias in the absorption aerosol optical depth, as well as smaller biases over the oceans in Ångström parameter and in the accumulation mode number concentration. The simulated geographical distribution of aerosol optical depth (AOD) is better correlated with the MODIS data, while the surface aerosol mass concentrations are very similar to those in the old version. The total aerosol water content in HAM2 is considerably closer to the multi-model average from Phase I of the AeroCom intercomparison project. Model deficiencies that require further efforts in the future include (i) positive biases in AOD over the ocean, (ii) negative biases in AOD and aerosol mass concentration in high-latitude regions, and (iii) negative biases in particle number concentration, especially that of the Aitken mode, in the lower troposphere in heavily polluted regions.
Abstract. Aerosols have important impacts on air quality and climate, but the processes affecting their removal from the atmosphere are not fully understood and are poorly constrained by observations. This makes modelled aerosol lifetimes uncertain. In this study, we make use of an observational constraint on aerosol lifetimes provided by radionuclide measurements and investigate the causes of differences within a set of global models. During the Fukushima Dai-Ichi nuclear power plant accident of March 2011, the radioactive isotopes cesium-137 (137Cs) and xenon-133 (133Xe) were released in large quantities. Cesium attached to particles in the ambient air, approximately according to their available aerosol surface area. 137Cs size distribution measurements taken close to the power plant suggested that accumulation-mode (AM) sulfate aerosols were the main carriers of cesium. Hence, 137Cs can be used as a proxy tracer for the AM sulfate aerosol's fate in the atmosphere. In contrast, the noble gas 133Xe behaves almost like a passive transport tracer. Global surface measurements of the two radioactive isotopes taken over several months after the release allow the derivation of a lifetime of the carrier aerosol. We compare this to the lifetimes simulated by 19 different atmospheric transport models initialized with identical emissions of 137Cs that were assigned to an aerosol tracer with each model's default properties of AM sulfate, and 133Xe emissions that were assigned to a passive tracer. We investigate to what extent the modelled sulfate tracer can reproduce the measurements, especially with respect to the observed loss of aerosol mass with time. Modelled 137Cs and 133Xe concentrations sampled at the same location and times as station measurements allow a direct comparison between measured and modelled aerosol lifetime. The e-folding lifetime τe, calculated from station measurement data taken between 2 and 9 weeks after the start of the emissions, is 14.3 days (95 % confidence interval 13.1–15.7 days). The equivalent modelled τe lifetimes have a large spread, varying between 4.8 and 26.7 days with a model median of 9.4 ± 2.3 days, indicating too fast a removal in most models. Because sufficient measurement data were only available from about 2 weeks after the release, the estimated lifetimes apply to aerosols that have undergone long-range transport, i.e. not for freshly emitted aerosol. However, modelled instantaneous lifetimes show that the initial removal in the first 2 weeks was quicker (lifetimes between 1 and 5 days) due to the emissions occurring at low altitudes and co-location of the fresh plume with strong precipitation. Deviations between measured and modelled aerosol lifetimes are largest for the northernmost stations and at later time periods, suggesting that models do not transport enough of the aerosol towards the Arctic. The models underestimate passive tracer (133Xe) concentrations in the Arctic as well but to a smaller extent than for the aerosol (137Cs) tracer. This indicates that in addition to too fast an aerosol removal in the models, errors in simulated atmospheric transport towards the Arctic in most models also contribute to the underestimation of the Arctic aerosol concentrations.
Abstract. The uncertainties associated with the wet removal of aerosols entrained above convective cloud bases are investigated in a global aerosol-climate model (ECHAM5-HAM) under a set of limiting assumptions for the wet removal of the entrained aerosols. The limiting assumptions for the wet removal of entrained aerosols are negligible scavenging and vigorous scavenging (either through activation, with size-dependent impaction scavenging, or with the prescribed fractions of the standard model). To facilitate this process-based study, an explicit representation of cloud-droplet-borne and ice-crystal-borne aerosol mass and number, for the purpose of wet removal, is introduced into the ECHAM5-HAM model. This replaces and is compared with the prescribed cloud-droplet-borne and ice-crystal-borne aerosol fraction scavenging scheme of the standard model. A 20% to 35% uncertainty in simulated global, annual mean aerosol mass burdens and optical depth (AOD) is attributed to different assumptions for the wet removal of aerosols entrained above convective cloud bases. Assumptions about the removal of aerosols entrained above convective cloud bases control modeled upper tropospheric aerosol concentrations by as much as one order of magnitude. Simulated aerosols entrained above convective cloud bases contribute 20% to 50% of modeled global, annual mean aerosol mass convective wet deposition (about 5% to 10% of the total dry and wet deposition), depending on the aerosol species, when including wet scavenging of those entrained aerosols (either by activation, size-dependent impaction, or with the prescribed fraction scheme). Among the simulations, the prescribed fraction and size-dependent impaction schemes yield the largest global, annual mean aerosol mass convective wet deposition (by about two-fold). However, the prescribed fraction scheme has more vigorous convective mixed-phase wet removal (by two to five-fold relative to the size-dependent impaction scheme) since nearly all entrained accumulation and coarse mode aerosols are assumed to be cloud-droplet borne or ice-crystal borne, and evaporation due to the Bergeron-Findeisen process is neglected. The simulated convective wet scavenging of entrained accumulation and coarse mode aerosols has feedbacks on new particle formation and the number of Aitken mode aerosols, which control stratiform and convective cloud droplet number concentrations and yield precipitation changes in the ECHAM5-HAM model. However, the geographic distribution of aerosol annual mean convective wet deposition change in the model is driven by changes to the assumptions regarding the scavenging of aerosols entrained above cloud bases rather than by precipitation changes, except for sea salt deposition in the tropics. Uncertainty in the seasonal, regional cycles of AOD due to assumptions about entrained aerosol wet scavenging is similar in magnitude to the estimated error in the AOD retrievals. The uncertainty in aerosol concentrations, burdens, and AOD attributed to different assumptions for the wet scavenging of aerosols entrained above convective cloud bases in a global model motivates the ongoing need to better understand and model the activation and impaction processes that aerosols undergo after entrainment into convective updrafts.
Abstract. In this paper, we investigate the coagulation of interstitial aerosol particles (particles too small to activate to cloud droplets) with cloud drops, a process often ignored in aerosol-climate models. We use the GEOS-Chem-TOMAS (Goddard Earth Observing System-Chemistry TwO-Moment Aerosol Sectional) global chemical transport model with aerosol microphysics to calculate the changes in the aerosol size distribution, cloud-albedo aerosol indirect effect, and direct aerosol effect due to the interstitial coagulation process. We find that inclusion of interstitial coagulation in clouds lowers total particle number concentrations by 15–21% globally, where the range is due to varying assumptions regarding activation diameter, cloud droplet size, and ice cloud physics. The interstitial coagulation process lowers the concentration of particles with dry diameters larger than 80 nm (a proxy for larger CCN) by 10–12%. These 80 nm particles are not directly removed by the interstitial coagulation but are reduced in concentration because fewer smaller particles grow to diameters larger than 80 nm. The global aerosol indirect effect of adding interstitial coagulation varies from +0.4 to +1.3 W m−2 where again the range depends on our cloud assumptions. Thus, the aerosol indirect effect of this process is significant, but the magnitude depends greatly on assumptions regarding activation diameter, cloud droplet size, and ice cloud physics. The aerosol direct effect of the interstitial coagulation process is minor (< 0.01 W m−2) due to the shift in the aerosol size distribution at sizes where scattering is most effective being small. We recommend that this interstitial scavenging process be considered in aerosol models when the size distribution and aerosol indirect effects are important.
Abstract. Wet deposition processes are highly efficient in the removal of aerosols from the atmosphere, and thus strongly influence global aerosol concentrations, and clouds, and their respective radiative forcings. In this study, physically detailed size-dependent below-cloud scavenging parameterizations for rain and snow are implemented in the ECHAM5-HAM global aerosol-climate model. Previously, below-cloud scavenging by rain in the ECHAM5-HAM was simply a function of the aerosol mode, and then scaled by the rainfall rate. The below-cloud scavenging by snow was a function of the snowfall rate alone. The global mean aerosol optical depth, and sea salt burden are sensitive to the below-cloud scavenging coefficients, with reductions near to 15% when the more vigorous size-dependent below-cloud scavenging by rain and snow is implemented. The inclusion of a prognostic rain scheme significantly reduces the fractional importance of below-cloud scavenging since there is higher evaporation in the lower troposphere, increasing the global mean sea salt burden by almost 15%. Thermophoretic effects are shown to produce increases in the global and annual mean number removal of Aitken size particles of near to 10%, but very small increases (near 1%) in the global mean below-cloud mass scavenging of carbonaceous and sulfate aerosols. Changes in the assumptions about the below-cloud scavenging by rain of particles with radius smaller than 10 nm do not cause any significant changes to the global and annual mean aerosol mass or number burdens, despite a change in the below-cloud number removal rate for nucleation mode particles by near to five-fold. Annual and zonal mean nucleation mode number concentrations are enhanced by up to 30% in the lower troposphere with the more vigourous size-dependent below-cloud scavenging. Closer agreement with different observations is found when the more physically detailed below-cloud scavenging parameterization is employed in the ECHAM5-HAM model.
