Natural aerosols explain seasonal and spatial patterns of Southern Ocean cloud albedo

See allHide authors and affiliations

Science Advances  17 Jul 2015:
Vol. 1, no. 6, e1500157
DOI: 10.1126/sciadv.1500157


Atmospheric aerosols, suspended solid and liquid particles, act as nucleation sites for cloud drop formation, affecting clouds and cloud properties—ultimately influencing the cloud dynamics, lifetime, water path, and areal extent that determine the reflectivity (albedo) of clouds. The concentration Nd of droplets in clouds that influences planetary albedo is sensitive to the availability of aerosol particles on which the droplets form. Natural aerosol concentrations affect not only cloud properties themselves but also modulate the sensitivity of clouds to changes in anthropogenic aerosols. It is shown that modeled natural aerosols, principally marine biogenic primary and secondary aerosol sources, explain more than half of the spatiotemporal variability in satellite-observed Nd. Enhanced Nd is spatially correlated with regions of high chlorophyll a, and the spatiotemporal variability in Nd is found to be driven primarily by high concentrations of sulfate aerosol at lower Southern Ocean latitudes (35o to 45oS) and by organic matter in sea spray aerosol at higher latitudes (45o to 55oS). Biogenic sources are estimated to increase the summertime mean reflected solar radiation in excess of 10 W m–2 over parts of the Southern Ocean, which is comparable to the annual mean increases expected from anthropogenic aerosols over heavily polluted regions of the Northern Hemisphere.

  • Atmospheric science
  • clouds
  • aerosols
  • weather
  • Climate
  • Southern Ocean
  • sea spray
  • albedo


The Southern Ocean (SO) is an expansive and dynamic ocean with rich ecosystems remote from most human influences. It is also the cloudiest region on Earth. These clouds influence the atmospheric and oceanic circulation of the entire Southern Hemisphere and beyond (1), and may help determine the Earth’s climate sensitivity (2). Its remoteness from anthropogenic and natural continental aerosol sources makes the SO a unique natural laboratory for our understanding of aerosol-cloud interactions. Aerosols influence clouds by acting as the cloud condensation nuclei (CCN) on which cloud droplets form, and the resulting concentration Nd of cloud droplets influences the amount of sunlight reflected by clouds (3). Aerosol processes remain a poorly understood influence on clouds (4). Processes regulating the concentration of naturally occurring aerosols, in particular, remain a major source of uncertainty, limiting our ability to quantify the magnitude of the human impact on climate from aerosols (5, 6).

The rich biological activity in the SO produces a range of biogenic aerosols and aerosol precursor gases (7). Marine biogenic emissions influence atmospheric aerosols in this region, both through primary emissions of organic matter in sea spray aerosol (SSA) (8) and through secondary aerosol formation processes, that is, the condensation of volatile sulfurous and organic compounds (79). The effect of marine biota on clouds and climate through biogenic sulfate aerosol has been the subject of intense research for several decades (10). In recent years, primary organic aerosol species have received considerable attention as well (8). Chlorophyll a (Chl-a), a proxy for the phytoplankton biomass believed to be a source for aerosol and aerosol precursors, has been found to be negatively correlated with the radius of cloud droplets over the seasonal cycle in selected biologically active regions of the SO, indicating that aerosol originating from marine microbial communities may be associated with enhanced concentrations of CCN (11, 12). Although previous investigations did much to probe the coupling between biological activity and clouds (8, 1014), they were limited in three ways. First, previous studies did not correct for seasonal biases in the remote sensing of cloud properties at high latitudes, which were unknown at the time. This conflated true seasonal variations in Nd due to CCN number with seasonal retrieval biases due to solar angle. Second, they did not examine the distributions of both primary aerosols and secondary sulfate aerosols from marine biota and their respective roles as mediators of the Chl-a–Nd correlation. Finally, most previous studies examined limited regions within the SO where strong seasonal cycles in Chl-a may correlate with underlying seasonal cycles in meteorological variables, confounding the detection of an unambiguous signal of phytoplankton effects on clouds.

The Moderate Resolution Imaging Spectroradiometer (MODIS) instrument on the NASA Terra satellite offers remotely sensed estimates of Nd with global coverage. In this work, we use a new MODIS Nd data set that eliminates certain retrieval biases (15), allowing for the first time the comprehensive study of the relationship between Nd and modeled marine biogenic aerosols over the entire SO. Biases due to seasonal-latitudinal variations in the solar zenith angle (SZA) affect the visible to near-infrared satellite retrievals of cloud optical thickness and droplet effective radius used by MODIS to estimate Nd. These biases are mitigated by removing all retrievals with SZA >65° (15).


With this new data set, we find that spring-summer-autumn (September to April) mean that Nd values are elevated in the regions of the SO with the most active ocean biology as indicated by Chl-a (Fig. 1), demonstrating that ocean biology drives variations in Nd across the entire SO.

Fig. 1 Elevated mean September to April cloud droplet concentrations over the SO are associated with regions of high Chl-a (indicating the presence of phytoplankton biomass).

Spring-summer-autumn (September to April) mean values of Nd are shown as a function of SeaWiFS Chl-a (for 5° latitude × 15° longitude boxes from 35° to 55°S). Three comparisons are shown: Nd from 2007 only using retrievals SZA <65° versus Chl-a from 2007 (red); Nd from 2007 (only retrievals with SZA <65°) versus a climatology of Chl-a from 2001 to 2009 (black); and Nd from each year from 2001 to 2009 (all retrievals) versus Chl-a from each year from 2001 to 2009 (brown). The Spearman rank correlations are 0.76, 0.63, and 0.57 ± 0.04, respectively, where the mean and SD of the latter is calculated from the individual correlations from each year from 2001 to 2009.

Correlations between Nd and candidate aerosols

