Geoengineering: Whiter skies?
Ben Kravitz,1 Douglas G. MacMartin,2 and Ken Caldeira1
Received 9 March 2012; revised 1 May 2012; accepted 2 May 2012; published 1 June 2012.
 One proposed side effect of geoengineering with stratospheric sulfate aerosols is sky whitening during the day and afterglows near sunset, as is seen after large volcanic eruptions. Sulfate aerosols in the stratosphere would increase diffuse light received at the surface, but with a non-uniform spectral distribution. We use a radiative transfer model to calculate spectral irradiance for idealized size distributions of sulfate aerosols. A 2% reduction in total irradiance, approximately enough to offset anthropogenic warming for a doubling of CO2 concentrations, brightens the sky (increase in diffuse light) by 3 to 5 times, depending on the aerosol size distribution. The relative increase is less when optically thin cirrus clouds are included in our simulations. Particles with small radii have little influence on the shape of the spectra. Particles of radius 0.5 mm preferentially increase diffuse irradiance in red wavelengths, whereas large particles (0.9 mm) preferentially increase diffuse irradiance in blue wavelengths. Spectra show little change in dominant wavelength, indicating little change in sky hue, but all particle size distributions produce an increase in white light relative to clear sky conditions. Diffuse sky spectra in our simulations of geoengineering with stratospheric aerosols are similar to those of average conditions in urban areas today.
Citation: Kravitz, B., D. G. MacMartin, and K. Caldeira (2012), Geoengineering: Whiter skies?, Geophys. Res. Lett., 39, L11801, doi:10.1029/2012GL051652.
 Robock [2008, 2011] argued that geoengineering with sulfate aerosols would cause whitening of the sky, as was seen after the eruption of Mount Pinatubo [Robock, 2000]. Geoengineering with stratospheric sulfate aerosols to offset anthropogenic warming [Crutzen, 2006] would result in direct irradiance being decreased by much more than total irradiance, which is compensated by an increase in diffuse irradiance. When the ratio of particle size to scattered wavelength is small, such as for air molecules scattering visible wavelengths (380–780 nm), scattering is in the Rayleigh regime, in which short wavelengths (blues) are scattered more efficiently than long wavelengths (reds), giving the sky its blue hue. When particles are similar in size to the scattered wavelength, as is the case for stratospheric sulfate aerosols from large volcanic eruptions, scattering is more uniform across the visible spectrum. Thus, in the presence of an aerosol layer, diffuse irradiance at longer wavelengths (reds) increases more than the increase in shorter wavelengths (blues). As scattering increases, the portion of diffuse light also increases, resulting in a brighter sky as seen from the ground, because sulfate aerosols preferentially forward scatter more radiation than they back scatter.
 In this paper, we use the term brightness of sky to refer to diffuse irradiance integrated over the visible band from an observer looking upward, and we are specifically excluding direct irradiance of the sun in reporting this quantity. We refer to white light as having a uniform irradiance spectrum, so a whiter sky would have a higher portion of irradiance in longer wavelengths than a blue sky. All references to solar irradiance reduction are in irradiance integrated over the visible band as measured at the surface.
1 Department of Global Ecology, Carnegie Institution for Science, Stanford, California, USA.
2 Department of Control and Dynamical Systems, California Institute of Technology, Pasadena, California, USA.
Corresponding author: B. Kravitz, Department of Global Ecology, Carnegie Institution for Science, 260 Panama St., Stanford, CA 94305, USA. (firstname.lastname@example.org)
 To determine the degree to which the sky will whiten and brighten as a result of idealized distributions of stratospheric sulfate aerosols, we use the radiative transfer model libRadtran [Mayer and Kylling, 2005] (see Text S1 in the auxiliary material) to simulate spectral irradiance.1 Because the scattering properties of sulfate aerosols are strongly dependent upon the size distribution, we compare the results for lognormal distributions with different median radii and standard deviations, the latter of which is a measure of distribution width.We investigate atmospheric loadings of aerosols that correspond to up to a 2% reduction in irradiance received at the surface, which is roughly the amount needed to offset the warming from a doubling of CO2 concentrations from the preindustrial era [Govindasamy and Caldeira, 2000].
