The scale dependency of variance of total water mixing ratio is explored by analyzing data from a general circulation model (GCM), a numerical weather prediction model (NWP), and large-eddy simulations (LESs). For clarification, direct numerical simulation (DNS) data are additionally included, but the focus is placed on defining a general scaling behavior for scales ranging from global down to cloud resolving. For this, appropriate power-law exponents are determined by calculating and approximating the power density spectrum. The large-scale models (GCM and NWP) show a consistent scaling with a power-law exponent of approximately 22. For the high-resolution LESs, the slope of the power density spectrum shows evidence of being somewhat steeper, although the estimates are more uncertain. Also the transition between resolved and parameterized scales in a current GCM is investigated. Neither a spectral gap nor a strong scale break is found, but a weak scale break at high wavenumbers cannot be excluded. The evaluation of the parameterized total water variance of a state-of-the-art statistical scheme shows that the scale dependency is underestimated by this parameterization. This study and the discovered general scaling behavior emphasize the need for a development of scale-dependent parameterizations.

This study explores the viability of parameter estimation in the comprehensive general circulation model ECHAM6 using ensemble Kalman filter data assimilation techniques. Four closure parameters of the cumulus-convection scheme are estimated using increasingly less idealized scenarios ranging from perfect-model experiments to the assimilation of conventional observations. Updated parameter values from experiments with real observations are used to assess the error of the model state on short 6 h forecasts and on climatological timescales. All parameters converge to their default values in single parameter perfect-model experiments. Estimating parameters simultaneously has a neutral effect on the success of the parameter estimation, but applying an imperfect model deteriorates the assimilation performance. With real observations, single parameter estimation generates the default parameter value in one case, converges to different parameter values in two cases, and diverges in the fourth case. The implementation of the two converging parameters influences the model state: Although the estimated parameter values lead to an overall error reduction on short timescales, the error of the model state increases on climatological timescales.