It is interesting that a significant correlation is found between regions of high Chl-a and enhanced Nd; however, a more physical understanding of what drives the spatiotemporal variability in the SO necessitates investigation of the SO CCN amount and provenance. Unfortunately, remote sensing alone cannot determine the detailed chemical composition of aerosol. In particular, it cannot distinguish between sea salt (SS), sulfate, and organic aerosol because they all primarily scatter rather than absorb incoming radiation. To determine which aerosol species are responsible for driving the observed Chl-a–Nd relationship over the SO, we combine Nd observations with candidate aerosol products derived from global modeling to develop a regression model to predict Nd over the SO. Sulfate and SS surface concentrations are taken from the multimodel median year 2000 global model simulations performed for phase 1 of the AeroCom project (16). Aerosol modeling is extremely complex and subject to a great deal of uncertainty and, consequently, intermodel disagreement (17). Because of this diversity, the use of the AeroCom multimodel median data set is likely to provide a more robust representation of the distribution of an aerosol species than any single model. The non-SS aerosol candidates selected here, marine organic matter and SO4, have been chosen on the basis of observational evidence for their efficacy as CCN. Sulfate aerosol is found to serve efficiently as a CCN (18, 19), and in situ measurements in marine air have identified a contribution of organic aerosols to CCN (20). In the SO region, dust may be important as a source of ice nuclei (21, 22) and oceanic micronutrients such as iron (23, 24). However, dust sources in the Southern Hemisphere are minimal, and the mean contribution to the total aerosol number and CCN in the SO is small, particularly when compared to sea spray and sulfate aerosol. This has been established by satellite and in situ observations and is reflected in global chemistry-climate models (17, 2531). We will now discuss the inclusion of organic enrichment of SSA as a predictor variable.

The mechanisms by which phytoplankton affect organic aerosol in the atmosphere are a subject of ongoing research and debate, and there is no consensus yet on the effect of marine phytoplankton and organic matter on the number, size, mixing state, and detailed organic chemistry of emitted particles. Model simulations of atmospheric transport and removal of aerosols, and of aerosol-cloud interactions, are extremely sensitive to the details of such assumptions in models, and differing choices may result in strong changes of the magnitude or even the sign of effects on CCN and Nd (12). Furthermore, although strong evidence supports the existence of a primary sea spray source of marine organic aerosol, the relative importance of condensational growth by secondary aerosol formation from marine precursor compounds is not well understood. However, a greater role for primary organics is suggested by evidence such as (i) the chemical similarity of aerosol generated at sea to ambient marine aerosol (3234); (ii) single-particle measurements of marine organic particles, which show that these particles frequently contain solid fragments and agglomerates that cannot be explained by condensation (3537); (iii) the correlation of wind speed with organic mass fraction (OMF) observed at sea in remote locations (38); and (iv) the observation that until now, most studies of which we are aware are able to demonstrate only a negligible to minor contribution to organic mass (OM) in naturally occurring marine aerosol from condensation of organic precursor gases [for example, oxalate (39, 40) or isoprene (41)] or from atmospheric oxidation processes (39, 42). Here, we use the OMF in emitted submicron SSA as a proxy for the effect of marine organic matter on primary SSA. OMF is calculated using OCEANFILMS, a physically based framework that models the enrichment of emitted primary sea spray particles with organic matter as being driven by adsorption of surface-active organics at the air-water interfaces of ocean bubbles, using fields from an ocean biogeochemical model to drive the analysis (43, 44). Observational evaluation of the OCEANFILMS model shows skill in predicting OMF (43) despite inherent limitations in our ability to fully represent the complexity of ocean biogeochemistry (45). Although the OMF of emitted particles is not a direct estimate of aerosol number concentration, it has the virtue of not requiring additional assumptions to be made about the emissions process, for which observational support is not yet available. In particular, the effect on emitted number is a critical parameter needed by models to simulate the impact of marine organics on CCN number in models (46). We are unaware of firm observational evidence from field or laboratory studies that can inform models as to whether the total number and mass of SSA emitted increase, decrease, or remain the same in the presence of marine organics, although limited evidence suggests that total emissions may increase under some circumstances (47).

The OCEANFILMS parameterization is the first mechanistic model of the enrichment of organic matter in SSA appropriate for use in a global atmospheric model. An important feature of this model that distinguishes it from previous parameterizations of OMF is that it allows for a contribution to organic enrichment of SSA by classes of molecules whose spatial distribution is decoupled from that of Chl-a, namely, the polysaccharide- and protein-like classes of semi-labile organic macromolecules. This is particularly important in the nutrient-poor regions of the ocean, where Chl-a and primary productivity are low, but longer-lived marine organic matter is imported by ocean currents from regions of high primary production. Observations show that for a given Chl-a concentration, marine aerosol organic enrichment in oligotrophic regions is higher than would be suggested by extrapolation from observations in mid-latitude bloom regions such as those in the North Atlantic (32, 4749).

The candidate aerosol products used as predictors for Nd are near-surface sulfate (SO4) concentration, near-surface SS concentration, and OMF of emitted SSA. Additional candidate predictors considered are the surface wind speed, sea surface temperature (SST), ocean dimethyl sulfide (DMS) concentration (50), and Chl-a.

The relationships between Nd and the predictor variables were determined by multiple linear regression on monthly mean data averaged over 15° longitude × 5° latitude bins for those months where Nd retrievals are available. By using geographically binned monthly mean data, we improve signal-to-noise ratio and mitigate potential sampling biases related to the greater availability of usable retrievals. The uncertainty range in the variability in Nd explained by each set of predictor variables (R2) was obtained using the 10-fold cross-validation (51) procedure repeated 10 times. The variability explained by each predictor is reported in Table 1. Individually, SO4 and OMF are the two variables that explain the most seasonal and geographical variability in Nd, explaining 48 ± 7% and 35 ± 10%, respectively. These relationships are shown in Fig. 1. An increase in Nd with OMF (Fig. 1A) is found for high, intermediate, and low SO4. Similarly, for a given concentration of SO4, Nd generally increases with OMF (Fig. 2B). The SO4 concentration and the OMF are only moderately correlated with each other (R2 = 0.39), indicating the appropriateness of including both as independent predictors here. Only 9 ± 4% of the Nd variance is explained by SS, indicating that sea salt mass loading is not an important contributor to the spatial and seasonal variability in Nd over the SO during September through April. The variances explained in the spatial structure of the September to April Nd by SO4 and OMF are 48 ± 24% and 33 ± 22%, respectively.

