1. Introduction
An oceanic thermal dome around the Republic of Seychelles in the southwestern Indian Ocean (Reverdin and Fieux 1987; Woodberry et al. 1989; McCreary et al. 1993) is referred to simply as Seychelles Dome (SD; Yokoi et al. 2008). Because the thermocline is very shallow over the SD, the region is known to be a center of strong ocean–atmosphere interaction. In fact, the sea surface temperature anomaly (SSTA) over the southwestern Indian Ocean including the SD region, which is induced by oceanic Rossby waves (Murtugudde et al. 2000; Xie et al. 2002; Rao et al. 2002; Huang and Kinter 2002; Rao and Behera 2005), influences the onset of the Indian summer monsoon (Annamalai et al. 2005; Izumo et al. 2008). When the dome is weak and the overlying SSTA is positive, a northward migration of the intertropical convergence zone (ITCZ) toward the Indian subcontinent is delayed. In addition, the SD is located close to the western pole of the Indian Ocean dipole (IOD; Saji et al. 1999; Webster et al. 1999; Yamagata et al. 2004; Chang et al. 2006), which plays an important role in the variability of the East African short rains (Black et al. 2003; Behera et al. 2005) and world climate (Guan and Yamagata 2003; Saji and Yamagata 2003). The SST variations in the dome region may even influence the Pacific–North America (PNA) pattern during El Niño years through atmospheric teleconnections (Annamalai et al. 2007b). Furthermore, the number of tropical cyclones is closely related to the strength of the SD (Xie et al. 2002). Thus, a successful simulation of the SD in coupled GCMs (CGCMs) is crucial to better climate predictability.
However, no study to date has compared the ability of different models to simulate the oceanic dome, even though it is of great importance to climate and ecology. Therefore, we believe that such an intercomparison study will enhance our understanding of the SD and the associated air–sea interaction and may contribute to the improvement of these coupled GCMs. In fact, it is a very challenging task for CGCMs to simulate the seasonal variation of the SD (Yokoi et al. 2008; Hermes and Reason 2008). Yokoi et al. (2008) demonstrated that it is dominated by a semiannual cycle due to the local Ekman upwelling. This semiannual nature is explained by unique contributions of two components of the Ekman pumping: one term proportional to the planetary beta and the zonal wind stress and the other term proportional to the wind stress curl. The former is determined by the seasonal change in the zonal component of the wind stress vector above the SD associated with the Indian monsoon and causes downwelling (upwelling) during boreal summer (boreal winter). On the other hand, the latter is due to the meridional gradient of the zonal wind stress and causes almost steady upwelling for 5 months from June to October even though the zonal wind stress itself varies significantly during this period. Thus, the semiannual variation, despite being remarkable, is due to a subtle balance between the seasonal variation of the zonal wind stress and the wind stress curl.
In this study, we have examined the ability of coupled GCMs to simulate the SD and its seasonal variation. The paper is organized as follows: A brief description of 22 coupled GCMs submitted to the World Climate Research Programme (WCRP) Coupled Model Intercomparison Project, phase 3, (CMIP3) multimodel dataset is given in section 2. In section 3, we analyze the SD simulated in each model and discuss causes of model biases. A summary and discussion are presented in the final section.
2. Description of CMIP3 models
We here examine ocean temperature, precipitation, and wind stress data from the “twentieth-century climate in coupled models” (20c3m) control runs of 22 coupled GCMs submitted to the CMIP3. Among 23 models submitted at the time of the present analysis, we could extract both ocean temperature and wind stress data for 22 models, listed in Table 1 (details of these models are described at http://www-pcmdi.llnl.gov/ipcc/model_documentation/ipcc_model_documentation.php). Using the available dataset, we have calculated the monthly mean climatology for the last 50 yr of the twentieth century to investigate the simulated seasonal cycle of the SD. For the Hadley Centre Global Environmental Model, version 1 (HadGEM1), only 45 yr of oceanic temperature data were available.
The model outputs are compared with various observational datasets including the World Ocean Atlas 1998 (WOA98; Antonov et al. 1998), the National Centers for Environmental Prediction–National Center for Atmospheric Research (NCEP–NCAR) reanalysis data (Kalnay et al. 1996), and the Climate Prediction Center (CPC) Merged Analysis of Precipitation (CMAP) data (Xie and Arkin 1996). The WOA98 provides a global hydrographic temperature data with a uniform horizontal resolution of 1° and 24 vertical levels. The NCEP–NCAR reanalysis data are used for the wind stress data.
3. Seasonal variation
Figure 1 shows the distribution of the annual mean depth of the 20°C isotherm (D20) in the tropical Indian Ocean. Except for four models [CGCM3.1 (T47), CGCM3.1 (T63), the Goddard Institute for Space Studies Atmosphere–Ocean Model (GISS-AOM), and the GISS Model E-R (GISS-ER)], all other models are successful in reproducing the SD; the D20 of GISS-AOM and GISS-ER is too deep (their D20 is deeper than 200 m), whereas the upwelling of CGCM3.1 (T47) and CGCM3.1 (T63) appears as coastal upwelling in the eastern Indian Ocean along Sumatra (Saji et al. 2006). However, the core of the dome in many models is located too far (more than 10° in longitude) to the east [the Commonwealth Scientific and Industrial Research Organisation Mark, version 3.0 (CSIRO-Mk3.0); the Geophysical Fluid Dynamics Laboratory Climate Model, version 2.0 (GFDL-CM2.0); the Model for Interdisciplinary Research on Climate 3.2 (MIROC3.2), in both the high-resolution (hires) and medium-resolution version (medres) versions; and the Community Climate System Model, version 3 (CCSM3)] or to the west (the Institute of Numerical Mathematics Coupled Model, version 3.0 (INM-CM3.0) and ECHAM and the global Hamburg Ocean Primitive Equation model (ECHO-G)] when compared with that of the WOA98, whose center is located at 62°E, 7.5°S (Fig. 2). Here, we have defined the core by the location of minimum D20 between 12° and 5°S in the Indian Ocean, excluding coastal regions. Only six models [the Centre National de Recherches Météorologiques Coupled Global Climate Model, version 3 (CNRM-CM3), the GFDL Climate Model, version 2.1 (GFDL-CM2.1), GISS Model E-H (GISS-EH), the Istituto Nazionale di Geofisica e Vulcanologia SXG (INGV-SXG), the Meteorological Research Institute Coupled General Circulation Model, version 2.3.2a (MRI-CGCM2.3.2), and HadGEM1] simulate the center within 5° of the observed position. In addition, the dome is too deep in MIROC3.2 (hires) and the third climate configuration of the Met Office Unified Model (HadCM3), whereas it is too shallow in the Parallel Climate Model (PCM).
Next, the seasonal variation of the D20 within the simulated SD and the amplitudes of the annual and semiannual harmonics are shown (Figs. 3 and 4). Here, we have calculated average D20 in a 16° (longitude) × 4° (latitude) box; it is shifted so that it can cover the simulated SD in each model. It turns out that most models except for CNRM-CM3, L’Institut Pierre-Simon Laplace Coupled Model, version 4 (IPSL-CM4), HadGEM1, and INM-CM3.0 underestimate the amplitude of the semiannual harmonic (Fig. 4b). In contrast, most models except for ECHO-G, the Flexible Global Ocean–Atmosphere–Land System Model, gridpoint version 1.0 (FGOALS-g1.0), and GISS-EH overestimate the amplitude of the annual harmonic (Fig. 4a). As a result, the ratio between the semiannual harmonic and the annual harmonic is much smaller in most coupled GCMs when compared with that of the WOA98 (Fig. 4c), suggesting that most CMIP3 models have difficulties in simulating the realistic seasonal variation of the SD.
To identify the origin of these model biases, we have examined the seasonal variation of the Ekman pumping in the SD region because it plays the dominant role in the seasonal variation of the SD (Yokoi et al. 2008). As expected, the ratio between the semiannual harmonic and the annual harmonic for the local Ekman pumping and that for D20 in the dome region have good correspondence (Fig. 5). The models with a stronger semiannual component in Ekman pumping over the dome region have a better skill in the semiannual cycle of the SD.
Because the seasonal variations in the curl and beta terms are closely related to the Indian monsoon and the meridional shift of the ITCZ, we have constructed time–latitude diagrams of the monthly climatology for precipitation, SST, and wind stress averaged over 60°–70°E (Fig. 8). The ITCZ remains to the north of the equator until December in some models (e.g., CSIRO-Mk3.0) even though that in the observation lies over the SD. One reason for this bias is the warmer Arabian Sea. As shown by Izumo et al. (2008), the strong Ekman upwelling in the western Arabian Sea off Oman associated with the summer monsoon leads to the cooling of SST there. It seems that this upwelling along the western boundary of the Arabian Sea is weaker in some models. In contrast, the ITCZ shifts too far south during boreal winter in some models (e.g., INGV-SXG and PCM). This results in erroneously strong westerly winds over the SD and leads to the positive skewness of the beta term.
4. Summary and discussion
Using outputs from 22 CMIP3 coupled GCMs, we have examined their skill in simulating the SD and its seasonal variation. Six models are successful in simulating the center of the dome within the accuracy of 5° in its location. However, almost all models underestimate the dominance of the semiannual harmonic when compared with the observation. This is because the simulated semiannual harmonic in the local Ekman pumping over the SD is weak. The main cause for this bias is the model’s inability to simulate the asymmetric nature of the monsoonal winds; the summer (winter) monsoon is shorter (longer) in its duration but stronger (weaker) in its magnitude.
It is interesting to note that Annamalai et al. (2007a) have shown that only GFDL-CM2.0, GFDL-CM2.1, ECHAM/Max Planck Institute Ocean Model (MPI-OM), MRI-CGCM2.3.2, PCM, and HadCM3 were successful in simulating the June–September rainfall climatology and seasonal evolution of the all-India rainfall. However, in our analysis, some of these models do not show high skill in simulating the semiannual cycle in the Ekman pumping over the SD, which is closely linked with the monsoon. The major reason for the discrepancy is the difference in the metrics used to evaluate the monsoon; those in Annamalai et al. (2007a) focused on the rainfall pattern correlations and root-mean-square differences during the summer monsoon, whereas ours take into account several factors including the monsoon transitions and the winter monsoon. It seems that models with a good skill in simulating the summer monsoon rainfall pattern do not necessarily succeed in simulating the winter monsoon.
The SST in the SD region may influence the interannual variation in the Indian summer monsoon, as discussed by Annamalai et al. (2005) and more recently by Izumo et al. (2008), and the tropical cyclone activity, as discussed by Xie et al. (2002). Because the realistic simulation of the SD is crucial for reproducing the climatic effects of the SD successfully, resolving the SD and its seasonal cycle may suitably provide a benchmark test for the performance of CGCMs.
Acknowledgments
We thank Dr. H. Annamalai and an anonymous reviewer for their valuable comments and suggestions. We acknowledge the modeling groups, the Program for Climate Model Diagnosis and Intercomparison (PCMDI), and the WCRP’s Working Group on Coupled Modelling (WGCM) for their roles in making available the WCRP CMIP3 multimodel dataset. Support of this dataset is provided by the Office of Science of the U.S. Department of Energy. The present research is supported by the Japan Society for Promotion of Science through Grant-in-Aid for Scientific Research (B) 20340125 for the senior author.
REFERENCES
An, S-I., and F-F. Jin, 2004: Nonlinearity and asymmetry of ENSO. J. Climate, 17 , 2399–2412.
Annamalai, H., P. Liu, and S-P. Xie, 2005: Southwest Indian Ocean SST variability: Its local effect and remote influence on Asian monsoons. J. Climate, 18 , 4150–4167.
Annamalai, H., K. Hamilton, and K. R. Sperber, 2007a: The South Asian summer monsoon and its relationship with ENSO in the IPCC AR4 simulations. J. Climate, 20 , 1071–1092.
Annamalai, H., H. Okajima, and M. Watanabe, 2007b: Possible impact of the Indian Ocean SST on the Northern Hemisphere circulation during El Niño. J. Climate, 20 , 3164–3189.
Antonov, J. I., S. Levitus, T. P. Boyer, M. E. Conkright, T. O’Brien, C. Stephens, and B. Trotsenko, 1998: Temperature of the Indian Ocean. Vol. 3, World Ocean Atlas 1998, NOAA Atlas NESDIS 29, 166 pp.
Behera, S. K., J-J. Luo, S. Masson, P. Delecluse, S. Gualdi, A. Navarra, and T. Yamagata, 2005: Paramount impact of the Indian Ocean dipole on the East African short rain: A CGCM study. J. Climate, 18 , 4514–4530.
Black, E., J. Slingo, and K. R. Sperber, 2003: An observational study of the relationship between excessively strong short rains in coastal East Africa and Indian Ocean SST. Mon. Wea. Rev., 131 , 74–94.
Chang, P., and Coauthors, 2006: Climate fluctuations of tropical coupled systems—The role of ocean dynamics. J. Climate, 19 , 5122–5174.
Guan, Z., and T. Yamagata, 2003: The unusual summer of 1994 in East Asia: IOD teleconnections. Geophys. Res. Lett., 30 , 1544. doi:10.1029/2002GL016831.
Hermes, J., and C. J. C. Reason, 2008: Annual cycle of the South Indian Ocean (Seychelles-Chagos) thermocline ridge in a regional ocean model. J. Geophys. Res., 113 , C04035. doi:10.1029/2007JC004363.
Huang, B., and J. L. Kinter, 2002: Interannual variability in the tropical Indian Ocean. J. Geophys. Res., 107 , 3199. doi:10.1029/2001JC001278.
Izumo, T., C. de Boyer Montégut, J-J. Luo, S. K. Behera, S. Masson, and T. Yamagata, 2008: The role of the western Arabian Sea upwelling in Indian monsoon rainfall variability. J. Climate, 21 , 5603–5623.
Kalnay, E., and Coauthors, 1996: The NCEP/NCAR 40-Year Reanalysis Project. Bull. Amer. Meteor. Soc., 77 , 437–471.
McCreary, J. P., P. K. Kundu, and R. L. Molinari, 1993: A numerical investigation of dynamics, thermodynamics, and mixed-layer processes in the Indian Ocean. Prog. Oceanogr., 31 , 181–244.
Murtugudde, R., J. P. McCreary, and A. J. Busalacchi, 2000: Oceanic processes associated with anomalous events in the Indian Ocean with relevance to 1997–1998. J. Geophys. Res., 105 , 3295–3306.
Rao, S. A., and S. K. Behera, 2005: Subsurface influence on SST in the tropical Indian Ocean: Structure and interannual variability. Dyn. Atmos. Oceans, 39 , 103–135.
Rao, S. A., S. K. Behera, Y. Masumoto, and T. Yamagata, 2002: Interannual subsurface variability in the tropical Indian Ocean with a special emphasis on the Indian Ocean dipole. Deep-Sea Res. II, 49 , 1549–1572.
Reverdin, G., and M. Fieux, 1987: Sections in the western Indian Ocean: Variability in the temperature structure. Deep-Sea Res., 34 , 601–626.
Saji, N. H., and T. Yamagata, 2003: Possible impacts of Indian Ocean dipole mode events on global climate. Climate Res., 25 , 151–169.
Saji, N. H., B. N. Goswami, P. N. Vinayachandran, and T. Yamagata, 1999: A dipole mode in the tropical Indian Ocean. Nature, 401 , 360–363.
Saji, N. H., S-P. Xie, and T. Yamagata, 2006: Tropical Indian Ocean variability in the IPCC twentieth-century climate simulations. J. Climate, 19 , 4397–4417.
Webster, P. J., A. Moore, J. Loschnigg, and M. Leben, 1999: Coupled ocean–atmosphere dynamics in the Indian Ocean during 1997–98. Nature, 401 , 356–360.
Woodberry, K. E., M. Luther, and J. J. O’Brien, 1989: The wind-driven seasonal circulation in the southern tropical Indian Ocean. J. Geophys. Res., 94 , 17985–18002.
Xie, P., and P. A. Arkin, 1996: Analyses of global monthly precipitation using gauge observations, satellite estimates, and numerical model predictions. J. Climate, 9 , 840–858.
Xie, S-P., H. Annamalai, F. A. Schott, and J. P. McCreary, 2002: Structure and mechanisms of south Indian Ocean climate variability. J. Climate, 15 , 864–878.
Yamagata, T., S. K. Behera, J-J. Luo, S. Masson, M. R. Jury, and S. A. Rao, 2004: Coupled ocean–atmosphere variability in the tropical Indian Ocean. Earth’s Climate: The Ocean–Atmosphere Interaction, Geophys. Monogr., Vol. 147, Amer. Geophys. Union, 189–211.
Yokoi, T., T. Tozuka, and T. Yamagata, 2008: Seasonal variation of the Seychelles Dome. J. Climate, 21 , 3740–3754.
Brief description of the CMIP3 models (ID is identifier).