ALBERT

Description: We present a sensitivity study on the effects of spatial and temporal resolution of atmospheric relative humidity (RH) on calculated aerosol optical thickness (AOT) and the aerosol direct radiative effects (DRE) in a global model. Using the same aerosol fields simulated in the Global Modeling Initiative (GMI) model, we find that, on a global average, the calculated AOT from RH in 1° latitude by 1.25° longitude spatial resolution is 11% higher than that in 2° by 2.5° resolution, and the corresponding DRE at the top of the atmosphere is 8–9% higher for total aerosols and 15% higher for only anthropogenic aerosols in the finer spatial resolution case. The difference is largest over surface escarpment regions (e.g. 〉200% over the Andes Mountains) where RH varies substantially with surface terrain. The largest zonal mean AOT difference occurs at 50–60°N (16–21%), where AOT is also relatively larger. A similar increase is also found when the time resolution of RH is increased. This increase of AOT and DRE with the increase of model resolution is due to the highly non-linear relationship between RH and the aerosol mass extinction efficiency (MEE) at high RH (〉80%). Our study suggests that caution should be taken in a multi-model comparison (e.g. AeroCom) since the comparison usually deals with results coming from different spatial/temporal resolutions.

Description: We examine the response of Arctic gas and aerosol concentrations to perturbations in pollutant emissions from Europe, East and South Asia, and North America using results from a coordinated model intercomparison. These sensitivities to regional emissions (mixing ratio change per unit emission) vary widely across models and species. Intermodel differences are systematic, however, so that the relative importance of different regions is robust. North America contributes the most to Arctic ozone pollution. For aerosols and CO, European emissions dominate at the Arctic surface but Asian emissions become progressively more important with altitude, and are dominant in the upper troposphere. Sensitivities show strong seasonality: surface sensitivities typically maximize during boreal winter for European and during spring for East Asian and North American emissions. Mid-tropospheric sensitivities, however, nearly always maximize during spring or summer for all regions. Deposition of black carbon (BC) onto Greenland is most sensitive to North American emissions. North America and Europe each contribute ~40% of total BC deposition to Greenland, with ~20% from East Asia. Elsewhere in the Arctic, both sensitivity and total BC deposition are dominated by European emissions. Model diversity for aerosols is especially large, resulting primarily from differences in aerosol physics and removal. Comparison of aerosols with observations indicates problems in either the models or interpretation of the measurements. For gas phase pollutants such as CO and O3, which are relatively well-simulated, the processes contributing most to uncertainties depend on the source region. Uncertainties in the Arctic surface CO response to emissions perturbations are dominated by emissions for East Asian sources, while uncertainties in transport, emissions, and oxidation are comparable for European and North American sources. At higher levels, model-to-model variations in transport and oxidation are most important. Differences in photochemistry appear to play the largest role in the intermodel variations in Arctic ozone sensitivity.

Description: Knowledge of the global distribution of tropospheric aerosols is important for studying the effects of aerosols on global climate. Chemical transport models rely on archived meteorological fields, accounting for aerosol sources, transport and removal processes can simulate the global distribution of atmospheric aerosols. However, the accuracy of global aerosol modeling is limited. Uncertainty in location and strength of aerosol emission sources is a major factor in limiting modeling accuracy. This paper describes an effort to develop an algorithm for retrieving global sources of aerosol from satellite observations by inverting the GOCART aerosol transport model. To optimize inversion algorithm performance, the inversion was formulated as a generalized multi-term least-squares-type fitting. This concept uses the principles of statistical optimization and unites diverse retrieval techniques into a single flexible inversion procedure. It is particularly useful for choosing and refining a priori constraints in the retrieval algorithm. For example, it is demonstrated that a priori limitations on the partial derivatives of retrieved characteristics, which are widely used in atmospheric remote sensing, can also be useful in inverse modeling for constraining time and space variability of the retrieved global aerosol emissions. The similarities and differences with the standard "Kalman filter" inverse modeling approach and the "Phillips-Tikhonov-Twomey" constrained inversion widely used in remote sensing are discussed. In order to retain the originally high space and time resolution of the global model in the inversion of a long record of observations, the algorithm was expressed using adjoint operators in a form convenient for practical development of the inversion from codes implementing forward model simulations. The inversion algorithm was implemented using the GOCART aerosol transport model. The numerical tests we conducted showed successful retrievals of global aerosol emissions with a 2°×2.5° resolution by inverting the GOCART output. For achieving satisfactory retrieval from satellite sensors such as MODIS, the emissions were assumed constant within the 24 h diurnal cycle and aerosol differences in chemical composition were neglected. Such additional assumptions were needed to constrain the inversion due to limitations of satellite temporal coverage and sensitivity to aerosol parameters. As a result, the algorithm was defined for the retrieval of emission sources of fine and coarse mode aerosols from the MODIS fine and coarse mode aerosol optical thickness data respectively. Numerical tests showed that such assumptions are justifiable, taking into account the accuracy of the model and observations and that it provides valuable retrievals of the location and the strength of the aerosol emissions. The algorithm was applied to MODIS observations during two weeks in August 2000. The global placement of fine mode aerosol sources retrieved from inverting MODIS observations was coherent with available independent knowledge. This was particularly encouraging since the inverse method did not use any a priori information about the sources and it was initialized under a "zero aerosol emission" assumption. The retrieval reproduced the instantaneous global MODIS observations with a standard deviation in fitting of aerosol optical thickness of ~0.04. The optical thickness during high aerosol loading events was reproduced with a standard deviation of ~48%. Applications of the algorithm for the retrieval of coarse mode aerosol emissions were less successful, mainly due to the currently existing lack of MODIS data over high reflectance desert dust sources. Possibilities for enhancing the global satellite data inversion by using diverse a priori constraints on the retrieval are demonstrated. The potential and limitations of applying our approach for the retrieval of global aerosol sources from aerosol remote sensing are discussed.

