Holocene precipitation change in different monsoon sub-regions (time-slices and transient data) simulated by different global climate models, with links to model results

The recently proposed global monsoon hypothesis interprets monsoon systems as part of one global-scale atmospheric overturning circulation, implying a connection between the regional monsoon systems and an in-phase behaviour of all northern hemispheric monsoons on annual timescales (Trenberth et al., 2000). Whether this concept can be applied to past climates and variability on longer timescales is still under debate, because the monsoon systems exhibit different regional characteristics such as different seasonality (i.e. onset, peak, and withdrawal). To investigate the interconnection of different monsoon systems during the pre-industrial Holocene, five transient global climate model simulations have been analysed with respect to the rainfall trend and variability in different sub-domains of the Afro-Asian monsoon region. Our analysis suggests that on millennial timescales with varying orbital forcing, the monsoons do not behave as a tightly connected global system. According to the models, the Indian and North African monsoons are coupled, showing similar rainfall trend and moderate correlation in rainfall variability in all models. The East Asian monsoon changes independently during the Holocene. The dissimilarities in the seasonality of the monsoon sub-systems lead to a stronger response of the North African and Indian monsoon systems to the Holocene insolation forcing than of the East Asian monsoon and affect the seasonal distribution of Holocene rainfall variations. Within the Indian and North African monsoon domain, precipitation solely changes during the summer months, showing a decreasing Holocene precipitation trend. In the East Asian monsoon region, the precipitation signal is determined by an increasing precipitation trend during spring and a decreasing precipitation change during summer, partly balancing each other. A synthesis of reconstructions and the model results do not reveal an impact of the different seasonality on the timing of the Holocene rainfall optimum in the different sub-monsoon systems. They rather indicate locally inhomogeneous rainfall changes and show, that single palaeo-records should not be used to characterise the rainfall change and monsoon evolution for entire monsoon sub-systems.

The physics of ENSO-monsoon connection

The physical mechanism through which Ei-Nino and Southern Oscillation (ENSO) tends to produce deficient precipitation over Indian continent is investigated using both observations as well as a general circulation model. Both analysis of observations and atmospheric general circulation model (AGCM) study show that the planetary scale response associated with ENSO primarily influences the equatorial Indian Ocean region. Through this interaction it tends to favour the equatorial heat source, enhance precipitation over the equatorial Indian Ocean and indirectly cause a decrease in continental precipitation through induced subsidence. This situation is further complicated by the fact the regional tropospheric quasi biennial oscillation (QBO) has a bimodal structure over this region with large amplitude over the Indian continent. While the ENSO response has a quasi-four year periodicity and tends peak during beginning of the calendar year, the QBO mode tends to peak during northern summer. Thus, the QBO mode exerts a stronger influence on the interannual variability of the monsoon. The strength of the Indian monsoon in a given year depends on the combined effect of the ENSO and the QBO mode. Sines the two oscillations have disparate time scales, exact phase information of the two modes during northern summer is important in determining the Indian summer monsoon. The physical mechanism of the interannual variations of the Indian monsoon precipitation associated with ENSO presented here is similar to the physical process that cause intraseasonal 'active', 'break' oscillations of the monsoon.

Ocean atmosphere coupling over monsoon regions

In monsoon regions, the seasonal migration of the intertropical convergence zone (ITCZ) is manifested as a seasonal reversal of winds. Most of the summer monsoon rainfall over India occurs owing to synoptic and large-scale convection associated with the continental ITCZ (Fig. 1). We have investigated the interaction between these large-scale convective systems and the ocean over which they are generated1â3, concentrating on the relationship between organized convection over the Indian Ocean and sea surface temperature (SST). We report here that on a monthly basis the degree of cloudiness correlates well with SST for the relatively colder oceans, but when SST is maintained above 28 °C it ceases to be an important factor in determining the variability of cloudiness. Over the major regions of convection east of 70°E, which are warm year after year, the observed cloudiness cannot be correlated with variations in SST.