Table 1 The variability in Nd explained by each predictor.

R2 is given in percent for both unweighted calculations and a robust calculation using bi-square weights, using a bias-corrected Nd data set. Spatial R2 values are calculated by applying the same regression procedure and using bi-square weights, but using time-averaged values of bias-corrected Nd for the period September through April. Additional notes following asterisks give the R2 values calculated using a 9-year monthly climatology of MODIS Nd including retrievals at all SZAs. Uncertainties are the SDs of estimates of R2 from 10-fold cross-validation, repeated 10 times.

View this table:
Fig. 2 Monthly mean Nd correlates with OMF and SO4 concentration.

Data points are given as gray dots. Box and whisker plots show binned means, boxes show interquartile range, whiskers show twice interquartile range, and outliers are shown as colored circles. Left: Nd as a function of OMF. Box and whisker plots are shown for three levels of SO4 concentration. Higher SO4 surface concentrations (orange) are associated with higher Nd for the same OMF. Right: Nd as a function of SO4 surface concentration. Box and whisker plots are shown for three different bins of OMF, where higher OMF (orange) is associated with higher Nd for a given SO4 concentration. All data are averaged monthly and over 5° × 15° boxes. R2 values are shown for linear regressions between Nd and the respective predictor variable.

The multiple linear fit to both SO4 and OMF explains 53 ± 6% of the variance in Nd, and the inclusion of additional aerosol-related predictors, including sea salt, did not significantly increase the explained variance. The total spatial variance in the September to April mean explained by OMF + SO4 is 53 ± 22%. We find that monthly mean Nd as predicted by this regression is within 20% of observed Nd 74% of the time. The regression model of Nd isEmbedded Image(1)where Nd has units of cm−3, the surface SO4 concentration has units of μg m−3, and OMF is dimensionless. The constant term can be considered as the expected Nd over an ocean devoid of biological activity and is similar to aircraft measurements of Nd made off the coast of Tasmania (52) during winter, when biogenic aerosol sources are minimal.

Although this analysis cannot explicitly determine the composition of the CCN that determines background values of Nd (≈36 cm−3), as estimated from the constant term in Eq. 1, it is likely that this background is largely attributable to SS. A detailed model study of CCN sources over the SO (7) found that SS contributed significantly, but the simulated number of CCN from SS was almost the same in summer and winter. As a result, SS explains little of the variance in Nd, but likely contributes to the constant term in the multiple linear regression. Under the assumption that the constant term is entirely attributable to SS, the summertime contribution of SS to Nd is about 55% and the wintertime contribution is near 80% in the 35° to 55°S region. As shown in Fig. 3, this is highly consistent with a previous model evaluation of the relative importance of SO4 and SS as CCN in the SO, with a January contribution from SS of 48% at 30° to 45°S, 58% at 45° to 60°S, and a wintertime contribution of >80% at 30° to 45°S and >90% farther south (7).

Fig. 3 Comparison of the fractional contributions to CCN from SS as estimated from this study and as modeled in the global aerosol model GLOMAP (7).

The percentage contribution of SS to sulfate and SS CCN is estimated in this study as the constant term in Eq. 1 divided by the sum of the sulfate and constant terms from Eq. 1. The winter and summertime contributions of SS to CCN modeled by GLOMAP are shown in the figure with triangles, and the dashed connecting lines are provided as guides to the eye. The SS contributions from this study (oceans only) are shown aggregated into summer, fall, winter, and spring seasons and into 35° to 45°S and 45° to 55°S.

Because the modeled sulfate concentration, the modeled OMF, and the satellite data do not represent the same years, our analysis cannot capture interannual variability in clouds or ocean biota. As seen in the observed correlation between Chl-a and Nd, the correlation decreases significantly when a climatology of Chl-a is used instead of the Chl-a observed in the same year as the Nd data (Fig. 1), suggesting that matching of the exact location, timing, and strength of phytoplankton blooms can improve the correlation between the presence of oceanic phytoplankton and cloud properties. This suggests that the predictive capability of the regression model in Eq. 1 might also be further improved by better representation of the interannual variability of simulated ocean biological activity and that the estimate of the explained variance due to OMF and SO4 offered in this work is conservative.

The spatial structures of observed and predicted September to April mean Nd agree well (Fig. 4, A and B). Contributions to Nd from SO4 (Fig. 4C) and OMF (Fig. 4D) determined from the regression model indicate that the impact of sulfate on summertime Nd is greater at lower SO latitudes (35o to 45oS), whereas the OMF contribution is more heterogeneous and spread over a broader range of latitudes with the major patches confined to higher latitudes (45° to 55°S). The seasonal cycle of Nd (Fig. 4F) shows a strong summertime peak in agreement with in situ measurements of CCN in the SO (53). Comparison of observed Nd (Fig. 4A), which shows a mean modeled Nd of 59 cm−3, with the background for a biologically inactive ocean Nd ≈ 36 cm−3, suggests that ocean biology increases Nd by about 60% on average during September to April (spring-summer-autumn).

Fig. 4 The model captures significant spatial and seasonal structure in the SO Nd.

(A) One year of Nd observations averaged from September through April (when observations are available for all months at all latitudes). (B) Predicted Nd from the regression model. (C and D) Contributions from OMF and SO4 (note the different color scale). The means and spatial SDs are shown above each panel. The correlation between modeled and observed September-April mean Nd is r = 0.71, and the correlation between the September-April mean OMF and SO4 contributions to Nd is r = 0.41. (E) Mean ocean Chl-a observed by SeaWiFS. (F) Time series of the observed Nd, modeled Nd, and the contributions of OMF, SO4, their sum, and the constant term to Nd. Time series data are monthly averages over all locations where Nd is observed.

