A New Two-Moment Bulk Stratiform Cloud Microphysics Scheme in the Community Atmosphere Model, Version 3 (CAM3). Part I: Description and Numerical Tests

2008 ◽  
Vol 21 (15) ◽  
pp. 3642-3659 ◽  
Author(s):  
Hugh Morrison ◽  
Andrew Gettelman

Abstract A new two-moment stratiform cloud microphysics scheme in a general circulation model is described. Prognostic variables include cloud droplet and cloud ice mass mixing ratios and number concentrations. The scheme treats several microphysical processes, including hydrometeor collection, condensation/evaporation, freezing, melting, and sedimentation. The activation of droplets on aerosol is physically based and coupled to a subgrid vertical velocity. Unique aspects of the scheme, relative to existing two-moment schemes developed for general circulation models, are the diagnostic treatment of rain and snow number concentration and mixing ratio and the explicit treatment of subgrid cloud water variability for calculation of the microphysical process rates. Numerical aspects of the scheme are described in detail using idealized one-dimensional offline tests of the microphysics. Sensitivity of the scheme to time step, vertical resolution, and numerical method for diagnostic precipitation is investigated over a range of conditions. It is found that, in general, two substeps are required for numerical stability and reasonably small time truncation errors using a time step of 20 min; however, substepping is only required for the precipitation microphysical processes rather than the entire scheme. A new numerical approach for the diagnostic rain and snow produces reasonable results compared to a benchmark simulation, especially at low vertical resolution. Part II of this study details results of the scheme in single-column and global simulations, including comparison with observations.

2014 ◽  
Vol 7 (6) ◽  
pp. 7559-7573
Author(s):  
S. Watanabe ◽  
K. Sato ◽  
Y. Kawatani ◽  
M. Takahashi

Abstract. The dependence of the gravity wave spectra of energy and momentum flux on the horizontal resolution and time step of atmospheric general circulation models (AGCMs) has been thoroughly investigated in the past. In contrast, much less attention has been given to the dependence of these gravity wave parameters on models' vertical resolutions. The present study demonstrates the dependence of gravity wave momentum flux in the stratosphere and mesosphere on the model's vertical resolution, which is evaluated using an AGCM with a horizontal resolution of about 0.56°. We performed a series of sensitivity test simulations changing only the model's vertical resolution above a height of 8 km, and found that inertial gravity waves with short vertical wavelengths simulated at higher vertical resolutions likely play an important role in determining the gravity wave momentum flux in the stratosphere and mesosphere.


Author(s):  
Daisuke Matsuoka ◽  
Fumiaki Araki ◽  
Hideharu Sasaki

Numerical study of ocean eddies has been carried out by using high-resolution ocean general circulation models. In order to understand ocean eddies from the large volume data produced by simulations, visualizing only eddy distribution at each time step is insufficient; time-variations in eddy events and phenomena must also be considered. However, existing methods cannot precisely find and track eddy events such as amalgamation and bifurcation. In this study, we propose an original approach for eddy detection, tracking, and event visualization based on an eddy classification system. The proposed method detects streams and currents as well as eddies, and it classifies discovered eddies into several categories using the additional stream and current information. By tracking how the classified eddies vary over time, detecting events such as eddy amalgamation and bifurcation as well as the interaction between eddies and ocean currents becomes achievable. We adopt the proposed method for two ocean areas in which strong ocean currents exist as case studies. We visualize the detected eddies and events in a time series of images, allowing us to acquire an intuitive understanding of a region of interest concealed in a high-resolution data set. Furthermore, our proposed method succeeded in clarifying the occurrence place and seasonality of each type of eddy event.


2014 ◽  
Vol 27 (12) ◽  
pp. 4391-4402 ◽  
Author(s):  
Dorian S. Abbot