Documenter la variabilité climatique récente est nécessaire à la compréhension des mécanismes en jeu, associés au rôle du bilan de masse de l’Antarctique pour l’élévation du niveau des mers globale. Les enregistrements issus des carottes peu profondes d’Antarctique sont des données précieuses, complémentaires aux observations instrumentales et satellitaires, pour couvrir en continu l’ensemble du continent. Mesurés le long de ces carottes de glace, les isotopes stables de l’eau sont traditionnellement utilisés pour quantifier les changements passés de la température locale.Cette thèse doctorale a été initiée dans le cadre du programme de l’Agence Nationale de la Recherche ASUMA (“Improving the Accurancy of SUrface Mass balance of Antarctica”), ayant pour objectif de reconstruire et identifier les processus contrôlant la variabilité spatio-temporelle du bilan de masse de surface (BMS) de la Terre Adélie. J’ai utilisé des données d’isotopes stables de l’eau enregistrées dans des carottes de névé, des simulations atmosphériques produites par le modèle atmosphérique de circulation générale de haute résolution ECHAM5-wiso équipé des isotopes stables de l’eau, des réanalyses atmosphériques, des rétro-trajectoires, ainsi que des observations instrumentales satellitaires et de surface.Dans une première partie, j’ai évalué les capacités du modèle ECHAM5-wiso à simuler les températures de l’Antarctique, le BMS, le δ18O et le d-excess (ci-après, d-excess), comme prérequis à l’exploitation du modèle pour interpréter les compositions isotopiques. J’ai développé des diagnostics pour les relations δ18O-température et d-excess- δ18O sur l’ensemble du continent de l’Antarctique, en montrant que les différences issues des pentes des relations δ18O-température spatiales, inter-annuelles et saisonnières. Au sein du groupe de travail international de PAGES (Past Global Changes) Antarctica2k, j’ai utilisé des calibrations établies issues du modèle ECHAM5-wiso pour reconstruire la température de 7 régions d’Antarctique à partir d’une synthèse d’enregistrements de δ18O issus de carottes de glace couvrant les 2 000 dernières années.Dans une seconde partie, de nouveaux enregistrements issus de deux carottes de névé extraites en Terre Adélie, la S1C1 et la TA192A, ont été exploités, couvrant respectivement les périodes 1947-2007 et 1998-2014. Les BMS reconstruits décrivent une grande variabilité spatiale (74,11 ± 14,1 cm w.e. y-1 et 21,8 ± 6,9 cm w.e. y-1 pour la TA192A et la S1C1 respectivement), cohérente avec les données de balise disponibles. En utilisant une base de données mise à jour des isotopes stables de l’eau de l’Antarctique, j’ai montré que les valeurs moyennes isotopiques de Terre Adélie appartiennent à l’intervalle des valeurs côtières de l’Antarctique. Des analyses statistiques montrent une absence de relation entre nos enregistrements avec la température de surface locale à l’échelle inter-annuelle, mais des relations significatives avec des rétro-trajectoires atmosphériques et des simulations isotopiques suggérant que les isotopes de l’eau de la Terre Adélie fournissent des indications de la variabilité de la dynamique atmosphérique et du transport d’humidité, aux échelles saisonnière et inter-annuelle.Les analyses de cette thèse ont été limitées par la quantité d’enregistrements isotopiques disponibles pour la Terre Adélie, ainsi que par le manque de compréhension des effets de dépôt et de post-dépôt. Il est donc nécessaire d’exploiter les nouvelles carottes de névé extraites au cours du programme ASUMA, et d’effectuer en continu des mesures de la composition isotopique des précipitations, de la vapeur d’eau et de la neige de surface de Terre Adélie, en combinaison avec des outils de simulations atmosphériques, tels que des rétro-trajectoires associées à un diagnostic des sources d’humidité, et des modèles atmosphériques de circulation générale et régionaux équipés des isotopes stables de l’eau. / Documenting recent Antarctic climate variability is needed in order to understand the mechanisms at play, associated with the role of Antarctic mass balance for global sea level rise. Proxy records from Antarctic shallow firn cores are precious data, which complement instrumental and remote sensing observations to continuously cover the whole continent. Within these ice cores, water stable isotopes are commonly used to quantify past changes in local temperature.This PhD thesis was initiated within the French Agence Nationale de la Recherche “Improving the Accurancy of SUrface Mass balance of Antarctica” (ASUMA) project, which aims to reconstruct and to identify the processes controlling the spatio-temporal variability of the surface mass balance (SMB) in Adélie Land. I used water stable isotopes records from recently drilled shallow firn cores, as well as atmospheric simulations performed with the high resolution atmospheric general circulation model ECHAM5-wiso model, equipped with water stable isotopes, atmospheric reanalyses and back-trajectories, instrumental and remote sensing climate observations.In a first part, I assessed the skills of the ECHAM5-wiso with respect to Antarctic temperature, SMB, δ18O and deuterium excess (hereafter d-excess), as a prerequisite for the exploitation of the model to interpret isotope compositions. I developed Antarctic-wide diagnostics of the δ18O-temperature and d-excess- δ18O relationships, showing differences in the spatial, seasonal and interannual δ18O-temperature slopes. Within the international working group of PAGES (Past Global Changes) Antarctica 2k, I used the calibrations inferred from ECHAM5-wiso to reconstruct temperatures over 7 Antarctic regions from a synthesis of ice core δ18O records spanning the past 2,000 years.In a second part, new water stable isotope records from two firn core drilled in Adélie Land, the S1C1 and the TA192A, were investigated, covering the periods 1947-2007 and 1998-2014 respectively. The reconstructed SMB display a high spatial variability (74.1 ± 14.1 cm w.e. y-1 and 21.8 ± 6.9 cm w.e. y-1 for the TA192A and S1C1 respectively), consistent with Adélie Land stake data. Using an updated database of Antarctic water stable isotope datasets, I showed that the mean isotopic values (δ18O and d-excess) in Adélie Land are in line with the range of Antarctic coastal values. Statistical analyses show no relationship between our records and local surface air temperature, at the inter-annual scale, but significant relationships with atmospheric back-trajectories and isotopic simulations, suggesting that water stable isotopes in Adélie Land provide fingerprints of the variability of atmospheric dynamics and moisture transport, at the seasonal and inter-annual scales.The analyses performed during this PhD thesis have been limited by the few available Adélie Land water stable isotope records, and by the lack of understanding of deposition and post-deposition processes. Further work is thus needed to exploit the new firn cores drilled within the ASUMA project, and to monitor continuously Adélie Land water stable isotopes in precipitation, surface water vapour and surface snow, in combination with tools of atmospheric simulations such as back-trajectory simulations provided with moisture sources diagnostics, as well as water stable isotopes-enabled atmospheric general and regional circulation models.

