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.

Sur les effets du changement climatique aux échelles globale et régionale. Il montre en particulierqu’aucun consensus ne peut être trouvé pour ce qui concerne l’évolution future de lapluviométrie — et de la dynamique atmosphérique associée — en région de mousson africaine.Ce mémoire revisite cette question à la lumière des nouvelles données disponibles et selon uneapproche évitant toute surreprésentation du nombre de simulations disponibles pour un type demodèle donné, tout en prenant en compte la diversité des modèles ainsi que leur évolution dansle temps : sorties de vingt modèles de circulation générale (MCGs) ayant participé aux exercicesCMIP3 (douze MCGs) et CMIP5 (huit MCGs) sous les scénarios d’émissions A1B et rcp4.5,respectivement. Les sorties sont analysées principalement sur deux fenêtres de quarante ans —périodes actuelle (1960-1999) et future (2031-2070) — et les résultats discutés au regard de leurvraisemblance selon une approche permettant à la fois de quantifier les différences futur moinsactuel, de mesurer les significativités et les robustesses statistiques et d’associer une probabilitémesurant le consensus des modèles en fonction des échelles et des variables considérées.Les analyses menées sur CMIP3 et CMIP5 montrent qu’un consensus sur l’effet du changementclimatique en Afrique de l’Ouest peut être obtenu si l’on ne fait pas de l’ensemble de labande sahélienne une entité homogène et qu’on raisonne à des échelles spatiales inférieures. Lesrésultats révèlent une évolution contrastée entre le centre et l’ouest du Sahel avec, pour le futur(i) une hausse des précipitations au centre s’expliquant surtout par une plus grande convergencedes flux dans les basses couches, ainsi qu’une pénétration plus au nord de la mousson ;(ii) une baisse des précipitations à l’ouest s’expliquant par le renforcement de la circulation detype Walker, du Jet d’Est Africain (JEA) et de la subsidence dans les couches moyennes. Parailleurs, on peut s’attendre à une modification du cycle annuel moyen avec un retrait retardé dela mousson. Ce retard est notamment lié aux apports supplémentaires d’humidité depuis l’Atlantique,dus au renforcement des contrastes thermiques et d’humidité entre océan et continent,mais aussi et surtout aux apports tardifs d’humidité depuis la Méditerranée et au renforcementdes flux de nord en septembre et octobre en direction du Sahel / The fourth IPCC report in 2007 established the synthesis of previously published work onthe effects of climate change on global and regional scales. It shows in particular that no consensuscan be found with regard to the future of rainfall — and atmospheric dynamics- associatedwith region — African monsoon. This dissertation revisits this issue in the light of new dataand using an approach avoiding over-representation of the number of simulations available forone type of model and taking into account the diversity of models and their evolution in time :twenty general circulation models (GCMs) participating in the exercises CMIP3 (twelve GCMs)and CMIP5 (eight GCMs) under the A1B emissions scenario and rcp4.5, respectively. Outputsare analyzed on two 40-year periods, — ‘Present’ (1960-1999) and ‘Future’ (2031-2070) — anddiscussed in terms of likelihood, through an approach allowing us to both quantify differences‘future’ minus ‘present’, measure robustness and statistical significances and associate a probabilitymeasuring the model consensus as a function of scales and variables.Analyzes conducted on CMIP3 and CMIP5 show that consensus on the effect of climatechange in West Africa can be achieved if we do not consider the Sahel as a whole and homogeneousentity but at lower scales. The results show contrasted responses over the centraland western Sahel, with for the future, (i) an increase in precipitation in the central regionexplained primarily by a greater convergence of flow in the lower layers and a most northerlymonsoon penetration over the continent, (ii) a rainfall decrease in the western Sahel explainedby increased Walker-type circulation, African easterly jet and mid-level subsidence. Moreover,we can expect a change in the mean annual cycle of the monsoon season with a delayed withdrawallinked to additional inputs of moisture from the Atlantic due to increasing thermal andmoisture contrasts between ocean and continent but also to a stronger contribution of moisturefluxes in September and October from the Mediterranean into the Sahel

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