Description: The effects of unified aerosol sources on global aerosol fields simulated by different models are examined in this paper. We compare results from two AeroCom experiments, one with different (ExpA) and one with unified emissions, injection heights, and particle sizes at the source (ExpB). Surprisingly, harmonization of aerosol sources has only a small impact on the simulated diversity for aerosol burden, and consequently optical properties, as the results are largely controlled by model-specific transport, removal, chemistry (leading to the formation of secondary aerosols) and parameterizations of aerosol microphysics (e.g. the split between deposition pathways) and to a lesser extent on the spatial and temporal distributions of the (precursor) emissions. The burdens of black carbon and especially sea salt become more coherent in ExpB only, because the large ExpA diversity for these two species was caused by few outliers. The experiment also indicated that despite prescribing emission fluxes and size distributions, ambiguities in the implementation in individual models can lead to substantial differences. These results indicate the need for a better understanding of aerosol life cycles at process level (including spatial dispersal and interaction with meteorological parameters) in order to obtain more reliable results from global aerosol simulations. This is particularly important as such model results are used to assess the consequences of specific air pollution abatement strategies.

Description: The AeroCom exercise diagnoses multi-component aerosol modules in global modeling. In an initial assessment simulated global distributions for mass and mid-visible aerosol optical thickness (aot) were compared among 20 different modules. Model diversity was also explored in the context of previous comparisons. For the component combined aot general agreement has improved for the annual global mean. At 0.11 to 0.14, simulated aot values are at the lower end of global averages suggested by remote sensing from ground (AERONET ca. 0.135) and space (satellite composite ca. 0.15). More detailed comparisons, however, reveal that larger differences in regional distribution and significant differences in compositional mixture remain. Of particular concern are large model diversities for contributions by dust and carbonaceous aerosol, because they lead to significant uncertainty in aerosol absorption (aab). Since aot and aab, both, influence the aerosol impact on the radiative energy-balance, the aerosol (direct) forcing uncertainty in modeling is larger than differences in aot might suggest. New diagnostic approaches are proposed to trace model differences in terms of aerosol processing and transport: These include the prescription of common input (e.g. amount, size and injection of aerosol component emissions) and the use of observational capabilities from ground (e.g. measurements networks) or space (e.g. correlations between aerosol and clouds).

