The fourth version of the Community Climate System Model (CCSM4) was recently completed and released to the climate community. This paper describes developments to all CCSM components, and documents fully coupled preindustrial control runs compared to the previous version, CCSM3. Using the standard atmosphere and land resolution of 18 results in the sea surface temperature biases in the major upwelling regions being comparable to the 1.48-resolution CCSM3. Two changes to the deep convection scheme in the atmosphere component result in CCSM4 producing El Niñ o-Southern Oscillation variability with a much more realistic frequency distribution than in CCSM3, although the amplitude is too large compared to observations. These changes also improve the Madden-Julian oscillation and the frequency distribution of tropical precipitation. A new overflow parameterization in the ocean component leads to an improved simulation of the Gulf Stream path and the North Atlantic Ocean meridional overturning circulation. Changes to the CCSM4 land component lead to a much improved annual cycle of water storage, especially in the tropics. The CCSM4 sea ice component uses much more realistic albedos than CCSM3, and for several reasons the Arctic sea ice concentration is improved in CCSM4. An ensemble of twentieth-century simulations produces a good match to the observed September Arctic sea ice extent from 1979 to 2005. The CCSM4 ensemble mean increase in globally averaged surface temperature between 1850 and 2005 is larger than the observed increase by about 0.48C. This is consistent with the fact that CCSM4 does not include a representation of the indirect effects of aerosols, although other factors may come into play. The CCSM4 still has significant biases, such as the mean precipitation distribution in the tropical Pacific Ocean, too much low cloud in the Arctic, and the latitudinal distributions of shortwave and longwave cloud forcings.

[1] We apply our Snow, Ice, and Aerosol Radiative (SNICAR) model, coupled to a general circulation model with prognostic carbon aerosol transport, to improve understanding of climate forcing and response from black carbon (BC) in snow. Building on two previous studies, we account for interannually varying biomass burning BC emissions, snow aging, and aerosol scavenging by snow meltwater. We assess uncertainty in forcing estimates from these factors, as well as BC optical properties and snow cover fraction. BC emissions are the largest source of uncertainty, followed by snow aging. The rate of snow aging determines snowpack effective radius (r e ), which directly controls snow reflectance and the magnitude of albedo change caused by BC. For a reasonable r e range, reflectance reduction from BC varies threefold. Inefficient meltwater scavenging keeps hydrophobic impurities near the surface during melt and enhances forcing. Applying biomass burning BC emission inventories for a strong (1998) and weak (2001) boreal fire year, we estimate global annual mean BC/snow surface radiative forcing from all sources (fossil fuel, biofuel, and biomass burning) of +0.054 (0.007-0.13) and +0.049 (0.007-0.12) W m À2 , respectively. Snow forcing from only fossil fuel + biofuel sources is +0.043 W m À2 (forcing from only fossil fuels is +0.033 W m À2 ), suggesting that the anthropogenic contribution to total forcing is at least 80%. The 1998 global land and sea-ice snowpack absorbed 0.60 and 0.23 W m À2 , respectively, because of direct BC/snow forcing. The forcing is maximum coincidentally with snowmelt onset, triggering strong snow-albedo feedback in local springtime. Consequently, the ''efficacy'' of BC/snow forcing is more than three times greater than forcing by CO 2 . The 1998 and 2001 land snowmelt rates north of 50°N are 28% and 19% greater in the month preceding maximum melt of control simulations without BC in snow. With climate feedbacks, global annual mean 2-meter air temperature warms 0.15 and 0.10°C, when BC is included in snow, whereas annual arctic warming is 1.61 and 0.50°C. Stronger highlatitude climate response in 1998 than 2001 is at least partially caused by boreal fires, which account for nearly all of the 35% biomass burning contribution to 1998 arctic forcing. Efficacy was anomalously large in this experiment, however, and more research is required to elucidate the role of boreal fires, which we suggest have maximum arctic BC/snow forcing potential during April-June. Model BC concentrations in snow agree reasonably well (r = 0.78) with a set of 23 observations from various locations, spanning nearly 4 orders of magnitude. We predict concentrations in excess of 1000 ng g À1 for snow in northeast China, enough to lower snow albedo by more than 0.13. The greatest instantaneous forcing is over the Tibetan Plateau, exceeding 20 W m À2 in some places during spring. These results indicate that snow darkening is an important component of carbon aerosol climate forcing.Citation: Flanner, M. G., C. S. Zender, J. T. Ra...

