Dynamics of ''internal'' interannual variability of the Indian summer monsoon in a GCM

1] The poor predictability of the Indian summer monsoon ( ISM) appears to be due to the fact that a large fraction of interannual variability (IAV) is governed by unpredictable "internal'' low frequency variations. Mechanisms responsible for the internal IAV of the monsoon have not been clearly identified. Here, an attempt has been made to gain insight regarding the origin of internal IAV of the seasonal ( June - September, JJAS) mean rainfall from "internal'' IAV of the ISM simulated by an atmospheric general circulation model (AGCM) driven by fixed annual cycle of sea surface temperature (SST). The underlying hypothesis that monsoon ISOs are responsible for internal IAV of the ISM is tested. The spatial and temporal characteristics of simulated summer intraseasonal oscillations ( ISOs) are found to be in good agreement with those observed. A long integration with the AGCM forced with observed SST, shows that ISO activity over the Asian monsoon region is not modulated by the observed SST variations. The internal IAV of ISM, therefore, appears to be decoupled from external IAV. Hence, insight gained from this study may be useful in understanding the observed internal IAV of ISM. The spatial structure of the ISOs has a significant projection on the spatial structure of the seasonal mean and a common spatial mode governs both intraseasonal and interannual variability. Statistical average of ISO anomalies over the season ( seasonal ISO bias) strengthens or weakens the seasonal mean. It is shown that interannual anomalies of seasonal mean are closely related to the seasonal mean of intraseasonal anomalies and explain about 50% of the IAV of the seasonal mean. The seasonal mean ISO bias arises partly due to the broad-band nature of the ISO spectrum allowing the time series to be aperiodic over the season and partly due to a non-linear process where the amplitude of ISO activity is proportional to the seasonal bias of ISO anomalies. The later relation is a manifestation of the binomial character of rainfall time series. The remaining 50% of the IAV may arise due to land-surface processes, interaction between high frequency variability and ISOs, etc.

Intraseasonal response of mixed layer temperature and salinity in the Bay of Bengal to heat and freshwater flux

Buoy and satellite data show pronounced subseasonal oscillations of sea surface temperature (SST) in the summertime Bay of Bengal. The SST oscillations are forced mainly by surface heat flux associated with the active break cycle of the south Asian summer monsoon. The input of freshwater (FW) from summer rain and rivers to the bay is large, but not much is known about subseasonal salinity variability. We use 2002-2007 observations from three Argo floats with 5 day repeat cycle to study the subseasonal response of temperature and salinity to surface heat and freshwater flux in the central Bay of Bengal. About 95% of Argo profiles show a shallow halocline, with substantial variability of mixed layer salinity. Estimates of surface heat and freshwater flux are based on daily satellite data sampled along the float trajectory. We find that intraseasonal variability of mixed layer temperature is mainly a response to net surface heat flux minus penetrative radiation during the summer monsoon season. In winter and spring, however, temperature variability appears to be mainly due to lateral advection rather than local heat flux. Variability of mixed layer freshwater content is generally independent of local surface flux (precipitation minus evaporation) in all seasons. There are occasions when intense monsoon rainfall leads to local freshening, but these are rare. Large fluctuations in FW appear to be due to advection, suggesting that freshwater from rivers and rain moves in eddies or filaments.

Climatic potential of delta O-18 of Abies spectabilis from the Nepal Himalaya

A 50-year tree-ring delta O-18 chronology of Abies spectabilis growing close to the tree line (3850 m asl) in the Nepal Himalaya is established to explore its dendroclimatic potential. Response function analysis with ambient climatic records revealed that tree-ring delta O-18 is primarily governed by rainfall during the monsoon season (June September), and the regression model accounts for 35% of the variance in rainfall. Extreme dry years identified in instrumental weather data are detected in the delta O-18 chronology. Further, tree-ring delta O-18 is much more sensitive to rainfall fluctuations than other tree-ring parameters such as width and density typically used in dendroclimatology. Correlation analyses with Nino 3.4 SST reveal time-dependent behavior of ENSO-monsoon relationships. (C) 2009 Elsevier GmbH. All rights reserved.