General Circulation Models (GCMs), which are mathematical models based on principles of fluid dynamics, thermodynamics and radiative transfer, are the most reliable tools available for projecting climate change. However, the spatial scale on which typical GCMs operate is very coarse as compared to that of a hydrologic process and hence, the output from a GCM cannot be directly used in hydrologic models. Statistical Downscaling (SD) derives a statistical or empirical relationship between the variables simulated by the GCM (predictors) and a point-scale meteorological series (predictand). In this work, a new downscaling model called CRF-downscaling model, is developed where the conditional distribution of the hydrologic predictand sequence, given atmospheric predictor variables, is represented as a conditional random field (CRF) to downscale the predictand in a probabilistic framework. Features defined in the downscaling model capture information about various factors influencing precipitation such as circulation patterns, temperature and pressure gradients and specific humidity levels. Uncertainty in prediction is addressed by projecting future cumulative distribution functions (CDFs) for a number of most likely precipitation sequences. Direct classification of dry/wet days as well as precipitation amount is achieved within a single modeling framework, and changes in the non-parametric distribution of precipitation and dry and wet spell lengths are projected. Application of the method is demonstrated with the case study of downscaling to daily precipitation in the Mahanadi basin in Orissa, with the A1B scenario of the MIROC3.2 GCM from the Center for Climate System Research (CCSR), Japan. An uncertainty modeling framework is presented in this work, which combines GCM, scenario and downscaling uncertainty using the Dempster-Shafer (D-S) evidence theory for representing and combining uncertainty. The methodology for combining uncertainties is applied to projections of hydrologic drought in terms of monsoon standardized streamflow index (SSFI-4) from streamflow projections for the Mahanadi river at Hirakud. The results from the work indicate an increasing probability of extreme, severe and moderate drought and decreasing probability of normal to wet conditions, as a result of a decrease in monsoon streamflow in the Mahanadi river due to climate change. In most studies to date, the nature of the downscaling relationship is assumed stationary, or remaining unchanged in a future climate. In this work, an uncertainty modeling framework is presented in which, in addition to GCM and scenario uncertainty, uncertainty in the downscaling relationship itself is explored by linking downscaling with changes in frequencies of modes of natural variability. Downscaling relationships are derived for each natural variability cluster and used for projections of hydrologic drought. Each projection is weighted with the future projected frequency of occurrence of that cluster, called ‘cluster-linking’, and scaled by the GCM performance with respect to the associated cluster for the present period, called ‘frequency scaling’. The uncertainty modeling framework is applied to a case study of projections of hydrologic drought or SSFI-4 classifications, using projected streamflows for the Mahanadi river at Hirakud. It is shown that a stationary downscaling relationship will either over- or under-predict downscaled hydrologic variable values and associated uncertainty. Results from the work show improved agreement between GCM predictions at the regional scale, which are validated for the 20th century, implying that frequency scaling and cluster-linking may indeed be a valid method for constraining uncertainty. To assess the impact of climate change on reservoir performance, in this study, a range of integrated hydrologic scenarios are projected for the future. The hydrologic scenarios incorporate increased irrigation demands; rule curves dictated by increased need for flood storage and downscaled projections of streamflow from an ensemble of GCMs and emission scenarios. The impact of climate change on multipurpose reservoir performance is quantified, using annual hydropower and RRV criteria, under GCM and scenario uncertainty. The ‘business-as-usual’ case using Standard Operating Policy (SOP) is studied initially for quantifying impacts. Adaptive Stochastic Dynamic Programming (SDP) policies are subsequently derived for the range of future hydrologic scenarios, with the objective of maximizing reliabilities with respect to multiple reservoir purposes of hydropower, irrigation and flood control. It is shown that the hydrologic impact of climate change is likely to result in decreases in performance criteria and annual hydropower generation for Hirakud reservoir. Adaptive policies show that a marginal reduction in irrigation and flood control reliability can achieve increased hydropower reliability in future. Hence, reservoir rules for flood control may have to be revised in the future.

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.

