River water quality response under hypothetical climate change scenarios in Tunga-Bhadra river, India

Analysis of climate change impacts on streamflow by perturbing the climate inputs has been a concern for many authors in the past few years, but there are few analyses for the impacts on water quality. To examine the impact of change in climate variables on the water quality parameters, the water quality input variables have to be perturbed. The primary input variables that can be considered for such an analysis are streamflow and water temperature, which are affected by changes in precipitation and air temperature, respectively. Using hypothetical scenarios to represent both greenhouse warming and streamflow changes, the sensitivity of the water quality parameters has been evaluated under conditions of altered river flow and river temperature in this article. Historical data analysis of hydroclimatic variables is carried out, which includes flow duration exceedance percentage (e.g. Q90), single low- flow indices (e.g. 7Q10, 30Q10) and relationships between climatic variables and surface variables. For the study region of Tunga-Bhadra river in India, low flows are found to be decreasing and water temperatures are found to be increasing. As a result, there is a reduction in dissolved oxygen (DO) levels found in recent years. Water quality responses of six hypothetical climate change scenarios were simulated by the water quality model, QUAL2K. A simple linear regression relation between air and water temperature is used to generate the scenarios for river water temperature. The results suggest that all the hypothetical climate change scenarios would cause impairment in water quality. It was found that there is a significant decrease in DO levels due to the impact of climate change on temperature and flows, even when the discharges were at safe permissible levels set by pollution control agencies (PCAs). The necessity to improve the standards of PCA and develop adaptation policies for the dischargers to account for climate change is examined through a fuzzy waste load allocation model developed earlier. Copyright (C) 2011 John Wiley & Sons, Ltd.

Climate change and crop yields in the Indian Cardamom Hills, 1978-2007 CE

In this study we analyzed climate and crop yields data from Indian cardamom hills for the period 1978-2007 to investigate whether there were significant changes in weather elements, and if such changes have had significant impact on the production of spices and plantation crops. Spatial and temporal variations in air temperatures (maximum and minimum), rainfall and relative humidity are evident across stations. The mean air temperature increased significantly during the last 30 years; the greatest increase and the largest significant upward trend was observed in the daily temperature. The highest increase in minimum temperature was registered for June (0.37A degrees C/18 years) at the Myladumpara station. December and January showed greater warming across the stations. Rainfall during the main monsoon months (June-September) showed a downward trend. Relative humidity showed increasing and decreasing trends, respectively, at the cardamom and tea growing tracts. The warming trend coupled with frequent wet and dry spells during the summer is likely to have a favorable effect on insect pests and disease causing organisms thereby pesticide consumption can go up both during excess rainfall and drought years. The incidence of many minor pest insects and disease pathogens has increased in the recent years of our study along with warming. Significant and slight increases in the yield of small cardamom (Elettaria cardamomum M.) and coffee (Coffea arabica), respectively, were noticed in the recent years.; however the improvement of yield in tea (Thea sinensis) and black pepper (Piper nigrum L.) has not been seen in our analysis.

Controls on intense silicate weathering in a tropical river, southwestern India

The Silicate Weathering Rate (SWR) and associated Carbon dioxide Consumption Rate (CCR) in tropical silicate terrain is assessed through a study of the major ion chemistry in a small west flowing river of Peninsular India, the Nethravati River. The specific features of the river basin are high mean annual rainfall and temperature, high runoff and a Precambrian basement composed of granitic-gneiss, charnockite and minor metasediments. The water samples (n = 56) were collected from three locations along the Nethravati River and from two of its tributaries over a period of twelve months. Chemical Weathering Rate (CWR) for the entire watershed is calculated by applying rainwater correction using river chloride as a tracer. Chemical Weathering Rate in the Nethravati watershed is estimated to 44 t.km(-2).y(-1) encompassing a SWR of 42 t.km(-2).y(-1) and a maximum carbonate contribution of 2 t.km(-2).y(-1). This SWR is among the highest reported for granito-gneissic terrains. The assessed CCR is 2.9 . 10(5) mol.km(-2).y(-1). The weathering index (Re). calculated from molecular ratios of dissolved cations and silica in the river, suggests an intense silicate weathering leading to kaolinite-gibbsite precipitation in the weathering covers. The intense SWR and CCR could be due to the combination of high runoff and temperature along with the thickness and nature of the weathering cover. The comparison of silicate weathering fluxes with other watersheds reveals that under similar morpho-climatic settings basalt weathering would be 2.5 times higher than the granite-gneissic rocks. (C) 2012 Elsevier B.V. All rights reserved.