Fast versus slow response in climate change: implications for the global hydrological cycle

Recent studies have shown that changes in global mean precipitation are larger for solar forcing than for CO2 forcing of similar magnitude.In this paper, we use an atmospheric general circulation model to show that the differences originate from differing fast responses of the climate system. We estimate the adjusted radiative forcing and fast response using Hansen's ``fixed-SST forcing'' method.Total climate system response is calculated using mixed layer simulations using the same model. Our analysis shows that the fast response is almost 40% of the total response for few key variables like precipitation and evaporation. We further demonstrate that the hydrologic sensitivity, defined as the change in global mean precipitation per unit warming, is the same for the two forcings when the fast responses are excluded from the definition of hydrologic sensitivity, suggesting that the slow response (feedback) of the hydrological cycle is independent of the forcing mechanism. Based on our results, we recommend that the fast and slow response be compared separately in multi-model intercomparisons to discover and understand robust responses in hydrologic cycle. The significance of this study to geoengineering is discussed.

Interannual Variations of Sea Surface Temperature over the Arabian Sea and the Indian Monsoon: A New Perspective

The interannual variation of surface fields over the Arabian Sea and Bay of Bengal are studied using data between 1900 and 1979. It is emphasized that the monthly mean sea surface temperature (SST) over the north Indian Ocean and monsoon rainfall are significantly affected by synoptic systems and other intraseasonal variations. To highlight the interannual signals it is important to remove the large-amplitude high-frequency noise and very low frequency long-term trends, if any. By suitable spatial and temporal averaging of the SST and the rainfall data and by removing the long-term trend from the SST data, we have been able to show that there exists a homogeneous region in the southeastern Arabian Sea over which the March�April (MA) SST anomalies are significantly correlated with the seasonal (June�September) rainfall over India. A potential of this premonsoon signal for predicting the seasonal rainfall over India is indicated. It is shown that the correlation between the SST and the seasonal monsoon rainfall goes through a change of sign from significantly positive with premonsoon SST to very small values with SST during the monsoon season and to significantly negative with SST during the post-monsoon months. For the first time, we have demonstrated that heavy or deficient rainfall years are associated with large-scale coherent changes in the SST (although perhaps of small amplitude) over the north Indian 0cean. We also indicate possible reasons for the apparent lack of persistence of the premonsoon SST anomalies.

Towards understanding the unusual Indian monsoon in 2009

The Indian summer monsoon season of 2009 commenced with a massive deficit in all-India rainfall of 48% of the average rainfall in June. The all-India rainfall in July was close to the normal but that in August was deficit by 27%. In this paper, we first focus on June 2009, elucidating the special features and attempting to identify the factors that could have led to the large deficit in rainfall. In June 2009, the phase of the two important modes, viz., El Nino and Southern Oscillation (ENSO) and the equatorial Indian Ocean Oscillation (EQUINOO) was unfavourable. Also, the eastern equatorial Indian Ocean (EEIO) was warmer than in other years and much warmer than the Bay. In almost all the years, the opposite is true, i.e., the Bay is warmer than EEIO in June. It appears that this SST gradient gave an edge to the tropical convergence zone over the eastern equatorial Indian Ocean, in competition with the organized convection over the Bay. Thus, convection was not sustained for more than three or four days over the Bay and no northward propagations occurred. We suggest that the reversal of the sea surface temperature (SST) gradient between the Bay of Bengal and EEIO, played a critical role in the rainfall deficit over the Bay and hence the Indian region. We also suggest that suppression of convection over EEIO in association with the El Nino led to a positive phase of EQUINOO in July and hence revival of the monsoon despite the El Nino. It appears that the transition to a negative phase of EQUINOO in August and the associated large deficit in monsoon rainfall can also be attributed to the El Nino.

