ALBERT

Description: Imperfect scaling in distributions of radar-derived rainfall fields Hydrology and Earth System Sciences Discussions, 10, 11385-11422, 2013 Author(s): M. J. van den Berg, L. Delobbe, and N. E. C. Verhoest Fine scale rainfall observations for modeling exercises are often not available, but rather coarser data derived from a variety of sources are used. Effectively using these data sources in models often requires the probability distribution of the data at the applicable scale. Although numerous models for scaling distributions exist, these are often based on theoretical developments, rather than on data. In this study, we develop a model based on the α-stable distribution of rainfall fields, and tested on 5 min radar data from a Belgian weather radar. We use these data to estimate functions that describe parameters of the distribution over various scales. Moreover, we study how the mean of the distribution and the intermittency change with scale, and validate and design functions to describe the shape parameter of the distribution. This information was combined into an effective model of the distribution. Finally, the model was fitted to data from numerous storms, and the resulting parameters were compared to investigate the change in scaling behavior through time.

Description: Spatially resolved information on karst conduit flow from in-cave dye-tracing Hydrology and Earth System Sciences Discussions, 10, 11311-11335, 2013 Author(s): U. Lauber, W. Ufrecht, and N. Goldscheider Artificial tracers are powerful tools to investigate karst systems. Tracers are commonly injected into sinking streams or dolines, while springs serve as monitoring sites. The obtained flow and transport parameters represent mixed information from the vadose, epiphreatic and phreatic zones, i.e., the aquifer remains a black box. Accessible active caves constitute valuable but underexploited natural laboratories to gain detailed insights into the hydrologic functioning of the aquifer. Two multi-tracer tests in the catchment of a major karst spring (Blautopf, Germany) with injections and monitoring in two associated water caves aimed at obtaining spatially and temporally resolved information on groundwater flow in different compartments of the system. Two tracers were injected in the caves to characterize the hydraulic connections between them and with the spring. Two injections at the land surface, far from the spring, aimed at resolving the aquifer's internal drainage structure. Tracer breakthrough curves were monitored by field fluorimeters in caves and at the spring. Results demonstrate the dendritic drainage structure of the aquifer. It was possible to obtain relevant flow and transport parameters for different sections of this system. The highest mean flow velocities (275 m h −1 ) were observed in the near-spring epiphreatic section (open-channel flow), while velocities in the phreatic zone (pressurized flow) were one order of magnitude lower. Determined conduit water volumes confirm results of water balances and hydrograph analyses. In conclusion, experiments and monitoring in caves can deliver spatially resolved information on karst aquifer heterogeneity and dynamics that cannot be obtained by traditional investigative methods.

Description: On the lack of robustness of hydrologic models regarding water balance simulation – a diagnostic approach on 20 mountainous catchments using three models of increasing complexity Hydrology and Earth System Sciences Discussions, 10, 11337-11383, 2013 Author(s): L. Coron, V. Andréassian, C. Perrin, M. Bourqui, and F. Hendrickx This paper investigates the robustness of rainfall–runoff models when their parameters are transferred in time. More specifically, we studied their ability to simulate water balance on periods with different hydroclimatic characteristics. The testing procedure consisted in a series of parameter transfers between 10-yr periods and the systematic analysis of mean-volume errors. This procedure was applied to three conceptual models of different structural complexity over 20 mountainous catchments in southern France. The results showed that robustness problems are common. Errors on 10-yr-mean flows were significant for all three models and calibration periods, even when the entire record was used for calibration. Various graphical and numerical tools were used to show strong similarities between the shapes of mean flow biases calculated on a 10-yr-long sliding window when various parameter sets are used. Unexpected behavioural similarities were observed between the three models tested, considering their large differences in structural complexity. While the actual causes for robustness problems in these models remain unclear, this work stresses the limited transferability in time of the water balance adjustments made through parameter optimization. Although absolute differences between simulations obtained with different calibrated parameter sets were sometimes substantial, relative differences in simulated mean flows between time periods remained similar regardless of the calibrated parameter sets.

Description: Review Article: On the relation between the seismic activity and the Hurst exponent of the geomagnetic field at the time of the 2000 Izu swarm Natural Hazards and Earth System Science, 13, 2189-2194, 2013 Author(s): F. Masci and J. N. Thomas Many papers document the observation of earthquake-related precursory signatures in geomagnetic field data. However, the significance of these findings is ambiguous because the authors did not adequately take into account that these signals could have been generated by other sources, and the seismogenic origin of these signals have not been validated by comparison with independent datasets. Thus, they are not reliable examples of magnetic disturbances induced by the seismic activity. Hayakawa et al. (2004) claim that at the time of the 2000 Izu swarm the Hurst exponent of the Ultra-Low-Frequency (ULF: 0.001–10 Hz) band of the geomagnetic field varied in accord with the energy released by the seismicity. The present paper demonstrates that the behaviour of the Hurst exponent was insufficiently investigated and also misinterpreted by the authors. We clearly show that during the Izu swarm the changes of the Hurst exponent were strongly related to the level of global geomagnetic activity and not to the increase of the local seismic activity.

Description: Forest fire danger rating in complex topography – results from a case study in the Bavarian Alps in autumn 2011 Natural Hazards and Earth System Science, 13, 2157-2167, 2013 Author(s): C. Schunk, C. Wastl, M. Leuchner, C. Schuster, and A. Menzel Forest fire danger rating based on sparse meteorological stations is known to be potentially misleading when assigned to larger areas of complex topography. This case study examines several fire danger indices based on data from two meteorological stations at different elevations during a major drought period. This drought was caused by a persistent high pressure system, inducing a pronounced temperature inversion and its associated thermal belt with much warmer, dryer conditions in intermediate elevations. Thus, a massive drying of fuels, leading to higher fire danger levels, and multiple fire occurrences at mid-slope positions were contrasted by moderate fire danger especially in the valleys. The ability of fire danger indices to resolve this situation was studied based on a comparison with the actual fire danger as determined from expert observations, fire occurrences and fuel moisture measurements. The results revealed that, during temperature inversion, differences in daily cycles of meteorological parameters influence fire danger and that these are not resolved by standard meteorological stations and fire danger indices (calculated on a once-a-day basis). Additional stations in higher locations or high-resolution meteorological models combined with fire danger indices accepting at least hourly input data may allow reasonable fire danger calculations under these circumstances.

Description: Preface "Landslide hazard and risk assessment at different scales" Natural Hazards and Earth System Science, 13, 2169-2171, 2013 Author(s): P. Reichenbach, A. Günther, and T. Glade

Description: Antecedent flow conditions and nitrate concentrations in the Mississippi River Basin Hydrology and Earth System Sciences Discussions, 10, 11451-11484, 2013 Author(s): J. C. Murphy, R. M. Hirsch, and L. A. Sprague The influence of antecedent flow conditions on nitrate concentrations was explored at eight sites in the Mississippi River Basin, USA. Antecedent moisture conditions have been shown to influence nutrient export from small, relatively homogenous basins, but this influence has not been observed at a regional or continental scale. Antecedent flow conditions were quantified as the ratio between the mean daily flow of the previous year and the mean daily flow from the period of record ( Q ratio), and the Q ratio was statistically related to nitrate anomalies (the unexplained variability in nitrate concentration after filtering out season, long-term trend, and contemporaneous flow effects) at each site. Nitrate anomaly and Q ratio were negatively related at three of the four major tributary sites and upstream in the Mississippi River, indicating that when the previous year was drier than average, at these sites, nitrate concentrations were higher than expected. The strength of these relationships increased when data were subdivided by contemporaneous flow conditions. Five of the eight sites had significant negative relationships ( p ≤ 0.05) at high or moderately high contemporaneous flows, suggesting nitrate that accumulates in these basins during a drought is flushed during subsequent storm events. At half of the sites, when flow during the previous year was 50% drier than average, nitrate concentration can be from 9 and 27% higher than nitrate concentrations that follow a year with average daily flow. Conversely, nitrate concentration can be from 8 and 21% lower than expected when the previous year was 50% wetter than average. These relationships between nitrate concentration and Q ratio serve as the basis for future studies that can better define specific hydrologic processes occurring during and after a drought, which influence nitrate concentration, such as the duration or magnitude of low flows, and the timing of low and high flows.

