ALBERT

Description: Global glacier changes: a revised assessment of committed mass losses and sampling uncertainties The Cryosphere, 7, 1565-1577, 2013 Author(s): S. H. Mernild, W. H. Lipscomb, D. B. Bahr, V. Radić, and M. Zemp Most glaciers and ice caps (GIC) are out of balance with the current climate. To return to equilibrium, GIC must thin and retreat, losing additional mass and raising sea level. Because glacier observations are sparse and geographically biased, there is an undersampling problem common to all global assessments. Here, we further develop an assessment approach based on accumulation-area ratios (AAR) to estimate committed mass losses and analyze the undersampling problem. We compiled all available AAR observations for 144 GIC from 1971 to 2010, and found that most glaciers and ice caps are farther from balance than previously believed. Accounting for regional and global undersampling errors, our model suggests that GIC are committed to additional losses of 32 ± 12% of their area and 38 ± 16% of their volume if the future climate resembles the climate of the past decade. These losses imply global mean sea-level rise of 163 ± 69 mm, assuming total glacier volume of 430 mm sea-level equivalent. To reduce the large uncertainties in these projections, more long-term glacier measurements are needed in poorly sampled regions.

Description: Early warning of snow-caused disasters in pastoral areas on the Tibetan Plateau Natural Hazards and Earth System Science, 13, 1411-1425, 2013 Author(s): W. Wang, T. Liang, X. Huang, Q. Feng, H. Xie, X. Liu, M. Chen, and X. Wang This study develops a model for early warning of snow-caused livestock disasters on a county basis and proposes a method of qualitative risk assessment of snow disasters at 500 m resolution for pastoral areas on the Tibetan Plateau (TP). Data used for the model development include remote sensing data, statistical data of weather, livestock, and social economy, and 45 typical snow disaster cases from 2000 to 2010. The principal component analysis (PCA) approach is used to choose 7 crucial factors that contribute over 85% of information for early warning snow disasters on the TP. They are mean annual probability of snow disaster, number of snow-covered days, livestock stocking rate, continual days of mean daily temperature below −10 °C, grassland burial index, rate of snow-covered grassland, and per livestock gross domestic product. The chosen 411 cases from 2008 to 2010 are used to validate the prediction results from the developed early warning model, with an overall accuracy of 85.64% in predicting snow disasters and no disasters. This suggests that the early warning approach developed in the study has operational potential for predicting snow disasters on the TP.

Description: How much does participatory flood management contribute to stakeholders' social capacity building? Empirical findings based on a triangulation of three evaluation approaches Natural Hazards and Earth System Science, 13, 1427-1444, 2013 Author(s): M. Buchecker, S. Menzel, and R. Home Recent literature suggests that dialogic forms of risk communication are more effective to build stakeholders' hazard-related social capacities. In spite of the high theoretical expectations, there is a lack of univocal empirical evidence on the relevance of these effects. This is mainly due to the methodological limitations of the existing evaluation approaches. In our paper we aim at eliciting the contribution of participatory river revitalisation projects on stakeholders' social capacity building by triangulating the findings of three evaluation studies that were based on different approaches: a field-experimental, a qualitative long-term ex-post and a cross-sectional household survey approach. The results revealed that social learning and avoiding the loss of trust were more relevant benefits of participatory flood management than acceptance building. The results suggest that stakeholder involvements should be more explicitly designed as tools for long-term social learning.

Description: Upscaling of evapotranspiration fluxes from instantaneous to daytime scales for thermal remote sensing applications Hydrology and Earth System Sciences Discussions, 10, 7325-7350, 2013 Author(s): C. Cammalleri, M. C. Anderson, and W. P. Kustas Four upscaling methods for estimating daytime evapotranspiration (ET) from single time-of-day snapshots, as commonly retrieved using remote sensing, were compared. These methods are based on the assumption of self-preservation of the ratio between ET and a given reference variable over the daytime hours. The analysis was performed using eddy covariance data collected at 12 AmeriFlux towers, sampling a fairly wide range in climatic and land cover conditions. The choice of energy budget closure method significantly impacted performance using different scaling methodologies. Therefore, a statistical evaluation approach was adopted to better account for the inherent uncertainty in ET fluxes using eddy covariance technique. Overall, this approach suggests that at-surface solar radiation is the most robust reference variable amongst those tested, due to high accuracy of upscaled fluxes and absence of systematic biases. Top-of-atmosphere irradiance was also tested and proved to be reliable under near clear-sky conditions, but tended to overestimate the observed daytime ET during cloudy days. Use of reference ET as a scaling flux did not perform as well as the solar radiation method, but similarly had errors with little seasonal dependency. Finally, the commonly-used evaporative fraction method yielded satisfactory results only in summer months, July and August, and tended to underestimate the observations in the fall/winter seasons from November to January at the flux sites studied.

Description: Errors in climate model daily precipitation and temperature output: time invariance and implications for bias correction Hydrology and Earth System Sciences, 17, 2147-2159, 2013 Author(s): E. P. Maurer, T. Das, and D. R. Cayan When correcting for biases in general circulation model (GCM) output, for example when statistically downscaling for regional and local impacts studies, a common assumption is that the GCM biases can be characterized by comparing model simulations and observations for a historical period. We demonstrate some complications in this assumption, with GCM biases varying between mean and extreme values and for different sets of historical years. Daily precipitation and maximum and minimum temperature from late 20th century simulations by four GCMs over the United States were compared to gridded observations. Using random years from the historical record we select a "base" set and a 10 yr independent "projected" set. We compare differences in biases between these sets at median and extreme percentiles. On average a base set with as few as 4 randomly-selected years is often adequate to characterize the biases in daily GCM precipitation and temperature, at both median and extreme values; 12 yr provided higher confidence that bias correction would be successful. This suggests that some of the GCM bias is time invariant. When characterizing bias with a set of consecutive years, the set must be long enough to accommodate regional low frequency variability, since the bias also exhibits this variability. Newer climate models included in the Intergovernmental Panel on Climate Change fifth assessment will allow extending this study for a longer observational period and to finer scales.

Description: Characterization of groundwater dynamics in landslides in varved clays Hydrology and Earth System Sciences, 17, 2171-2183, 2013 Author(s): J. E. van der Spek, T. A. Bogaard, and M. Bakker Groundwater dynamics may play a significant role in landslides. A detailed model is developed of the groundwater dynamics in landslides in varved clays in the Trièves area in the French Alps. The varved clays consist of a sequence of alternating silt and clay layers, covered by a colluvium layer and intersected by fissures. The hydraulic conductivity of the clay layers is negligible compared to the silt layers. It is conceptualized that fissures form a hydraulic connection between the colluvium and the varved clays. Groundwater recharge flows through the colluvium into the fissures, where water is exchanged horizontally between the fissure and the silt layers of the varved clays. Groundwater flow in the colluvium is simulated with the Boussinesq equation, while flow in the silt layers of the varved clays is simulated with the Richards equation. Longitudinal outflow from the fissure is simulated with a linear-reservoir model. Scattered data of relatively short monitoring periods is available for several landslides in the region. A good similarity between observed and simulated heads is obtained, especially when considering the lack of important physical parameters such as the fissure width and the distance between the monitoring point and the fissure. A simulation for the period 1959–2004 showed some correlation between peaks in the simulated heads and the recorded occurrence of landslides, while the bottom of the varved clays remained saturated during the entire simulation period.

Description: The physical basis for gas transport through polar firn: a case study at Summit, Greenland The Cryosphere Discussions, 7, 2455-2487, 2013 Author(s): A. C. Adolph and M. R. Albert Compared to other natural porous materials, relatively little is known about the physical nature of polar firn. This intricate network of ice and pore space that comprises the top 60–100 m of the polar ice sheets is the framework that forms the natural archive of past climate information. Despite the many implications for ice core interpretation, direct measurements of physical properties throughout the firn column are limited. Models of gas transport through firn are used to interpret in-situ chemical data which is retrieved to analyze past atmospheric composition. These traditional models treat the firn as a "black box," with gas transport parameters tuned to match gas concentrations with depth to known atmospheric histories. Though this method has been largely successful and provided very useful insights, there are still many questions and uncertainties to be addressed. This work seeks to understand the impact of firn structure on gas transport in firn from a first principles standpoint through direct measurements of permeability, gas diffusivity and microstructure. The relationships between gas transport properties and microstructure will be characterized and compared to existing relationships for general porous media. Direct measurements of gas diffusivity are compared to diffusivities deduced from models based on firn air chemical sampling. Our comparison illuminates the primary importance of including microstructural parameters, beyond just porosity or density, in mass transport modeling, and it provides insights about the nature of gas transport throughout the firn column. Guidance is provided for development of next-generation firn air transport models.

Description: Assessing the hydrological effect of the check dams in the Loess Plateau, China, by model simulations Hydrology and Earth System Sciences, 17, 2185-2193, 2013 Author(s): Y. D. Xu, B. J. Fu, and C. S. He Check dams are commonly used for soil conservation. In the Loess Plateau of China, check dams have been widely constructed as the principal means to retain floodwater and intercept soil sediments since the 1970s. For instance, there are more than 6572 check dams in the Yanhe watershed with an area of 7725 km 2 in the Loess Plateau. However, little research has been done to quantify the hydrological effects of the check dams. In this research, the SWAT model (Soil and Water Assessment Tool) was applied to simulate the runoff and sediment in the Yanhe watershed. We treated the 1950s to 1960s as the reference period since there were very few check dams during the period. The model was firstly calibrated and validated in the reference period. The calibrated model was then used in the later periods to simulate the hydrological effects of the check dams. The results showed that the check dams had a regulation effect on runoff and a retention effect on sediment. From 1984 to 1987, the runoff in rainy season (from May to October) decreased by 1.54 m 3 s −1 (14.7%) to 3.13 m 3 s −1 (25.9%) due to the check dams; while in dry season (from November to the following April), runoff increased by 1.46 m 3 s −1 (60.5%) to 1.95 m 3 s −1 (101.2%); the sediment in rainy season decreased by 2.49 × 10 6 ton (34.6%) to 4.35 × 10 6 ton (48.0%). From 2006 to 2008, the runoff in rainy season decreased by 0.79 m 3 s −1 (15.5%) to 1.75 m 3 s −1 (28.9%), and the runoff in dry season increased by 0.51 m 3 s −1 (20.1%) to 0.97 m 3 s −1 (46.4%); the sediment in rainy season decreased by 2.03 × 10 6 ton (79.4%) to 3.12 × 10 6 ton (85.5%). Construction of the large number of check dams in the Loess Plateau has enhanced the region's capacity to control the runoff and sediment. In the Yanhe watershed, the annual runoff was reduced by less than 14.3% due to the check dams; and the sediment in rainy season was blocked by up to 85.5%. Thus, check dams are effective measures for soil erosion control in the Loess Plateau.

Description: Global isoscapes for δ 18 O and δ 2 H in precipitation: improved prediction using regionalized climatic regression models Hydrology and Earth System Sciences Discussions, 10, 7351-7393, 2013 Author(s): S. Terzer, L. I. Wassenaar, L. J. Araguás-Araguás, and P. K. Aggarwal A Regionalized Climatic Water Isotope Prediction (RCWIP) approach, based on the Global Network for Isotopes in Precipitation (GNIP), was demonstrated for the purposes of predicting point- and large-scale spatiotemporal patterns of the stable isotope compositions of water (δ 2 H, δ 18 O) in precipitation around the world. Unlike earlier global domain and fixed regressor models, RCWIP pre-defined thirty-six climatic cluster domains, and tested all model combinations from an array of climatic and spatial regressor variables to obtain the best predictive approach to each cluster domain, as indicated by RMSE and variogram analysis. Fuzzy membership fractions were thereafter used as the weights to seamlessly amalgamate results of the optimized climatic zone prediction models into a single predictive mapping product, such as global or regional amount-weighted mean annual, mean monthly or growing-season δ 18 O/δ 2 H in precipitation. Comparative tests revealed the RCWIP approach outperformed classical global-fixed regression-interpolation based models more than 67% of the time, and significantly improved upon predictive accuracy and precision. All RCWIP isotope mapping products are available as gridded GeoTIFF files from the IAEA website ( www.iaea.org/water ) and are for use in hydrology, climatology, food authenticity, ecology, and forensics.

Description: Stakeholders' issues for action during the warning process and the interpretation of forecasts' uncertainties Natural Hazards and Earth System Science, 13, 1469-1479, 2013 Author(s): L. Créton-Cazanave and C. Lutoff This article considers the socio-geographical approach carried out as part of the MedUp program. It presents a study of the ways that the French "actors" manage forecast uncertainties during a flash flood warning process. In order to better understand the role of forecasts' uncertainties in decision making, we focused on the actions people took and how what they say explains their actions. The practices of actors involved in warnings for the Vidourle watershed (Gard, France), in particular, are analyzed using a practice-based approach. A set of categories of the "actors" was developed based on their descriptions of the problems they faced during the flash flood warning, independent of their socio-professional status and position in the warning chain. Five actor profiles result from this: Translators, Managers, Committed, Navigators and Vulnerable. For each profile, specific action contexts are defined, determining how each deals with uncertainty.

Description: Augmenting Austrian flood management practices through geospatial predictive analytics: a study in Carinthia Natural Hazards and Earth System Science, 13, 1445-1455, 2013 Author(s): S. M. Ward and G. Paulus The Danube River basin has long been the location of significant flooding problems across central Europe. The last decade has seen a sharp increase in the frequency, duration and intensity of these flood events, unveiling a dire need for enhanced flood management policy and tools in the region. Located in the southern portion of Austria, the state of Carinthia has experienced a significant volume of intense flood impacts over the last decade. Although the Austrian government has acknowledged these issues, their remedial actions have been primarily structural to date. Continued focus on controlling the natural environment through infrastructure while disregarding the need to consider alternative forms of assessing flood exposure will only act as a provisional solution to this inescapable risk. In an attempt to remedy this flaw, this paper highlights the application of geospatial predictive analytics and spatial recovery index as a proxy for community resilience, as well as the cultural challenges associated with the application of foreign models within an Austrian environment.

Description: Water displacement by sewer infrastructure in the Grote Nete catchment, Belgium, and its hydrological regime effects Hydrology and Earth System Sciences Discussions, 10, 7425-7467, 2013 Author(s): D. Vrebos, T. Vansteenkiste, J. Staes, P. Willems, and P. Meire Urbanization and especially impervious areas, in combination with wastewater treatment infrastructure, can exert several pressures on the hydrological cycle. These pressures were studied for the Grote Nete catchment in Belgium (8.18% impervious area and 3.89% effective impervious area), based on a combination of empirical and model-based approaches. The effective impervious area, combined with the extent of the wastewater collection regions which do not coincide with the natural catchment boundaries, was used as an indicator for the urbanization pressure. Our study revealed changes in the total upstream areas of the subcatchments between −16% and +3%, and in upstream impervious areas between −99% and +64%. These changes lead to important inter-catchment water transfers. Based on simulations with a physically-based and spatially-distributed hydrological catchment model, profound impacts of effective impervious area on infiltration and runoff were found. The model results show that the changes in impervious areas and related water displacements in and between catchments due to the installation of the wastewater treatment infrastructure severely impacted low flows, peak flows and seasonal trends. They moreover show that it is difficult, but of utmost importance, to incorporate these pressures and artificial processes in an accurate way during the development of hydrological models for urbanized catchments.

Description: Recent changes in spring snowmelt timing in the Yukon River basin detected by passive microwave satellite data The Cryosphere, 7, 905-916, 2013 Author(s): K. A. Semmens and J. M. Ramage Spring melt is a significant feature of high latitude snowmelt dominated drainage basins influencing hydrological and ecological processes such as snowmelt runoff and green-up. Melt duration, defined as the transition period from snowmelt onset until the end of the melt refreeze, is characterized by high diurnal amplitude variations (DAV) where the snowpack is melting during the day and refreezing at night, after which the snowpack melts constantly until depletion. Determining trends for this critical period is necessary for understanding how the Arctic is changing with rising temperatures and provides a baseline from which to assess future change. To study this dynamic period, brightness temperature ( T b ) data from the Special Sensor Microwave Imager (SSM/I) 37 V-GHz frequency from 1988 to 2010 were used to assess snowmelt timing trends for the Yukon River basin, Alaska/Canada. Annual T b and DAV for 1434 Equal-Area Scalable Earth (EASE)-Grid pixels (25 km resolution) were processed to determine melt onset and melt refreeze dates from T b and DAV thresholds previously established in the region. Temporal and spatial trends in the timing of melt onset and melt refreeze, and the duration of melt were analyzed for the 13 sub-basins of the Yukon River basin with three different time interval approaches. Results show a lengthening of the melt period for the majority of the sub-basins with a significant trend toward later end of melt refreeze after which the snowpack melts day and night leading to snow clearance, peak discharge, and green-up. Earlier melt onset trends were also found in the higher elevations and northernmost sub-basins (Porcupine, Chandalar, and Koyukuk rivers). Latitude and elevation displayed the dominant controls on melt timing variability and spring solar flux was highly correlated with melt timing in middle (∼600–1600 m) elevations.

Description: Solving Richards Equation for snow improves snowpack meltwater runoff estimations The Cryosphere Discussions, 7, 2373-2412, 2013 Author(s): N. Wever, C. Fierz, C. Mitterer, H. Hirashima, and M. Lehning The runoff from the snow cover during spring snow melt or rain-on-snow events is an important factor in the hydrological cycle. In this study, water transport schemes for a 1-dimensional physical based snowpack model are compared to 14 yr of lysimeter measurements at a high alpine site. The schemes include a simple bucket-type approach, an approximation of Richards Equation (RE), and the full RE. The results show that daily sums of runoff are strongly related to a positive energy balance of the snow cover and therefore, all water transport schemes show very similar performance in terms of Nash–Sutcliffe efficiency (NSE) coefficients (around 0.59) and r 2 values (around 0.77). Timing of the arrival of meltwater in spring at the bottom of the snowpack showed differences between the schemes, where especially in the bucket-type and approximated RE approach, meltwater release is slower than in the measurements. Overall, solving RE for the snow cover yields the best agreement between modelled and measured runoff. On sub-daily time scales, the water transport schemes behave very differently. Also here, solving RE provides the highest agreement between modelled and measured runoff in terms of NSE coefficient (0.48), where other water transport schemes loose any predictive power. This appears to be mainly due to bad timing of meltwater release during the day. Accordingly, solving RE for the snow cover improves several aspects of modelling snow cover runoff. The additional computational cost was found to be in the order of a factor of 1.5.