On the Maximum Cloud Zone and the ITCZ over Indian, Longitudes during the Southwest Monsoon

An investigation is presented of the daily variation of the maximum cloud zone (MCZ) and the 7W mb trough in the Northern Hemisphere over the Indian longitudes 70–90°E during April–October for 1973–77. It is found that during June–September there are two favorable locations for a MCZ over these longitudes–on a majority of days the MCZ is present in the monsoon zone north of 15°N, and often a secondary MCZ occurs in the equatorial region (0–10°N). The monsoon MCZ gets established by northward movement of the MCZ occurring over the equatorial Indian ocean in April and May. The secondary MCZ appears intermittently, and is characterized by long spells of persistence only when the monsoon MCZ is absent. In each of the seasons studied, the MCZ temporarily disappeared from the mean summer monsoon location (15–28°N) about four weeks after it was established near the beginning of July. It is reestablished by the northward movement of the secondary MCZ, which becomes active during the absence of the monsoon MCZ, in a manner strikingly similar to that observed in the spring to summer transition. A break in monsoon conditions prevails just prior to the temporary disappearance of the monsoon MCZ. Thus we conclude that the monsoon MCZ cannot survive for longer than a month without reestablishment by the secondary MCZ. Possible underlying mechanisms are also discussed.

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.

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.

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.

On breaks of the Indian monsoon

For over a century, the term break has been used for spells in which the rainfall over the Indian monsoon zone is interrupted. The phenomenon of 'break monsoon' is of great interest because long intense breaks are often associated with poor monsoon seasons. Such breaks have distinct circulation characteristics (heat trough type circulation) and have a large impact on rainfed agriculture. Although interruption of the monsoon rainfall is considered to be the most important feature of the break monsoon, traditionally breaks have been identified on the basis of the surface pressure and wind patterns over the Indian region. We have defined breaks (and active spells) on the basis of rainfall over the monsoon zone. The rainfall criteria are chosen so as to ensure a large overlap with the traditional breaks documented by Ramamurthy (1969) and De et al (1998). We have identified these rainbreaks for 1901-89. We have also identified active spells on the basis of rainfall over the Indian monsoon zone. We have shown that the all-India summer monsoon rainfall is significantly negatively correlated with the number of rainbreak days (correlation coefficient -0.56) and significantly positively correlated with the number of active days (correlation coefficient 0.47). Thus the interannual variation of the all-India summer monsoon rainfall is shown to be related to the number of days of rainbreaks and active spells identified here. There have been several studies of breaks (and also active spells in several cases) identified on the basis of different criteria over regions differing in spatial scales (e.g., Webster et al 1998; Krishnan et al 2000; Goswami and Mohan 2000; and Annamalai and Slingo 2001). We find that there is considerable overlap between the rainbreaks we have identified and breaks based on the traditional definition. There is some overlap with the breaks identified by Krishnan et al (2000) but little overlap with breaks identified by Webster et al (1998). Further, there are three or four active-break cycles in a season according to Webster et al (1998) which implies a time scale of about 40 days for which Goswami and Mohan (2000), and Annamalai and Slingo'(2001) have studied breaks and active minus break fluctuations. On the other hand, neither the traditional breaks (Ramamurthy 1969; and De et al 1998) nor the rainbreaks occur every year. This suggests that the 'breaks' in these studies axe weak spells of the intraseasonal variation of the monsoon, which occur every year. We have derived the OLR and circulation patterns associated with rainbreaks and active spells and compared them with the patterns associated with breaks/active minus break spells from these studies. Inspite of differences in the patterns over the Indian region, there is one feature which is seen in the OLR anomaly patterns of breaks identified on the basis of different criteria as well as the rainbreaks identified in this paper viz., a quadrapole over the Asia-west Pacific region arising from anomalies opposite (same) in sign to those over the Indian region occurring over the equatorial Indian Ocean and northern tropical (equatorial) parts of the west Pacific. Thus it appears that this quadrapole is a basic feature of weak spells of the intraseasonal,variation over the Asia-west Pacific region. Since the rainbreaks are intense weak spells, this basic feature is also seen in the composite patterns of these breaks. We find that rainbreaks (active spells) are also associated with negative (positive) anomalies over a part of the cast Pacific suggesting that the convection over the Indian region is linked to that over the east Pacific not only on the interannual scale (as evinced by the link between the Indian summer monsoon rainfall and ENSO) but on the intraseasonal scale as well.