Description: Contribution of snow and glacier melt to discharge for highly glacierised catchments in Norway Hydrology and Earth System Sciences Discussions, 10, 11485-11517, 2013 Author(s): M. Engelhardt, T. V. Schuler, and L. M. Andreassen Glacierised catchments significantly alter the streamflow regime due to snow and glacier meltwater contribution to discharge. In this study, we modelled the mass balance and discharge rates for three highly glacierised catchments (〉50% glacier cover) in western Norway over the period 1961–2012. The spatial pattern of the catchments follows a gradient in climate continentality from west to east. The model uses gridded temperature and precipitation values from seNorge ( http://senorge.no ) as input which are available at a daily resolution. It accounts for accumulation of snow, transformation of snow to firn and ice, evaporation and melt. The model was calibrated for each catchment based on measurements of seasonal glacier mass-balances and daily discharge rates. For validation, daily melt rates were compared with measurements from sonic rangers located in the ablation zones of two of the glaciers and an uncertainty analysis was performed for the third catchment. The discharge contributions from snowmelt, glacier melt and rain were analysed with respect to spatial variations and temporal evolution. The model simulations reveal an increase of the relative contribution from glacier melt for the three catchments from less than 10% in the early 1990s to 15–30% in the late 2000s. The decline in precipitation by 10–20% in the same period was therefore overcompensated resulting in an increase of the annual discharge by 5–20%. Annual discharge sums and annual glacier melt are strongest correlated with annual and winter precipitation at the most maritime glacier and, with increased climate continentality, variations in both glacier melt contribution and annual discharge are becoming stronger correlated with variations in summer temperatures.

Description: Bias correction can modify climate model-simulated precipitation changes without adverse affect on the ensemble mean Hydrology and Earth System Sciences Discussions, 10, 11585-11611, 2013 Author(s): E. P. Maurer and D. W. Pierce When applied to remove climate model biases in precipitation, quantile mapping can in some settings modify the simulated trends. This has important implications when the precipitation will be used to drive an impacts model that is sensitive to changes in precipitation. We use daily precipitation output from 12 general circulation models (GCMs) over the conterminous United States interpolated to a common 1° grid, and gridded observations aggregated to the same scale, to compare precipitation differences before and after quantile mapping bias correction. The change in seasonal mean (winter, DJF, and summer, JJA) precipitation between different 30-yr historical periods is compared to examine (1) the consensus among GCMs as to whether the bias correction tends to amplify or diminish their simulated precipitation trends, and (2) whether the modification of the change in precipitation tends to improve or degrade the correspondence to observed changes in precipitation for the same periods. In some cases, for a particular GCM, the trend modification can be as large as the original simulated change, though the areas where this occurs varies among GCMs so the ensemble median shows smaller trend modification. In specific locations and seasons the trend modification by quantile mapping improves correspondence with observed trends, and in others it degrades it. In the majority of the domain the ensemble median is for little effect on the correspondence of simulated precipitation trends with observed. This highlights the need to use an ensemble of GCMs rather than relying on a small number of models to estimate impacts.

Description: Linking ICT and society in early warning and adaptation to hydrological extremes in mountains Natural Hazards and Earth System Science, 13, 2253-2270, 2013 Author(s): C. de Jong The assessment of the societal impact of hydrological extremes is particularly important in mountain regions, since mountains can be considered both as the generators and victims of extreme events. ICT can provide a powerful tool for transmitting hydro-meteorological information to predict, prepare and adapt to such events. However, in remote regions, such as mountains, the poles, deserts and islands, preventive and adaptive measures are often restricted by data availability and lack and/or incoherence of data networks. This paper distinguishes between early warning of floods and droughts, emphasising the latter in particular in mountains and explores the possibilities of enhancing the role of society in data collection, the identification, activation and application of stakeholder knowledge and transferral of data from gauged to ungauged catchments.

Description: The environmental impact of the Puyehue–Cordon Caulle 2011 volcanic eruption on Buenos Aires Natural Hazards and Earth System Science, 13, 2319-2330, 2013 Author(s): G. B. Raga, D. Baumgardner, A. G. Ulke, M. Torres Brizuela, and B. Kucienska On 4 June 2011, the volcanic complex Puyehue–Cordon Caulle located in the Chilean Andes erupted, producing a plume of gases and particles that eventually circled the Southern Hemisphere, disrupting air travel and depositing ash in large quantities. On eight occasions, the plume passed over the city of Buenos Aires, Argentina, leading local authorities to close the two international airports. The eruption occurred during an on-going field campaign when measurements of the properties of atmospheric aerosol particles were being made in Buenos Aires as part of a year-long study of the concentration and optical properties of aerosol at one site in the city. The suite of instruments deployed in Buenos Aires were not tailored to measurements of volcanic ash, but were designed to characterize urban conditions. Nevertheless, these measurements were analysed for periods when vertical profiles of aerosol backscatter, made with a ceilometer, clearly showed the presence of the volcano plume over the research site and resulted in airport closure. Aerosol optical thickness derived from AERONET, MODIS and a ceilometer at our research site, all show enhanced values clearly indicating that the three platforms identified the volcanic plume simultaneously. However, a quantitative comparison of the different estimates proves difficult, suggesting large spatial and temporal variability of the plume. Our results indicate that the number concentration of condensation nuclei (CN), the mass concentration of particle-bound polycyclic aromatic hydrocarbons (PPAH) and the light absorption coefficient exceeded the average background values by more than one standard deviation during the events of volcanic plume. The anomalous concentrations of CN suggest new particle formation, presumably from the conversion of SO 2 , while the anomalous concentrations of PPAH may come from the uptake of PAHs on the plume particles or from chemical reactions on the surface of plume particles. The anomalous absorption coefficients indicate that plume particles may contain certain compounds that can absorb radiation at 550 nm. Another possible explanation consistent with the observations is the scavenging of black carbon from urban sources as the plume descends through the boundary layer to the surface. In addition, the volcanic plume influenced the local meteorology resulting in a decrease of the temperature when compared to the average temperature during days with no plume present.

Description: Contribution of land use changes to future flood damage along the river Meuse in the Walloon region Natural Hazards and Earth System Science, 13, 2301-2318, 2013 Author(s): A. Beckers, B. Dewals, S. Erpicum, S. Dujardin, S. Detrembleur, J. Teller, M. Pirotton, and P. Archambeau Managing flood risk in Europe is a critical issue because climate change is expected to increase flood hazard in many european countries. Beside climate change, land use evolution is also a key factor influencing future flood risk. The core contribution of this paper is a new methodology to model residential land use evolution. Based on two climate scenarios ("dry" and "wet"), the method is applied to study the evolution of flood damage by 2100 along the river Meuse. Nine urbanization scenarios were developed: three of them assume a "current trend" land use evolution, leading to a significant urban sprawl, while six others assume a dense urban development, characterized by a higher density and a higher diversity of urban functions in the urbanized areas. Using damage curves, the damage estimation was performed by combining inundation maps for the present and future 100 yr flood with present and future land use maps and specific prices. According to the dry scenario, the flood discharge is expected not to increase. In this case, land use changes increase flood damages by 1–40%, to €334–462 million in 2100. In the wet scenario, the relative increase in flood damage is 540–630%, corresponding to total damages of €2.1–2.4 billion. In this extreme scenario, the influence of climate on the overall damage is 3–8 times higher than the effect of land use change. However, for seven municipalities along the river Meuse, these two factors have a comparable influence. Consequently, in the "wet" scenario and at the level of the whole Meuse valley in the Walloon region, careful spatial planning would reduce the increase in flood damage by no more than 11–23%; but, at the level of several municipalities, more sustainable spatial planning would reduce future flood damage to a much greater degree.

Description: The effect of training image and secondary data integration with multiple-point geostatistics in groundwater modeling Hydrology and Earth System Sciences Discussions, 10, 11829-11860, 2013 Author(s): X. He, T. O. Sonnenborg, F. Jørgensen, and K. H. Jensen Multiple-point geostatistic simulation (MPS) has recently become popular in stochastic hydrogeology, primarily because of its capability to derive multivariate distributions from the training image (TI). However, its application in three dimensional simulations has been constrained by the difficulty of constructing 3-D TI. The object-based TiGenerator may be a useful tool in this regard; yet the sensitivity of model predictions to the training image has not been documented. Another issue in MPS is the integration of multiple geophysical data. The best way to retrieve and incorporate information from high resolution geophysical data is still under discussion. This work shows that TI from TiGenerator delivers acceptable results when used for groundwater modeling, although the TI directly converted from high resolution geophysical data leads to better simulation. The model results also indicate that soft conditioning in MPS is a convenient and efficient way of integrating secondary data such as 3-D airborne electromagnetic data, but over conditioning has to be avoided.

Description: Multi-decadal river flows variations in France Hydrology and Earth System Sciences Discussions, 10, 11861-11900, 2013 Author(s): J. Boé and F. Habets In this article, multi-decadal variations in French hydroclimate are investigated, with a specific focus on river flows. Based on long observed series, it is shown that river flows in France generally exhibit large multi-decadal variations on the historical period, especially in spring. Differences of means between two 21 yr periods of the 20th century as large as 40% are indeed found for many gauging stations. Multi-decadal spring river flows variations are associated with variations in spring precipitation and temperature. These multi-decadal variations in precipitation are themselves found to be driven by large-scale atmospheric circulation, more precisely by a multi-decadal oscillation in a sea level pressure dipole between western Europe and the East Atlantic. It is suggested that the Atlantic Multidecadal Variability, the main mode of decadal variability in the North Atlantic/Europe sector, controls those variations in large-scale circulation and is therefore the main ultimate driver of multi-decadal variations in spring river flows. Multi-decadal variations in river flows in other seasons, and in particular summer, are also noted. As they are not associated with significant surface climate anomalies (i.e. temperature, precipitation) in summer, other mechanisms are investigated based on hydrological simulations. The impact of climate variations in spring on summer soil moisture, and the impact of soil moisture in summer on the runoff to precipitation ratio, could potentially play a role in multi-decadal summer river flows variations. The large amplitude of the multi-decadal variations in French river flows suggests that internal variability may play a very important role in the evolution of river flows during the next decades, potentially temporarily limiting, reversing or seriously aggravating the long-term impacts of anthropogenic climate change.

Description: Estimating soil suction from electrical resistivity Natural Hazards and Earth System Science, 13, 2369-2379, 2013 Author(s): E. Piegari and R. Di Maio Soil suction and resistivity strongly depend on the degree of soil saturation and, therefore, both are used for estimating water content variations. The main difference between them is that soil suction is measured using tensiometers, which give point information, while resistivity is obtained by tomography surveys, which provide distributions of resistivity values in large volumes, although with less accuracy. In this paper, we have related soil suction to electrical resistivity with the aim of obtaining information about soil suction changes in large volumes, and not only for small areas around soil suction probes. We derived analytical relationships between soil matric suction and electrical resistivity by combining the empirical laws of van Genuchten and Archie. The obtained relationships were used to evaluate maps of soil suction values in different ashy layers originating in the explosive activity of the Mt Somma-Vesuvius volcano (southern Italy). Our findings provided a further example of the high potential of geophysical methods in contributing to more effective monitoring of soil stress conditions; this is of primary importance in areas where rainfall-induced landslides occur periodically.

Description: Teleconnection analysis of runoff and soil moisture over the Pearl River basin in South China Hydrology and Earth System Sciences Discussions, 10, 11943-11982, 2013 Author(s): J. Niu, J. Chen, and B. Sivakumar This study explores the teleconnection of two climatic patterns, namely the El Niño-Southern Oscillation (ENSO) and the Indian Ocean Dipole (IOD), with hydrological processes over the Pearl River basin in South China. The Variable Infiltration Capacity (VIC) model is used to simulate the daily hydrological processes over the basin for the study period 1952–2000, and then, using the simulation results, the time series of the monthly runoff and soil moisture anomalies for its ten sub-basins are aggregated. Wavelet analysis is performed to explore the variability properties of these time series at 49 timescales ranging from 2 months to 9 yr. Use of wavelet coherence and rank correlation method reveals that the dominant variabilities of the time series of runoff and soil moisture are basically correlated with IOD. The influences of ENSO on the terrestrial hydrological processes are mainly found in the eastern sub-basins. The teleconnections between climatic patterns and hydrological variability also serve as a reference basis for inferences on the occurrence of extreme hydrological events (e.g. floods and droughts).

Description: Sulfur dioxide emissions from Papandayan and Bromo, two Indonesian volcanoes Natural Hazards and Earth System Science, 13, 2399-2407, 2013 Author(s): P. Bani, Surono, M. Hendrasto, H. Gunawan, and S. Primulyana Indonesia hosts 79 active volcanoes, representing 14% of all active volcanoes worldwide. However, little is known about their SO 2 contribution into the atmosphere, due to isolation and access difficulties. Existing SO 2 emission budgets for the Indonesian archipelago are based on extrapolations and inferences as there is a considerable lack of field assessments of degassing. Here, we present the first SO 2 flux measurements using differential optical absorption spectroscopy (DOAS) for Papandayan and Bromo, two of the most active volcanoes in Indonesia. Results indicate mean SO 2 emission rates of 1.4 t d −1 from the fumarolic activity of Papandayan and more than 22–32 t d −1 of SO 2 released by Bromo during a declining eruptive phase. These DOAS results are very encouraging and pave the way for a better evaluation of Indonesian volcanic emissions.

Description: Hydrological functions of sinkholes and characteristics of point recharge in groundwater basins Hydrology and Earth System Sciences Discussions, 10, 11423-11449, 2013 Author(s): N. Somaratne, K. Smettem, J. Lawson, K. Nguyen, and J. Frizenschaf Karstic limestone aquifers are hydrologically and hydrochemically extremely heterogeneous and point source recharge via sinkholes and fissures is a common feature. We studied three groundwater systems in karstic settings dominated by point source recharge in order to assess the relative contributions to total recharge from point sources using chloride and δ 18 O relations. Preferential groundwater flows were observed through an inter-connected network of highly conductive zones with groundwater mixing along flow paths. Measurements of salinity and chloride indicated that fresh water pockets exist at point recharge locations. A measurable fresh water plume develops only when a large quantity of surface water enters the aquifer as a point recharge source. The difference in chloride concentrations in diffuse and point recharge zones decreases as aquifer saturated thickness increases and the plumes become diluted through mixing. The chloride concentration in point recharge fluxes crossing the watertable plane can remain at or near surface runoff chloride concentrations, rather than in equilibrium with groundwater chloride. In such circumstances the conventional chloride mass balance method that assumes equilibrium of recharge water chloride with groundwater requires modification to include both point and diffuse recharge mechanisms.

Description: Interviewing insights regarding the fatalities inflicted by the 2011 Great East Japan Earthquake Natural Hazards and Earth System Science, 13, 2173-2187, 2013 Author(s): M. Ando, M. Ishida, Y. Hayashi, C. Mizuki, Y. Nishikawa, and Y. Tu One hundred fifty survivors of the 11 March 2011 Great East Japan Earthquake (Tohoku-oki earthquake) ( M w = 9.0) were interviewed to study the causes of deaths from the associated tsunami in coastal areas of Tohoku. The first official tsunami warning underestimated the height of the tsunami and 40% of the interviewees did not obtain this warning due to immediate blackouts and a lack of communication after the earthquake. Many chose to remain in dangerous locations based on the underestimated warning and their experiences with previous smaller tsunamis and/or due to misunderstanding the mitigating effects of nearby breakwaters in blocking incoming tsunamis. Some delayed their evacuation to perform family safety checks, and in many situations, the people affected misunderstood the risks involved in tsunamis. In this area, three large tsunamis have struck in the 115 yr preceding the 2011 tsunami. These tsunamis remained in the collective memory of communities, and numerous measures against future tsunami damage, such as breakwaters and tsunami evacuation drills, had been implemented. Despite these preparedness efforts, approximately 18 500 deaths and cases of missing persons occurred. The death rate with the age of 65 and above was particularly high, four times higher than that with other age groups. These interviews indicate that deaths resulted from a variety of reasons, but if residents had taken immediate action after the major ground motion stopped, most residents might have been saved. Education about the science behind earthquakes and tsunamis could help save more lives in the future.

Description: Climate information based streamflow and rainfall forecasts for Huai River Basin using Hierarchical Bayesian Modeling Hydrology and Earth System Sciences Discussions, 10, 11559-11584, 2013 Author(s): X. Chen, Z. Hao, N. Devineni, and U. Lall A Hierarchal Bayesian model for forecasting regional summer rainfall and streamflow season-ahead using exogenous climate variables for East Central China is presented. The model provides estimates of the posterior forecasted probability distribution for 12 rainfall and 2 streamflow stations considering parameter uncertainty, and cross-site correlation. The model has a multilevel structure with regression coefficients modeled from a common multivariate normal distribution results in partial-pooling of information across multiple stations and better representation of parameter and posterior distribution uncertainty. Covariance structure of the residuals across stations is explicitly modeled. Model performance is tested under leave-10-out cross-validation. Frequentist and Bayesian performance metrics used include Receiver Operating Characteristic, Reduction of Error, Coefficient of Efficiency, Rank Probability Skill Scores, and coverage by posterior credible intervals. The ability of the model to reliably forecast regional summer rainfall and streamflow season-ahead offers potential for developing adaptive water risk management strategies.

Description: Review Article: Potential geomorphic consequences of a future great ( M w = 8.0+) Alpine Fault earthquake, South Island, New Zealand Natural Hazards and Earth System Science, 13, 2279-2299, 2013 Author(s): T. R. Robinson and T. R. H. Davies The Alpine Fault in New Zealand's South Island has not sustained a large magnitude earthquake since ca. AD 1717. The time since this rupture is close to the average inferred recurrence interval of the fault (~300 yr). The Alpine Fault is therefore expected to generate a large magnitude earthquake in the near future. Previous ruptures of this fault are inferred to have generated M w = 8.0 or greater earthquakes and to have resulted in, amongst other geomorphic hazards, large-scale landslides and landslide dams throughout the Southern Alps. There is currently 85% probability that the Alpine Fault will cause a M w = 8.0+ earthquake within the next 100 yr. While the seismic hazard is fairly well understood, that of the consequential geomorphic activity is less well studied, and these consequences are explored herein. They are expected to include landsliding, landslide damming, dam-break flooding, debris flows, river aggradation, liquefaction, and landslide-generated lake/fiord tsunami. Using evidence from previous events within New Zealand as well as analogous international examples, we develop first-order estimates of the likely magnitude and possible locations of the geomorphic effects associated with earthquakes. Landsliding is expected to affect an area 〉 30 000 km 2 and involve 〉 1billion m 3 of material. Some tens of landslide dams are expected to occur in narrow, steep-sided gorges in the affected region. Debris flows will be generated in the first long-duration rainfall after the earthquake and will continue to occur for several years as rainfall (re)mobilises landslide material. In total more than 1000 debris flows are likely to be generated at some time after the earthquake. Aggradation of up to 3 m will cover an area 〉 125 km 2 and is likely to occur on many West Coast alluvial fans and floodplains. The impact of these effects will be felt across the entire South Island and is likely to continue for several decades.

Description: A spatial bootstrap technique for parameter estimation of rainfall annual maxima distribution Hydrology and Earth System Sciences Discussions, 10, 11755-11794, 2013 Author(s): F. Uboldi, A. N. Sulis, C. Lussana, M. Cislaghi, and M. Russo Estimation of extreme event distributions and depth-duration-frequency (DDF) curves is achieved at any target site by repeated sampling among all available raingauge data in the surrounding area. The estimate is computed over a gridded domain in Northern Italy, using precipitation time series from 1929 to 2011, including data from historical analog stations and from the present-day automatic observational network. The presented local regionalisation naturally overcomes traditional station-point methods, with their demand of long historical series and their sensitivity to very rare events occurring at very few stations, possibly causing unrealistic spatial gradients in DDF relations. At the same time, the presented approach allows for spatial dependence, necessary in a geographical domain such as Lombardy, complex for both its topography and its climatology. The bootstrap technique enables evaluating uncertainty maps for all estimated parameters and for rainfall depths at assigned return periods.

Description: Overview of the first HyMeX Special Observation Period over Italy: observations and model results Hydrology and Earth System Sciences Discussions, 10, 11643-11710, 2013 Author(s): R. Ferretti, E. Pichelli, S. Gentile, I. Maiello, D. Cimini, S. Davolio, M. M. Miglietta, G. Panegrossi, L. Baldini, F. Pasi, F. S. Marzano, A. Zinzi, S. Mariani, M. Casaioli, G. Bartolini, N. Loglisci, A. Montani, C. Marsigli, A. Manzato, A. Pucillo, M. E. Ferrario, V. Colaiuda, and R. Rotunno During the first Hymex campaign (5 September–6 November 2012) referred to as Special Observation Period (SOP-1), dedicated to heavy precipitation events and flash floods in Western Mediterranean, three Italian hydro-meteorological monitoring sites were activated: Liguria-Tuscany, North-Eastern Italy and Central Italy. The extraordinary deployment of advanced instrumentation, including instrumented aircrafts, and the use of several different operational weather forecast models has allowed an unprecedented monitoring and analysis of high impact weather events around the Italian hydro-meteorological sites. This activity has seen the strict collaboration between the Italian scientific and operational communities. In this paper, an overview of the Italian organization during the SOP-1 is provided, and selected Intensive Observation Periods (IOPs) are described. A significant event for each Italian target area is chosen for this analysis: IOP2 (12–13 September 2012) in North-Eastern Italy, IOP13 (15–16 October 2012) in Central Italy and IOP19 (3–5 November 2012) in Liguria and Tuscany. For each IOP the meteorological characteristics, together with special observations and weather forecasts, are analyzed with the aim of highlighting strengths and weaknesses of the forecast modeling systems. Moreover, using one of the three events, the usefulness of different operational chains is highlighted.

Description: Correcting basin-scale snowfall in a mountainous basin using a distributed snowmelt model and remote sensing data Hydrology and Earth System Sciences Discussions, 10, 11711-11753, 2013 Author(s): M. Shrestha, L. Wang, T. Koike, H. Tsutsui, Y. Xue, and Y. Hirabayashi Adequate estimation of the spatial distribution of snowfall is critical in hydrologic modeling. However, this is a well-known problem in estimating basin-scale snowfall, especially in mountainous basins with data scarcity. This study focuses on correction and estimation of this spatial distribution, which considers topographic effects within the basin. A method is proposed that optimizes an altitude-based snowfall correction factor ( C fsnow ). This is done through multi-objective calibration of a spatially distributed, multilayer energy and water balance-based snowmelt model (WEB-DHM-S) with observed discharge and remotely sensed snow cover data from the Moderate Resolution Imaging Spectroradiometer (MODIS). The Shuffled Complex Evolution – University of Arizona automatic search algorithm is used to obtain the optimal value of C fsnow for minimum cumulative error in discharge and snow cover simulations. Discharge error is quantified by Nash–Sutcliffe efficiency and relative volume deviation, and snow cover error was estimated by pixel-by-pixel analysis. The study region is the heavily snow-fed Yagisawa Basin of the Upper Tone River in northeast Japan. First, the system was applied to one snow season (2002–2003), obtaining an optimized C fsnow of 0.0007 m −1 . For validation purposes, the optimized C fsnow was implemented to correct snowfall in 2004, 2002 and 2001. Overall, the system was effective, implying improvements in correlation of simulated vs. observed discharge and snow cover. The 4 yr mean of basin-average snowfall for the corrected spatial snowfall distribution was 1160 mm (780 mm before correction). Execution of sensitivity runs against other model input and parameters indicated that C fsnow could be affected by uncertainty in shortwave radiation and setting of the threshold air temperature parameter. Our approach is suitable to correct snowfall and estimate its distribution in poorly-gauged basins, where elevation dependence of snowfall amount is strong.

Description: Attribution of hydrologic forecast uncertainty within scalable forecast windows Hydrology and Earth System Sciences Discussions, 10, 11795-11828, 2013 Author(s): L. Yang, F. Tian, Y. Sun, X. Yuan, and H. Hu Hindcasts based on the Extended Streamflow Prediction (ESP) approach are carried out in a typical rainfall-dominated basin in China, aiming to examine the roles of initial condition (IC), future atmospheric forcing (FC) and hydrologic model uncertainty (MU) in the streamflow forecast skill. The combined effects of IC and FC are explored within the framework of a forecast window. By implementing virtual numerical simulations without the consideration of MU, it is found that the dominance of IC could last up to 90 days in dry season, while its impact gives way to FC for lead times exceeding 30 days in the wet season. The combined effects of IC and FC on the forecast skill are further investigated by proposing a dimensionless parameter ( β ) that represents the ratio of the total amount of initial water storage and the incoming rainfall. The forecast skill increases exponentially with β , and varies greatly in different forecast windows. Moreover, the influence of MU on forecast skill is examined by focusing on the uncertainty of model parameters. Two different hydrologic model calibration strategies are carried out. The results indicate that the uncertainty of model parameters exhibits a more significant influence on the forecast skill in the dry season than in the wet season. The ESP approach is more skillful in monthly streamflow forecast during the transition period from wet to dry than otherwise. For the transition period from dry to wet, the low skill of the forecasts could be attributed to the combined effects of IC and FC, but less to the biases in the hydrologic model parameters. For the forecasting in dry season, the usefulness of the ESP approach is heavily dependent on the strategy of the model calibration.

Description: Integrating spatial, temporal, and size probabilities for the annual landslide hazard maps in the Shihmen watershed, Taiwan Natural Hazards and Earth System Science, 13, 2353-2367, 2013 Author(s): C. Y. Wu and S. C. Chen Landslide spatial, temporal, and size probabilities were used to perform a landslide hazard assessment in this study. Eleven intrinsic geomorphological, and two extrinsic rainfall factors were evaluated as landslide susceptibility related factors as they related to the success rate curves, landslide ratio plots, frequency distributions of landslide and non-landslide groups, as well as probability–probability plots. Data on landslides caused by Typhoon Aere in the Shihmen watershed were selected to train the susceptibility model. The landslide area probability, based on the power law relationship between the landslide area and a noncumulative number, was analyzed using the Pearson type 5 probability density function. The exceedance probabilities of rainfall with various recurrence intervals, including 2, 5, 10, 20, 50, 100 and 200 yr, were used to determine the temporal probabilities of the events. The study was conducted in the Shihmen watershed, which has an area of 760 km 2 and is one of the main water sources for northern Taiwan. The validation result of Typhoon Krosa demonstrated that this landslide hazard model could be used to predict the landslide probabilities. The results suggested that integration of spatial, area, and exceedance probabilities to estimate the annual probability of each slope unit is feasible. The advantage of this annual landslide probability model lies in its ability to estimate the annual landslide risk, instead of a scenario-based risk.

Description: On the puzzling feature of the silence of precursory electromagnetic emissions Natural Hazards and Earth System Science, 13, 2381-2397, 2013 Author(s): K. Eftaxias, S. M. Potirakis, and T. Chelidze It has been suggested that fracture-induced MHz–kHz electromagnetic emissions (EME), which emerge from a few days up to a few hours before the main seismic shock occurrence permit a real-time monitoring of the damage process during the last stages of earthquake preparation, as it happens at the laboratory scale. Despite fairly abundant evidence, electromagnetic (EM) precursors have not been adequately accepted as credible physical phenomena. These negative views are enhanced by the fact that certain "puzzling features" are repetitively observed in candidate fracture-induced pre-seismic EME. More precisely, EM silence in all frequency bands appears before the main seismic shock occurrence, as well as during the aftershock period. Actually, the view that "acceptance of "precursive" EM signals without convincing co-seismic signals should not be expected" seems to be reasonable. In this work we focus on this point. We examine whether the aforementioned features of EM silence are really puzzling ones or, instead, reflect well-documented characteristic features of the fracture process, in terms of universal structural patterns of the fracture process, recent laboratory experiments, numerical and theoretical studies of fracture dynamics, critical phenomena, percolation theory, and micromechanics of granular materials. Our analysis shows that these features should not be considered puzzling.

Description: A conceptual model of check dam hydraulics for gully control Hydrology and Earth System Sciences Discussions, 10, 11901-11941, 2013 Author(s): C. Castillo, R. Pérez, and J. A. Gómez There is little information in scientific literature regarding the modifications induced by check dam systems in flow regimes in restored gully reaches, despite it being a crucial issue for the design of conservation measures. Here, we develop a conceptual model to classify flow regimes in straight rectangular channels for initial and dam-filling conditions as well as a method of estimating efficiency in order to provide guidelines for optimal design. The model integrates several previous mathematical approaches for assessing the main processes involved (hydraulic jump HJ, impact flow, gradually varied flows). Its performance was compared with the simulations obtained from IBER, a bi-dimensional hydrodynamic model. The impact of check dam spacing (defined by the geometric factor of influence c ) on efficiency was explored. Eleven main classifications of flow regimes were identified depending on the element and level of influence. The model produced similar results when compared with IBER, but led to higher estimations of HJ and impact lengths. Total influence guaranteed maximum efficiency and HJ control defining the location of the optimal c . Geometric total influence ( c = 1) was a valid criterion for the different stages of the structures in a wide range of situations provided that hydraulic roughness conditions remained high within the gully, e.g. through revegetation. Our total influence criterion involved shorter spacing than that habitually recommended in technical manuals for restoration, but was in line with those values found in spontaneous and stable step-pools systems, which might serve as a reference for man-made interventions.

Description: Senstitivity of water balance components to environmental changes in a mountainous watershed: uncertainty assessment based on models comparison Hydrology and Earth System Sciences Discussions, 10, 11983-12026, 2013 Author(s): E. Morán-Tejeda, J. Zabalza, K. Rahman, A. Gago-Silva, J. I. López-Moreno, S. Vicente-Serrano, A. Lehmann, C. L. Tague, and M. Beniston This paper evaluates the response of stream flow and other components of the water balance to changes in climate and land-use in a Pyrenean watershed. It further provides a measure of uncertainty in water resources forecasts by comparing the performance of two hydrological models: Soil and Water Assessment Tool (SWAT) and Regional Hydro-Ecological Simulation System (RHESSys). Regional Climate Model outputs for the 2021–2050 time-frame, and hypothetical (but plausible) land-use scenarios considering re-vegetation and wildfire processes were used as inputs to the models. Results indicate an overall decrease in river flows when the scenarios are considered, except for the post-fire vegetation scenario, in which stream flows are simulated to increase. However the magnitude of these projections varies between the two models used, as SWAT tends to produce larger hydrological changes under climate change scenarios, and RHESSys shows more sensitivity to changes in land-cover. The final prediction will therefore depend largely on the combination of the land-use and climate scenarios, and on the model utilized.

Description: Hydrologic impact of climate change on Murray Hotham catchment of Western Australia: a projection of rainfall-runoff for future water resources planning Hydrology and Earth System Sciences Discussions, 10, 12027-12076, 2013 Author(s): S. A. Islam, M. A. Bari, and A. H. M. F. Anwar Reduction of rainfall and runoff in recent years across South West Western Australia (SWWA) has drawn attention about climate change impact on water resources and its availability in this region. In this paper, hydrologic impact of climate change on Murray Hotham catchment in SWWA is investigated using multi-model ensemble approach. The Land Use Change Incorporated Catchment (LUCICAT) model was used for hydrologic modelling. Model calibration was performed using (5 km) grid rainfall data from Australian Water Availability Project (AWAP). Downscaled and bias corrected rainfall data from 11 General Circulation Models (GCMs) for Intergovernmental Panel on Climate Change (IPCC) emission scenarios A2 and B1 was used in LUCICAT model to derive rainfall and runoff scenarios for 2046–2065 (mid this century) and 2081–2100 (late this century). The results of climate scenarios were compared with observed past (1961–1980) climate. The mean annual rainfall averaged over the catchment during recent time (1981–2000) was reduced by 2.3% with respect to observed past (1961–1980) and resulting runoff reduction was found 14%. Compared to the past, the mean annual rainfall reductions, averaged over 11 ensembles and over the period for the catchment for A2 scenario are 13.6 and 23.6% for mid and late this century respectively while the corresponding runoff reductions are 36 and 74%. For B1 scenario, the rainfall reductions were 11.9 and 11.6% for mid and late this century and corresponding runoff reductions were 31 and 38%. Spatial distribution of rainfall and runoff changes showed that the rate of changes were higher in high rainfall part compared to the low rainfall part. Temporal distribution of rainfall and runoff indicate that high rainfall in the catchment reduced significantly and further reductions are projected resulting significant runoff reductions. A catchment scenario map has been developed through plotting decadal runoff reduction against corresponding rainfall reduction at four gauging stations for observed and projected period. This could be useful for planning future water resources in the catchment. Projection of rainfall and runoff made based on the GCMs varied significantly for the time periods and emission scenarios. Hence, considerable uncertainty involved in this study though ensemble mean was used to explain the findings.

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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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.

Description: Wind waves in the Black Sea: results of a hindcast study Natural Hazards and Earth System Science, 14, 2883-2897, 2014 Author(s): V. S. Arkhipkin, F. N. Gippius, K. P. Koltermann, and G. V. Surkova In this study we describe the wind wave fields in the Black Sea. The general aims of the work were the estimation of statistical wave parameters and the assessment of interannual and seasonal wave parameter variability. The domain of this study was the entire Black Sea. Wave parameters were calculated by means of the SWAN wave model on a 5 × 5 km rectangular grid. Initial conditions (wind speed and direction) for the period between 1949 and 2010 were derived from the NCEP/NCAR reanalysis. According to our calculations the average significant wave height on the Black Sea does not exceed 0.7 m. Areas of most significant heavy sea are the southwestern and the northeastern parts of the sea as expressed in the spatial distribution of significant wave heights, wave lengths and periods. Besides, long-term annual variations of wave parameters were estimated. Thus, linear trends of the annual total duration of storms and of their quantity are nearly stable over the hindcast period. However, an intensification of storm activity is observed in the 1960s–1970s.

Description: Complex networks, streamflow, and hydrometric monitoring system design Hydrology and Earth System Sciences Discussions, 11, 13663-13710, 2014 Author(s): M. Halverson and S. Fleming Network theory is applied to an array of streamflow gauges located in the Coast Mountains of British Columbia and Yukon, Canada. The goal of the analysis is to assess whether insights from this branch of mathematical graph theory can be meaningfully applied to hydrometric data, and more specifically, whether it may help guide decisions concerning stream gauge placement so that the full complexity of the regional hydrology is efficiently captured. The streamflow data, when represented as a complex network, has a global clustering coefficient and average shortest path length consistent with small-world networks, which are a class of stable and efficient networks common in nature, but the results did not clearly suggest a scale-free network. Stability helps ensure that the network is robust to the loss of nodes; in the context of a streamflow network, stability is interpreted as insensitivity to station removal at random. Community structure is also evident in the streamflow network. A community detection algorithm identified 10 separate communities, each of which appears to be defined by the combination of its median seasonal flow regime (pluvial, nival, hybrid, or glacial, which in this region in turn mainly reflects basin elevation) and geographic proximity to other communities (reflecting shared or different daily meteorological forcing). Betweenness analyses additionally suggest a handful of key stations which serve as bridges between communities and might therefore be highly valued. We propose that an idealized sampling network should sample high-betweenness stations, as well as small-membership communities which are by definition rare or undersampled relative to other communities, while retaining some degree of redundancy to maintain network robustness.

Description: Factors affecting flood insurance purchase in residential properties in Johor, Malaysia Natural Hazards and Earth System Science, 14, 3297-3310, 2014 Author(s): U. G. Aliagha, T. E. Jin, W. W. Choong, M. Nadzri Jaafar, and H. M. Ali High-impact floods have become a virtually annual experience in Malaysia, yet flood insurance has remained a grossly neglected part of comprehensive integrated flood risk management. Using discriminant analysis, this study seeks to identify the demand-side variables that best predict flood insurance purchase and risk aversion between two groups of residential homeowners in three districts of Johor State, Malaysia: those who purchased flood insurance and those who did not. Our results revealed an overall 34% purchase rate, with Kota Tinggi district having the highest (44%) and thus the highest degree of flood risk aversion. The Wilks' lambda F test for equality of group means, standardised discriminant function coefficients, structure correlation, and canonical correlation has clearly shown that there are strong significant attribute differences between the two groups of homeowners, based on the measures of objective flood risk exposure, subjective risk perception, and socio-economic cum demographic variables. However, the measures of subjective risk perception were found to be more predictive of flood insurance purchase and flood risk aversion.

Description: Monitoring and modelling of soil–plant interactions: the joint use of ERT, sap flow and Eddy Covariance data to characterize the volume of an orange tree root zone Hydrology and Earth System Sciences Discussions, 11, 13353-13384, 2014 Author(s): G. Cassiani, J. Boaga, D. Vanella, M. T. Perri, and S. Consoli Mass and energy exchanges between soil, plants and atmosphere control a number of key environmental processes involving hydrology, biota and climate. The understanding of these exchanges also play a critical role for practical purposes e.g. in precision agriculture. In this paper we present a methodology based on coupling innovative data collection and models in order to obtain quantitative estimates of the key parameters of such complex flow system. In particular we propose the use of hydro-geophysical monitoring via 4-D Electrical Resistivity Tomography (ERT) in conjunction with measurements of plant transpiration via sap flow and evapotranspiration from Eddy Covariance (EC). This abundance of data is fed to a spatially distributed soil model in order to characterize the distribution of active roots. We conducted experiments in an orange orchard in Eastern Sicily (Italy), characterized by the typical Mediterranean semi-arid climate. The subsoil dynamics, particularly influenced by irrigation and root uptake, were characterized mainly by the ERT setup, consisting of 48 buried electrodes on 4 instrumented micro boreholes (about 1.2 m deep) placed at the corners of a square (about 1.3 m in side) surrounding the orange tree, plus 24 mini-electrodes on the surface spaced 0.1 m on a square grid. During the monitoring, we collected repeated ERT and TDR soil moisture measurements, soil water samples, sap flow measurements from the orange tree and EC data. We conducted a laboratory calibration of the soil electrical properties as a function of moisture content and pore water electrical conductivity. Irrigation, precipitation, sap flow and ET data are available allowing knowledge of the system's long term forcing conditions on the system. This information was used to calibrate a 1-D Richards' equation model representing the dynamics of the volume monitored via 3-D ERT. Information on the soil hydraulic properties was collected from laboratory and field experiments. The successful results of the calibrated modeling exercise allow the quantification of the soil volume interested by root water uptake. This volume is much smaller (with a surface area less than 2 m 2 , and about 40 cm thickness) than expected and assumed in the design of classical drip irrigation schemes that prove to be losing at least half of the irrigated water that is not uptaken by the plants.

Description: Kinematics of a mass movement constrained by sparse and inhomogeneous data Natural Hazards and Earth System Science, 11, 1609-1618, 2011 Author(s): M. Karbon, E. Brückl, E. Hegedüs, and A. Preh On 12 February 2008, a landslide occurred along a 50 m high bank of the Danube river near Dunaszekcsö, Hungary. The initial state is only incompletely documented and the geodetic data acquired after the mass movement are sparse. A generalized 3-D topographic model of the landslide and its surrounding area was assembled and a representative longitudinal profile extracted. The reconstruction of the original surface is based on an orthophoto as well as on morphological considerations. Recorded observations include the locations of the outcrops of basal sliding surfaces, displacements at the main scarp and in the lower part of the slide, and a value to describe the total mass transport. Such sparse and inhomogeneous data were insufficient to derive a comprehensive documentation of the landslide or obtain adequate constraints for an accurate numerical analysis. Therefore, slider block models were fitted to the field data, which have only a small number of free parameters. A general view on the morphology of the mass movement justifies its classification as a rotational slide. A double slider block model fits all observational parameters within their error margin and supplies valuable information on the geometry of the slide. Estimates of the residual friction angles were derived and the question of reactivation was addressed. Finite Difference (FD) modelling and the application of conventional stability analysis support the geometry of the slider blocks and the computed average residual friction angles. Generally, the results are assumed to represent preliminary information, which could only be attained by the combination of the thinly distributed geodetic data with qualitative morphological observations and the implementation of a model. This type of information can be gained quickly and may be valuable for preliminary hazard mitigation measures or the planning of a comprehensive exploration and monitoring program.

Description: Looking for evidence of climate change impacts in the eastern Irish Sea Natural Hazards and Earth System Science, 11, 1641-1656, 2011 Author(s): L. S. Esteves, J. J. Williams, and J. M. Brown Although storminess is often cited as a driver of long-term coastal erosion, a lack of suitable datasets has only allowed objective assessment of this claim in a handful of case studies. This reduces our ability to understand and predict how the coastline may respond to an increase in "storminess" as suggested by global and regional climate models. With focus on 16 km of the Sefton coastline bordering the eastern Irish Sea (UK), this paper analyses available measured datasets of water level, surge level, wave height, wind speed and barometric pressure with the objective of finding trends in metocean climate that are consistent with predictions. The paper then examines rates of change in shoreline position over the period 1894 to 2005 with the aim of establishing relationships with climatic variability using a range of measured and modelled metocean parameters (with time spans varying from two to eight decades). With the exception of the mean monthly wind speed, available metocean data do not indicate any statistically significant changes outside seasonal and decadal cycles. No clear relationship was found between changes in metocean conditions and rates of shoreline change along the Sefton coast. High interannual variability and the lack of long-term measurements make unambiguous correlations between climate change and shoreline evolution problematic. However, comparison between the North Atlantic Oscillation winter index (NAOw) and coastline changes suggest increased erosion at times of decreasing NAOw values and reduced erosion at times of increasing NAOw values. Erosion tends to be more pronounced when decreasing NAOw values lead to a strong negative NAO phase. At present, anthropogenic changes in the local sediment budget and the short-term impact of extreme events are still the largest threat likely to affect coastal flooding and erosion risk in the short- and medium-term. Nevertheless, the potential impacts of climate change in the long-term should not be ignored.

Description: Comparing historical-hydrogeomorphological reconstitution and hydrological-hydraulic modelling in the estimation of flood-prone areas – a case study in Central Portugal Natural Hazards and Earth System Science, 11, 1669-1681, 2011 Author(s): P. P. Santos, A. O. Tavares, and A. I. A. S. S. Andrade The Arunca River basin in Central Portugal has a historical record of hazardous events related to floods, causing widespread disturbance. This article describes the application of two approaches based on well-known methods for the estimation of flood-prone areas: (i) historical-hydrogeomorphological reconstitution, applied to the entire Arunca River basin, and (ii) hydrological-hydraulic modelling, applied to four sections selected from different (upper, middle and lower) sectors of the basin and including urban and rural areas along the Arunca River. The mapping of the flood-prone areas obtained by these two methods was compared in order to identify the main differences and similarities. Human interventions (river channel and floodplain morphological changes) were found to be the main factor explaining the differences and similarities between the results obtained by both methods. The application of hydrological-hydraulic modelling proved important in reinforcing the results of the historical-hydrogeomorphological method; it also helped in complementing the results produced by the latter method in urban areas and in areas with insufficient historical records. The application of the historical-hydrogeomorphological method, in turn, allowed for the size of the flood-prone areas to be determined where the primary data (e.g. geometry, roughness and flow) was not accurate enough for hydrological-hydraulic modelling. The methodological approach adopted demonstrates the strong complementary relationship between the different existing methods for estimating flood-prone areas, and may be reproduced for other drainage basins.

Description: Domestic wells have high probability of pumping septic tank leachate Hydrology and Earth System Sciences Discussions, 8, 5701-5732, 2011 Author(s): J. E. Horn and T. Harter Onsite wastewater treatment systems such as septic systems are common in rural and semi-rural areas around the world; in the US, about 25–30 % of households are served by a septic system and a private drinking water well. Site-specific conditions and local groundwater flow are often ignored when installing septic systems and wells. Particularly in areas with small lots, thus a high septic system density, these typically shallow wells are prone to contamination by septic system leachate. Typically, mass balance approaches are used to determine a maximum septic system density that would prevent contamination of the aquifer. In this study, we estimate the probability of a well pumping partially septic system leachate. A detailed groundwater and transport model is used to calculate the capture zone of a typical drinking water well. A spatial probability analysis is performed to assess the probability that a capture zone overlaps with a septic system drainfield depending on aquifer properties, lot and drainfield size. We show that a high septic system density poses a high probability of pumping septic system leachate. The hydraulic conductivity of the aquifer has a strong influence on the intersection probability. We conclude that mass balances calculations applied on a regional scale underestimate the contamination risk of individual drinking water wells by septic systems. This is particularly relevant for contaminants released at high concentrations, for substances which experience limited attenuation, and those being harmful even in low concentrations.

Description: The contribution of engineering surveys by means of GPS to the determination of crustal movements in Istanbul Natural Hazards and Earth System Science, 11, 1705-1713, 2011 Author(s): M. Özyaşar and M. T. Özlüdemir Global Navigation Satellite Systems (GNSS) are space based positioning techniques and widely used in geodetic applications. Geodetic networking accomplished by engineering surveys constitutes one of these tasks. Geodetic networks are used as the base of all kinds of geodetic implementations, Co from the cadastral plans to the relevant surveying processes during the realization of engineering applications. Geodetic networks consist of control points positioned in a defined reference frame. In fact, such positional information could be useful for other studies as well. One of such fields is geodynamic studies that use the changes of positions of control stations within a network in a certain time period to understand the characteristics of tectonic movements. In Turkey, which is located in tectonically active zones and struck by major earthquakes quite frequently, the positional information obtained in engineering surveys could be very useful for earthquake related studies. For this purpose, a GPS (Global Positioning System) network of 650 stations distributed over Istanbul (Istanbul GPS Triangulation Network; abbreviated IGNA) covering the northern part of the North Anatolian Fault Zone (NAFZ) was established in 1997 and measured in 1999. From 1998 to 2004, the IGNA network was extended to 1888 stations covering an area of about 6000 km 2 , the whole administration area of Istanbul. All 1888 stations within the IGNA network were remeasured in 2005. In these two campaigns there existed 452 common points, and between these two campaigns two major earthquakes took place, on 17 August and 12 November 1999 with a Richter scale magnitude of 7.4 and 7.2, respectively. Several studies conducted for estimating the horizontal and vertical displacements as a result of these earthquakes on NAFZ are discussed in this paper. In geodynamic projects carried out before the earthquakes in 1999, an annual average velocity of 2–2.5 cm for the stations along the NAFZ were estimated. Studies carried out using GPS observations in the same area after these earthquakes indicated that point displacements vary depending on their distance to the epicentres of the earthquakes. But the directions of point displacements are similar. The results obtained through the analysis of the IGNA network also show that there is a common trend in the directions of point displacements in the study area. In this paper, the past studies about the tectonics of Marmara region are summarised and the results of the displacement analysis on the IGNA network are discussed.

Description: Use of ENVISAT ASAR Global Monitoring Mode to complement optical data in the mapping of rapid broad-scale flooding in Pakistan Hydrology and Earth System Sciences Discussions, 8, 5769-5809, 2011 Author(s): D. O'Grady, M. Leblanc, and D. Gillieson Envisat ASAR Global Monitoring Mode (GM) data are used to produce maps of the extent of the flooding in Pakistan which are made available to the rapid response effort within 24 h of acquisition. The high temporal frequency and independence of the data from cloud-free skies makes GM data a viable tool for mapping flood waters during those periods where optical satellite data is unavailable, which may be crucial to rapid response disaster planning, where thousands of lives are affected. Image differencing techniques are used, with pre-flood baseline image backscatter values being deducted from target values to eliminate regions with a permanent flood-like radar response due to volume scattering and attenuation, and to highlight the low response caused by specular reflection by open flood water. The effect of local incidence angle on the received signal is mitigated by ensuring that the deducted image is acquired from the same orbit track as the target image. Poor separability of the water class with land in areas beyond the river channels is tackled using a region-growing algorithm which seeks threshold-conformance from seed pixels at the center of the river channels. The resultant mapped extents are tested against MODIS SWIR data where available, with encouraging results.

Description: Soil weathering rates in 21 catchments of the Canadian Shield Hydrology and Earth System Sciences Discussions, 8, 5743-5768, 2011 Author(s): D. Houle, P. Lamoureux, N. Bélanger, M. Bouchard, C. Gagnon, S. Couture, and A. Bouffard Soil mineral weathering represents an essential source of nutrient base cation (Ca, Mg and K) for forest growth in addition to provide a buffering power against precipitation acidity for soils and surface waters. Weathering rates of base cations were obtained for 21 catchments located within the temperate and the boreal forest of the Canadian Shield with the geochemical model PROFILE. Weathering rates ranged from 0.58 to 4.46 kmol c ha −1 yr −1 and their spatial variation within the studied area was mostly in agreement with spatial variations in soil mineralogy. Weathering rates of Ca and Mg were significantly correlated ( r = 0.80 and 0.64) with their respective lake concentrations. Weathering rates of K and Na did not correlate with lake concentrations of K and Na. The modeled weathering rates for each catchment were also compared with estimations of net catchment exportations. The result show that modeled weathering rates of Ca were not significantly different than the net catchment exportations while modeled weathering rates of Mg were higher by 51 %. Larger differences were observed for K and Na weathering rates that were significantly different than net catchment exportations being 6.9 and 2.2 times higher than net exportations, respectively. The results for K were expected given its high reactivity with biotic compartments and suggest that most of the K produced by weathering reactions was retained within soil catchments and/or above ground biomass. This explanation does not apply to Na, however, which is a conservative element in forest ecosystems because of the insignificant needs of Na for soil microorganisms and above ground vegetations. It raises concern about the liability of the PROFILE model to provide reliable values of Na weathering rates. Overall, we concluded that the PROFILE model is powerful enough to reproduce spatial geographical gradients in weathering rates for relatively large areas as well as adequately predict absolute weathering rates values for the sum of base cations, Ca and Mg.

Description: Quantitative estimation of landslide risk from rapid debris slides on natural slopes in the Nilgiri hills, India Natural Hazards and Earth System Science, 11, 1723-1743, 2011 Author(s): P. Jaiswal, C. J. van Westen, and V. Jetten A quantitative procedure for estimating landslide risk to life and property is presented and applied in a mountainous area in the Nilgiri hills of southern India. Risk is estimated for elements at risk located in both initiation zones and run-out paths of potential landslides. Loss of life is expressed as individual risk and as societal risk using F-N curves, whereas the direct loss of properties is expressed in monetary terms. An inventory of 1084 landslides was prepared from historical records available for the period between 1987 and 2009. A substantially complete inventory was obtained for landslides on cut slopes (1042 landslides), while for natural slopes information on only 42 landslides was available. Most landslides were shallow translational debris slides and debris flowslides triggered by rainfall. On natural slopes most landslides occurred as first-time failures. For landslide hazard assessment the following information was derived: (1) landslides on natural slopes grouped into three landslide magnitude classes, based on landslide volumes, (2) the number of future landslides on natural slopes, obtained by establishing a relationship between the number of landslides on natural slopes and cut slopes for different return periods using a Gumbel distribution model, (3) landslide susceptible zones, obtained using a logistic regression model, and (4) distribution of landslides in the susceptible zones, obtained from the model fitting performance (success rate curve). The run-out distance of landslides was assessed empirically using landslide volumes, and the vulnerability of elements at risk was subjectively assessed based on limited historic incidents. Direct specific risk was estimated individually for tea/coffee and horticulture plantations, transport infrastructures, buildings, and people both in initiation and run-out areas. Risks were calculated by considering the minimum, average, and maximum landslide volumes in each magnitude class and the corresponding minimum, average, and maximum run-out distances and vulnerability values, thus obtaining a range of risk values per return period. The results indicate that the total annual minimum, average, and maximum losses are about US$ 44 000, US$ 136 000 and US$ 268 000, respectively. The maximum risk to population varies from 2.1 × 10 −1 for one or more lives lost to 6.0 × 10 −2 yr −1 for 100 or more lives lost. The obtained results will provide a basis for planning risk reduction strategies in the Nilgiri area.

Description: The role of catchment classification in rainfall-runoff modeling Hydrology and Earth System Sciences Discussions, 8, 6113-6153, 2011 Author(s): Y. He, A. Bárdossy, and E. Zehe A sound catchment classification scheme is a fundamental step towards improved catchment hydrology science and prediction in ungauged basins. Two categories of catchment classification methods are presented in the paper. The first one is based directly on physiographic properties and climatic conditions over a catchment and regarded as a Linnaean type or natural classification scheme. The second one is based on numerical clustering and regionalization methods and considered as a statistical or arbitrary classification scheme. This paper reviews each category including what has been done since recognition of the intrinsic value of catchment classification, what is being done in the current research, as well as what is to be done in the future.

Description: Infiltration-soil moisture redistribution under natural conditions: experimental evidence as a guideline for realizing simulation models Hydrology and Earth System Sciences Discussions, 8, 6199-6225, 2011 Author(s): R. Morbidelli, C. Corradini, C. Saltalippi, A. Flammini, and E. Rossi The evolution in time, t , of the experimental soil moisture vertical profile under natural conditions is investigated in order to address the corresponding simulation modelling. The measurements were conducted in a plot with a bare silty loam soil. The soil water content, θ, was continuously monitored at different depths, z , using a Time Domain Reflectometry (TDR) system. For each profile four buriable three-rod waveguides were inserted horizontally at different depths (5, 15, 25 and 35 cm). In addition, we used sensors of air temperature and relative humidity, wind speed, solar radiation, evaporation and rain as supports for the application of selected simulation models, as well as for the detection of elements leading to their improvement. The results indicate that, under natural conditions, very different trends of the θ( z , t ) function can be observed in the given fine-textured soil, where the formation of a sealing layer over the parent soil requires an adjustment of the simulation modelling commonly used for hydrological applications. In particular, because of the considerable variations in the shape of the moisture content vertical profile as a function of time, a generalization of the existing models should incorporate a representation of the variability in time of the saturated hydraulic conductivity of the uppermost soil. This conclusion is supported by the fact that the observed shape of θ( z ) can be appropriately reproduced by adopting this approach, however the observed rainfall rate and the occurrence of freeze-thaw cycles with high soil moisture contents have to be explicitly incorporated.

Description: Trends in climatic variables and future reference evapotranspiration in Duero Valley (Spain) Natural Hazards and Earth System Science, 11, 1795-1805, 2011 Author(s): R. Moratiel, R. L. Snyder, J. M. Durán, and A. M. Tarquis The impact of climate change and its relation with evapotranspiration was evaluated in the Duero River Basin (Spain). The study shows possible future situations 50 yr from now from the reference evapotranspiration (ET o ). The maximum temperature ( T max ), minimum temperature ( T min ), dew point ( T d ), wind speed ( U ) and net radiation ( R n ) trends during the 1980–2009 period were obtained and extrapolated with the FAO-56 Penman-Montheith equation to estimate ET o . Changes in stomatal resistance in response to increases in CO 2 were also considered. Four scenarios were done, taking the concentration of CO 2 and the period analyzed (annual or monthly) into consideration. The scenarios studied showed the changes in ET o as a consequence of the annual and monthly trends in the variables T max , T min , T d , U and R n with current and future CO 2 concentrations (372 ppm and 550 ppm). The future ET o showed increases between 118 mm (11 %) and 55 mm (5 %) with respect to the current situation of the river basin at 1042 mm. The months most affected by climate change are May, June, July, August and September, which also coincide with the maximum water needs of the basin's crops.

Description: Possible interrelation between the lead time of precursory seismic electric signals (SES) and geodynamics in Aegean Sea Natural Hazards and Earth System Science, 11, 1599-1603, 2011 Author(s): E. Dologlou The seismicity of the last 15 years in the Aegean Sea revealed that earthquakes ( M w 〉 5) with epicentres falling within the Sporades basin and the confined area north of Samos island were preceded by electric seismic signals (SES) with a remarkably long lead time. A possible explanation of this behaviour by means of specific tectonics and geodynamics which characterise these two regions, such as a significant small crustal thickness and a high heat flow rate, has been attempted. New data seem to strengthen the above hypothesis.

Description: Sand box experiments to evaluate the influence of subsurface temperature probe design on temperature based water flux calculation Hydrology and Earth System Sciences Discussions, 8, 6155-6197, 2011 Author(s): M. Munz, S. E. Oswald, and C. Schmidt Quantification of subsurface water fluxes based on the one dimensional solution to the heat transport equation depends on the accuracy of measured subsurface temperatures. The influence of temperature probe setup on the accuracy of vertical water flux calculation was systematically evaluated in this experimental study. Four temperature probe setups were installed into a sand box experiment to measure temporal highly resolved vertical temperature profiles under controlled water fluxes in the range of ±1.3 m d −1 . Pass band filtered time series provided amplitude and phase of the diurnal temperature signal varying with depth depending on water flux. Amplitude ratios of setups directly installed into the saturated sediment significantly varied with sand box hydraulic gradients. Amplitude ratios provided an accurate basis for the analytical calculation of water flow velocities, which matched measured flow velocities. Calculated flow velocities were sensitive to thermal properties of saturated sediment and to probe distance, but insensitive to thermal dispersivity equal to solute dispersivity. Amplitude ratios of temperature probe setups indirectly installed into piezometer pipes were influenced by thermal exchange processes within the pipes and significantly varied with water flux direction only. Temperature time lags of small probe distances of all setups were found to be insensitive to vertical water flux.

Description: Daytime variations of foE connected to earthquakes Natural Hazards and Earth System Science, 11, 1807-1812, 2011 Author(s): E. V. Liperovskaya, V. V. Bogdanov, P.-F. Biagi, C.-V. Meister, and V. A. Liperovsky In the present work it is shown that, in accordance with the observations of the vertical sounding station "Tashkent", the critical foE-frequency of the daytime E-layer increases about one day before winter-earthquakes with magnitudes M 〉 5 and depths of the epicentre of h 〈 60 km, which appeared at distances of R 〈 2000 km from the station. The reliability of the result is larger than 99 %. The phenomenon is not observed for summer-earthquakes. It seems to be determined by the atmospheric wind system. Further, the variations of the foE-frequency are compared with possible simultaneous variations of the critical frequency foF2 of the F2-layer. First results show that only very large changes of the ionisation density in the E-layer influence the ionisation density in the F-region. Therefore, no synchronous growth of the foE- and foF2-frequencies 1–2 days before seismic shocks could be observed.

Description: The use of LIDAR as a data source for digital elevation models – a study of the relationship between the accuracy of digital elevation models and topographical attributes in northern peatlands Hydrology and Earth System Sciences Discussions, 8, 5497-5522, 2011 Author(s): A. Hasan, P. Pilesjö, and A. Persson It is important to study the factors affecting estimates of wetness since wetness is crucial in climate change studies. The availability of digital elevation models (DEMs) generated with high resolution data is increasing, and their use is expanding. LIDAR earth elevation data have been used to create several DEMs with different resolutions, using various interpolation parameters, in order to compare the models with collected surface data. The aim is to study the accuracy of DEMs in relation to topographical attributes such as slope and drainage area, which are normally used to estimate the wetness in terms of topographic wetness indices. Evaluation points were chosen from the high-resolution LIDAR dataset at a maximum distance of 10 mm from the cell center for each DEM resolution studied, 0.5, 1, 5, 10, 30 and 90 m. The interpolation method used was inverse distance weighting method with four search radii: 1, 2, 5 and 10 m. The DEM was evaluated using a quantile-quantile test and the normalized median absolute deviation. The accuracy of the estimated elevation for different slopes was tested using the DEM with 0.5 m resolution. Drainage areas were investigated at three resolutions, with coinciding evaluation points. The ability of the model to generate the drainage area at each resolution was obtained by pairwise comparison of three data subsets. The results show that the accuracy of the elevations obtained with the DEM model are the same for different resolutions, but vary with search radius. The accuracy of the values (NMAD of errors) varies from 29.7 mm to 88.9 mm, being higher for flatter areas. It was also found that the accuracy of the drainage area is highly dependent on DEM resolution. Coarse resolution yielded larger estimates of the drainage area but lower slope values. This may lead to overestimation of wetness values when using a coarse resolution DEM.

Description: Multivariate design via Copulas Hydrology and Earth System Sciences Discussions, 8, 5523-5558, 2011 Author(s): G. Salvadori, C. De Michele, and F. Durante Calculating return periods and design quantiles in a multivariate framework is a difficult problem: essentially, this is due to the lack of a natural total order in multi-dimensional Euclidean spaces. This paper tries to make the issue clear. First, we outline a possible way to introduce a coherent notion of multivariate total order, and discuss its consequences on the calculation of multivariate return period: in particular, the latter is based on Copulas and the Kendall's measure, which provides a consistent notion of multivariate quantile. Secondly, we introduce several approaches for the identification of critical design events: these latter quantities are of utmost importance in practical applications, but their calculation is yet limited, due to the lack of a suitable theoretical setting where to embed the problem. Throughout the paper, a case study involving the behavior of a dam is used to illustrate the new concepts outlined in this work.

Description: Influence of initial heterogeneities and recharge limitations on the evolution of aperture distributions in carbonate aquifers Hydrology and Earth System Sciences Discussions, 8, 5631-5666, 2011 Author(s): B. Hubinger and S. Birk Karst aquifers evolve where the dissolution of soluble rocks causes the enlargement of discrete pathways along fractures or bedding planes, thus creating highly conductive solution conduits. To identify general interrelations between hydrogeological conditions and the properties of the evolving conduit systems the aperture-size frequency distributions resulting from generic models of conduit evolution are analysed. For this purpose, a process-based numerical model coupling flow and rock dissolution is employed. Initial protoconduits are represented by tubes with log-normally distributed aperture sizes with a mean of 0.5 mm. Apertures are spatially uncorrelated and widen up to the metre range due to dissolution by chemically aggressive waters. Several examples of conduit development are examined focussing on influences of the initial heterogeneity and the available amount of recharge. If the available recharge is sufficiently high the evolving conduits compete for flow and those with large apertures and high hydraulic gradients attract more and more water. As a consequence, the positive feedback between increasing flow and dissolution causes the breakthrough of a conduit pathway connecting the recharge and discharge sides of the modelling domain. Under these competitive flow conditions dynamically stable bimodal aperture distributions are found to evolve, i.e. a certain percentage of tubes continues to be enlarged while the remaining tubes stay small-sized. The percentage of strongly widened tubes is found to be independent of the breakthrough time and decreases with increasing heterogeneity of the initial apertures and decreasing amount of available water. If the competition for flow is suppressed because the availability of water is strongly limited breakthrough of a conduit pathway is inhibited and the conduit pathways widen very slowly. The resulting aperture distributions are found to be unimodal covering some orders of magnitudes in size. Under these suppressed flow conditions the entire range of apertures continues to be enlarged. Hence, the number of tubes reaching aperture sizes in the order of centimetres or decimetres continues to increase with time and in the long term may exceed the number of large-sized tubes evolving under competitive flow conditions. This suggests that conduit development under suppressed flow conditions may significantly enhance the permeability of the formation e.g. in deep-seated carbonate settings.

Description: Skewness as measure of the invariance of instantaneous renormalized drop diameter distributions – Part 1: Convective vs. stratiform precipitation Hydrology and Earth System Sciences Discussions, 8, 5605-5629, 2011 Author(s): M. Ignaccolo and C. De Michele We investigate the variability of the instantaneous distribution shape of the renormalized drop diameter making use of the third order central moment: the skewness . Disdrometer data, collected at Darwin Australia, are considered either as whole or as divided in convective and stratiform precipitation intervals. We show that in all cases the distribution of the skewness is strongly peaked around 0.64. This allows to identify a most common distribution of renormalized drop diameters and two main variations, one with larger and one with smaller skewness. The distributions' shapes are independent from the stratiform vs. convective classification.

Description: The sensitivity of land emissivity estimates from AMSR-E at C and X bands to surface properties Hydrology and Earth System Sciences Discussions, 8, 5667-5699, 2011 Author(s): H. Norouzi, M. Temimi, W. B. Rossow, C. Pearl, M. Azarderakhsh, and R. Khanbilvardi Microwave observations at low frequencies exhibit more sensitivity to surface and subsurface properties with little interference from the atmosphere. The objective of this study is to develop a global land emissivity product using passive microwave observations from the Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) and to investigate its sensitivity to land surface properties. The developed product complements existing land emissivity products from SSM/I and AMSU by adding land emissivity estimates at two lower frequencies, 6.9 and 10.65 GHz (C- and X-band, respectively). Observations at these low frequencies penetrate deeper into the soil layer. Ancillary data used in the analysis, such as surface skin temperature and cloud mask, are obtained from International Satellite Cloud Climatology Project (ISCCP). Atmospheric properties are obtained from the TIROS Operational Vertical Sounder (TOVS) observations to determine the small upwelling and downwelling atmospheric emissions as well as the atmospheric transmission. A sensitivity test confirms the small effect of the atmosphere but shows that skin temperature accuracy can significantly affect emissivity estimates. Retrieved emissivities at C- and X-bands and their polarization differences exhibit similar patterns of variation with changes in land cover type, soil moisture, and vegetation density as seen at SSM/I-like frequencies (Ka and Ku bands). The emissivity maps from AMSR-E at these higher frequencies agree reasonably well with the existing SSM/I-based product. The inherent but small discrepancy introduced by the difference between SSM/I and AMSR-E frequencies and incidence angles has been examined and found to be small. Large differences between emissivity estimates from ascending and descending overpasses were found at the lower frequencies due to the inconsistency between the thermal IR skin temperatures and passive microwave brightness temperatures which can come from below the surface. This issue must be addressed in future studies to improve the accuracy of the emissivity estimates at lower frequencies.