Description: An upper-bound estimate for the accuracy of volume-area scaling The Cryosphere Discussions, 7, 2293-2331, 2013 Author(s): D. Farinotti and M. Huss Volume-area scaling is the most popular method for estimating the ice volume of large glacier samples. Here, a series of resampling experiments based on different sets of synthetic data are presented in order to derive an upper-bound estimate (i.e. a level achieved only with ideal conditions) for the accuracy of its application. We also quantify the maximum accuracy expected when scaling is used for determining the glacier volume change, and area change of a given glacier population. A comprehensive set of measured glacier areas, volumes, area and volume changes is evaluated to investigate the impact of real-world data quality on the so assessed accuracies. For populations larger than a few thousand glaciers, the total ice volume can be recovered within 30% if all measurements available worldwide are used for estimating the scaling coefficients. Assuming no systematic biases in ice volume measurements, their uncertainty is of secondary importance. Knowing the individual areas of a glacier sample for two points in time allows recovering the corresponding ice volume change within 40% for populations larger than a few hundred glaciers, both for steady-state and transient geometries. If ice volume changes can be estimated without bias, glacier area changes derived from volume-area scaling show similar uncertainties as for the volume changes. This paper does not aim at making a final judgement about the suitability of volume-area scaling, but provides the means for assessing the accuracy expected from its application.

Description: Evaluation of the snow regime in dynamic vegetation land surface models using field measurements The Cryosphere Discussions, 7, 2333-2372, 2013 Author(s): E. Kantzas, M. Lomas, S. Quegan, and E. Zakharova An increasing number of studies have demonstrated the significant climatic and ecological changes occurring in the northern latitudes over the past decades. As coupled, earth-system models attempt to describe and simulate the dynamics and complex feedbacks of the Arctic environment, it is important to reduce their uncertainties in short-term predictions by improving the description of both the systems processes and its initial state. This study focuses on snow-related variables and extensively utilizes a historical data set (1966–1996) of field snow measurements acquired across the extend of the Former Soviet Union (FSU) to evaluate a range of simulated snow metrics produced by a variety of land surface models, most of them embedded in IPCC-standard climate models. We reveal model-specific issues in simulating snow dynamics such as magnitude and timings of SWE as well as evolution of snow density. We further employ the field snow measurements alongside novel and model-independent methodologies to extract for the first time (i) a fresh snow density value (57–117 kg m –3 ) for the region and (ii) mean monthly snowpack sublimation estimates across a grassland-dominated western (November–February) [9.2, 6.1, 9.15, 15.25] mm and forested eastern sub-sector (November–March) [1.53, 1.52, 3.05, 3.80, 12.20] mm; we subsequently use the retrieved values to assess relevant model outputs. The discussion session consists of two parts. The first describes a sensitivity study where field data of snow depth and snow density are forced directly into the surface heat exchange formulation of a land surface model to evaluate how inaccuracies in simulating snow metrics affect important modeled variables and carbon fluxes such as soil temperature, thaw depth and soil carbon decomposition. The second part showcases how the field data can be assimilated with ready-available optimization techniques to pinpoint model issues and improve their performance.

Description: Virtual water trade and development in Africa Hydrology and Earth System Sciences Discussions, 10, 7291-7324, 2013 Author(s): M. Konar and K. Caylor A debate has long existed on the relationships between human population, natural resources, and development. Recent research has expanded this debate to include the impacts of trade; specifically, virtual water trade, or the water footprint of traded commodities. We conduct an empirical analysis of the relationships between virtual water trade, population, and development in Africa. We find that increases in virtual water imports do not lead to increases in population growth nor do they diminish human welfare. We establish a new index of virtual water trade openness and show that levels of undernourishment tend to fall with increased values of virtual water trade openness. Countries with small dam storage capacity obtain a higher fraction of their agricultural water requirements from external sources, which may indicate implicit "infrastructure sharing" across nations. Globally, increased crop exports tends to correlate with increased crop water use efficiency, though this relationship does not hold for Africa. However, internal African trade is much more efficient in terms of embodied water resources than any other region in the world. Thus, internal African trade patterns may be compensating for poor internal production systems.

Description: Seasonal controls on snow distribution and aerial ablation at the snow-patch and landscape scales, McMurdo Dry Valleys, Antarctica The Cryosphere, 7, 917-931, 2013 Author(s): J. W. Eveland, M. N. Gooseff, D. J. Lampkin, J. E. Barrett, and C. D. Takacs-Vesbach Accumulated snow in the McMurdo Dry Valleys, while limited, has great ecological significance to subnivian soil environments. Though sublimation dominates the ablation process in this region, measurable increases in soil moisture and insulation from temperature extremes provide more favorable conditions with respect to subnivian soil communities. While precipitation is not substantial, significant amounts of snow can accumulate, via wind transport, in topographic lees along the valley bottoms, forming thousands of discontinuous snow patches. These patches have the potential to act as significant sources of local meltwater, controlling biogeochemical cycling and the landscape distribution of microbial communities. Therefore, determining the spatial and temporal dynamics of snow at multiple scales is imperative to understanding the broader ecological role of snow in this region. High-resolution satellite imagery acquired during the 2009–2010 and 2010–2011 austral summers was used to quantify the distribution of snow across Taylor and Wright valleys. Extracted snow-covered area from the imagery was used as the basis for assessing inter-annual variability and seasonal controls on accumulation and ablation of snow at multiple scales. In addition to landscape analyses, fifteen 1 km 2 plots (3 in each of 5 study regions) were selected to assess the prevalence of snow cover at finer spatial scales, referred to herein as the snow-patch scale. Results confirm that snow patches tend to form in the same locations each year with some minor deviations observed. At the snow-patch scale, neighboring patches often exhibit considerable differences in aerial ablation rates, and particular snow patches do not reflect trends for snow-covered area observed at the landscape scale. These differences are presumably related to microtopographic influences acting on individual snow patches, such as wind sheltering and differences in snow depth due to the underlying topography. This highlights the importance of both the landscape and snow-patch scales in assessing the effects of snow cover on biogeochemical cycling and microbial communities.

Description: An interdisciplinary swat ecohydrological model to define catchment-scale hydrologic partitioning Hydrology and Earth System Sciences Discussions, 10, 7235-7290, 2013 Author(s): C. L. Shope, G. R. Maharjan, J. Tenhunen, B. Seo, K. Kim, J. Riley, S. Arnhold, T. Koellner, Y. S. Ok, S. Peiffer, B. Kim, J.-H. Park, and B. Huwe Land use and climate change have long been implicated in modifying ecosystem services, such as water quality and water yield, biodiversity, and agricultural production. To account for future effects on ecosystem services, the integration of physical, biological, economic, and social data over several scales must be implemented to assess the effects on natural resource availability and use. Our objective is to assess the capability of the SWAT model to capture short-duration monsoonal rainfall-runoff processes in complex mountainous terrain under rapid, event-driven processes in a monsoonal environment. To accomplish this, we developed a unique quality-control gap-filling algorithm for interpolation of high frequency meteorological data. We used a novel multi-location, multi-optimization calibration technique to improve estimations of catchment-wide hydrologic partitioning. We calibrated the interdisciplinary model to a combination of statistical, hydrologic, and plant growth metrics. In addition, we used multiple locations of different drainage area, aspect, elevation, and geologic substrata distributed throughout the catchment. Results indicate scale-dependent sensitivity of hydrologic partitioning and substantial influence of engineered features. While our model accurately reproduced observed discharge variability, the addition of hydrologic and plant growth objective functions identified the importance of culverts in catchment-wide flow distribution. The results of this study provide a valuable resource to describe landscape controls and their implication on discharge, sediment transport, and nutrient loading. This study also shows the challenges of applying the SWAT model to complex terrain and extreme environments. By incorporating anthropogenic features into modeling scenarios, we can greatly enhance our understanding of the hydroecological impacts on ecosystem services.

Description: Determining spatial variability of dry spells: a Markov-based method, applied to the Makanya catchment, Tanzania Hydrology and Earth System Sciences, 17, 2161-2170, 2013 Author(s): B. M. C. Fischer, M. L. Mul, and H. H. G. Savenije With a growing world population and a trend towards more resource-intensive diets, pressure on land and water resources for food production will continue to increase in the coming decades. Large parts of the world rely on rainfed agriculture for their food security. In Africa, 90% of the food production is from rainfed agriculture, generally with low yields and a high risk of crop failure. One of the main reasons for crop failure is the occurrence of dry spells during the growing season. Key indicators are the critical dry spell duration and the probability of dry spell occurrence. In this paper a new Markov-based framework is presented to spatially map the length of dry spells for fixed probabilities of non-exceedance. The framework makes use of spatially varying Markov coefficients that are correlated to readily available spatial information such as elevation and distance to the sea. The dry spell map thus obtained is compared to the spatially variable critical dry spell duration, based on soil properties and crop water requirements, to assess the probability of crop failure in different locations. The results show that in the Makanya catchment the length of dry spell occurrence is highly variable in space, even over relatively short distances. In certain areas the probability of crop failure reaches levels that make rainfed agricultural unsustainable, even close to areas where currently rainfed agriculture is successfully being practised. This method can be used to identify regions that are vulnerable to dry spells and, subsequently, to develop strategies for supplementary irrigation or rainwater harvesting.

Description: Spatial debris-cover effect on the maritime glaciers of Mount Gongga, south-eastern Tibetan Plateau The Cryosphere Discussions, 7, 2413-2453, 2013 Author(s): Y. Zhang, Y. Hirabayashi, K. Fujita, S. Liu, and Q. Liu The Tibetan Plateau and surroundings contain a large number of debris-covered glaciers, on which debris cover affects glacier response to climate change by altering ice melting rates and spatial patterns of mass loss. Insufficient spatial distribution of debris thickness data makes it difficult to analyze regional debris-cover effects. Mount Gongga glaciers, maritime glaciers in the south-eastern Tibetan Plateau, are characterized by a substantial reduction in glacier length and ice mass in recent decades. Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER)-derived thermal property of the debris layer reveals that 68% of the glaciers have extensive mantles of supraglacial debris in their ablation zones, in which the proportion of debris cover to total glacier area varies from 1.74% to 53.0%. Using a surface energy-mass balance model accounting for the debris-cover effect applied at a regional scale, we find that although the presence of supraglacial debris has a significant insulating effect on heavily debris-covered glaciers, it accelerates ice melting on ~ 10.2% of the total ablation area and produces rapid wastage of ~ 25% of the debris-covered glaciers, resulting in the similar mass losses between debris-covered and debris-free glaciers. Widespread debris cover also facilitates the development of active terminus regions. Regional differences in the debris-cover effect are apparent, highlighting the importance of debris cover for understanding glacier status and hydrology in both the Tibetan Plateau and other mountain ranges around the world.

Description: Decadal changes from a multi-temporal glacier inventory of Svalbard The Cryosphere Discussions, 7, 2489-2532, 2013 Author(s): C. Nuth, J. Kohler, M. König, A. von Deschwanden, J. O. Hagen, A. Kääb, G. Moholdt, and R. Pettersson We present a multi-temporal digital inventory of Svalbard glaciers with the most recent from the late 2000s containing 33 775 km 2 of glaciers, or 57% of the total land area of the archipelago. At present, 68% of the glaciated area of Svalbard drains through tidewater glaciers that have a summed terminus width of ~ 740 km. The glaciated area over the entire archipelago has decreased by an average of 80 km 2 a −1 over the past ~ 30 yr, representing a reduction of 7%. For a sample of ~ 400 glaciers (10 000 km 2 ) in the south and west of Spitsbergen, three digital inventories are available from 1930/60s, 1990 and 2007 from which we calculate average changes during 2 epochs. In the more recent epoch, the terminus retreat was larger than in the earlier epoch while area shrinkage was smaller. The contrasting pattern may be explained by the decreased lateral wastage of the glacier tongues. Temporal retreat rates for individual glaciers show a mix of accelerating and decelerating trends, reflecting the large spatial variability of glacier types and climatic/dynamic response times in Svalbard. Last, retreat rates estimated by dividing glacier area changes by the tongue width are larger than centerline retreat due to a more encompassing frontal change estimate with inclusion of lateral area loss.

Description: Technical note: Method of Morris effectively reduces the computational demands of global sensitivity analysis for distributed watershed models Hydrology and Earth System Sciences Discussions, 10, 4275-4299, 2013 Author(s): J. D. Herman, J. B. Kollat, P. M. Reed, and T. Wagener The increase in spatially distributed hydrologic modeling warrants a corresponding increase in diagnostic methods capable of analyzing complex models with large numbers of parameters. Sobol ' sensitivity analysis has proven to be a valuable tool for diagnostic analyses of hydrologic models. However, for many spatially distributed models, the Sobol ' method requires a prohibitive number of model evaluations to reliably decompose output variance across the full set of parameters. We investigate the potential of the method of Morris, a screening-based sensitivity approach, to provide results sufficiently similar to those of the Sobol ' method at a greatly reduced computational expense. The methods are benchmarked on the Hydrology Laboratory Research Distributed Hydrologic Model (HL-RDHM) model over a six-month period in the Blue River Watershed, Oklahoma, USA. The Sobol ' method required over six million model evaluations to ensure reliable sensitivity indices, corresponding to more than 30 000 computing hours and roughly 180 gigabytes of storage space. We find that the method of Morris is able to correctly identify sensitive and insensitive parameters with 300 times fewer model evaluations, requiring only 100 computing hours and 1 gigabyte of storage space. Method of Morris proves to be a promising diagnostic approach for global sensitivity analysis of highly parameterized, spatially distributed hydrologic models.

Description: Estimating actual, potential, reference crop and pan evaporation using standard meteorological data: a pragmatic synthesis Hydrology and Earth System Sciences, 17, 1331-1363, 2013 Author(s): T. A. McMahon, M. C. Peel, L. Lowe, R. Srikanthan, and T. R. McVicar This guide to estimating daily and monthly actual, potential, reference crop and pan evaporation covers topics that are of interest to researchers, consulting hydrologists and practicing engineers. Topics include estimating actual evaporation from deep lakes and from farm dams and for catchment water balance studies, estimating potential evaporation as input to rainfall-runoff models, and reference crop evapotranspiration for small irrigation areas, and for irrigation within large irrigation districts. Inspiration for this guide arose in response to the authors' experiences in reviewing research papers and consulting reports where estimation of the actual evaporation component in catchment and water balance studies was often inadequately handled. Practical guides using consistent terminology that cover both theory and practice are not readily available. Here we provide such a guide, which is divided into three parts. The first part provides background theory and an outline of the conceptual models of potential evaporation of Penman, Penman–Monteith and Priestley–Taylor, as well as discussions of reference crop evapotranspiration and Class-A pan evaporation. The last two sub-sections in this first part include techniques to estimate actual evaporation from (i) open-surface water and (ii) landscapes and catchments (Morton and the advection-aridity models). The second part addresses topics confronting a practicing hydrologist, e.g. estimating actual evaporation for deep lakes, shallow lakes and farm dams, lakes covered with vegetation, catchments, irrigation areas and bare soil. The third part addresses six related issues: (i) automatic (hard wired) calculation of evaporation estimates in commercial weather stations, (ii) evaporation estimates without wind data, (iii) at-site meteorological data, (iv) dealing with evaporation in a climate change environment, (v) 24 h versus day-light hour estimation of meteorological variables, and (vi) uncertainty in evaporation estimates. This paper is supported by a Supplement that includes 21 sections enhancing the material in the text, worked examples of many procedures discussed in the paper, a program listing (Fortran 90) of Morton's WREVAP evaporation models along with tables of monthly Class-A pan coefficients for 68 locations across Australia and other information.

Description: The effect of black carbon on reflectance of snow in the accumulation area of glaciers in the Baspa basin, Himachal Pradesh, India The Cryosphere Discussions, 7, 1359-1382, 2013 Author(s): A. V. Kulkarni, G. Vinay Kumar, H. S. Negi, J. Srinivasan, and S. K. Satheesh Himalayan glaciers are being extensively debated in scientific and public forums, as changes in their distribution can significantly affect the availability of water in many rivers originating in the region. The distribution of glaciers can be influenced by mass balance, and most of the glaciers located in the Pir Panjal and Greater Himalayan mountain ranges are losing mass at the rate of almost a meter per year. The Equilibrium Line Altitude (ELA) has also shifted upward by 400 m in the last two decades. This upward migration of ELA and the loss in mass could have been influenced by changes in temperature, precipitation and by the deposition of black carbon in the accumulation area of glaciers. The deposition of black carbon can reduce the albedo of snow in the accumulation area leading to faster melting of snow and causing more negative mass balance. In this investigation, a change in reflectance in the accumulation area of the Baspa basin is analysed for the year 2009, as the region has experienced extensive forest fires along with northern Indian biomass burning. The investigation has shown that: (1) The number of forest fires in the summer of 2009 was substantially higher than in any other year between 2001 and 2010; (2) the drop in reflectance in the visible region from April to May in the accumulation area was significantly higher in the year 2009 than in any other year from 2000 to 2012; (3) the temperature of the region was substantially lower than the freezing point during the active fire period of 2009, indicating the small influence of liquid water and grain size; (4) the drop in reflectance was observed only in the visible region, indicating role of contamination; (5) in the visible region, a mean drop in reflectance of 21± 5% was observed during the active fire period in the accumulation area. At some places, the drop was as high as 50 ± 5%. This can only be explained by the deposition of black carbon. The study suggests that a change in snow albedo in the accumulation area due to the deposition of black carbon from anthropogenic and natural causes can influence the mass balance of the glaciers in the Baspa basin, Himachal Pradesh, India.

Description: Socio-hydrology: conceptualising human-flood interactions Hydrology and Earth System Sciences Discussions, 10, 4515-4536, 2013 Author(s): G. Di Baldassarre, A. Viglione, G. Carr, L. Kuil, J. L. Salinas, and G. Blöschl Over history, humankind has tended to settle near streams because of the role of rivers as transportation corridors and the fertility of riparian areas. However, human settlements in floodplains have been threatened by the risk of flooding. Possible responses have been to resettle away and/or modify the river system by building flood control structures. This has led to a complex web of interactions and feedback mechanisms between hydrological and social processes in settled floodplains. This paper is an attempt to conceptualise these interplays for hypothetical human-flood systems. We develop a simple, dynamic model to represent the interactions and feedback loops between hydrological and social processes. The model is then used to explore the dynamics of the human-flood system and the effect of changing individual characteristics, including external forcing such as technological development. The results show that the conceptual model is able to reproduce reciprocal effects between floods and people as well as the emergence of typical patterns. For instance, when levees are built or raised to protect floodplain areas, their presence not only reduces the frequency of flooding, but also exacerbates high water levels. Then, because of this exacerbation, higher flood protection levels are required by the society. As a result, more and more flooding events are avoided, but rare and catastrophic events take place.