General Circulation Models (GCMs) are tools designed to simulate time series of climate variables globally, accounting for effects of greenhouse gases in the atmosphere. They attempt to represent the physical processes in the atmosphere, ocean, cryosphere and land surface. They are currently the most credible tools available for simulating the response of the global climate system to increasing greenhouse gas concentrations, and to provide estimates of climate variables (e.g. air temperature, precipitation, wind speed, pressure etc.) on a global scale. GCMs demonstrate a significant skill at the continental and hemispheric spatial scales and incorporate a large proportion of the complexity of the global system; they are, however, inherently unable to represent local subgrid-scale features and dynamics. The spatial scale on which a GCM can operate (e.g., 3.75° longitude x 3.75° latitude for Coupled Global Climate Model, CGCM2) is very coarse compared to that of a hydrologic process (e.g., precipitation in a region, streamflow in a river etc.) of interest in the climate change impact assessment studies. Moreover, accuracy of GCMs, in general, decreases from climate related variables, such as wind, temperature, humidity and air pressure to hydrologic variables such as precipitation, evapotranspiration, runoff and soil moisture, which are also simulated by GCMs. These limitations of the GCMs restrict the direct use of their output in hydrology. This thesis deals with developing statistical downscaling models to assess climate change impacts and methodologies to address GCM and scenario uncertainties in assessing climate change impacts on hydrology. Downscaling, in the context of hydrology, is a method to project the hydrologic variables (e.g., rainfall and streamflow) at a smaller scale based on large scale climatological variables (e.g., mean sea level pressure) simulated by a GCM. A statistical downscaling model is first developed in the thesis to predict the rainfall over Orissa meteorological subdivision from GCM output of large scale Mean Sea Level Pressure (MSLP). Gridded monthly MSLP data for the period 1948 to 2002, are obtained from the National Center for Environmental Prediction/ National Center for Atmospheric Research (NCEP/NCAR) reanalysis project for a region spanning 150 N -250 N in latitude and 800 E -900 E in longitude that encapsulates the study region. The downscaling model comprises of Principal Component Analysis (PCA), Fuzzy Clustering and Linear Regression. PCA is carried out to reduce the dimensionality of the larger scale MSLP and also to convert the correlated variables to uncorrelated variables. Fuzzy clustering is performed to derive the membership of the principal components in each of the clusters and the memberships obtained are used in regression to statistically relate MSLP and rainfall. The statistical relationship thus obtained is used to predict the rainfall from GCM output. The rainfall predicted with the GCM developed by CCSR/NIES with B2 scenario presents a decreasing trend for non-monsoon period, for the case study. Climate change impact assessment models developed based on downscaled GCM output are subjected to a range of uncertainties due to both ‘incomplete knowledge’ and ‘unknowable future scenario’ (New and Hulme, 2000). ‘Incomplete knowledge’ mainly arises from inadequate information and understanding about the underlying geophysical process of global change, leading to limitations in the accuracy of GCMs. This is also termed as GCM uncertainty. Uncertainty due to ‘unknowable future scenario’ is associated with the unpredictability in the forecast of socio-economic and human behavior resulting in future Green House Gas (GHG) emission scenarios, and can also be termed as scenario uncertainty. Downscaled outputs of a single GCM with a single climate change scenario represent a single trajectory among a number of realizations derived using various GCMs and scenarios. Such a single trajectory alone can not represent a future hydrologic scenario, and will not be useful in assessing hydrologic impacts due to climate change. Nonparametric methods are developed in the thesis to model GCM and scenario uncertainty for prediction of drought scenario with Orissa meteorological subdivision as a case study. Using the downscaling technique described in the previous paragraph, future rainfall scenarios are obtained for all available GCMs and scenarios. After correcting for bias, equiprobability transformation is used to convert the precipitation into Standardized Precipitation Index-12 (SPI-12), an annual drought indicator, based on which a drought may be classified as a severe drought, mild drought etc. Disagreements are observed between different predictions of SPI-12, resulting from different GCMs and scenarios. Assuming SPI-12 to be a random variable at every time step, nonparametric methods based on kernel density estimation and orthonormal series are used to determine the nonparametric probability density function (pdf) of SPI-12. Probabilities for different categories of drought are computed from the estimated pdf. It is observed that there is an increasing trend in the probability of extreme drought and a decreasing trend in the probability of near normal conditions, in the Orissa meteorological subdivision. The single valued Cumulative Distribution Functions (CDFs) obtained from nonparametric methods suffer from limitations due to the following: (a) simulations for all scenarios are not available for all the GCMs, thus leading to a possibility that incorporation of these missing climate experiments may result in a different CDF, (b) the method may simply overfit to a multimodal distribution from a relatively small sample of GCMs with a limited number of scenarios, and (c) the set of all scenarios may not fully compose the universal sample space, and thus, the precise single valued probability distribution may not be representative enough for applications. To overcome these limitations, an interval regression is performed to fit an imprecise normal distribution to the SPI-12 to provide a band of CDFs instead of a single valued CDF. Such a band of CDFs represents the incomplete nature of knowledge, thus reflecting the extent of what is ignored in the climate change impact assessment. From imprecise CDFs, the imprecise probabilities of different categories of drought are computed. These results also show an increasing trend of the bounds of the probability of extreme drought and decreasing trend of the bounds of the probability of near normal conditions, in the Orissa meteorological subdivision. Water resources planning requires the information about future streamflow scenarios in a river basin to combat hydrologic extremes resulting from climate change. It is therefore necessary to downscale GCM projections for streamflow prediction at river basin scales. A statistical downscaling model based on PCA, fuzzy clustering and Relevance Vector Machine (RVM) is developed to predict the monsoon streamflow of Mahanadi river at Hirakud reservoir, from GCM projections of large scale climatological data. Surface air temperature at 2m, Mean Sea Level Pressure (MSLP), geopotential height at a pressure level of 500 hecto Pascal (hPa) and surface specific humidity are considered as the predictors for modeling Mahanadi streamflow in monsoon season. PCA is used to reduce the dimensionality of the predictor dataset and also to convert the correlated variables to uncorrelated variables. Fuzzy clustering is carried out to derive the membership of the principal components in each of the clusters and the memberships thus obtained are used in RVM regression model. RVM involves fewer number of relevant vectors and the chance of overfitting is less than that of Support Vector Machine (SVM). Different kernel functions are used for comparison purpose and it is concluded that heavy tailed Radial Basis Function (RBF) performs best for streamflow prediction with GCM output for the case considered. The GCM CCSR/NIES with B2 scenario projects a decreasing trend in future monsoon streamflow of Mahanadi which is likely to be due to high surface warming. A possibilistic approach is developed next, for modeling GCM and scenario uncertainty in projection of monsoon streamflow of Mahanadi river. Three GCMs, Center for Climate System Research/ National Institute for Environmental Studies (CCSR/NIES), Hadley Climate Model 3 (HadCM3) and Coupled Global Climate Model 2 (CGCM2) with two scenarios A2 and B2 are used for the purpose. Possibilities are assigned to GCMs and scenarios based on their system performance measure in predicting the streamflow during years 1991-2005, when signals of climate forcing are visible. The possibilities are used as weights for deriving the possibilistic mean CDF for the three standard time slices, 2020s, 2050s and 2080s. It is observed that the value of streamflow at which the possibilistic mean CDF reaches the value of 1 reduces with time, which shows reduction in probability of occurrence of extreme high flow events in future and therefore there is likely to be a decreasing trend in the monthly peak flow. One possible reason for such a decreasing trend may be the significant increase in temperature due to climate warming. Simultaneous occurrence of reduction in Mahandai streamflow and increase in extreme drought in Orissa meteorological subdivision is likely to pose a challenge for water resources engineers in meeting water demands in future.

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.