Description: Simulation results of global aerosol models have been assembled in the framework of the AeroCom intercomparison exercise. In this paper, we analyze the life cycles of dust, sea salt, sulfate, black carbon and particulate organic matter as simulated by sixteen global aerosol models. The differences among the results (model diversities) for sources and sinks, burdens, particle sizes, water uptakes, and spatial dispersals have been established. These diversities have large consequences for the calculated radiative forcing and the aerosol concentrations at the surface. Processes and parameters are identified which deserve further research. The AeroCom all-models-average emissions are dominated by the mass of sea salt (SS), followed by dust (DU), sulfate (SO4), particulate organic matter (POM), and finally black carbon (BC). Interactive parameterizations of the emissions and contrasting particles sizes of SS and DU lead generally to higher diversities of these species, and for total aerosol. The lower diversity of the emissions of the fine aerosols, BC, POM, and SO4, is due to the use of similar emission inventories, and does therefore not necessarily indicate a better understanding of their sources. The diversity of SO4-sources is mainly caused by the disagreement on depositional loss of precursor gases and on chemical production. The diversities of the emissions are passed on to the burdens, but the latter are also strongly affected by the model-specific treatments of transport and aerosol processes. The burdens of dry masses decrease from largest to smallest: DU, SS, SO4, POM, and BC. The all-models-average residence time is shortest for SS with about half a day, followed by SO4 and DU with four days, and POM and BC with six and seven days, respectively. The wet deposition rate is controlled by the solubility and increases from DU, BC, POM to SO4 and SS. It is the dominant sink for SO4, BC, and POM, and contributes about one third to the total removal of SS and DU species. For SS and DU we find high diversities for the removal rate coefficients and deposition pathways. Models do neither agree on the split between wet and dry deposition, nor on that between sedimentation and other dry deposition processes. We diagnose an extremely high diversity for the uptake of ambient water vapor that influences the particle size and thus the sink rate coefficients. Furthermore, we find little agreement among the model results for the partitioning of wet removal into scavenging by convective and stratiform rain. Large differences exist for aerosol dispersal both in the vertical and in the horizontal direction. In some models, a minimum of total aerosol concentration is simulated at the surface. Aerosol dispersal is most pronounced for SO4 and BC and lowest for SS. Diversities are higher for meridional than for vertical dispersal, they are similar for the individual species and highest for SS and DU. For these two components we do not find a correlation between vertical and meridional aerosol dispersal. In addition the degree of dispersals of SS and DU is not related to their residence times. SO4, BC, and POM, however, show increased meridional dispersal in models with larger vertical dispersal, and dispersal is larger for longer simulated residence times.

Description: Aerosols affect the Earth's energy budget directly by scattering and absorbing radiation and indirectly by acting as cloud condensation nuclei and, thereby, affecting cloud properties. However, large uncertainties exist in current estimates of aerosol forcing because of incomplete knowledge concerning the distribution and the physical and chemical properties of aerosols as well as aerosol-cloud interactions. In recent years, a great deal of effort has gone into improving measurements and datasets. It is thus feasible to shift the estimates of aerosol forcing from largely model-based to increasingly measurement-based. Our goal is to assess current observational capabilities and identify uncertainties in the aerosol direct forcing through comparisons of different methods with independent sources of uncertainties. Here we assess the aerosol optical depth (τ), direct radiative effect (DRE) by natural and anthropogenic aerosols, and direct climate forcing (DCF) by anthropogenic aerosols, focusing on satellite and ground-based measurements supplemented by global chemical transport model (CTM) simulations. The multi-spectral MODIS measures global distributions of aerosol optical depth (τ) on a daily scale, with a high accuracy of ±0.03±0.05τ over ocean. The annual average τ is about 0.14 over global ocean, of which about 21%±7% is contributed by human activities, as estimated by MODIS fine-mode fraction. The multi-angle MISR derives an annual average AOD of 0.23 over global land with an uncertainty of ~20% or ±0.05. These high-accuracy aerosol products and broadband flux measurements from CERES make it feasible to obtain observational constraints for the aerosol direct effect, especially over global the ocean. A number of measurement-based approaches estimate the clear-sky DRE (on solar radiation) at the top-of-atmosphere (TOA) to be about -5.5±0.2 Wm-2 (median ± standard error from various methods) over the global ocean. Accounting for thin cirrus contamination of the satellite derived aerosol field will reduce the TOA DRE to -5.0 Wm-2. Because of a lack of measurements of aerosol absorption and difficulty in characterizing land surface reflection, estimates of DRE over land and at the ocean surface are currently realized through a combination of satellite retrievals, surface measurements, and model simulations, and are less constrained. Over the oceans the surface DRE is estimated to be -8.8±0.7 Wm-2. Over land, an integration of satellite retrievals and model simulations derives a DRE of -4.9±0.7 Wm-2 and -11.8±1.9 Wm-2 at the TOA and surface, respectively. CTM simulations derive a wide range of DRE estimates that on average are smaller than the measurement-based DRE by about 30-40%, even after accounting for thin cirrus and cloud contamination. A number of issues remain. Current estimates of the aerosol direct effect over land are poorly constrained. Uncertainties of DRE estimates are also larger on regional scales than on a global scale and large discrepancies exist between different approaches. The characterization of aerosol absorption and vertical distribution remains challenging. The aerosol direct effect in the thermal infrared range and in cloudy conditions remains relatively unexplored and quite uncertain, because of a lack of global systematic aerosol vertical profile measurements. A coordinated research strategy needs to be developed for integration and assimilation of satellite measurements into models to constrain model simulations. Enhanced measurement capabilities in the next few years and high-level scientific cooperation will further advance our knowledge.

