Search

The effective radiative forcing of anthropogenic aerosols (ERFaer) is an important measure of the anthropogenic aerosol effects simulated by a global climate model. Here we analyze ERFaer simulated by the E3SM version 1 (E3SMv1) atmospheric model using both century-long free-running atmosphere–land simulations and short nudged simulations. We relate the simulated ERFaer to characteristics of the aerosol composition and optical properties, and we evaluate the relationships between key aerosol and cloud properties. In terms of historical changes from the year 1870 to 2014, our results show that the global mean anthropogenic aerosol burden and optical depth increase during the simulation period as expected, but the regional averages show large differences in the temporal evolution. The largest regional differences are found in the emission-induced evolution of the burden and optical depth of the sulfate aerosol: a strong decreasing trend is seen in the Northern Hemisphere high-latitude region after around 1970, while a continued increase is simulated in the tropics. The relationships between key aerosol and cloud properties (relative changes between pre-industrial and present-day conditions) also show evident changes after 1970, diverging from the linear relationships exhibited for the period of 1870–1969. In addition to the regional differences in the simulated relationships, a reduced sensitivity in cloud droplet number and other cloud properties to aerosol perturbations is seen when the aerosol perturbation is large. Consequently, the global annual mean ERFaer magnitude does not increase after around 1970. The ERFaer in E3SMv1 is relatively large compared to the recently published multi-model estimates; the primary reason is the large indirect aerosol effect (i.e., through aerosol–cloud interactions). Compared to other models, E3SMv1 features large relative changes in the cloud droplet effective radius in response to aerosol perturbations. Large sensitivity is also seen in the liquid cloud optical depth, which is determined by changes in both the effective radius and liquid water path. Aerosol-induced changes in liquid and ice cloud properties in E3SMv1 are found to have a strong correlation, as the evolution of anthropogenic sulfate aerosols affects both the liquid cloud formation and the homogeneous ice nucleation in cirrus clouds (that causes a large effect on longwave ERFaer). As suggested by a previous study, the large ERFaer appears to be one of the reasons why the model cannot reproduce the observed global mean temperature evolution in the second half of the 20th century. Sensitivity simulations are performed to understand which parameterization and/or parameter changes have a large impact on the simulated ERFaer. The ERFaer estimates in E3SMv1 for the shortwave and longwave components are sensitive to the parameterization changes in both liquid and ice cloud processes. When the parameterization of ice cloud processes is modified, the top-of-model forcing changes in the shortwave and longwave components largely offset each other, so the net effect is negligible. This suggests that, to reduce the magnitude of the net ERFaer, it would be more effective to reduce the anthropogenic aerosol effect through liquid or mixed-phase clouds.

An advanced aerosol treatment, with a focus on semivolatile nitrate formation, is introduced into the Community Atmosphere Model version 5 with interactive chemistry (CAM5‐chem) by coupling the Model for Simulating Aerosol Interactions and Chemistry (MOSAIC) with the 7‐mode Modal Aerosol Module (MAM7). An important feature of MOSAIC is dynamic partitioning of all condensable gases to the different fine and coarse mode aerosols, as governed by mode‐resolved thermodynamics and heterogeneous chemical reactions. Applied in the free‐running mode from 1995 to 2005 with prescribed historical climatological conditions, the model simulates global distributions of sulfate, nitrate, and ammonium in good agreement with observations and previous studies. Inclusion of nitrate resulted in ∼10% higher global average accumulation mode number concentrations, indicating enhanced growth of Aitken mode aerosols from nitrate formation. While the simulated accumulation mode nitrate burdens are high over the anthropogenic source regions, the sea‐salt and dust modes respectively constitute about 74% and 17% of the annual global average nitrate burden. Regional clear‐sky shortwave radiative cooling of up to −5 W m−2 due to nitrate is seen, with a much smaller global average cooling of −0.05 W m−2. Significant enhancements in regional cloud condensation nuclei (at 0.1% supersaturation) and cloud droplet number concentrations are also attributed to nitrate, causing an additional global average shortwave cooling of −0.8 W m−2. Taking into consideration of changes in both longwave and shortwave radiation under all‐sky conditions, the net change in the top of the atmosphere radiative fluxes induced by including nitrate aerosol is −0.7 W m−2., Key Points A modal version of the advanced aerosol chemistry module MOSAIC is developed and introduced in a climate model to simulate nitrate aerosolMOSAIC provides an accurate and efficient treatment for dynamically partitioning semivolatile gases over to entire aerosol size distributionThe modeled global distribution of nitrate is in good agreement with observations and its impact on the radiative effects is quantified

With increasing computational capabilities, cumulus parameterizations that are adaptable to the smaller grid spacing and temporal interval for high-resolution climate model simulations are needed. In this study, we propose a method to improve the resolution adaptability of the Zhang-McFarlane (ZM) scheme, by implementing spatial and temporal averaging to the CAPE tendency. This method allows for a more consistent application of the quasi-equilibrium (QE) hypothesis at high spatial and temporal resolutions. The resolution adaptability of the original ZM scheme, the scheme with spatial averaging, and the scheme with spatiotemporal averaging at 4–32 km grid spacings are assessed using the Weather Research and Forecasting (WRF) model by comparing to cloud resolving model (CRM) simulation results coarse-grained to these same grid spacings. We show the original ZM scheme has poor resolution adaptability, with spatiotemporally averaged subgrid convective transport and convective precipitation increasing significantly as the resolution increases. The spatial averaging method improves the resolution adaptability of the ZM scheme and better conserves total transport and total precipitation. Temporal averaging further improves the resolution adaptability of the scheme. With better constrained (although smoothed) convective transport and precipitation, both the spatial distribution and time series of total precipitation at 4 and 8 km grid spacings are improved with the averaging methods. The results could help develop resolution adaptability for other cumulus parameterizations that are based on the QE assumption.

We evaluate the Community Atmosphere Model Version 5 (CAM5) with a higher-order turbulence closure scheme, named Cloud Layers Unified By Binomials (CLUBB), and a Multiscale Modeling Framework, referred as the “super-parameterization” (SP) with two different microphysics configurations to investigate their influences on rainfall simulations over Southern Amazonia. The two different microphysics configurations in SP are the one-moment cloud microphysics without aerosol treatment (SP1) and two-moment cloud microphysics coupled with aerosol treatment (SP2). Results show that both SP2 and CLUBB effectively reduce the low biases of rainfall, mainly during the wet season, and reduce low biases of humidity in the lower troposphere with further reduced shallow clouds and increased surface solar flux. These changes increase moist static energy in the lower atmosphere, contribute to stronger convection and more rainfall. SP2 appears to realistically capture the observed increase of relative humidity prior to deep convection and it significantly increases rainfall in the afternoon; CLUBB significantly delays the afternoon peak rainfall and produces more precipitation in the early morning, due to more gradual transition between shallow and deep convection. In CAM5 and CAM5 with CLUBB, occurrence of more deep convection appears to be a result of stronger heating rather than higher relative humidity.

We used two 150-year pre-industrial simulations of the Community Earth System Model (CESM), one with interactive dust and the other with prescribed dust, to quantify the impacts of changes in wind during East Asian Winter Monsoon (EAWM) season on dust emissions, and the resulting consequences for interannual variations of temperature and precipitation over East Asia. The simulated December-January-February dust column burden and dust optical depth are lower over northern China in the strongest EAWM years than those of the weakest years by 38.3% and 37.2%, respectively. The decrease in dust over the dust source regions and the downwind region leads to an increase in direct radiative forcing (RF) at the surface by up to 1.5 Wm-2. The effects of EAWM-related variations in surface winds, precipitation and their effects on dust emissions and wet removal contribute 67% to the total dust-induced variations of direct RF at the surface and partly offset the cooling that occurs with the EAWM strengthening by heating the surface. The variations of surface air temperature induced by the changes in wind and dust emissions between the strongest and weakest EAWM years (strongest minus weakest) decrease by 0.4-0.6 K from eastern coastal China to Japan, which weakens the impact of EAWM on surface air temperature by 3–18% in these regions. The warming results from the combined effects of changes in direct RF, turbulent heat flux at the surface, and northwesterly wind anomalies that bring cold and dry air from Siberia to these regions. Over eastern coastal China, the variations of large-scale precipitation induced by the feedback of EAWM-related changes in wind on dust emissions decrease by 10-30% in winter because of the reducing changes in surface air temperature and the anomalous circulation.