Abstract Recent general circulation model (GCM) simulations have challenged the idea that a snowball Earth would be nearly entirely cloudless. This is important because clouds would provide a strong warming to a high-albedo snowball Earth. GCM results suggest that clouds could lower the threshold CO2 needed to deglaciate a snowball by a factor of 10–100, enough to allow consistency with geochemical data. Here a cloud-resolving model is used to investigate cloud and convection behavior in a snowball Earth climate. The model produces convection that extends vertically to a similar temperature as modern tropical convection. This convection produces clouds that resemble stratocumulus clouds under an inversion on modern Earth, which slowly dissipate by sedimentation of cloud ice. There is enough cloud ice for the clouds to be optically thick in the longwave, and the resulting cloud radiative forcing is similar to that produced in GCMs run in snowball conditions. This result is robust to large changes in the cloud microphysics scheme because the cloud longwave forcing, which dominates the total forcing, is relatively insensitive to cloud amount and particle size. The cloud-resolving model results are therefore consistent with the idea that clouds would provide a large warming to a snowball Earth, helping to allow snowball deglaciation.


2010 ◽  
Vol 10 (16) ◽  
pp. 8037-8064 ◽  
Author(s):  
M. Salzmann ◽  
Y. Ming ◽  
J.-C. Golaz ◽  
P. A. Ginoux ◽  
H. Morrison ◽  
...  

Abstract. A new stratiform cloud scheme including a two-moment bulk microphysics module, a cloud cover parameterization allowing ice supersaturation, and an ice nucleation parameterization has been implemented into the recently developed GFDL AM3 general circulation model (GCM) as part of an effort to treat aerosol-cloud-radiation interactions more realistically. Unlike the original scheme, the new scheme facilitates the study of cloud-ice-aerosol interactions via influences of dust and sulfate on ice nucleation. While liquid and cloud ice water path associated with stratiform clouds are similar for the new and the original scheme, column integrated droplet numbers and global frequency distributions (PDFs) of droplet effective radii differ significantly. This difference is in part due to a difference in the implementation of the Wegener-Bergeron-Findeisen (WBF) mechanism, which leads to a larger contribution from super-cooled droplets in the original scheme. Clouds are more likely to be either completely glaciated or liquid due to the WBF mechanism in the new scheme. Super-saturations over ice simulated with the new scheme are in qualitative agreement with observations, and PDFs of ice numbers and effective radii appear reasonable in the light of observations. Especially, the temperature dependence of ice numbers qualitatively agrees with in-situ observations. The global average long-wave cloud forcing decreases in comparison to the original scheme as expected when super-saturation over ice is allowed. Anthropogenic aerosols lead to a larger decrease in short-wave absorption (SWABS) in the new model setup, but outgoing long-wave radiation (OLR) decreases as well, so that the net effect of including anthropogenic aerosols on the net radiation at the top of the atmosphere (netradTOA = SWABS-OLR) is of similar magnitude for the new and the original scheme.


2014 ◽  
Vol 7 (2) ◽  
pp. 2173-2216
Author(s):  
H. Wan ◽  
P. J. Rasch ◽  
K. Zhang ◽  
Y. Qian ◽  
H. Yan ◽  
...  

Abstract. This paper explores the feasibility of an experimentation strategy for investigating sensitivities in fast components of atmospheric general circulation models. The basic idea is to replace the traditional serial-in-time long-term climate integrations by representative ensembles of shorter simulations. The key advantage of the proposed method lies in its efficiency: since fewer days of simulation are needed, the computational cost is less, and because individual realizations are independent and can be integrated simultaneously, the new dimension of parallelism can dramatically reduce the turnaround time in benchmark tests, sensitivities studies, and model tuning exercises. The strategy is not appropriate for exploring sensitivity of all model features, but it is very effective in many situations. Two examples are presented using the Community Atmosphere Model version 5. The first example demonstrates that the method is capable of characterizing the model cloud and precipitation sensitivity to time step length. A nudging technique is also applied to an additional set of simulations to help understand the contribution of physics-dynamics interaction to the detected time step sensitivity. In the second example, multiple empirical parameters related to cloud microphysics and aerosol lifecycle are perturbed simultaneously in order to explore which parameters have the largest impact on the simulated global mean top-of-atmosphere radiation balance. Results show that in both examples, short ensembles are able to correctly reproduce the main signals of model sensitivities revealed by traditional long-term climate simulations for fast processes in the climate system. The efficiency of the ensemble method makes it particularly useful for the development of high-resolution, costly and complex climate models.