Description: CDRD and PNPR satellite passive microwave precipitation retrieval algorithms: EuroTRMM/EURAINSAT origins and H-SAF operations Natural Hazards and Earth System Science, 13, 887-912, 2013 Author(s): A. Mugnai, E. A. Smith, G. J. Tripoli, B. Bizzarri, D. Casella, S. Dietrich, F. Di Paola, G. Panegrossi, and P. Sanò Satellite Application Facility on Support to Operational Hydrology and Water Management (H-SAF) is a EUMETSAT (European Organisation for the Exploitation of Meteorological Satellites) program, designed to deliver satellite products of hydrological interest (precipitation, soil moisture and snow parameters) over the European and Mediterranean region to research and operations users worldwide. Six satellite precipitation algorithms and concomitant precipitation products are the responsibility of various agencies in Italy. Two of these algorithms have been designed for maximum accuracy by restricting their inputs to measurements from conical and cross-track scanning passive microwave (PMW) radiometers mounted on various low Earth orbiting satellites. They have been developed at the Italian National Research Council/Institute of Atmospheric Sciences and Climate in Rome (CNR/ISAC-Rome), and are providing operational retrievals of surface rain rate and its phase properties. Each of these algorithms is physically based, however, the first of these, referred to as the Cloud Dynamics and Radiation Database (CDRD) algorithm, uses a Bayesian-based solution solver, while the second, referred to as the PMW Neural-net Precipitation Retrieval (PNPR) algorithm, uses a neural network-based solution solver. Herein we first provide an overview of the two initial EU research and applications programs that motivated their initial development, EuroTRMM and EURAINSAT (European Satellite Rainfall Analysis and Monitoring at the Geostationary Scale), and the current H-SAF program that provides the framework for their operational use and continued development. We stress the relevance of the CDRD and PNPR algorithms and their precipitation products in helping secure the goals of H-SAF's scientific and operations agenda, the former helpful as a secondary calibration reference to other algorithms in H-SAF's complete mix of algorithms. Descriptions of the algorithms' designs are provided including a few examples of their performance. This aspect of the development of the two algorithms is placed in the context of what we refer to as the TRMM era, which is the era denoting the active and ongoing period of the Tropical Rainfall Measuring Mission (TRMM) that helped inspire their original development. In 2015, the ISAC-Rome precipitation algorithms will undergo a transformation beginning with the upcoming Global Precipitation Measurement (GPM) mission, particularly the GPM Core Satellite technologies. A few years afterward, the first pair of imaging and sounding Meteosat Third Generation (MTG) satellites will be launched, providing additional technological advances. Various of the opportunities presented by the GPM Core and MTG satellites for improving the current CDRD and PNPR precipitation retrieval algorithms, as well as extending their product capability, are discussed.

Description: Integration of remote sensing, RUSLE and GIS to model potential soil loss and sediment yield (SY) Hydrology and Earth System Sciences Discussions, 10, 4567-4596, 2013 Author(s): H. Kamaludin, T. Lihan, Z. Ali Rahman, M. A. Mustapha, W. M. R. Idris, and S. A. Rahim Land use activities within a basin serve as one of the contributing factors which cause deterioration of river water quality through its potential effect on erosion. Sediment yield in the form of suspended solid in the river water body which is transported to the coastal area occurs as a sign of lowering of the water quality. Hence, the aim of this study was to determine potential soil loss using the Revised Universal Soil Loss Equation (RUSLE) model and the sediment yield, in the Geographical Information Systems (GIS) environment within selected sub-catchments of Pahang River Basin. RUSLE was used to estimate potential soil losses and sediment yield by utilizing information on rainfall erosivity ( R ) using interpolation of rainfall data, soil erodibility ( K ) using field measurement and soil map, vegetation cover ( C ) using satellite images, topography (LS) using DEM and conservation practices ( P ) using satellite images. The results indicated that the rate of potential soil loss in these sub-catchments ranged from very low to extremely high. The area covered by very low to low potential soil loss was about 99%, whereas moderate to extremely high soil loss potential covered only about 1% of the study area. Sediment yield represented only 1% of the potential soil loss. The sediment yield (SY) value in Pahang River turned out to be higher closer to the river mouth because of the topographic character, climate, vegetation type and density, and land use within the drainage basin.

Description: Probability distributions for explaining hydrological losses in South Australian catchments Hydrology and Earth System Sciences Discussions, 10, 4597-4626, 2013 Author(s): S. H. P. W. Gamage, G. A. Hewa, and S. Beecham The wide variability of hydrological losses in catchments is due to multiple variables that affect the rainfall-runoff process. Accurate estimation of hydrological losses is required for making vital decisions in design applications that are based on design rainfall models and rainfall-runoff models. Using representative single values of losses, despite their wide variability, is common practice, especially in Australian studies. This practice leads to issues such as over or under estimation of design floods. Probability distributions can be used as a better representation of losses. In particular, using joint probability approaches (JPA), probability distributions can be incorporated into hydrological loss parameters in design models. However, lack of understanding of loss distributions limits the benefit of using JPA. The aim of this paper is to identify a probability distribution function that can successfully describe hydrological losses in South Australian (SA) catchments. This paper describes suitable parametric and non-parametric distributions that can successfully describe observed loss data. The goodness-of-fit of the fitted distributions and quantification of the errors associated with quantile estimation are also discussed a two-parameter Gamma distribution was identified as one that successfully described initial loss (IL) data of the selected catchments. Also, a non-parametric standardised distribution of losses that describes both IL and continuing loss (CL) data were identified. The results obtained for the non-parametric methods were compared with similar studies carried out in other parts of Australia and a remarkable degree of consistency was observed. The results will be helpful in improving design flood applications.

Description: Ice-shelf buttressing and the stability of marine ice sheets The Cryosphere, 7, 647-655, 2013 Author(s): G. H. Gudmundsson Ice-shelf buttressing and the stability of marine-type ice sheets are investigated numerically. Buttressing effects are analysed for a situation where a stable grounding line is located on a bed sloping upwards in the direction of flow. Such grounding-line positions are known to be unconditionally unstable in the absence of transverse flow variations. It is shown that ice-shelf buttressing can restore stability under these conditions. Ice flux at the grounding line is, in general, not a monotonically increasing function of ice thickness. This, possibly at first somewhat counterintuitive result, is found to be fully consistent with recent theoretical work. Grounding lines on retrograde slopes are conditionally stable, and the stability regime is a non-trivial function of bed and ice-shelf geometry. The stability of grounding lines cannot be assessed from considerations of local bed slope only.

Description: Evidence and analysis of 2012 Greenland records from spaceborne observations, a regional climate model and reanalysis data The Cryosphere, 7, 615-630, 2013 Author(s): M. Tedesco, X. Fettweis, T. Mote, J. Wahr, P. Alexander, J. E. Box, and B. Wouters A combined analysis of remote sensing observations, regional climate model (RCM) outputs and reanalysis data over the Greenland ice sheet provides evidence that multiple records were set during summer 2012. Melt extent was the largest in the satellite era (extending up to ∼97% of the ice sheet) and melting lasted up to ∼2 months longer than the 1979–2011 mean. Model results indicate that near surface temperature was ∼3 standard deviations (σ) above the 1958–2011 mean, while surface mass balance (SMB) was ∼3σ below the mean and runoff was 3.9σ above the mean over the same period. Albedo, exposure of bare ice and surface mass balance also set new records, as did the total mass balance with summer and annual mass changes of, respectively, −627 Gt and −574 Gt, 2σ below the 2003–2012 mean. We identify persistent anticyclonic conditions over Greenland associated with anomalies in the North Atlantic Oscillation (NAO), changes in surface conditions (e.g., albedo, surface temperature) and preconditioning of surface properties from recent extreme melting as major driving mechanisms for the 2012 records. Less positive if not increasingly negative SMB will likely occur should these characteristics persist.

Description: Corrigendum to "Calibration of aerodynamic roughness over the Tibetan Plateau with Ensemble Kalman Filter analysed heat flux" published in Hydrol. Earth Syst. Sci., 16, 4291–4302, 2012 Hydrology and Earth System Sciences, 17, 1309-1310, 2013 Author(s): J. H. Lee, J. Timmermans, Z. Su, and M. Mancini No abstract available.

Description: Technical Note: Using wavelet analyses on water depth time series to detect glacial influence in high-mountain hydrosystems Hydrology and Earth System Sciences Discussions, 10, 4369-4395, 2013 Author(s): S. Cauvy-Fraunié, T. Condom, A. Rabatel, M. Villacis, D. Jacobsen, and O. Dangles Worldwide, the rapid shrinking of glaciers in response to ongoing climate change is currently modifying the glacial meltwater contribution to hydrosystems in glacierized catchments. Assessing the contribution of glacier run-off to stream discharge is therefore of critical importance to evaluate potential impact of glacier retreat on water quality and aquatic biota. This task has challenged both glacier hydrologists and ecologists over the last 20 yr due to both structural and functional complexity of the glacier-stream system interface. Here we propose a new methodological approach based on wavelet analyses on water depth time series to determine the glacial influence in glacierized catchments. We performed water depth measurement using water pressure loggers over ten months in 15 stream sites in two glacier-fed catchments in the Ecuadorian Andes (〉 4000 m). We determined the global wavelet spectrum of each time series and defined the Wavelet Glacier Signal (WGS) as the ratio between the global wavelet power spectrum value at a 24 h-scale and its corresponding significance value. To test the relevance of the WGS we compared it with the percentage of the glacier cover in the catchments, a metric of glacier influence often used in the literature. We then tested whether one month data could be sufficient to reliably determine the glacial influence. As expected we found that the WGS of glacier-fed streams decreased downstream with the increasing of non-glacial tributaries. We also found that the WGS and the percentage of the glacier cover in the catchment were significantly positively correlated and that one month data was sufficient to identify and compare the glacial influence between two sites, provided that the water level time series were acquired over the same period. Furthermore, we found that our method permits to detect glacial signal in supposedly non-glacial sites, thereby evidencing glacial meltwater infiltrations. While we specifically focused on the tropical Andes in this paper, our approach to determine glacier influence would be applicable to temperate and arctic glacierized catchments. The WGS therefore appears as a powerful and cost effective tool to better understand the hydrological links between glaciers and hydrosystems and assess the consequences of rapid glacier melting.

Description: Spatio-temporal heterogeneity of riparian soil morphology in a restored floodplain Hydrology and Earth System Sciences Discussions, 10, 4337-4367, 2013 Author(s): B. Fournier, C. Guenat, G. Bullinger-Weber, and E. A. D. Mitchell Floodplains have been intensively altered in industrialized countries, but are now increasingly being restored and it is therefore important to assess the effect of these restoration projects on the aquatic and terrestrial components of ecosystems. Soils are a functionally crucial component of terrestrial ecosystems but are generally overlooked in floodplain restoration assessment. We studied the spatio-temporal heterogeneity of soil morphology in a restored (riverbed widening) river reach along River Thur (Switzerland) using three criteria (soil diversity, dynamism and typicality) and their associated indicators. We hypothesized that these criteria would correctly discriminate the post-restoration changes in soil morphology within the study site, and that these changes correspond to patterns of vascular plant diversity. Soil diversity and dynamism increased five years after the restoration, but typical soils of braided rivers were still missing. Soil typicality and dynamism correlated to vegetation changes. These results suggest a limited success of the project in agreement with evaluations carried out at the same site using other, more resource demanding methods (e.g. soil fauna, fish, ecosystem functioning). Soil morphology provides structural and functional information on floodplain ecosystems and allows predicting broad changes in plant diversity. The spatio-temporal heterogeneity of soil morphology represents a cost-efficient ecological indicator that could easily be integrated into rapid assessment protocols of floodplain and river restoration projects.

Description: Space-time kriging extension of precipitation variability at 12 km spacing from tree-ring chronologies and its implications for drought analysis Hydrology and Earth System Sciences Discussions, 10, 4301-4335, 2013 Author(s): F. Biondi Understanding and preparing for future hydroclimatic variability greatly benefits from long (i.e., multi-century) records at seasonal to annual time steps that have been gridded at km-scale spatial intervals over a geographic region. Kriging is a geostatistical technique commonly used for optimal interpolation of environmental data, and space-time geostatistical models can improve kriging estimates when long temporal sequences of observations exist at relatively few points on the landscape. Here I present how a network of 22 tree-ring chronologies from single-leaf pinyon ( Pinus monophylla ) in the central Great Basin of North America was used to extend hydroclimatic records both temporally and spatially. First, the Line of Organic Correlation (LOC) method was used to reconstruct October–May total precipitation anomalies at each tree-ring site, as these ecotonal environments at the lower forest border are typically moisture limited. Individual site reconstructions were then combined using a hierarchical model of spatio-temporal kriging that produced annual anomaly maps on a 12 × 12 km grid during the period in common among all chronologies (1650–1976). Hydro-climatic episodes were numerically identified and modeled using their duration, magnitude, and peak. Spatial patterns were more variable during wet years than during dry years, and the evolution of drought episodes over space and time could be visualized and quantified. The most remarkable episode in the entire reconstruction was the early 1900s pluvial, followed by the late 1800s drought. The 1930s "Dust Bowl" drought was among the top ten hydroclimatic episodes in the past few centuries. These results directly address the needs of water and natural resource managers with respect to planning for "worst case" scenarios of drought duration and magnitude at the watershed level. For instance, it is possible to analyze which geographical areas are more likely to be impacted by severe and sustained droughts at annual or multiannual timescales and at spatial resolutions commonly used by regional climate models.

Description: Structural break or long memory: an empirical survey on daily rainfall data sets across Malaysia Hydrology and Earth System Sciences, 17, 1311-1318, 2013 Author(s): F. Yusof, I. L. Kane, and Z. Yusop A short memory process that encounters occasional structural breaks in mean can show a slower rate of decay in the autocorrelation function and other properties of fractional integrated I (d) processes. In this paper we employed a procedure for estimating the fractional differencing parameter in semiparametric contexts proposed by Geweke and Porter-Hudak (1983) to analyse nine daily rainfall data sets across Malaysia. The results indicate that all the data sets exhibit long memory. Furthermore, an empirical fluctuation process using the ordinary least square (OLS)-based cumulative sum (CUSUM) test for the break date was applied. Break dates were detected in all data sets. The data sets were partitioned according to their respective break date, and a further test for long memory was applied for all subseries. Results show that all subseries follows the same pattern as the original series. The estimate of the fractional parameters d 1 and d 2 on the subseries obtained by splitting the original series at the break date confirms that there is a long memory in the data generating process (DGP). Therefore this evidence shows a true long memory not due to structural break.

Description: Snow glacier melt estimation in tropical Andean glaciers using artificial neural networks Hydrology and Earth System Sciences, 17, 1265-1280, 2013 Author(s): V. Moya Quiroga, A. Mano, Y. Asaoka, S. Kure, K. Udo, and J. Mendoza Snow and glacier melt (SGM) estimation plays an important role in water resources management. Although melting process can be modelled by energy balance methods, such studies require detailed data, which is rarely available. Hence, new and simpler approaches are needed for SGM estimations. The present study aims at developing an artificial neural networks (ANN) based technique for estimating the energy available for melt (EAM) and SGM rates using available and easy to obtain data such as temperature, short-wave radiation and relative humidity. Several ANN and multiple linear regression models (MLR) were developed to represent the energy fluxes and estimate the EAM. The models were trained using measured data from the Zongo glacier located in the outer tropics and validated against measured data from the Antizana glacier located in the inner tropics. It was found that ANN models provide a better generalisation when applied to other data sets. The performance of the models was improved by including Antizana data into the training set, as it was proved to provide better results than other techniques like the use of a prior logarithmic transformation. The final model was validated against measured data from the Alpine glaciers Argentière and Saint-Sorlin. Then, the models were applied for the estimation of SGM at Condoriri glacier. The estimated SGM was compared with SGM estimated by an enhanced temperature method and proved to have the same behaviour considering temperature sensibility. Moreover, the ANN models have the advantage of direct application, while the temperature method requires calibration of empirical coefficients.

Description: Optimising predictor domains for spatially coherent precipitation downscaling Hydrology and Earth System Sciences Discussions, 10, 4015-4061, 2013 Author(s): S. Radanovics, J.-P. Vidal, E. Sauquet, A. Ben Daoud, and G. Bontron Statistical downscaling is widely used to overcome the scale gap between predictors from Numerical Weather Prediction (NWP) models or General Circulation Models (GCMs) and predictands like local precipitation, required for example for medium-term operational forecasts or climate change impact studies. The predictors are considered over a given spatial domain which is rarely optimised with respect to the target predictand location. In this study the geopotential predictor domains used by an analogue downscaling method are optimised for 608 target zones covering France. An extended version of the growing rectangular domain algorithm provides an ensemble of five near-optimum domains for each target zone. All five near-optimum domains are consistently equally skillful based on the Continuous Rank Probability Score. Relevance maps calculated for selected target zones first reveal high skill geopotential regions with specific shapes for locations in south-eastern France compared to the rest of the country. In all cases, the optimised domains tend to include the most relevant area on the relevance maps. The domain centers of the optimised domains are mainly distributed following the geographical location of the target location, but there are apparent differences between the windward and the lee side of mountain ridges. Moreover, domains for target zones located in south-eastern France are centered more east and south than the ones for target locations on the same longitude. The size of the optimised domains tends to be larger in the southeastern part of the country, while domains with a very small meridional extent can be found in a east-west band around 47° N. Sensitivity tests on the archive length for the analogue method show a general robustness except for zones with high interannual variability like in the Cévennes area. Moreover, results appear to be rather unsensitive to the starting point of the optimisation algorithm except for zones located in the transition area north of the zones having optimized domains with a small meridional extent. This study paves the way for defining regions with homogeneous geopotential predictor domains for precipitation downscaling over France.

