Abstract
We have explored the direct and indirect radiative effects on climate of organic particles likely to have been present on early Earth by measuring their hygroscopicity and cloud nucleating ability. The early Earth analog aerosol particles were generated via ultraviolet photolysis of an early Earth analog gas mixture, which was designed to mimic possible atmospheric conditions before the rise of oxygen. An analog aerosol for the present-day atmosphere of Saturn's moon Titan was tested for comparison. We exposed the early Earth aerosol to a range of relative humidities (RHs). Water uptake onto the aerosol was observed to occur over the entire RH range tested (RH = 80–87%). To translate our measurements of hygroscopicity over a specific range of RHs into their water uptake ability at any RH < 100% and into their ability to act as cloud condensation nuclei (CCN) at RH > 100%, we relied on the hygroscopicity parameter κ, developed by Petters and Kreidenweis. We retrieved κ = 0.22 ± 0.12 for the early Earth aerosol, which indicates that the humidified aerosol (RH < 100 %) could have contributed to a larger antigreenhouse effect on the early Earth atmosphere than previously modeled with dry aerosol. Such effects would have been of significance in regions where the humidity was larger than 50%, because such high humidities are needed for significant amounts of water to be on the aerosol. Additionally, Earth organic aerosol particles could have activated into CCN at reasonable—and even low—water-vapor supersaturations (RH > 100%). In regions where the haze was dominant, it is expected that low particle concentrations, once activated into cloud droplets, would have created short-lived, optically thin clouds. Such clouds, if predominant on early Earth, would have had a lower albedo than clouds today, thereby warming the planet relative to current-day clouds. Key Words: Archean—Early Earth—Organic haze—Tholin. Astrobiology 11, 135–149.
1. Introduction
A
A simple radiative balance calculation indicates that the limited amount of solar radiation received by early Earth would have resulted in a frozen planet if the atmosphere had the same greenhouse effect as it has today (Sagan and Mullen, 1972; Sagan and Chyba, 1997).Yet studies of uranium, lead, and oxygen isotopes in ancient zircons (Mojzsis et al., 2001; Wilde et al., 2001) have indicated the presence of ubiquitous liquid water over 4 billion years ago. There is also ample fossil evidence that life, which is presumably reliant on liquid water, existed throughout the Archean (Knoll and Barghoorn, 1977; Walsh, 1992; Schopf et al., 2007). One solution to this so-called “faint young Sun paradox” is the presence of larger than current-day concentrations of greenhouse gases. From paleosol record constraints, it is believed that 2.8 billion years ago carbon dioxide (CO2) was present in concentrations approximately 20 to no more than 100 times higher than concentrations in the present-day atmosphere (Rye et al., 1995; Sheldon, 2006). However, recent studies of banded iron formations suggest that CO2 concentrations were no more than 3 times higher than those of the present (Rosing et al., 2010). Despite the discrepancies in the geological data regarding Archean CO2 concentrations, one fact remains clear: according to models, if CO2 were the only greenhouse gas, even the upper limit of CO2 concentrations suggested by paleosol data would not have provided enough warming to offset the faint young Sun (Pavlov et al., 2000; Haqq-Misra et al., 2008). Increased concentration of atmospheric methane (CH4) has been suggested as an additional source of greenhouse warming. During the Archean, CH4 may have existed in higher concentrations than those of the present day due to methanogenic bacteria on early Earth and the relatively long lifetime for CH4 in an anoxic atmosphere (Haqq-Misra et al., 2008). Models of early Earth 2.8 billion years ago have indicated that <0.03 bar CO2 and a mixing ratio of at least 10−3 CH4 in a nitrogen (N2)-based atmosphere could have provided enough warming to boost the global average surface temperature above 273 K (Haqq-Misra et al., 2008). When the radiative effects of spherical organic particles are included in the models, however, the warming is offset by the antigreenhouse effect of the haze, and the paradox persists (Haqq-Misra et al., 2008). More recently, two new solutions to the problem have been suggested. Rosing et al. (2010) suggested that, despite their new greenhouse-restrictive constraint on CO2 levels, Earth could have stayed warm because the surface albedo would have been less than it is at present due to the smaller continental area. They also suggested that clouds would have had lower albedo because, with fewer cloud condensation nuclei (CCN), the droplets would have been larger. These two effects could have maintained the Earth above freezing. Alternatively, the fractal organic aerosols suggested by Wolf and Toon (2010) may have provided a UV shield, which would have allowed NH3 to build up to 1–10 ppmv. Such levels could also keep Earth above freezing even with CO2 and CH4 at only 1000 ppmv levels.
Titan, a moon of Saturn, may act as an analogue to early Earth because its stratosphere is mostly N2 with 1.5% CH4 (Flasar et al., 2005; Niemann et al., 2005). A ubiquitous organic haze is formed from the photolysis and electron-dissociation of CH4 and N2. When CO2 is added to analog atmospheres—either in the laboratory or in photochemical models—a haze forms from a gas mixture with a ratio of CH4/CO2 = 0.1 or higher (Domagal-Goldman et al., 2008; DeWitt et al., 2009). On Titan, the direct effect of the haze, which is to scatter and absorb sunlight, lowers the surface temperature on Titan by approximately 10% or 9 K, compared to a scenario without the haze (McKay et al., 1991). A ubiquitous haze on early Earth would have had profound consequences on climate in terms of both its direct and indirect effects. Several laboratory and modeling studies have investigated the direct effect of possible early Earth aerosol (Khare et al., 1984; Haqq-Misra et al., 2008; Wolf and Toon, 2010), though none have investigated how humidified aerosols would have altered the direct effect on early Earth.
The effects of early Earth aerosol on cloud radiative properties, termed the indirect effect of the aerosol, have not been previously investigated in any laboratory study nor explicitly explored by models. Yet, there are compelling reasons to investigate the cloud-forming ability of these aerosol particles. First, the haze layer formed by these aerosol particles would likely have been nearly global in extent (Wolf and Toon, 2010). Second, many of the sources of CCN in our present-day atmosphere would not have existed, or would have been greatly reduced, in the early Earth atmosphere. Most CCN on current-day Earth are sulfates, which should have been rare on early Earth, or sea salt aerosol, which dominate in the marine boundary layer but usually exist no more than 2 km in altitude (Seinfeld and Pandis, 1998). Therefore, an organic aerosol haze, generated from UV photolysis of CH4, may have dominated the CCN in some regions of the Archean atmosphere. Specifically, the lack of other CCN sources would have been due to the absence of oxygen, which precludes most current-day natural biogenic and anthropogenic aerosols, and the lack of landmass, which supplies dust aerosol to the current-day atmosphere. Typical CCN concentrations in our current atmosphere are on the order of 100–1000 cm−3, depending upon whether the air mass is maritime or continental (Rogers and Yau, 1996). Estimates of non-anthropogenic or “natural” concentrations over the current-day oceans range from 50–200 CCN cm−3, due largely to products of biological activity and inorganics (Andreae and Rosenfeld, 2008). In the most remote areas of the present world, where anthropogenic and oxygen-reliant biogenic aerosol influences are minimized—perhaps reminiscent of aerosol concentrations on early Earth—typical concentrations are only a few tens per cubic centimeter (Andreae and Rosenfeld, 2008). Kump and Pollard (2008) showed that fewer CCN during the Cretaceous may have decreased cloud coverage and, hence, decreased Earth's effective albedo, enough to warm the planet significantly. Recently, Rosing et al. (2010) modeled the effects of fewer CCN during the Archean than at present and showed that the resultant planetary albedo would be low enough to provide enough warming to offset the faint young Sun. Rondanelli and Lindzen (2010) showed that increased cirrus clouds in the tropics could have offset a faint young Sun but mentioned no mechanism for their formation. Recent work by Goldblatt and Zahnle (2010) pointed out that a better understanding of the mechanisms by which early Earth clouds formed is needed before the role clouds played in radiative forcing on early Earth can be properly constrained.
Investigations into the hygroscopicity and ability of early Earth analog aerosols to serve as CCN are especially compelling projects, given the important role humidified aerosols and clouds must have played in both the radiative balance and hydrological cycle of our early planet. The present study is the first to report hygroscopicity measurements of an early Earth analog aerosol. From these measurements, we have calculated how the aerosol will take up water at any relative humidity (RH) (above or below 100%), which thus allowed us to estimate the humidified aerosol's inferred direct (i.e., antigreenhouse effect) and indirect (i.e., CCN ability) effects. We measured the hygroscopicity of a Titan analog aerosol for comparison.
2. Experimental Methods
A tandem cavity ringdown aerosol extinction spectrometer (CRD-AES) was used to measure the optical growth factor, f(RH), of early Earth and Titan analog aerosol for a range of particle diameters. The optical growth is defined as
where

