Aerosol climate effects and air quality impacts from 1980 to 2030

We investigate aerosol effects on climate for 1980, 1995 (meant to reflect present day) and 2030 using the NASA Goddard Institute for Space Studies climate model coupled to an on-line aerosol source and transport model with interactive oxidant and aerosol chemistry. Aerosols simulated include sulfates, organic matter (OM), black carbon (BC), sea-salt and dust and, additionally, the amount of tropospheric ozone is calculated, allowing us to estimate both changes to air quality and climate for different time periods and emission amounts. We include both the direct aerosol effect and indirect aerosol effects for liquid-phase clouds. Future changes for the 2030 A1B scenario are examined, focusing on the Arctic and Asia, since changes are pronounced in these regions. Our results for the different time periods include both emission changes and physical climate changes. We find that the aerosol indirect effect (AIE) has a large impact on photochemical processing, decreasing ozone amount and ozone forcing, especially for the future (2030–1995). Ozone forcings increase from 0 to 0.12 W m−2 and the total aerosol forcing decreases from −0.10 to −0.94 W m−2 (AIE decreases from −0.13 to −0.68 W m−2) for 1995–1980 versus 2030–1995. Over the Arctic we find that compared to ozone and the direct aerosol effect, the AIE contributes the most to net radiative flux changes. The AIE, calculated for 1995–1980, is positive (1.0 W m−2), but the magnitude decreases (−0.3 W m−2) considerably for the future scenario. Over Asia, we evaluate the role of biofuel- and transportation-based emissions (for BC and OM) via a scenario (2030A) that includes a projected increase (factor of 2) in biofuel- and transport-based emissions for 2030 A1B over Asia. Projected changes from present day due to the 2030A emissions versus 2030 A1B are a factor of 4 decrease in summertime precipitation in Asia. Our results are sensitive to emissions used. Uncertainty in present-day emissions suggests that future climate projections warrant particular scrutiny.