Description: Assessment of shallow subsurface characterisation with non-invasive geophysical methods at the intermediate hill-slope scale Hydrology and Earth System Sciences, 17, 1297-1307, 2013 Author(s): S. Popp, D. Altdorff, and P. Dietrich Hill-slopes of several hectares in size represent a difficult scale for subsurface characterisation, as these landscape units are well beyond the scope of traditional point-scale techniques. By means of electromagnetic induction (EMI) and gamma-ray spectroscopy, spatially distributed soil proxy data were collected from a heterogeneous hill-slope site. Results of EMI mapping using the EM38DD showed that soil electrical conductivity (ECa) is highly variable at both temporal and spatial scales. Calibration of the integral ECa signal to a specific target like soil moisture is hampered by the ambiguous response of EMI to the clay-rich hill-slope underground. Gamma-ray results were obtained during a single survey, along with EMI measurements and selected soil sampling. In contrast to ECa, a noticeable correlation between Total Count and K emission data and soil-water content seemed to be present. Relevant proxy variables from both methods were used for k means clustering in order to distinguish between hill-slope areas with different soil conditions. As a result, we obtained a suitable partition of hill-slope that was comparable with a previously obtained zonation model based on ecological factors.

Description: Multivariate return periods in hydrology: a critical and practical review focusing on synthetic design hydrograph estimation Hydrology and Earth System Sciences, 17, 1281-1296, 2013 Author(s): B. Gräler, M. J. van den Berg, S. Vandenberghe, A. Petroselli, S. Grimaldi, B. De Baets, and N. E. C. Verhoest Most of the hydrological and hydraulic studies refer to the notion of a return period to quantify design variables. When dealing with multiple design variables, the well-known univariate statistical analysis is no longer satisfactory, and several issues challenge the practitioner. How should one incorporate the dependence between variables? How should a multivariate return period be defined and applied in order to yield a proper design event? In this study an overview of the state of the art for estimating multivariate design events is given and the different approaches are compared. The construction of multivariate distribution functions is done through the use of copulas, given their practicality in multivariate frequency analyses and their ability to model numerous types of dependence structures in a flexible way. A synthetic case study is used to generate a large data set of simulated discharges that is used for illustrating the effect of different modelling choices on the design events. Based on different uni- and multivariate approaches, the design hydrograph characteristics of a 3-D phenomenon composed of annual maximum peak discharge, its volume, and duration are derived. These approaches are based on regression analysis, bivariate conditional distributions, bivariate joint distributions and Kendall distribution functions, highlighting theoretical and practical issues of multivariate frequency analysis. Also an ensemble-based approach is presented. For a given design return period, the approach chosen clearly affects the calculated design event, and much attention should be given to the choice of the approach used as this depends on the real-world problem at hand.

Description: Opportunities and challenges for the use of scintillometer-based catchment-averaged evapotranspiration estimates as model forcing Hydrology and Earth System Sciences Discussions, 10, 3973-4013, 2013 Author(s): B. Samain and V. R. N. Pauwels To date, lumped rainfall-runoff models rely on rough estimates of catchment-averaged potential evapotranspiration (ET p ) rates as meteorological forcing. A model parameter converts this ET p input into actual evapotranspiration (ET act ) estimates. This paper examines the potential use of scintillometer-based ET act rates for rainfall-runoff modeling. It has been found that the reservoir-structure of the rainfall-runoff model functions as a low-pass filter for the ET p input. If the long-term volume of the ET p used in the model simulations is consistent with the data set used for calibration, a good match of the seasonal pattern, using temporally constant ET p data, is sufficient to obtain adequate discharge simulations. However, these results are then obtained with strongly erroneous evapotranspiration estimates. A better match of the diurnal cycle does not lead to better model results. Replacing the ET p inputs by scintillometer-based ET act estimates does not lead to better model predictions. Small underestimations of ET act under stable conditions, which occur at night and during the Winter, and which accumulate to significant amounts, are the cause of this problem. Consistent with other studies, the scintillometer-based ET act estimates can be considered reliable and realistic under unstable conditions. These values can thus be used as forcing for rainfall-runoff models.

Description: Temporal stability of soil moisture patterns measured by proximal ground-penetrating radar Hydrology and Earth System Sciences Discussions, 10, 4063-4097, 2013 Author(s): J. Minet, N. E. C. Verhoest, S. Lambot, and M. Vanclooster We analyzed the temporal stability of soil moisture patterns acquired using a proximal ground-penetrating radar (GPR) in a 2.5 ha agricultural field at five different dates over three weeks. The GPR system was mounted on a mobile platform, allowing for real-time mapping of soil moisture with a high spatial resolution (2–5 m). The spatio-temporal soil moisture patterns were in accordance with the meteorological data and with soil moisture measurements from soil core sampling. Time-stable areas showing the field-average moisture could be revealed by two methods: (1) by the computation of temporal stability indicators based on relative differences of soil moisture to the field-average and (2) by the spatial intersection of the areas showing the field-average. Locations where the mean relative difference was below 0.02 m 3 m −3 extended up to 10% of the field area whereas the intersection of areas showing the field-average within a tolerance of 0.02 m 3 m −3 covered 5% of the field area. Compared to most of the previous studies about temporal stability of soil moisture, time-stable areas and their spatial patterns could be revealed instead of single point locations, owing to the advanced GPR method for real-time mapping. It is believed that determining spatially coherent time-stable areas is more informative rather than determining time-stable points. Other acquisitions over larger time periods would be necessary to assert the robustness of the time-stable areas.