The highest SO4 concentrations occur primarily at lower latitudes, whereas high OMF is found at all latitudes. Sunlight is required to oxidize DMS, and continental sources of natural sulfate aerosol are also stronger at lower latitudes (7). The sulfate sources affecting the northern SO (35° to 45°S) potentially have a significant volcanic component (54), in addition to the DMS sources of sulfate which dominate at higher latitudes (5456). Anthropogenic emissions contribute less than a third of the SO4 concentration on the annual mean over the study region, with a peak contribution to SO4 of 35% at 35°S rapidly declining to 10% at 45°S; this is consistent with estimates of the DMS contributed SO4 in this region provided by dedicated aerosol tagging studies (56). The Nd pattern is consistent with attributing much of the variability in Nd in the 35° to 45°S latitude band to sulfate (Fig. 4D). Organics in SSA contribute strongly in regions of high Chl-a (such as off the eastern coast of South America).

Increases in reflected shortwave radiation

These results imply that natural aerosols have the potential to significantly affect cloud droplet concentration over the SO. Observational estimates have shown that cloud microphysical processes can significantly affect the strength of the enhancement in cloud albedo with warming in the SO (57). Thus, accurate representation of natural, non-SS aerosols is likely to be important not only for accurate representation of the SO albedo in the current climate by global climate models (GCMs), which is generally too low (2), but also to help determine the strength of the optical depth feedback as the climate warms. Here, we quantify the effect of sulfate and primary organic aerosols on reflected shortwave (RSW) radiation, estimated in two ways: (i) from observed cloud properties and radiation (11), and (ii) using an offline radiative transfer calculation performed in Community Atmosphere Model version 5.3.38 (CAM5) using Portable Offline Radiative Transfer (PORT) (58). These calculations consider only the impact of changing Nd on RSW (the cloud albedo effect). Other aerosol effects, such as direct forcing, the cloud lifetime effect, and the semidirect effect, are not considered. We find that the increase in RSW from biogenic aerosol effects on Nd (Fig. 5A) exceeds 10 W m−2 during summer. Although the latitude range of the observed Nd data is limited by the range of SZAs at which a measurement can be made, the OMF and SO4 data are modeled globally, and thus we may use them to extend the geographic range of our estimate of increased RSW to 70°S. This is especially interesting because of the increase in phytoplankton, and consequently OMF and RSW, at high latitudes during the summer as insolation increases and sea ice recedes (Fig. 5B). Although the sulfate contribution to the summertime increase in RSW is larger than that from OMF over the SO as a whole, both are found to be important. Another interesting finding is that enhanced RSW from OMF is often stronger than the contribution from SO4 at high latitudes (Figs. 4 and 5). This suggests that detailed in situ observations of aerosol physical and chemical properties in combination with measurements of aerosol precursors and cloud properties may be able to distinguish the contributions from these two important species.

Fig. 5 OMF and SO4 significantly affect SO cloud albedo.

(A to C) Contours show the zonal mean change in top-of-atmosphere (TOA) net shortwave radiation (RSW) due to modification of Nd by (A) SO4 and OMF, (B) OMF alone, and (C) SO4 alone. (D) The annual mean contribution from OMF, SO4, and the combination of OMF and SO4. In (A) to (C), filled contours show the analytical estimate of the change in RSW based on observed cloud properties, and dashed white lines show the change in RSW calculated by an offline radiative transfer model using simulated cloud properties. Contour intervals are 1 W m−2 for both. Hatched white lines show the observed zonal mean sea ice extent (in percent). (D) Annual mean changes in RSW from the analytical estimate (solid lines) and the offline radiative transfer model (crosses).

Previous studies have estimated the effect on annual mean RSW of sulfates from DMS emissions at 16 W m−2 over the SO in summer (59) and the effects of the observed seasonal cycle of effective radius, hypothesized to be due to biogenic aerosols, at 6 to 8 W m−2 during summer in the 40° to 60°S region (60). Here, the annual mean contribution to reflected shortwave (SW) is estimated to be between 4 and 6 W m−2 at 40°S with decreasing contributions poleward, and the summertime contribution is 6 to 10 W m−2 (Fig. 5D). This is comparable to the annual mean radiative forcing from aerosol-cloud interaction downwind of highly polluted regions as simulated in GCMs (55, 61). It is also interesting to note that the bias in absorbed SW in the SO between 30° and 50°S in CMIP3 (2) is similar to the effect of SO4 and OMF on RSW as estimated here. Given the significant enhancement in RSW estimated here, a detectable signature in cloud albedo should be evident from remote sensing. Recent studies have revealed a pattern of cloud albedo enhancement that is strikingly similar to the pattern of enhanced Nd described by this study, indicating the robustness of our results (62).

Possible pathways of enriched sea spray influence on Nd

While in situ studies have shown that organic aerosols can affect CCN (20) and our results indicate that organic matter in SSA over regions of high Chl-a affects Nd, there is no broad consensus on the mechanisms that generate these relationships. Some studies suggest that the presence of marine surfactants enhances the number of particles emitted in the range of diameters (roughly 50 to 200 nm) that are activated to form the majority of cloud droplets in the marine boundary layer (47, 63). Aerosol chemistry is likely to be altered by biological activity, both through primary aerosol emissions and through condensation of biogenic volatile organic compounds, and experimental studies show that increases in the organic fraction in sea spray particles with diameters less than 180 nm correlated with increased biological activity (63, 64). Another possibility is that secondary organic material condenses onto the primary spray particles, growing the smallest of them into CCN active sizes (11). Other observations suggest that although SSAs containing high fractions of organic material exhibit low hygroscopic growth at subsaturated relative humidities, they may nevertheless serve effectively as CCN under supersaturated conditions, possibly due to the surface activity of the organic matter (65). Surface-active organic material can reduce the aerosol surface tension, leading to a lowering of the critical activation radius and enhancing activation. Even comparatively weakly surface-active organic substances can facilitate droplet activation when they make up a large fraction of aerosol mass (66), as has been observed downwind of strong phytoplankton blooms (67, 68). A fifth possibility is that the presence of a condensed film kinetically inhibits aerosol water uptake and growth, thereby delaying droplet formation (69). In the absence of a sufficient concentration of more efficient CCN, high supersaturations may be reached, ultimately allowing the activation of a larger proportion of the aerosol population. Finally, the mixing state of the aerosol may also evolve under the influence of biological changes to ocean chemistry, affecting CCN number (70).


