Indian summer monsoon (ISM) is an important component of the tropical climate system, known for its regular seasonality and abundance of rainfall over the country. The droughts and floods associated with the year-to-year variation of the average seasonal rainfall have devastating effect on people, agriculture and economy of this region. The demand for prediction of seasonal monsoon rainfall, therefore, is overwhelming. A number of attempts to predict the seasonal mean monsoon have been made over a century, but neither dynamical nor empirical models provide skillful forecasts of the extremes of the monsoon such as the unprecedented drought of 2002. This study investigates the problems and prospects of extended range monsoon prediction. An evaluation of the potential predictability of the ISM with the aid of an ensemble of Atmospheric General Circulation Model (AGCM) simulations indicates that the interannual variability (IAV) of ISM is contributed equally by the slow boundary forcing (‘externally’ forced variability) and the inherent climate noise (‘internal’ variability) in the atmosphere. Success in predicting the ISM would depend on our ability to extract the predictable signal from a background of noise of comparable amplitude. This would be possible only if the ‘external’ variability is separable from the ‘internal’ variability. A serious effort has been made to understand and isolate the sea surface temperature (SST) forced component of ISM variability that is not strongly influenced by the ‘internal’ variability. In addition, we have investigated to unravel the mechanism of generation of ‘internal’ IAV so that the method of isolating it from forced variability may be found. Since the primary forcing mechanism of the monsoon is the large-scale meridional gradient of deep tropospheric heat sources, large-scale changes in tropospheric temperature (TT) due to the boundary forcing can induce interannual variations of the timing and duration of the monsoon season. The concept of interannually varying monsoon season is introduced here, with the onset and withdrawal of monsoon definitions based on the reversal of meridional gradient of TT between north and south. This large scale definition of the monsoon season is representative of the planetary scale influence of the El Ni˜no Southern Oscillation (ENSO) on monsoon through the modification of TT and the cross equatorial pressure gradient over the ISM region. A sig- nificant relationship between ENSO and monsoon, that has remained steady over the decades, is discovered by which an El Ni˜no (La Ni˜na) delays (advances) the onset, advances (delays) the withdrawal and suppresses (enhances) the strength of the monsoon. The integral effect of the meridional gradient of TT from the onset to withdrawal proves to be a useful index of seasonal monsoon which isolates the boundary forced signal from the influence of internal variations that has remained steady even in the recent decades. However, consistent with the estimates of potential predictability, the boundary forced variability isolated with the above definitions explains only about 50% of the total interannual variability of ISM. Detailed diagnostics of the onset and withdrawal processes are performed to understand how the ENSO forcing modifies the onset and withdrawal, and thus the seasonal mean monsoon. It is found that during an El Ni˜no, the onset is delayed due to the enhanced adiabatic subsidence that inhibits vertical mixing of sensible heating from the warm landmass during pre-monsoon months, and the withdrawal is advanced due to the horizontal advective cooling. This link between ENSO and monsoon is realized through the advective processes associated with the stationary waves in the upper troposphere set up by the tropical ENSO heating. The remaining 50% of the monsoon IAV is governed by internal processes. To unravel the mechanism of the generation of internal IAV, we perform another set of AGCM simulations, forced with climatological monthly mean SSTs, to extract the pure internal IAV. We find that the spatial structure of the intraseasonal oscillations (ISOs) in these simulations has 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 broadband nature of the ISO spectrum, allowing the intraseasonal 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 relationship is a manifestation of the binomial character of the rainfall time series. The remaining part of IAV may arise due to the complex land-surface processes, scale interactions, etc. We also find that the ISOs over the ISM region are not significantly modulated by the Pacific and Indian Ocean SST variations. Thus, even with a perfect prediction of SST, only about 50% of the observed IAV of ISM could be predicted with the best model in forced mode. Even so, prediction of all India rainfall (AIR) representing the average conditions of the whole country and the season may not always serve the purposes of monsoon forecasting. One reason is the large inhomogeneities in the rainfall distribution during a normal seasonal monsoon. Agriculture and hydrological sector could benefit more if provided with regional scale forecasts of active/break spells 2-3 weeks ahead. Therefore, we advocate an alternative strategy to the seasonal prediction. Here, we present a method to estimate the potential predictability of active and break conditions from daily rainfall and circulation from observations for the recent 24 years. We discover that transitions from break to active conditions are much more chaotic than those from active to break, a fundamental property of the monsoon ISOs. The potential predictability limit of monsoon breaks (∼20 days) is significantly higher than that of the active conditions (∼10 days). An empirical real- time forecasting strategy to predict the sub-seasonal variations of monsoon up to 4 pentads (20 days) in advance is developed. The method is physically based, with the consideration that the large-scale spatial patterns and slow evolution of monsoon intraseasonal variations possess some similarity in their evolutions from one event to the other. This analog method is applied on NOAA outgoing longwave radiation (OLR) pentad mean data which is available on a near real time basis. The elimination of high frequency variability and the use of spatial and temporal analogs produces high and useful skill of predictions over the central and northern Indian region for a lead-time of 4-5 pentads. An important feature of this method is that, unlike other empirical methods to forecast monsoon ISOs, this uses minimal time filtering to avoid any possible end-point effects, and hence it has immense potential for real-time applications.