One limitation of most global climate models (GCMs) is that with the horizontal resolutions they typically employ, they cannot resolve the sub-grid variability (SGV) of clouds and aerosols, adding extra uncertainties to the aerosol radiative forcing estimation. To inform the development of an aerosol sub-grid variability parameterization, here we analyze the aerosol SGV over the southern Pacific Ocean simulated by the high-resolution Weather Research and Forecasting model coupled to Chemistry (WRF-Chem). We find that within a typical GCM grid, the aerosol mass sub-grid standard deviation is 15% of the grid-box mean mass near the surface on a 1-month mean basis. The fraction can increase to 50% in the free troposphere. The relationships between the sea-salt mass concentration, meteorological variables, and sea-salt emission rate are investigated in both the clear and cloudy portion. Under clear-sky conditions, marine aerosol sub-grid standard deviation is highly correlated with the standard deviations of vertical velocity, cloud water mixing ratio, and sea-salt emission rates near the surface. It is also strongly connected to the grid box mean aerosol in the free troposphere (between 2 km to 4 km). In the cloudy area, interstitial sea-salt aerosol mass concentrations are smaller, but higher correlation is found between the sub-grid standard deviations of aerosol mass and vertical velocity. Additionally, we find that decreasing the model grid resolution can reduce the marine aerosol SGV but strengthen the correlations between the aerosol SGV and the total water mixing ratio (sum of water vapor, cloud liquid, and cloud ice mixing ratios).

From April 2009 to December 2010, the Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) program carried out an observational field campaign on Graciosa Island, targeting the marine boundary layer (MBL) clouds over the Azores region. In this paper, we present an inter-comparison of the MBL cloud properties, namely, cloud liquid water path (LWP), cloud optical thickness (COT) and cloud-droplet effective radius (CER), among retrievals from the ARM mobile facility (AMF) and two Moderate Resolution Spectroradiometer (MODIS) cloud products (GSFC-MODIS and CERES-MODIS). A total of 63 daytime single-layer MBL cloud cases are selected for inter-comparison. Comparison of collocated retrievals indicates that the two MODIS cloud products agree well on both COT and CER retrievals, with the correlation coefficient R>0.95. despite their significant difference in spatial sampling. In both MODIS products, the CER retrievals based on the 2.1 µm band (CER(2.1)) is significantly smaller than that based on the 3.7 µm band (CER(3.7)). The GSFC-MODIS cloud product is collocated and compared with ground-based ARM observations at several temporal-spatial scales. In general, the correlation increases with more precise collocation. For the 63 selected MBL cloud cases, the GSFC-MODIS LWP and COT retrievals agree reasonably well with the ground-based observations with no apparent bias and correlation coefficient R around 0.85 and 0.70, respectively. However, GSFC-MODIS CER(3.7) and CER(2.1) retrievals have a lower correlation (R~0.5) with the ground-based retrievals. For the 63 selected cases, they are on average larger than ground observations by about 1.5 µm and 3.0 µm, respectively. Taking into account that the MODIS CER retrievals are only sensitive to cloud top reduces the bias only by 0.5 µm.

Aerosol-cloud interaction processes can be represented more physically with bin cloud microphysics relative to bulk microphysical parameterizations. However, due to computational power limitations in the past, bin cloud microphysics was often run with very simple aerosol treatments. The purpose of this study is to represent better aerosol-cloud interaction processes in the Chemistry version of Weather Research and Forecast model (WRF-Chem) at convection-permitting scales by coupling spectral-bin cloud microphysics (SBM) with the MOSAIC sectional aerosol model. A flexible interface is built that exchanges cloud and aerosol information between them. The interface contains a new bin aerosol activation approach, which replaces the treatments in the original SBM. It also includes the modified aerosol resuspension and in-cloud wet removal processes with the droplet loss tendencies and precipitation fluxes from SBM. The newly coupled system is evaluated for two marine stratocumulus cases over the Southeast Pacific Ocean with either a simplified aerosol setup or full-chemistry. We compare the aerosol activation process in the newly coupled SBM-MOSAIC against the SBM simulation without chemistry using a simplified aerosol setup, and the results show consistent activation rates. A longer time simulation reinforces that aerosol resuspension through cloud drop evaporation plays an important role in replenishing aerosols and impacts cloud and precipitation in marine stratocumulus clouds. Evaluation of the coupled SBM-MOSAIC with full-chemistry using aircraft measurements suggests that the new model works realistically for the marine stratocumulus clouds, and improves the simulation of cloud microphysical properties compared to a simulation using MOSAIC coupled with the Morrison two-moment microphysics.

We use 150 year preindustrial simulations of the Community Earth System Model to quantify the impacts of the East Asian Monsoon strength on interannual variations of springtime dust concentrations over China. The simulated interannual variations in March-April-May (MAM) dust column concentrations range between 20–40% and 10–60% over eastern and western China, respectively. The dust concentrations over eastern China correlate negatively with the East Asian Monsoon (EAM) index, which represents the strength of monsoon, with a regionally averaged correlation coefficient of −0.64. Relative to the strongest EAM years, MAM dust concentrations in the weakest EAM years are higher over China, with regional relative differences of 55.6%, 29.6%, and 13.9% in the run with emissions calculated interactively and of 33.8%, 10.3%, and 8.2% over eastern, central, and western China, respectively, in the run with prescribed emissions. Both interactive run and prescribed emission run show the similar pattern of climate change between the weakest and strongest EAM years. Strong anomalous northwesterly and westerly winds over the Gobi and Taklamakan deserts during the weakest EAM years result in larger transport fluxes, and thereby increase the dust concentrations over China. These differences in dust concentrations between the weakest and strongest EAM years (weakest-strongest) lead to the change in the net radiative forcing by up to −8 and −3 W m−2 at the surface, compared to −2.4 and +1.2 W m−2 at the top of the atmosphere over eastern and western China, respectively.

Efficient simulation strategies are crucial for the development and evaluation of high-resolution climate models. This paper evaluates simulations with constrained meteorology for the quantification of parametric sensitivities in the Community Atmosphere Model version 5 (CAM5). Two parameters are perturbed as illustrating examples: the convection relaxation time scale (TAU), and the threshold relative humidity for the formation of low-level stratiform clouds (rhminl). Results suggest that the fidelity of the constrained simulations depends on the detailed implementation of nudging and the mechanism through which the perturbed parameter affects precipitation and cloud. The relative computational costs of nudged and free-running simulations are determined by the magnitude of internal variability in the physical quantities of interest, as well as the magnitude of the parameter perturbation. In the case of a strong perturbation in convection, temperature and/or wind nudging with a 6-hour relaxation time scale leads to non-negligible side effects due to the distorted interactions between resolved dynamics and parameterized convection, while 1-year free running simulations can satisfactorily capture the annual mean precipitation and cloud forcing sensitivities. In the case of a relatively weak perturbation in the large-scale condensation scheme, results from 1-year free-running simulations are strongly affected by natural noise, while nudging winds effectively reduces the noise, and reasonably reproduces the sensitivities. These results indicate that caution is needed when using nudged simulations to assess precipitation and cloud forcing sensitivities to parameter changes in general circulation models. We also demonstrate that ensembles of short simulations are useful for understanding the evolution of model sensitivities. This article is protected by copyright. All rights reserved.