Understanding climate change and variability at urban scale is essential for water resource management, land use planning, and development of adaption plans. However, there are serious challenges to meet these goals due to unavailability of observed and / or simulated high resolution spatial and temporal climate data. Recent efforts made possible the availability of high resolution climate data from non-hydrostatic regional climate model (RCM) and statistically downscaled General Circulation Models (GCMs). This study investigates trends in climate and urbanization and their impact on surface water supply for the city of Addis Ababa, Ethiopia. The methodology presented in this study focused on the observed and projected NIMRHadGEM2- AO model and Special Report on Emissions Scenarios (SRES) of B2 and A2 of HadCM3 model are also employed for rainfall, maximum temperature and minimum temperature data using for climate analysis. Water Evaluation and Planning (WEAP) modeling system was used for determination of climate and urbanization impacts on water. Land-Sat images were analyzed using Normalized Differencing Vegetation Index (NDVI). Statistical downscaling model (SDSM) was employed to investigate the major changes and intensity of the urban heat island (UHI). The result indicates monthly rainfall anomalies with respect to the baseline mean showing wet anomaly in summer (kiremt) during 2030s and 2050s, and a dry anomaly in the 2080s under A2 and B2 scenarios with exception of a wet anomaly in September over the city. The maximum temperature anomalies under Representative Concentration Pathways (RCPs) also show warming during near, mid and end terms. The mean monthly minimum temperature anomalies under A2 and B2 scenarios are warm but the anomalies are much lower than RCPs. The climate under the RCP 8.5 and high population growth (3.3 %) scenario will lead to the unmet demand of 462.77 million m3 by 2039. Future projection of urban heat island under emission pathway of A2 and B2 scenario shows that, the nocturnal UHI will be intense in winter or dry season episodes in the city. Under A2 scenario the highest urban warming will occur during October to December (2.5 ºC to 3.2 ºC). Under RCP 8.5 scenario the highest urban warming will occur during October to December (0.5 ºC to 1.0 °C) in the 2050s and 2080s. Future management and adaptation strategies are to expand water supply to meet future demand and to implement demand side water management systems of the city and UHI / Environmental Sciences / Ph. D. (Environmental Management)

