On the basin-scale detection and attribution of human-induced climate change in monsoon precipitation and streamflow

Detecting and quantifying the presence of human-induced climate change in regional hydrology is important for studying the impacts of such changes on the water resources systems as well as for reliable future projections and policy making for adaptation. In this article a formal fingerprint-based detection and attribution analysis has been attempted to study the changes in the observed monsoon precipitation and streamflow in the rain-fed Mahanadi River Basin in India, considering the variability across different climate models. This is achieved through the use of observations, several climate model runs, a principal component analysis and regression based statistical downscaling technique, and a Genetic Programming based rainfall-runoff model. It is found that the decreases in observed hydrological variables across the second half of the 20th century lie outside the range that is expected from natural internal variability of climate alone at 95% statistical confidence level, for most of the climate models considered. For several climate models, such changes are consistent with those expected from anthropogenic emissions of greenhouse gases. However, unequivocal attribution to human-induced climate change cannot be claimed across all the climate models and uncertainties in our detection procedure, arising out of various sources including the use of models, cannot be ruled out. Changes in solar irradiance and volcanic activities are considered as other plausible natural external causes of climate change. Time evolution of the anthropogenic climate change ``signal'' in the hydrological observations, above the natural internal climate variability ``noise'' shows that the detection of the signal is achieved earlier in streamflow as compared to precipitation for most of the climate models, suggesting larger impacts of human-induced climate change on streamflow than precipitation at the river basin scale.

Physically Based, Hydrologic Model Results Based on Three Precipitation Products

The main objective of the study is to examine the accuracy of and differences among simulated streamflows driven by rainfall estimates from a network of 22 rain gauges spread over a 2,170 km2 watershed, NEXRAD Stage III radar data, and Tropical Rainfall Measuring Mission (TRMM) 3B42 satellite data. The Gridded Surface Subsurface Hydrologic Analysis (GSSHA), a physically based, distributed parameter, grid-structured, hydrologic model, was used to simulate the June-2002 flooding event in the Upper Guadalupe River watershed in south central Texas. There were significant differences between the rainfall fields estimated by the three types of measurement technologies. These differences resulted in even larger differences in the simulated hydrologic response of the watershed. In general, simulations driven by radar rainfall yielded better results than those driven by satellite or rain-gauge estimates. This study also presents an overview of effects of land cover changes on runoff and stream discharge. The results demonstrate that, for major rainfall events similar to the 2002 event, the effect of urbanization on the watershed in the past two decades would not have made any significant effect on the hydrologic response. The effect of urbanization on the hydrologic response increases as the size of the rainfall event decreases.

Assessment of hydrologic impacts of climate change in Tunga-Bhadra river basin, India with HEC-HMS and SDSM

Climate change would significantly affect many hydrologic systems, which in turn would affect the water availability, runoff, and the flow in rivers. This study evaluates the impacts of possible future climate change scenarios on the hydrology of the catchment area of the TungaBhadra River, upstream of the Tungabhadra dam. The Hydrologic Engineering Center's Hydrologic Modeling System version 3.4 (HEC-HMS 3.4) is used for the hydrological modelling of the study area. Linear-regression-based Statistical DownScaling Model version 4.2 (SDSM 4.2) is used to downscale the daily maximum and minimum temperature, and daily precipitation in the four sub-basins of the study area. The large-scale climate variables for the A2 and B2 scenarios obtained from the Hadley Centre Coupled Model version 3 are used. After model calibration and testing of the downscaling procedure, the hydrological model is run for the three future periods: 20112040, 20412070, and 20712099. The impacts of climate change on the basin hydrology are assessed by comparing the present and future streamflow and the evapotranspiration estimates. Results of the water balance study suggest increasing precipitation and runoff and decreasing actual evapotranspiration losses over the sub-basins in the study area.

Modeling the vulnerability of an urban groundwater system due to the combined impacts of climate change and management scenarios