Insights to urban dynamics through landscape spatial pattern analysis

Urbanisation is a dynamic complex phenomenon involving large scale changes in the land uses at local levels. Analyses of changes in land uses in urban environments provide a historical perspective of land use and give an opportunity to assess the spatial patterns, correlation, trends, rate and impacts of the change, which would help in better regional planning and good governance of the region. Main objective of this research is to quantify the urban dynamics using temporal remote sensing data with the help of well-established landscape metrics. Bangalore being one of the rapidly urbanising landscapes in India has been chosen for this investigation. Complex process of urban sprawl was modelled using spatio temporal analysis. Land use analyses show 584% growth in built-up area during the last four decades with the decline of vegetation by 66% and water bodies by 74%. Analyses of the temporal data reveals an increase in urban built up area of 342.83% (during 1973-1992), 129.56% (during 1992-1999), 106.7% (1999-2002), 114.51% (2002-2006) and 126.19% from 2006 to 2010. The Study area was divided into four zones and each zone is further divided into 17 concentric circles of 1 km incrementing radius to understand the patterns and extent of the urbanisation at local levels. The urban density gradient illustrates radial pattern of urbanisation for the period 1973-2010. Bangalore grew radially from 1973 to 2010 indicating that the urbanisation is intensifying from the central core and has reached the periphery of the Greater Bangalore. Shannon's entropy, alpha and beta population densities were computed to understand the level of urbanisation at local levels. Shannon's entropy values of recent time confirms dispersed haphazard urban growth in the city, particularly in the outskirts of the city. This also illustrates the extent of influence of drivers of urbanisation in various directions. Landscape metrics provided in depth knowledge about the sprawl. Principal component analysis helped in prioritizing the metrics for detailed analyses. The results clearly indicates that whole landscape is aggregating to a large patch in 2010 as compared to earlier years which was dominated by several small patches. The large scale conversion of small patches to large single patch can be seen from 2006 to 2010. In the year 2010 patches are maximally aggregated indicating that the city is becoming more compact and more urbanised in recent years. Bangalore was the most sought after destination for its climatic condition and the availability of various facilities (land availability, economy, political factors) compared to other cities. The growth into a single urban patch can be attributed to rapid urbanisation coupled with the industrialisation. Monitoring of growth through landscape metrics helps to maintain and manage the natural resources. (C) 2012 Elsevier B.V. All rights reserved.

Causality analysis of groundwater dynamics based on a Vector Autoregressive model in the semi-arid basin of Gundal (South India)

Causal relationships existing between observed levels of groundwater in a semi-arid sub-basin of the Kabini River basin (Karnataka state, India) are investigated in this study. A Vector Auto Regressive model is used for this purpose. Its structure is built on an upstream/downstream interaction network based on observed hydro-physical properties. Exogenous climatic forcing is used as an input based on cumulated rainfall departure. Optimal models are obtained thanks to a trial approach and are used as a proxy of the dynamics to derive causal networks. It appears to be an interesting tool for analysing the causal relationships existing inside the basin. The causal network reveals 3 main regions: the Northeastern part of the Gundal basin is closely coupled to the outlet dynamics. The Northwestern part is mainly controlled by the climatic forcing and only marginally linked to the outlet dynamic. Finally, the upper part of the basin plays as a forcing rather than a coupling with the lower part of the basin allowing for a separate analysis of this local behaviour. The analysis also reveals differential time scales at work inside the basin when comparing upstream oriented with downstream oriented causalities. In the upper part of the basin, time delays are close to 2 months in the upward direction and lower than 1 month in the downward direction. These time scales are likely to be good indicators of the hydraulic response time of the basin which is a parameter usually difficult to estimate practically. This suggests that, at the sub-basin scale, intra-annual time scales would be more relevant scales for analysing or modelling tropical basin dynamics in hard rock (granitic and gneissic) aquifers ubiquitous in south India. (c) 2012 Elsevier B.V. All rights reserved.

Parameter estimation of a two-horizon soil profile by combining crop canopy and surface soil moisture observations using GLUE