Climatological winds over the equatorial Indian Ocean (EqlO) are westerly most of the year. Twice a year, in April-May ("spring") and October-December ("fall"), strong, sustained westerly winds generate eastward equatorial jets in the ocean. There are several unresolved issues related to the equatorial jets. They accelerate rapidly to speeds over lms"1 when westerly wind stress increases to about 0.7 dyne cm"2 in spring and fall, but decelerate while the wind stress continues to be westerly; each jet is followed by westward flow in the upper ocean lasting a month or longer. In addition to the semi-annual cycle, the equatorial winds and currents have strong in-traseasonal fluctuations. Observations show strong 30-60 day variability of zonal flow, and suggest that there might be variability with periods shorter than 20 days in the central EqlO. Observations from moored current meter arrays along 80.5°E south of Sri Lanka showed a distinct 15 day oscillation of equatorial meridional velocity (v) and off-equatorial zonal velocity (u). Recent observations from current meter moorings at the equator in the eastern EqlO show continuous 10-20 day, or biweekly, oscillations of v. The main motivation for the present study is to understand the dynamics of intraseasonal variability in the Indian Ocean that has been documented in the observational literature. What physical processes are responsible for the peculiar behavior of the equatorial jets? What are the relative roles of wind stress and large scale ocean dynamics? Does intraseasonal variability of wind stress force intraseasonal jets? What is the structure and origin of the biweekly variability? The intraseasonal and longer timescale variability of the equatorial Indian Ocean circulation is studied using an ocean general circulation model (OGCM) and recent in Abstract ii situ observations. The OGCM simulations are validated against other available observations. In this thesis, we document the space-time structure of the variability of equatorial Indian Ocean circulation, and attempt to find answers to some of the questions raised above. The main results are based on OGCM simulations forced by high frequency reanalysis and satellite scatterometer (QuikSCAT) winds. Several model experiments with idealized winds are used to interpret the results of the simulations. In addition to the OGCM simulations, the origin of observed intraseasonal anomalies of sea surface temperature (SST) in the eastern EqlO and Bay of Bengal, and related air-sea interaction, are investigated using validated satellite data. The main findings of the thesis can be summarized as: • High frequency accurate winds are required for accurate simulation of equatorial Indian Ocean currents, which have strong variability on intraseasonal to interannual time scales. • The variability in the equatorial waveguide is mainly driven by variability of the winds; there is some intraseasonal variability near the western boundary and in the equatorial waveguide due to dynamic instability of seasonal "mean" flows. • The fall equatorial jet is generally stronger and longer lived than the spring jet; the fall jet is modulated on intraseasonal time scales. Westerly wind bursts can drive strong intraseasonal equatorial jets in the eastern EqlO during the summer monsoon. • Eastward equatorial jets create a westward zonal pressure gradient force by raising sea level, and deepening the thermocline, in the east relative to the west. The zonal pressure force relaxes via Rossby wave radiation from the eastern boundary. • The zonal pressure force exerts strong control on the evolution of zonal flow; the decel eration of the eastward jets, and the subsequent westward flow in the upper ocean in the presence of westerly wind stress, is due to the zonal pressure force. • Neither westward currents in the upper ocean nor subsurface eastward flow (the ob served spring and summer "undercurrent") requires easterly winds; they can be gener ated by equatorial adjustment due to Kelvin (Rossby) waves generated at the western (eastern) boundary. • The biweekly variability in the EqlO is associated with forced mixed Rossby-gravity (MRG) waves generated by intraseasonal variability of winds. The biweekly MRG wave in has westward and upward phase propagation, zonal wavelength of 3000-4500 km and phase speed of 4 m s"1; it is associated with deep off equatorial upwelling/downwelling. Intraseasonal SST anomalies are forced mainly by net heat flux anomalies in the central and eastern EqlO; the large northward propagating SST anomalies in summer in the Bay of Bengal are due to net heat flux anomalies associated with the monsoon active-break cycle. Coherent variability in the atmosphere and ocean suggests air-sea interaction.