Climate change impact on a groundwater-dependent small urban town has been investigated in the semiarid hard rock aquifer in southern India. A distributed groundwater model was used to simulate the groundwater levels in the study region for the projected future rainfall (2012-32) obtained from a general circulation model (GCM) to estimate the impacts of climate change and management practices on groundwater system. Management practices were based on the human-induced changes on the urban infrastructure such as reduced recharge from the lakes, reduced recharge from water and wastewater utility due to an operational and functioning underground drainage system, and additional water extracted by the water utility for domestic purposes. An assessment of impacts on the groundwater levels was carried out by calibrating a groundwater model using comprehensive data gathered during the period 2008-11 and then simulating the future groundwater level changes using rainfall from six GCMs Institute of Numerical Mathematics Coupled Model, version 3.0 (INM-CM. 3.0); L'Institut Pierre-Simon Laplace Coupled Model, version 4 (IPSL-CM4); Model for Interdisciplinary Research on Climate, version 3.2 (MIROC3.2); ECHAM and the global Hamburg Ocean Primitive Equation (ECHO-G); Hadley Centre Coupled Model, version 3 (HadCM3); and Hadley Centre Global Environment Model, version 1 (HadGEM1)] that were found to show good correlation to the historical rainfall in the study area. The model results for the present condition indicate that the annual average discharge (sum of pumping and natural groundwater outflow) was marginally or moderately higher at various locations than the recharge and further the recharge is aided from the recharge from the lakes. Model simulations showed that groundwater levels were vulnerable to the GCM rainfall and a scenario of moderate reduction in recharge from lakes. Hence, it is important to sustain the induced recharge from lakes by ensuring that sufficient runoff water flows to these lakes.

Development of a regional model for the North Indian Ocean

We have developed a one-way nested Indian Ocean regional model. The model combines the National Oceanic and Atmospheric Administration (NOAA) Geophysical Fluid Dynamics Laboratory's (GFDL) Modular Ocean Model (MOM4p1) at global climate model resolution (nominally one degree), and a regional Indian Ocean MOM4p1 configuration with 25 km horizontal resolution and 1 m vertical resolution near the surface. Inter-annual global simulations with Coordinated Ocean-Ice Reference Experiments (CORE-II) surface forcing over years 1992-2005 provide surface boundary conditions. We show that relative to the global simulation, (i) biases in upper ocean temperature, salinity and mixed layer depth are reduced, (ii) sea surface height and upper ocean circulation are closer to observations, and (iii) improvements in model simulation can be attributed to refined resolution, more realistic topography and inclusion of seasonal river runoff. Notably, the surface salinity bias is reduced to less than 0.1 psu over the Bay of Bengal using relatively weak restoring to observations, and the model simulates the strong, shallow halocline often observed in the North Bay of Bengal. There is marked improvement in subsurface salinity and temperature, as well as mixed layer depth in the Bay of Bengal. Major seasonal signatures in observed sea surface height anomaly in the tropical Indian Ocean, including the coastal waveguide around the Indian peninsula, are simulated with great fidelity. The use of realistic topography and seasonal river runoff brings the three dimensional structure of the East India Coastal Current and West India Coastal Current much closer to observations. As a result, the incursion of low salinity Bay of Bengal water into the southeastern Arabian Sea is more realistic. (C) 2013 Elsevier Ltd. All rights reserved.

Hydrographical characteristics and oxygen isotopic signatures of water in a coastal environment (Mangalore) along the southeastern Arabian Sea

Coastal marine environments are important links between the continents and the open ocean. The coast off Mangalore forms part of the upwelling zone along the southeastern Arabian Sea. The temperature, salinity, density, dissolved oxygen and stable oxygen isotope ratio (delta O-18) of surface waters as well as those of bottom waters off coastal Mangalore were studied every month from October 2010 to May 2011. The coastal waters were stratified in October and November due to precipitation and runoff. The region was characterised by upwelled bottom waters in October, whereas the region exhibited a temperature inversion in November. The surface and bottom waters presented almost uniform properties from December until April. The coastal waters were observed to be most dense in January and May. Comparatively cold and poorly oxygenated bottom waters during the May sampling indicated the onset of upwelling along the region. delta O-18 of the coastal waters successfully documented the observed variations in the hydrographical characteristics of the Mangalore coast during the monthly sampling period. We also noted that the monthly variability in the properties of the coastal waters of Mangalore was related to the hydrographical characteristics of the adjacent open ocean inferred from satellite-derived surface winds, sea surface height anomaly data and sea surface temperatures.

Vegetation impact on stream chemical fluxes: Mule Hole watershed (South India)