Simulating the daily discharge of the Mandovi River, west coast of India

A hydrological modelling framework was assembled to simulate the daily discharge of the Mandovi River on the Indian west coast. Approximately 90% of the west-coast rainfall, and therefore discharge, occurs during the summer monsoon (June-September), with a peak during July-August. The modelling framework consisted of a digital elevation model (DEM) called GLOBE, a hydrological routing algorithm, the Terrestrial Hydrological Model with Biogeochemistry (THMB), an algorithm to map the rainfall recorded by sparse rain-gauges to the model grid, and a modified Soil Conservation Service Curve Number (SCS-CN) method. A series of discharge simulations (with and without the SCS method) was carried out. The best simulation was obtained after incorporating spatio-temporal variability in the SCS parameters, which was achieved by an objective division of the season into five regimes: the lean season, monsoon onset, peak monsoon, end-monsoon, and post-monsoon. A novel attempt was made to incorporate objectively the different regimes encountered before, during and after the Indian monsoon, into a hydrological modelling framework. The strength of our method lies in the low demand it makes on hydrological data. Apart from information on the average soil type in a region, the entire parameterization is built on the basis of the rainfall that is used to force the model. That the model does not need to be calibrated separately for each river is important, because most of the Indian west-coast basins are ungauged. Hence, even though the model has been validated only for the Mandovi basin, its potential region of application is considerable. In the context of the Prediction in Ungauged Basins (PUB) framework, the potential of the proposed approach is significant, because the discharge of these (ungauged) rivers into the eastern Arabian Sea is not small, making them an important element of the local climate system.

Monsoon sensitivity to aerosol direct radiative forcing in the community atmosphere model

Aerosol forcing remains a dominant uncertainty in climate studies. The impact of aerosol direct radiative forcing on Indian monsoon is extremely complex and is strongly dependent on the model, aerosol distribution and characteristics specified in the model, modelling strategy employed as well as on spatial and temporal scales. The present study investigates (i) the aerosol direct radiative forcing impact on mean Indian summer monsoon when a combination of quasi-realistic mean annual cycles of scattering and absorbing aerosols derived from an aerosol transport model constrained with satellite observed Aerosol Optical Depth (AOD) is prescribed, (ii) the dominant feedback mechanism behind the simulated impact of all-aerosol direct radiative forcing on monsoon and (iii) the relative impacts of absorbing and scattering aerosols on mean Indian summer monsoon. We have used CAM3, an atmospheric GCM (AGCM) that has a comprehensive treatment of the aerosol-radiation interaction. This AGCM has been used to perform climate simulations with three different representations of aerosol direct radiative forcing due to the total, scattering aerosols and black carbon aerosols. We have also conducted experiments without any aerosol forcing. Aerosol direct impact due to scattering aerosols causes significant reduction in summer monsoon precipitation over India with a tendency for southward shift of Tropical Convergence Zones (TCZs) over the Indian region. Aerosol forcing reduces surface solar absorption over the primary rainbelt region of India and reduces the surface and lower tropospheric temperatures. Concurrent warming of the lower atmosphere over the warm oceanic region in the south reduces the land-ocean temperature contrast and weakens the monsoon overturning circulation and the advection of moisture into the landmass. This increases atmospheric convective stability, and decreases convection, clouds, precipitation and associated latent heat release. Our analysis reveals a defining negative moisture-advection feedback that acts as an internal damping mechanism spinning down the regional hydrological cycle and leading to significant circulation changes in response to external radiative forcing perturbations. When total aerosol loading (both absorbing and scattering aerosols) is prescribed, dust and black carbon aerosols are found to cause significant atmospheric heating over the monsoon region but the aerosol-induced weakening of meridional lower tropospheric temperature gradient (leading to weaker summer monsoon rainfall) more than offsets the increase in summer-time rainfall resulting from the atmospheric heating effect of absorbing aerosols, leading to a net decrease of summer monsoon rainfall. Further, we have carried out climate simulations with globally constant AODs and also with the constant AODs over the extended Indian region replaced by realistic AODs. Regional aerosol radiative forcing perturbations over the Indian region is found to have impact not only over the region of loading but over remote tropical regions as well. This warrants the need to prescribe realistic aerosol properties in strategic regions such as India in order to accurately assess the aerosol impact.