Schematic of aerosol generation and the humidified tandem cavity ringdown aerosol extinction system. Abbreviations are as follows: CPC, condensation particle counter; CRD, cavity ringdown aerosol extinction spectrometer; DMA, differential mobility analyzer; RH, relative humidity; PMT, photomultiplier tube.
2.1. Description of aerosol generation
The hygroscopicity of both early Earth and Titan analog aerosols were studied. The early Earth analog aerosol was produced by UV photolysis of an early Earth atmospheric gas mixture of 0.1% ultrahigh-purity CH4 and 0.1% research-grade CO2 with an ultrahigh-purity N2 background. This mixture has a C/O ratio of 1 and is considered to be a plausible early Earth analog atmosphere, in terms of both the absolute concentrations of CH4 and CO2 and C/O ratio (Pavlov et al., 2000; Haqq-Misra et al., 2008). The CO2 concentrations, however, are a thousand times larger than the new constraints based on banded iron formations (Rosing et al., 2010).
The Titan analog aerosol was generated by the photolysis of a mixture of 0.1% ultrahigh-purity CH4 in an ultrahigh-purity N2 background. Past work has shown that photolysis of this gas mixture results in aerosols of similar composition to photolysis of 2% CH4 in N2, similar to conditions on Titan (Trainer et al., 2006). Also, because the Titan analogue has an identical mixing ratio of CH4 to our terrestrial analogue, it can provide an interesting comparison to the early Earth analogue. From previous studies in our laboratory (Trainer et al., 2006), we know that the early Earth analog aerosol generated from the early Earth gas mixture has several oxygenated species incorporated into it, which the Titan analogue lacks. Using analog gas mixtures with identical mixing ratios of CH4 illustrates the impact on hygroscopicity of including an oxygenated species in an analog gas mixture. Both analog gas mixtures used in this study were chosen for optimal aerosol production with regard to the CH4/N2 ratio for the Titan aerosol and the C/O ratio for the early Earth aerosol. These analog gas mixtures, therefore, provide the maximum possible signal for our experiments.
Images of the generated aerosols were obtained through tapping-mode atomic force microscopy on a hydrophilic silica substrate (SiOx/Si(100)) and are shown in Fig. 2. The aerosols in both panels range from tens of nanometers in diameter to roughly 100 nm. This particle size is consistent with both previous work in this laboratory and observations of Titan haze (Tomasko et al., 2005, 2008). In both images in Fig. 2, the analog aerosols appear spherical, which is in agreement with previous work that used other techniques to evaluate sphericity of our analog aerosols (Trainer et al., 2006). The spherical shape of the aerosols allows them to be treated as spherical Mie scatterers. The complex refractive indices of both these aerosols, which is necessary for analyzing the optical growth results at λ = 532 nm, were determined previously in Hasenkopf et al. (2010). Our method of aerosol generation was described in greater detail by Hasenkopf et al. (2010) and Trainer et al. (2006).