Submitted by Marco Antonio Oliveira da Silva (maosilva@ipen.br) on 2016-08-26T12:15:03Z No. of bitstreams: 0 / Made available in DSpace on 2016-08-26T12:15:03Z (GMT). No. of bitstreams: 0 / O escoamento bifásico de gás-líquido é encontrado em muitos circuitos fechados que utilizam circulação natural para fins de resfriamento. O fenômeno da circulação natural é importante nos recentes projetos de centrais nucleares para a remoção de calor. O circuito de circulação natural (Circuito de Circulação Natural - CCN), instalado no Instituto de Pesquisas Energéticas e Nucleares, IPEN / CNEN, é um circuito experimento concebido para fornecer dados termo-hidráulicos relacionados com escoamento monofásico ou bifásico em condições de circulação natural. A estimativa de transferência de calor tem sido melhorada com base em modelos que requerem uma previsão precisa de transições de padrão de escoamento. Este trabalho apresenta testes experimentais desenvolvidos no CCN para a visualização dos fenômenos de instabilidade em ciclos de circulação natural básica e classificar os padrões de escoamento bifásico associados aos transientes e instabilidades estáticas de escoamento. As imagens são comparadas e agrupadas utilizando mapas auto-organizáveis de Kohonen (SOM), aplicados em diferentes características da imagem digital. Coeficientes da Transformada Discreta de Cossenos de Quadro Completo (FFDCT) foram utilizados como entrada para a tarefa de classificação, levando a bons resultados. Os protótipos de FFDCT obtidos podem ser associados a cada padrão de escoamento possibilitando uma melhor compreensão da instabilidade observada. Uma metodologia sistemática foi utilizada para verificar a robustez do método. / Dissertação (Mestrado em Tecnologia Nuclear) / IPEN/D / Instituto de Pesquisas Energéticas e Nucleares - IPEN-CNEN/SP