Estimation of soil parameters by inverse modeling using observations on either surface soil moisture or crop variables has been successfully attempted in many studies, but difficulties to estimate root zone properties arise when heterogeneous layered soils are considered. The objective of this study was to explore the potential of combining observations on surface soil moisture and crop variables - leaf area index (LAI) and above-ground biomass for estimating soil parameters (water holding capacity and soil depth) in a two-layered soil system using inversion of the crop model STICS. This was performed using GLUE method on a synthetic data set on varying soil types and on a data set from a field experiment carried out in two maize plots in South India. The main results were (i) combination of surface soil moisture and above-ground biomass provided consistently good estimates with small uncertainity of soil properties for the two soil layers, for a wide range of soil paramater values, both in the synthetic and the field experiment, (ii) above-ground biomass was found to give relatively better estimates and lower uncertainty than LAI when combined with surface soil moisture, especially for estimation of soil depth, (iii) surface soil moisture data, either alone or combined with crop variables, provided a very good estimate of the water holding capacity of the upper soil layer with very small uncertainty whereas using the surface soil moisture alone gave very poor estimates of the soil properties of the deeper layer, and (iv) using crop variables alone (else above-ground biomass or LAI) provided reasonable estimates of the deeper layer properties depending on the soil type but provided poor estimates of the first layer properties. The robustness of combining observations of the surface soil moisture and the above-ground biomass for estimating two layer soil properties, which was demonstrated using both synthetic and field experiments in this study, needs now to be tested for a broader range of climatic conditions and crop types, to assess its potential for spatial applications. (C) 2012 Elsevier B.V. All rights reserved.

A comparative study of deglaciation in two neighbouring basins (Warwan and Bhut) of Western Himalaya

Glaciers of the Himalaya contribute significantly in the processes linking atmosphere, biosphere and hydrosphere, thus need to be monitored in view of the climatic variations. In this direction, many studies have been carried out during the last two decades and satellite-based multispectral data have been used extensively for this purpose throughout the world. The present study is aimed at mapping of glaciers in two adjacent basins (Warwan and Bhut) of the Western Himalaya with almost similar altitude and latitude and comparing the changes in the two time-frames with respect to three parameters, i.e. area, debris cover and area altitude distribution of glaciers. The two time-frames are topographical maps of 1962 and IRS LISS III images of 2001/02. Deglaciation was observed in both the basins with 19% and 9% loss in the glaciated area in Warwan and Bhut respectively. This difference may be due to: (i) the smaller size of the glaciers of the Warwan Basin (e.g. 164 glaciers having <1 sq. km area in comparison to 101 glaciers in the Bhut Basin), (ii) lower percentage of moraine cover in Warwan (18) than in the Bhut Basin (30) and (iii) higher percentage of glaciated area lying below 5100 m (80) in Warwan than in the Bhut Basin (70).

Climate response to changes in atmospheric carbon dioxide and solar irradiance on the time scale of days to weeks

Recent studies show that fast climate response on time scales of less than a month can have important implications for long-term climate change. In this study, we investigate climate response on the time scale of days to weeks to a step-function quadrupling of atmospheric CO2 and contrast this with the response to a 4% increase in solar irradiance. Our simulations show that significant climate effects occur within days of a stepwise increase in both atmospheric CO2 content and solar irradiance. Over ocean, increased atmospheric CO2 warms the lower troposphere more than the surface, increasing atmospheric stability, moistening the boundary layer, and suppressing evaporation and precipitation. In contrast, over ocean, increased solar irradiance warms the lower troposphere to a much lesser extent, causing a much smaller change in evaporation and precipitation. Over land, both increased CO2 and increased solar irradiance cause rapid surface warming that tends to increase both evaporation and precipitation. However, the physiological effect of increased atmospheric CO2 on plant stomata reduces plant transpiration, drying the boundary layer and decreasing precipitation. This effect does not occur with increased solar irradiance. Therefore, differences in climatic effects from CO2 versus solar forcing are manifested within days after the forcing is imposed.

Longitudinal localization of tropical intraseasonal variability