Introduction
To address the impacts of future climate, projections in the concentrations of both greenhouse gases (GHGs) and non-GHGs and their potential climate effects are important. The counteracting influence to GHG warming from the direct and indirect aerosol effects remains uncertain due in part to the complexity in treating aerosol-cloud interactions (Intergovernmental Panel on Climate Change 2007). Apart from the climate impacts of aerosols, changes to cloud and radiation fields from aerosols also affect the formation of ozone, which depends on the photo-oxidation of precursor emissions. Thus short-lived species such as aerosols and ozone are coupled, are both important for their climate and air quality impacts and have the potential to be important in near-term mitigation efforts.
The radiative forcing associated with ozone is ∼0.4 W m −2 and for aerosols, forcings range from 0.3 (for black carbon (BC)) to −2.2 W m −2 (for reflecting aerosols such as sulfates, organic matter (OM), nitrates and aerosol-cloud changes) based on emission changes from 1750 to present day (PD) (Intergovernmental Panel on Climate Change 2007). Absorbing aerosols are thought to inhibit surface ozone formation by reducing the photolysis rates in polluted areas; whereas reflecting aerosols are thought to increase ozone formation (He and Carmichael 1999). Thus, future forcings depend on emission scenarios and relative changes in the distribution of precursors and reflecting versus absorbing aerosols (BC).
Here, we examine ozone and aerosol-cloud climate effects for three time periods: 1980, 1995 and 2030 using the NASA Goddard Institute for Space Studies (GISS) climate model (ModelE) coupled to an on-line aerosol chemistry/transport model. Aerosols simulated include sulfates, OM, BC, seasalt and dust. The amount of tropospheric ozone is also calculated, thus allowing us to estimate both changes to air quality and climate for different time periods and emission amounts. Changes to future ozone forcing from coupling with an interactive sulfate chemistry model have been included in prior studies (Unger et al 2006b(Unger et al , 2006a, however the addition of carbonaceous aerosols and their cloud impacts on future ozone forcing have not been included. Also, most prior studies of the aerosol indirect effects (AIE) focus on PD to preindustrial (PI) time eras, and very few studies (e.g. Kloster et al (2008)) account for changes in the AIE for future scenarios in simulations that include full couplings between the ozone and sulfate chemistry.
Since BC aerosols are more absorbing in the visible and due to their climate warming potential (Menon and Del Genio 2007), we single out their climate effects. We thus examine an alternative scenario with twice the 2030 A1B emissions from biofuel and transportation sources for OM and BC over Asia, due to the presence of a large proportion of these emissions over Asia (Fernandes et al 2007) and a factor of two uncertainty in current residential emissions. Over India a large proportion of BC is from biofuel (Venkataraman et al 2005); and the transportation sector over Asia is expected to grow with rapid economic expansion (Streets 2005). Furthermore, the transport of Asian emissions to Arctic snow surfaces (Koch and Hansen 2005); and implied association between mid-latitude pollution and longwave cloud emissivity in the Arctic Zhao 2006, Lubin andVogelmann 2006) suggest that these relevant issues should be examined for future changes. In section 2 we describe the methodology and model simulations; section 3 includes the results from the simulations for both climate and air quality effects and the regional impacts of particulate matter (aerosols) and ozone for the Arctic and Asia. Finally in section 4 the discussion and conclusions are presented.

Methodology
We use the GISS GCM (ModelE)  (4 • × 5 • horizontal resolution and 20 vertical layers from 974 to 0.2 hPa) coupled to an on-line aerosol chemistry/transport model (Koch and Hansen 2005, Koch et al 2007a, Unger et al 2006b that includes schemes to treat the AIE (Menon et al 2008, Menon andDel Genio 2007) and thus feedbacks between aerosols and clouds. We simulate the sulfur cycle and include source terms for OM and BC aerosols that are allowed to be transported and subject to the same physical processes (wet/dry deposition) as sulfates, but with reduced solubility. Both sea-salt and dust are treated as natural emissions. We use emissions for 1995 meant to reflect PD simulations and emissions for 1980 are used mainly to examine impacts over time periods where emission trends differ considerably across the developed and developing world. For future predictions, we confine our analysis to the 2030 A1B scenario since it represents a mid-range estimate of future growth in fossil fuel emissions (other A1B aerosol projections may differ substantially ). In this scenario, the atmospheric column burden for sulfates increase by 40% and the carbonaceous aerosols decrease by 16%. For the 1980, PD and 2030 simulations we use SST and sea-ice climatologies for 1975-1984, 1990-1999 and 2030-2039, respectively. Thus, both emission changes and climate effects are included and we do not separate between them. Unger et al (2006a) found that climate and air pollution effects of emissions are more determined by emission changes for the 2030 A1B scenario than physical climate changes. Further details of the simulations and emissions used are in appendix A.
We examine three sets of simulations for 1980, PD and 2030, shown in table 1, listed as Sim 1980, Sim PD, Sim 2030, respectively, that include ozone-sulfate couplings, carbonaceous aerosols and the direct and indirect aerosol effects. Additional sensitivity tests include similar sets of simulations with only sulfates (OS) and without the AIE (listed as Sim 1980OS, Sim PDOS, and Sim 2030OS), and a simulation where biofuel and transportation sources of OM and BC are multiplied by 2 for Sim 2030 over Asia (10 • S-40 • N, 58 • -122 • E), described as Sim 2030A.

Results
Here, we present our analysis of changes in climate and air pollution due to the effects of ozone, aerosols and aerosolcloud interactions as well as associated physical climate changes over the different time intervals we examine.

Aerosol and ozone distribution
To evaluate climate change for 1980-2030, we analyze the aerosol and ozone distributions for differences between PD and 1980 (referred to as 'past') and for differences between  2030A1B and PD (referred to as 'future'). Figure 1 shows annual surface concentrations for ozone (ppbv) and sulfate (µg m −3 ). Past ozone concentrations increased overall over most parts of Asia and over Africa but decreased slightly over the US and Europe (similar to observed changes) in response to precursor emissions reductions legislation in these regions. Similar changes were found for sulfates described in more detail in Unger et al (2006b) that performed similar analysis but without carbonaceous aerosols and the AIE. Figures 2 and 3 show OM and BC surface concentrations (µg m −3 ) for the past and future from both fossil, biofuel and biomass emissions. Similar to sulfate concentrations, the US and western Europe indicate a decrease in OM and BC emissions from fossil and biofuel sources for the past, whereas emissions from these sources increased over parts of S America, Africa, India and China. For biomass sources, increases (>2 µg m −3 , especially for OM) over Amazonia, southern parts of Africa, south-east Asia and Australia are noted for the past. For the past case the all-sky total (sulfate, BC and OM) aerosol optical depth (AOD) is 0.00. Both aerosol optical properties and mass distributions have been extensively evaluated with measurements, where available, and are described in Yu et al (2006) and Koch et al (2007b) in more detail. Additionally, AOD for Sim PD has been evaluated with satellite data for the Atlantic Ocean region (Menon et al 2008).
In general, the model tends to underestimate AOD, probably due to the assumed particle sizes and the underestimation of natural aerosols in southern tropical oceans (Koch et al 2007b).
For the future case, mainly sulfates-, ozone-and biomassbased carbonaceous aerosols are projected to increase over Asia, while the Middle East shows an increase in fossil and biofuel BC emissions (due to increased oil and gas extraction under the A1B scenario). The AOD for the future case is 0.02 (17% increase). While most emissions for sulfates, OM and BC decrease over the US and Europe, small ozone increases of the order of 2 ppbv are noted everywhere, which have been shown in a previous study to be largely due to future increases in methane since other short-lived precursor emissions decrease across the US and Europe for this scenario (Unger et al 2006b). Differences between the ozone and sulfate concentrations for Sim 2030 versus Sim 2030OS are shown in figure 4. Globally, although the decreases for both ozone and  sulfate are small, their distributions differ substantially. Ozone concentrations increase by a small amount over Europe and parts of China, eastern Australia and the Amazon, whereas, sulfate concentrations increase over tropical ocean regions and SE Asia by over 10%. Asia, Europe and eastern US all show an increase of 5-8%. Physical reasons for the differences are further discussed in the next sections where we consider cloud and radiation changes due to aerosols that could affect ozone formation and distribution, and sulfate precursors.

Aerosol and ozone forcings
For the aerosol and ozone distributions described above, we calculate the annual average instantaneous radiative forcings at the tropopause for aerosols and ozone for the past and future that are given in table 2. Note that the differences between forcings for the different time periods also include changes from SSTs/sea-ice and GHGs in addition to differences in aerosol emissions. Given almost neutral changes in carbonaceous aerosols for the past case and the slight decrease Table 2. Annual average instantaneous shortwave radiative forcings (W m −2 ) at the tropopause for the listed aerosol species for the three simulations. Values for ozone and the aerosol indirect effect include longwave forcings (W m −2 ). OM refers to organic matter, BCB to black carbon from biomass sources and BCF to black carbon from fossil and biofuel sources.

Species
Sim 1980   in sulfate levels, the direct effect for the past case is close to 0.0 W m −2 , and the larger increase in sulfate concentrations for the future results in higher negative forcings (−0.26 W m −2 ) since carbonaceous aerosols are almost constant. Table 3 indicates the forcing efficiency (radiative forcing/column burden) for the different aerosols we simulate that may be used to estimate the relative role that mitigation efforts may play in modulating aerosol distributions. BC from fossil and biofuel sources clearly dominate forcing efficiency values and are almost a factor of 4 higher than sulfates and a factor of 10 higher than OM. Thus, for similar emission reductions, climate impacts from changes to BC distributions may play a larger role in mitigation studies than would impacts from changes to sulfates. For ozone, we obtain forcings of 0.00 and 0.12 W m −2 for the past and future cases. For the future case, forcings of Table 3. Radiative forcing efficiencies (W g −1 ) calculated from the ratio of direct forcing to column burdens for the listed aerosol species. OM refers to organic matter, BCB to black carbon from biomass sources and BCF to black carbon from fossil and biofuel sources.

Species
Sim 1980  were found to be larger for the future case since the sulfate increase and the indirect effect is much more pronounced in the future scenario. The AIE may influence ozone chemistry through a number of pathways that change photolysis rates and wet deposition of precursors and thus ozone forcing. In general, for all of the three time periods analyzed, key photolysis rates increase by around 20% over the oceans and decrease by a few per cent over land when carbonaceous aerosols and the AIE are included in the model versus the OS simulations (not shown). The increases in J(O 1 D) (photolysis of ozone to form O 1 D) drive the reductions in surface ozone (figure 4) resulting in enhanced OH levels. OH increases by at least 5% across most mid-latitude regions and up to 20-30% at higher latitudes. The higher OH levels drive increases in gas-phase production of sulfate that result in higher surface concentrations across the northeastern US and Europe (figure 4). The sulfate increases in tropical regions shown in figure 4 are caused by increased aqueous-phase production for the simulation with the AIE and carbonaceous aerosols versus OS.

Aerosol-cloud interactions
The values for the AIE (calculated from changes to net cloud radiative forcing) for the past and future cases are −0.13 and −0.68 W m −2 , respectively, as shown in table 2. Figure 5 shows the spatial distribution of cloud droplet number (CDNC), cloud optical depth (COD) and the AIE for the past and future cases. For past changes, the indirect effect is negative over land areas, but is positive over the South Asian subcontinent and SE Asia, parts of Africa and the northern US. These changes should be related to CDNC such that an increase in CDNC (and thereby a reduction in cloud droplet effective radii (R eff )) should generally coincide with areas of negative indirect effects, unless cloud liquid water and meteorology change. This is observed for the northern US and African regions. The increase in CDNC for the SE Asian region is due to the increase in carbonaceous aerosols, especially from biomass sources. However, the indirect effect was positive and the change in the sign of the AIE may be related to the decrease in cloud cover (CC) (not shown). Although CC (both total and low) increases globally for the past case (0.58%), these changes are not a result of the indirect effect as similar increases were also observed for simulations OS and may be related to climate feedbacks from imposed SSTs/sea-ice and GHG changes. For the future case, similar to the increases observed for COD, the liquid water path (LWP) (not shown) indicates an increase of 2.36 g m −2 . Thus, the enhanced AIE for the future case is mainly due to the increase in LWP (at the expense of suppressed precipitation), and COD rather than CC. (Although physical climate changes may cause some of the increase in LWP or COD, compared to Sim (2030-PD), the increase in COD and LWP for Sim (2030OS-PDOS) are smaller by a factor of 1.6 and 20%, respectively, thus indicating an indirect effect signal). In terms of the distribution of the AIE impacts, for future changes the average AIE over the NH land areas is 0.61 W m −2 compared to −1.49 W m −2 found for the SH land areas, indicative of the effects of decreased aerosol emissions (especially sulfates) over most mid-latitude land areas. AIE values over land and ocean are 0.03 and −0.89 W m −2 , respectively, indicating where the response is likely to shift with a change in emissions.

Impacts of the aerosol indirect effect on ozone forcing
The impact of the AIE and carbonaceous aerosols on the annual average (J(O 1 D)) and (J(NO 2 )) photolysis rates at the surface for the future case given as Sim (2030-PD)-Sim (2030OS-PDOS) is shown in figure 6. Photochemistry is suppressed in subtropical and tropical regions where COD has increased due to the increased aerosol loading between 2030 and 1995 (figure 5) resulting in a dampening of the global ozone increase for the future case. The inclusion of the AIE may also impact other processes that affect ozone, especially wet deposition.  other wet deposited species. It is evident that the spatial pattern of the changes follows the changes in CDNC and LWP/COD for the time period (e.g. figure 5). Where precipitation decreases in subtropical regions due to increases in aerosol loading, the wet deposition of ozone precursors is reduced relative to the OS simulations. Conversely, in some NH midlatitude regions where aerosols have been reduced, especially the northeastern US, the opposite effect (not desirable for reducing acid deposition) is apparent. A decreased in wet deposition of ozone precursors would tend to increase ozone amounts in the troposphere. However, the decrease in surface ozone found in figure 4 suggest that the influence of aerosolcloud interactions on wet deposition must be outweighed by the influence on photolysis rates.

Regional impacts
Due to the heterogeneity in aerosol emissions, we evaluate regional changes over two areas of interest: the Arctic (60 • -90 • N) and Asia (10 • S-40 • N, 58 • -122 • E). These regions are chosen to analyze the influence of short-lived pollutants (aerosols and ozone) on the Arctic surface radiation (given the sensitivity of high-latitude regions to future climate change) and changes to precipitation and air quality from aerosols and ozone over the Asian monsoon region.

Aerosol and ozone impacts on the Arctic climate.
To understand the relative influence of aerosol emissions on the distant Arctic climate we examine the annual values of radiative fluxes, cloud cover and temperature over the Arctic region and include seasonality since Arctic summer versus winter clouds have different distributions of liquid versus solidphase clouds. Note that our simulations include physical climate changes in addition to aerosol changes and thus radiative forcing differences are actually flux differences. Our simulations do not include impacts of aerosols on cold cloud processes, via changes to ice nuclei. Thus changes we identify are a result of feedbacks associated with aerosol direct radiative effects and aerosol impacts on liquid-phase clouds, whenever they exist. Table 4 lists the annual, summer and winter averages for several variables (radiation and cloud fields) for the three time periods. We find an increase in surface radiation, TOA net radiation and a decrease in low and total CC for the future case compared to the past. The indirect effect decreases (more negative) for the future case. As expected during winter LWCRF dominates the net cloud radiative forcing giving us a net positive indirect effect though the magnitude gets weaker for the future. Although the low and total cloud cover decrease, liquid-phase cloud cover increases for the future case (at the expense of a decrease in the ice-phase clouds) and surface air temperature (SAT) decreases. Projected changes in ozone forcing for the Arctic are small and decreases for the future (goes from 0.03 to −0.02 W m −2 ). The aerosol direct forcing values in table 4 are also rather small. Given such a large change in the AIE, from 1.0 to −0.30 W m −2 (from the past to the future), we compare some of the simulated cloud microphysical values for Sim PD with available measurements over the Arctic (Barrow region), to constrain the contribution of the indirect effect to Arctic surface radiative changes, in appendix B.
As shown in appendix B, while R eff and COD are within retrieval uncertainties as given in Zhao and Garrett (2008), larger/smaller values are simulated for CDNC and LWP that are not within retrieval uncertainties. However, radiative fluxes are more sensitive to both R eff and COD and thus due to compensating biases in CDNC and LWP, simulated radiative fluxes may be a useful indicator of future changes. Since downwelling longwave fluxes (DLW) are expected to contribute most to the Arctic surface radiation balance (Francis andHunter 2007, Garrett andZhao 2006), we examine DLW fluxes obtained from the television and infrared observation satellite operational vertical sounder (TOVS) instrument (Francis 1997) and compare them to simulated DLW fluxes. These DLW fluxes are sensitive to cloud fraction, cloud droplet size, and cloud water paths and thus provide a measure of the relative impacts of clouds on the surface Table 4. Annual, June-July-August (JJA) and December-January-February (DJF) average values for TOA net radiation (NR-TOA), surface net radiation (NR-Sfc), shortwave and longwave cloud radiative forcing (SWCRF, LWCRF), aerosol direct forcing (DF) and ozone forcing (OF) at the tropopause, surface air temperature (SAT), water cloud cover (WCC), low cloud cover (LCC), total cloud cover (TCC) for the Arctic region (60 • -90 • N) for the three simulations.  Shindell (2007) finds that SAT response may be more coupled to remote forcings than local forcings and can show an opposite response to local forcings. Thus, SAT response may be more dominated by non-aerosol forcings.

Climate impacts from increase in transportation-and biofuel-based carbonaceous emissions over Asia.
For the future scenario considered here, carbonaceous aerosols are projected to decline over most parts of the globe, including Asia, where concentrations are generally high. Koch et al (2007b) finds that PD carbonaceous aerosols over Asia may be lower than observed amounts which may in part be due to emission sources used (Cao et al 2006). Thus, it is quite likely that the sensitivity test (Sim 2030A) we carry out, where carbonaceous aerosols (OM and BC) from biofuel and transportation sources for 2030 A1B increase by 2× over Asia (10 • S-40 • N and 58 • -122 • E), may not be unrealistic.  Table 5. Annual average differences in ozone amount, aerosol surface mass, TOA net radiation (NR-TOA), surface net radiation (NR-Sfc), ozone forcing and aerosol forcing (aerosol direct +indirect effect) for the Asian region (10 • S-40 • N; 58 • to 122 • E) for the future cases.

Simulation
Ozone Aerosol mass NR-TOA NR-Sfc Ozone forcing Aerosol forcing Sim 2030-PD 8. We examine the differences in aerosol and ozone forcings for this region for the future case as shown in table 5. Globally, average values for the ozone and the total aerosol forcing (direct + indirect) are 0.12 and −0.94 (−0.26 + −0.68) W m −2 , respectively, while over the Asian region considered here, average values for ozone and the total aerosol forcing are 0.23 and −0.86 (−0.50 + −0.36) W m −2 , respectively. Thus, Asian averages are about twice the global averages for ozone and the direct effect and almost half the global value for the AIE. For Sim 2030A, ozone forcing increases by 9% and the total aerosol forcing decreases by 75% mainly due to the decrease in the indirect effect (decrease in low CC and LWP probably due to the heating effects of BC aerosols), as shown in table 5. The indirect effect thus increases from −0.36 to 0.27 W m −2 . With the increased carbonaceous emissions, apart from the large increase in submicron aerosol mass, net TOA radiation decreases by 30%, implying a less reflective atmosphere and a warmer surface, with an increase in surface ozone and particulate mass.
Precipitation changes for June-August are shown in figure 9 for the past and future cases. The increase in summer precipitation over southern China, similar to that found in Menon et al (2002), and large parts of SE Asia for the past case is reduced for the future cases. Similar changes are observed for the eastern and western parts of India. These changes in precipitation are due to both emission and physical climate changes. The impact of the increase in carbonaceous aerosol emissions on precipitation for the Asian region studied here are a factor of 4 decrease for Sim (2030A-PD) versus Sim (2030-PD), shown in figure 9. Thus, increases in biofuel-and transport-based carbonaceous emissions result in a reduction in average precipitation, with some regions (e.g. northern China and eastern India) indicating an increase in precipitation. Similar spatial redistributions in precipitation due to the heating effects of Asian haze were found by Chung and Ramanathan (2006) and Ramanathan et al (2005) for the South Asian region.
Since Koch and Hansen (2005) found that Asian emissions may impact the distant Arctic we evaluate similar diagnostics as shown in table 4 for Sim 2030A-PD, but for annual differences only. We find that values for the aerosol direct forcing (−0.03 W m −2 ) and ozone forcing (0.00 W m −2 ) are negligible, though are more positive than those obtained for Sim 2030-PD. The indirect effect (−0.51 W m −2 ) and values for net radiation both at the surface (1.39 W m −2 ) and TOA (0.37 W m −2 ) are larger than that found for Sim 2030-PD, mainly due to changes in LCC (−0.95) and WCC (−0.45). Thus, the bigger impact from increased carbonaceous emissions over Asia is in the aerosol indirect effect (mainly from cloud changes).

Discussion and conclusions
We use a series of model simulations, for 1980for , 1995 and 2030, to estimate changes in aerosol and ozone forcing, both globally and across two regions of interest-the Arctic and Asia. Simulated climate include effects from both changing emissions and a changing climate (from changes in SSTs/seaice and GHGs).
With a projected increase in sulfate emissions and a slight reduction in carbonaceous aerosols for 2030, the aerosol forcing (both the direct and indirect effects) decreases from −0.10 (for PD versus 1980) to −0.94 (for 2030 versus PD) W m −2 . Similar changes for the ozone forcing are an increase from 0 to 0.12 W m −2 due to an increase in ozone precursors. Differences between simulations with carbonaceous aerosols and the AIE versus simulations with just sulfate aerosol direct effects are a small decrease in surface sulfate and ozone with larger changes in surface sulfate over tropical regions. Globally, ozone forcing decreases between 2030 and 1995 relative to simulations with just sulfates. The ozone reductions appear to be driven by a reduction in photolysis due to increased cloud optical depth in the subtropics as a result of the enhanced aerosol loading there.
Over the Arctic, the AIE dominates the ozone and direct aerosol forcing, with positive values of the indirect effect of +1 W m −2 , calculated for PD-1980, changing to negative values of −0.3 W m −2 for the future case (2030-PD). Since the AIE was only applied to liquid-phase clouds, model simulated cloud microphysical properties and radiative fluxes for PD were compared to available measurements for the summer season when shallow liquid-phase stratus clouds persist in the Arctic. Simulated cloud droplet size and optical thickness are comparable to retrieved values, however, cloud droplet number and LWP were over and under-predicted. Since radiative fluxes are more sensitive to cloud optical depth and droplet size, simulated downwelling LW fluxes were comparable to retrieved values. Both total and low level clouds decrease in the future scenario, and thus the downwelling LW flux and net surface radiation increase. However, the surface temperature decreases, which may be subject more to large-scale dynamical changes than to the aerosol forcing.
Over Asia, both ozone and the aerosol direct effect are twice the global averages, while the AIE is half of the global average. With the expected increase in transportation over Asia and the larger proportion of global biofuel used in Asia, we carried out a sensitivity test (Sim 2030A) with twice the biofuel-and transport-based emissions from carbonaceous aerosols for Asia for the 2030 A1B scenario. Resulting changes are a 9% increase in ozone forcing, a 75% reduction in the aerosol forcing and a factor of four decrease in precipitation for 2030A-PD versus 2030-PD. Impacts on the distant Arctic due to the increased Asian emissions are mostly evident in the indirect aerosol effect (from cloud changes).
Our results are sensitive to the change in physical climate from GHG and SSTs/sea-ice changes and projected emissions of reflecting versus absorbing aerosols. Since prior work by Unger et al (2006a) indicated that climate effects are more sensitive to emission changes than physical climate changes for the future scenario considered here, uncertainties in aerosol emissions or precursors would necessarily result in uncertain climate projections. Given the large aerosol forcing values we obtain and the uncertainty in indirect effect simulations, the PD simulations used in this letter were compared to satellite retrievals and meteorological data for June to August 2002 over the Atlantic Ocean to constrain the magnitude of the indirect effect (Menon et al 2008). Results suggest that our PD simulations may slightly over-predict the indirect effect over the Atlantic Ocean region and thus our AIE estimate may be considered an upper limit for expected changes. Relatively few studies account for the impacts of the AIE on tropospheric ozone. Our analysis suggests large potential impacts on photochemical processing and thus the aerosolcloud interactions need to be included in simulations of past and future ozone changes. In future work we will explore the impacts of the AIE on ozone production and loss mechanisms in different regions for past and future changes. are from Streets et al (2006). Biomass burning for the past and future are based on projected changes (Streets et al 2004(Streets et al , 2006 applied to the GFED emissions v1, 1997-2001, also described in Koch et al (2007b) (from e.g. van der Werf et al (2004)) together with the carbonaceous aerosol emission factors from Andreae and Merlet (2001).
Aerosols are externally mixed and to convert aerosol mass to aerosol number we use assumed log-normal size distributions as specified in Menon and Rotstayn (2006). Although the model resolution is coarse the tracers are advected with a highly non-diffusive quadratic upstream scheme that increases the effective resolution of the tracer fields. The aerosol-cloud interactions we simulate for liquidphase stratus and cumulus clouds, include (a) the first AIE: increased (decreased) cloud reflectivity due to an increase (decrease) in aerosols and cloud droplet number concentrations (CDNC) and reduced (increased) droplet sizes (R eff ) for similar cloud liquid water content (LWC) (Twomey 1991); and (b) the second AIE: change in cloud cover (CC), cloud liquid water path (LWP) and precipitation due to smaller droplet sizes that inhibit precipitation processes, thereby increasing CC and LWP (and thus cloud optical depths (COD)) (Albrecht 1989). Aerosols that can impact cloud processes via CDNC include sulfates, OM, BC and sea-salt.
Aerosol radiative forcing is obtained from the difference in the radiative flux with and without aerosols calculated at each model timestep at the tropopause. Our aerosol indirect effect is calculated from differences in the total and clear sky net TOA radiative flux (SW + LW) at each model timestep. This is thus the cloudy sky radiative flux and differences between sets of simulations includes the cloudy sky forcing changes and other feedbacks. Since we allow for cloud changes, the definition does not follow that of the Intergovernmental Panel on Climate Change (2007) as also stated in Kloster et al (2008). Model simulations are run for 6 years (equilibrium simulations) and averages for the last 5 years are given, allowing for sufficient spin-up time. Although statistical significance of the results are not shown, in general, most results, above the minimum difference indicated by the white bar in the figures, are significant at the 95% level based on a student's t-test, as also shown in Bell et al (2005), and based on similar simulations with a similar version of the model.

Appendix B. Comparing simulated cloud microphysical properties with observed data from the Arctic
Ground-based remote sensing measurements for single-layer Arctic stratus clouds using data from the Atmospheric Radiation Measurement Program's site at the North Slope of Alaska (Barrow region) are available for 2000-2003(Zhao and Garrett 2008 for single-layer clouds with bases below 4 km, COD 8 and R eff 25 µm. Climatological ranges of cloud properties for June to August (maximum occurrence of liquid-phase clouds) for all single-layer thin clouds of thickness between 300 and 1200 m are 10-12 µm and 20-40 cm −3 for R eff and CDNC, respectively. Values of 3.5-4.5 and ∼30 g m −2 were obtained for COD and LWP, respectively. Retrieval uncertainties for the 95% confidence level for R eff , CDNC, COD and LWP are 7%, 30%, 15% and 10%, respectively. Although ground-based measurements cannot easily be compared with values from a coarse grid model, the general characteristics may be compared to evaluate basic features simulated by Sim PD. We use daily data for similar time periods (June to August), sampled once a day (around 10:30 am) (as in Menon et al (2008)) for the 60 • -70 • N and 145 • -165 • W region for Sim PD and restrict our analyses to single-layer clouds with bases below 4 km and for similar optical depth and cloud droplet sizes as retrievals. Simulated average values for cloud thickness ranges from 100 to 900 m (average of 310 m) for CDNC, R eff , LWP and COD are 60 cm −3 , 9.2 µm, 18 g m −2 , and 3.8, respectively. Although measurements were available at different time periods, and simulations were for a specific time period, differences between day and night measurements were relatively small compared to seasonal differences (Zhao 2007