Research leading to this thesis aims to assess the policy relevant mitigation potential of Indian forests as well as aims to assess the impact of climate change on carbon stocks, vegetation boundary shifts, Net Primary Productivity (NPP) and the mitigation potential of Indian forests. To project the impact of climate change we chose a dynamic global vegetation model ‘Integrated Biosphere Simulator’ (IBIS V.2.6b4). We selected A2 and B2 scenarios for projecting the impacts. Mitigation potential was assessed using the ‘Generalized Comprehensive Mitigation Assessment Process’ (GCOMAP) model. We assess the mitigation potential of Indian forests in the light of India’s long-term policy objective of bringing 33% of its total geographical area under forest cover. We analyzed the mitigation potential of this policy objective under two scenarios: the first comprising of rapid afforestation scenario with the target to achieving the goal by 2020 and the second a moderate afforestation scenario in which this goal is achieved by 2030. We estimate that afforestation could offset about 9% of India’s average national emissions over the 2010-2030 period, while about 6.7% could be mitigated under the moderate afforestation scenario over the same period. We analyze the impact of climate change on the four key attributes of Indian forests, i.e. impact on vegetation distribution, impact on forest productivity (NPP), impact on soil carbon (SOC) and impact on biomass carbon. IBIS simulations suggest that approximately 39% and 34% of forest grids are projected to experience change in vegetation type under A2 and B2 climate scenarios, respectively over the period 2070¬2100. Simulations further indicate that NPP is projected to increase by an average of 66% under the A2 scenario and 49% under the B2 scenario. The increase is higher in the northeastern part of India. However, in the central and western Indian forests NPP remains stable or increases only moderately, and in some places even decreases. Our assessment of the impact of climate change on Soil Organic Carbon (SOC) suggests a trend similar to NPP distribution, which is to be expected as increased NPP is the primary driver of higher litter input to the soil. However, the quantum of increase in this case is lower, around 37% and 30%, for the A2 and B2 scenario respectively (averaged over India). The biomass carbon is also projected to increase all over India on the lines similar to NPP gains. However, projected gains in biomass, NPP and SOC should be viewed with caution as IBIS tends to simulate a fairly strong CO2 fertilization effect that may not necessarily be realized under conditions of nutrient and water limitations and under conditions of increased pest and fire outbreaks. Further we analyzed the impact of climate change on the mitigation potential of Indian forests by linking impact assessment models to mitigation potential assessment model GCOMAP. Two impact assessment models BIOME4 and IBIS are used for simulating the impact of climate change. IBIS is a dynamic vegetation model while BIOME4 is an equilibrium model. Our assessment suggests that with the BIOME4 simulations the cumulative mitigation potential increases by up to 21% under the A2 scenario over the period 2008 to 2108, whereas, under the B2 scenario the mitigation potential increases only by 14% over the same period. However cumulative mitigation potential estimates obtained from the IBIS simulations suggest much smaller gains, where mitigation potential increases by only 6% and 5% over the period 2008 to 2108, under A2 and B2 scenarios respectively. To enable effective policy analysis and to build a synergy between the mitigation and adaptation efforts in the Indian forest sector, a vulnerability index for the forested grids is constructed. The vulnerability index is based on the premise that forests in India are already subjected to multiple stresses including over extraction, insect outbreaks, live¬stock grazing, forest fires and other anthropogenic pressures -with climate change being an additional stress. The forest vulnerability index suggests that nearly 39% of the forest grids in India are projected to be vulnerable to the impacts of climate change under the A2 scenario, while 34% of the forest grids are projected to be vulnerable under the B2 scenario. The vulnerability index suggests that forests in the central part of India, a significant part of the western Himalayan forests and northern and central parts of the Western Ghats are particularly vulnerable to the impacts of climate change. Forests in the northeastern part of India are seemingly resilient to the impacts of climate change. It also suggests that given the high deforestation rate in northeast, this region be prioritized for reducing deforestation and forest degradation (REDD) projects under the United Nations Framework Convention on Climate Change (UNFCCC) mechanisms.

