Interannual variation of Indian summer (June-September: JJAS) monsoon rainfall (ISMR) depends on its relative intensity during early (June-July: JJ; contribution 52%) and late (August-September: AS; contribution 49%) phases. Apart from variations in sea surface temperature (SST), the primary reasons behind the variability during JJ and AS can be very different due to change in climatic conditions on account of post-onset processes. Here, using a high resolution general circulation model with seasonally varying climatological SST, mechanisms those govern the intensity of rainfall during JJ and AS are investigated. There is no significant relation-ship between intensity of precipitation over Indian region in JJ and AS. Moreover, the factors determining early monsoon (JJ) precipitation are different than that for late monsoon (AS). In absence of interannual SST variation, pre-monsoon soil moisture do not play a significant role for the interannual variation of monsoon precipitation over India. A large scale oscillation of the ITCZ is noticed on interannual time scale spanning from around 60◦E to 150◦E that brings spatially coherent flood and drought over this region. Early monsoon precipitation has a larger dependency on spring snow depth over Eurasia and phase of the upper tropospheric Rossby wave in May. However, late monsoon precipitation over India is mainly governed by the intensity and time scale of the intraseasonally varying convective cloud bands. This study suggests that early monsoon (JJ) precipitation over Indian region is more correlated with pre-monsoon signatures of land-atmosphere parameters. However, in later parts after the onset (AS), the monsoon intensity is primarily driven by its internal dynamics and characteristics of intraseasonal oscillation.

This thesis examines the structure of the Tropical Easterly Jet (TEJ) in a General Circulation Model (GCM). The TEJ is observed only during the Indian summer monsoon period and is strongest during July and August. The jet structure simulated by an atmospheric GCM (CAM-3.1) in July has been compared with reanalysis data. The simulated TEJ was displaced westward by ~ 25◦ when compared to observations. The removal of orography had no impact on the jet structure. This demonstrated that the Tibetan Plateau did not play an important role in the location and structure of the jet. The changes in cumulus scheme in the GCM had a large influence on the location of the jet maxima. To examine the factors which control the location and structure of the jet, a series of experiments were conducted using an aqua-planet version of the model. The impact of different sea surface temperature (SST) profiles was studied. The rainfall in the GCM was primarily in the regions where the SST attained a maximum. By altering the location of SST maximum (and hence the rainfall maximum), the impact of location of rainfall maximum on the location and structure of the jet was studied. When the rainfall maximum was located close to the equator, it did not generate a strong jet but had an influence on the vertical structure of the jet. A large number of simulations were conducted with multiple rainfall maxima and the need for these was demonstrated since only then was the observed jet structure well simulated. Based on the simulations, it was concluded that the simulation of the TEJ by CAM-3.1 was unrealistic because of large unrealistic rainfall over Saudi Arabia in this GCM. Equatorial heating has been shown to be important to simulate proper jet structure. The zonal structure of the jet was also influenced by rainfall in the Pacific Ocean. Although the aqua-planet configuration of the CAM-3.1 GCM provided several useful insights, the simulation was not perfect on account of errors in the simulation of the temperature profile in the lower troposphere. An ideal-physics configuration of the GCM was used. This removed the cumulus physics and instead imposed the observed heating pro-files. Both upper tropospheric friction and radiative-convective atmospheric temperatures were required to simulate the TEJ. The problems with the simulation of structure in the jet exit region was corrected by using radiative-convective atmospheric temperatures that were qualitatively similar to those observed in northern hemisphere summer time. The ideal-physics configuration reconfirmed that the Saudi Arabian rainfall was responsible for the westward shift of the TEJ in the simulations. The ideal-physics simulations showed that the simple analytical model proposed by Gillin1980 was not suitable for the simulation of TEJ. The above the simulations indicate that a shift in the location of the jet is related to a shift in the rainfall pattern. Based on this insight one would expect that the jet location will be different in good and bad monsoon periods. This is indeed the case. In July 2002 the Indian monsoon failed after beginning well in June. In June the TEJ is consequently located west ward compared to July. The same situation prevails even in good and poor monsoon years. In a good monsoon year (July 1988) the jet maximum is located westward when compared to a bad monsoon year (July 2002). In this thesis we have clearly demonstrated the role of anomalous rainfall on the location of the TEJ. This thesis has shown that an accurate simulation of the TEJ depends upon the accurate simulation of various rainfall centers that act as multiple heat sources in the atmosphere. The rainfall in the equatorial region does not influence the strength of the TEJ but alters the vertical structure of the jet. The strength the jet is dependent on the intensity of rainfall and the latitudinal distance from the equator. The complex vertical structure of the jet is not simulated by simple analytical models of the jet.