Implementing a nonlinear gravity wave (GW) parameterization into a mechanistic middle and upper atmosphere model, which extends to the lower thermosphere (160 km), we study the response of the atmosphere in terms of the circulation patterns, temperature distribution, and migrating terdiurnal solar tide activity to the upward propagating smallscale internal GWs originating in the lower atmosphere. We perform three test simulations for the Northern Hemisphere winter conditions in order to assess the effects of variations in the initial GWspectrum on the climatology and tidal patterns of the mesosphere and lower thermosphere. We find that the overall strength of the source level momentum flux has a relatively small impact on the zonal mean climatology. The tails of the GW source level spectrum, however, are crucial for the lower thermosphere climatology. With respect to the terdiurnal tide, we find a strong dependence of tidal amplitude on the induced GW drag, generally being larger when GW drag is increased.

I present herein an ensemble of ROCKE-3D aqua-planet simulations which I compare with the simulations presented in the work of Yang et al. (2019) and other similar works. The focus was on contrasting differences in the greenhouse effect between the models. In contrast to their work, I examined simulations with a dynamic ocean instead of a slab ocean, as well as the inclusion of dynamic sea ice for 2 out of 4 of them. A subset of the simulations examined prevented the formation of sea ice to make them more comparable to Yang et al. (2019), but they never reached radiative equilibrium and this made it difficult to utilize their results. When contrasting the sea ice simulations of ROCKE-3D with the CAM4_Wolf/ExoCAM simulation of Komacek & Abbot(2019), I found that the inclusion of ocean heat transport through a dynamic ocean increases the ice-free region around the sub-tropics for the rapidly rotating aqua-planet around a G-star, thus, resulting in a lower Bond albedo and more surface warming. Supporting previous intercomparisons (Sergeevet al. 2021), ROCKE-3D produces less low- to midlevel clouds toward the equator/substellar point, compared to other models. Consequently, this leads to less cooling through the shortwave cloud radiative forcing. Lastly, I looked at the specific humidity. ROCKE-3D produced the highest stratospheric water vapor content in the M-star scenario, which suggests that ROCKE-3D is closer to the moist greenhouse limit of Kasting et al. (1993); although, the model is still far off. / I detta arbete så presenterar jag vattenplanet simulationer producerade med hjälp av ROCKE-3D, som jag sedan jämför med simulationerna som presenteras i Yang et al. (2019). Fokuset för jämförelsen låg på att jämföra skillnader gällande den producerade växthuseffekten. Alla simulationer utnyttjade ett dynamiskt hav i stället för ett enklare "platt hav", varav 2 utav 4 av dom simulationer som presenteras tillåter havsis att formas. De simulationer som inte tillät is uppnådde aldrig termisk jämvikt vilket gjorde det svårt att antyda något utifrån dom. Vid jämförelse av is-simulationerna som producerades av ROCKE-3D med de is-simulationer producerade med hjälp av CAM4/ExoCAM i Komacek & Abbot (2019), så visades det sig att den is-fria regionen runtomkring de subtropiska områderna vart större för den snabbt roterande vattenplaneten runt en G-klassad stjärna om man inkluderar ett dynamiskt hav i stället för ett "platt hav". Vidare, så stödjer detta arbete dom resultat presenterade i Sergeev et al. (2021), där ROCKE-3D producerar mindre låg- och medelhöga molnformationer runtom ekvatorn/substellära regionen, jämfört med andra modeller. Vidare, så leder detta till en mindre kylningseffekt genom molnens reflektion av kortvågsstrålning. Sist, så undersökte jag den specifika fuktigheten, där ROCKE-3D visar på ett högre värde av stratosfärisk vattenånga i fallet av en tidsvattenlåst havsplanet runt en M-stjärna. Detta tyder på att ROCKE-3D är närmare den fuktiga växthusgränsen som presenteras i Kasting et al. (1993). Dock, så är den fortfarande långt ifrån att uppnås.