Description: Balanced conditions or slight mass gain of glaciers in the Lahaul and Spiti region (northern India, Himalaya) during the nineties preceded recent mass loss The Cryosphere, 7, 569-582, 2013 Author(s): C. Vincent, Al. Ramanathan, P. Wagnon, D. P. Dobhal, A. Linda, E. Berthier, P. Sharma, Y. Arnaud, M. F. Azam, P. G. Jose, and J. Gardelle The volume change of the Chhota Shigri Glacier (India, 32° 20 N, 77° 30' E) between 1988 and 2010 has been determined using in situ geodetic measurements. This glacier has experienced only a slight mass loss between 1988 and 2010 (−3.8 ± 2.0 m w.e. (water equivalent) corresponding to −0.17 ± 0.09 m w.e. yr −1 ). Using satellite digital elevation models (DEM) differencing and field measurements, we measure a negative mass balance (MB) between 1999 and 2010 (−4.8 ± 1.8 m w.e. corresponding to −0.44 ± 0.16 m w.e. yr −1 ). Thus, we deduce a slightly positive or near-zero MB between 1988 and 1999 (+1.0 ± 2.7 m w.e. corresponding to +0.09 ± 0.24 m w.e. yr −1 ). Furthermore, satellite DEM differencing reveals that the MB of the Chhota Shigri Glacier (−0.39 ± 0.15 m w.e. yr −1 ) has been only slightly less negative than the MB of a 2110 km 2 glaciarized area in the Lahaul and Spiti region (−0.44 ± 0.09 m w.e. yr −1 ) during 1999−2011. Hence, we conclude that the ice wastage is probably moderate in this region over the last 22 yr, with near equilibrium conditions during the nineties, and an ice mass loss after. The turning point from balanced to negative mass budget is not known but lies probably in the late nineties and at the latest in 1999. This positive or near-zero MB for Chhota Shigri Glacier (and probably for the surrounding glaciers of the Lahaul and Spiti region) during at least part of the 1990s contrasts with a recent compilation of MB data in the Himalayan range that indicated ice wastage since 1975. However, in agreement with this compilation, we confirm more negative balances since the beginning of the 21st century.

Description: Supercooled interfacial water in fine grained soils probed by dielectric spectroscopy The Cryosphere Discussions, 7, 1441-1493, 2013 Author(s): A. Lorek and N. Wagner Water as thermodynamic state parameter affects nearly all physical, chemical and biological processes on the earth. Recent Mars observations as well as laboratory investigations suggest that water is also a key factor of current physical and chemical processes on the martian surface, e.g. rheological phenomena. Therefore it is of particular interest to get information about the liquid like state of water on martian analog soils in the temperature range below 0 °C. In this context, a parallel plate capacitor has been developed to obtain isothermal dielectric spectra of fine grained soils in the frequency range from 10 Hz to 1.1 MHz at martian like temperatures down to −70 °C. Two martian analogue soils have been investigated: a Ca-Bentonite (specific surface of 237 m 2 g −1 , up to 9.4% w/w gravimetric water content) and JSC Mars 1, a volcanic ash (specific surface of 146 m 2 g −1 , up to 7.4% w/w ). Three soil-specific relaxation processes are observed in the investigated frequency-temperature range: two weak high frequency processes (bound or hydrated water as well as ice) and a strong low frequency process due to counter ion relaxation and the Maxwell–Wagner effect. To characterize the dielectric relaxation behavior, a generalized fractional dielectric relaxation model is applied assuming three active relaxation processes with relaxation time of the i th process according to an Eyring equation. The real part of effective complex soil permittivity at 350 kHz was used to determine ice and liquid like water content by means of the Birchak or CRIM equation. There are evidence that Bentonite down to −70 °C has a liquid like water content of 1.17 monolayers and JSC Mars 1 a liquid like water content of 1.96 mono layers.

Description: Impacts of tropical cyclones on hydrochemistry of a subtropical forest Hydrology and Earth System Sciences Discussions, 10, 4537-4566, 2013 Author(s): C. T. Chang, S. P. Hamburg, J. L. Hwong, N. H. Lin, M. L. Hsueh, M. C. Chen, and T. C. Lin Tropical cyclones (typhoons/hurricanes) have major impacts on the biogeochemistry of forest ecosystems, but the stochastic nature and the long intervals between storms means that there are limited data on their effects. We characterized the impacts of 14 typhoons over six years on hydrochemistry of a subtropical forest plantation in Taiwan, a region experiencing frequent typhoons. Typhoons contributed 1/3 of annual rainfall on average, but ranged from 4% to 55%. The stochastic nature of annual typhoon related precipitation poses a challenge with respect to managing the impacts of these extreme events. This challenge is exacerbated by the fact that typhoon-related rainfall is not significantly correlated with wind velocity, the current focus of weather forecasts. Thus little advance warning is provided for the hydrological impacts of these storms. The typhoons we studied contributed approximately one third of the annual input and output of most nutrients (except nitrogen) during an average 9.5d yr −1 period, resulting in nutrient input/output rates an order of magnitude greater than during non-typhoon period. Nitrate output balanced input during the non-typhoon period, but during the typhoon period an average of 10 kg ha −1 yr −1 nitrate was lost. Streamwater chemistry exhibited similarly high variability during typhoon and non-typhoon periods and returned to pre-typhoon levels one to three weeks following each typhoon. The streamwater chemistry appears to be very resilient in response to typhoons, resulting in minimal loss of nutrients.

Description: Derivation of RCM-driven potential evapotranspiration for hydrological climate change impact analysis in Great Britain: a comparison of methods and associated uncertainty in future projections Hydrology and Earth System Sciences, 17, 1365-1377, 2013 Author(s): C. Prudhomme and J. Williamson Potential evapotranspiration (PET) is the water that would be lost by plants through evaporation and transpiration if water was not limited in the soil, and it is commonly used in conceptual hydrological modelling in the calculation of runoff production and hence river discharge. Future changes of PET are likely to be as important as changes in precipitation patterns in determining changes in river flows. However PET is not calculated routinely by climate models so it must be derived independently when the impact of climate change on river flow is to be assessed. This paper compares PET estimates from 12 equations of different complexity, driven by the Hadley Centre's HadRM3-Q0 model outputs representative of 1961–1990, with MORECS PET, a product used as reference PET in Great Britain. The results show that the FAO56 version of the Penman–Monteith equations reproduces best the spatial and seasonal variability of MORECS PET across GB when driven by HadRM3-Q0 estimates of relative humidity, total cloud, wind speed and linearly bias-corrected mean surface temperature. This suggests that potential biases in HadRM3-Q0 climate do not result in significant biases when the physically based FAO56 equations are used. Percentage changes in PET between the 1961–1990 and 2041–2070 time slices were also calculated for each of the 12 PET equations from HadRM3-Q0. Results show a large variation in the magnitude (and sometimes direction) of changes estimated from different PET equations, with Turc, Jensen–Haise and calibrated Blaney–Criddle methods systematically projecting the largest increases across GB for all months and Priestley–Taylor, Makkink, and Thornthwaite showing the smallest changes. We recommend the use of the FAO56 equation as, when driven by HadRM3-Q0 climate data, this best reproduces the reference MORECS PET across Great Britain for the reference period of 1961–1990. Further, the future changes of PET estimated by FAO56 are within the range of uncertainty defined by the ensemble of 12 PET equations. The changes show a clear northwest–southeast gradient of PET increase with largest (smallest) changes in the northwest in January (July and October) respectively. However, the range in magnitude of PET changes due to the choice of PET method shown in this study for Great Britain suggests that PET uncertainty is a challenge facing the assessment of climate change impact on hydrology mostly ignored up to now.

Description: The geomorphological effect of cornice fall avalanches in the Longyeardalen valley, Svalbard The Cryosphere, 7, 1361-1374, 2013 Author(s): M. Eckerstorfer, H. H. Christiansen, L. Rubensdotter, and S. Vogel The study of snow avalanches and their geomorphological effect in the periglacial parts of the cryosphere is important for enhanced geomorphological process understanding as well as hazard-related studies. Only a few field studies, and particularly few in the High Arctic, have quantified avalanche sedimentation. Snow avalanches are traditionally ranked behind rockfall in terms of their significance for mass-wasting processes of rockslopes. Cornice fall avalanches are at present the most dominant snow avalanche type at two slope systems, called Nybyen and Larsbreen, in the valley Longyeardalen in central Svalbard. Both slope systems are on northwest-facing lee slopes underneath a large summit plateau, with annual cornices forming on the top. High-frequency and magnitude cornice fall avalanching is observed by daily automatic time-lapse photography. In addition, rock debris sedimentation by cornice fall avalanches was measured directly in permanent sediment traps or by snow inventories. The results from a maximum of seven years of measurements in a total of 13 catchments show maximum mean rock debris sedimentation rates ranging from 8.2 to 38.7 kg m −2 at Nybyen, and from 0.8 to 55.4 kg m −2 at Larsbreen. Correspondingly, avalanche fan surfaces accreted from 2.6 to 8.8 mm yr −1 at Nybyen, and from 0.2 to 13.9 mm yr −1 at Larsbreen. This comparably efficient rockslope mass wasting is due to collapsing cornices producing cornice fall avalanches containing large amounts of rock debris throughout the entire winter. The rock debris of different origin stems from the plateau crests, the adjacent free rock face and the transport pathway, accumulating distinct avalanche fans at both slope systems. Cornice fall avalanche sedimentation also contributed to the development of a rock glacier at the Larsbreen site during the Holocene. We have recorded present maximum rockwall retreat rates of 0.9 mm yr −1 at Nybyen, but as much as 6.7 mm yr −1 at Larsbreen, while average Holocene rockwall retreat rates of 1.1 mm yr −1 at Nybyen have been determined earlier. As cornice fall avalanches are the dominant type of avalanche in central Svalbard, the related geomorphological effect is assumed to be of significance at periglacial landscape scale. A climate-induced shift in prevailing winter wind direction could change the rockslope sedimentation effectively by changing the snow avalanche activity.

Description: Influence of management of variables, sampling zones and land units on LR analysis for landslide spatial prevision Natural Hazards and Earth System Science, 13, 2209-2221, 2013 Author(s): R. Greco and M. Sorriso-Valvo Several authors, according to different methodological approaches, have employed logistic Regression (LR), a multivariate statistical analysis adopted to assess the spatial probability of landslide, even though its fundamental principles have remained unaltered. This study aims at assessing the influence of some of these methodological approaches on the performance of LR, through a series of sensitivity analyses developed over a test area of about 300 km 2 in Calabria (southern Italy). In particular, four types of sampling (1 – the whole study area; 2 – transects running parallel to the general slope direction of the study area with a total surface of about 1/3 of the whole study area; 3 – buffers surrounding the phenomena with a 1/1 ratio between the stable and the unstable area; 4 – buffers surrounding the phenomena with a 1/2 ratio between the stable and the unstable area), two variable coding modes (1 – grouped variables; 2 – binary variables), and two types of elementary land (1 – cells units; 2 – slope units) units have been tested. The obtained results must be considered as statistically relevant in all cases (Aroc values 〉 70%), thus confirming the soundness of the LR analysis which maintains high predictive capacities notwithstanding the features of input data. As for the area under investigation, the best performing methodological choices are the following: (i) transects produced the best results (0 〈 P(y) ≤ 93.4%; Aroc = 79.5%); (ii) as for sampling modalities, binary variables (0 〈 P(y) ≤ 98.3%; Aroc = 80.7%) provide better performance than ordinated variables; (iii) as for the choice of elementary land units, slope units (0 〈 P(y) ≤ 100%; Aroc = 84.2%) have obtained better results than cells matrix.

Description: Risk evolution: how can changes in the built environment influence the potential loss of natural hazards? Natural Hazards and Earth System Science, 13, 2195-2207, 2013 Author(s): B. Schwendtner, M. Papathoma-Köhle, and T. Glade Alpine areas often suffer significant loss and damage due to a range of natural processes such as landslides, debris flows, snow avalanches or floods. Sealing of the soil surface, settling in endangered areas and enhanced human intervention in the natural settings, as well as socio-economic changes, increase the risk and susceptibility of built environments to natural hazards and the costs of the consequences in a spatio-temporal context. The present study examines the loss estimation of a particular debris flow event for different points in time. The event occurred in August 1987, affected the municipality Martell in South Tyrol, Italy, and resulted in a total cost of € 25 million. The approach presented in this paper focuses on the changes of the land use and settlement expansion in the area since 1954 and attempts to assess the monetary impact of a similar event, which could have happened before (1954, 1985) or following the actual event (1992, 1999, 2006). The method applied is based on the use of a vulnerability curve which was developed for the specific area, based on the documentation of the damage of the 1987 event. Based on this curve, a loss estimation was carried out in order to visualise the risk evolution in a period of 52 yr (1954 to 2006). The results show a significant increase in the extent of the built environment (number, size and value of buildings) which consequently reflect an increase of the potential overall loss through the years. The method can be used in order to assess the potential loss for future scenarios based on different spatial patterns of the built environment.

Description: Wind shear over the Nice Côte d'Azur airport: case studies Natural Hazards and Earth System Science, 13, 2223-2238, 2013 Author(s): A. Boilley and J.-F. Mahfouf The Nice Côte d'Azur international airport is subject to horizontal low-level wind shears. Detecting and predicting these hazards is a major concern for aircraft security. A measurement campaign took place over the Nice airport in 2009 including 4 anemometers, 1 wind lidar and 1 wind profiler. Two wind shear events were observed during this measurement campaign. Numerical simulations were carried out with Meso-NH in a configuration compatible with near-real time applications to determine the ability of the numerical model to predict these events and to study the meteorological situations generating an horizontal wind shear. A comparison between numerical simulation and the observation dataset is conducted in this paper.

Description: Impact of long-term drainage on summer groundwater flow patterns in the Mer Bleue peatland, Ontario, Canada Hydrology and Earth System Sciences, 17, 3485-3498, 2013 Author(s): B. J. Kopp, J. H. Fleckenstein, N. T. Roulet, E. Humphreys, J. Talbot, and C. Blodau Long-term impacts of drier conditions on the hydrology of northern peatlands are poorly understood. We used long-term drainage near a historic drainage ditch, separating an area from the main peatland, as an analogue for long-term drying in a northern temperate bog. The objective was to identify the impact of drier conditions on ecohydrological processes and groundwater flow patterns in an area now forested and an area that maintained a bog-like character. Groundwater flow patterns alternated between mostly downward flow and occasionally upward flow in the bog area and were mostly upward-orientated in the forested area, which suggested that there the flow pattern had shifted from bog- to fen-like conditions. Flow patterns were in agreement with changes in post-drainage hydraulic conductivities, storage capacity of the peat and water table levels. Compared to the bog, hydraulic conductivities in the forested area were one to three orders of magnitude lower in the uppermost 0.75 m of peat (paired t test, p 〈 0.05). Bulk density had increased and the water table level was lower and more strongly fluctuating in the forested area. Our findings suggest hydraulic gradients and flow patterns have changed due to increased evapotranspiration and interception with the emergence of a tree cover. The smaller size of the now-forested area relative to the remaining bog area appeared to be important for the hydrological change. With the main Mer Bleue bog as hinterland, enhanced runoff to the drainage channel had little effect on hydrologic and vegetation patterns. In the cut-off, smaller, now forested area pervasive changes in vegetation and hydrologic processes occurred. The difference in response to local drainage raises questions about tipping points with respect to the impact of drying on peatland ecosystems that need to be addressed in future research.

Description: Hurricane impacts on a pair of coastal forested watersheds: implications of selective hurricane damage to forest structure and streamflow dynamics Hydrology and Earth System Sciences Discussions, 10, 11519-11557, 2013 Author(s): A. D. Jayakaran, T. M. Williams, H. Ssegane, D. M. Amatya, B. Song, and C. C. Trettin Hurricanes are infrequent but influential disruptors of ecosystem processes in the southeastern Atlantic and Gulf coasts. Every southeastern forested wetland has the potential to be struck by a tropical cyclone. We examined the impact of Hurricane Hugo on two paired coastal watersheds in South Carolina in terms of stream flow and vegetation dynamics, both before and after the hurricane's passage in 1989. The study objectives were to quantify the magnitude and timing of changes including a reversal in relative streamflow-difference between two paired watersheds, and to examine the selective impacts of a hurricane on the vegetative composition of the forest. We related these impacts to their potential contribution to change watershed hydrology through altered evapotranspiration processes. Using over thirty years of monthly rainfall and streamflow data we showed that there was a significant transformation in the hydrologic character of the two watersheds – a transformation that occurred soon after the hurricane's passage. We linked the change in the rainfall-runoff relationship to a catastrophic shift in forest vegetation due to selective hurricane damage. While both watersheds were located in the path of the hurricane, extant forest structure varied between the two watersheds as a function of experimental forest management techniques on the treatment watershed. We showed that the primary damage was to older pines, and to some extent larger hardwood trees. We believe that lowered vegetative water use impacted both watersheds with increased outflows on both watersheds due to loss of trees following hurricane impact. However, one watershed was able to recover to pre hurricane levels of canopy transpiration at a quicker rate due to the greater abundance of pine seedlings and saplings in that watershed.

Description: Explosive development of winter storm Xynthia over the subtropical North Atlantic Ocean Natural Hazards and Earth System Science, 13, 2239-2251, 2013 Author(s): M. L. R. Liberato, J. G. Pinto, R. M. Trigo, P. Ludwig, P. Ordóñez, D. Yuen, and I. F. Trigo In winter of 2009–2010 south-western Europe was hit by several destructive windstorms. The most important was Xynthia (26–28 February 2010), which caused 64 reported casualties and was classified as the 2nd most expensive natural hazard event for 2010 in terms of economic losses. In this work we assess the synoptic evolution, dynamical characteristics and the main impacts of storm Xynthia, whose genesis, development and path were very uncommon. Wind speed gusts observed at more than 500 stations across Europe are evaluated as well as the wind gust field obtained with a regional climate model simulation for the entire North Atlantic and European area. Storm Xynthia was first identified on 25 February around 30° N, 50° W over the subtropical North Atlantic Ocean. Its genesis occurred on a region characterized by warm and moist air under the influence of a strong upper level wave embedded in the westerlies. Xynthia followed an unusual SW–NE path towards Iberia, France and central Europe. The role of moist air masses on the explosive development of Xynthia is analysed by considering the evaporative sources. A lagrangian model is used to identify the moisture sources, sinks and moisture transport associated with the cyclone during its development phase. The main supply of moisture is located over an elongated region of the subtropical North Atlantic Ocean with anomalously high SST, confirming that the explosive development of storm Xynthia had a significant contribution from the subtropics.

Description: Potential and limitations of multidecadal satellite soil moisture observations for selected climate model evaluation studies Hydrology and Earth System Sciences, 17, 3523-3542, 2013 Author(s): A. Loew, T. Stacke, W. Dorigo, R. de Jeu, and S. Hagemann Soil moisture is an essential climate variable (ECV) of major importance for land–atmosphere interactions and global hydrology. An appropriate representation of soil moisture dynamics in global climate models is therefore important. Recently, a first multidecadal, observation-based soil moisture dataset has become available that provides information on soil moisture dynamics from satellite observations (ECVSM, essential climate variable soil moisture). The present study investigates the potential and limitations of this new dataset for several applications in climate model evaluation. We compare soil moisture data from satellite observations, reanalysis and simulations from a state-of-the-art land surface model and analyze relationships between soil moisture and precipitation anomalies in the different dataset. Other potential applications like model parameter optimization or model initialization are not investigated in the present study. In a detailed regional study, we show that ECVSM is capable to capture well the interannual and intraannual soil moisture and precipitation dynamics in the Sahelian region. Current deficits of the new dataset are critically discussed and summarized at the end of the paper to provide guidance for an appropriate usage of the ECVSM dataset for climate studies.

Description: Simultaneous estimation of model state variables and observation and forecast biases using a two-stage hybrid Kalman filter Hydrology and Earth System Sciences, 17, 3499-3521, 2013 Author(s): V. R. N. Pauwels, G. J. M. De Lannoy, H.-J. Hendricks Franssen, and H. Vereecken In this paper, we present a two-stage hybrid Kalman filter to estimate both observation and forecast bias in hydrologic models, in addition to state variables. The biases are estimated using the discrete Kalman filter, and the state variables using the ensemble Kalman filter. A key issue in this multi-component assimilation scheme is the exact partitioning of the difference between observation and forecasts into state, forecast bias and observation bias updates. Here, the error covariances of the forecast bias and the unbiased states are calculated as constant fractions of the biased state error covariance, and the observation bias error covariance is a function of the observation prediction error covariance. In a series of synthetic experiments, focusing on the assimilation of discharge into a rainfall-runoff model, it is shown that both static and dynamic observation and forecast biases can be successfully estimated. The results indicate a strong improvement in the estimation of the state variables and resulting discharge as opposed to the use of a bias-unaware ensemble Kalman filter. Furthermore, minimal code modification in existing data assimilation software is needed to implement the method. The results suggest that a better performance of data assimilation methods should be possible if both forecast and observation biases are taken into account.

Description: Seismic behavior of buried pipelines constructed by design criteria and construction specifications of both Korea and the US Natural Hazards and Earth System Science, 13, 2271-2278, 2013 Author(s): S.-S. Jeon Earthquake loss estimation systems in the US, for example HAZUS (Hazard in US), have been established based on sufficient damage records for the purpose of prevention and efficient response to earthquake hazards; however, in Korea, insufficient data sets of earthquakes and damage records are currently available. In this study, the earthquake damages to pipelines in Korea using the pipeline repair rate (RR) recommended in HAZUS was reevaluated with the degree of confidence when RR is used without modification for the damage estimation of pipelines in Korea. The numerical analyses using a commercial finite element model, ABAQUS, were carried out to compare stresses and strains mobilized in both brittle and ductile pipelines constructed by the design criteria and construction specifications of both Korea and the US. These pipelines were embedded in dense sand overlying three different in situ soils (clay, sand, and gravel) subjected to earthquake excitations with peak ground accelerations (PGAs) of 0.2 to 1.2 g and 1994 Northridge and 1999 Chi-Chi earthquake loadings. The numerical results show that differences in the stress and strain rates are less than 10%. This implies that RR in HAZUS can be used for earthquake damage estimation of pipelines with a 90% confidence level in Korea.

Description: Heterogeneity in glacier response in the upper Shyok valley, northeast Karakoram The Cryosphere, 7, 1385-1398, 2013 Author(s): R. Bhambri, T. Bolch, P. Kawishwar, D. P. Dobhal, D. Srivastava, and B. Pratap Glaciers in the Karakoram show long-term irregular behaviour with comparatively frequent and sudden advances. A glacier inventory of the upper Shyok valley situated in northeast Karakoram has been generated for the year 2002 using Landsat ETM+ and SRTM3 DEM as baseline data for the investigations and subsequent change analysis. The upper Shyok valley contained 2123 glaciers (larger than 0.02 km 2 in size) with an area of 2977.9 ± 95.3 km 2 in 2002. Out of these, 18 glaciers with an area of 1004.1 ± 32.1 km 2 showed surge-type behaviour. Change analysis based on Hexagon KH-9 (years 1973 and 1974) and Landsat TM/ETM+ (years 1989, 2002 and 2011) images had to be restricted to a subset of 136 glaciers (covering an area of 1609.7 ± 51.5 km 2 in 2002) due to adverse snow conditions. The area of the investigated glaciers, including the 18 surge-type glaciers identified, showed no significant changes during all studied periods. However, the analysis provides a hint that the overall glacier area slightly decreased until about 1989 (area 1973: 1613.6 ± 43.6 km 2 ; area 1989: 1602.0 ± 33.6 km 2 ) followed by an increase (area 2002: 1609.7 ± 51.5; area 2011: 1615.8 ± 35.5 km 2 ). Although the overall change in area is insignificant, advances in glacier tongues since the end of the 1980s are clearly visible. Detailed estimations of length changes for individual glaciers since the 1970s and for Central Rimo Glacier since the 1930s confirm the irregular retreat and advance.

Description: Influence of snow depth distribution on surface roughness in alpine terrain: a multi-scale approach The Cryosphere Discussions, 7, 4633-4680, 2013 Author(s): J. Veitinger, B. Sovilla, and R. S. Purves In alpine terrain, the snow covered winter surface deviates from its underlying summer terrain due to the progressive smoothing caused by snow accumulation. Terrain smoothing is believed to be an important factor in avalanche formation, avalanche dynamics and affects surface heat transfer, energy balance as well as snow depth distribution. To characterize the effect of snow on terrain we use the concept of roughness. Roughness is calculated for several snow surfaces and its corresponding underlying terrain for three alpine basins in the Swiss Alps characterized by low medium and high terrain roughness. To this end, elevation models of winter and summer terrain are derived from high-resolution (1 m) measurements performed by airborne and terrestrial LIDAR. We showed that on basin scale terrain smoothing not only depends on mean snow depth in the basin but also on its variability. Terrain smoothing can be modelled in function of mean snow depth and its standard deviation using a power law. However, a relationship between terrain smoothing and snow depth does not exist on a pixel scale. Further we demonstrated the high persistence of snow surface roughness even in between winter seasons. Those persistent patterns might be very useful to improve the representation of a winter terrain without modelling of the snow cover distribution. This can potentially improve avalanche release area definition and in the long term natural hazard management strategies.

Description: Land use change effects on runoff generation in a humid tropical montane cloud forest region Hydrology and Earth System Sciences, 17, 3543-3560, 2013 Author(s): L. E. Muñoz-Villers and J. J. McDonnell While tropical montane cloud forests (TMCF) provide critical hydrological services to downstream regions throughout much of the humid tropics, catchment hydrology and impacts associated with forest conversion in these ecosystems remain poorly understood. Here, we compare the annual, seasonal and event-scale streamflow patterns and runoff generation processes of three neighbouring headwater catchments in central Veracruz (eastern Mexico) with similar pedological and geological characteristics, but different land cover: old-growth TMCF, 20 yr-old naturally regenerating TMCF and a heavily grazed pasture. We used a 2 yr record of high resolution rainfall and stream flow data (2008–2010) in combination with stable isotope and chemical tracer data collected for a series of storms during a 6-week period of increasing antecedent wetness (wetting-up cycle). Our results showed that annual and seasonal streamflow patterns in the mature and secondary forest were similar. In contrast, the pasture showed a 10% higher mean annual streamflow, most likely because of a lower rainfall interception. During the wetting-up cycle, storm runoff ratios increased at all three catchments (from 11 to 54% for the mature forest, 7 to 52% for the secondary forest and 3 to 59% for the pasture). With the increasing antecedent wetness, hydrograph separation analysis showed progressive increases of pre-event water contributions to total stormflow (from 35 to 99% in the mature forest, 26 to 92% in the secondary forest and 64 to 97% in the pasture). At all three sites, rainfall-runoff responses were dominated by subsurface flow generation processes for the majority of storms. However, for the largest and most intense storm (typically occurring once every 2 yr), sampled under wet antecedent conditions, the event water contribution in the pasture (34% on average) was much higher than in the forests (5% on average), indicating that rainfall infiltration capacity of the pasture was exceeded. This result suggests that despite the high permeability of the volcanic soils and underlying substrate in this TMCF environment, the conversion of forest to pasture may lead to important changes in runoff generation processes during large and high intensity storms. On the other hand, our results also showed that 20 yr of natural regeneration may be enough to largely restore the original hydrological conditions of this TMCF.

Description: Recent evolution of China's virtual water trade: analysis of selected crops and considerations for policy Hydrology and Earth System Sciences Discussions, 10, 11613-11641, 2013 Author(s): J. Shi, J. Liu, and L. Pinter China has dramatically increased its virtual water import unconsciously for recent years. Many studies have focused on the quantity of traded virtual water but very few go into analysing geographic distribution and the properties of China's virtual water trade network. This paper provides a calculation and analysis of the crop-related virtual water trade network of China based on 27 major primary crops between 1986 and 2009. The results show that China is a net importer of virtual water from water-abundant areas of North and South America, and a net virtual water exporter to water-stressed areas of Asia, Africa, and Europe. Virtual water import is far larger than virtual water export and in both import and export a small number of trade partners control the supply chain. Grain crops are the major contributors to virtual water trade, and among grain crops soybeans, mostly imported from the US, Brazil and Argentina are the most significant. As crop yield and crop water productivity in North and South America are generally higher than those in Asia and Africa, the effect of China's crop-related virtual water trade positively contributes to optimizing crop water use efficiency at the global scale. In order to mitigate water scarcity and secure the food supply, virtual water should be actively incorporated into national water management strategies. From the national perspective, China should reduce the export and increase the import of water-intensive crops. But the sources of virtual water import need to be further diversified to reduce supply chain risks and increase resilience.

Description: Weak layer fracture: facets and depth hoar The Cryosphere, 7, 1447-1453, 2013 Author(s): I. Reiweger and J. Schweizer Understanding failure initiation within weak snow layers is essential for modeling and predicting dry-snow slab avalanches. We therefore performed laboratory experiments with snow samples containing a weak layer consisting of either faceted crystals or depth hoar. During these experiments the samples were loaded with different loading rates and at various tilt angles until fracture. The strength of the samples decreased with increasing loading rate and increasing tilt angle. Additionally, we took pictures of the side of four samples with a high-speed video camera and calculated the displacement using a particle image velocimetry (PIV) algorithm. The fracture process within the weak layer could thus be observed in detail. Catastrophic failure started due to a shear fracture just above the interface between the depth hoar layer and the underlying crust.

Description: Scatter of mass changes estimates at basin scale for Greenland and Antarctica The Cryosphere, 7, 1411-1432, 2013 Author(s): V. R. Barletta, L. S. Sørensen, and R. Forsberg During the last decade, the GRACE mission has provided valuable data for determining the mass changes of the Greenland and Antarctic ice sheets. Yet, discrepancies still exist in the published mass balance results, and comprehensive analyses on the sources of errors and discrepancies are lacking. Here, we present monthly mass changes together with trends derived from GRACE data at basin scale for both the Greenland and Antarctic ice sheets, and we assess the variability and errors for each of the possible sources of discrepancies, and we do this in an unprecedented systematic way, taking into account mass inference methods, data sets and background models. We find a very good agreement between the monthly mass change results derived from two independent methods, which represents a cross validation. For the monthly solutions, we find that most of the scatter is caused by the use of the two different data sets rather than the two different methods applied. Besides the well-known GIA trend uncertainty, we find that the geocenter motion and the recent de-aliasing corrections significantly impact the trends, with contributions of +13.2 Gt yr −1 and −20 Gt yr −1 , respectively, for Antarctica, which is more affected by these than Greenland. We show differences between the use of release RL04 and the new RL05 and confirm a lower noise content in the new release. The overall scatter of the solutions well exceeds the uncertainties propagated from the data errors and the leakage (as done in the past); hence we calculate new sound total errors for the monthly solutions and the trends. We find that the scatter in the monthly solutions caused by applying different estimates of geocenter motion time series (degree-1 corrections) is significant – contributing with up to 40% of the total error. For the whole GRACE period (2003–2011) our trend estimate for Greenland is −234 ± 20 Gt yr −1 and −83 ± 36 Gt yr −1 for Antarctica (−111 ± 15 Gt yr −1 in the western part). We also find a clear (with respect to our errors) increase of mass loss in the last four years.

Description: Seismic zones for Azores based on statistical criteria Natural Hazards and Earth System Science, 13, 2337-2351, 2013 Author(s): M. C. M. Rodrigues and C. S. Oliveira The objective of this paper is to define seismic zones in the Azores based on statistical criteria. These seismic zones will likely be used in seismic simulations of occurrences in the Azores Archipelago. The data used in this work cover the time period from 1915 to 2011. The Azores region was divided into 1° × 1° area units, for which the seismicity and the maximum magnitudes of events were calculated. The seismicity, the largest earthquakes recorded and the geological characteristics of the region were used to group these area units because similar seismic zones must delineate areas with homogeneous seismic characteristics. We have identified seven seismic zones. To verify that the defined areas differ statistically, we considered the following dissimilarity measures (variables): time, size and seismic conditions – the number of seismic events with specific characteristics. Statistical tests, particularly goodness-of-fit tests, allowed us to conclude that, considering these three variables, the seven earthquake zones defined here are statistically distinct.

Description: Sub-ionospheric very low frequency perturbations associated with the 12 May 2008 M = 7.9 Wenchuan earthquake Natural Hazards and Earth System Science, 13, 2331-2336, 2013 Author(s): A. K. Maurya, R. Singh, B. Veenadhari, S. Kumar, and A. K. Singh The present study reports the VLF (very low frequency) sub-ionospheric perturbations observed on transmitter JJI (22.1 kHz), Japan, received at the Indian low-latitude station, Allahabad ( geographic lat. 25.41° N, long 81.93° E), due to Wenchuan earthquake (EQ) that occurred on 12 May 2008 with the magnitude 7.9 and at the depth of 19 km in Sichuan province of Southwest China, located at 31.0° N, 103.4° E. The nighttime amplitude fluctuation analysis gives a significant increase in fluctuation and dispersion two days before EQ, when it crosses 2σ criterion. However, there was no significant change observed in the amplitude trend. The diurnal amplitude variation shows a significant increase in the amplitude of JJI signal on 11 and 12 May 2008. The gravity wave channel and changes in the electric field associated with this EQ seem to be the potential factors of the observed nighttime amplitude fluctuation, dispersion, and significant increase in the signal strength.

Description: Using the nonlinear aquifer storage–discharge relationship to simulate the base flow of glacier- and snowmelt-dominated basins in northwest China Hydrology and Earth System Sciences, 17, 3577-3586, 2013 Author(s): R. Gan and Y. Luo Base flow is an important component in hydrological modeling. This process is usually modeled by using the linear aquifer storage–discharge relation approach, although the outflow from groundwater aquifers is nonlinear. To identify the accuracy of base flow estimates in rivers dominated by snowmelt and/or glacier melt in arid and cold northwestern China, a nonlinear storage–discharge relationship for use in SWAT (Soil Water Assessment Tool) modeling was developed and applied to the Manas River basin in the Tian Shan Mountains. Linear reservoir models and a digital filter program were used for comparisons. Meanwhile, numerical analysis of recession curves from 78 river gauge stations revealed variation in the parameters of the nonlinear relationship. It was found that the nonlinear reservoir model can improve the streamflow simulation, especially for low-flow period. The higher Nash–Sutcliffe efficiency, logarithmic efficiency, and volumetric efficiency, and lower percent bias were obtained when compared to the one-linear reservoir approach. The parameter b of the aquifer storage–discharge function varied mostly between 0.0 and 0.1, which is much smaller than the suggested value of 0.5. The coefficient a of the function is related to catchment properties, primarily the basin and glacier areas.

Description: Technical Note: Reducing the spin-up time of integrated surface water–groundwater models Hydrology and Earth System Sciences, 18, 5169-5179, 2014 Author(s): H. Ajami, J. P. Evans, M. F. McCabe, and S. Stisen One of the main challenges in the application of coupled or integrated hydrologic models is specifying a catchment's initial conditions in terms of soil moisture and depth-to-water table (DTWT) distributions. One approach to reducing uncertainty in model initialization is to run the model recursively using either a single year or multiple years of forcing data until the system equilibrates with respect to state and diagnostic variables. However, such "spin-up" approaches often require many years of simulations, making them computationally intensive. In this study, a new hybrid approach was developed to reduce the computational burden of the spin-up procedure by using a combination of model simulations and an empirical DTWT function. The methodology is examined across two distinct catchments located in a temperate region of Denmark and a semi-arid region of Australia. Our results illustrate that the hybrid approach reduced the spin-up period required for an integrated groundwater–surface water–land surface model (ParFlow.CLM) by up to 50%. To generalize results to different climate and catchment conditions, we outline a methodology that is applicable to other coupled or integrated modeling frameworks when initialization from an equilibrium state is required.

Description: Extending periodic eddy covariance latent heat fluxes through tree sapflow measurements to estimate long-term total evaporation in a peat swamp forest Hydrology and Earth System Sciences Discussions, 11, 13607-13661, 2014 Author(s): A. D. Clulow, C. S. Everson, M. G. Mengistu, J. S. Price, A. Nickless, and G. P. W. Jewitt A combination of measurement and modelling was used to find a pragmatic solution to estimate the annual total evaporation (ET) from the rare and indigenous Nkazana Peat Swamp Forest (PSF) on the east coast of Southern Africa to improve the water balance estimates within the area. Total evaporation was measured during three window periods (between seven and nine days each) using an eddy covariance (EC) system on a telescopic mast above the forest canopy. Sapflow of an understory and an emergent tree was measured using a low maintenance heat pulse velocity system for an entire hydrological year (October 2009 to September 2010). An empirical model was derived, describing the relationship between the observed ET of the Nkazana PSF measured during two of the window periods ( R 2 = 0.92 and 0.90) which, overlapped with sapflow measurements, thereby providing hourly estimates of predicted ET of the Nkazana PSF for a year, totalling 1125 mm (while rainfall was 650 mm). In building the empirical model, it was found that including the understory tree sapflow provided no benefit to the model performance. In addition, the observed emergent tree sapflow relationship with observed ET between the two field campaigns was consistent and could be represented by a single empirical model ( R 2 = 0.90; RMSE = 0.08 mm). During the window periods of EC measurement, no single meteorological variable was found to describe the Nkazana PSF ET satisfactorily. However, in terms of evaporation models, the hourly FAO56 Penman–Monteith equation best described the observed ET from EC during the August 2009 ( R 2 = 0.75), November 2009 ( R 2 = 0.85) and March 2010 ( R 2 = 0.76) field campaigns, compared to the Priestley–Taylor model ( R 2 = 0.54, 0.74 and 0.62 during the respective field campaigns). From the empirical model of ET and the FAO56 Penman–Monteith equation, a monthly crop factor ( K c ) was derived for the Nkazana PSF providing a method of estimating long-term swamp forest ET from meteorological data. The monthly crop factor indicated two distinct periods. From February to May, it was between 1.2 and 1.4 compared with June to January, when the crop factor was 0.8 to 1.0. The derived monthly K c values were verified as accurate (to one significant digit) using historical data measured at the same site, also using EC, from a~previous study. The measurements provided insights into the microclimate within a subtropical peat swamp forest and the contrasting sapflow of emergent and understory trees. They showed that expensive, high maintenance equipment can be used during manageable window periods in conjunction with low maintenance systems, dedicated to individual trees, to derive a model to estimate long-term ET over remote heterogeneous forests. In addition, the contrast in ET and rainfall emphasises the reliance of the Nkazana PSF on groundwater.

Description: Long-term variability of storm surge frequency in the Venice Lagoon: an update thanks to eighteenth century sea level observations Natural Hazards and Earth System Sciences Discussions, 2, 7465-7486, 2014 Author(s): F. Raicich Sea level observations made in the Venice Lagoon between 1751 and 1792 have been recovered, consisting of two time series of daily data on high and low waters at Venice and Chioggia. From comparisons with modern observations the quality of the 18th century data appears to be good enough to allow a useful analysis. A composite time series of daily mean sea level is obtained by merging the 18th century data and 1872–2004 observations at Venice Punta della Salute. The absence of reliable information on vertical references prevents the connection of the two 18th century time series with each other and with modern observations. However, daily sea level anomalies relative to the mean sea level enable to recognize storm surge events, that appear to occur more frequently in the second half of the 18th century than in the late 19th and 20th centuries, particularly during the 1751–1769 period. The record-breaking storm surge of 4 November 1966 turns out to be a remarkable event also in comparison with the events extracted from the 18th century time series. Further work is required to fill the gap between the old and modern observations.

Description: Ground-penetrating radar observations for estimating the vertical displacement of rotational landslides Natural Hazards and Earth System Sciences Discussions, 2, 7487-7506, 2014 Author(s): C. Lissak, O. Maquaire, J.-P. Malet, F. Lavigne, C. Virmoux, C. Gomez, and R. Davidson The objective of this paper is to demonstrate the applicability of Ground Penetrating Radar (GPR) for monitoring the displacement of slow-moving landslides. GPR data is used to estimate the vertical movement of rotational slides in combination with other surveying techniques. The experimental site is located along the Normandy coast (North East France) here several rotational landslides are continuously affected by a seasonal kinematic pattern (low displacement rates of 0.01 to 0.10 m yr −1 ) and periodically by major acceleration events (high displacement of 1.0 to 7.0 m per event).

Description: Continental Portuguese Territory Flood Susceptibility Index – contribution for a vulnerability index Natural Hazards and Earth System Sciences Discussions, 2, 7521-7552, 2014 Author(s): R. Jacinto, N. Grosso, E. Reis, L. Dias, F. D. Santos, and P. Garrett This work defines a national flood susceptibility index for the Portuguese continental territory, by proposing the aggregation of different variables which represent natural conditions for permeability, runoff and accumulation. This index is part of the national vulnerability index developed in the scope of Flood Maps in Climate Change Scenarios (CIRAC) project, supported by the Portuguese Association of Insurers (APS). This approach expands on previous works by trying to bridge the gap between different floods mechanisms (e.g. progressive and flash floods) occurring at different spatial scales in the Portuguese territory through: (a) selecting homogeneously processed datasets, (b) aggregating their values to better translate the spatially continuous and cumulative influence in floods at multiple spatial scales. Results show a good ability to capture, in the higher susceptibility classes, different flood types: progressive floods and flash floods. Lower values are usually related to: mountainous areas, low water accumulation potential and more permeable soils. Validation with independent flood datasets confirmed these index characteristics, although some overestimation can be seen in the southern region of Alentejo where, due to a dense hydrographic network and an overall low slope, floods are not as frequent as a result of lower precipitation mean values. Future work will focus on: (i) including extreme precipitation datasets to represent the triggering factor, (ii) improving representation of smaller and stepper basins, (iii) optimizing variable weight definition process, (iii) developing more robust independent flood validation datasets.

Description: Dye tracing for investigating flow and transport properties of hydrocarbon-polluted Rabots glaciär, Kebnekaise, Sweden Hydrology and Earth System Sciences Discussions, 11, 13711-13744, 2014 Author(s): C. C. Clason, C. Coch, J. Jarsjö, K. Brugger, P. Jansson, and G. Rosqvist Over 11 000 L of hydrocarbon pollution was deposited on the surface of Rabots glaciär on the Kebnekaise Massif, northern Sweden, following the crash of a Royal Norwegian Air Force aircraft in March 2012. An environmental monitoring programme was subsequently commissioned, including water, snow and ice sampling. The scientific programme further included a series of dye tracing experiments during the 2013 melt season, conducted to investigate flow pathways for pollutants through the glacier hydrological system, and to gain new insight to the internal hydrological system of Rabots glaciär. Results of dye tracing reveal a degree of homogeneity in the topology of the drainage system throughout July and August, with an increase in efficiency as the season progresses, as reflected by decreasing temporary storage and dispersivity. Early onset of melting likely led to formation of an efficient, discrete drainage system early in the melt season, subject to decreasing sinuosity and braiding as the season progressed. Analysis of turbidity-discharge hysteresis further supports the formation of discrete, efficient drainage, with clockwise diurnal hysteresis suggesting easy mobilisation of readily-available sediments in channels. Dye injection immediately downstream of the pollution source zone revealed prolonged storage of dye followed by fast, efficient release. Twinned with a low dye recovery, and supported by sporadic detection of hydrocarbons in the proglacial river, we suggest that meltwater, and thus pollutants in solution, may be released periodically from this zone of the glacier hydrological system. The here identified dynamics of dye storage, dispersion and breakthrough indicate that the ultimate fate and permanence of pollutants in the glacier system is likely to be governed by storage of pollutants in the firn layer and ice mass, or within the internal hydrological system, where it may refreeze. This shows that future studies on the fate of hydrocarbons in pristine, glaciated mountain environments should address the extent to which pollutants in solution act like water molecules or whether they are more susceptible to, for example, refreezing into the surrounding ice, becoming stuck in micro-fractures and pore spaces, or sorption onto subglacial sediments.

Description: Air temperature variability over three glaciers in the Ortles-Cevedale (Italian Alps): effects of glacier disintegration, intercomparison of calculation methods, and impacts on mass balance modeling The Cryosphere Discussions, 8, 6147-6192, 2014 Author(s): L. Carturan, F. Cazorzi, F. De Blasi, and G. Dalla Fontana Glacier mass balance models rely on accurate spatial calculation of input data, in particular air temperature. Lower temperatures (the so-called glacier cooling effect), and lower temperature variability (the so-called glacier damping effect) generally occur over glaciers, compared to ambient conditions. These effects, which depend on the geometric characteristics of glaciers and display a high spatial and temporal variability, have been mostly investigated on medium- to large-size glaciers so far, while observations on smaller ice bodies are scarce. Using a dataset from 8 on-glacier and 4 off-glacier weather stations, collected in summer 2010 and 2011, we analyzed the air temperature variability and wind regime over three different glaciers in the Ortles-Cevedale. The magnitude of the cooling effect and the occurrence of katabatic boundary layer (KBL) processes showed remarkable differences among the three ice bodies, suggesting the likely existence of important reinforcing mechanisms during glacier decay and disintegration. None of the methods proposed in the literature for calculating on-glacier temperature from off-glacier data fully reproduced our observations. Among them, the more physically-based procedure of Greuell and Böhm (1998) provided the best overall results where the KBL prevail, but it was not effective elsewhere (i.e. on smaller ice bodies and close to the glacier margins). The accuracy of air temperature estimations strongly impacted the results from a mass balance model which was applied to the three investigated glaciers. Most importantly, even small temperature deviations caused distortions in parameter calibration, thus compromising the model generalizability.

Description: Are seasonal calving dynamics forced by buttressing from ice mélange or undercutting by melting? Outcomes from full-Stokes simulations of Store Glacier, West Greenland The Cryosphere, 8, 2353-2365, 2014 Author(s): J. Todd and P. Christoffersen We use a full-Stokes 2-D model (Elmer/Ice) to investigate the flow and calving dynamics of Store Glacier, a fast-flowing outlet glacier in West Greenland. Based on a new, subgrid-scale implementation of the crevasse depth calving criterion, we perform two sets of simulations: one to identify the primary forcing mechanisms and another to constrain future stability. We find that the mixture of icebergs and sea ice, known as ice mélange or sikussak, is principally responsible for the observed seasonal advance of the ice front. On the other hand, the effect of submarine melting on the calving rate of Store Glacier appears to be limited. Sensitivity analysis demonstrates that the glacier's calving dynamics are sensitive to seasonal perturbation, but are stable on interannual timescales due to the strong topographic control on the flow regime. Our results shed light on the dynamics of calving glaciers and may help explain why neighbouring glaciers do not necessarily respond synchronously to changes in atmospheric and oceanic forcing.

Description: Deglaciation of the Caucasus Mountains, Russia/Georgia, in the 21st century observed with ASTER satellite imagery and aerial photography The Cryosphere, 8, 2367-2379, 2014 Author(s): M. Shahgedanova, G. Nosenko, S. Kutuzov, O. Rototaeva, and T. Khromova Changes in the map area of 498 glaciers located on the Main Caucasus ridge (MCR) and on Mt. Elbrus in the Greater Caucasus Mountains (Russia and Georgia) were assessed using multispectral ASTER and panchromatic Landsat imagery with 15 m spatial resolution in 1999/2001 and 2010/2012. Changes in recession rates of glacier snouts between 1987–2001 and 2001–2010 were investigated using aerial photography and ASTER imagery for a sub-sample of 44 glaciers. In total, glacier area decreased by 4.7 ± 2.1% or 19.2 ± 8.7 km 2 from 407.3 ± 5.4 km 2 to 388.1 ± 5.2 km 2 . Glaciers located in the central and western MCR lost 13.4 ± 7.3 km 2 (4.7 ± 2.5%) in total or 8.5 km 2 (5.0 ± 2.4%) and 4.9 km 2 (4.1 ± 2.7%) respectively. Glaciers on Mt. Elbrus, although located at higher elevations, lost 5.8 ± 1.4 km 2 (4.9 ± 1.2%) of their total area. The recession rates of valley glacier termini increased between 1987–2000/01 and 2000/01–2010 (2000 for the western MCR and 2001 for the central MCR and Mt.~Elbrus) from 3.8 ± 0.8, 3.2 ± 0.9 and 8.3 ± 0.8 m yr −1 to 11.9 ± 1.1, 8.7 ± 1.1 and 14.1 ± 1.1 m yr −1 in the central and western MCR and on Mt. Elbrus respectively. The highest rate of increase in glacier termini retreat was registered on the southern slope of the central MCR where it has tripled. A positive trend in summer temperatures forced glacier recession, and strong positive temperature anomalies in 1998, 2006, and 2010 contributed to the enhanced loss of ice. An increase in accumulation season precipitation observed in the northern MCR since the mid-1980s has not compensated for the effects of summer warming while the negative precipitation anomalies, observed on the southern slope of the central MCR in the 1990s, resulted in stronger glacier wastage.

Description: Inversion method for initial tsunami waveform reconstruction Natural Hazards and Earth System Sciences Discussions, 2, 7735-7772, 2014 Author(s): V. V. Voronin, T. A. Voronina, and V. A. Tcheverda This paper deals with the application of r-solution method to recover the initial tsunami waveform in a tsunami source area by remote water-level measurements. Wave propagation is considered within the scope of a linear shallow-water theory. An ill-posed inverse problem is regularized by means of least square inversion using a truncated SVD approach. The properties of obtained solution are determined to a large extent by the properties of an inverse operator, which were numerically investigated. The method presented allows one to control instability of the numerical solution and to obtain an acceptable result in spite of ill-posedness of the problem. It is shown that the accuracy of tsunami source reconstruction strongly depends on the signal-to-noise ratio, the azimuthal coverage of recording stations with respect to the source area and bathymetric features along the wave path. The numerical experiments were carried out with synthetic data and various computational domains including a real bathymetry. The method proposed allows us to make a preliminary prediction of the efficiency of the inversion with a given set of the recording stations and to find out the most informative part of the existing observation system. This essential property of the method can prove to be useful in designing a monitoring system for tsunamis.

Description: Glacier-like forms on Mars The Cryosphere, 8, 2047-2061, 2014 Author(s): B. Hubbard, C. Souness, and S. Brough More than 1300 glacier-like forms (GLFs) are located in Mars' mid-latitudes. These GLFs are predominantly composed of ice–dust mixtures and are visually similar to terrestrial valley glaciers, showing signs of downhill viscous deformation and an expanded former extent. However, several fundamental aspects of their behavior are virtually unknown, including temporal and spatial variations in mass balance, ice motion, landscape erosion and deposition, and hydrology. Here, we investigate the physical glaciology of martian GLFs. We use satellite images of specific examples and case studies to build on existing knowledge relating to (i) GLF current and former extent, exemplified via a GLF located in Phlegra Montes; (ii) indicators of GLF motion, focusing on the presence of surface crevasses on several GLFs; (iii) processes of GLF debris transfer, focusing on mapping and interpreting boulder trains on one GLF located in Protonilus Mensae, the analysis of which suggests a best-estimate mean GLF flow speed of 7.5 mm a −1 ; and (iv) GLF hydrology, focusing on supra-GLF gulley networks. On the basis of this information, we summarize the current state of knowledge of the glaciology of martian GLFs and identify future research avenues.

Description: Strategy of valid 14 C dates choice in syngenetic permafrost The Cryosphere Discussions, 8, 5589-5621, 2014 Author(s): Y. K. Vasil'chuk and A. C. Vasil'chuk The main problem of radiocarbon dating within permafrost is the uncertain reliability of the 14 C dates. Syngenetic sediments contain allochthonous organic deposit that originated at a distance from its present position. Due to the very good preservation of organic materials in permafrost conditions and numerous re-burials of the fossils from ancient deposits into younger ones the dates could be both younger and older than the true age of dated material. The strategy for the most authentic radiocarbon date selection for dating of syncryogenic sediments is considered taking into account the fluvial origin of the syngenetic sediments. The re-deposition of organic material is discussed in terms of cyclic syncryogenic sedimentation and also the possible re-deposition of organic material in subaerial-subaqueous conditions. The advantages and the complications of dating organic micro-inclusions from ice wedges by the accelerator mass spectrometry (AMS) method are discussed applying to true age of dated material search. Radiocarbon dates of different organic materials from the same samples are compared. The younger age of the yedoma from cross-sections of Duvanny Yar in Kolyma River and Mamontova Khayata in the mouth of Lena River is substantiated due to the principle of the choice of the youngest 14 C date from the set.

Description: Modeling the elastic transmission of tidal stresses to great distances inland in channelized ice streams The Cryosphere, 8, 2007-2029, 2014 Author(s): J. Thompson, M. Simons, and V. C. Tsai Geodetic surveys suggest that ocean tides can modulate the motion of Antarctic ice streams, even at stations many tens of kilometers inland from the grounding line. These surveys suggest that ocean tidal stresses can perturb ice stream motion at distances about an order of magnitude farther inland than tidal flexure of the ice stream alone. Recent models exploring the role of tidal perturbations in basal shear stress are primarily one- or two-dimensional, with the impact of the ice stream margins either ignored or parameterized. Here, we use two- and three-dimensional finite-element modeling to investigate transmission of tidal stresses in ice streams and the impact of considering more realistic, three-dimensional ice stream geometries. Using Rutford Ice Stream as a real-world comparison, we demonstrate that the assumption that elastic tidal stresses in ice streams propagate large distances inland fails for channelized glaciers due to an intrinsic, exponential decay in the stress caused by resistance at the ice stream margins. This behavior is independent of basal conditions beneath the ice stream and cannot be fit to observations using either elastic or nonlinear viscoelastic rheologies without nearly complete decoupling of the ice stream from its lateral margins. Our results suggest that a mechanism external to the ice stream is necessary to explain the tidal modulation of stresses far upstream of the grounding line for narrow ice streams. We propose a hydrologic model based on time-dependent variability in till strength to explain transmission of tidal stresses inland of the grounding line. This conceptual model can reproduce observations from Rutford Ice Stream.

Description: Simulating the Antarctic ice sheet in the Late-Pliocene warm period: PLISMIP-ANT, an ice-sheet model intercomparison project The Cryosphere Discussions, 8, 5539-5588, 2014 Author(s): B. de Boer, A. M. Dolan, J. Bernales, E. Gasson, H. Goelzer, N. R. Golledge, J. Sutter, P. Huybrechts, G. Lohmann, I. Rogozhina, A. Abe-Ouchi, F. Saito, and R. S. W. van de Wal In the context of future climate change, understanding the nature and behaviour of ice sheets during warm intervals in Earth history is of fundamental importance. The Late-Pliocene warm period (also known as the PRISM interval: 3.264 to 3.025 million years before present) can serve as a potential analogue for projected future climates. Although Pliocene ice locations and extents are still poorly constrained, a significant contribution to sea-level rise should be expected from both the Greenland ice sheet and the West and East Antarctic ice sheets based on palaeo sea-level reconstructions. Here, we present results from simulations of the Antarctic ice sheet by means of an international Pliocene Ice Sheet Modeling Intercomparison Project (PLISMIP-ANT). For the experiments, ice-sheet models including the shallow ice and shelf approximations have been used to simulate the complete Antarctic domain (including grounded and floating ice). We compare the performance of six existing numerical ice-sheet models in simulating modern control and Pliocene ice sheets by a suite of four sensitivity experiments. Ice-sheet model forcing fields are taken from the HadCM3 atmosphere–ocean climate model runs for the pre-industrial and the Pliocene. We include an overview of the different ice-sheet models used and how specific model configurations influence the resulting Pliocene Antarctic ice sheet. The six ice-sheet models simulate a comparable present-day ice sheet, although the models are setup with their own parameter settings. For the Pliocene simulations using the Bedmap1 bedrock topography, some models show a small retreat of the East Antarctic ice sheet, which is thought to have happened during the Pliocene for the Wilkes and Aurora basins. This can be ascribed to either the surface mass balance, as the HadCM3 Pliocene climate shows a significant increase over the Wilkes and Aurora basin, or the initial bedrock topography. For the latter, our simulations with the recently published Bedmap2 bedrock topography indicate a significantly larger contribution to Pliocene sea-level rise from the East Antarctic ice sheet for all six models relative to the simulations with Bedmap1.

Description: Fluctuations of a Greenlandic tidewater glacier driven by changes in atmospheric forcing: observations and modelling of Kangiata Nunaata Sermia, 1859–present The Cryosphere, 8, 2031-2045, 2014 Author(s): J. M. Lea, D. W. F. Mair, F. M. Nick, B. R. Rea, D. van As, M. Morlighem, P. W. Nienow, and A. Weidick Many tidewater glaciers in Greenland are known to have undergone significant retreat during the last century following their Little Ice Age maxima. Where it is possible to reconstruct glacier change over this period, they provide excellent records for comparison to climate records, as well as calibration/validation for numerical models. These glacier change records therefore allow for tests of numerical models that seek to simulate tidewater glacier behaviour over multi-decadal to centennial timescales. Here we present a detailed record of behaviour from Kangiata Nunaata Sermia (KNS), SW Greenland, between 1859 and 2012, and compare it against available oceanographic and atmospheric temperature data between 1871 and 2012. We also use these records to evaluate the ability of a well-established one-dimensional flow-band model to replicate behaviour for the observation period. The record of terminus change demonstrates that KNS has advanced/retreated in phase with atmosphere and ocean climate anomalies averaged over multi-annual to decadal timescales. Results from an ensemble of model runs demonstrate that observed dynamics can be replicated. Model runs that provide a reasonable match to observations always require a significant atmospheric forcing component, but do not necessarily require an oceanic forcing component. Although the importance of oceanic forcing cannot be discounted, these results demonstrate that changes in atmospheric forcing are likely to be a primary driver of the terminus fluctuations of KNS from 1859 to 2012. We propose that the detail and length of the record presented makes KNS an ideal site for model validation exercises investigating links between climate, calving rates, and tidewater glacier dynamics.

Description: Hydrogeology of an Alpine rockfall aquifer system and its role in flood attenuation and maintaining baseflow Hydrology and Earth System Sciences, 18, 4437-4452, 2014 Author(s): U. Lauber, P. Kotyla, D. Morche, and N. Goldscheider The frequency and intensity of extreme hydrological events in Alpine regions is projected to increase with climate change. The goal of this study is to better understand the functioning of aquifers composed of complex alluvial and rockfall deposits in Alpine valleys and to quantify the role of these natural storage spaces in flood attenuation and baseflow maintenance. Geomorphological and hydrogeological mapping, tracer tests, and continuous flow measurements were conducted in the Reintal (German Alps), where runoff from a karst spring infiltrates a series of postglacial alluvial/rockfall aquifers. During high-flow conditions, groundwater velocities of 30 m h −1 were determined along 500 m; hydrograph analyses revealed short lag times (5 h) between discharge peaks upstream and downstream from the aquifer series; the maximum discharge ratio downstream (22) and the peak recession coefficient (0.196 d −1 ) are low compared with other Alpine catchments. During low-flow conditions, the underground flow path length increased to 2 km and groundwater velocities decreased to 13 m h −1 . Downstream hydrographs revealed a delayed discharge response after 101 h and peaks damped by a factor of 1.5. These results indicate that alluvial/rockfall aquifers might play an important role in the flow regime and attenuation of floods in Alpine regions.

Description: Observed groundwater temperature response to recent climate change Hydrology and Earth System Sciences, 18, 4453-4466, 2014 Author(s): K. Menberg, P. Blum, B. L. Kurylyk, and P. Bayer Climate change is known to have a considerable influence on many components of the hydrological cycle. Yet, the implications for groundwater temperature, as an important driver for groundwater quality, thermal use and storage, are not yet comprehensively understood. Furthermore, few studies have examined the implications of climate-change-induced groundwater temperature rise for groundwater-dependent ecosystems. Here, we examine the coupling of atmospheric and groundwater warming by employing stochastic and deterministic models. Firstly, several decades of temperature time series are statistically analyzed with regard to climate regime shifts (CRSs) in the long-term mean. The observed increases in shallow groundwater temperatures can be associated with preceding positive shifts in regional surface air temperatures, which are in turn linked to global air temperature changes. The temperature data are also analyzed with an analytical solution to the conduction–advection heat transfer equation to investigate how subsurface heat transfer processes control the propagation of the surface temperature signals into the subsurface. In three of the four monitoring wells, the predicted groundwater temperature increases driven by the regime shifts at the surface boundary condition generally concur with the observed groundwater temperature trends. Due to complex interactions at the ground surface and the heat capacity of the unsaturated zone, the thermal signals from distinct changes in air temperature are damped and delayed in the subsurface, causing a more gradual increase in groundwater temperatures. These signals can have a significant impact on large-scale groundwater temperatures in shallow and economically important aquifers. These findings demonstrate that shallow groundwater temperatures have responded rapidly to recent climate change and thus provide insight into the vulnerability of aquifers and groundwater-dependent ecosystems to future climate change.

Description: Snow mass decrease in the Northern Hemisphere (1979/80–2010/11) The Cryosphere Discussions, 8, 5623-5644, 2014 Author(s): Z. Li, J. Liu, L. Huang, N. Wang, B. Tian, J. Zhou, Q. Chen, and P. Zhang Snow cover has a key effect on climate change and hydrological cycling, as well as water supply to a sixth of the world's population across the Northern Hemisphere. However, reliable data on trends in snow cover in the Northern Hemisphere is lacking. Snow water equivalent (SWE) is a common measure of the amount of equivalent water of the snow pack. Here we verify the accuracy of three existing global SWE products and merge the most accurate aspects of them to generate a new SWE product covering the last 32 years (1979/80–2010/11). Using this new SWE product, we show that there has been a significant decreasing trend in the total mass of snow in the Northern Hemisphere. The most notable changes in total snow mass are −16.45 ± 6.68 and −13.55 ± 7.80 Gt year −1 in January and February, respectively. These are followed by March and December, which have trends of −12.58 ± 6.88 and −10.70 ± 5.62 Gt year −1 , respectively, from 1979/80 to 2010/11. During the same period, the temperature in the study area raised 0.17 °C decade −1 , which is thought to be the main reason of SWE decline.

Description: Factors controlling Slope Environmental Lapse Rate (SELR) of temperature in the monsoon and cold-arid glacio-hydrological regimes of the Himalaya The Cryosphere Discussions, 8, 5645-5686, 2014 Author(s): R. J. Thayyen and A. P. Dimri Moisture, temperature and precipitation interplay forced through the orographic processes sustains the Himalayan cryospheric system. However, factors controlling the Slope Environmental Lapse Rate (SELR) of temperature along the higher Himalayan mountain slopes across various glacio-hydrologic regimes remain as a key knowledge gap. Present study dwells on the orographic processes driving the moisture–temperature interplay in the monsoon and cold-arid glacio-hydrological regimes of the Himalaya. Systematic data collection at three altitudes between 2540 and 3763 m a.s.l. in the Garhwal Himalaya (hereafter called monsoon regime) and between 3500 and 5600 m a.s.l. in the Ladakh Himalaya (herefater called cold-arid regime) revealed moistrue control on temperature distribution at temporal and spatial scales. Observed daily SELR of temperature ranges between 9.0 to 1.9 °C km −1 and 17.0 to 2.8 °C km −1 in the monsoon and cold-arid regimes respectively highlighting strong regional variability. Moisture influx to the region, either from Indian summer monsoon (ISM) or from Indian winter monsoon (IWM) forced lowering of SELR. This phenophena of "monsoon lowering" of SELR is due to the release latent heat of condensation from orographically focred lifted air parcel. Seasonal response of SELR in the monsoon regime is found to be closly linked with the variations in the local lifting condensation levels (LCL). Contrary to this, cold-arid system is characterised by the extremely high values of daily SELR upto 17 °C km −1 signifying the extremely arid conditions prevailing in summer. Distinctly lower SELR devoid of monsoon lowering at higher altitude sections of monsoon and cold-arid regimes suggests sustained wetter high altitude regimes. We have proposed a SELR model for both glacio-hydrological regimes demostrating with two sections each using a derivative of the Clausius–Clapeyron relationship by deriving monthly SELR indices. It has been proposed that the manifestations of presence or absence of moisture is the single most important factor determining the temperature distribution along the higher Himalayan slopes driven by the orographic forcings. This work also suggests that the arbitary use of temperature lapse rate to extrapolate temperature to the higher Himalaya is extremely untenable.

Description: Evaluation of the satellite-based Global Flood Detection System for measuring river discharge: influence of local factors Hydrology and Earth System Sciences, 18, 4467-4484, 2014 Author(s): B. Revilla-Romero, J. Thielen, P. Salamon, T. De Groeve, and G. R. Brakenridge One of the main challenges for global hydrological modelling is the limited availability of observational data for calibration and model verification. This is particularly the case for real-time applications. This problem could potentially be overcome if discharge measurements based on satellite data were sufficiently accurate to substitute for ground-based measurements. The aim of this study is to test the potentials and constraints of the remote sensing signal of the Global Flood Detection System for converting the flood detection signal into river discharge values. The study uses data for 322 river measurement locations in Africa, Asia, Europe, North America and South America. Satellite discharge measurements were calibrated for these sites and a validation analysis with in situ discharge was performed. The locations with very good performance will be used in a future project where satellite discharge measurements are obtained on a daily basis to fill the gaps where real-time ground observations are not available. These include several international river locations in Africa: the Niger, Volta and Zambezi rivers. Analysis of the potential factors affecting the satellite signal was based on a classification decision tree (random forest) and showed that mean discharge, climatic region, land cover and upstream catchment area are the dominant variables which determine good or poor performance of the measure\-ment sites. In general terms, higher skill scores were obtained for locations with one or more of the following characteristics: a river width higher than 1km; a large floodplain area and in flooded forest, a potential flooded area greater than 40%; sparse vegetation, croplands or grasslands and closed to open and open forest; leaf area index 〉 2; tropical climatic area; and without hydraulic infrastructures. Also, locations where river ice cover is seasonally present obtained higher skill scores. This work provides guidance on the best locations and limitations for estimating discharge values from these daily satellite signals.

Description: Sampling frequency trade-offs in the assessment of mean transit times of tropical montane catchment waters under semi-steady-state conditions Hydrology and Earth System Sciences Discussions, 11, 12443-12488, 2014 Author(s): E. Timbe, D. Windhorst, R. Celleri, L. Timbe, P. Crespo, H.-G. Frede, J. Feyen, and L. Breuer Stream and soil waters were collected on a weekly basis in a tropical montane cloud forest catchment for two years and analyzed for stable water isotopes in order to infer transit time distribution functions and to define the mean transit times. Depending on the water type (stream or soil water), lumped distribution functions such as Exponential-Piston flow, Linear-Piston flow and Gamma models using temporal isotopic variations of precipitation event samples as input, were fitted. Samples were aggregated to daily, weekly, biweekly, monthly and bimonthly time scales in order to check the sensitivity of temporal sampling on model predictions. The study reveals that the effect of decreasing sampling frequency depends on the water type. For soil waters with transit times in the order of weeks to months, there was a clear trend of over prediction. In contrast, the trend of prediction for stream waters, with a dampened isotopic signal and mean transit times in the order of 2 to 4 years, was less clear and depending on the type of model used. The trade-off to coarse data resolutions could potentially lead to misleading conclusions on how water actually moves through the catchment, while at the same time predictions can reach better fitting efficiencies, lesser uncertainties, errors and biases. For both water types an optimal sampling frequency seems to be one or at most two weeks. The results of our analyses provide information for the planning (in particular in terms of cost-benefit and time requirements) of future fieldwork in similar Andean or other catchments.

Description: A model study of Abrahamsenbreen, a surging glacier in northern Spitsbergen The Cryosphere Discussions, 8, 5687-5726, 2014 Author(s): J. Oerlemans and W. J. J. van Pelt The climate sensitivity of Abrahamsenbreen, a 20 km long surge-type glacier in northern Spitsbergen, is studied with a simple glacier model. A scheme to describe the surges is included, which makes it possible to account for the effect of surges on the total mass budget of the glacier. A climate reconstruction back to AD 1300, based on ice-core data from Lomonosovfonna and climate records from Longyearbyen, is used to drive the model. The model is calibrated by requesting that it produces the correct Little Ice Age maximum glacier length and simulates the observed magnitude of the 1978-surge. Abrahamsenbreen is strongly out of balance with the current climate. If climatic conditions will remain as they were for the period 1989–2010, the glacier will ultimately shrink to a length of about 4 km (but this will take hundreds of years). For a climate change scenario involving a 2 m yr −1 rise of the equilibrium line from now onwards, we predict that in the year 2100 Abrahamsenbreen will be about 12 km long. The main effect of a surge is to lower the mean surface elevation and to increase the ablation area, thereby causing a negative perturbation of the mass budget. We found that the occurrence of surges leads to a somewhat stronger retreat of the glacier in a warming climate. Because of the very small bed slope, Abrahamsenbreen is sensitive to small perturbations in the equilibrium-line altitude E . For a decrease of E of only 160 m, the glacier would steadily grow into the Woodfjorddalen until after 2000 years it would reach the Woodfjord and calving could slow down the advance.

Description: Technical Note: On the Matt–Shuttleworth approach to estimate crop water requirements Hydrology and Earth System Sciences, 18, 4341-4348, 2014 Author(s): J. P. Lhomme, N. Boudhina, and M. M. Masmoudi The Matt–Shuttleworth method provides a way to make a one-step estimate of crop water requirements with the Penman–Monteith equation by translating the crop coefficients, commonly available in United Nations Food and Agriculture Organization (FAO) publications, into equivalent surface resistances. The methodology is based upon the theoretical relationship linking crop surface resistance to a crop coefficient and involves the simplifying assumption that the reference crop evapotranspiration (ET 0 ) is equal to the Priestley–Taylor estimate with a fixed coefficient of 1.26. This assumption, used to eliminate the dependence of surface resistance on certain weather variables, is questionable; numerical simulations show that it can lead to substantial differences between the true value of surface resistance and its estimate. Consequently, the basic relationship between surface resistance and crop coefficient, without any assumption, appears to be more appropriate for inferring crop surface resistance, despite the interference of weather variables.

Description: Variational assimilation of remotely sensed flood extents using a 2-D flood model Hydrology and Earth System Sciences, 18, 4325-4339, 2014 Author(s): X. Lai, Q. Liang, H. Yesou, and S. Daillet A variational data assimilation (4D-Var) method is proposed to directly assimilate flood extents into a 2-D dynamic flood model to explore a novel way of utilizing the rich source of remotely sensed data available from satellite imagery for better analyzing or predicting flood routing processes. For this purpose, a new cost function is specially defined to effectively fuse the hydraulic information that is implicitly indicated in flood extents. The potential of using remotely sensed flood extents for improving the analysis of flood routing processes is demonstrated by applying the present new data assimilation approach to both idealized and realistic numerical experiments.

Description: Quantitative high-resolution observations of soil water dynamics in a complicated architecture with time-lapse Ground-Penetrating Radar Hydrology and Earth System Sciences Discussions, 11, 12365-12404, 2014 Author(s): P. Klenk, S. Jaumann, and K. Roth High-resolution time-lapse Ground-Penetrating Radar (GPR) observations of advancing and retreating water tables can yield a wealth of information about near-surface water content dynamics. In this study, we present and analyze a series of imbibition, drainage and infiltration experiments which have been carried out at our artificial ASSESS test site and observed with surface based GPR. The test site features a complicated but known subsurface architecture constructed with three different kinds of sand. It allows studying soil water dynamics with GPR under a wide range of different conditions. Here, we assess in particular (i) the accurate determination of soil water dynamics averaged over the whole vertical extent by evaluating the bottom reflection and (ii) the feasibility of monitoring the dynamic shape of the capillary fringe reflection. The phenomenology of the GPR response of a dynamically changing capillary fringe is developed from a soil physical point of view. We then explain experimentally observed phenomena based on numerical simulations of both the water content dynamics and the expected GPR response.

Description: Translating aboveground cosmic-ray neutron intensity to high-frequency soil moisture profiles at sub-kilometer scale Hydrology and Earth System Sciences, 18, 4363-4379, 2014 Author(s): R. Rosolem, T. Hoar, A. Arellano, J. L. Anderson, W. J. Shuttleworth, X. Zeng, and T. E. Franz Above-ground cosmic-ray neutron measurements provide an opportunity to infer soil moisture at the sub-kilometer scale. Initial efforts to assimilate those measurements have shown promise. This study expands such analysis by investigating (1) how the information from aboveground cosmic-ray neutrons can constrain the soil moisture at distinct depths simulated by a land surface model, and (2) how changes in data availability (in terms of retrieval frequency) impact the dynamics of simulated soil moisture profiles. We employ ensemble data assimilation techniques in a "nearly-identical twin" experiment applied at semi-arid shrubland, rainfed agricultural field, and mixed forest biomes in the USA. The performance of the Noah land surface model is compared with and without assimilation of observations at hourly intervals, as well as every 2 days. Synthetic observations of aboveground cosmic-ray neutrons better constrain the soil moisture simulated by Noah in root-zone soil layers (0–100cm), despite the limited measurement depth of the sensor (estimated to be 12–20cm). The ability of Noah to reproduce a "true" soil moisture profile is remarkably good, regardless of the frequency of observations at the semi-arid site. However, soil moisture profiles are better constrained when assimilating synthetic cosmic-ray neutron observations hourly rather than every 2 days at the cropland and mixed forest sites. This indicates potential benefits for hydrometeorological modeling when soil moisture measurements are available at a relatively high frequency. Moreover, differences in summertime meteorological forcing between the semi-arid site and the other two sites may indicate a possible controlling factor to soil moisture dynamics in addition to differences in soil and vegetation properties.

Description: Large-scale 3-D modeling by integration of resistivity models and borehole data through inversion Hydrology and Earth System Sciences, 18, 4349-4362, 2014 Author(s): N. Foged, P. A. Marker, A. V. Christansen, P. Bauer-Gottwein, F. Jørgensen, A.-S. Høyer, and E. Auken We present an automatic method for parameterization of a 3-D model of the subsurface, integrating lithological information from boreholes with resistivity models through an inverse optimization, with the objective of further detailing of geological models, or as direct input into groundwater models. The parameter of interest is the clay fraction, expressed as the relative length of clay units in a depth interval. The clay fraction is obtained from lithological logs and the clay fraction from the resistivity is obtained by establishing a simple petrophysical relationship, a translator function, between resistivity and the clay fraction. Through inversion we use the lithological data and the resistivity data to determine the optimum spatially distributed translator function. Applying the translator function we get a 3-D clay fraction model, which holds information from the resistivity data set and the borehole data set in one variable. Finally, we use k -means clustering to generate a 3-D model of the subsurface structures. We apply the procedure to the Norsminde survey in Denmark, integrating approximately 700 boreholes and more than 100 000 resistivity models from an airborne survey in the parameterization of the 3-D model covering 156 km 2 . The final five-cluster 3-D model differentiates between clay materials and different high-resistivity materials from information held in the resistivity model and borehole observations, respectively.

Description: Simulating long-term past changes in the balance between water demand and availability and assessing their main drivers at the river basin management scale Hydrology and Earth System Sciences Discussions, 11, 12315-12364, 2014 Author(s): J. Fabre, D. Ruelland, A. Dezetter, and B. Grouillet The aim of this study was to assess the balance between water demand and availability and its spatial and temporal variability from 1971 to 2009 in the Herault (2500 km 2 , France) and the Ebro (85 000 km 2 , Spain) catchments. Natural streamflow was evaluated using a conceptual hydrological model. The regulation of river flow was accounted for through a widely applicable demand-driven reservoir management model applied to the largest dam in the Herault basin and to 11 major dams in the Ebro basin. Urban water demand was estimated from population and monthly unit water consumption data. Water demand for irrigation was computed from irrigated area, crop and soil data, and climatic forcing. Finally, a series of indicators comparing water supply and water demand at strategic resource and demand nodes were computed at a 10 day time step. Variations in water stress in each catchment over the past 40 years were successfully modeled, taking into account climatic and anthropogenic pressures and changes in water management strategies over time. Observed changes in discharge were explained by separating human and hydro-climatic pressures on water resources: respectively 20 and 3% of the decrease in the Ebro and the Herault discharges were linked to human-induced changes. Although key areas of the Herault basin were shown to be highly sensitive to hydro-climatic variability, the balance between water uses and availability in the Ebro basin appears to be more critical, owing to high agricultural pressure on water resources. The proposed modeling framework is currently being used to assess water stress under climatic and socio-economic prospective scenarios. Further research will investigate the effectiveness of adaptation policies aimed at maintaining the balance between water use and availability.

Description: Time-series analysis of the long-term hydrologic impacts of afforestation in the Águeda watershed of North-Central Portugal Hydrology and Earth System Sciences Discussions, 11, 12223-12256, 2014 Author(s): D. Hawtree, J. P. Nunes, J. J. Keizer, R. Jacinto, J. Santos, M. E. Rial-Rivas, A.-K. Boulet, F. Tavares-Wahren, and K.-H. Feger The north-central region of Portugal has undergone significant afforestation of the species Pinus pinaster and Eucalyptus globulus since the early 1900s; however, the long-term hydrologic impacts of this land cover change are not fully understood. To contribute to a better understanding of the potential hydrologic impacts of this land cover change, this study examines the temporal trends in 7 years of data from the Águeda watershed (part of the Vouga Basin) over the period of 1936 to 2010. Meteorological and hydrological records were analysed using a combined Thiel–Sen/Mann–Kendall trend testing approach, to assess the magnitude and significance of patterns in the observed data. These trend tests indicated that there had been no significant reduction in streamflow yield over either the entire test period, or during sub-record periods, despite the large-scale afforestation which had taken place. This lack of change is attributed to both the characteristics of the watershed and the nature of the land cover change. By contrast, a number of significant trends were found for baseflow index, which showed positive trends in the early data record (primarily during Pinus pinaster afforestation), followed by a reversal to negative trends later in the data record (primarily during Eucalyptus globulus afforestation). These changes are attributed to vegetation impacts on streamflow generating processes, both due to the species differences and to alterations in soil properties (i.e. promoting water repellency of the topsoil). These results highlight the importance of considering both vegetation types/dynamics and watershed characteristic when assessing hydrologic impacts, in particular with respect to soil properties.

Description: Reply to D. L. Peters' comment on "Streamflow input to Lake Athabasca, Canada" by Rasouli et al. (2013) Hydrology and Earth System Sciences Discussions, 11, 12257-12270, 2014 Author(s): K. Rasouli, M. A. Hernández-Henríquez, and S. J. Déry This paper provides a reply to a comment from Peters (2014) on our recent effort focused on evaluating changes in streamflow input to Lake Athabasca, Canada. Lake Athabasca experienced a 21.2% decline in streamflow input between 1960 and 2010 that has led to a marked decline in its water levels in recent decades. A reassessment of trends in naturalized Lake Athabasca water levels shows insignificant changes from our previous findings reported in Rasouli et al. (2013), and hence our previous conclusions remain unchanged. The reply closes with recommendations for future research to minimize uncertainties in historical assessments of trends in Lake Athabasca water levels and to better project its future water levels driven by climate change and anthropogenic activities in the Athabasca Lake Basin.

Description: Is sinuosity a function of slope and bankfull discharge? – A case study of the meandering rivers in the Pannonian Basin Hydrology and Earth System Sciences Discussions, 11, 12271-12290, 2014 Author(s): J. Petrovszki, G. Timár, and G. Molnár Pre-regulation channel sinuosities of the meandering rivers of the Pannonian Basin are analysed in order to define a mathematical model to estimate the influence of the bankfull discharge and the channel slope on them. As a primary database, data triplets of slope, discharge and sinuosity values were extracted from historical and modern datasets and pre-regulation historical topographic maps. Channel slope values were systematically modified to estimate figures valid before the river regulation works. The bankfull discharges were estimated from the average discharges using a robust yet complex method. The "classical" graphs of Leopold and Wolman (1957), Ackers and Charlton (1970b) and Schumm and Khan (1972) were compiled to a set up a theoretical surface, whose parameters are estimated by the real values of the above database, containing characteristics of the Pannonian Basin rivers. As a result it occurred that there is a two-dimensional function of the bankfull discharges, which provides a good estimation of the most probable sinuosity values of the rivers with the given slope and discharge characteristics. The average RMS error of this estimation is around 15% on this dataset and believed to be the effect of the non-analysed changes in the sediment discharge and size distribution.

Description: Technical Note: Field experiences using UV/VIS sensors for high-resolution monitoring of nitrate in groundwater Hydrology and Earth System Sciences Discussions, 11, 12291-12314, 2014 Author(s): M. Huebsch, F. Grimmeisen, M. Zemann, O. Fenton, K. G. Richards, P. Jordan, A. Sawarieh, P. Blum, and N. Goldscheider Two different in-situ spectrophotometers are compared that were used in the field to determine nitrate-nitrogen (NO 3 -N) concentrations at two distinct spring discharge sites. One sensor was a double wavelength spectrophotometer (DWS) and the other a multiple wavelength spectrophotometer (MWS). The objective of the study was to review the hardware options, determine ease of calibration, accuracy, influence of additional substances and to assess positive and negative aspects of the two sensors as well as troubleshooting and trade-offs. Both sensors are sufficient to monitor highly time-resolved NO 3 -N concentrations in emergent groundwater. However, the chosen path length of the sensors had a significant influence on the sensitivity and the range of detectable NO 3 -N. The accuracy of the calculated NO 3 -N concentrations of the sensors can be affected, if the content of additional substances such as turbidity, organic matter, nitrite or hydrogen carbonate significantly varies after the sensors have been calibrated to a particular water matrix. The MWS offers more possibilities for calibration and error detection, but requires more expertise compared with the DWS.

Description: A catalog of high-impact windstorms in Switzerland since 1859 Natural Hazards and Earth System Science, 14, 2867-2882, 2014 Author(s): P. Stucki, S. Brönnimann, O. Martius, C. Welker, M. Imhof, N. von Wattenwyl, and N. Philipp In recent decades, extremely hazardous windstorms have caused enormous losses to buildings, infrastructure and forests in Switzerland. This has increased societal and scientific interest in the intensity and frequency of historical high-impact storms. However, high-resolution wind data and damage statistics mostly span recent decades only. For this study, we collected quantitative (e.g., volumes of windfall timber, losses relating to buildings) and descriptive (e.g., forestry or insurance reports) information on the impact of historical windstorms. To define windstorm severity, normalized and declustered quantitative data were processed by extreme value statistics. Descriptive information was classified using a conceptual guideline. Validation with independent damage information, as well as comparison with wind measurements and a reanalysis, indicates that the most hazardous winter storms are captured, while too few moderate windstorms are detected. Strong storms in the wind measurements and reanalysis are thus added to the catalog. The final catalog encompasses approximately 240 high-impact windstorms in Switzerland since 1859. It features three robust severity classes and contains eight extreme windstorms. Evidence of high winter storm activity in the early and late 20th century compared to the mid-20th century in both damage and wind data indicates a co-variability of hazard and related damage on decadal timescales.

Description: Generalized drought assessment in Dongliao river basin based on water resources system Natural Hazards and Earth System Sciences Discussions, 2, 6703-6746, 2014 Author(s): B. S. Weng, D. H. Yan, H. Wang, T. L. Qin, and J. Yin Drought is firstly a resource issue, and with its development it transforms into a disaster issue. The occurrences of drought events usually feature determinacy and randomness. Drought issue has become one of the major factors to affect sustainable economic and social development. In this paper, we propose the generalized drought assessment index (GDAI) based on water resources system for assessing drought events. The GDAI considers water supply and water demand using a distributed hydrological model. We demonstrate the use of the proposed index in the Dongliao river basin (DRB) in the northeast China. The results simulated by the GDAI are then compared to observed drought disaster records in DRB. As second, the temporal distribution of drought events and the spatial distribution of drought frequency from the GDAI are compared with the traditional approach (i.e. the SPI, the PDSI, and the RWD). Then, generalized drought times (GDT), generalized drought duration (GDD), and generalized drought severity (GDS) were calculated by theory of runs. Application of the GDT, the GDD, and the GDS of various drought levels (i.e. mild drought, moderate drought, severe drought, and extreme drought) to the period 1960–2010 shows that the centers of gravity of them are all distributed in the middle reached of DRB, and change with time. The proposed methodology helps water managers in water-stressed regions to quantify the impact of drought, consequently, to make decisions regarding coping with drought issue.

Description: A GIS-based model to estimate flood consequences and the degree of accessibility and operability of strategic emergency response structures in urban areas Natural Hazards and Earth System Science, 14, 2847-2865, 2014 Author(s): R. Albano, A. Sole, J. Adamowski, and L. Mancusi Efficient decision-making regarding flood risk reduction has become a priority for authorities and stakeholders in many European countries. Risk analysis methods and techniques are a useful tool for evaluating costs and benefits of possible interventions. Within this context, a methodology to estimate flood consequences was developed in this paper that is based on GIS, and integrated with a model that estimates the degree of accessibility and operability of strategic emergency response structures in an urban area. The majority of the currently available approaches do not properly analyse road network connections and dependencies within systems, and as such a loss of roads could cause significant damages and problems to emergency services in cases of flooding. The proposed model is unique in that it provides a maximum-impact estimation of flood consequences on the basis of the operability of the strategic emergency structures in an urban area, their accessibility, and connection within the urban system of a city (i.e. connection between aid centres and buildings at risk), in the emergency phase. The results of a case study in the Puglia region in southern Italy are described to illustrate the practical applications of this newly proposed approach. The main advantage of the proposed approach is that it allows for defining a hierarchy between different infrastructure in the urban area through the identification of particular components whose operation and efficiency are critical for emergency management. This information can be used by decision-makers to prioritize risk reduction interventions in flood emergencies in urban areas, given limited financial resources.