In this thesis, we investigated modes of intraseasonal variability (ISV) observed in the Indian monsoon rainfall and how these modes modulate rainfall over India. We identified a decreasing trend in the intensity of low-frequency intraseasonal mode with increasing strength in synoptic variability over India. We also made an attempt to understand the reason for these observed trends using numerical simulations. In the first part of the thesis, satellite rainfall estimates are used to understand the spatiotem-poral structures of convection in the intraseasonal timescale and their intensity during boreal sum-mer over south Asia. Two dominant modes of variability with periodicities of 10–20-days (high-frequency) and 20–60-days (low-frequency) are found, with the latter strongly modulated by sea surface temperature. The 20–60-day mode shows northward propagation from the equatorial In-dian Ocean linked with eastward propagating modes of convective systems over the tropics. The 10–20-day mode shows a complex space-time structure with a northwestward propagating anoma-lous pattern emanating from the Indonesian coast. This pattern is found to be interacting with a structure emerging from higher latitudes propagating southeastwards. This could be related to ver-tical shear of zonal wind over northern India. The two modes exhibit variability in their intensity on the interannual time scale and contribute a significant amount to the daily rainfall variability in a season. The intensities of the 20–60-day and 10–20-day modes show significantly strong inverse and direct relationship, respectively, with the all-India June–September rainfall. This study also establishes that the probability of occurrence of substantial rainfall over central India increases significantly if the two intraseasonal modes simultaneously exhibit positive anomalies over the region. There also exists a phase-locking between the two modes. In the second part of the thesis, we investigated the changing nature of these intraseasonal modes over Indian region, and their association with extreme rainfall events using ground based observed rainfall. We found that the relative strength of the northward propagating 20–60-day mode has a significant decreasing trend during the past six decades, possibly attributed to the weakening of large-scale circulation in the region during monsoon. This reduction is compensated by a gain in synoptic-scale (3–9 days) variability. The decrease in the low-frequency ISV 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 a significant increasing trend in the percentage of extreme events in break phase. We also find a significant rise in occurrence of extremes during early- and late-monsoon months, mainly over the eastern coastal regions of India. We do not observe any significant trend in the high-frequency ISV. In the last part of the thesis, we used numerical simulations to understand the observed changes in the ISV features. Using the atmospheric component of a global climate model (GCM), we have performed two experiments: control experiment (CE) and heating experiment (HE). The CE is the default simulation for 10 years. In HE, we prescribed heating in the atmosphere in such a way that it mimics the conditions for extreme rainfall events as observed over central India during June– September. Heating is prescribed primarily during the break phase of the 20–60-day mode. This basically increases the number of extremes, majority of which are in break phase. The design of the experiment reflects the observed current scenario of increased extreme events during breaks. We found that the increased extreme events in the HE decreased the intensity of the 20–60-day mode over the Indian region. This reduction is associated with a reduction of rainfall in active phase and increase in the length of break phase. A reduction in the seasonal mean over India is also observed. The reduction of active phase rainfall is linked with an increased stability of the atmosphere over central India. Lastly, we propose a possible mechanism for the reduction of rainfall in active phase. We found that there is a significant reduction in the strength of the vertical easterly shear over the northern Indian region during break–active transition phase. This basically weakens the conditions for the growth of Rossby wave instability, thereby elongating break phase and reducing the rainfall intensity in the following active phase. This study highlights the redistribution of rainfall intensity among periodic (low-frequency) and non-periodic (extreme) modes in a changing climate scenario, which is further tested in a modeling study. The results presented in this thesis will provide a pathway to understand, using observations and numerical model simulations, the ISV and its relative contribution to the Indian summer monsoon. It can also be used for model evaluation.

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.

Smallholder farmers often have great entrepreneurial qualities that build on generations of experience. However, many farm management practices are poorly adapted to current climate change conditions. In order for farmers to understand the risks they are undertaking by following certain farming practices and to adapt accordingly, a decision support tool is being developed by researchers at TU Delft. The tool runs a socio-hydrological model, created in Python, in the back-end and provides farmer specific investment and profit data for different crops in the front-end. The aim of this study is to develop a risk assessment process that integrates hydro-climatic variability in the decision support tool, and to identify ways of communicating risk to smallholder farmers in Maharashtra, India. Two sources of variability were characterised based on a literature review of Indian farmers’ own risk perceptions; the untimely onset of the Indian Summer Monsoon and the frequency of dry spells. A sensitivity analysis was then carried out to investigate their respective effects on the farmers’ crop yields. The method proposed to evaluate these risks used a single variable, precipitation data, and a two-dimensional risk matrix to compound the two risk factors, over a time span of 14 years (2003-2016). However, the results indicate that it might be more beneficial to define dry spells in terms of crop water stress, instead of a precipitation threshold. This study also proposed a method for translating a cumulative distribution curve into a risk representation that is adapted for low-literacy users by combining numbers and text with graphics, color and voice descriptions. Ultimately, however, the usability of the tool cannot be determined solely through literature, but must involve the end-users in its design. / Småskaliga jordbrukare är goda entreprenörer som samlat på sig kunskaper och erfarenheter över flera generationer. Däremot är vissa metoder som jordbrukarna använder sig av idag för att förvalta sitt jordbruk inte anpassade till nutida klimatförändringar. För att jordbrukarna ska förstå riskerna som de åtar sig vid valet av dessa metoder försöker forskare vid TU Delft nu ta fram ett verktyg för att underlätta jordbrukares förmåga att ta självständiga men välgrundade beslut om sitt jordbruk. Verktyget är baserat på en socio-hydrologisk modell som är framtagen i Python och som förser specifika investerings- och inkomstdata för enskilda jordbrukare. Syftet med detta kandidatarbete är att bidra till verktyget genom att undersöka de hydroklimatiska risker som uppstår till följd av föränderliga och osäkra klimatologiska förhållanden för jordbrukare i delstaten Maharashtra, Indien. Två riskfaktorer karakteriserades baserat på en litteraturstudie om indiska jordbrukares riskuppfattningar: avvikelser i starten på den indiska sommarmonsunen och antal torrperioder under monsunsäsongen. Dessutom utfördes en känslighetsanalys för att undersöka om och hur den existerande modellens utdata av skörd påverkades av de valda riskfaktorerna. Monsunstarten och torrperioderna togs fram genom metoder som enbart använde historiska nederbördsdata över tidsperioden 2003-2016 och kombinerades sedan med hjälp av en tvådimensionell riskmatris. Resultaten visade att det fanns anledning att ifrågasätta hur torrperioderna definierades och att det kan vara mer fördelaktigt att undersöka vattenbrist för grödan, snarare än att enbart förlita sig på nederbördsdata. Vidare föreslog denna studie en metod för att översätta en kumulativ fördelningsfunktion till en grafisk riskframställning som är anpassad till användare med låg läskunnighet genom att kombinera siffror med text, grafik, färg och ljudförklaringar. I slutändan kan dock inte användbarheten av verktyget enbart avgöras utifrån litteratur, utan måste även inkludera återkoppling från slutanvändarna.