The proportion of chemical elements passing through vegetation prior to being exported in a stream was quantified for a forested tropical watershed(Mule Hole, South India) using an extensive hydrological and geochemical monitoring at several scales. First, a solute annual mass balance was established at the scale of the soil-plant profile for assessing the contribution of canopy interaction and litter decay to the solute fluxes of soil inputs (overland flow) and soil outputs (pore water flow as seepages). Second, based on the respective contributions of overland flow and seepages to the stream flow as estimated by a hydrological lumped model, we assigned the proportion of chemical elements in the stream that transited through the vegetation at both flood event (End Member Mixing Analysis) and seasonal scales. At the scale of the 1D soil-plant profile, leaching from the canopy constituted the main source of K above the ground surface. Litter decay was the main source of Si, whereas alkalinity, Ca and Mg originated in the same proportions from both sources. The contribution of vegetation was negligible for Na. Within the soil, all elements but Na were removed from the pore water in proportions varying from 20% for Cl to 95% for K: The soil output fluxes corresponded to a residual fraction of the infiltration fluxes. The behavior of K, Cl, Ca and Mg in the soil-plant profile can be explained by internal cycling, as their soil output fluxes were similar to the atmospheric inputs. Na was released from soils as a result of Na-plagioclase weathering and accompanied by additional release of Si. Concentration of soil pore water by evapotranspiration might limit the chemical weathering in the soil. Overall, the solute K, Ca, Mg, alkalinity and Si fluxes associated with the vegetation turnover within the small experimental watershed represented 10-15 times the solute fluxes exported by the stream, of which 83-97% transited through the vegetation. One important finding is that alkalinity and Si fluxes at the outlet were not linked to the ``current weathering'' of silicates in this watershed. These results highlight the dual effect of the vegetation cover on the solute fluxes exported from the watershed: On one hand the runoff was limited by evapotranspiration and represented only 10% of the annual rainfall, while on the other hand, 80-90% of the overall solute flux exported by the stream transited through the vegetation. The approach combining geochemical monitoring and accurate knowledge of the watershed hydrological budget provided detailed understanding of several effects of vegetation on stream fluxes: (1) evapotranspiration (limiting), (2) vertical transfer through vegetation from vadose zone to ground surface (enhancing) and (3) redistribution by throughfalls and litter decay. It provides a good basis for calibrating geochemical models and more precisely assessing the role of vegetation on soil processes. (C) 2014 Elsevier Ltd. All rights reserved.

Combining data sets of satellite-retrieved products for basin-scale water balance study: 2. Evaluation on the Mississippi Basin and closure correction model

In this study, we applied the integration methodology developed in the companion paper by Aires (2014) by using real satellite observations over the Mississippi Basin. The methodology provides basin-scale estimates of the four water budget components (precipitation P, evapotranspiration E, water storage change Delta S, and runoff R) in a two-step process: the Simple Weighting (SW) integration and a Postprocessing Filtering (PF) that imposes the water budget closure. A comparison with in situ observations of P and E demonstrated that PF improved the estimation of both components. A Closure Correction Model (CCM) has been derived from the integrated product (SW+PF) that allows to correct each observation data set independently, unlike the SW+PF method which requires simultaneous estimates of the four components. The CCM allows to standardize the various data sets for each component and highly decrease the budget residual (P - E - Delta S - R). As a direct application, the CCM was combined with the water budget equation to reconstruct missing values in any component. Results of a Monte Carlo experiment with synthetic gaps demonstrated the good performances of the method, except for the runoff data that has a variability of the same order of magnitude as the budget residual. Similarly, we proposed a reconstruction of Delta S between 1990 and 2002 where no Gravity Recovery and Climate Experiment data are available. Unlike most of the studies dealing with the water budget closure at the basin scale, only satellite observations and in situ runoff measurements are used. Consequently, the integrated data sets are model independent and can be used for model calibration or validation.

Impact of compost, vermicompost and biochar on soil fertility, maize yield and soil erosion in Northern Vietnam: A three year mesocosm experiment

Compost, vermicompost and biochar amendments are thought to improve soil quality and plant yield. However, little is known about their long-term impact on crop yield and the environment in tropical agro-ecosystems. In this study we investigated the effect of organic amendments (buffalo manure, compost and verrnicompost) and biochar (applied alone or with vermicompost) on plant yield, soil fertility, soil erosion and water dynamics in a degraded Acrisol in Vietnam. Maize growth and yield, as well as weed growth, were examined for three years in terrestrial mesocosms under natural rainfall. Maize yield and growth showed high inter-annual variability depending on the organic amendment. Vermicompost improved maize growth and yield but its effect was rather small and was only significant when water availability was limited (year 2). This suggests that vermicompost could be a promising substrate for improving the resistance of agrosystems to water stress. When the vermicompost biochar mixture was applied, further growth and yield improvements were recorded in some cases. When applied alone, biochar had a positive influence on maize yield and growth, thus confirming its interest for improving long-term soil productivity. All organic amendments reduced water runoff, soil detachment and NH4+ and NO3- transfer to water. These effects were more significant with vermicompost than with buffalo manure and compost, highlighting that the beneficial influence of vermicompost is not limited to its influence on plant yield. In addition, this study showed for the first time that the combination of vermicompost and biochar may not only improve plant productivity but also reduce the negative impact of agriculture on water quality. (C) 2015 Elsevier B.V. All rights reserved.

Evaluating Methods to Predict Streamflow at Ungauged Sites using Regional Flow Duration Curves: A Case Study