Indian Ocean Dipole: Processes and impacts

Equatorial Indian Ocean is warmer in the east, has a deeper thermocline and mixed layer, and supports a more convective atmosphere than in the west. During certain years, the eastern Indian Ocean becomes unusually cold, anomalous winds blow from east to west along the equator and southeastward off the coast of Sumatra, thermocline and mixed layer lift up and the atmospheric convection gets suppressed. At the same time, western Indian Ocean becomes warmer and enhances atmospheric convection. This coupled ocean-atmospheric phenomenon in which convection, winds, sea surface temperature (SST) and thermocline take part actively is known as the Indian Ocean Dipole (IOD). Propagation of baroclinic Kelvin and Rossby waves excited by anomalous winds, play an important role in the development of SST anomalies associated with the IOD. Since mean thermocline in the Indian Ocean is deep compared to the Pacific, it was believed for a long time that the Indian Ocean is passive and merely responds to the atmospheric forcing. Discovery of the IOD and studies that followed demonstrate that the Indian Ocean can sustain its own intrinsic coupled ocean-atmosphere processes. About 50% percent of the IOD events in the past 100 years have co-occurred with El Nino Southern Oscillation (ENSO) and the other half independently. Coupled models have been able to reproduce IOD events and process experiments by such models – switching ENSO on and off – support the hypothesis based on observations that IOD events develop either in the presence or absence of ENSO. There is a general consensus among different coupled models as well as analysis of data that IOD events co-occurring during the ENSO are forced by a zonal shift in the descending branch of Walker cell over to the eastern Indian Ocean. Processes that initiate the IOD in the absence of ENSO are not clear, although several studies suggest that anomalies of Hadley circulation are the most probable forcing function. Impact of the IOD is felt in the vicinity of Indian Ocean as well as in remote regions. During IOD events, biological productivity of the eastern Indian Ocean increases and this in turn leads to death of corals over a large area.Moreover, the IOD affects rainfall over the maritime continent, Indian subcontinent, Australia and eastern Africa. The maritime continent and Australia suffer from deficit rainfall whereas India and east Africa receive excess. Despite the successful hindcast of the 2006 IOD by a coupled model, forecasting IOD events and their implications to rainfall variability remains a major challenge as understanding reasons behind an increase in frequency of IOD events in recent decades.

Aperiodic Variability in the Cane-Zebiak Model: A Diagnostic Study

The finite predictability of the coupled ocean-atmosphere system is determined by its aperiodic variability. To gain insight regarding the predictability of such a system, a series of diagnostic studies has been carried out to investigate the role of convergence feedback in producing the aperiodic behavior of the standard version of the Cane-Zebiak model. In this model, an increase in sea surface temperature (SST) increases atmospheric heating by enhancing local evaporation (SST anomaly feedback) and low-level convergence (convergence feedback). The convergence feedback is a nonlinear function of the background mean convergence field. For the set of standard parameters used in the model, it is shown that the convergence feedback contributes importantly to the aperiodic behaviour of the model. As the strength of the convergence feedback is increased from zero to its standard value, the model variability goes from a periodic regime to an aperiodic regime through a broadening of the frequency spectrum around the basic periodicity of about 4 years. Examination of the forcing associated with the convergence feedback reveals that it is intermittent, with relatively large amplitude only during 2 or 3 months in the early part of the calendar year. This seasonality in the efficiency of the convergence feedback is related to the strong seasonality of the mean convergence over the eastern Pacific. It is shown that if the mean convergence field is fixed at its March value, aperiodic behavior is produced even in the absence of annual cycles in the other mean fields. On the, other hand, if the mean convergence field is fixed at its September value, the coupled model evolution remains close to periodic, even in the presence of the annual cycle in the other fields. The role of convergence feedback on the aperiodic variability of the model for other parameter regimes is also examined. It is shown that a range exists in the strength of the SST anomaly feedback for which the model variability is aperiodic even without the convergence feedback. It appears that in the absence of convergence feedback, enhancement of the strength of the air-sea coupling in the model through other physical processes also results in aperiodicity in the model.