Space-Time Scales of Northward Propagation of Convection during Boreal Summer

This study uses precipitation estimates from the Tropical Rainfall Measuring Mission to quantify the spatial and temporal scales of northward propagation of convection over the Indian monsoon region during boreal summer. Propagating modes of convective systems in the intraseasonal time scales such as the Madden-Julian oscillation can interact with the intertropical convergence zone and bring active and break spells of the Indian summer monsoon. Wavelet analysis was used to quantify the spatial extent (scale) and center of these propagating convective bands, as well as the time period associated with different spatial scales. Results presented here suggest that during a good monsoon year the spatial scale of this oscillation is about 30 degrees centered around 10 degrees N. During weak monsoon years, the scale of propagation decreases and the center shifts farther south closer to the equator. A strong linear relationship is obtained between the center/scale of convective wave bands and intensity of monsoon precipitation over Indian land on the interannual time scale. Moreover, the spatial scale and its center during the break monsoon were found to be similar to an overall weak monsoon year. Based on this analysis, a new index is proposed to quantify the spatial scales associated with propagating convective bands. This automated wavelet-based technique developed here can be used to study meridional propagation of convection in a large volume of datasets from observations and model simulations. The information so obtained can be related to the interannual and intraseasonal variation of Indian monsoon precipitation.

Monsoon circulation interaction with Western Ghats orography under changing climate - Projection by a 20-km mesh AGCM

In this study, the authors have investigated the likely future changes in the summer monsoon over the Western Ghats (WG) orographic region of India in response to global warming, using time-slice simulations of an ultra high-resolution global climate model and climate datasets of recent past. The model with approximately 20-km mesh horizontal resolution resolves orographic features on finer spatial scales leading to a quasi-realistic simulation of the spatial distribution of the present-day summer monsoon rainfall over India and trends in monsoon rainfall over the west coast of India. As a result, a higher degree of confidence appears to emerge in many aspects of the 20-km model simulation, and therefore, we can have better confidence in the validity of the model prediction of future changes in the climate over WG mountains. Our analysis suggests that the summer mean rainfall and the vertical velocities over the orographic regions of Western Ghats have significantly weakened during the recent past and the model simulates these features realistically in the present-day climate simulation. Under future climate scenario, by the end of the twenty-first century, the model projects reduced orographic precipitation over the narrow Western Ghats south of 16A degrees N that is found to be associated with drastic reduction in the southwesterly winds and moisture transport into the region, weakening of the summer mean meridional circulation and diminished vertical velocities. We show that this is due to larger upper tropospheric warming relative to the surface and lower levels, which decreases the lapse rate causing an increase in vertical moist static stability (which in turn inhibits vertical ascent) in response to global warming. Increased stability that weakens vertical velocities leads to reduction in large-scale precipitation which is found to be the major contributor to summer mean rainfall over WG orographic region. This is further corroborated by a significant decrease in the frequency of moderate-to-heavy rainfall days over WG which is a typical manifestation of the decrease in large-scale precipitation over this region. Thus, the drastic reduction of vertical ascent and weakening of circulation due to `upper tropospheric warming effect' predominates over the `moisture build-up effect' in reducing the rainfall over this narrow orographic region. This analysis illustrates that monsoon rainfall over mountainous regions is strongly controlled by processes and parameterized physics which need to be resolved with adequately high resolution for accurate assessment of local and regional-scale climate change.