2008 ◽  
Vol 8 (11) ◽  
pp. 2949-2963 ◽  
Author(s):  
R. Posselt ◽  
U. Lohmann

Abstract. Prognostic equations for the rain mass mixing ratio and the rain drop number concentration are introduced into the large-scale cloud microphysics parameterization of the ECHAM5 general circulation model (ECHAM5-PROG). To this end, a rain flux from one level to the next with the appropriate fall speed is introduced. This maintains rain water in the atmosphere to be available for the next time step. Rain formation in ECHAM5-PROG is, therefore, less dependent on the autoconversion rate than the standard ECHAM5 but shifts the emphasis towards the accretion rates in accordance with observations. ECHAM5-PROG is tested and evaluated with Single Column Model (SCM) simulations for two cases: the marine stratocumulus study EPIC (October 2001) and the continental mid-latitude ARM Cloud IOP (shallow frontal cloud case – March 2000). In case of heavy precipitation events, the prognostic equations for rain hardly affect the amount and timing of precipitation at the surface in different SCM simulations because heavy rain depends mainly on the large-scale forcing. In case of thin, drizzling clouds (i.e., stratocumulus), surface precipitation is sensitive to the number of sub-time steps used in the prognostic rain scheme. Cloud microphysical quantities, such as cloud liquid and rain water within the atmosphere, are sensitive to the number of sub-time steps in both considered cases. This results from the decreasing autoconversion rate and increasing accretion rate.


2010 ◽  
Vol 10 (3) ◽  
pp. 6375-6446 ◽  
Author(s):  
M. Salzmann ◽  
Y. Ming ◽  
J.-C. Golaz ◽  
P. A. Ginoux ◽  
H. Morrison ◽  
...  

Abstract. A new stratiform cloud scheme including a two-moment bulk microphysics module, a cloud cover parameterization allowing ice supersaturation, and an ice nucleation parameterization has been implemented into the recently developed GFDL AM3 general circulation model (GCM) as part of an effort to treat aerosol-cloud-radiation interactions more realistically. Unlike the original scheme, the new scheme facilitates the study of cloud-ice-aerosol interactions via influences of dust and sulfate on ice nucleation. While liquid and cloud ice water path associated with stratiform clouds are similar for the new and the original scheme, column integrated droplet numbers and global frequency distributions (PDFs) of droplet effective radii differ significantly. This difference is in part due to a difference in the implementation of the Wegener-Bergeron-Findeisen (WBF) mechanism, which leads to a larger contribution from super-cooled droplets in the original scheme. Clouds are more likely to be either completely glaciated or liquid due to the WBF mechanism in the new scheme. Super-saturations over ice simulated with the new scheme are in qualitative agreement with observations, and PDFs of ice numbers and effective radii appear reasonable in the light of observations. Especially, the temperature dependence of ice numbers qualitatively agrees with in-situ observations. The global average long-wave cloud forcing decreases in comparison to the original scheme as expected when super-saturation over ice is allowed. Anthropogenic aerosols lead to a larger decrease in short-wave absorption (SWABS) in the new model setup, but outgoing long-wave radiation (OLR) decreases as well, so that the net effect of including anthropogenic aerosols on the net radiation at the top of the atmosphere (netradTOA = SWABS-OLR) is of similar magnitude for the new and the original scheme.


2018 ◽  
Vol 11 (4) ◽  
pp. 1557-1576 ◽  
Author(s):  
Remo Dietlicher ◽  
David Neubauer ◽  
Ulrike Lohmann