Implementing a nonlinear gravity wave (GW) parameterization into a mechanistic middle and upper atmosphere model, which extends to the lower thermosphere (160 km), we study the response of the atmosphere in terms of the circulation patterns, temperature distribution, and migrating terdiurnal solar tide activity to the upward propagating small scale internal GWs originating in the lower atmosphere. We perform three test simulations for the Northern Hemisphere winter conditions in order to assess the effects of variations in the initial GW spectrum on the climatology and tidal patterns of the mesosphere and lower thermosphere. We find that the overall strength of the source level momentum flux has a relatively small impact on the zonal mean climatology. The tails of the GW source level spectrum, however, are crucial for the lower thermosphere climatology. With respect to the terdiurnal tide, we find a strong dependence of tidal amplitude on the induced GW drag, generally being larger when GW drag is increased.

In first part of the study we discuss impact of dynamical core in simulation of tropical rainfall. Over years many new dynamical cores have been developed for atmospheric models to increase efficiency and reduce numerical errors. CAM3.1 gives an opportunity to study the impact of the dynamical core on simulations with its three dynamical cores namely Eulerian spectral(EUL) , Semilagrangian dynamics(SLD) and Finite volume(FV) coupled to a single parametrization package. A past study has compared dynamical cores of CAM3 in terms on tracer transport and has showed advantages using FV in terms of tracer transport. In this study we compare the dynamical cores in climate simulations and at their optimal configuration, which is the intended use of the model. The model is forced with AMIP type SST and rainfall over seasonal, interannual scales is compared. The significant differences in simulation of seasonal mean exist over tropics and over monsoon regions with observations and among dynamical cores. The differences among EUL and SLD, which use spectral transform methods are lesser compared that of with FV clearly indicating role of numerics in differences. There exist major errors in simulation of seasonal cycle in all dynamical cores and errors in simulation of seasonal means over many regions are associated with errors in simulation of seasonal cycle such as over south china sea. Seasonal cycle in FV is weaker compared to SLD and EUL. The dynamical cores exhibit different interannual variability of rainfall over Indian monsoon region, the period of maximum power corresponding to a dynamical core differs substantially with another. From this study there seems no superiority associated with FV dynamical core over all climate scales as seen in tracer transport. The next part of the study deals with impact of diurnal ocean atmosphere coupling in an AGCM,CAM3.1. Due to relatively low magnitude of diurnal cycle of SST and lack of SST observations over diurnal scales current atmospheric models are forced with SSTs of periods grater than a day. CAM 3.1 standalone model is forced with monthly SSTs but the interpolation is linear to every time step between any two months and this linear interpolation implies a linear diurnal and intraseasonal variation of SST which is not true in nature. To test the sensitivity of CAM3.1 to coupling of SST on diurnal scales, we prescribed over tropics(20S20N) a diurnal cycle of SST over daily mean interpolated SST of different magnitudes and phase comparable to observations. This idea of using a diurnal cycle of SST retaining seasonal mean SST in an atmospheric model is novel and provides an interesting frame work to test sensitivity of model to interpolations used in coupling of boundary conditions. Our analysis shows a high impact of using diurnal cycle of SST on simulation of mean rainfall over tropics. The impact in a case where diurnal cycle of SST is fixed and retained to daily mean SST implies that changes associated with a coupled model are to some extent due to change in representation of diurnal cycle of SST. A decrease of excess rainfall over western coast of Bay of Bengal and an increase of rainfall over northern bay of Bengal in such case is similar to the improvement due to coupling atmospheric model to a slab ocean model. This also implies that problems with current AMIP models in simulation of seasonal mean Indian monsoon rainfall could be due to erroneous representation of diurnal cycle of SST in models over this region where the diurnal cycle of SST is high in observations. The high spatial variability of the impact in various cases over tropics implies that a similar spatial variation of diurnal cycle could be important for accurate simulation of rainfall over tropics. Preliminary analysis shows that impact on rainfall was due to changes in moisture convergence. We also hypothesized that diurnal cycle of SST could trigger convection over regions such as northern Bay of Bengal and rainfall convergence feedback sustains it. The impact was also found on simulation of internal interannual variability of rainfall

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.