Here, we present observational data indicating a significant spatial correlation between regions of elevated Chl-a and Nd across the SO, and show that modeled OMF + SO4 explains 53 ± 22% of the spatial variability in observed Nd. Further, a spatial and temporal correlation is shown relating observed Nd to modeled SO4 and SSA organic enrichment, where the SO4 is likely to be significantly biogenic. The regression model captures a large fraction, 53 ± 6%, of the spatial and seasonal variability in Nd, demonstrating that distributions of natural aerosols simulated by present-day climate models have sufficient fidelity to be useful in predicting climatological cloud properties over the SO. Because of the complexity of chemical and physical interactions affecting sea spray production and processing in the atmosphere, experimental and observational studies have been unable to unambiguously determine whether marine organic aerosol derived from phytoplankton is associated with an increase or with a decrease in CCN and Nd in marine clouds. The results of this study will inform future modeling of these processes by providing the first top-down estimate of the sign and magnitude of the impact of marine biogenic aerosol on Nd across the entire SO.

Our results suggest that ocean biology augments the cloud droplet concentration over the biologically active SO by about 60% in the annual mean and likely doubles it during summer when the incident solar radiation is strongest. The enhancement in upwelling SW due to these changes in cloud droplet number concentration is between 4 and 6 Wm−2 annually averaged and exceeds 10 W m−2 in summer. Accurate representation of the biogenic contribution to CCN and Nd over the oceans is therefore necessary for the accurate representation of climate forcing from aerosol-cloud interactions. To represent these processes in a changing climate will require models capable of representing the dynamic response of ocean biology to the warming and acidification of the oceans and to the retreat of sea ice and predicting the subsequent effects on natural aerosol sources.


Calculation of upwelling SW changes due to biogenic aerosol-cloud interactions

Meskhidze and Nenes method. The enhancement in upwelling shortwave radiation (RSW) caused by increases in Nd was estimated using the following equation, consistent with previous studies (11),Embedded Imagewhere ΔF is the change in RSW, Fin is the downwelling SW at cloud top, Ac is the monthly MODIS cloud fraction, Rc is the monthly MODIS cloud albedo, ΔlnNdb is the change in Nd versus the background concentration ΔlnNdb = (Nd − Nb)/Nd, and Nd and Nb are the perturbed and background Nd, respectively.

Monthly, gridded observations of cloud properties, sea ice, and radiation from 2007 were used in the calculation of changes in RSW, excluding land surfaces. Cloud albedo was estimated from the MODIS cloud optical depth using the cloud albedo approximation of Lacis and Hansen (71). Cloud fraction was taken from MODIS, and downwelling SW at cloud top was approximated by the Clouds and Earth’s Radiant Energy Systems (CERES) Energy Balanced and Filled (EBAF) 2.8 downwelling clear-sky SW at the surface. Sea ice distributions were taken from HADISST (Hadley Centre Sea Ice and Sea Surface Temperature) as included in the NOAA (National Oceanic & Atmospheric Administration) Optimum Interpolation (OI) SST, and the cloudy fraction of each grid box was assumed not to contribute to changes in RSW if clouds were over sea ice (72). Maximum overlap between sea ice and MODIS cloud fraction was assumed, resulting in a conservative estimate of the change in RSW in the presence of ice.

The Nd resulting from the presence of natural aerosols was calculated using the AeroCom sulfate concentration and B14 OMF data in combination with the regression detailed in the article (Eq. 1). The impacts of SO4 and OMF + SO4 on the TOA SW forcing were estimated as the change in RSW when Nd was increased from the background term in Eq. 1 to include the SO4 term, OMF term, and their combination.