There is ample growth in scientific evidence about climate change. Since, hydrometeorological processes are sensitive to climate variability and changes, ascertaining the linkages and feedbacks between the climate and the hydrometeorological processes becomes critical for environmental quality, economic development, social well-being etc. As the river basin integrates some of the important systems like ecological and socio-economic systems, the knowledge of plausible implications of climate change on hydrometeorology of a river basin will not only increase the awareness of how the hydrological systems may change over the coming century, but also prepare us for adapting to the impacts of climate changes on water resources for sustainable management and development. In general, quantitative climate impact studies are based on several meteorological variables and possible future climate scenarios. Among the meteorological variables, sic “cardinal” variables are identified as the most commonly used in impact studies (IPCC, 2001). These are maximum and minimum temperatures, precipitation, solar radiation, relative humidity and wind speed. The climate scenarios refer to plausible future climates, which have been constructed for explicit use for investigating the potential consequences of anthropogenic climate alterations, in addition to the natural climate variability. Among the climate scenarios adapted in impact assessments, General circulation model(GCM) projections based on marker scenarios given in Intergovernmental Panel on Climate Change’s (IPCC’s) Special Report on Emissions Scenarios(SRES) have become the standard scenarios. The GCMs are run at coarse resolutions and therefore the output climate variables for the various scenarios of these models cannot be used directly for impact assessment on a local(river basin)scale. Hence in the past, several methodologies such as downscaling and disaggregation have been developed to transfer information of atmospheric variables from the GCM scale to that of surface meteorological variables at local scale. The most commonly used downscaling approaches are based on transfer functions to represent the statistical relationships between the large scale atmospheric variables(predictors) and the local surface variables(predictands). Recently Support vector machine (SVM) is proposed, and is theoretically proved to have advantages over other techniques in use such as transfer functions. The SVM implements the structural risk minimization principle, which guarantees the global optimum solution. Further, for SVMs, the learning algorithm automatically decides the model architecture. These advantages make SVM a plausible choice for use in downscaling hydrometeorological variables. The literature review on use of transfer function for downscaling revealed that though a diverse range of transfer functions has been adopted for downscaling, only a few studies have evaluated the sensitivity of such downscaling models. Further, no studies have so far been carried out in India for downscaling hydrometeorological variables to a river basin scale, nor there was any prior work aimed at downscaling CGCM3 simulations to these variables at river basin scale for various IPCC SRES emission scenarios. The research presented in the thesis is motivated to assess the impact of climate change on streamflow at river basin scale for the various IPCC SRES scenarios (A1B, A2, B1 and COMMIT), by integrating implications of climate change on all the six cardinal variables. The catchment of Malaprabha river (upstream of Malaprabha reservoir) in India is chosen as the study area to demonstrate the effectiveness of the developed models, as it is considered to be a climatically sensitive region, because though the river originates in a region having high rainfall it feeds arid and semi-arid regions downstream. The data of the National Centers for Environmental Prediction (NCEP), the third generation Canadian Global Climate Model (CGCM3) of the Canadian Center for Climate Modeling and Analysis (CCCma), observed hydrometeorological variables, Digital Elevation model (DEM), land use/land cover map, and soil map prepared based on PAN and LISS III merged, satellite images are considered for use in the developed models. The thesis is broadly divided into four parts. The first part comprises of general introduction, data, techniques and tools used. The second part describes the process of assessment of the implications of climate change on monthly values of each of the six cardinal variables in the study region using SVM downscaling models and k-nearest neighbor (k-NN) disaggregation technique. Further, the sensitivity of the SVM downscaling models to the choice of predictors, predictand, calibration period, season and location is evaluated. The third part describes the impact assessment of climate change on streamflow in the study region using the SWAT hydrologic model, and SVM downscaling models. The fourth part presents summary of the work presented in the thesis, conclusions draws, and the scope for future research. The development of SVM downscaling model begins with the selection of probable predictors (large scale atmospheric variables). For this purpose, the cross-correlations are computed between the probable predictor variables in NCEP and GCM data sets, and the probable predictor variables in NCEP data set and the predictand. A pool of potential predictors is then stratified (which is optional and variable dependant) based on season and or location by specifying threshold values for the computed cross-correlations. The data on potential predictors are first standardized for a baseline period to reduce systemic bias (if any) in the mean and variance of predictors in GCM data, relative to those of the same in NCEP reanalysis data. The standardized NCEP predictor variables are then processed using principal component analysis (PCA) to extract principal components (PCs) which are orthogonal and which preserve more than 98% of the variance originally present in them. A feature vector is formed for each month using the PCs. The feature vector forms the input to the SVM model, and the contemporaneous value of predictand is its output. Finally, the downscaling model is calibrated to capture the relationship between NCEP data on potential predictors (i.e feature vectors) and the predictand. Grid search procedure is used to find the optimum range for each of the parameters. Subsequently, the optimum values of parameters are obtained from the selected ranges, using the stochastic search technique of genetic algorithm. The SVM model is subsequently validated, and then used to obtain projections of predictand for simulations of CGCM3. Results show that precipitation, maximum and minimum temperature, relative humidity and cloud cover are projected to increase in future for A1B, A2 and B1 scenarios, whereas no trend is discerned with theCOMMIT. The projected increase in predictands is high for A2 scenario and is least for B1 scenario. The wind speed is not projected to change in future for the study region for all the aforementioned scenarios. The solar radiation is projected to decrease in future for A1B, A2 and B1 scenarios, whereas no trend is discerned with the COMMIT. To assess the monthly streamflow responses to climate change, two methodologies are considered in this study namely (i) downscaling and disaggregating the meteorological variables for use as inputs in SWAT and (ii) directly downscaling streamflow using SVM. SWAT is a physically based, distributed, continuous time hydrological model that operates on a daily time scale. The hydrometeorologic variables obtained using SVM downscaling models are disaggregated to daily scale by using k-nearest neighbor method developed in this study. The other inputs to SWAT are DEM, land use/land cover map, soil map, which are considered to be the same for the present and future scenarios. The SWAT model has projected an increase in future streamflows for A1B, A2 andB1 scenarios, whereas no trend is discerned with the COMMIT. The monthly projections of streamflow at river basin scale are also obtained using two SVM based downscaling models. The first SVM model (called one-stage SVM model) considered feature vectors prepared based on monthly values of large scale atmospheric variables as inputs, whereas the second SVM model (called two-stage SVM model) considered feature vectors prepared from the monthly projections of cardinal variables as inputs. The trend in streamflows projected using two-stage SVM model is found to be similar to that projected by SWAT for each of the scenarios considered. The streamflow is not projected to change for any of the scenarios considered with the one-stage SVM downscaling model. The relative performance of the SWAT and the two SVM downscaling models in simulating observed streamflows is evaluated. In general, all the three models are able to simulate the streamflows well. Nevertheless, the performance of SWAT model is better. Further, among the two SVM models, the performance of one-stage streamflow downscaling model is marginally better than that of the two-stage streamflow downscaling model.