The ability of 11 models in simulating the aerosol vertical distribution from regional to global scales, as part of the second phase of the AeroCom model intercomparison initiative (AeroCom II), is assessed and compared to results of the first phase. The evaluation is performed using a global monthly gridded data set of aerosol extinction profiles built for this purpose from the CALIOP (Cloud-Aerosol Lidar with Orthogonal Polarization) Layer Product 3.01. Results over 12 subcontinental regions show that five models improved, whereas three degraded in reproducing the interregional variability in Zα0–6 km, the mean extinction height diagnostic, as computed from the CALIOP aerosol profiles over the 0–6 km altitude range for each studied region and season. While the models' performance remains highly variable, the simulation of the timing of the Zα0–6 km peak season has also improved for all but two models from AeroCom Phase I to Phase II. The biases in Zα0–6 km are smaller in all regions except Central Atlantic, East Asia, and North and South Africa. Most of the models now underestimate Zα0–6 km over land, notably in the dust and biomass burning regions in Asia and Africa. At global scale, the AeroCom II models better reproduce the Zα0–6 km latitudinal variability over ocean than over land. Hypotheses for the performance and evolution of the individual models and for the intermodel diversity are discussed. We also provide an analysis of the CALIOP limitations and uncertainties contributing to the differences between the simulations and observations.

PUBLICATIONS Journal of Geophysical Research: Atmospheres RESEARCH ARTICLE 10.1002/2015JD024503 Key Points: • Interannual variability in cloud radiative effects is driven by midlevel and high clouds • Wind-related feedbacks on natural aerosol emissions enhance this variability by 3 to 5% • Variations in natural aerosol concentrations enhance interannual variability by 1 to 3% Supporting Information: • Supporting Information S1 Correspondence to: L. M. Russell, lmrussell@ucsd.edu Citation: Yang, Y., et al. (2016), Impacts of ENSO events on cloud radiative effects in preindustrial conditions: Changes in cloud fraction and their dependence on interactive aerosol emissions and concentrations, J. Geophys. Res. Atmos., 121, 6321–6335, doi:10.1002/ 2015JD024503. Received 13 NOV 2015 Accepted 16 MAY 2016 Accepted article online 19 MAY 2016 Published online 2 JUN 2016 Impacts of ENSO events on cloud radiative effects in preindustrial conditions: Changes in cloud fraction and their dependence on interactive aerosol emissions and concentrations Yang Yang 1,2 , Lynn M. Russell 1 , Li Xu 1 , Sijia Lou 1 , Maryam A. Lamjiri 1 , Richard C. J. Somerville 1 , Arthur J. Miller 1 , Daniel R. Cayan 1 , Michael J. DeFlorio 1 , Steven J. Ghan 3 , Ying Liu 3 , Balwinder Singh 3 , Hailong Wang 3 , Jin-Ho Yoon 4 , and Philip J. Rasch 3 Scripps Institution of Oceanography, University of California, San Diego, La Jolla, California, USA, 2 Now at Atmospheric Science and Global Change Division, Pacific Northwest National Laboratory, Richland, Washington, USA, 3 Atmospheric Science and Global Change Division, Pacific Northwest National Laboratory, Richland, Washington, USA, 4 Gwangju Institute of Science and Technology, Gwangju, South Korea Abstract We use three 150 year preindustrial simulations of the Community Earth System Model to quantify the impacts of El Nino–Southern Oscillation (ENSO) events on shortwave and longwave cloud radiative effects (CRE SW and CRE LW ). Compared to recent observations from the Clouds and the Earth’s Radiant Energy System data set, the model simulation successfully reproduces larger variations of CRE SW and CRE LW over the tropics. The ENSO cycle is found to dominate interannual variations of cloud radiative effects. Simulated cooling (warming) effects from CRE SW (CRE LW ) are strongest over the tropical western and central Pacific Ocean during warm ENSO events, with the largest difference between 20 and 60 W m 2 , with weaker effects of 10–40 W m 2 over Indonesian regions and the subtropical Pacific Ocean. Sensitivity tests show that variations of cloud radiative effects are mainly driven by ENSO-related changes in cloud fraction. The variations in midlevel and high cloud fractions each account for approximately 20–50% of the interannual variations of CRE SW over the tropics and almost all of the variations of CRE LW between 60°S and 60°N. The variation of low cloud fraction contributes to most of the variations of CRE SW over the midlatitude oceans. Variations in natural aerosol concentrations explained 10–30% of the variations of both CRE SW and CRE LW over the tropical Pacific, Indonesian regions, and the tropical Indian Ocean. Changes in natural aerosol emissions and concentrations enhance 3–5% and 1–3% of the variations of cloud radiative effects averaged over the tropics. 1. Introduction Clouds strongly influence the Earth’s radiation balance. They reflect incoming solar radiation back to space, which enhances the reflected solar flux by 47.5 ± 3 W m 2 globally, and absorb outgoing infrared radiation, which reduces the outgoing longwave flux relative to clear sky by approximately 26.4 ± 4 W m 2 . Overall, clouds exert a net cooling effect of about 21.1 ± 5 W m 2 at the top of atmosphere (TOA) [Stephens et al., 2012], which is 6 times larger than that from doubling CO 2 concentration [Ramanathan et al., 1989; Loeb et al., 2009]. Any changes in cloud properties such as cloud fraction, cloud top height, and microphysical fea- tures would perturb cloud radiative forcing and greatly modulate the radiative balance of the Earth system [Slingo, 1990; Wielicki et al., 1998; Curry et al., 2000; Stephens, 2005]. The Intergovernmental Panel on Climate Change reported that simulations of clouds and their radiative feedbacks are still one of the largest uncertainties in the fifth-generation climate models [Boucher et al., 2013]. ©2016. American Geophysical Union. All Rights Reserved. YANG ET AL. On interannual time scales, many regional changes in the global climate system are associated with the El Nino–Southern Oscillation (ENSO). ENSO is characterized by anomalous sea surface temperatures (SSTs) in the equatorial Pacific and has far-reaching impacts on global and regional temperature, precipitation, and circulation. Using cloud data from the Extended Edited Cloud Reports Archive (EECRA) from year 1954 to recent years, Park and Leovy [2004] and Eastman et al. [2011] both showed that interannual variations of cloud cover in the tropics have strong correlations to the ENSO index. For example, warmer central and eastern tropical Pacific SST (warm ENSO phase, i.e., El Nino) is associated with increased cloud cover in the tropical central Pacific Ocean and reduced cloud cover over the Indonesian and eastern Pacific regions, and vice versa for cool ENSO phase (La Nina) events. IMPACT OF ENSO ON CLOUD RADIATIVE EFFECT

Accurate representation of global stratospheric aerosols from volcanic and non-volcanic sulfur emissions is key to understanding the cooling effects and ozone-losses that may be linked to volcanic activity. Attribution of climate variability to volcanic activity is of particular interest in relation to the post-2000 slowing in the rate of global average temperature increases. We have compiled a database of volcanic SO2 emissions and plume altitudes for eruptions from 1990 to 2014, and developed a new prognostic capability for simulating stratospheric sulfate aerosols in the Community Earth System Model (CESM). We used these combined with other non-volcanic emissions of sulfur sources to reconstruct global aerosol properties from 1990 to 2014. Our calculations show remarkable agreement with ground-based lidar observations of stratospheric aerosol optical depth (SAOD), and with in situ measurements of stratospheric aerosol surface area density (SAD). These properties are key parameters in calculating the radiative and chemical effects of stratospheric aerosols. Our SAOD calculations represent a clear improvement over available satellite-based analyses, which generally ignore aerosol extinction below 15 km, a region that can contain the vast majority of stratospheric aerosol extinction at mid- and high-latitudes. Our SAD calculations greatly improve on that provided for the Chemistry-Climate Model Initiative, which misses about 60% of the SAD measured in situ on average during both volcanically active and volcanically quiescent periods.