Atomic force microscope images of (
2.2. Description of f(RH) experiments
After the aerosol exited the UV formation chamber, the aerosol flow rate was increased to a total of 300 standard cubic centimeters per minute (sccm) by adding a dilution flow of pre-purified N2. The flow of dry, polydisperse aerosol then entered a differential mobility analyzer (TSI Electrostatic Classifier Model 3080, DMA), where it was surrounded by a sheath flow of 1500 sccm of dry N2. The DMA was used to size select the particles by their electrostatic mobility. The aerosol exiting the DMA at a given voltage was dominated by a monodisperse population of single-sized particles but also included monodisperse populations of more highly charged particles of larger sizes. For our experimental setup, the number of doubly charged particles was small but significant enough to affect the extinction coefficient measured by the CRD-AES. The number of particles at higher charges than +2 is negligible.
To account for the doubly charged particles' effect on extinction, an optical effective diameter D
eff was employed, as has been done in previous optical growth studies (Garland et al., 2007; Beaver et al., 2008). The effective diameter D
eff is the equivalent diameter of identically sized particles that yield the same extinction coefficient b
ext for the same concentration C (particles cm−3) and have the same refractive indices as the polydisperse population. To calculate D
eff for a given selected particle size, the effective extinction cross section
After the aerosol particles were size-selected in the DMA, they entered the tandem cavity ringdown aerosol extinction spectrometers with humidity control. The details of the CRD-AES have been described previously in Baynard et al. (2007), Garland et al. (2007), and Beaver et al. (2008) and thus are only briefly described here. The body of each CRD-AES cell is a 0.85 m stainless steel tube. Highly reflective (reflectivity >99.998%) mirrors (Advanced Thin Films) are mounted on both ends of the tube. The second harmonic, λ = 532 nm, from a neodymium-doped yttrium aluminum garnet laser (Nd:YAG, Big Sky Laser Technologies, Quantel USA), is pulsed into the cavity and reflects back and forth many times on the mirrors. In a cell with only the background gas (devoid of particles) and over a characteristic timescale τ
0, the intensity of the light falls to 1/e of its initial intensity. A typical τ
0 for our system was 100 μs, which corresponds to an effective pathlength of 30 km. When the cell was filled with aerosol, the ringdown time was decreased to a new ringdown time, τ. By comparing the two ringdown times, a measure of the extinction coefficient b
ext can be made (Pettersson et al., 2004):
where R L is the ratio of the total cavity length (0.85 m) to the cavity length occupied by aerosol sample (0.75 m) and c is the speed of light.
In the CRD-AES, the aerosol particles first travel into a CRD-AES cell, in which the extinction coefficient of the dry aerosol (RH ≤ 13%),
Before each experiment was conducted, it was ensured that there were no particle losses between the cells and that the cells responded identically. Particles at a chosen diameter were flowed through both cells, where the RH was held constant and low (RH < 13%). The cells are considered to have identical response when their extinctions match within error of the extinction measurements due to fluctuations in particle concentrations, which is generally within 2% of each other. The similarity of the extinction between the cells indicates that negligible particle loss has occurred and that the cells have a similar response to the particles.
3. Results
3.1. Experimental validation using ammonium sulfate
To verify the system, we measured the f(RH) of ammonium sulfate particles at nine diameters between 150 and 500 nm at RH = 80%. The particles were generated via atomization of a 10% ammonium sulfate-HPLC water solution. A syringe pump (Harvard Apparatus Model 70-2208) injected the solution into an atomizer (TSI Model 3076) with pre-purified N2 as the carrier gas. The aqueous particles were then flowed through several driers to ensure that dry particles entered the CRD-AES. Ammonium sulfate aerosol generation and characterization is described in more detail in Hasenkopf et al. (2010).
Figure 3 shows the results from the test case. Each f(RH) value at a given diameter is an average of at least three data points, each data point itself an average of extinction measurements sampled every 10.3 s for 2 min. Each f(RH) value was calculated from data obtained on at least two different days for 79% < RH < 81%. The error bars in Fig. 3 indicate the uncertainty in RH of ± 3%. Figure 3 also compares the test case values with previous ammonium sulfate f(RH) studies at RH = 80% (Garland et al., 2007; Beaver et al., 2008) and theoretical predictions calculated from the extended Aerosol Inorganic Model (e-AIM) (

f(RH = 80%) of ammonium sulfate for a range of particle diameters from different studies [this study = black circles; Beaver et al. (2008) = gray triangles; and Garland et al. (2007) = gray squares]. The solid black line represents the f(RH) predicted for ammonium sulfate using e-AIM and Mie theory. The error bars represent error based on the uncertainty in the measured RH.
The geometric growth factor (G
f) of ammonium sulfate at RH = 80% can also be compared with other measurements in the literature. The G
f is the ratio of the humidified particle diameter over the dry particle diameter. To convert the measured f(RH) values to G
f, Mie theory is applied. The humidified particle's extinction coefficient
where D
wet is the diameter of the wet particle, N is the concentration of particles (cm−3), and
where
Table 1 shows a comparison of G f values at RH = 80% determined from this study with several others. As with the f(RH) measurements shown in Fig. 3, the G f calculated from this study is comparable to others that have used systems similar to the one used in this work (Garland et al., 2007; Beaver et al., 2008) and also to those that have used different techniques (Petters and Kreidenweis, 2007). Calculating G f for 77% < RH < 83% using the e-AIM, which corresponds to the error at RH = 80% of the RH probes, we find that our measured G f falls within the calculated range.
The italicized values were calculated from the boldface values in the same row.
This work.
Beaver et al. (2008)—growth factor is calculated from f(RH = 80%) measurements.
Garland et al. (2007)—growth factor is calculated from f(RH = 80%) measurements.
Tang (1996).
Petters and Kreidenweis (2007)—upper and lower κ-values determined from growth factor and CCN data.
e-AIM calculation (
Brooks et al. (2004)—growth factor is measured from humidified tandem DMA data.
Svenningsson et al. (2006)—growth factor is measured by humidified tandem DMA data.
3.2. f(RH) and G f at RH = 80% for early Earth and Titan analog aerosols
Figure 4 shows the f(RH = 80%) of early Earth and Titan analog aerosols. An RH of 80% was chosen to provide comparison with previous optical growth studies of ammonium sulfate and slightly soluble organics (namely, phthalic and pyromellitic acids) also performed at RH = 80% (Garland et al., 2007; Beaver et al., 2008). In addition, within the detection limits of the system, it is difficult to obtain enough difference in the dry versus the humidified signal of these aerosols for studies conducted much lower than RH = 80%. The size range of the diameters tested was dependent upon available signal, which is affected by both aerosol concentration and the particle's optical properties. Like the ammonium sulfate data, each f(RH) value at a given diameter is an average of at least three sets of data points, each data point itself an average of extinction measurements sampled every 10.3 s for 2 min. As in the ammonium sulfate test case, each f(RH) value was calculated from data points obtained on at least two different days for 79% < RH < 81%. The error bars indicate the standard deviation among the data points used.

Optical growth f(RH) of early Earth and Titan analog aerosols of various particle diameters at RH = 80% (this work) compared with other slightly soluble organics (Beaver et al., 2008). The dotted line at f(RH) = 1 indicates no growth due to water uptake on the particles. The error bars represent the standard deviation of repeated measurements.
Both analog aerosols have f(RH) > 1, exhibiting water uptake at RH = 80%. The f(RH) of the early Earth analogue is similar to phthalic and pyromellitic acids, both of which are aromatic dicarboxylic acids. Though the refractive indices of phthalic and pyromellitic acids are not reported in the literature, they presumably have very similar real and imaginary refractive indices, since they have similar chemical structures. The chemical structure and measured G f at RH = 80% of phthalic acid [assuming n = 1.60 and k = 0.000, as done in Beaver et al. (2008)] is shown in Table 1. The similarity between phthalic and pyromellitic acids and the early Earth analog aerosol in terms of hygroscopicity is qualitatively consistent with previous chemical composition studies of the early Earth analog aerosol that identified dicarboxylic acids and aromatics in aerosol mass spectra of analog aerosol generated in an identical manner to this study (Trainer et al., 2006).
The f(RH) versus diameter curves in Fig. 4 depend upon G
f and
For a particle of a given composition over the range of RHs studied, G f is a strong function of RH and a very weak function of dry particle diameter for particles in the size range studied in this work. The weak dependence of G f on dry particle diameter is due to the negligible importance of the Kelvin effect for particles greater than D dry ∼ 70 nm (Lewis, 2006; Gao et al., 2007). The error in G f due to treating it as a constant with particle diameter is approximately 2% for the smallest sizes studied (D dry = 125 nm), which is smaller than experimental uncertainties (Lewis, 2006).
The deviation of f(RH) from a straight horizontal line in Eq. 5 is due to the
The G f at RH = 80% is reported in Table 1 for early Earth and Titan analog aerosols. The G f reported is the average G f measured at all sizes, and the error reported is the standard deviation. The conversion from f(RH) to G f was performed in the same manner as G f calculations for ammonium sulfate in Section 3.1. The G f at RH = 80% for early Earth and Titan are 1.17 ± 0.06 and 1.08 ± 0.01, respectively. From chemical considerations, it is not surprising that the early Earth aerosol particles are slightly more hygroscopic than the Titan aerosol. Trainer et al. (2006) found that the early Earth analog aerosol contains significant amounts of oxygen. It is expected from theoretical considerations, laboratory studies (Kotzick and Niessner, 1999), and field campaigns (Jimenez et al., 2009) that a more oxidized version of a given organic will be more hygroscopic.
3.3. Hygroscopicity parameter κ at RH = 80% for early Earth and Titan analog aerosols
The geometric growth factors can be used to calculate the hygroscopicity parameter κ. This parameter describes a particle's behavior in both sub- and supersaturated conditions based on its chemical composition. The hygroscopicity parameter enables investigators to take our measurements, performed in subsaturated conditions at one RH or over a range of RHs, and calculate an analog aerosol's direct effect on the radiative balance of the early Earth atmosphere at any RH and also its ability to act as CCN. This parameter was first developed by Petters and Kreidenweis (2007) and can be used to express what they term as “κ-Kohler theory,” which describes the RH (i.e., the saturation ratio S over an aqueous solution droplet:
where D is the wet droplet diameter, D dry is the dry diameter of the dry aerosol, σ s/a is the surface tension of the solution/air interface, M w is the molecular weight of water, R is the universal gas constant, T is the temperature, and ρ w is the density of water. As Petters and Kreidenweis (2007) described, the use of the hygroscopicity parameter κ circumvents the need to define aerosol properties such as density, molecular weight, or dissociation constants. For reference, pure water has a κ-value of 0, and an extremely hygroscopic compound such as sodium chloride has a κ-value of 1.4 (Petters and Kreidenweis, 2007). For our calculations, we assume σ s/a = 0.072 J m−2, which is the surface tension of water. Though the actual σ s/a for the organic analog aerosol-water mixture is unknown, it has been shown that for subsaturated regimes of RH < 95%, hygroscopic growth is “almost insensitive to σ s/a” (Wex et al., 2008), based on fundamental properties of the Kohler equation. For a given f(RH) measurement, the RH is controlled, D dry is known, D is calculated as described in Sections 3.1 and 3.2, and the temperature at which all experiments were conducted was 293 K. The only unknown in Eq. 6 is κ, which can then be calculated.
In Table 1, we have listed the calculated κ-values for the early Earth and Titan analog aerosols and compared them to other compounds. Also, Table 1 includes κ calculated for our ammonium sulfate f(RH) measurements, which, like the f(RH) and G f measurements, agree with other experimental studies and theory. For all compounds, the italicized values are calculated from the boldface values in the same row. For the early Earth analog aerosol, we calculate κ = 0.15 ± 0.08 at RH = 80%, and for Titan κ = 0.06 ± 0.01 at RH = 80%. As expected from the growth factor values reported in Section 3.2, the early Earth κ-value is higher than the Titan value. The early Earth analog aerosol has a comparable κ-value to levoglucosan, an oxidized organic that is a product of wood combustion. Interestingly, κ ∼ 0.1 is a typical range of values found for secondary organic aerosol generated in the laboratory that is designed to represent aerosol found in the current-day atmosphere (Prenni et al., 2007; Wex et al., 2009).
3.4. f(RH) and κ of early Earth aerosols over a range of relative humidities
Petters and Kreidenweis (2007) showed experimentally that κ derived from hygroscopic data obtained in a subsaturated regime for many organic and inorganic aerosols are within experimental uncertainties of those derived from κ obtained from a supersaturated regime. However, several other studies have found that κ is not constant in present-day Earth secondary organic aerosol particles with increasing RH in both the subsaturated and supersaturated regimes (Hartz et al., 2005; VanReken et al., 2005; King et al., 2007; Prenni et al., 2007; Duplissy et al., 2008; Engelhart et al., 2008; Wex et al., 2009). We assessed the constancy of κ for our early Earth analog aerosol by measuring κ over a range of subsaturated conditions. Assessing the constancy of κ allowed us to assess the CCN ability of the early Earth analog aerosol in the Archean atmosphere based on our results obtained in subsaturated conditions.
Figure 5 shows the 74 individual f(RH) measurements for various particle sizes from RH = 78–87%. Each f(RH) value at a given diameter is one data point that represents an average of extinction measurements sampled every 10.3 s for 2 min. As expected, f(RH) decreases with increasing particle size and increases at a given size as RH increases. At each of the 74 points in Fig. 5, a κ-value can be calculated. The average κ for the entire range and all early Earth analog diameters is 0.22 ± 0.12.

Optical growth f(RH) of early Earth analog aerosols of various particle diameters over a range of relative humidities (RH = 78–87%).
Figure 6 shows a plot of κ versus RH. From Fig. 6, there appears to be a dependence of κ on RH. The RH dependence of κ is quantitatively assessed by calculating the statistical significance of the correlation with Spearman's rank correlation coefficient ρ. A value of ρ = +0.69 is found, which, for 74 data points, indicates a very significant correlation (>99.9% confidence) between κ and RH. Again, this dependence and the direction of this dependence—increasing κ with RH—appears to be common in the literature. A study on three types of secondary organic aerosol by Wex et al. (2009) showed an order of magnitude increase in κ from RH = 80% to RH = 99.6%. Also, the error in the average κ-value for the early Earth analog aerosol over the humidified range is larger than the error obtained at RH = 80%. This difference can be explained by this dependence of κ on RH. Given the trend in increasing κ with RH and similar results in secondary organic aerosol in other studies, our reported early Earth κ-value likely represents an underestimate of the true CCN ability of this early Earth analog aerosol. Additionally, the statistical significance of the correlation between particle diameter and RH was tested. The correlation between these two variables was found to be statistically insignificant (<95% confidence) with a ρ = −0.10 for 74 data points.

κ as a function of RH for early Earth analog aerosols. A positive and statistically significant (>99.9% confidence) correlation between κ and RH is found using a Spearman's rank correlation coefficient analysis.
4. Implications
4.1. Direct effect of humidified aerosol
We estimated the radiative impact of the humidified organic aerosol by using the calculated hygroscopicity parameter κ to compare the globally averaged direct aerosol radiative forcing of the humidified aerosol ΔF
hum to that of the dry aerosol ΔF
dry (both in units of W m−2). The globally averaged solar direct radiative forcing ΔF of an absorbing aerosol layer is proportional to
where a is the albedo of the underlying surface, β is the fraction of the radiation scattered into the upper hemisphere, and τ
sca and τ
abs are the optical depths of the aerosol layer due to scattering and absorption, respectively (Chylek and Wong, 1995). To calculate the absolute magnitude of ΔF, a prefactor including the solar constant, transmittance of the atmosphere above the aerosol layer, and fraction of the sky covered by clouds would need to be included (Chylek and Wong, 1995). Because all these values in the prefactor are not well constrained and, for the same sky conditions, do not explicitly change for aerosol of different optical properties (corresponding to aerosol exposed to different RHs), we did not include them in our calculations. To compare ΔF
hum and ΔF
dry for two otherwise identical aerosol layers of particle concentration C and thickness l, we took the ratio of the two and substituted τ = σCl, where σ is the scattering or absorption cross section of the particle (cm2 particle−1), as indicated by subscript:
For the calculations that follow, a is assumed to be 0.22, the planetary albedo estimated for an atmosphere of 10−3 bar CO2 and 0.1% CH4 (Haqq-Misra et al., 2008). The upscattering fraction β is approximated as β = (1 – g/2)/2 (Sagan and Pollack, 1967; Chylek and Wong, 1995), where g is the asymmetry parameter. For the size range of particles considered, β corresponds to the upscattered fraction for an average 60° solar zenith angle (Wiscombe and Grams, 1976). The optical cross sections σ sca and σ abs are functions of RH and therefore of κ, which dictates the humidified aerosol's size and volume-weighted refractive index. The humidified aerosol size (or G f) is determined by inputting our measured κ-value into Eq. 6. Once G f is calculated, the volume-weighted refractive index of the humidified aerosol can be input into a Mie code assuming spherical particles, in which the optical cross sections and g (and therefore β) are determined as a function of RH.
Figure 7 (a) shows the term proportional to ΔF (given in Eq. 7) over a range of particle diameters and RHs. This calculation assumes the hygroscopicity parameter measured in this work and the optical properties of the dry early Earth analog aerosol measured by Hasenkopf et al. (2010). For dry (RH = 0%) particles in the size range of 115 nm > D > 665 nm, a net cooling effect is induced by the aerosol. Increasing the RH of the particles in this size range further increases their cooling effect. Because of the significant imaginary refractive index of the early Earth analog aerosols, at dry diameters larger than 665 nm, a net warming effect is induced by the aerosol. This heating effect is mitigated as particles are exposed to higher RHs.

(
Figure 7(b) shows the ratio, RΔF , of the ΔF term of humidified haze aerosol to the ΔF term of dry haze aerosol as a function of RH for particles of diameters 200–500 nm. Even at an RH as low as 50%, substantial cooling induced by the humidified aerosol is observed relative to the dry. This increased cooling effect is increased at higher RHs. In the particle diameter range of 200–500 nm, cooling can be increased by a factor of 1.5–3.5 at RH = 90%.
A similar calculation for fractal aggregates would be very interesting yet much more complicated because the amount of exposed surface area would be dependent upon the contact angle between connecting monomers. Also, it is unclear at what RH collapse into a spherical droplet would occur; however, this would be interesting to determine because fractal aggregate collapse into a sphere at a given RH might have significant consequences for the effectiveness of a fractal aggregate UV shield. Regardless, it is clear that the degree of hygroscopicity exhibited by the early Earth analog aerosols indicates that the humidified aerosol, spherical, or fractal could have a significant cooling or decreased heating effect, depending on the size of the aerosol and the specific RH.
4.2. Clouds on early Earth
Photochemical models of haze formation indicate that aerosol present at cloud-relevant altitudes (20 km and below) would be large in size (D ∼ 1 μm) and sparse in concentration (0.1–10 particles cm−3). The number concentration is inversely proportional to the particle radius cubed for a fixed mass so that with a mass production rate similar to today's sulfate aerosols, the concentration of these 1 μm aerosol particles would likely be about 100 times less than current 0.2 μm sulfate CCN (Pavlov et al., 2001a; Wolf and Toon, 2010). It is uncertain whether these large aerosol particles would be spherical or fractal aggregates composed of spherical monomers. Observations on Titan indicate that its organic aerosol particles are fractal aggregates composed of possibly several hundred to several thousand spherical monomer units, each unit roughly 100 nm in diameter (Tomasko et al., 2005, 2008; Bellucci et al., 2009). However, as aerosol particles descend lower in Titan's atmosphere, there is some evidence that the morphology or composition, or both, change possibly due to CH4 condensing on the aerosol particles (Tomasko et al., 2008). On current-day Earth, collapsing and consequent compaction of soot particles in laboratory simulations to a more spherical shape due to ambient water RH has been observed (Zhang et al., 2008). It is unclear what shape the early Earth organic haze aerosols would have been when exposed to a given RH. In the following calculations of CCN ability of the early Earth analog aerosol, we treat the analog aerosols as perfect spheres with 1 μm diameters. Because of the many assumptions the calculations would otherwise require, we neglect any fractal shape. This may underestimate the CCN ability of the aerosol, since theoretical studies have shown that a pair of joined spherical particles will activate at supersaturations much lower than the particles would separately due to the negative curvature of the space at which they are joined and subsequent enhanced capillary action (Crouzet and Marlow, 1995).
The critical supersaturation, s c, defined as the RH above 100% for a given dry aerosol diameter at which an aerosol will grow rapidly and be activated into a cloud droplet, can be calculated by finding the maximum of S in Eq. 6 for a given κ-value and dry diameter. Figure 8 shows a plot of s c versus dry diameter for a range of constant κ-values, including the early Earth and Titan analog aerosols from this study and other reference materials. According to Fig. 8, a 1 μm early Earth analog particle would activate into a cloud droplet well below even a modest supersaturation of 0.1% (equivalently, RH = 100.1%). In fact, even a relatively small early Earth analog aerosol 200 nm in diameter would activate at such a low supersaturation. For comparison, the CCN concentration subset of a population of aerosols is typically defined as what fraction of aerosol is CCN active at s c = 1% (Seinfeld and Pandis, 1998). Figure 8 shows that the early Earth analog aerosol would be CCN active at reasonable—and, in fact, low—supersaturations. Because these early Earth analog aerosols are good CCN and are predicted in models to have been present in low concentrations on Archean Earth, they would grow rapidly due to lack of competition for water vapor, fall out as precipitation, and therefore have a short lifetime anywhere critical supersaturation is reached. This “rain without clouds” scenario is similar to one on Titan, which was first suggested by Toon et al. (1988) and reported in observations on Titan (Tokano et al., 2006). On Titan, the optically thin clouds are thought to perhaps provide a constant source of liquid methane drizzle to the surface, which has a significant effect on the surface geology (Tokano et al., 2006). The degree to which such a constant drizzle would have been present on early Earth is unclear and would need to be resolved with modeling.

Critical supersaturation s c as a function of dry diameter for the early Earth and Titan analog aerosols and other organic and inorganic materials. The top axis shows the corresponding scale for the hygroscopicity parameter κ. The values of κ for the early Earth and Titan analog aerosols are inferred from f(RH) measurements in a subsaturated regime.
If there were indeed a dearth of other types of aerosols on most regions of early Earth, the organic haze may have been a main source of CCN. The degree to which the clouds formed from organic haze aerosol would have been optically thin can be estimated by a simple analytical expression from Schwartz and Slingo (1996), which assumes the two-stream approximation and relates the albedo A
c of a non-absorbing, horizontally homogenous cloud optical depth τ
c:
Where g
c is the asymmetry parameter of the cloud droplets and is approximately 0.85 for particles with radii much smaller than the incident light, a reasonable assumption for cloud droplets. The optical depth τ
c can, in turn, be estimated by an expression from Schwartz and Slingo (1996) that is a function of cloud droplet concentration N
c:
where z is the physical cloud thickness and L is the liquid water volume fraction. By combining Eqs. 9 and 10, the albedo of clouds can be calculated with different droplet concentrations but identical cloud thicknesses and liquid water volume fractions. First, using Eq. 10 and assuming a typical liquid water volume fraction of a cloud, 0.3 cm3 m−3, we calculate that a cloud with 10 droplets cm−3, perhaps typical of an Archean cloud, has only a third of the optical depth of a cloud with 500 droplets cm−3, more typical of a current-day cloud. By extending this calculation to cloud albedo using Eq. 9, the albedo of the early Earth cloud is 33–80% that of a current-day cloud over a range of cloud thicknesses, 50–1500 m. The much lower albedo of these Archean clouds could substantially warm the planet. However, this calculation relies on the Wolf and Toon (2010) and Pavlov et al. (2001a) aerosol models for particle concentrations. Both models must make assumptions about wet deposition of aerosols whose hygroscopic properties have not been previously measured. Wolf and Toon (2010) used predictions of cloud formation and rainfall from the model to simulate removal. However, these predictions are based on current Earth clouds, which may differ from those on ancient Earth. For a given temperature and ocean coverage, the rainfall rate must be about the same as on present Earth because the evaporation rate only depends on temperature and water availability. However, the rainfall does not need to be distributed in the same way spatially. Pavlov et al. (2001a) prescribed a lifetime of 5 days to organic haze aerosol, which is the same as that of highly soluble gases and that previously assigned to sulfur haze aerosols (Kasting et al., 1989). Such a short lifetime may result in an overestimate of the solubility of an organic haze, underestimate atmospheric particle concentrations, and reduce the warming effect due to optically thin clouds. To resolve how the CCN ability of these aerosols affect climate, the hygroscopicity of early Earth analog aerosol would need to be included in early Earth climate models.
5. Conclusion
We measured the hygroscopicity parameter κ = 0.22 ± 0.12 for early Earth analog aerosol, which enables one to calculate the direct radiative effects at any RH and their indirect radiative effects based on their CCN ability. Figure 9 shows a schematic that summarizes the implications for the early Earth climate from this work. The degree of hygroscopicity indicates that the humidified aerosol, compared to the same aerosol that is dry, would likely have a significant cooling or decreased heating effect, depending on the size of the aerosol and the specific RH. For instance, at an RH of 85%, the optical growth of the early Earth analog is 2, which means the optical depth of the aerosol layer is doubled. Thus, the growth of particles due to humidification substantially increases the antigreenhouse effect of the haze that is located in regions with humidities greater than about 50%.

Summary of the implications of a humidified early Earth haze exposed to humidification and supersaturation for early Earth climate based on the hygroscopicity of the early Earth analog aerosol. The hygroscopic nature of the early Earth analog aerosol from this study indicates that the humidified aerosol, compared to the same aerosol that is dry, would likely have a significant cooling or decreased heating effect, depending on the size of the aerosol and the specific RH. For instance, at an RH of 85%, the f(RH) of the early Earth analogue is 2, doubling the optical depth (OD) of an aerosol layer. Assuming a typical liquid water volume fraction of a cloud and using Eq. 10, a typical Archean cloud would have only a third of the OD of a cloud with more current-day cloud droplet concentration. Extending this calculation to cloud albedo using Eq. 9, we find that, over a range of cloud thicknesses, 50–1500 m, the albedo of the early Earth cloud is 33–80% that of a current-day cloud. By including the hygroscopic nature of the aerosol and therefore higher rainout rate in aerosol models, the albedo could be even more substantially reduced. The much lower albedo of a typical Archean cloud could substantially warm the planet.
Because of the biogenic nature of the haze, Domagal-Goldman et al. (2008) and Haqq-Misra et al. (2008) suggested that, regardless of the optical properties of the haze, the degree of cooling it caused on early Earth would have been regulated by a corresponding decrease in CH4 production, which in turn would have reduced haze and restored surface temperatures. However, this negative feedback cycle could have been at least partially offset by a positive feedback cycle involving a UV shield that may have been provided by the haze. If CH4 production was reduced because the antigreenhouse effect of haze was too large, then a thinner haze would have resulted. The thinner haze, while inducing warming, could also have potentially destroyed the organisms' UV shield. The optical properties of the aerosol—both dry and humidified—would have dictated the balance between these two feedback cycles on early Earth.
In contrast, warming relative to current-day clouds is expected due to the indirect effect of the organic haze. The value of the hygroscopicity parameter κ measured indicates that the early Earth aerosol could have activated into cloud droplets at reasonable—in fact, low—supersaturations. Photochemical models of the haze predict very low concentrations, compared to current-day Earth standards, of micron-sized organic aerosol located 20 km and lower in altitude. If these aerosol particles were the only particles present, they could activate and grow unimpeded by water-vapor competition of other aerosols, further lowering the particle concentrations than current Archean aerosol models predict. The lifetime of the droplets in these optically thin clouds would be short. If globally present on early Earth, these optically thin clouds would lower the planetary albedo in comparison to today's value and therefore warm the planet. This “rain without clouds” scenario would mean a hydrological cycle very different from the present one.
Footnotes
Acknowledgments
This material is based on work supported by NASA Exobiology Grant NNXOAAV55G issued through the Office of Space Science. C.A.H. was supported with a National Science Foundation Graduate Research Fellowship. M.A.F. acknowledges support from the NOAA Climate and Global Change Postdoctoral Fellowship Program administered by the University Corporation for Atmospheric Research. M.R.B. was supported by an EPA-STAR fellowship.
Disclosure Statement
No competing financial interests exist.
Abbreviations
CCN, cloud condensation nuclei; CRD-AES, cavity ringdown aerosol extinction spectrometer; DMA, differential mobility analyzer; e-AIM, extended Aerosol Inorganic Model; RH, relative humidity.