Seasonal prediction of the Indian monsoon

Under the project `Seasonal Prediction of the Indian Monsoon' (SPIM), the prediction of Indian summer monsoon rainfall by five atmospheric general circulation models (AGCMs) during 1985-2004 was assessed. The project was a collaborative effort of the coordinators and scientists from the different modelling groups across the country. All the runs were made at the Centre for Development of Advanced Computing (CDAC) at Bangalore on the PARAM Padma supercomputing system. Two sets of simulations were made for this purpose. In the first set, the AGCMs were forced by the observed sea surface temperature (SST) for May-September during 1985-2004. In the second set, runs were made for 1987, 1988, 1994, 1997 and 2002 forced by SST which was obtained by assuming that the April anomalies persist during May-September. The results of the first set of runs show, as expected from earlier studies, that none of the models were able to simulate the correct sign of the anomaly of the Indian summer monsoon rainfall for all the years. However, among the five models, one simulated the correct sign in the largest number of years and the second model showed maximum skill in the simulation of the extremes (i.e. droughts or excess rainfall years). The first set of runs showed some common bias which could arise either from an excessive sensitivity of the models to El Nino Southern Oscillation (ENSO) or an inability of the models to simulate the link of the Indian monsoon rainfall to Equatorial Indian Ocean Oscillation (EQUINOO), or both. Analysis of the second set of runs showed that with a weaker ENSO forcing, some models could simulate the link with EQUINOO, suggesting that the errors in the monsoon simulations with observed SST by these models could be attributed to unrealistically high sensitivity to ENSO.

Tropical deep convection, convective available potential energy and sea surface temperature

The convective available potential energy (CAFE) based on monthly mean sounding has been shown to be relevant to deep convection in the tropics. The variation of CAFE with SST has been found to be similar to the variation of the frequency of deep convection at one station each in the tropical Atlantic and W. Pacific oceans. This suggests a strong link between the frequency of tropical convection and CAFE. It has been shown that CAFE so derived can be interpreted as the work potential of the atmosphere above the boundary layer with ascent in the convective region and subsidence in the surrounding cloud-free region.

Interannual variations of Indian summer monsoon in a GCM: External conditions versus internal feedbacks

The potential predictability of the Indian summer monsoon due to slowly varying sea surface temperature (SST) forcing is examined. Factors responsible for limiting the predictability are also investigated. Three multiyear simulations with the R30 version of the Geophysical Fluid Dynamics Laboratory's climate model are carried out for this purpose, The mean monsoon simulated by this model is realistic including the mean summer precipitation over the Indian continent. The interannual variability of the large-scale component of the monsoon such as the "monsoon shear index" and its teleconnection with Pacific SST is well simulated by the model in a 15-yr integration with observed SST as boundary condition. On regional scales, the skill in simulating the interannual variability of precipitation over the Indian continent by the model is rather modest and its simultaneous correlation with eastern Pacific SST is negative but poor as observed. The poor predictability of precipitation over the Indian region in the model is related to the fact that contribution to the interannual variability over this region due to slow SST variations [El Nino-Southern Oscillation (ENSO) related] is comparable to those due to regional-scale fluctuations unrelated to ENSO SST. The physical mechanism through which ENSO SST tend to produce reduction in precipitation over the Indian continent is also elucidated. A measure of internal variability of the model summer monsoon is obtained from a 20-yr integration of the same model with fixed annual cycle SST as boundary conditions but with predicted soil moisture and snow cover. A comparison of summer monsoon indexes between this run and the observed SST run shows that the internal oscillations can account for a large fraction of the simulated monsoon variability. The regional-scale oscillations in the observed SST run seems to arise from these internal oscillations. It is discovered that most of the interannual internal variability is due to an internal quasi-biennial oscillation (QBO) of the model atmosphere. Such a QBO is also found in the author's third 18-yr simulation in which fixed annual cycle of SST as well as soil moisture and snow cover are prescribed. This shows that the model QBO is not due to land-surface-atmosphere interaction. It is proposed that the model QBO arises due to an interaction between nonlinear intraseasonal oscillations and the annual cycle. Spatial structure of the QBO and its role in limiting the predictability of the Indian summer monsoon is discussed.

Indian Ocean dipole mode events in an ocean general circulation model

The evolution of the dipole mode (DM) events in the Indian Ocean is examined using an ocean model that is driven by the NCEP fluxes for the period 1975-1998. The positive DM events during 1997, 1994 and 1982 and negative DM events during 1996 and 1984-1985 are captured by the model and it reproduces both the surface and subsurface features associated with these events. In its positive phase, the DM is characterized by warmer than normal SST in the western Indian Ocean and cooler than normal SST in the eastern Indian Ocean. The DM events are accompanied by easterly wind anomalies along the equatorial Indian Ocean and upwelling-favorable alongshore wind anomalies along the coast of Sumatra. The Wyrtki jets are weak during positive DM events, and the thermocline is shallower than normal in the eastern Indian Ocean and deeper in the west. This anomaly pattern reverses during negative DM events. During the positive phase of the DM easterly wind anomalies excite an upwelling equatorial Kelvin wave. This Kelvin wave reflects from the eastern boundary as an upwelling Rossby wave which propagates westward across the equatorial Indian Ocean. The anomalies in the eastern Indian Ocean weaken after the Rossby wave passes. A similar process excites a downwelling Rossby wave during the negative phase. This Rossby wave is much weaker but wind forcing in the central equatorial Indian Ocean amplifies the downwelling and increases its westward phase speed. This Rossby wave initiates the deepening of the thermocline in the western Indian Ocean during the following positive phase of the DM. Rossby wave generated in the southern tropical Indian Ocean by Ekman pumping contributes to this warming. Concurrently, the temperature equation of the model shows upwelling and downwelling to be the most important mechanism during both positive events of 1994 and 1997. (C) 2002 Elsevier Science Ltd. All rights reserved.

Numerical Modelling of the Observed Diurnal Cycle of Indian Monsoon Rainfall

We have analysed the diurnal cycle of rainfall over the Indian region (10S-35N, 60E-100E) using both satellite and in-situ data, and found many interesting features associated with this fundamental, yet under-explored, mode of variability. Since there is a distinct and strong diurnal mode of variability associated with the Indian summer monsoon rainfall, we evaluate the ability of the Weather Research and Forecasting Model (WRF) to simulate the observed diurnal rainfall characteristics. The model (at 54km grid-spacing) is integrated for the month of July, 2006, since this period was particularly favourable for the study of diurnal cycle. We first evaluate the sensitivity of the model to the prescribed sea surface temperature (SST), by using two different SST datasets, namely, Final Analyses (FNL) and Real-time Global (RTG). It was found that with RTG SST the rainfall simulation over central India (CI) was significantly better than that with FNL. On the other hand, over the Bay of Bengal (BoB), rainfall simulated with FNL was marginally better than with RTG. However, the overall performance of RTG SST was found to be better than FNL, and hence it was used for further model simulations. Next, we investigated the role of the convective parameterization scheme on the simulation of diurnal cycle of rainfall. We found that the Kain-Fritsch (KF) scheme performs significantly better than Betts-Miller-Janjić (BMJ) and Grell-Devenyi schemes. We also studied the impact of other physical parameterizations, namely, microphysics, boundary layer, land surface, and the radiation parameterization, on the simulation of diurnal cycle of rainfall, and identified the “best” model configuration. We used this configuration of the “best” model to perform a sensitivity study on the role of various convective components used in the KF scheme. In particular, we studied the role of convective downdrafts, convective timescale, and feedback fraction, on the simulated diurnal cycle of rainfall. The “best” model simulations, in general, show a good agreement with observations. Specifically, (i) Over CI, the simulated diurnal rainfall peak is at 1430 IST, in comparison to the observed 1430-1730 IST peak; (ii) Over Western Ghats and Burmese mountains, the model simulates a diurnal rainfall peak at 1430 IST, as opposed to the observed peak of 1430-1730 IST; (iii) Over Sumatra, both model and observations show a diurnal peak at 1730 IST; (iv) The observed southward propagating diurnal rainfall bands over BoB are weakly simulated by WRF. Besides the diurnal cycle of rainfall, the mean spatial pattern of total rainfall and its partitioning between the convective and stratiform components, are also well simulated. The “best” model configuration was used to conduct two nested simulations with one-way, three-level nesting (54-18-6km) over CI and BoB. While, the 54km and 18km simulations were conducted for the whole of July, 2006, the 6km simulation was carried out for the period 18 - 24 July, 2006. The results of our coarse- and fine-scale numerical simulations of the diurnal cycle of monsoon rainfall will be discussed.

Some salient features of the atmosphere observed over the north Bay of Bengal during BOBMEX

This paper describes the near surface characteristics and vertical variations based on the observations made at 17.5degreesN and 89degreesE from ORV Sagar Kanya in the north Bay of Bengal during the Bay of Bengal Monsoon Experiment (BOBMEX) carried out in July-August 1999. BOBMEX captured both the active and weak phases of convection. SST remained above the convection threshold throughout the BOBMEX. While the response of the SST to atmospheric forcing was clearly observed, the response of the atmosphere to SST changes was not clear. SST decreased during periods of large scale precipitation, and increased during a weak phase of convection. It is shown that the latent heat flux at comparable wind speeds was about 25-50% lower over the Bay during BOBMEX compared to that over the Indian Ocean during other seasons and tropical west Pacific. On the other hand, the largest variations in the surface daily net heat flux are observed over the Bay during BOBMEX. SST predicted using observed surface fluxes showed that 1-D heat balance model works sometime but not always, and horizontal advection is important. The high resolution Vaisala radiosondes launched during BOBMEX could clearly bring out the changes in the vertical structure of the atmosphere between active and weak phases of convection. Convective Available Potential Energy of the surface air decreased,by 2-3 kJ kg(-1) following convection, and recovered in a time period of one or two days. The mid tropospheric relative humidity and water vapor content, and wind direction show the major changes between the active and weak phases of convection.

Mechanisms of summer intraseasonal sea surface temperature oscillations in the Bay of Bengal

Intraseasonal variations (ISV) of sea surface temperature (SST) in the Bay of Bengal (BoB) is highest in its northwestern part. An Indian Ocean model forced by QuikSCAT winds and climatological river discharge (QR run) reproduces ISV of SST, albeit with weaker magnitude. Air-sea fluxes, in the presence of a shallow mixed layer, efficiently effect intraseasonal SST fluctuations. Warming during intraseasonal events is smaller (<1°C) for June - July period and larger (1.5° to 2°C) during September, the latter due to a thinner mixed layer. To examine the effect of salinity on ISV, the model was run by artificially increasing the salinity (NORR run) and by decreasing it (MAHA10 run). In NORR, both rainfall and river discharge were switched off and in MAHA10 the discharge by river Mahanadi was increased tenfold. The spatial pattern of ISV as well as its periodicity was similar in QR, NORR and MAHA10. The ISV was stronger in NORR and weaker in MAHA10, compared to QR. In NORR, both intraseasonal warming and cooling were higher than in QR, the former due to reduced air-sea heat loss as the mean SST was lower, and the latter due to enhanced subsurface processes resulting from weaker stratification. In MAHA10, both warming and cooling were lower than in QR, the former due to higher air-sea heat loss owing to higher mean SST, and the latter due to weak subsurface processes resulting from stronger stratification. These model experiments suggest that salinity effects are crucial in determining amplitudes of intraseasonal SST variations in the BoB.