PORT method. Offline radiative transfer calculations were performed using the PORT tool (58) within the stand-alone CAM, which is the atmosphere component of the Community Earth System Model (CESM) (73). Details can be found online ( Model simulations were performed using the CAM5 physics package (74). The model was compiled in an atmosphere-only configuration with prescribed ocean and ice boundary conditions with 1.9° × 2.5° grid spacing. CAM5 uses the Rapid Radiative Transfer Model for General circulation models (RRTMG) (75) to compute the radiative forcing of aerosols, clouds, and gases.

To obtain samples of the atmospheric state, instantaneous values of atmospheric state variables required for calculation of radiative transfer were recorded during a 2-year simulation (beginning 1 January 2000), of which the first year was discarded, allowing ample time for model spin-up. Atmospheric state variables were sampled every 73 model time steps, which corresponds to 1.5 days + 1 time step. This sampling frequency has been found to produce a representative sample of diurnal and seasonal cycles while reducing the computation cost of the calculation (58).

After the time samples were produced, the cloud drop size distribution of the time samples was then modified to correspond to each of three regression cases: background, background + SO4, and background + SO4 + OMF (described in point 3 below) for the subsequent PORT calculation. All other atmospheric state variables remained unchanged; therefore, changes in radiative forcing reflect changes in the cloud drop number concentration and related changes in cloud drop size distribution only, with no feedbacks onto cloud lifetime, atmospheric dynamics, or other processes.

The modification of the cloud drop size distribution was performed by the following procedure:

1. Calculate Ndinst_cldtop, the cloud drop number concentration at cloud top (in each column), using the value of the in-cloud droplet number concentration that is passed to the radiation code. Cloud top was defined as occurring in the highest model layer in which the stratus cloud liquid water mixing ratio (icwmrst) was greater than 1 × 10−6 kg/kg.Embedded Imagewhere ρ is the atmospheric density.

2. Calculate the climatological monthly mean of the cloud top droplet concentration (Ndinst,cldtop) over the entire 5-year time series:Embedded Image

3. Calculate Ndreg, the cloud top droplet number concentration as predicted from the regression model (Eq. 1), for each of three cases: (i) FULL: with all terms; (ii) SO4_BKGD: with coefficient of OMF term set to zero; (iii) BKGD_ONLY: with coefficients of SO4 and OMF terms set to zero (constant term only).

The following steps are performed separately for each regression case.

4. Update the cloud drop number concentration in each column by multiplying by the ratio of the cloud top values of Nd from the regression and from the monthly mean:Embedded Image

5. Constrain Ndnew to be in the range [min[Ndinst], max(Ndinst)] on a per-model-level basis:Embedded Image

Similarly, constrain Ndnew so that the grid mean value (as opposed to the in-cloud value) of Nd is within the range of values in the original simulation (on a per-level basis).

6. Calculate new gamma parameters (radμ, radλ) using Ndnew and instantaneous values of air density (RHO) and stratus cloud fraction (AST). This calculation is identical to the calculation in the microphysics (76). Write out the updated values of radμ, radλ to the PORT input files.

Finally, PORT was run using the modified CAM history files to provide inputs to the offline radiative transfer code. This produces a new set of instantaneous history files that include the TOA radiation terms (upward, downward, and net, for SW and longwave radiation).

Estimation of anthropogenic component of SO4 in the remote SO

Although the relative contributions of volcanoes and DMS are difficult to quantify, models indicate that near-surface sulfate over the SO is primarily from natural sources, with a peak biogenic contribution in the summer (5456, 77). Here, SO4 is represented by the present-day surface concentration from the AeroCom A phase 1 multimodel median product. This includes SO4 from volcanic, biogenic, and anthropogenic sources. Because simulated aerosols are not routinely tagged according to their sources in standard atmospheric model simulations, estimation of the volcanic component of the SO4 would require a dedicated simulation (56). The volcanic contribution cannot be determined post hoc for the AeroCom product; however, previous studies have indicated that biogenic sulfate is the source of the majority of SO4 over the SO poleward of 45°S, with increasing volcanic contributions northward of this latitude (55, 56).

AeroCom phase 1 included three experiments: A (models are run with their native, present-day emissions), B (all models use the same, present-day emissions), and PRE (all models use the same, preindustrial emissions). Here, AeroCom A has been used to represent the present-day aerosols. Although the DMS emission inventories used in the AeroCom median are relatively similar, which leads to a lower model diversity in the case of SO4 emissions (78), it incorporates a wider range of aerosol emissions, and thus, its median is hoped to represent the most diversity. On the other hand, this diversity makes it difficult to estimate anthropogenic increases in sulfate because the preindustrial and scenario A experiments do not use consistent emissions, so differences between the two may be either anthropogenic or natural (28).

Scenarios B and PRE have identical natural emissions, and their anthropogenic emissions are scaled with population size, crop production, and wood consumption (28). Thus, differences between scenario B and PRE are purely anthropogenic and provide an approximation of the total present-day anthropogenic aerosol contribution (because preindustrial anthropogenic sources are small). Changes in SO4 column density and surface concentration are shown in figs. S1 and S2, respectively. It can be see that in the eight models within the AeroCom median (all except MPI-HAM), for which SO4 column density data were available, the anthropogenic contribution to loading is relatively minor and ranges between 10 and 20%; however, column density is not the variable examined in this study, and it is expected that near-surface SO4 will be significantly less affected by anthropogenic emissions (77). The four models considered in the AeroCom phase 1 median data set (LSCE, GISS, LOA, and UIO-GCM) for which near-surface SO4 concentration was available in the preindustrial data show that the increase in SO4 from preindustrial to present day is relatively small over the SO (fig. S3). The mean fractional contribution of industrial-era anthropogenic SO4 is at most 35% (at 35°S) and is 10% or less at latitudes south of 45°S (fig. S4), indicating that the majority of sulfate in the remote SO is natural, rather than anthropogenic.

Description of OMB14 (OCEANFILMS) parameterization and contributions of different components of marine OM to total OM

The recently developed OCEANFILMS (43) parameterization is the first mechanistic model of the enrichment of organic matter in SSA appropriate for use in a global atmospheric model. Enrichment is modeled as competitive physical adsorption of several classes of marine organic matter onto bubble surfaces, including a lipid-like class (LIPID) that is strongly associated with fresh primary production, and polysaccharide-like (POLY) and protein-like (PROT) classes that are associated with semi-labile dissolved organic matter with a modeled chemical lifetime of 100 days.

Marine biota and organic matter concentrations are simulated by the Biogeochemical Ecosystems Cycling (BEC) model within the Parallel Ocean Program (POP) (43). The BEC model has recently been optimized and evaluated against an extensive set of marine organic matter measurements (44). The distributions of macromolecular classes used in OCEANFILMS are derived by association with biogeochemical variables as calculated by BEC, including phytoplankton and zooplankton concentrations (LIPID), and semi-labile dissolved organic carbon (DOC) (POLY and PROT). The chemical properties of these classes and their propensity to adsorb to the air-water interface are estimated from the properties of well-studied proxy molecules. A more detailed evaluation of the ocean concentrations of the macromolecular classes is in progress (Ogunro et al., 2015, in revision). In an initial evaluation, the parameterization has been shown to agree reasonably well with observations of the OMF of SSA in clean marine air from both shipboard and coastal clean-marine measurements, as well as with spray generated at sea. The OMF distributions used here are based on the 33rd year of a free-running simulation initialized in the year 2000.

A unique feature of the OCEANFILMS model is its treatment of marine organic matter as a sum of several broad classes of marine macromolecules. The spatial variability in the model’s total OMF is almost entirely determined by two components: a “lipid-like” labile portion that is strongly associated with regions of primary productivity and highly correlated with Chl-a (LIPID), and a semi-labile portion with a model lifetime of 100 days (“polysaccharide-like”/POLY and “protein-like”/PROT classes). The R2 values for regressions of Nd against individual classes (LIPID and POLY) are given in table S1. Breaking down the total organic fraction into these two components reveals differences in the spatial patterns of their contribution to Nd (fig. S5). The labile LIPID group contributes much more in the immediate vicinity of phytoplankton blooms (as indicated by Chl-a) than in oligotrophic regions. The total semi-labiles, POLY + PROT, contribute more in the northern part of the domain (35° to 45°S), and their contribution is less variable in both space (fig. S5) and time (not shown). The inclusion of the POLY + PROT fraction therefore reduces the overall spatiotemporal variability of OMF and of its effect on Nd in the regression model. Additional details about the OMF simulation can be found in the paper describing the parameterization (43).

Comparison of OMB14 to observationally derived parameterizations of OMF

The work of Burrows et al. (43) (B14) is different from previous methods of parameterizing the submicron SSA OMF in that it is a physically based framework as opposed to an empirical relationship derived from observations of Chl-a and OMF.

To assess whether parameterizations of OMF based on Chl-a can better predict Nd, B14 is compared to two such parameterizations (49, 79). The latter study [Gantt et al. (79)] includes empirical fits to data from both the Mace Head Atmospheric Research Station (53.33°N, 9.9°W) and Point Reyes National Sea Shore (38.12°N, 122.91°W). To evaluate if B14 significantly improves on empirical parameterizations, OMF consistent with each empirical fit was calculated using the POP/BEC ocean biogeochemistry-simulated chlorophyll rather than observed chlorophyll, to provide a fair comparison with B14. Because ocean biology will respond to future changes in ocean temperature, acidity, and sea ice extent, it is important to evaluate how well parameterizations can perform within the context of simulated ocean biology, which can dynamically change in response to ocean conditions. The correlations of Nd with the OMF calculated using each of these parameterizations or the POP Chl-a fields is low, with less than 20% of the variance explained by any of the parameterizations. By contrast, B14 explains 37 ± 10% of the variance in Nd as a single predictor (fig. S6). This is because B14 includes components driven by two categories of ocean biological variables with differing spatial and seasonal distributions: a “lipid-like” component that is strongly correlated with Chl-a and a “polysaccharide-like” component that is tied to the longer-lived semi-labile DOC, whose distribution is simulated by POP/BEC.

The prediction of biogeochemistry is challenging, and the BEC ocean model predicts a different distribution of Chl-a than is observed by SeaWiFS (Sea-Viewing Wide Field-of-View Sensor) (45). It is likely that the simulation of OMF by B14 is limited by the oceanic biogeochemistry model’s accuracy. To estimate the degree to which the B14 OMF might be improved by a more accurate simulation of oceanic biology, a relationship between Chl-a and the B14 lipid mass fraction has been derived using POP/BEC Chl-a (fig. S7). This relationship is used to calculate a distribution of the lipid-like class that is consistent with SeaWiFS-observed Chl-a.

The variance explained by OMF is calculated for several common empirical functions of Chl-a (49, 79). This is compared to the variance explained by OMF calculated by B14. Two cases are considered: Chl-a as observed by SeaWiFS and as simulated by the POP/BEC model. In both cases, the OMF consistent with the empirical parameterizations is calculated on the basis of Chl-a. In the former case, the B14 analogous OMF is calculated using the combination of polysaccharide-like OMF calculated from the POP/BEC model and lipid-like OMF calculated using SeaWiFS Chl-a and the function shown in fig. S7. In the latter case, OMF from B14 using the full simulation of ocean biogeochemistry is used. The variance explained by the OMF consistent with SeaWiFS Chl-a is shown in fig. S8, and the variance explained by OMF consistent with POP/BEC is shown in fig. S4. When observed distributions of Chl-a are used, B14-derived OMF performs as well as the upper end of empirical parameterizations of OMF. When POP/BEC Chl-a is used in combination with the empirical parameterizations to predict OMF, it explains significantly less variance than the B14-simulated OMF. OMF consistent with B14 explains roughly the same variance whether observed Chl-a is used as an input or the full oceanic biogeochemistry model is used. This suggests that the ability of the physically based framework used by B14 to simulate organic matter in SSA offers additional insight into the enrichment of SSA and will continue to improve as the ocean biogeochemistry model’s ability to accurately simulate Chl-a improves.

Statistical analysis

Regression modeling. The relationships between predictor variables and Nd are explored using multivariate linear regressions. Here, we present the coefficients for several additional regressions. In addition to the OMF from B14 (OMF) and the sulfate concentration from the AeroCom median (SO4), we also examine the regression based on the separation of SSA organic matter into lipid-like organic matter (LIPID) and polysaccharide-like organic matter (POLY), which have distinct seasonal and geographic patterns. Details of this partitioning are explained in section 5 (figs. S4 and S8) and in B14.

Multivariate fits to the Nd from Grosvenor and Wood (15) for the year 2007 with smoothed spatial interpolation (80, 81) (described in greater detail below) yield the following coefficients. Uncertainties are SDs of the coefficients obtained from 10-fold cross-validation (51), repeated 10 times. Model results indicate that the relationship between CCN and Nd should be approximately linear in the SO (82), allowing the use of simple linear regression in this region.

Nd = (115.6971 ± 2.7742 cm−3) OMF + 40.6347 ± 0.4087 cm−3

Nd = (176.6927 ± 2.4004 μg−1 m3 cm−3) SO4 + 40.4987 ± 0.2209 cm−3

Nd = (137.1134 ± 2.5679 μg−1 m3 cm−3) SO4 +

(58.2682 ± 1.9953 cm−3) OMF + 35.6441 ± 0.2233 cm−3

Nd = (139.1144 ± 2.9050 μg−1 m3 cm−3) SO4 + (58.0173 ± 1.6236 cm−3) LIPID + 35.9051 ± 0.2432 cm−3 + (300.0160 ± 27.1975 cm−3) POLY

Several other predictors were investigated to rule out possible confounding variables. AeroCom SS concentration and NOAA OI SST (72) each predict only a small fraction of the variability in Nd (fig. S6). A linear fit to monthly mean wind speed predicted about 20 ± 10% of the variability in Nd; however, a multivariate fit including wind speed, OMF, and SO4 concentration did not significantly increase the variance predicted by OMF and SO4 concentration alone (fig. S9). A climatology of ocean DMS concentration (50) explained about 25 ± 10% of the variance in Nd (fig. S6). Although not negligible, the variance explained by each of these variables is relatively slight compared to OMF and SO4, indicating that they are not the underlying drivers of the spatial and temporal structure of Nd in the SO, and none of these variables significantly increase predictability in a regression that also includes SO4 and OMF.

Details of data interpolation

Regressions were performed using monthly mean Nd data averaged over 5° latitude × 15° longitude bins. Data were interpolated spatially at its native resolution and then averaged. Nd data were left blank in bins where reliable retrievals could not be made for the whole bin. Only data from oceans not covered by sea ice were used as estimated by the HADISST sea ice coverage as included in the NOAA OI SST (72) data.

Observations of ocean and cloud properties from space are not always possible, leading to gaps in observational data sets. SeaWiFS requires a cloud-free view of the ocean, and robust retrievals of Nd are only possible at low SZA (15). Resulting data gaps in Nd and observed predictor variables (Chl-a, wind speed, SST, etc.) were filled using a smoothed interpolation procedure described in (80, 81). The interpolated data were then averaged in 5° latitude × 15° longitude bins. Analysis was repeated for 5° latitude × 5° longitude bins, and the variance explained was found to be fairly similar (fig. S10).


Supplementary material for this article is available at

Fig. S1. The AeroCom scenario B and PRE SO4 column density (black) across eight models (MPI-HAM excluded due to unavailability of SO4 column density) and difference in SO4 column density (μg/m2) between preindustrial and present day (red).

Fig. S2. The AeroCom scenario B and PRE SO4 surface concentration from LSCE, GISS, LOA, and UIO-GCM.

Fig. S3. The fractional change in SO4 surface concentration between the present-day and preindustrial period from AeroCom scenario B and PRE SO4 surface concentration from LSCE, GISS, LOA, and UIO-GCM.

Fig. S4. Tenfold cross-validation performed using only POP/BEC-modeled chlorophyll-consistent quantities.

Fig. S5. Contribution of the labile “lipid-like” group and total contribution of the semi-labile “polysaccharide-like” and “protein-like” groups to the total Nd from the regression model.

Fig. S6. The 10-fold cross-validation repeated using wind speed from the Multi-Platform Ocean Surface Wind Velocity L3.5, SST from the NOAA OI data product, SS aerosol concentration from AeroCom, and DMS (50).

Fig. S7. Scatter plot of lipid-like OMF in SSA as calculated by B14 as a function of POP chlorophyll.

Fig. S8. Tenfold cross-validation performed using only SeaWiFS consistent Chl-a.

Fig. S9. Tenfold cross-validation repeated using wind and SS concentration in addition to OMB14 and SO4 concentration.

Fig. S10. As in fig. S1, but averaged over 5° latitude × 5° longitude bins.

Table S1. Correlation coefficients in space and time using labile and semi-labile fields from B14.

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: We wish to thank the AeroCom intercomparison project for making the sea salt and sulfate model median data used in this study available, and the modeling centers whose models constitute the median data set: UMI, UIO-CTM, LOA, LSCE, MPI-HAM, GISS, UIO-GCM, SPRINTARS, and ULAQ. CERES data were obtained from the NASA Langley Research Center CERES ordering tool. MODIS data were obtained from the NASA Goddard Land Processes data archive. We thank A. Conley and F. Vitt for assistance setting up the PORT offline radiative transfer calculations; M. Maltrud for providing the POP/BEC ocean biogeochemistry simulation data underlying the B14 parameterization; and B. Singh for technical support at PNNL. We also wish to thank the three anonymous reviewers, R. Charlson, C. Gaston, A. McCoy, and J. McCoy, for useful discussion and advice. POP/BEC, PORT, and CAM5 are components of the Community Earth System Model (CESM). The CESM project is supported by the National Science Foundation and the Office of Science of the U.S. Department of Energy. The PORT radiative transfer calculations were performed using EMSL, a DOE Office of Science User Facility sponsored by the Office of Biological and Environmental Research (BER) and located at Pacific Northwest National Laboratory. Funding: We wish to acknowledge the following funding sources: D.L.H. and D.T.M. were supported under NASA Grant NNX14AG26G, and D.T.M. acknowledges government support awarded by DoD, Air Force Office of Scientific Research, National Defense Science and Engineering Graduate (NDSEG) Fellowship, 32 CFR 168a. S.M.B., P.-L.M., and P.J.R. were supported by the U.S. Department of Energy, Office of Science, BER as part of the Earth System Modeling program. D.P.G. was supported by the School of Earth and Environment at the University of Leeds. R.W. was supported by NASA Grant NNX13AQ35G. S.M.E. was supported by several DOE SciDAC and SFA earth system modeling projects (Scientific Discovery through Advanced Computing, Science Focus Areas). Author contributions: D.T.M. and S.M.B. performed analysis and prepared figures and the manuscript. R.W., D.P.G., S.M.E., P.-L.M., P.J.R., and D.L.H. assisted with the analyses and preparation of the manuscript. All authors proofread and corrected the final version. Competing interests: The authors declare that they have no competing interests.
View Abstract

Stay Connected to Science Advances

Navigate This Article