3. Spectral Irradiance and Mie Calculations
 Figure 1 shows perturbations to the visible light spectrum from our simulations of stratospheric sulfate aerosols (unimodal lognormal, median radius rg = 0.5 mm, s = 0.1) with an aerosol optical depth that reduces total integrated visible band surface irradiance by 2%. Relative to a clear sky case (no clouds or aerosols), a 2% reduction in total irradiance with this aerosol size distribution reduces daytime direct irradiance by 22% and increases diffuse irradiance by over four times. The increase in diffuse irradiance is more prominent in red wavelengths than blue wavelengths.
 For our daytime simulations, we use a solar zenith angle of 0. We examined the impacts of changing solar zenith angle on the aerosol effects on irradiance spectra. Except for large solar zenith angles (i.e., near sunrise or sunset), increasing the angle reduces irradiance by a factor
that is nearly constant at all wavelengths, so daytime sky hue changes from stratospheric aerosols are relatively insensitive to solar zenith angle (not shown). At large solar zenith angles, the sun appears redder than in daytime, as the amount of atmosphere through which light must travel to reach an observer increases. This increases the amount of Rayleigh scattering by air molecules, which preferentially scatters shorter wavelengths, and consequently longer wavelengths (reds and oranges) reach the observer as a higher portion of light [Corfidi, 1996]. Therefore, direct solar radiation at sunset appears to be redder, whereas the diffuse sky retains its blue hue (Figure 1, bottom). However, in the presence of a stratospheric aerosol layer, solar radiation is scattered by the aerosol layer, resulting in an afterglow that is characteristic of volcanic sunsets [Corfidi, 1996]. Indeed, our simulations for the eruption of Pinatubo (described in Section 3 of Text S1) show an increase in diffuse light by approximately a factor of two at long wavelengths; direct irradiance decreases by up to 75% at long wavelengths. Similar to the eruption of Pinatubo, sulfate aerosols under the specifications we simulated (rg = 0.5 mm, s = 0.1) would cause sunsets similar in appearance to volcanic sunsets.
 Changes in irradiance are dependent upon the size distribution of the particles, both in terms of median radius and distribution width. Figure 2 shows spectral irradiance and asymmetry parameters for a range of median radii and distribution standard deviations. The asymmetry parameter g indicates the average direction of light scattering. Pure backscattering corresponds to g = 1, pure forward scattering to g = 1, and isotropic scattering to g = 0. Note that g varies with wavelength in ways that depend on the particle size distribution. All cases show at least a factor of 2 increase in diffuse irradiance over the clear sky case, indicating a sky brightening.
 The distribution with median radius rg = 0.1 mm shows an asymmetry parameter that is almost uniform in wavelength, meaning the same approximate percentage of diffuse light will be forward scattered at all wavelengths. Because the solar spectrum has more irradiance in shorter wavelengths than longer ones, the diffuse spectrum resulting from stratospheric sulfate aerosols with this particular size distribution has a larger increase in irradiance in blue wavelengths than red wavelengths. The results for rg = 0.1 and 0.3 mm are similar. For rg = 0.5 mm, the asymmetry parameter favors forward scattering of longer wavelengths, resulting in a flattening of the diffuse spectrum. This whitens the sky, where white light is defined as having uniform irradiance for all visible wavelengths. rg = 0.7 mm has an asymmetry parameter that favors forward scattering of both short and long visible wavelengths, and rg = 0.9 mm favors forward scattering of short wavelengths. The latter results in the smallest change in the overall shape of the diffuse spectrum.
 The asymmetry parameter depends more strongly on wavelength for narrower size distributions. For a distribution median radius of rg = 0.5 mm, as the distribution width increases, the asymmetry parameter becomes more uniform, slightly favoring forward scattering of shorter wavelengths. The asymmetry parameter also increases (more forward scattering) for all wavelengths when s changes from 0.5 to 1.0, which skews diffuse irradiance toward shorter wavelengths.
 Table 1 shows increases in diffuse sky brightness (diffuse irradiance integrated over the visible band) for all simulations we performed. Larger reductions in solar irradiance are achieved by larger optical depths (Table 2 in Text S1), so large reductions have more sky brightening. Increasing s from 0.1 to 0.25 decreases brightness for small distribution median radii, and increasing s from 0.25 to 0.5 or 1.0 increases brightness. This can be explained by integrating under the irradiance curves in Figure 2 (top right). For large median radii, the increase in brightness is monotonic with s.
 The dependence of brightness on the median radius of the distribution strongly depends upon the distribution width (Table 1). For narrow distributions (s = 0.1), the radius that gives the largest increase in brightness is rg = 0.3 mm, and the smallest increase is for rg = 0.8 mm. For s ≥ 0.25, the largest brightness increase is for the smallest radius of rg = 0.1 mm, and brightness increases diminish almost monotonically with increasing radius.
 We also performed simulations of thin cirrus clouds to determine the effects of these clouds on our simulations (see Section 4 of Text S1). A cloud fraction of 20% [Mace et al., 2001] diminishes brightness increases by approximately 40% for all size distributions (Table 3 in Text S1). Effects of the geoengineering case rg = 0.5 mm, s = 0.1 show nearly negligible changes in total irradiance, up to a 5% decrease in direct irradiance, and over a 20% increase in diffuse irradiance (Figure 1 in Text S1). The enhancement of diffuse light by cirrus clouds is by an approximately constant factor across all wavelengths, whereas stratospheric sulfate aerosols preferentially enhance longer wavelengths over shorter ones.