Submitted by Cássia Santos (cassia.bcufg@gmail.com) on 2016-05-31T09:35:51Z No. of bitstreams: 2 Tese - Frederico Augusto Martins Valtuille Faleiro - 2016.pdf: 7096330 bytes, checksum: 04cfce04ef128c5bd6e99ce18bb7f650 (MD5) license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) / Approved for entry into archive by Luciana Ferreira (lucgeral@gmail.com) on 2016-05-31T10:52:51Z (GMT) No. of bitstreams: 2 Tese - Frederico Augusto Martins Valtuille Faleiro - 2016.pdf: 7096330 bytes, checksum: 04cfce04ef128c5bd6e99ce18bb7f650 (MD5) license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) / Made available in DSpace on 2016-05-31T10:52:51Z (GMT). No. of bitstreams: 2 Tese - Frederico Augusto Martins Valtuille Faleiro - 2016.pdf: 7096330 bytes, checksum: 04cfce04ef128c5bd6e99ce18bb7f650 (MD5) license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) Previous issue date: 2016-03-07 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - CAPES / The climate changes are one of the major threats to the biodiversity and it is expected to increase its impact along the 21st century. The climate change affect all levels of the biodiversity from individuals to biomes, reducing the ecosystem services. Despite of this, the prediction of climate change impacts on biodiversity is still a challenge. Overcoming these issues depends on improvements in different aspects of science that support predictions of climate change impact on biodiversity. The common practice to predict the climate change impact consists in formulate ecological niche models based in the current climate and project the changes based in the future climate predicted by the climate models. However, there are some recognized limitations both in the formulation of the ecological niche model and in the use of predictions from the climate models that need to be analyzed. Here, in the first chapter we review the science behind the climate models in order to reduce the knowledge gap between the scientific community that formulate the climate models and the community that use the predictions of these models. We showed that there is not consensus about evaluate the climate models, obtain regional models with higher spatial resolution and define consensual models. However, we gave some guidelines for use the predictions of the climate models. In the second chapter, we tested if the predictions of correlative ecological niche models fitted with presence-absence match the predictions of models fitted with abundance data on the metrics of climate change impact on orchid bees in the Atlantic Forest. We found that the presence-absence models were a partial proxy of change in abundance when the output of the models was continuous, but the same was not true when the predictions were converted to binary. The orchid bees in general will decrease the abundance in the future, but will retain a good amount of suitable sites in the future and the distance to gained climatic suitable areas can be very close, despite of great variation. The change in the species richness and turnover will be mainly in the western and some regions of southern of the Atlantic Forest. In the third chapter, we discussed the drawbacks in using the estimations of realized niche instead the fundamental niche, such as overpredicting the effect of climate change on species’ extinction risk. We proposed a framework based on phylogenetic comparative and missing data methods to predict the dimensions of the fundamental niche of species with missing data. Moreover, we explore sources of uncertainty in predictions of fundamental niche and highlight future directions to overcome current limitations of phylogenetic comparative and missing data methods to improve predictions. We conclude that it is possible to make better use of the current knowledge about species’ fundamental niche with phylogenetic information and auxiliary traits to predict the fundamental niche of poorly-studied species. In the fourth chapter, we used the framework of the chapter three to test the performance of two recent phylogenetic modeling methods to predict the thermal niche of mammals. We showed that PhyloPars had better performance than Phylogenetic Eigenvector Maps in predict the thermal niche. Moreover, the error and bias had similar phylogenetic pattern for both margins of the thermal niche while they had differences in the geographic pattern. The variance in the performance was explained by taxonomic differences and not by methodological aspects. Finally, our models better predicted the upper margin than the lower margin of the thermal niche. This is a good news for predicting the effect of climate change on species without physiological data. We hope our finds can be used to improve the predictions of climate change effect on the biodiversity in future studies and support the political decisions on minimizing the effects of climate change on biodiversity. / As mudanças climáticas são uma das principais ameaças à biodiversidade e é esperado que aumente seu impacto ao longo do século XXI. As mudanças climáticas afetam todos os níveis de biodiversidade, de indivíduos à biomas, reduzindo os serviços ecossistêmicos. Apesar disso, as predições dos impactos das mudanças climáticas na biodiversidade é ainda um desafio. A superação dessas questões depende de melhorias em diferentes aspectos da ciência que dá suporte para predizer o impacto das mudanças climáticas na biodiversidade. A prática comum para predizer o impacto das mudanças climáticas consiste em formular modelos de nicho ecológico baseado no clima atual e projetar as mudanças baseadas no clima futuro predito pelos modelos climáticos. No entanto, existem algumas limitações reconhecidas na formulação do modelo de nicho ecológico e no uso das predições dos modelos climáticos que precisam ser analisadas. Aqui, no primeiro capítulo nós revisamos a ciência por detrás dos modelos climáticos com o intuito de reduzir a lacuna de conhecimentos entre a comunidade científica que formula os modelos climáticos e a comunidade que usa as predições dos modelos. Nós mostramos que não existe consenso sobre avaliar os modelos climáticos, obter modelos regionais com maior resolução espacial e definir modelos consensuais. No entanto, nós damos algumas orientações para usar as predições dos modelos climáticos. No segundo capítulo, nós testamos se as predições dos modelos correlativos de nicho ecológicos ajustados com presença-ausência são congruentes com aqueles ajustados com dados de abundância nas medidas de impacto das mudanças climáticas em abelhas de orquídeas da Mata Atlântica. Nós encontramos que os modelos com presença-ausência foram substitutos parciais das mudanças na abundância quando o resultado dos modelos foi contínuo (adequabilidade), mas o mesmo não ocorreu quando as predições foram convertidas para binárias. As espécies de abelhas, de modo geral, irão diminuir em abundância no futuro, mas reterão uma boa quantidade de locais adequados no futuro e a distância para áreas climáticas adequadas ganhadas podem estar bem próximo, apesar da grande variação. A mudança na riqueza e na substituição de espécies ocorrerá principalmente no Oeste e algumas regiões no sul da Mata Atlântica. No terceiro capítulo, nós discutimos as desvantagens no uso de estimativas do nicho realizado ao invés do nicho fundamental, como superestimar o efeito das mudanças climáticas no risco de extinção das espécies. Nós propomos um esquema geral baseado em métodos filogenéticos comparativos e métodos de dados faltantes para predizer as dimensões do nicho fundamental das espécies com dados faltantes. Além disso, nós exploramos as fontes de incerteza nas predições do nicho fundamental e destacamos direções futuras para superar as limitações atuais dos métodos comparativos filogenéticas e métodos de dados faltantes para melhorar as predições. Nós concluímos que é possível fazer melhor uso do conhecimento atual sobre o nicho fundamental das espécies com informação filogenética e caracteres auxiliares para predizer o nicho fundamental de espécies pouco estudadas. No quarto capítulo, nós usamos o esquema geral do capítulo três para testar a performance de dois novos métodos de modelagem filogenética para predizer o nicho térmico dos mamíferos. Nós mostramos que o “PhyloPars” teve uma melhor performance que o “Phylogenetic Eigenvector Maps” em predizer o nicho térmico. Além disso, o erro e o viés tiveram um padrão filogenético similar para ambas as margens do nicho térmico, enquanto eles apresentaram diferentes padrões espaciais. A variância na performance foi explicada pelas diferenças taxonômicas e não pelas diferenças em aspectos metodológicos. Finalmente, nossos modelos melhor predizem a margem superior do que a margem inferior do nicho térmico. Essa é uma boa notícia para predizer o efeito das mudanças climáticas em espécies sem dados fisiológicos. Nós esperamos que nossos resultados possam ser usados para melhorar as predições do efeito das mudanças climáticas na biodiversidade em estudos futuros e dar suporte para decisões políticas para minimização dos efeitos das mudanças climáticas na biodiversidade.