As a result of increase in global average surface temperature, abnormalities in different hydroclimatic components such as evapotranspiration, stream flow and precipitation have been experienced. So investigation has to be carried out to assess the hidden abnormality subsisting in the hydroclimatological time series in the form of trend. This thesis broadly consists of following four parts. The first part comprises of a detailed review of various trend detection approaches. Approaches incorporating the effect of serial correlation for trend detection and interesting developments concerning various non parametric approaches are focused explicitly. Recent trends in annual, monthly, and seasonl (winter, pre-monsoon, monsoon and post-monsoon) Tmax and Tmin have been analyzed considering three time slots viz. 1901-2003, 1948-2003 and 1970-2003. For this purpose, time series of Tmax and Tmin of India as a whole and for seven homogeneous regions, viz. Western Himalaya (WH), Northwest (NW), Northeast (NE), North Central (NC), East coast (EC), West coast (WC) and Interior Peninsula (IP) were originally considered. During the last three decades significant upward trend in Tmin is found to be present in all regions considered either at annual or seasonal level. Sequential Mann Kendall test revealed that most of the significant upward trends both in Tmax and Tmin began after 1970. The second part discusses about numerous climate models from both Coupled Model Inter comparison Project-5 and 3 (i.e. CMIP5, CMIP3) and their skills in simulating Indian climate and assessing their performance using various evaluation measures. Performances of climate models were evaluated for whole of India and over all the individual grid points covering India. The newly defined metric symbolized as Skill_All is an intersection of the three metrics i.e. Skill_r, Skill_s and Skill_rmse, is used for overall model evaluation analysis. A notable enhancement of Skill_All for CMIP5 over CMIP3 was found. After overall model evaluation study, Compromise Programming, a distance based decision making technique, was employed to rank the GCMs gridwise. Entropy method was employed to obtain weights of the chosen indicators. Group decision making methodology was used to arrive at a consensus based on the ranking pattern obtained by individual grid points. In the third part, a detailed detection and attribution (D&A) analysis is performed to determine the causes of changes in seasonal Tmax and Tmin during the period 1950-2005. This formal D&A exercise helps in providing better insight (than trend detection analysis) into the nature of the observed seasonal temperature changes. It was noticed that the emergence of observed trend was more pronounced in Tmin compared to Tmax. Although observed changes were not solely associated with one specific causative factor, most of the changes in Tmin are above the bounds of natural internal climate variability. Finally in the fourth part, to understand the climate change impact on the hydrological cycle, a spatiotemporal change detection study of potential evapotranspiration (PET) along with Tmax and Tmin over India has been performed. Climatology patterns for PET confirmed a greater PET rate during the month of March, April, May and June. A significant increasing trend in both Tmax and Tmin (Tmin being more) was observed in more number of grid points compared to PET. Significant positive trends in Tmax, Tmin and PET were observed over most of the grid points in the IP region. Heterogeneities existed in the spatiotemporal variability of PET over all India. This spatio-temporal change detection study would be helpful for present and future water resources management.