Description: Influence of soil parameters on the skewness coefficient of the annual maximum flood peaks Hydrology and Earth System Sciences Discussions, 8, 5559-5604, 2011 Author(s): A. Gioia, V. Iacobellis, S. Manfreda, and M. Fiorentino Understanding the spatial variability of key parameters of flood probability distributions represents a strategy to provide insights on hydrologic similarity and building probabilistic models able to reduce the uncertainty in flood prediction in ungauged basins. In this work, we exploited the theoretically derived distribution of floods TCIF (Gioia et al., 2008), based on two different threshold mechanisms associated respectively to ordinary and extraordinary events. The model is based on the hypotheses that ordinary floods are generally due to rainfall events exceeding a threshold infiltration rate in a small source area, while the so-called outlier events, responsible of the high skewness of flood distributions, are triggered when severe rainfalls exceed a storage threshold over a large portion of the basin. Within this scheme, a sensitivity analysis was performed in order to analyze the effects of climatic and geomorphologic parameters on the skewness coefficient. In particular, the analysis was conducted investigating the influence on flood distribution of physical factors such as rainfall intensity, soil infiltration capacity, and basin area, in order to provide insights in catchment classification and process conceptualization.

Description: Shallow circulation groundwater – the main type of water containing hazardous radon concentration Natural Hazards and Earth System Science, 11, 1695-1703, 2011 Author(s): T. A. Przylibski The main factors affecting the value of 222 Rn activity concentration in groundwater are the emanation coefficient of reservoir rocks ( K em ), the content of parent 226 Ra in these rocks ( q ), changes in the volume and flow velocity as well as the mixing of various groundwater components in the circulation system. The highest values of 222 Rn activity concentration are recorded in groundwaters flowing towards an intake through strongly cracked reservoir rocks undergoing weathering processes. Because of these facts, waters with hazardous radon concentration levels, i.e. containing more than 100 Bq dm −3 222 Rn, could be characterised in the way that follows. They are classified as radon waters, high-radon waters and extreme-radon waters. They belong to shallow circulation systems (at less than a few dozen metres below ground level) and are contemporary infiltration waters, i.e. their underground flow time ranges from several fortnights to a few decades. Because of this, these are usually poorly mineralised waters (often below 0.2–0.5 g dm −3 ). Their resources are renewable, but also vulnerable to contamination. Waters of this type are usually drawn from private intakes, supplying water to one or at most a few households. Due to an increased risk of developing lung tumours, radon should be removed from such waters when still in the intake. To achieve this aim, appropriate legislation should be introduced in many countries.