Intraseasonal time-scales play an important role in tropical variability. Two modes that contribute significantly to tropical intraseasonal variability (ISV) are the eastward-propagating MaddenJulian Oscillation (MJO), and westward-moving moist equatorial Rossby waves. This note reports on a correspondence between the longitudinal gradient of mean tropical precipitable water (PW), and the geographical regions of genesis, and convective activity, of both these large-scale tropical systems. Our finding is based on an analysis of PW from the MERRA reanalysis product. The data indicate that the mean tropical PW has a dominant wavenumber two (three) structure in longitude in the Northern (Southern) Hemisphere. Departures from a longitudinally homogeneous state are attributed to the influence of subtropical anticyclones, and are accentuated by monsoonal regions of both hemispheres. This mean structure results in a sharply localized longitudinal gradient of PW. Remarkably, regions with positive gradients (such as the Northern and Southern Hemisphere western Indian Ocean), i.e. they have larger PW to the east, are the very zones that are implicated in the formation, and show high levels of convective activity, of the eastward-moving MJO. On the other hand, regions with negative gradients (such as the Southern Hemisphere central Pacific) are the very regions where genesis, and maxima in variance, of westward-moving moist equatorial Rossby waves are known to occur. Apart from providing a first-order longitudinal footprint of the convective phase of these systems, this correspondence reinforces the role of the mean climatic state in tropical ISV. Copyright (c) 2012 Royal Meteorological Society

Assessing severe Drought and wet events over India in a future climate using a nested bias-correction approach

General circulation models (GCMs) are routinely used to simulate future climatic conditions. However, rainfall outputs from GCMs are highly uncertain in preserving temporal correlations, frequencies, and intensity distributions, which limits their direct application for downscaling and hydrological modeling studies. To address these limitations, raw outputs of GCMs or regional climate models are often bias corrected using past observations. In this paper, a methodology is presented for using a nested bias-correction approach to predict the frequencies and occurrences of severe droughts and wet conditions across India for a 48-year period (2050-2099) centered at 2075. Specifically, monthly time series of rainfall from 17 GCMs are used to draw conclusions for extreme events. An increasing trend in the frequencies of droughts and wet events is observed. The northern part of India and coastal regions show maximum increase in the frequency of wet events. Drought events are expected to increase in the west central, peninsular, and central northeast regions of India. (C) 2013 American Society of Civil Engineers.

A note on new indices for the equatorial Indian Ocean oscillation

It is now well known that there is a strong association of the extremes of the Indian summer monsoon rainfall (ISMR) with the El Nio and southern oscillation (ENSO) and the Equatorial Indian Ocean Oscillation (EQUINOO), later being an east-west oscillation in convection anomaly over the equatorial Indian Ocean. So far, the index used for EQUINOO is EQWIN, which is based on the surface zonal wind over the central equatorial Indian Ocean. Since the most important attribute of EQUINOO is the oscillation in convection/precipitation, we believe that the indices based on convection or precipitation would be more appropriate. Continuous and reliable data on outgoing longwave radiation (OLR), and satellite derived precipitation are now available from 1979 onwards. Hence, in this paper, we introduce new indices for EQUINOO, based on the difference in the anomaly of OLR/precipitation between eastern and western parts of the equatorial Indian Ocean. We show that the strong association of extremes of the Indian summer monsoon with ENSO and EQUINOO is also seen when the new indices are used to represent EQUINOO.

Flying between Sky Islands: The Effect of Naturally Fragmented Habitat on Butterfly Population Structure

High elevation montane areas are called ``sky islands'' when they occur as a series of high mountains separated by lowland valleys. Different climatic conditions at high elevations makes sky islands a specialized type of habitat, rendering them naturally fragmented compared to more continuous habitat at lower elevations. Species in sky islands face unsuitable climate in the intervening valleys when moving from one montane area to another. The high elevation shola-grassland mosaic in the Western Ghats of southern India form one such sky island complex. The fragmented patches make this area ideal to study the effect of the spatial orientation of suitable habitat patches on population genetic structure of species found in these areas. Past studies have suggested that sky islands tend to have genetically structured populations, possibly due to reduced gene flow between montane areas. To test this hypothesis, we adopted the comparative approach. Using Amplified Fragment Length Polymorphisms, we compared population genetic structures of two closely related, similar sized butterfly species: Heteropsis oculus, a high elevation shola-grassland specialist restricted to the southern Western Ghats, and Mycalesis patnia, found more continuously distributed in lower elevations. In all analyses, as per expectation the sky island specialist H. oculus exhibited a greater degree of population genetic structure than M. patnia, implying a difference in geneflow. This difference in geneflow in turn appears to be due to the natural fragmentation of the sky island complexes. Detailed analysis of a subset of H. oculus samples from one sky island complex (the Anamalais) showed a surprising genetic break. A possible reason for this break could be unsuitable conditions of higher temperature and lower rainfall in the intervening valley region. Thus, sky island species are not only restricted by lack of habitat continuity between montane areas, but also by the nature of the intervening habitat.

The influence of debris cover and glacial lakes on the recession of glaciers in Sikkim Himalaya, India