Abstract. A new scheme for stratiform cloud microphysics has been implemented in the ECHAM6-HAM2 general circulation model. It features a widely used description of cloud water with two categories for cloud droplets and raindrops. The unique aspect of the new scheme is the break with the traditional approach to describe cloud ice analogously. Here we parameterize cloud ice by a single category that predicts bulk particle properties (P3). This method has already been applied in a regional model and most recently also in the Community Atmosphere Model 5 (CAM5). A single cloud ice category does not rely on heuristic conversion rates from one category to another. Therefore, it is conceptually easier and closer to first principles. This work shows that a single category is a viable approach to describe cloud ice in climate models. Prognostic representation of sedimentation is achieved by a nested approach for sub-stepping the cloud microphysics scheme. This yields good results in terms of accuracy and performance as compared to simulations with high temporal resolution. Furthermore, the new scheme allows for a competition between various cloud processes and is thus able to unbiasedly represent the ice formation pathway from nucleation to growth by vapor deposition and collisions to sedimentation. Specific aspects of the P3 method are evaluated. We could not produce a purely stratiform cloud where rime growth dominates growth by vapor deposition and conclude that the lack of appropriate conditions renders the prognostic parameters associated with the rime properties unnecessary. Limitations inherent in a single category are examined.


2016 ◽  
Author(s):  
Kristofer Döös ◽  
Bror Jönsson ◽  
Joakim Kjellsson

Abstract. Two different trajectory schemes for oceanic and atmospheric general circulation models are compared in two different experiments. The theories of the two trajectory schemes are presented showing the differential equations they solve and why they are mass conserving. One scheme assumes that the velocity fields are stationary for a limited period of time and solves the trajectory path from a differential equation only as a function of space, i.e. "stepwise stationary". The second scheme uses a continuous linear interpolation of the fields in time and solves the trajectory path from a differential equation as a function of both space and time, i.e. "time-dependent". A special case of the "stepwise-stationary" scheme, when velocities are assumed constant between GCM outputs, is also considered, named "fixed GCM time step". The trajectory schemes are tested "off-line", i.e. using the already integrated and stored velocity fields from a GCM. The first comparison of the schemes uses trajectories calculated using the velocity fields from an eddy-resolving ocean general circulation model in the Agulhas region. The second comparison uses trajectories calculated using the wind fields from an atmospheric reanalysis. The study shows that using the "time-dependent" scheme over the "stepwise-stationary" scheme greatly improves accuracy with only a small increase in computational time. It is also found that with decreasing time steps the "stepwise-stationary" scheme becomes more accurate but at increased computational cost. The "time-dependent" scheme is therefore preferred over the "stepwise-stationary" scheme. However, when averaging over large ensembles of trajectories the two schemes are comparable, as intrinsic variability dominates over numerical errors. The "fixed GCM time step" is found to be less accurate than the "stepwise-stationary" scheme, even when considering averages over large ensembles.


2008 ◽  
Vol 21 (15) ◽  
pp. 3660-3679 ◽  
Author(s):  
A. Gettelman ◽  
H. Morrison ◽  
S. J. Ghan

Abstract The global performance of a new two-moment cloud microphysics scheme for a general circulation model (GCM) is presented and evaluated relative to observations. The scheme produces reasonable representations of cloud particle size and number concentration when compared to observations, and it represents expected and observed spatial variations in cloud microphysical quantities. The scheme has smaller particles and higher number concentrations over land than the standard bulk microphysics in the GCM and is able to balance the top-of-atmosphere radiation budget with 60% the liquid water of the standard scheme, in better agreement with retrieved values. The new scheme diagnostically treats both the mixing ratio and number concentration of rain and snow, and it is therefore able to differentiate the two key regimes, consisting of drizzle in shallow, warm clouds and larger rain drops in deeper cloud systems. The modeled rain and snow size distributions are consistent with observations.


Sign in / Sign up

Export Citation Format

Share Document