Precise information on streamflows is of major importance for planning and monitoring of water resources schemes related to hydro power, water supply, irrigation, flood control, and for maintaining ecosystem. Engineers encounter challenges when streamflow data are either unavailable or inadequate at target locations. To address these challenges, there have been efforts to develop methodologies that facilitate prediction of streamflow at ungauged sites. Conventionally, time intensive and data exhaustive rainfall-runoff models are used to arrive at streamflow at ungauged sites. Most recent studies show improved methods based on regionalization using Flow Duration Curves (FDCs). A FDC is a graphical representation of streamflow variability, which is a plot between streamflow values and their corresponding exceedance probabilities that are determined using a plotting position formula. It provides information on the percentage of time any specified magnitude of streamflow is equaled or exceeded. The present study assesses the effectiveness of two methods to predict streamflow at ungauged sites by application to catchments in Mahanadi river basin, India. The methods considered are (i) Regional flow duration curve method, and (ii) Area Ratio method. The first method involves (a) the development of regression relationships between percentile flows and attributes of catchments in the study area, (b) use of the relationships to construct regional FDC for the ungauged site, and (c) use of a spatial interpolation technique to decode information in FDC to construct streamflow time series for the ungauged site. Area ratio method is conventionally used to transfer streamflow related information from gauged sites to ungauged sites. Attributes that have been considered for the analysis include variables representing hydrology, climatology, topography, land-use/land- cover and soil properties corresponding to catchments in the study area. Effectiveness of the presented methods is assessed using jack knife cross-validation. Conclusions based on the study are presented and discussed. (C) 2015 The Authors. Published by Elsevier B.V.

Reliability Assessment of a Storm Water Drain Network

Most of the cities in India are undergoing rapid development in recent decades, and many rural localities are undergoing transformation to urban hotspots. These developments have associated land use/land cover (LULC) change that effects runoff response from catchments, which is often evident in the form of increase in runoff peaks, volume and velocity in drain network. Often most of the existing storm water drains are in dilapidated stage owing to improper maintenance or inadequate design. The drains are conventionally designed using procedures that are based on some anticipated future conditions. Further, values of parameters/variables associated with design of the network are traditionally considered to be deterministic. However, in reality, the parameters/variables have uncertainty due to natural and/or inherent randomness. There is a need to consider the uncertainties for designing a storm water drain network that can effectively convey the discharge. The present study evaluates performance of an existing storm water drain network in Bangalore, India, through reliability analysis by Advance First Order Second Moment (AFOSM) method. In the reliability analysis, parameters that are considered to be random variables are roughness coefficient, slope and conduit dimensions. Performance of the existing network is evaluated considering three failure modes. The first failure mode occurs when runoff exceeds capacity of the storm water drain network, while the second failure mode occurs when the actual flow velocity in the storm water drain network exceeds the maximum allowable velocity for erosion control, whereas the third failure mode occurs when the minimum flow velocity is less than the minimum allowable velocity for deposition control. In the analysis, runoff generated from subcatchments of the study area and flow velocity in storm water drains are estimated using Storm Water Management Model (SWMM). Results from the study are presented and discussed. The reliability values are low under the three failure modes, indicating a need to redesign several of the conduits to improve their reliability. This study finds use in devising plans for expansion of the Bangalore storm water drain system. (C) 2015 The Authors. Published by Elsevier B.V.

Fast and slow climate responses to CO2 and solar forcing: A linear multivariate regression model characterizing transient climate change

Climate change in response to a change in external forcing can be understood in terms of fast response to the imposed forcing and slow feedback associated with surface temperature change. Previous studies have investigated the characteristics of fast response and slow feedback for different forcing agents. Here we examine to what extent that fast response and slow feedback derived from time-mean results of climate model simulations can be used to infer total climate change. To achieve this goal, we develop a multivariate regression model of climate change, in which the change in a climate variable is represented by a linear combination of its sensitivity to CO2 forcing, solar forcing, and change in global mean surface temperature. We derive the parameters of the regression model using time-mean results from a set of HadCM3L climate model step-forcing simulations, and then use the regression model to emulate HadCM3L-simulated transient climate change. Our results show that the regression model emulates well HadCM3L-simulated temporal evolution and spatial distribution of climate change, including surface temperature, precipitation, runoff, soil moisture, cloudiness, and radiative fluxes under transient CO2 and/or solar forcing scenarios. Our findings suggest that temporal and spatial patterns of total change for the climate variables considered here can be represented well by the sum of fast response and slow feedback. Furthermore, by using a simple 1-D heat-diffusion climate model, we show that the temporal and spatial characteristics of climate change under transient forcing scenarios can be emulated well using information from step-forcing simulations alone.