Role of orography in modulating space-time scales of convection over South Asia

Propagation of convective systems in the meridional direction during boreal summer is responsible for active and break phases of monsoon over south Asia. This region is unique in the world in its characteristics of monsoon variability and is in close proximity of mountains like the Himalayas. Here, using an atmospheric general circulation model, we try to understand the role of orography in determining spatial and temporal scales of these convective systems. Absence of orography (noGlOrog) decreased the simulated seasonal mean precipitation over India by 23 % due to delay in onset by about a month vis-a-vis the full-mountain case. In noGlOrog, poleward propagations were absent during the delayed period prior to onset. Post-onset, both simulations had similar patterns of poleward propagations. The spatial and temporal scales of propagating clouds bands were determined using wavelet analysis. These scales were found to be different in full-mountain and no-mountain experiments in June-July. However, after the onset of monsoon in noGlOrog, these scales become similar to that with orography. Simulations with two different sets of convection schemes confirmed this result. Further analysis shows that the absence (presence) of meridional propagations during early (late) phase of summer monsoon in noGlOrog was associated with weaker (stronger) vertical shear of zonal wind over south Asia. Our study shows that orography plays a major role in determining the time of onset over the Indian region. However, after onset, basic characteristics of propagating convective systems and therefore the monthly precipitation over India, are less sensitive to the presence of orography and are modulated by moist convective processes.

Local and remote impacts of direct aerosol forcing on Asian monsoon

The impact of heating by black carbon aerosols on Indian summer monsoon has remained inconclusive. Some investigators have predicted that black carbon aerosols reduce monsoon rainfall while others have argued that it will increase monsoon rainfall. These conclusions have been based on local influence of aerosols on the radiative fluxes. The impact of aerosol-like heating in one region on the rainfall in a remote region has not been examined in detail. Here, using an atmospheric general circulation model, it has been shown that remote influence of aerosol-like heating can be as important as local influence on Indian summer monsoon. Precipitation in northern Arabian Sea and north-west Indian region increased by 16% in June to July when aerosol-like heating were present globally. The corresponding increase in precipitation due to presence of aerosol-like heating only over South Asia (local impact) and East Asia (remote impact) were 28 and 13%, respectively. This enhancement in precipitation was due to destabilization of the atmosphere in pre-monsoon season that affected subsequent convection. Moreover, pre-monsoon heating of the lower troposphere changed the circulation substantially that enabled influx of more moisture over certain regions and reduced the moist static stability of the atmosphere. It has been shown that regional aerosol heating can have large impact on the phase of upper tropospheric Rossby wave in pre-monsoon season, which acts as a primary mechanism behind teleconnection and leads to the change in precipitation during monsoon season. These results demonstrate that changes in aerosol in one region can influence the precipitation in a remote region through changes in circulation.

Decreasing intensity of monsoon low-frequency intraseasonal variability over India

Understanding the changing nature of the intraseasonal oscillatory (ISO) modes of Indian summer monsoon manifested by active and break phase, and their association with extreme rainfall events are necessary for probabilistic estimation of flood-related risks in a warming climate. Here, using ground-based observed rainfall, we define an index to measure the strength of monsoon ISOs and show that the relative strength of the northward-propagating low-frequency ISO (20-60 days) modes have had a significant decreasing trend during the past six decades, possibly attributed to the weakening of large-scale circulation in the region during monsoon season. This reduction is compensated by a gain in synoptic-scale (3-9 days) variability. The decrease in low-frequency ISO variability is associated with a significant decreasing trend in the percentage of extreme events during the active phase of the monsoon. However, this decrease is balanced by significant increasing trends in the percentage of extreme events in the break and transition phases. We also find a significant rise in the occurrence of extremes during early and late monsoon months, mainly over eastern coastal regions. Our study highlights the redistribution of rainfall intensity among periodic (low-frequency) and non-periodic (extreme) modes in a changing climate scenario.