A large number of processes are involved in the chain from emissions of aerosol precursor gases and primary particles to impacts on cloud radiative forcing. Those processes are manifest in a number of relationships that can be expressed as factors dlnX/dlnY driving aerosol effects on cloud radiative forcing. These factors include the relationships between cloud condensation nuclei (CCN) concentration and emissions, droplet number and CCN concentration, cloud fraction and droplet number, cloud optical depth and droplet number, and cloud radiative forcing and cloud optical depth. The relationship between cloud optical depth and droplet number can be further decomposed into the sum of two terms involving the relationship of droplet effective radius and cloud liquid water path with droplet number. These relationships can be constrained using observations of recent spatial and temporal variability of these quantities. However, we are most interested in the radiative forcing since the preindustrial era. Because few relevant measurements are available from that era, relationships from recent variability have been assumed to be applicable to the preindustrial to present-day change. Our analysis of Aerosol Comparisons between Observations and Models (AeroCom) model simulations suggests that estimates of relationships from recent variability are poor constraints on relationships from anthropogenic change for some terms, with even the sign of some relationships differing in many regions. Proxies connecting recent spatial/temporal variability to anthropogenic change, or sustained measurements in regions where emissions have changed, are needed to constrain estimates of anthropogenic aerosol impacts on cloud radiative forcing.

This paper evaluates the impact of sub-grid variability of surface wind on sea salt and dust emissions in the Community Atmosphere Model version 5 (CAM5). The basic strategy is to calculate emission fluxes multiple times, using different wind speed samples of a Weibull probability distribution derived from model-predicted grid-box mean quantities. In order to derive the Weibull distribution, the sub-grid standard deviation of surface wind speed is estimated by taking into account four mechanisms: turbulence under neutral and stable conditions, dry convective eddies, moist convective eddies over the ocean, and air motions induced by mesoscale systems and fine-scale topography over land. The contributions of turbulence and dry convective eddy are parameterized using schemes from the literature. Wind variabilities caused by moist convective eddies and fine-scale topography are estimated using empirical relationships derived from an operational weather analysis data set at 15 km resolution. The estimated sub-grid standard deviations of surface wind speed agree well with reference results derived from 1 year of global weather analysis at 15 km resolution and from two regional model simulations with 3 km grid spacing.The wind-distribution-based emission calculations are implemented in CAM5. In terms of computational cost, the increase in total simulation time turns out to be less than 3 %. Simulations at 2° resolution indicate that sub-grid wind variability has relatively small impacts (about 7 % increase) on the global annual mean emission of sea salt aerosols, but considerable influence on the emission of dust. Among the considered mechanisms, dry convective eddies and mesoscale flows associated with topography are major causes of dust emission enhancement. With all the four mechanisms included and without additional adjustment of uncertain parameters in the model, the simulated global and annual mean dust emission increase by about 50 % compared to the default model. By tuning the globally constant dust emission scale factor, the global annual mean dust emission, aerosol optical depth, and top-of-atmosphere radiative fluxes can be adjusted to the level of the default model, but the frequency distribution of dust emission changes, with more contribution from weaker wind events and less contribution from stronger wind events. In Africa and Asia, the overall frequencies of occurrence of dust emissions increase, and the seasonal variations are enhanced, while the geographical patterns of the emission frequency show little change.