Abstract.A modal aerosol module (MAM) has been developed for the Community Atmosphere Model version 5 (CAM5), the atmospheric component of the Community Earth System Model version 1 (CESM1). MAM is capable of simulating the aerosol size distribution and both internal and external mixing between aerosol components, treating numerous complicated aerosol processes and aerosol physical, chemical and optical properties in a physically-based manner. Two MAM versions were developed: a more complete version with seven lognormal modes (MAM7), and a version with three lognormal modes (MAM3) for the purpose of long-term (decades to centuries) simulations. In this paper a description and evaluation of the aerosol module and its two representations are provided. Sensitivity of the aerosol lifecycle to simplifications in the representation of aerosol is discussed.Simulated sulfate and secondary organic aerosol (SOA) mass concentrations are remarkably similar between MAM3 and MAM7. Differences in primary organic matter (POM) and black carbon (BC) concentrations between MAM3 and MAM7 are also small (mostly within 10 %). The mineral dust global burden differs by 10 % and sea salt burden by 30-40 % between MAM3 and MAM7, mainly due to the different size ranges for dust and sea salt modes and different standard deviations of the log-normal size distribution for sea salt modes between MAM3 and MAM7. The model is able to qualitatively capture the observed geographical and temporal variations of aerosol mass and number concentrations, size distributions, and aerosol optical properties. However, there are noticeable biases; e.g., simulated BC concentrations are significantly lower than measurements in the Arctic. There is a low bias in modeled aerosol optical depth on the global scale, especially in the developing countries. These biases in aerosol simulations clearly indicate the need for improvements of aerosol processes (e.g., emission fluxes of anthropogenic aerosols and precursor gases in developing countries, boundary layer nucleation) and properties (e.g., primary aerosol emission size, POM hygroscopicity). In addition, the critical role of cloud properties (e.g., liquid water content, cloud fraction) responsible for the wet scavenging of aerosol is highlighted.

No abstract