In the recent scenario of climate change, the natural variability and uncertainty associated with the hydrologic variables is of great concern to the community. This thesis opens up a new area of multi-disciplinary research. It is a promising field of research in hydrology and water resources that uses the information from the field of atmospheric science. A new way to identify and capture the variability and uncertainty associated with the hydrologic variables is established through this thesis. Assessment of hydroclimatic teleconnection for Indian subcontinent and its use in basin-scale hydrologic time series analysis and forecasting is the broad aim of this PhD thesis. The initial part of the thesis is devoted to investigate and establish the dependence of Indian summer monsoon rainfall (ISMR) on large-scale Oceanic-atmospheric circulation phenomena from tropical Pacific Ocean and Indian Ocean regions. El Niño-Southern Oscillation (ENSO) is the well established coupled Ocean-atmosphere mode of tropical Pacific Ocean whereas Indian Ocean Dipole (IOD) mode is the recently identified coupled Ocean-atmosphere mode of tropical Indian Ocean. Equatorial Indian Ocean Oscillation (EQUINOO) is known as the atmospheric component of IOD mode. The potential of ENSO and EQUINOO for predicting ISMR is investigated by Bayesian dynamic linear model (BDLM). A major advantage of this method is that, it is able to capture the dynamic nature of the cause-effect relationship between large-scale circulation information and hydrologic variables, which is quite expected in the climate change scenario. Another new method, proposed to capture the dependence between the teleconnected hydroclimatic variables is based on the theory of copula, which itself is quite new to the field of hydrology. The dependence of ISMR on ENSO and EQUINOO is captured and investigated for its potential use to predict the monthly variation of ISMR using the proposed method. The association of monthly variation of ISMR with the combined information of ENSO and EQUINOO, denoted by monthly composite index (MCI), is also investigated and established. The spatial variability of such association is also investigated. It is observed that MCI is significantly associated with monthly rainfall variation all over India, except over North-East (NE) India, where it is poor. Having established the hydroclimatic teleconnection at a comparatively larger scale, the hydroclimatic teleconnection for basin-scale hydrologic variables is then investigated and established. The association of large-scale atmospheric circulation with inflow during monsoon season into Hirakud reservoir, located in the state of Orissa in India, has been investigated. The strong predictive potential of the composite index of ENSO and EQUINOO is established for extreme inflow conditions. So the methodology of inflow prediction using the information of hydroclimatic teleconnection would be very suitable even for ungauged or poorly gauged watersheds as this approach does not use any information about the rainfall in the catchment. Recognizing the basin-scale hydroclimatic association with both ENSO and EQUINOO at seasonal scale, the information of hydroclimatic teleconnection is used for streamflow forecasting for the Mahanadi River basin in the state of Orissa, India, both at seasonal and monthly scale. It is established that the basin-scale streamflow is influenced by the large-scale atmospheric circulation phenomena. Information of streamflow from previous month(s) alone, as used in most of the traditional modeling approaches, is shown to be inadequate. It is successfully established that incorporation of large-scale atmospheric circulation information significantly improves the performance of prediction at monthly scale. Again, the prevailing conditions/characteristics of watershed are also important. Thus, consideration of both the information of previous streamflow and large-scale atmospheric circulations are important for basin-scale streamflow prediction at monthly time-scale. Adopting the developed approach of using the information of hydroclimatic teleconnection, hydrologic variables can be predicted with better accuracy which will be a very useful input for better management of water resources.