Description: Precipitation change in Southern Italy linked to global scale oscillation indexes Natural Hazards and Earth System Science, 11, 1683-1694, 2011 Author(s): T. Caloiero, R. Coscarelli, E. Ferrari, and M. Mancini This study investigates precipitation variability in five regions of Southern Italy (Campania, Apulia, Basilicata, Calabria and Sicily) using a homogeneous database of about 70 rain gauges with more than 50 years of observation. First, a statistical analysis was performed through the Mann-Kendall non-parametric test in order to determine rainfall the trend on both yearly and seasonal scales. Then, the relationship between the rainfall and some teleconnection pattern indexes was investigated using Spearman's test. The results show remarkable statistically significant negative trends for annual and winter aggregations in most part of the series. Moreover, a strong correlation has emerged between the teleconnection patterns and precipitation in Southern Italy, particularly in winter and on the Tyrrhenian side of the study area.

Description: Technical note: Towards a continuous classification of climate using bivariate colour mapping Hydrology and Earth System Sciences Discussions, 8, 5733-5742, 2011 Author(s): A. J. Teuling Climate is often defined in terms of discrete classes. Here I use bivariate colour mapping to show that the global distribution of Köppen-Geiger climate classes can largely be reproduced by combining the simple means of two key states of the climate system (i.e., air temperature and relative humidity). This allows for a classification that is not only continuous in space, but can be applied at and transferred between timescales ranging from minutes to decades.