Using a 50-year (1951-2000) gridded (1-degree) daily rainfall data set over the Indian land region, we study two main aspects of the Indian monsoon. The first aspect deals with the changing nature of the rainfall annual cycle. This, to our knowledge, is the first attempt at studying the changing behaviour of the Indian monsoon rainfall annual cycle in a systematic way. The annual cycle is defined as a combination of the first few Fourier harmonics of daily rainfall. We then identify five attributes of the annual cycle for each year and location (grid): (a) the day of maximum intensity (peak day); (b) maximum intensity (peak value); (c) beginning; (d) end; and (e) duration of the annual cycle. An extensive statistical analysis of these five attributes over the central Indian region (16.5 – 26.5N; 74.5 – 86.5E) shows that the probability distributions of all attributes, barring the peak value, show a significant change in the last 25 years (1976-2000) compared to the first 25 years (1951-1975). The second issue addressed in this thesis deals with the behaviour of the intraseasonal oscillations in flood and drought years. Previous studies on this issue have been limited to only specific flood or drought years. Our analysis confirms earlier findings such as the northwestward propagation of the 10-20 day ISO. However, we also find, for the first time, based on 9 flood and 9 drought years, that the 20-60 day has an eastward propagation during drought years and remains stationary in flood years. The analysis is primarily statistical in nature, and providing a physical explanation for some of our findings is beyond the scope of our work. Finally, it is worth noting here that without the long-term gridded data, it would have been difficult to assess coherent changes over a large region and long time-period.

This thesis documents a collection of experimental investigations in which a differentially-heated annulus was used to investigate the effects of topography on the atmospheric and oceanic circulation. To this end a number of experiments were devised, each using a different topographic base to study a different aspect of the impact of topography, motivated by the most notable outstanding questions found in a review of the literature, namely exploring the effects of topographic resonance, blocking via partial barriers, and azimuthally differential-heating via thermal topography. First of all, whilst employing sinusoidal wavenumber-3 topography to extend the experimental parameter space of a similar study, namely Read and Risch (2011), a new regime within a region of structural vacillation was encountered. Denoted as the ‘stationary-transition’ regime, it featured periodic oscillations between a dominant stationary wavenumber-3 flow and axisymmetric or chaotic flow. An investigation into topographic resonance followed, keeping the wavenumber-3 base, but with a sloped lid to add a beta effect to the annulus. This acted to increase the occurrence of stationary waves, along with the ‘stationary-transition’ regime, which was discovered to be a near-resonant region where nonlinear topographic resonant instability led to a 23 to 42 ‘day’ oscillatory structure. The base was then replaced with an isolated ridge, forming a partial barrier to study the difference between blocked and unblocked flow. The topography was found to impact the circulation at a level much higher than its own peak, causing a unique flow structure when the drifting flow and the topography interacted in the form of an ‘interference’ regime at low Taylor Numbers, as well as forming an erratic ‘irregular’ regime at higher Taylor Numbers. Lastly, this isolated ridge was replaced by flat heating elements covering the same azimuthal extent, in order to observe whether thermal topography could be comparable to mechanical topography. These azimuthally-varying heating experiments produced much the same results as the partial barriers study, despite the lack of a physical peak or bottom-trapped waves, suggesting that blocking is independent of these activities. Evidence of resonant wave-triads was noted in all experiments, though the component wavenumbers of the wave-triads and their impact on the flow was found to depend on the topography in question.

In this thesis we examine intra-seasonal oscillations (ISO) in the aqua-planet setup of the Community Atmospheric Model (CAM) version 5.1, mainly based on July and January climatological sea surface temperature (SST). We investigate mainly two questions -what should be the SST distribution for the existence of (a) northward moving ISO in summer, and (b) eastward moving MJO-like modes in winter. In the first part of the thesis we discuss the northward propagation. A series of experiments were performed with zonally symmetric and asymmetric SST distributions. The basic lower boundary condition is specified from zonally averaged observed July and January SST. The zonally symmetric July SST experiment produced an inter tropical convergence zone (ITCZ) on both sides of the equator. Poleward movement is not clear, and it is confined to the region between the double ITCZ. In July, the Bay of Bengal (BOB) and West Pacific SST is high compared to the rest of the northern tropics. When we impose a zonally asymmetric SST structure with warm SST spanning about 80 of longitude, the model shows a monsoon-like circulation, and some northward propagating convective events. Analysis of these events shows that two adjacent cells with cyclonic and anticyclonic vorticity are created over the warm SST anomaly and to the west. The propagation occurs due to the convective region drawn north in the convergence zone between these vortices. Zonally propagating Madden-Julian oscillations (MJO) are discussed in the second part of the thesis. All the experiments in this part are based on the zonally symmetric SST. The zonally symmetric January SST configuration gives an MJO-like mode, with zonal wave number 1 and a period of 40-90 days. The SST structure has a nearly meridionally symmetric structure, with local SST maxima on either side of the equator, and a small dip in the equatorial region. If we replace this dip with an SST maximum, the time-scale of MJO becomes significantly smaller (20-40 days). The implication is that an SST maximum in the equatorial region reduces the strength of MJO, and a flat SST profile in the equatorial region is required for more energetic of MJO. This result was tested and found to be valid in a series of further experiments.

