Inclusion of ice microphysics in the NCAR community atmospheric model version 3 (CAM3)

A prognostic equation for ice crystal number concentration together with an ice nucleation scheme are implemented in the National Center for Atmospheric Research (NCAR) Community Atmospheric Model version 3 (CAM3) with the aim of studying the indirect effect of aerosols on cold clouds. The effective radius of ice crystals, which is used in the radiation and gravitational settlement calculations, is now calculated from model-predicted mass and number of ice crystals rather than diagnosed as a function of temperature. A water vapor deposition scheme is added to replace the condensation and evaporation (C-E) in the standard CAM3 for ice clouds. The repartitioning of total water into liquid and ice in mixed-phase clouds as a function of temperature is removed, and ice supersaturation is allowed. The predicted ice water content in the modified CAM3 is in better agreement with the Aura Microwave Limb Sounder (MLS) data than that in the standard CAM3. The cirrus cloud fraction near the tropical tropopause, which is underestimated in the standard CAM3 as revealed through comparison with the Stratospheric Aerosol and Gas Experiment II (SAGE II) data, is increased by 20%-30%, and the cold temperature bias there is reduced by 1-2 K. However, an increase in the cloud fraction in polar regions makes the underestimation ( by similar to 20 W m(-2)) of downwelling shortwave radiation in the standard CAM3 even worse. A sensitivity test reducing the threshold relative humidity with respect to ice (RHi) for heterogeneous ice nucleation from 120% to 105% (representing nearly perfect ice nuclei) increases the global cloud cover by 1.4%, temperature near the tropical tropopause by 4-5 K, and water vapor in the stratosphere by 50%-80%.