Description: Reconstruction of Atlantic historical winter coastal storms in the Spanish coasts of the Gulf of Cadiz, 1929–2005 Natural Hazards and Earth System Science, 11, 1715-1722, 2011 Author(s): P. Ribera, D. Gallego, C. Pena-Ortiz, L. Del Rio, T. A. Plomaritis, and J. Benavente This paper presents the reconstruction of a climatological series of winter coastal storms on the northern coasts of the Gulf of Cadiz. This series has been put together using information extracted from regional and local Spanish newspapers. It includes all the storms coming from the Atlantic sector that have been detected during the winter season, from October to March, between 1929 and 2005. In order to validate this historical storm series, it has been compared with storms series identified from quasi-observational data and using different wave heights as thresholds to decide what is to be considered as a coastal storm. Nearly 2.6 reports per year about coastal storms are published in the press which correspond to waves of 3.6 m high or more and to prevailing winds from a direction ranging between SSW and WNW. A long- term positive trend has been detected for the complete storm series. If only the instrumental period is analysed, no significant trend is detected. It is suggested that this difference might be associated with the impact of the North Atlantic Oscillation over the occurrence of storms in this area.

Description: The neutral temperature in the ionospheric dynamo region and the ionospheric F region density during Wenchuan and Pingtung Doublet earthquakes Natural Hazards and Earth System Science, 11, 1759-1768, 2011 Author(s): Y. Y. Sun, K.-I. Oyama, J. Y. Liu, H. K. Jhuang, and C. Z. Cheng One of the possible candidates which modifies the ionosphere before large earthquake is electric field. We presume that the electric field associated with large earthquakes is generated in the ionosphere dynamo region (100–120 km). This paper tries to identify the evidence of the contribution of the neutral atmosphere in the dynamo region. The relationship between the critical frequency at the F2 peak (foF2) and the height profile of the neutral atmosphere temperature was studied for two large earthquakes: Wenchuan, 2008 and Pingtung Doublet, 2006. It is found that the wave amplitude of the vertical wavelength 20–30 km which is usually superposed on the height profile of the neutral atmosphere temperature enhances when the foF2 increases. The correlation between the wave amplitude and foF2 is found better along a longitudinal direction than along latitude direction.