In this thesis, we have presented a systematic analysis of the change of cloud properties due to variation in aerosol concentration over Indian region using satellite observations, and Weather Research and Forecasting Model coupled with Chemistry (WRF-Chem) simulations. The Tropical Rainfall Measurement Mission (TRMM) based Microwave Imager (TMI) estimates (2A12) have been used to compare and contrast the characteristics of cloud liquid water and ice over the Indian land region and the surrounding oceans, during the pre-monsoon (May) and monsoon (June–September) seasons. Based on the spatial homogeneity of rainfall, we have selected five regions for our study (three over ocean, two over land). In general, we find that the mean cloud liquid water and cloud ice content of land and oceanic regions are different, with the ocean regions showing higher amount of CLW. A comparison across the ocean regions suggests that the cloud liquid water over the or graphically influenced Arabian Sea (close to the Indian west coast) behaves differently from the cloud liquid water over a trapped ocean (Bay of Bengal) or an open ocean (Equatorial Indian Ocean). Specifically, the Arabian Sea region shows higher liquid water for a lower range of rainfall, whereas the Bay of Bengal and the Equatorial Indian Ocean show higher liquid water for a higher range of rainfall. Apart from geographic differences, we also documented seasonal differences by comparing cloud liquid water profiles between monsoon and pre-monsoon periods, as well as between early and peak phases of the monsoon. We find that the cloud liquid water during the lean periods of rainfall (May or June) is higher than during the peak and late monsoon season (July-September) for raining clouds over central India. However, this is not true over the ocean. As active and break phases are important signatures of the monsoon progression, we also analyzed the differences in cloud liquid water during various phases of the monsoon, namely, active, break, active-to-break (a2b) and break-to-active (b2a) transition phases. We find that the cloud liquid water content during the b2a transition phase is significantly higher than that during the a2b transition phase over central India. We speculate that this could be attributed to higher amount of aerosol loading over this region during the break phase. We lend credence to this aerosol-liquid water/rain association by comparing the central Indian cloud liquid water with Southeast Asia (where the aerosol loading is significantly smaller) and find that in the latter region, there are no significant differences in cloud liquid water during the different phases of their monsoon. The second part of our study involves evaluating the ability of the Weather Research and Forecasting Model coupled with Chemistry (WRF-Chem) to simulate the observed variation of cloud liquid water and rain efficiency. We have used no chemistry option, and the model was run with constant aerosol concentration. The model simulations (at 4.5 km resolution) are done for the month of June–July 2004 since this period was particularly favorable for the study of an active–break cycle of the monsoon. We first evaluate the sensitivity of the model to different parameterizations (microphysical, boundary layer, land surface) on the simulation of rain over central India and Bay of Bengal. This is done to identify an “optimal” combination of parameterizations which reproduces the best correlation with observed rain over these regions. In this default configuration (control run), where the aerosol concentration is kept constant throughout the simulation period, the model is not able to reproduce the observed variations of cloud liquid water during the different phases of an active-break cycle. To this end, we proceeded to modify the model by developing an aerosol-rain relation, using Aerosol Robotic Network (AERONET) and TRMM 3B42 data that realistically captures the variation of aerosol with rain. It is worth highlighting here that our goal was to primarily isolate the indirect effect of aerosols in determining the observed changes in cloud liquid water (CLW) during the active-break phases of the Indian monsoon, without getting into the complexity of a full chemistry model such as that incorporated in WRF-Chem. Moreover, the proposed modification (modified run) is necessitated by the lack of realistic emission estimates over the Indian region as well as the presence of inherent biases in monsoon simulation in WRF. The main differences we find between the modified and control simulations is in the mean as well as spatial variability of CLW. We find that the proposed modification (i.e., rate of change of aerosol concentration as a function of rain rate) leads to a realistic variation in the CLW during the active-break cycle of Indian monsoon. Specifically, the peak value of CLW in the b2a (a2b) phase is larger (smaller) in the modified as compared to the control run. These results indicate a stronger change in CLW amount in the upper levels between the two transition phases in the modified scheme as compared to the control simulation. More significantly, we also observe a change in sign at the lower levels of the atmosphere, i.e., from a strong positive difference in the control run to a negative difference in the modified simulation, similar to that observed. Additionally, we investigated the impact of the proposed modification, via CLW changes, on cloud coverage, size of clouds and their spatial variability. We find that the transformation of optically thin clouds to thick clouds during the break phase was associated with larger cloud size in modified compared to the control simulation. Moreover, the higher rate of decay of the spatial variability of CLW with grid resolution, using the modified scheme, suggests that clusters of larger clouds are more in the modified compared to control simulation. Taken together, the interactive aerosol loading proposed in this thesis yields model simulations that better mimic the observed CLW variability between the transition phases.