Description: We evaluate black carbon (BC) model predictions from the AeroCom model intercomparison project by considering the diversity among year 2000 model simulations and comparing model predictions with available measurements. These model-measurement intercomparisons include BC surface and aircraft concentrations, aerosol absorption optical depth (AAOD) retrievals from AERONET and Ozone Monitoring Instrument (OMI) and BC column estimations based on AERONET. In regions other than Asia, most models are biased high compared to surface concentration measurements. However compared with (column) AAOD or BC burden retreivals, the models are generally biased low. The average ratio of model to retrieved AAOD is less than 0.7 in South American and 0.6 in African biomass burning regions; both of these regions lack surface concentration measurements. In Asia the average model to observed ratio is 0.7 for AAOD and 0.5 for BC surface concentrations. Compared with aircraft measurements over the Americas at latitudes between 0 and 50N, the average model is a factor of 8 larger than observed, and most models exceed the measured BC standard deviation in the mid to upper troposphere. At higher latitudes the average model to aircraft BC ratio is 0.4 and models underestimate the observed BC loading in the lower and middle troposphere associated with springtime Arctic haze. Low model bias for AAOD but overestimation of surface and upper atmospheric BC concentrations at lower latitudes suggests that most models are underestimating BC absorption and should improve estimates for refractive index, particle size, and optical effects of BC coating. Retrieval uncertainties and/or differences with model diagnostic treatment may also contribute to the model-measurement disparity. Largest AeroCom model diversity occurred in northern Eurasia and the remote Arctic, regions influenced by anthropogenic sources. Changing emissions, aging, removal, or optical properties within a single model generated a smaller change in model predictions than the range represented by the full set of AeroCom models. Upper tropospheric concentrations of BC mass from the aircraft measurements are suggested to provide a unique new benchmark to test scavenging and vertical dispersion of BC in global models.

Description: We examine the response of Arctic gas and aerosol concentrations to perturbations in pollutant emissions from Europe, East and South Asia, and North America using results from a coordinated model intercomparison. These sensitivities to regional emissions (mixing ratio change per unit emission) vary widely across models and species. Intermodel differences are systematic, however, so that the relative importance of different regions is robust. North America contributes the most to Arctic ozone pollution. For aerosols and CO, European emissions dominate at the Arctic surface but East Asian emissions become progressively more important with altitude, and are dominant in the upper troposphere. Sensitivities show strong seasonality: surface sensitivities typically maximize during boreal winter for European and during spring for East Asian and North American emissions. Mid-tropospheric sensitivities, however, nearly always maximize during spring or summer for all regions. Deposition of black carbon (BC) onto Greenland is most sensitive to North American emissions. North America and Europe each contribute ~40% of total BC deposition to Greenland, with ~20% from East Asia. Elsewhere in the Arctic, both sensitivity and total BC deposition are dominated by European emissions. Model diversity for aerosols is especially large, resulting primarily from differences in aerosol physical and chemical processing (including removal). Comparison of modeled aerosol concentrations with observations indicates problems in the models, and perhaps, interpretation of the measurements. For gas phase pollutants such as CO and O3, which are relatively well-simulated, the processes contributing most to uncertainties depend on the source region and altitude examined. Uncertainties in the Arctic surface CO response to emissions perturbations are dominated by emissions for East Asian sources, while uncertainties in transport, emissions, and oxidation are comparable for European and North American sources. At higher levels, model-to-model variations in transport and oxidation are most important. Differences in photochemistry appear to play the largest role in the intermodel variations in Arctic ozone sensitivity, though transport also contributes substantially in the mid-troposphere.

Description: The effects of unified aerosol sources on global aerosol fields simulated by different models are examined in this paper. We compare results from two AeroCom experiments, one with different (ExpA) and one with unified emissions, injection heights, and particle sizes at the source (ExpB). Surprisingly, harmonization of aerosol sources has only a small impact on the simulated inter-model diversity of the global aerosol burden, and consequently global optical properties, as the results are largely controlled by model-specific transport, removal, chemistry (leading to the formation of secondary aerosols) and parameterizations of aerosol microphysics (e.g., the split between deposition pathways) and to a lesser extent by the spatial and temporal distributions of the (precursor) emissions. The burdens of black carbon and especially sea salt become more coherent in ExpB only, because the large ExpA diversities for these two species were caused by a few outliers. The experiment also showed that despite prescribing emission fluxes and size distributions, ambiguities in the implementation in individual models can lead to substantial differences. These results indicate the need for a better understanding of aerosol life cycles at process level (including spatial dispersal and interaction with meteorological parameters) in order to obtain more reliable results from global aerosol simulations. This is particularly important as such model results are used to assess the consequences of specific air pollution abatement strategies.