We investigated area changes in glaciers covering an area of similar to 200 km(2) in the Tista basin, Sikkim, Eastern Indian Himalaya, between similar to 1990 and 2010 using Landsat Thematic Mapper (TM) and Indian Remote-sensing Satellite (IRS) images and related the changes to debris cover, supraglacial lakes and moraine-dam lakes. The glaciers lost an area of 3.3 +/- 0.8% between 1989/90 and 2010. More detailed analysis revealed an area loss of 2.00 +/- 0.82, 2.56 +/- 0.61 and 2.28 +/- 2.01 km(2) for the periods 1989-97, 1997-2004/05 and 2004-2009/10, respectively. This indicates an accelerated retreat of glaciers after 1997. On further analysis, we observed (1) the formation and expansion of supraglacial lakes on many debris-covered glaciers and (2) the merging of these lakes over time, leading to the development of large moraine-dam lakes. We also observed that debris-covered glaciers with lakes lose a greater area than debris-covered glaciers without lakes and debris-free glaciers. The climatic data for 24 years (1987-2011), measured at the Gangtok meteorological station (1812 m a.s.l.), showed that the region experienced a 1.0 degrees C rise in the summer minimum temperature and a 2.0 degrees C rise in the winter minimum temperature, indicating hotter summers and warmer winters. There was no significant trend in the total annual precipitation. We find that glacier retreat is caused mainly by a temperature increase and that debris-covered glaciers can retreat at a faster rate than debris-free glaciers, if associated with lakes.

A simple model for spatial disaggregation of evaporative fraction: Comparative study with thermal sharpened land surface temperature data over India

[1] Evaporative fraction (EF) is a measure of the amount of available energy at the earth surface that is partitioned into latent heat flux. The currently operational thermal sensors like the Moderate Resolution Imaging Spectroradiometer (MODIS) on satellite platforms provide data only at 1000 m, which constraints the spatial resolution of EF estimates. A simple model (disaggregation of evaporative fraction (DEFrac)) based on the observed relationship between EF and the normalized difference vegetation index is proposed to spatially disaggregate EF. The DEFrac model was tested with EF estimated from the triangle method using 113 clear sky data sets from the MODIS sensor aboard Terra and Aqua satellites. Validation was done using the data at four micrometeorological tower sites across varied agro-climatic zones possessing different land cover conditions in India using Bowen ratio energy balance method. The root-mean-square error (RMSE) of EF estimated at 1000 m resolution using the triangle method was 0.09 for all the four sites put together. The RMSE of DEFrac disaggregated EF was 0.09 for 250 m resolution. Two models of input disaggregation were also tried with thermal data sharpened using two thermal sharpening models DisTrad and TsHARP. The RMSE of disaggregated EF was 0.14 for both the input disaggregation models for 250 m resolution. Moreover, spatial analysis of disaggregation was performed using Landsat-7 (Enhanced Thematic Mapper) ETM+ data over four grids in India for contrasted seasons. It was observed that the DEFrac model performed better than the input disaggregation models under cropped conditions while they were marginally similar under non-cropped conditions.

Observed changes in Himalayan glaciers

In the Himalaya, large areas are covered by glaciers and seasonal snow. They are an important source of water for the Himalayan rivers. In this article, observed changes in glacial extent and mass balance have been discussed. Various studies suggest that most of the Himalayan glaciers are retreating though the rate of retreat varies from glacier to glacier, ranging from a few meters to almost 61 m/year, depending upon the terrain and meteorological parameters. In addition, mapping of almost 11,000 out of 40,000 sq. km of glaciated area, distributed in all major climatic zones of the Himalaya, suggests an almost 13% loss in area in the last 4-5 decades. The glacier mass balance observations and estimates made using methods like field, AAR, ELA and geodetic measurements, suggest a significant increase in mass wastage of Himalayan glaciers in the last 3-4 decades. In the last four decades loss in glacial ice has been estimated at 19 +/- 7 m. This suggests loss of 443 +/- 136 Gt of glacial mass out of a total 3600-4400 Gt of glacial stored water in the Indian Himalaya. This study has also shown that mean loss in glacier mass in the Indian Himalaya is accelerated from -9 +/- 4 to -20 +/- 4 Gt/year between the periods 1975-85 and 2000-2010. The estimate of glacial stored water in the Indian Himalaya is based on glacier inventory on a 1 : 250,000 scale and scaling methods; therefore, we assume uncertainties to be large.