[1] We have developed a global three-dimensional chemical transport model called Model of Ozone and Related Chemical Tracers (MOZART), version 2. This model, which will be made available to the community, is built on the framework of the National Center for Atmospheric Research (NCAR) Model of Atmospheric Transport and Chemistry (MATCH) and can easily be driven with various meteorological inputs and model resolutions. In this work, we describe the standard configuration of the model, in which the model is driven by meteorological inputs every 3 hours from the middle atmosphere version of the NCAR Community Climate Model (MACCM3) and uses a 20-min time step and a horizontal resolution of 2.8°latitude  2.8°longitude with 34 vertical levels extending up to approximately 40 km. The model includes a detailed chemistry scheme for tropospheric ozone, nitrogen oxides, and hydrocarbon chemistry, with 63 chemical species. Tracer advection is performed using a flux-form semi-Lagrangian scheme with a pressure fixer. Subgrid-scale convective and boundary layer parameterizations are included in the model. Surface emissions include sources from fossil fuel combustion, biofuel and biomass burning, biogenic and soil emissions, and oceanic emissions. Parameterizations of dry and wet deposition are included. Stratospheric concentrations of several long-lived species (including ozone) are constrained by relaxation toward climatological values. The distribution of tropospheric ozone is well simulated in the model, including seasonality and horizontal and vertical gradients. However, the model tends to overestimate ozone near the tropopause at high northern latitudes. Concentrations of nitrogen oxides (NO x ) and nitric acid (HNO 3 ) agree well with observed values, but peroxyacetylnitrate (PAN) is overestimated by the model in the upper troposphere at several locations. Carbon monoxide (CO) is simulated well at most locations, but the seasonal cycle is underestimated at some sites in the Northern Hemisphere. We find that in situ photochemical production and loss dominate the tropospheric ozone budget, over input from the stratosphere and dry deposition. Approximately 75% of the tropospheric production and loss of ozone occurs within the tropics, with large net production in the tropical upper troposphere. Tropospheric production and loss of ozone are three to four times greater in the northern extratropics than the southern extratropics. The global sources of CO consist of photochemical production (55%) and direct emissions (45%). The tropics dominate the chemistry of CO, accounting for about 75% of the tropospheric production and loss. The global budgets of tropospheric ozone and CO are generally consistent with the range found in recent studies. The lifetime of methane (9.5 years) and methylchloroform (5.7 years) versus oxidation by tropospheric hydroxyl radical (OH), two useful measures of the global abundance of OH, agree well with recent estimates. Concentrations of nonmethane hydrocarbons and oxygenated intermediates (carbony...

We report on the AeroCom Phase II direct aerosol effect (DAE) experiment where 16 detailed global aerosol models have been used to simulate the changes in the aerosol distribution over the industrial era. All 16 models have estimated the radiative forcing (RF) of the anthropogenic DAE, and have taken into account anthropogenic sulphate, black carbon (BC) and organic aerosols (OA) from fossil fuel, biofuel, and biomass burning emissions. In addition several models have simulated the DAE of anthropogenic nitrate and anthropogenic influenced secondary organic aerosols (SOA). The model simulated all-sky RF of the DAE from total anthropogenic aerosols has a range from −0.58 to −0.02 Wm−2, with a mean of −0.27 Wm−2 for the 16 models. Several models did not include nitrate or SOA and modifying the estimate by accounting for this with information from the other AeroCom models reduces the range and slightly strengthens the mean. Modifying the model estimates for missing aerosol components and for the time period 1750 to 2010 results in a mean RF for the DAE of −0.35 Wm−2. Compared to AeroCom Phase I (Schulz et al., 2006) we find very similar spreads in both total DAE and aerosol component RF. However, the RF of the total DAE is stronger negative and RF from BC from fossil fuel and biofuel emissions are stronger positive in the present study than in the previous AeroCom study. We find a tendency for models having a strong (positive) BC RF to also have strong (negative) sulphate or OA RF. This relationship leads to smaller uncertainty in the total RF of the DAE compared to the RF of the sum of the individual aerosol components. The spread in results for the individual aerosol components is substantial, and can be divided into diversities in burden, mass extinction coefficient (MEC), and normalized RF with respect to AOD. We find that these three factors give similar contributions to the spread in results

A new version of the Community Atmosphere Model (CAM) has been developed and released to the climate community. CAM Version 3 (CAM3) is an atmospheric general circulation model that includes the Community Land Model (CLM3), an optional slab ocean model, and a thermodynamic sea ice model. The dynamics and physics in CAM3 have been changed substantially compared to implementations in previous versions. CAM3 includes options for Eulerian spectral, semi-Lagrangian, and finite-volume formulations of the dynamical equations. It supports coupled simulations using either finite-volume or Eulerian dynamics through an explicit set of adjustable parameters governing the model time step, cloud parameterizations, and condensation processes. The model includes major modifications to the parameterizations of moist processes, radiation processes, and aerosols. These changes have improved several aspects of the simulated climate, including more realistic tropical tropopause temperatures, boreal winter land surface temperatures, surface insolation, and clear-sky surface radiation in polar regions. The variation of cloud radiative forcing during ENSO events exhibits much better agreement with satellite observations. Despite these improvements, several systematic biases reduce the fidelity of the simulations. These biases include underestimation of tropical variability, errors in tropical oceanic surface fluxes, underestimation of implied ocean heat transport in the Southern Hemisphere, excessive surface stress in the storm tracks, and offsets in the 500-mb height field and the Aleutian low.

No abstract