The model used here is a coupled ocean-atmosphere model that consists of the CCSR/NIES atmospheric GCM, the CCSR ocean GCM, a thermodynamic sea-ice model, and a river routing model (Abe-Ouchi et al., 1996). The spatial resolution is T21 spectral truncation (roughly 5.6 degrees latitude/longitude) and 20 vertical levels for the atmospheric part, and roughly 2.8 degrees horizontal grid and 17 vertical levels for the oceanic part. Flux adjustment for atmosphere-ocean heat and water exchange is applied to prevent a drift of the modelled climate. The atmospheric model adopts a radiation scheme based on the k-distribution, two-stream discrete ordinate method (DOM) ...