4. Perceived Color
 Determining whether the sky would be whiter under geoengineering is not straightforward. The main difficulty in determining this effect is that color is not an intrinsic property of light; rather, it is perceived and interpreted. Therefore, the relevance of showing color swatches of sky color in red-green-blue color space (or another equivalent) is not immediately apparent, as it is difficult to interpret these color swatches out of their normal perceptual context. Sky color temperature varies with sunlight and meteorological conditions, and many different color matching algorithms have been proposed, all of which can radically affect color perception. (See Section 7 of Text S1 for a discussion and examples.) Furthermore, there is no guarantee the colors will be displayed with fidelity, as color appearance depends on the means of viewing, e.g., the computer monitor and printer.
 Despite these difficulties, one can obtain some measure of sky color change. A metameric match is the monochromatic wavelength which will elicit the same response in the retinas as a specified irradiance spectrum. Using an average profile of cone sensitivities [Stockman and Sharpe, 2000] and matching functions [Smith, 2005] (see Section 6 of Text S1), we determined the metameric matches to spectra produced by our simulations. Additionally, the perceived sky can be approximated as a mixture of this monochromatic light with white light (light with a uniform irradiance
 Figure 3 illustrates these two concepts for our simuations, assuming a 2% reduction in visible integrated solar irradiance and a narrow size distribution to highlight the effects of radius. The black lines in Figure 3 show the case for a clear sky (no aerosols). Small radii show the largest wavelength departure, reaching a maximum difference of 8 nm. rg = 0.8 mm shows the smallest departure by only 1 nm. All of the wavelengths in Figure 3 fall distinctly within the blue range [Boynton et al., 1964], so none of these departures is indicative of a radical change in sky color. However, all simulations pictured in Figure 3 show a whitening of the sky. The white to monochromatic light ratio increases by over a factor of 2 for all cases, and by a full order of magnitude for rg = 0.5 and 0.6 mm. This figure shows that aerosols in the middle of the range we considered would likely cause the largest changes in perceived color. For reference, we have also included a calculation for the 1991 eruption of Mount Pinatubo (red asterisks, see Section 3 of Text S1 for specifications of the simulation). Nearly all of our geoengineering simulations that reduce total irradiance by 2% show more sky whitening than was seen in March after the eruption of Pinatubo.
5. Comparisons: Urban Pollution
 To provide context for our results, we compare the spectra from our stratospheric sulfate aerosol geoengineering simulations with spectra characteristic of a range of urban air pollution conditions. Since pollution in urban areas varies on a daily basis, depending upon meteorological conditions, as well as the particular urban area, we compared our results to an average of aerosol measurements (see Table 1 in Text S1) taken in three urban areas (Greenbelt, MD; Paris, France; and Mexico City, Mexico), including a low aerosol case (relatively clear day), a high aerosol case (relatively polluted day), and an average/typical case [Dubovik et al., 2002] (see Section 5 of Text S1). Stratospheric aerosols with a size
distribution specified by rg = 0.5 mm and s = 1.0 and with a reduction of visible integrated solar irradiance by 1.5% have a similar spectrum to an average pollution day in the urban areas considered (Figure 4). A day with heavy pollution would correspond to a level of stratospheric sulfate aerosols that would produce a solar irradiance reduction far exceeding 2%. For reference, we have also included a simulation of the 1991 eruption of Mount Pinatubo (specifications in Section 3 of Text S1).
 For stratospheric sulfate aerosol distributions with rg = 0.5 mm and s = 1.0, an optical depth of 0.15 results in a solar irradiance reduction of 1%, and reduction in irradiance scales linearly with optical depth (Table 2 in Text S1). Therefore, assuming the average urban pollution case has an optical depth of t550 = 0.21, an additional solar reduction of 1% from stratospheric aerosols would, on average, cause the sky over formerly pristine areas to look similar to the sky over urban areas. The variation of optical depth in urban areas that already occurs is much larger than the added optical depth from our simulations.
 According to our simulations, stratospheric geoengineering with sulfate aerosols at the levels considered here likely would not cause profound changes in sky hue but would whiten the sky, potentially exceeding the amount of whitening observed after the 1991 eruption of Mount Pinatubo. These levels of geoengineering would also cause sunsets similar to those seen after large volcanic eruptions.
 The results of this study have multiple implications, although some are not easily quantified. One of the motivations for Robock  is the concern over potential psychological implications. We have provided the physics-based predictions of sky changes from stratospheric sulfate geoengineering, and our results can provide a basis for psychological research. Additionally, public perception of the effects discussed in this paper could be investigated, as the whitening of skies could contribute to opposition to geoengineering.
 One of the more quantifiable impacts is a potential enhancement of the land carbon sink, as was seen after the 1991 eruption of Mount Pinatubo [Mercado et al., 2009]. Diffuse light more easily penetrates canopies, and it can also more easily strike more leaf surfaces, which can increase plant photosynthetic activity. However, the degree to which this effect could occur depends on the size distribution of the sulfate aerosols, as our results show.
 Additionally, solar power generating capacity could be affected by geoengineering [Murphy, 2009]. Generation of solar power by concentrating solar thermal plants is largely proportional to the amount of direct sunlight received. Therefore, the reductions in direct irradiance that our results show could reduce the effectiveness of this particular method of obtaining low-carbon-emission energy, which could lead to increased use of other energy sources, including fossil fuel combustion.
 Although we do not discuss photochemical impacts, Tilmes et al.  have shown that sulfate aerosol geoengineering has the potential to decrease stratospheric ozone concentrations. Possible impacts of sulfate geoengineering on chemistry warrant further study.
 Different means of generating the aerosol particles would produce different aerosol size distributions, affecting development of engineering approaches to geoengineering. If the aerosols are formed from photochemical conversion of SO2 into sulfate aerosols, the expected particle size distribution width is large [Heckendorn et al., 2009], meaning the median particle radius has diminished importance in determining the amount of sky whitening, but sky brightening will be larger than for narrower distributions. However, if a monodisperse size distribution of aerosols is created, which could result from direct condensation of gaseous sulfuric acid [Pierce et al., 2010], median radius becomes important in determining sky whiteness. Our results also have implications for geoengineering with engineered particles, in that the sky color and brightness will depend very strongly on the asymmetry parameter of the particle.
 To achieve a minimal amount of sky whitening under geoengineering with stratospheric sulfate aerosols, there are three options: (1) choose a larger particle size so the particles preferentially forward scatter more in blue wavelengths than red wavelengths (the shape of g will be similar to that
of r = 0.9 mm in Figure 2); (2) Use particles that have less forward scattering in all wavelengths (smaller value of g); or (3) Limit the amount of sulfate geoengineering or avoid such geoengineering entirely. The first option is problematic, as large aerosols have a higher fall speed than smaller aerosols, meaning a higher rate of injection is needed for these large aerosols than for smaller ones to achieve the same amount of atmospheric mass loading. The second is problematic, as can be seen from Figure 2, because all size distributions of sulfate aerosols considered scatter prominently in the forward direction. Thus, an obvious path to avoiding sky whitening is to greatly limit or avoid entirely the intentional introduction of sulfate particles into the stratosphere.
 Acknowledgments. We thank Alan Robock for suggesting this topic and for helpful comments throughout the process of drafting this paper, as well as the reviewers for their comments.
 The Editor thanks the anonymous reviewers for assisting with the evaluation of this paper.
Boynton, R. M., W. Schafer, and M. E. Neun (1964), Hue-wavelength relation measured by color-naming method for three retinal locations, Science, 146, 666–668.
Corfidi, S. F. (1996), The colors of twilight, Weatherwise, 49(3), 14–19, doi:10.1080/00431672.1996.9925403.
Crutzen, P. J. (2006), Albedo enhancement by stratospheric sulfur injections: A contribution to resolve a policy dilemma? Climatic Change, 77, 211–219, doi:10.1007/s10584-006-9101-y.
Dubovik, O., et al. (2002), Variability of absorption and optical properties of key aerosol types observed in worldwide locations, J. Atmos. Sci., 59, 590–608.
Govindasamy, B., and K. Caldeira (2000), Geoengineering Earth’s radiation balance to mitigate CO2-induced climate change, Geophys. Res. Lett., 27, 2141–2144, doi:10.1029/1999GL006086.
Heckendorn, P., et al. (2009), The impact of geoengineering aerosols on stratospheric temperature and ozone, Environ. Res. Lett., 4, 045108, doi:10.1088/1748-9326/4/4/045108.
Mace, G. G., E. E. Clothiaux, and T. P. Ackerman (2001), The composite characteristics of cirrus clouds: Bulk properties revealed by one year of continuous cloud data, J. Clim., 14, 2185–2203.
Mayer, B., and A. Kylling (2005), Technical note: The libRadtran software package for radiative transfer calculations—Description and examples of use, Atmos. Chem. Phys., 5, 1855–1857, doi:10.5194/acp-5-1855-2005.
Mercado, L. M., N. Bellouin, S. Sitch, O. Boucher, C. Huntingford, M. Wild, and P. M. Cox (2009), Impact of changes in diffuse radiation on the global land carbon sink, Nature, 458, 1014–1017, doi:10.1038/ nature07949.
Murphy, D. M. (2009), Effect of stratospheric aerosols on direct sunlight and implications for concentrating solar power, Environ. Sci. Technol., 48(8), 2784–2786, doi:10.1021/es802206b.
Pierce, J. R., D. K. Weisenstein, P. Heckendorn, T. Peter, and D. W. Keith (2010), Efficient formation of stratospheric aerosol for climate engineering by emission of condensible vapor from aircraft, Geophys. Res. Lett., 37, L18805, doi:10.1029/2010GL043975.
Robock, A. (2000), Volcanic eruptions and climate, Rev. Geophys., 38, 191–219, doi:10.1029/1998RG000054.
Robock, A. (2008), 20 reasons why geoengineering may be a bad idea, Bull. At. Sci., 64, 14–18, doi:10.2968/064002006.
Robock, A. (2011), Bubble, bubble, toil and trouble: An editorial comment, Clim. Change, 105, 383–385, doi:10.1007/s10584-010-0017-1.
Smith, G. S. (2005), Human color vision and the unsaturated blue color of the daytime sky, Am. J. Phys., 73, 590–597, doi:10.1119/1.1858479.
Stockman, A., and L. T. Sharpe (2000), The spectral sensitivites of the middle and long-wavelength-sensitive cones derived from measurements in observers of known genotype, Vision Res., 40, 1711–1737.
Tilmes, S., R. Müller, and R. Salawitch (2008), The sensitivity of polar ozone depletion to proposed geoengineering schemes, Science, 320,1201–1204, doi:10.1126/science.1153966.
Der Geo-Engineering Himmel heute..