Wet scavenging of aerosols by continental deep convective clouds is studied for a supercell storm complex observed over Oklahoma during the Deep Convective Clouds and Chemistry campaign. A new passive-tracer-based transport analysis framework is developed to characterize convective transport using vertical profiles of several passive trace gases. For this case, the analysis estimates that observed passive gas mixing ratios in the upper troposphere convective outflow consist of 47% low level (

In this study, we investigate the sensitivity of simulated low clouds to 14 selected tunable parameters of Cloud Layers Unified By Binormals (CLUBB), a higher order closure (HOC) scheme, and 4 parameters of the Zhang-McFarlane (ZM) deep convection scheme in the Community Atmosphere Model version 5 (CAM5). A quasi-Monte Carlo (QMC) sampling approach is adopted to effectively explore the high-dimensional parameter space and a generalized linear model is applied to study the responses of simulated cloud fields to tunable parameters. Our results show that the variance in simulated low-cloud properties (cloud fraction and liquid water path) can be explained by the selected tunable parameters in two different ways: macrophysics itself and its interaction with microphysics. First, the parameters related to dynamic and thermodynamic turbulent structure and double Gaussians closure are found to be the most influential parameters for simulating low clouds. The spatial distributions of the parameter contributions show clear cloud-regime dependence. Second, because of the coupling between cloud macrophysics and cloud microphysics, the coefficient of the dissipation term in the total water variance equation is influential. This parameter affects the variance of in-cloud cloud water, which further influences microphysical process rates, such as autoconversion, and eventually low-cloud fraction. Thismore » study improves understanding of HOC behavior associated with parameter uncertainties and provides valuable insights for the interaction of macrophysics and microphysics.« less

The Community Atmosphere Model Version 5 is run at horizontal grid spacing of 2, 1, 0.5, and 0.25°, with the meteorology nudged toward the Year Of Tropical Convection analysis, and cloud simulators and the collocated A-Train satellite observations are used to explore the resolution dependence of aerosol-cloud interactions. The higher-resolution model produces results that agree better with observations, showing an increase of susceptibility of cloud droplet size, indicating a stronger first aerosol indirect forcing (AIF), and a decrease of susceptibility of precipitation probability, suggesting a weaker second AIF. The resolution sensitivities of AIF are attributed to those of droplet nucleation and precipitation parameterizations. The annual average AIF in the Northern Hemisphere midlatitudes (where most anthropogenic emissions occur) in the 0.25° model is reduced by about 1 W m−2 (−30%) compared to the 2° model, leading to a 0.26 W m−2 reduction (−15%) in the global annual average AIF.

Secondary organic aerosols (SOA) are large contributors to fine-particle loadings and radiative forcing but are often represented crudely in global models. We have implemented three new detailed SOA treatments within the Community Atmosphere Model version 5 (CAM5) that allow us to compare the semivolatile versus nonvolatile SOA treatments (based on some of the latest experimental findings) and to investigate the effects of gas-phase fragmentation reactions. The new treatments also track SOA from biomass burning and biofuel, fossil fuel, and biogenic sources. For semivolatile SOA treatments, fragmentation reactions decrease the simulated annual global SOA burden from 7.5 Tg to 1.8 Tg. For the nonvolatile SOA treatment with fragmentation, the burden is 3.1 Tg. Larger differences between nonvolatile and semivolatile SOA (up to a factor of 5) exist in areas of continental outflow over the oceans. According to comparisons with observations from global surface Aerosol Mass Spectrometer measurements and the U.S. Interagency Monitoring of Protected Visual Environments (IMPROVE) network measurements, the FragNVSOA treatment, which treats SOA as nonvolatile and includes gas-phase fragmentation reactions, agrees best at rural locations. Urban SOA is underpredicted, but this may be due to the coarse model resolution. All three revised treatments show much better agreement with aircraft measurements of organic aerosols (OA) over the North American Arctic and sub-Arctic in spring and summer, compared to the standard CAM5 formulation. This is mainly due to the oxidation of SOA precursor gases from biomass burning, not included in standard CAM5, and long-range transport of biomass burning OA at high altitudes. The revised model configurations that include fragmentation (both semivolatile and nonvolatile SOA) show much better agreement with MODerate resolution Imaging Spectrometers (MODIS) aerosol optical depth data over regions dominated by biomass burning during the summer compared to standard CAM5, and predict biomass burning and biofuel as the largest global source of OA, followed by biogenic and fossil fuel sources. The large contribution of biomass burning OA in the revised treatments is supported by these measurements, but the emissions and aging of SOA precursors and POA are uncertain, and need further investigation. The nonvolatile and semivolatile configurations with fragmentation predict the direct radiative forcing of SOA as −0.5 W m−2 and −0.26 W m−2 respectively, at top of the atmosphere, which are higher than previously estimated by most models, but in reasonable agreement with a recent constrained modeling study. This study highlights the importance of improving process-level representation of SOA in global models.

Following Part I, in which 3-D cloud-resolving model (CRM) simulations of a squall line and mesoscale convective complex in the mid-latitude continental and the tropical regions are conducted and evaluated, we examine the scale-dependence of eddy transport of water vapor, evaluate different eddy transport formulations, and improve the representation of convective transport across all scales by proposing a new formulation that more accurately represents the CRM-calculated eddy flux. CRM results show that there are strong grid-spacing dependencies of updraft and downdraft fractions regardless of altitudes, cloud life stage, and geographical location. As for the eddy transport of water vapor, updraft eddy flux is a major contributor to total eddy flux in the lower and middle troposphere. However, downdraft eddy transport can be as large as updraft eddy transport in the lower atmosphere especially at the mature stage of 38 mid-latitude continental convection. We show that the single updraft approach significantly underestimates updraft eddy transport of water vapor because it fails to account for the large internal variability of updrafts, while a single downdraft represents the downdraft eddy transport of water vapor well. We find that using as few as 3 updrafts can account for the internal variability of updrafts well.more » Based on evaluation with the CRM simulated data, we recommend a simplified eddy transport formulation that considers three updrafts and one downdraft. Such formulation is similar to the conventional one but much more accurately represents CRM-simulated eddy flux across all grid scales.« less

The ultimate goal of this study is to improve the representation of convective transport by cumulus parameterization for mesoscale and climate models. As Part 1 of the study, we perform extensive evaluations of cloud-resolving simulations of a squall line and mesoscale convective complexes in midlatitude continent and tropical regions using the Weather Research and Forecasting model with spectral bin microphysics (SBM) and with two double-moment bulk microphysics schemes: a modified Morrison (MOR) and Milbrandt and Yau (MY2). Compared to observations, in general, SBM gives better simulations of precipitation and vertical velocity of convective cores than MOR and MY2 and therefore will be used for analysis of scale dependence of eddy transport in Part 2. The common features of the simulations for all convective systems are (1) the model tends to overestimate convection intensity in the middle and upper troposphere, but SBM can alleviate much of the overestimation and reproduce the observed convection intensity well; (2) the model greatly overestimates Ze in convective cores, especially for the weak updraft velocity; and (3) the model performs better for midlatitude convective systems than the tropical system. The modeled mass fluxes of the midlatitude systems are not sensitive to microphysics schemes but are very sensitive for the tropical case indicating strong microphysics modification to convection. Cloud microphysical measurements of rain, snow, and graupel in convective cores will be critically important to further elucidate issues within cloud microphysics schemes.

In this study, a higher-order turbulence closure scheme, called Cloud Layers Unified By Binormals (CLUBB), is implemented into a Multiscale Modeling Framework (MMF) model to improve low-cloud simulations. The performance of CLUBB in MMF simulations with two different microphysics configurations (one-moment cloud microphysics without aerosol treatment and two-moment cloud microphysics coupled with aerosol treatment) is evaluated against observations and further compared with results from the Community Atmosphere Model, Version 5 (CAM5) with conventional cloud parameterizations. CLUBB is found to improve low-cloud simulations in the MMF, and the improvement is particularly evident in the stratocumulus-to-cumulus transition regions. Compared to the single-moment cloud microphysics, CLUBB with two-moment microphysics produces clouds that are closer to the coast and agrees better with observations. In the stratocumulus-to-cumulus transition regions, CLUBB with two-moment cloud microphysics produces short-wave cloud forcing in better agreement with observations, while CLUBB with single-moment cloud microphysics overestimates short-wave cloud forcing. CLUBB is further found to produce quantitatively similar improvements in the MMF and CAM5, with slightly better performance in the MMF simulations (e.g., MMF with CLUBB generally produces low clouds that are closer to the coast than CAM5 with CLUBB). Improved low-cloud simulations in MMF make it an even more attractive tool for studying aerosol-cloud-precipitation interactions.

Dust variability in the climate system has been studied for several decades, yet there remains an incomplete understanding of the dynamical mechanisms controlling interannual and decadal variations in dust transport. The sparseness of multi-year observational datasets has limited our understanding of the relationship between climate variations and atmospheric dust. We use available in situ and satellite observations of dust and a century-length fully coupled Community Earth System Model (CESM) simulation to show that the El Nino-Southern Oscillation (ENSO) exerts a control on North African dust transport during boreal summer. In CESM, this relationship is stronger over the dusty tropical North Atlantic than near Barbados, one of the few sites having a multi-decadal observed record. During strong La Nina summers in CESM, a statistically significant increase in lower tropospheric easterly wind is associated with an increase in North African dust transport over the Atlantic. Barbados dust and Pacific SST variability are only weakly correlated in both observations and CESM, suggesting that other processes are controlling the cross-basin variability of dust. We also use our CESM simulation to show that the relationship between downstream North African dust transport and ENSO fluctuates on multidecadal timescales and is associated with a phase shift in the North Atlantic Oscillation. Our findings indicate that existing observations of dust over the tropical North Atlantic are not extensive enough to completely describe the variability of dust and dust transport, and demonstrate the importance of global models to supplement and interpret observational records.

Eastern China has experienced severe and persistent winter haze episodes in recent years due to intensification of aerosol pollution. In addition to anthropogenic emissions, the winter aerosol pollution over eastern China is associated with unusual meteorological conditions, including weaker wind speeds. Here we show, based on model simulations, that during years with decreased wind speed, large decreases in dust emissions (29%) moderate the wintertime land–sea surface air temperature difference and further decrease winds by −0.06 (±0.05) m s−1 averaged over eastern China. The dust-induced lower winds enhance stagnation of air and account for about 13% of increasing aerosol concentrations over eastern China. Although recent increases in anthropogenic emissions are the main factor causing haze over eastern China, we conclude that natural emissions also exert a significant influence on the increases in wintertime aerosol concentrations, with important implications that need to be taken into account by air quality studies., Anthropogenic aerosol and calm conditions give rise to winter haze episodes in eastern China. Yang et al. show that these weak winds also decrease natural dust emissions, reducing the land–ocean temperature difference and associated winds, enhancing air stagnation and pollution in this region.

In this study, we investigate the sensitivity of simulated shallow cumulus and stratocumulus to selected tunable parameters of Cloud Layers Unified by Binormals (CLUBB) in the single-column version of Community Atmosphere Model version 5 (SCAM5). A quasi-Monte Carlo (QMC) sampling approach is adopted to effectively explore the high-dimensional parameter space and a generalized linear model is adopted to study the responses of simulated cloud fields to tunable parameters. One stratocumulus and two shallow cumulus cases are configured at both coarse and fine vertical resolutions in this study. Our results show that most of the variance in simulated cloud fields can be explained by a small number of tunable parameters. The parameters related to Newtonian and buoyancy-damping terms of total water flux are found to be the most influential parameters for stratocumulus. For shallow cumulus, the most influential parameters are those related to skewness of vertical velocity, reflecting the strong coupling between cloud properties and dynamics in this regime. The influential parameters in the stratocumulus case are sensitive to the vertical resolution while little sensitivity is found for the shallow cumulus cases, as eddy mixing length (or dissipation time scale) plays a more important role and depends more strongly on the vertical resolution in stratocumulus than in shallow convections. The influential parameters remain almost unchanged when the number of tunable parameters increases from 16 to 35. This study improves understanding of the CLUBB behavior associated with parameter uncertainties and provides valuable insights for other high-order turbulence closure schemes.

This study uses a century length preindustrial climate simulation by the Community Earth System Model (CESM 1.0) to explore statistical relationships between dust, clouds, and atmospheric circulation and to suggest a semidirect dynamical mechanism linking subtropical North Atlantic lower tropospheric cloud cover with North African dust transport. The length of the run allows us to account for interannual variability of North African dust emissions and transport in the model. CESM's monthly climatology of both aerosol optical depth and surface dust concentration at Cape Verde and Barbados, respectively, agree well with available observations, as does the aerosol size distribution at Cape Verde. In addition, CESM shows strong seasonal cycles of dust burden and lower tropospheric cloud fraction, with maximum values occurring during boreal summer, when a strong correlation between these two variables exists over the subtropical North Atlantic. Calculations of Estimated Inversion Strength (EIS) and composites of EIS on high and low downstream North African dust months during boreal summer reveal that dust is likely increasing inversion strength over this region due to both solar absorption and reflection. We find no evidence for a microphysical link between dust and lower tropospheric clouds in this region. These results yield new insight over an extensive period of time into the complex relationship between North African dust and North Atlantic lower tropospheric clouds, which has previously been hindered by spatiotemporal constraints of observations. Our findings lay a framework for future analyses using different climate models and submonthly data over regions with different underlying dynamics.

Large-eddy simulations of mixed-phase Arctic clouds by 11 different models are analyzed with the goal of improving understanding and model representation of processes controlling the evolution of these clouds. In a case based on observations from the Indirect and Semi-Direct Aerosol Campaign (ISDAC), it is found that ice number concentration, Ni, exerts significant influence on the cloud structure. Increasing Ni leads to a substantial reduction in liquid water path (LWP), in agreement with earlier studies. In contrast to previous intercomparison studies, all models here use the same ice particle properties (i.e., mass-size, mass-fall speed, and mass-capacitance relationships) and a common radiation parameterization. The constrained setup exposes the importance of ice particle size distributions (PSDs) in influencing cloud evolution. A clear separation in LWP and IWP predicted by models with bin and bulk microphysical treatments is documented and attributed primarily to the assumed shape of ice PSD used in bulk schemes. Compared to the bin schemes that explicitly predict the PSD, schemes assuming exponential ice PSD underestimate ice growth by vapor deposition and overestimate mass-weighted fall speed leading to an underprediction of IWP by a factor of two in the considered case. Sensitivity tests indicate LWP and IWP are much closer to the bin model simulations when a modified shape factor which is similar to that predicted by bin model simulation is used in bulk scheme. These results demonstrate the importance of representation of ice PSD in determining the partitioning of liquid and ice and the longevity of mixed-phase clouds.

Most estimates of the global mean indirect effect of anthropogenic aerosol on the Earth's energy balance are from simulations by global models of the aerosol lifecycle coupled with global models of clouds and the hydrologic cycle. Extremely simple models have been developed for integrated assessment models, but lack the flexibility to distinguish between primary and secondary sources of aerosol. Here a simple but more physically based model expresses the aerosol indirect effect (AIE) using analytic representations of cloud and aerosol distributions and processes. Although the simple model is able to produce estimates of AIEs that are comparable to those from some global aerosol models using the same global mean aerosol properties, the estimates by the simple model are sensitive to preindustrial cloud condensation nuclei concentration, preindustrial accumulation mode radius, width of the accumulation mode, size of primary particles, cloud thickness, primary and secondary anthropogenic emissions, the fraction of the secondary anthropogenic emissions that accumulates on the coarse mode, the fraction of the secondary mass that forms new particles, and the sensitivity of liquid water path to droplet number concentration. Estimates of present-day AIEs as low as 5 W/sq m and as high as 0.3 W/sq m are obtained for plausible sets of parameter values. Estimates are surprisingly linear in emissions. The estimates depend on parameter values in ways that are consistent with results from detailed global aerosol-climate simulation models, which adds to understanding of the dependence on AIE uncertainty on uncertainty in parameter values.

Black carbon aerosol plays a unique and important role in Earth's climate system. Black carbon is a type of carbonaceous material with a unique combination of physical properties. This assessment provides an evaluation of black-carbon climate forcing that is comprehensive in its inclusion of all known and relevant processes and that is quantitative in providing best estimates and uncertainties of the main forcing terms: direct solar absorption; influence on liquid, mixed phase, and ice clouds; and deposition on snow and ice. These effects are calculated with climate models, but when possible, they are evaluated with both microphysical measurements and field observations. Predominant sources are combustion related, namely, fossil fuels for transportation, solid fuels for industrial and residential uses, and open burning of biomass. Total global emissions of black carbon using bottom-up inventory methods are 7500 Gg yr−1 in the year 2000 with an uncertainty range of 2000 to 29000. However, global atmospheric absorption attributable to black carbon is too low in many models and should be increased by a factor of almost 3. After this scaling, the best estimate for the industrial-era (1750 to 2005) direct radiative forcing of atmospheric black carbon is +0.71 W m−2 with 90% uncertainty bounds of (+0.08, +1.27) W m−2. Total direct forcing by all black carbon sources, without subtracting the preindustrial background, is estimated as +0.88 (+0.17, +1.48) W m−2. Direct radiative forcing alone does not capture important rapid adjustment mechanisms. A framework is described and used for quantifying climate forcings, including rapid adjustments. The best estimate of industrial-era climate forcing of black carbon through all forcing mechanisms, including clouds and cryosphere forcing, is +1.1 W m−2 with 90% uncertainty bounds of +0.17 to +2.1 W m−2. Thus, there is a very high probability that black carbon emissions, independent of co-emitted species, have a positive forcing and warm the climate. We estimate that black carbon, with a total climate forcing of +1.1 W m−2, is the second most important human emission in terms of its climate forcing in the present-day atmosphere; only carbon dioxide is estimated to have a greater forcing. Sources that emit black carbon also emit other short-lived species that may either cool or warm climate. Climate forcings from co-emitted species are estimated and used in the framework described herein. When the principal effects of short-lived co-emissions, including cooling agents such as sulfur dioxide, are included in net forcing, energy-related sources (fossil fuel and biofuel) have an industrial-era climate forcing of +0.22 (−0.50 to +1.08) W m−2 during the first year after emission. For a few of these sources, such as diesel engines and possibly residential biofuels, warming is strong enough that eliminating all short-lived emissions from these sources would reduce net climate forcing (i.e., produce cooling). When open burning emissions, which emit high levels of organic matter, are included in the total, the best estimate of net industrial-era climate forcing by all short-lived species from black-carbon-rich sources becomes slightly negative (−0.06 W m−2 with 90% uncertainty bounds of −1.45 to +1.29 W m−2). The uncertainties in net climate forcing from black-carbon-rich sources are substantial, largely due to lack of knowledge about cloud interactions with both black carbon and co-emitted organic carbon. In prioritizing potential black-carbon mitigation actions, non-science factors, such as technical feasibility, costs, policy design, and implementation feasibility play important roles. The major sources of black carbon are presently in different stages with regard to the feasibility for near-term mitigation. This assessment, by evaluating the large number and complexity of the associated physical and radiative processes in black-carbon climate forcing, sets a baseline from which to improve future climate forcing estimates.

Many global aerosol and climate models, including the widely used Community Atmosphere Model version 5 (CAM5), have large biases in predicting aerosols in remote regions such as upper troposphere and high latitudes. In this study, we conduct CAM5 sensitivity simulations to understand the role of key processes associated with aerosol transformation and wet removal affecting the vertical and horizontal long-range transport of aerosols to the remote regions. Improvements are made to processes that are currently not well represented in CAM5, which are guided by surface and aircraft measurements together with results from a multi-scale aerosol-climate model (PNNL-MMF) that explicitly represents convection and aerosol-cloud interactions at cloud-resolving scales. We pay particular attention to black carbon (BC) due to its importance in the Earth system and the availability of measurements. We introduce into CAM5 a new unified scheme for convective transport and aerosol wet removal with explicit aerosol activation above convective cloud base. This new implementation reduces the excessive BC aloft to better simulate observed BC profiles that show decreasing mixing ratios in the mid- to upper-troposphere. After implementing this new unified convective scheme, we examine wet removal of submicron aerosols that occurs primarily through cloud processes. The wet removal depends strongly on the sub-grid scale liquid cloud fraction and the rate of conversion of liquid water to precipitation. These processes lead to very strong wet removal of BC and other aerosols over mid- to high latitudes during winter months. With our improvements, the Arctic BC burden has a10-fold (5-fold) increase in the winter (summer) months, resulting in a much better simulation of the BC seasonal cycle as well. Arctic sulphate and other aerosol species also increase but to a lesser extent. An explicit treatment of BC aging with slower aging assumptions produces an additional 30-fold (5-fold) increase in the Arctic winter (summer) BC burden. This BC aging treatment, however, has minimal effect on other under-predicted species. Interestingly, our modifications to CAM5 that aim at improving prediction of high-latitude and upper tropospheric aerosols also produce much better aerosol optical depth over various other regions globally when compared to multi-year AERONET retrievals. The improved aerosol distributions have impacts on other aspects of CAM5, improving the simulation of global mean liquid water path and cloud forcing.

The objective of this study is to document and evaluate recent changes and updates to the module for aerosols and aerosol–cloud–radiation interactions in the atmospheric module CAM4-Oslo of the core version of the Norwegian Earth System Model (NorESM), NorESM1-M. Particular attention is paid to the role of natural organics, sea salt, and mineral dust in determining the gross aerosol properties as well as the anthropogenic contribution to these properties and the associated direct and indirect radiative forcing. The aerosol module is extended from earlier versions that have been published, and includes life-cycling of sea salt, mineral dust, particulate sulphate, black carbon, and primary and secondary organics. The impacts of most of the numerous changes since previous versions are thoroughly explored by sensitivity experiments. The most important changes are: modified prognostic sea salt emissions; updated treatment of precipitation scavenging and gravitational settling; inclusion of biogenic primary organics and methane sulphonic acid (MSA) from oceans; almost doubled production of land-based biogenic secondary organic aerosols (SOA); and increased ratio of organic matter to organic carbon (OM/OC) for biomass burning aerosols from 1.4 to 2.6. Compared with in situ measurements and remotely sensed data, the new treatments of sea salt and dust aerosols give smaller biases in near-surface mass concentrations and aerosol optical depth than in the earlier model version. The model biases for mass concentrations are approximately unchanged for sulphate and BC. The enhanced levels of modeled OM yield improved overall statistics, even though OM is still underestimated in Europe and overestimated in North America. The global anthropogenic aerosol direct radiative forcing (DRF) at the top of the atmosphere has changed from a small positive value to −0.08 W m−2 in CAM4-Oslo. The sensitivity tests suggest that this change can be attributed to the new treatment of biomass burning aerosols and gravitational settling. Although it has not been a goal in this study, the new DRF estimate is closer both to the median model estimate from the AeroCom intercomparison and the best estimate in IPCC AR4. Estimated DRF at the ground surface has increased by ca. 60%, to −1.89 W m−2. We show that this can be explained by new emission data and omitted mixing of constituents between updrafts and downdrafts in convective clouds. The increased abundance of natural OM and the introduction of a cloud droplet spectral dispersion formulation are the most important contributions to a considerably decreased estimate of the indirect radiative forcing (IndRF). The IndRF is also found to be sensitive to assumptions about the coating of insoluble aerosols by sulphate and OM. The IndRF of −1.2 W m−2, which is closer to the IPCC AR4 estimates than the previous estimate of −1.9 W m−2, has thus been obtained without imposing unrealistic artificial lower bounds on cloud droplet number concentrations.

Two 150-yr preindustrial simulations with and without interactive sea salt emissions from the Community Earth System Model (CESM) are performed to quantify the interactions between sea salt emissions and El Niño–Southern Oscillation (ENSO). Variations in sea salt emissions over the tropical Pacific Ocean are affected by changing wind speed associated with ENSO variability. ENSO-induced interannual variations in sea salt emissions result in decreasing (increasing) aerosol optical depth (AOD) by 0.03 over the equatorial central-eastern (western) Pacific Ocean during El Niño events compared to those during La Niña events. These changes in AOD further increase (decrease) radiative fluxes into the atmosphere by +0.2 (−0.4) W m−2 over the tropical eastern (western) Pacific. Thereby, sea surface temperature increases (decreases) by 0.2–0.4 K over the tropical eastern (western) Pacific Ocean during El Niño compared to La Niña events and enhances ENSO variability by 10%. The increase in ENSO amplitude is a result of systematic heating (cooling) during the warm (cold) phase of ENSO in the eastern Pacific. Interannual variations in sea salt emissions then produce the anomalous ascent (subsidence) over the equatorial eastern (western) Pacific between El Niño and La Niña events, which is a result of heating anomalies. Owing to variations in sea salt emissions, the convective precipitation is enhanced by 0.6–1.2 mm day−1 over the tropical central-eastern Pacific Ocean and weakened by 0.9–1.5 mm day−1 over the Maritime Continent during El Niño compared to La Niña events, enhancing the precipitation variability over the tropical Pacific.

We examined the multi-year mean and variability of dimethyl sulfide (DMS) and its relationship to sulfate aerosols, as well as cloud microphysical and radiative properties. We conducted a 150-year simulation using pre-industrial conditions produced by the Community Earth System Model embedded with a dynamic DMS module. The model simulated the mean spatial distribution of DMS emissions and burden, as well as sulfur budgets associated with DMS, SO2, H2SO4, and sulfate that were generally similar to available observations and inventories for a variety of regions. Changes in simulated sea-to-air DMS emissions and associated atmospheric abundance, along with associated aerosols and cloud and radiative properties, were consistently dominated by the El Nino-Southern Oscillation (ENSO) cycle in the tropical Pacific region. Simulated DMS, aerosols, and clouds showed a weak positive feedback on sea surface temperature. This feedback suggests a link among DMS, aerosols, clouds, and climate on interannual timescales. The variability of DMS emissions associated with ENSO was primarily caused by a higher variation in wind speed during La Nina events. The simulation results also suggest that variations in DMS emissions increase the frequency of La Nina events but do not alter the ENSO variability in terms of the standard deviation of Nino 3 SST anomalies.

Seventeen experts from communities focused on improving understanding and 6 modeling of processes controlling marine low clouds met to discuss future research 7 directions.

A large number of processes are involved in the chain from emissions of aerosol precursor gases and primary particles to impacts on cloud radiative forcing. Those processes are manifest in a number of relationships that can be expressed as factors dlnX/dlnY driving aerosol effects on cloud radiative forcing. These factors include the relationships between cloud condensation nuclei (CCN) concentration and emissions, droplet number and CCN concentration, cloud fraction and droplet number, cloud optical depth and droplet number, and cloud radiative forcing and cloud optical depth. The relationship between cloud optical depth and droplet number can be further decomposed into the sum of two terms involving the relationship of droplet effective radius and cloud liquid water path with droplet number. These relationships can be constrained using observations of recent spatial and temporal variability of these quantities. However, we are most interested in the radiative forcing since the preindustrial era. Because few relevant measurements are available from that era, relationships from recent variability have been assumed to be applicable to the preindustrial to present-day change. Our analysis of Aerosol Comparisons between Observations and Models (AeroCom) model simulations suggests that estimates of relationships from recent variability are poor constraints on relationships from anthropogenic change for some terms, with even the sign of some relationships differing in many regions. Proxies connecting recent spatial/temporal variability to anthropogenic change, or sustained measurements in regions where emissions have changed, are needed to constrain estimates of anthropogenic aerosol impacts on cloud radiative forcing.

Cloud-system resolving simulations with the chemistry version of the Weather Research and Forecasting (WRF-Chem) model are used to quantify the relative impacts of regional anthropogenic and oceanic emissions on changes in aerosol properties, cloud macro- and microphysics, and cloud radiative forcing over the Southeast Pacific (SEP) during the VAMOS Ocean-Cloud-Atmosphere-Land Study Regional Experiment (VOCALS-REx) (15 October–16 November 2008). Two distinct regions are identified. The near-coast polluted region is characterized by low surface precipitation rates, the strong suppression of non-sea-salt particle activation due to sea-salt particles, a predominant albedo effect in aerosol indirect effects, and limited impact of aerosols associated with anthropogenic emissions on clouds. Opposite sensitivities to natural marine and anthropogenic aerosol perturbations are seen in cloud properties (e.g., cloud optical depth and cloud-top and cloud-base heights), precipitation, and the top-of-atmosphere and surface shortwave fluxes over this region. The relatively clean remote region is characterized by large contributions of aerosols from non-regional sources (lateral boundaries) and much stronger drizzle at the surface. Under a scenario of five-fold increase in regional anthropogenic emissions, this relatively clean region shows large cloud responses, for example, a 13% increase in cloud-top height and a 9% increase in albedo in response to a moderate increase (25% of the reference case) in cloud condensation nuclei (CCN) concentration. The reduction of precipitation due to this increase in anthropogenic aerosols more than doubles the aerosol lifetime in the clean marine boundary layer. Therefore, the aerosol impacts on precipitation are amplified by the positive feedback of precipitation on aerosol, which ultimately alters the cloud micro- and macro-physical properties, leading to strong aerosol-cloud-precipitation interactions. The high sensitivity is also related to an increase in cloud-top entrainment rate (by 16% at night) due to the increased anthropogenic aerosols. The simulated aerosol-cloud-precipitation interactions due to the increased anthropogenic aerosols have a stronger diurnal cycle over the clean region compared to the near-coast region with stronger interactions at night. During the day, solar heating results in more frequent decoupling of the cloud and sub-cloud layers, thinner clouds, reduced precipitation, and reduced sensitivity to the increase in anthropogenic emissions. This study shows the importance of natural aerosols in accurately quantifying anthropogenic forcing within a regional modeling framework. The results of this study also imply that the energy balance perturbations from increased anthropogenic emissions are larger in the more susceptible clean environment than in already polluted environment and are larger than possible from the first indirect effect alone.

The authors have decomposed the anthropogenic aerosol radiative forcing into direct contributions from each aerosol species to the planetary energy balance through absorption and scattering of solar radiation, indirect effects of anthropogenic aerosol on solar and infrared radiation through droplet and crystal nucleation on aerosol, and semidirect effects through the influence of solar absorption on the distribution of clouds. A three-mode representation of the aerosol in version 5.1 of the Community Atmosphere Model (CAM5.1) yields global annual mean radiative forcing estimates for each of these forcing mechanisms that are within 0.1 W m−2 of estimates using a more complex seven-mode representation that distinguishes between fresh and aged black carbon and primary organic matter. Simulating fresh black carbon particles separately from internally mixed accumulation mode particles is found to be important only near fossil fuel sources. In addition to the usual large indirect effect on solar radiation, this study finds an unexpectedly large positive longwave indirect effect (because of enhanced cirrus produced by homogenous nucleation of ice crystals on anthropogenic sulfate), small shortwave and longwave semidirect effects, and a small direct effect (because of cancelation and interactions of direct effects of black carbon and sulfate). Differences between the three-mode and seven-mode versions are significantly larger (up to 0.2 W m−2) when the hygroscopicity of primary organic matter is decreased from 0.1 to 0 and transfer of the primary carbonaceous aerosol to the accumulation mode in the seven-mode version requires more hygroscopic material coating the primary particles. Radiative forcing by cloudborne anthropogenic black carbon is only −0.07 W m−2.

Many present-day numerical weather prediction (NWP) models are run at resolutions that permit deep convection. In these models, however, the boundary layer turbulence and boundary layer cloud features are still grossly underresolved. Underresolution is also present in climate models that use a multiscale modeling framework (MMF), in which a convection-permitting model is run in each grid column of a global general circulation model. To better represent boundary layer clouds and turbulence in convection-permitting models, a parameterization was developed that models the joint probability density function (PDF) of vertical velocity, heat, and moisture. Although PDF-based parameterizations are more complex and computationally expensive than many other parameterizations, in principle PDF parameterizations have several advantages. For instance, they ensure consistency of liquid (cloud) water and cloud fraction; they avoid using separate parameterizations for different cloud types such as cumulus and stratocumulus; and they have an appropriate formulation in the “terra incognita” in which updrafts are marginally resolved. In this paper, an implementation of a PDF parameterization is tested to see whether it improves the simulations of a state-of-the-art convection-permitting model. The PDF parameterization used is the Cloud Layers Unified By Binormals (CLUBB) parameterization. The host cloud-resolving model used is the System for Atmospheric Modeling (SAM). SAM is run both with and without CLUBB implemented in it. Simulations of two shallow cumulus (Cu) cases and two shallow stratocumulus (Sc) cases are run in a 3D configuration at 2-, 4-, and 16-km horizontal grid spacings. Including CLUBB in the simulations improves some of the simulated fields—such as vertical velocity variance, horizontal wind fields, cloud water content, and drizzle water content—especially in the two Cu cases. Implementing CLUBB in SAM improves the simulations slightly at 2-km horizontal grid spacing, significantly at 4-km grid spacing, and greatly at 16-km grid spacing. Furthermore, the simulations that include CLUBB exhibit a reduced sensitivity to horizontal grid spacing.

Desert dust plays an important role in the climate system through its impact on Earth¿s radiative budget and its role in the biogeochemical cycle as a source of iron in highnutrient- low-chlorophyll regions. A large degree of diversity exists between the many global models that simulate the dust cycle to estimate its impact on climate. We present the results of a broad intercomparison of a total of 15 global aerosol models within the AeroCom project. Each model is compared to observations focusing on variables responsible for the uncertainties in estimating the direct radiative effect and the dust impact on the biogeochemical cycle, i.e., aerosol optical depth (AOD) and dust deposi10 tion. Additional comparisons to Angstro¨m Exponent (AE), coarse mode AOD and dust surface concentration are included to extend the assessment of model performance. These datasets form a benchmark data set which is proposed for model inspection and future dust model developments. In general, models perform better in simulating climatology of vertically averaged integrated parameters (AOD and AE) in dusty sites 15 than they do with total deposition and surface concentration. Almost all models overestimate deposition fluxes over Europe, the Indian Ocean, the Atlantic Ocean and ice core data. Differences among the models arise when simulating deposition at remote sites with low fluxes over the Pacific and the Southern Atlantic Ocean. This study also highlights important differences in models ability to reproduce the deposition flux over Antarctica. The cause of this discrepancy could not be identified but different dust regimes at each site and issues with data quality should be considered. Models generally simulate better surface concentration at stations downwind of the main sources than at remote ones. Likewise, they simulate better surface concentration at stations affected by Saharan dust than at stations affected by Asian dust. Most models simulate the gradient in AOD and AE between the different dusty regions, however the seasonality and magnitude of both variables is better simulated at African stations than Middle East ones. The models also reproduce the dust transport across the Atlantic in terms of both AOD and AE; they simulate the offshore transport of West Africa throughout the year and limit the transport across the Atlantic to the summer months, yet overestimating the AOD and transporting too fine particles. However, most of the models do not reproduce the southward displacement of the dust cloud during the winter responsible of the transport of dust into South America. Based on the dependency of AOD on aerosol 5 burden and size distribution we use model data bias with respect to AOD and AE and infer on the over/under estimation of the dust emissions. According to this we suggest the emissions in the Sahara be between 792 and 2271 Tg/yr and the one in the Middle East between 376 and 526 Tg/yr., JRC.DDG.H.2-Climate change and air quality