The most important component of the climate system over the Indian Ocean region is the southwest monsoon, which dictates the life and economy of billions of people in the tropics. Being a phenomena that involves interaction between atmosphere, ocean and land, the southwest monsoon is strongly influenced by upper ocean, primarily through warm sea surface temperature (SST). This is particularly true about the southeastern Arabian Sea (SEAS) and the onset of southwest monsoon over the peninsular India. A localized patch of warm water, known as the Arabian Sea mini warm pool (ASMWP), forms in the SEAS during February–March. It remain as the warmest spot in the northern Indian Ocean till early April. A large region, surrounding the SEAS, attains SST exceeding 30°C during April–May, with often the ASMWP as its core. The ASMWP is believed to have a critical impact on the air-sea interaction during the onset phase of southwest monsoon and on the formation of the onset vortex, during late May or early June. This thesis addresses the formation mechanisms of ASMWP, using a high-resolution Ocean General Circulation Model (OGCM) of the Indian Ocean. In addition to the formation of ASMWP, the SEAS is characterized by several features in its hydrography and circulation, which have been invoked in the past to explain the preferential warming of this oceanic region. During November–January, the prevailing surface currents transport low-salinity water from the Bay of Bengal into the SEAS and leads to strong haline stratification in the upper layer and formation of barrier layer (layer between mixed layer and isothermal layer). The vertical distribution of temperature in the SEAS exhibit inversions (higher subsurface temperature than that at surface) during December–February. A high in sea level and anticyclonic eddies develop in the SEAS during December and they propagate westward. These eddies modify the hydrography through downwelling and play an important role in the redistribution of advected low-salinity water within the SEAS. The seasonally reversing coastal and equatorial currents present in and around SEAS also have a major contribution in setting up the hydrography, through the advection and redistribution of cooler low-salinity water. These features make the SEAS a unique oceanographic region. The first hypothesis on the formation of ASMWP, which has been suggested by diagnostic studies, is based on the barrier layer mechanism. The barrier layer, caused by the influx of low-salinity water at surface, is argued to maintain a shallow mixed layer which can warm more efficiently. In addition, presence of barrier layer can prevent mixed layer cooling, by cutting off the interaction of mixed layer with cooler thermocline water below. However, a coupled model study have shown that there is no significant impact on the ASMWP formation from barrier layer, but only a weak warming effect during it mature phase during April. The second hypothesis, which is based on an OGCM study, has suggested that the temperature inversions present within the barrier layer can heat the mixed layer through turbulent entrainment and in turn lead to the formation of ASMWP during February–March. Both hypotheses rule out the possibility of air-sea heat fluxes being the primary reason in its formation. The strong salinity stratification in the SEAS during December–March is central to the hypotheses about formation of the ASMWP. Observational studies have only limited success in assessing the contribution from barrier layer and temperature inversions, as the ASMWP always form in their presence. OGCMs offer a better alternative. However, modelling processes in the northern Indian Ocean, especially that in the SEAS, is a challenging problem. Previous Indian Ocean models have had serious difficulties in simulating the low-salinity water in the Bay of Bengal and its intrusion into the SEAS. The northward advection of low-salinity water in the SEAS, along the west coast of India, is used to be absent in model simulations. Moreover, the coarse resolution inhibited those models from simulating faster surface currents and vigorous eddies as seen in the observations. In this thesis, we use an OGCM of the Indian Ocean, based on the recent version of Modular Ocean Model (MOM4p0), to study the ASMWP. The model has high resolutions in the horizontal (1/4o x 1/4o) and vertical (40 levels, with 5 m spacing in upper 60 m), and has been forced with daily values momentum, heat and freshwater fluxes. The turbulent (latent and sensible) and long wave heat fluxes have been calculated as a function of model SST. The freshwater forcing consists of precipitation, evaporation and river runoff, and there are no surface restoring or flux adjustments. The river runoff has been distributed over several grid points about the river mouth instead of discharging into a singe grid point, which has resulted in remarkable improvements in salinity simulation. The model simulates the Indian Ocean temperature, salinity and circulation remarkably well. The pattern of model temperature distribution and evolution matches very well with that in the observations. Significant improvements have been made in the salinity simulation, including the Bay of Bengal freshwater plume and intrusion of low-salinity water from the bay into the SEAS. The salinity distribution within the SEAS is also well represented in the model. The use of appropriate horizontal friction parameters has resulted in the simulation of realistic currents. The observed features in the SEAS, including the life cycle of the ASMWP, low-salinity water, barrier layer, temperature inversions, eddies and currents are well represented in the model. Present study has unraveled the processes involved in the life cycle of barrier layer and temperature inversions in the SEAS. Presence of low-salinity water is necessary for their formation. Barrier layer develops in the SEAS during November, after the intrusion of low-salinity water from the Bay of Bengal. The barrier layer is thickest during January–February, and it dissipates during March–April. The variations and peak of barrier layer thickness is controlled by variations in isothermal layer depth, which in turn is dominated by the downwelling effects of anticyclonic eddies. The intense solar heating during March–April leads to the formation of shallow isothermal layer and results in the dissipation of barrier layer. Temperature inversions starts developing in the SEAS during December, reaches its peak during January–February and dissipates in the following months. Advection of cooler low-salinity water over warmer salty water and penetrating shortwave radiation is found to cause temperature inversions within the SEAS, whereas winter cooling is also important to the north and south of the SEAS. There is significant variation in the magnitude, depth of occurrence and formation mechanisms of temperature inversions within the SEAS. Analysis of model mixed layer heat budget has shown that the SEAS SST is mainly controlled by atmospheric forcing, including the life cycle of ASMWP. It has also shown that the heating from temperature inversions do not contribute to the formation of ASMWP. In an experiment in which a constant salinity of 35 psu was maintained over the entire model domain, the ASMWP evolved very similar to that in the standard run, suggesting that the salinity effects are not necessary for the formation of ASMWP. Examination of wind field show that the winds over the SEAS during November–February are low due to the blocking of northeasterly winds by Western Ghats. Several process experiments by modifying the wind and turbulent heat fluxforcing fields have shown that these low winds lead to the formation of ASMWP in the SEAS during February–March. The low winds reduce latent heat loss, resulting in net heat gain by the ocean. This helps the SEAS to keep warmer SST while the surrounding region experience intense cooling under the strong dry northeasterly winds. As the winds are weak over the SEAS, the mixed layer is not able to feel the stratification beneath and the mixed layer depth is determined by solar heating, with or without salinity effects. In addition, the weak winds are not able to entrain the temperature inversions present in the barrier layer. The winds are weak during March–April too, and the air-sea heat fluxes dictate the SST evolution during this period. Therefore, during November–April, the SEAS acts as a low wind heat-dominated regime, where the evolution of sea surface temperature is solely determined by atmospheric forcing. We show that, in such regions, the evolution of surface layer temperature is not dependent on the characteristics of subsurface ocean, including the presence of barrier layer and temperature inversions.