Future climate change can be due to internal and external variability  . In a recent study,  noted that the observed trends of climate change are very unlikely attributed to internal variability even if the contemporary climate models are found to underestimate it. The future projection of climate change obtained by state-of-art models may be uncertain mainly because of the difficulty in the determination of quantitative and spatial distribution of future greenhouse gases and other radiatively important external parameters. Also the knowledge how the climate system responds to these external influences is still incomplete. To complicate further it is uncertain how useful can a model be for climate change studies if its numerous parameters are tuned to present climate  . Since seasonal changes display climate changes, the ability of a model to reproduce the seasonal cycles is a positive indication that the model can also be used for the investigation of climate change  . Thus models can predict future climate accurately only if they can reproduce seasonal cycles accurately. Also the seasonal changes for example in surface air temperature, are much larger than the changes in the long periods, even centuries and it is unwise to ignore large ones compared to small climate shifts.
In the extra-tropics surface air temperature determines well the climate and climate change, while in the tropics rainfall variations are of greater significance.  determined the changes the phase and amplitude of the annual cycle of surface air temperature in extra tropics and found an advance of 1.7 days in phase for the land area between 1954 and 2007. They also found significant changes in the amplitude. None of the IPCC 5 models reproduce the shift towards earlier seasons. However,  examined only the annual cycle. But the seasonal shift is given not only by the first harmonic of 12 months period but the other higher harmonics also contribute to the seasonal shift.
Billions of people around the world in the tropics depend crucially on the rainfall for their survival. The seasonal cycles in precipitation of the tropics are known as monsoons, the most energetic part of tropical climate system. Here we find the veracity of IPCC models to reproduce the observed seasonal cycles in precipitation in the tropics, thus finding how accurately the models can predict climate change. For this purpose we calculated the amplitude and phase of the first four harmonics (not only the first) in the rainfall over land area from 10 IPCC models: Miroc3.2  , CCSM3  , HadCM3  , ECHAM5  , Mk3.0  , CGCM2.3.2  , CM4  , BCM2.0  , CGCM3  and CM3.0  . Such an analysis also gives an idea how accurately the present state-of-art models reproduce the global monsoons. We determine the accuracy with which a suite of ten coupled atmospheric general circulation models reproduce the observed seasonal cycles in rainfall. This gives an idea how well these models can predict climate change.
2. Data Sources
For the land area of the tropics (50˚N - 50˚S) the rainfall data were obtained for the period 1960-1990 from the University of Delaware (1950-1999)  . We calculated the amplitude and phase of the first four harmonics in the rainfall over land area from 10 IPCC models (http://www.ipcc-data.org/ar4/gcm_data.html) (see Table 1). The first four harmonics were obtained from 30 years mean monthly rainfall using the Fourier series. The calculation of amplitude and phase using the Fourier series ensures their determination by least square method  .
Table 1. Ten IPCC models used in the study.
3. Results and Discussion
The left hand side of Figure 1 shows the observed annual, January and July rainfall. January and July represent the typical boreal winter and summer and austral summer and winter months respectively. Thus the change in rainfall between these two months gives the seasonal change. The middle figure shows the corresponding rainfall in multi model ensemble (MME) of the 10 Models. The right hand side shows the standard deviation amongst the 10 models.
In observations, (top of Figure 1) the principle regions of rainfalls the Amazon basin, Congo basin and Indonesia can be seen in the annual mean. In the MME these three regions are reproduced well. But, there are quantitative differences both in magnitude and spatial extension, particularly over western Amazon region. Top right side figure shows variations amongst the ten models, again over western Amazon and Indonesia, varying from 1.5 - 2.5 per day. The models seem to have difficulty in obtaining the correct rainfall near mountainous region such as Andes and Himalaya. The generally high (low) rainfall in summer (winter) months of northern and southern hemispheres can be seen in the middle and bottom figures, thus simulating the seasonal variations. The monsoonal
Figure 1. Left hand side upper figure: Observed annual rainfall in mm/day; Middle figure: Observed January rainfall in mm/day; Lower figure: Observed July rainfall in mm/day.
rainfall variations are reproduced well in January and July by MME, for example, the east Asian Indian and South American monsoon regions show high (low) rainfall in summer (winter).
Figure 2 shows the variance explained by the first two (annual and semiannual) harmonics both in observations and MME. The lower figure shows the sum of the variance of the first two harmonics. It is clear from the figure that at least the first two harmonics are necessary to explain high percentage of variance of the annual variation of rainfall. Some studies of global monsoons   included only the annual (first) harmonic.
Figure 3(a) and Figure 3(b) show respectively the amplitude (h1 and h2) and phase (φ1 and φ2) of the first and second harmonics in the observations and MME and the difference between the two. In the lower part of the Figure 3(a) (left figure) MME seem to reproduce the amplitude values reasonably well over subtropical regions of South America, Africa and Australia and also over the tropical South America and India, while there are phase differences (right figure) of 10 days or more over some regions. In Figure 3(b), it can be seen that the amplitude of the second harmonic is substantial.  noted the important of the semi-annual component over Indian monsoon region. The MME showed the high amplitude in semi-annual component over western Amazonia and Congo basin. The lower part of the Figure 3(b) shows in some specific regions there are some differences between the observations and MME both in amplitude and phase of semi-annual component.
The root mean square (RMS) error, absolute error (ABSE) and the bias the models make with respect to observations in determining the amplitude and phase of the first four harmonics and mean annual rainfall are given Table 2, for the 10 models and MME. Also given are the spatial correlation coefficients (CC) between the observed and model amplitudes and phases. In general errors are
Figure 2. Observed and MME variance explained by first and second (annual and semiannual) harmonics.
Figure 3. (a) Upper figure: Observed and MME amplitude of first (annual) harmonic; Lower figure: left hand side: the difference between the amplitude of first (annual) harmonic in MME and the observations; right hand side: difference in phase between the observations for the first (annual) harmonic; (b) Same as the (a) except for the second harmonic.
Table 2. Bias, ABSE and RMSE and correlation in the spatial distribution of rainfall between the observed and the ten models. Last rows: grid size of models.
relatively less and CC are higher for the MME. High CC of 0.88 is found for the MME with mean annual rainfall. In the case of the annual harmonic the CC is 0.86. Even in the case of other harmonics CCs are high for MME. This shows that the MME is able to reproduce global monsoons well. Even in the case of some individual models the CCs are high. All the models seem to do well in reproducing mean annual rainfall. Although the models have some difficulty in reproducing the amplitude of the four harmonics, they seem to do well in reproducing phases. Note that the CCs are high for the phases. This means that all the models and MME do well in reproducing the phase of global monsoons such as onset and withdrawal.
Now we examine how these models and MME perform for the regional monsoons over India and South America.
3.1. India Monsoon
Figure 4. Same as Figure 1, but for the Indian monsoon.
figures shows that the MME is able represent the general patterns of variations of annual rainfall, although there are regional differences. The Western Peninsular maximum is reproduced well. But the MME over estimates rainfall in southeastern peninsula. Top right hand side figure shows that there are variations in the models with maximum variations around 30˚N around 90˚E.
Middle figure shows that the MME has some difficulty in reproducing the January rainfall quantitatively although there is qualitative agreement. Except over the Gujarat state and around the standard deviation amongst the models is small.
The lower figure shows that in the principal summer monsoon month, July there is over all agreement between the observations and the MME in the rainfall although there are some regional differences. The lower right hand side figure shows that in central India there are large differences amongst the models.
Figure 5 shows the variance explained by the first and second harmonics for
Figure 5. Same as Figure 2, but for the Indian monsoon.
the annual rainfall. It is very clear from the figure that in the observations the second harmonic explains large variance. Note that the lower figure shows the sum of the variance explained by first and second harmonics. Right hand side figure shows that the model is able to capture the importance of the second harmonic over the Indian monsoon region. Earlier studies,   did not include the second harmonic.
Figure 6(a) and Figure 6(b) show that the MME is able to capture very well the amplitude of first and second harmonics. There are some differences only in the northern part. In the case of phase there are differences only in the east and west side with good agreement in the middle.
The RMS and ABSE and bias the models make with respect to observations in determining the amplitude and phase of first four harmonics and mean annual rainfall for the Indian monsoon are given in Table 3. Also given are the spatial CCs between the observed and model amplitudes and phases. The CCs are high for MME and in many cases they are more than 0.8. High CC of 0.92 is found for the amplitude of first harmonic. All the models and MME seem to have same difficulty in reproducing higher harmonics. Just like in the case of global monsoons, the CCs for the phases are high in individual models and MME. This shows models and MME do well in representing the phase of the Indian monsoon such the onset and withdrawal.
3.2. South American Monsoon
Figure 7 is again similar to Figure 1. A comparison of top left hand side and middle figures shows that the MME is able to reproduce well the general distribution of rainfall. It is interesting to note that the northeast semi-arid region is reproduced well although there are quantitative differences. The northwestern maximum rainfall over northeast Amazonia is not reproduced. Top right hand side figure show there are large variations amongst the models are over western South America and also over the Sao Paulo, Rio de Janeiro states.
The middle figure shows that the rainfall in the principal summer monsoon month January is reproduced well but there are quantitative differences. The MME has difficulty in reproducing the high rainfall over northwest region of South America. Further it falsely produces a dry region over northern region. The dry region over northeast Brazil is also not reproduced well. The right hand side figure in the middle shows that the all the models falsely reproduce the dry region. There are large differences amongst the models over western South America.
The lower figure shows very good agreement between MME and the observations. The right hand side figure shows all the models perform well in reproducing July rainfall in austral winter. The maximum rainfall noted in the observations over south Brazil is not reproduced in MME.
Figure 8 shows the variance explained by the first harmonic and the total variance explained by the first two harmonics. In the observations one can note the
Figure 6. Same as Figure 3(a) and Figure 3(b), but for the Indian monsoon.
Figure 7. Same as Figure 1, but for the South American monsoon.
Figure 8. Same as Figure 2, but for the South American monsoon.
importance for the second harmonic. MME is able to reproduce the variance explained by the first two harmonic except over the northwest and southeast South America.
Figure 9(a) and Figure 9(b) show that the MME able to capture the importance of the first harmonic in central Brazil. There are large differences in the MME and observations in middle of South America. But the differences are very less over northern and southern parts in the annual component. In the semi-annual component also the differences are large over eastern and western South America. Regarding the phase except over the region of Mato Grosso and Minas Gerais states the differences are less for the annual harmonic. The differences are very low over northeast Brazil. The differences are large over Argentina for the annual harmonic but the differences are low in the semi-annual
Figure 9. Same as Figure 3(a) and Figure 3(b), but for the South American monsoon.
Table 4 gives similar to Table 2 and Table 3, the ABSE, RMSE, bias and the correlation. The CCs are small for the amplitudes unlike those of Table 1 and Table 2. But the CCs are very large for the phases just as in the case of global and Indian monsoons. This shows that the MME and individual models are able to capture the phase of South American monsoon such as the onset and withdrawal.
In the above we have seen that the individual models and MME reproduced reasonably well the global monsoons and Indian and South American monsoons. That is they reproduce well the seasonal cycles in rainfall there by suggesting that they are useful in prediction of future climate change and the mitigation adverse climate change.
To our knowledge, a systematic attempt has not yet been made to assess quantitatively the capability of models in reproduction seasonal cycles in tropical
precipitation. Since seasonal cycles exhibit climate change, models can predict climate change accurately only when they reproduce seasonal cycles accurately. Seasonal cycles in tropical precipitation are known as monsoons.
In this paper, we use a suite of 10 state-of-art models and observations to evaluate how accurately the models reproduce seasonal cycles in tropical precipitation or monsoons. Earlier studies of monsoons use only the first harmonic of
yearly (12 months) rainfall. Seasonal cycles can be reproduced only with at least first two harmonics. We calculated the root mean square error, absolute error and bias the models make with respect to observations in reproducing first four harmonics of 12 months cycle. We find that the individual models and also the multi model ensemble reproduce reasonably well the seasonal cycles in rainfall thus suggesting their utility in climate prediction and simulation of global and regional monsoons.
 Collins, W.D., Bitz, C.M., Blackmon, M.L., Bonan, G.B., Bretherton, C.S., Carton, J.A., Chang, P., Doney, S.C., Hack, J.J., Henderson, T.B., Kiehl, J.T., Large, W.G., McKenna, D.S., Santer, B.D. and Smith, R.D. (2006) The Community Climate System Model Version 3 (CCSM3). Journal of Climate, 19, 2122-2143.
 Jones, C.D., Gregory, J.M., Thorpe, R.B., Cox, P.M., Murphy, J.M., Sexton, D.M.H. and Valdes, P. (2004) Systematic Optimisation and Climate Simulation of FAMOUS, a Fast Version of HadCM3. Hadley Centre Technical Note 60, Hadley Centre, Exeter.
 Roeckner, E., Bäuml, G., Bonaventura, L., Brokopf, R., Esch, M., Giorgetta, M., Hagemann, S., Kirchner, I., Kornblueh, L., Manzini, E., Rhodin, A., Schlese, U., Schulzweida, U. and Tompkins, A. (2003) The Atmospheric General Circulation Model ECHAM 5. PART I: Model Description. Max Planck Institute for Meteorology Rep. 349, 127 p.
 Gordon, H.B., Rotstayn, I.D., McGregor, J.I., Dix, M.R., Kowalezyk, O’Farrell, S.P., Waterman, I.J., Hirst, A.C., Wilson, S.G., Collier, M.A., Watterson, I.G. and Elliot, T.I. (2002) The CSIRO Mk3 Climate System Model. Asoendale: CSIRO Atmospheric Research Technical Paper No. 60.
 Yukimoto, S., et al. (2001) The New Meteorological Research Institute Coupled GCM (MRI-CGM2)-Model Climate and Variability. Meteorology and Geophysics, 52, 47-88.
 Marti, O., Braconnot, P., Bellier, J., Benshila, R., Bony, S., Brockmann, P., Cadule, P., Caubel, A., Denvil, S., Dufresne, J.L., Fairhead, L., Filiberti, M.-A., Foujols, M.-A., Fichefet, T., Friedlingstein, P., Goosse, H., Grandpeix, J.Y., Hourdin, F., Krinner, G., Lévy, C., Madec, G., Musat, I., de Noblet, N., Polcher, J. and Talandier, C. (2005) The New IPSL Climate System Model: IPSL-CM4. Note du Pôle de Modélisation, 26, 1-86.
 Furevik, T., Bentsen, M., Drange, H., Kindem, I.K.T., Kvansto, N.G. and Sorteberg, A. (2003) Description and Evaluation of the Bergen Climate Model: ARPEGE Coupled with MICOM. Climate Dynamics, 21, 27-51.
 Lee, J.L., Wang, B., Kang, I.S., Shukla, J., Kumar, A., Kuj, J.S., Scheman, J.K.E., Luo, J.J., Yamagata, T., Fu, X., Alles, O., Stezn, B., Rosati, T. and Park, C.K. (2010) How Are Seasonal Prediction Skills Related to Models’ Performance on Mean State and Annual Cycle? Climate Dynamics, 35, 267-283.