ALBERT

Description: Land use change as conversion pasture to forest produces several changes on hydrological cycle. In this paper we analyze the effects on stream discharge of afforestation of a small watershed devoted to pasture using the HBV hydrological model. Streamflow data obtained over the first ten years after planting were employed to evaluate the capacity of HBV model to simulate hydrological behavior of catchment after afforestation. Obtained results indicate that the estimation of streamflow was accurate as reflected by statistics (R 2  = 0.90, NSC = 0.89 and PBIAS = 0.34). Afterwards, streamflow under pasture land use (if afforestation had not occurred) was simulated using hydro-meteorological data collected during the period of study and model parameters optimized previously, together with two parameters, pcorr and cevpfo , that were adjusted for pasture conditions. HBV model results indicate that afforestation produced a water yield reduction around 2000 mm (22% of total stream discharge) during the first ten years of planting growth. The differences between forest and pasture land cover are increasing in all seasons year by year. The greatest streamflow reduction was observed in wet period (autumn and winter) with 76% of total reduction. In summer, streamflow reduction represents only 3% of total, however, represents 24.7% of discharge in this season. Streamflow reduction was related to increase of rainfall interception (mainly in wet periods) and the increase of evapotranspiration by plantation in dry periods. This article is protected by copyright. All rights reserved.

Description: Air flows from the atmosphere into an unconfined aquifer when the water table falls during pumping tests. Pumping test results in unconfined aquifers may be significantly affected by low-permeability zones (LPZs) near the initial water table position because they restrict the downward movement of air. A transient, three-dimensional (3D) air-water two-phase flow model is employed to investigate numerically the effects of local heterogeneity on pumping test results in unconfined aquifers. Two cases of local heterogeneities are considered herein: a LPZ around the pumping well and on one side of the pumping well. Results show that the drawdown with the LPZ is significantly greater than that of the homogeneous aquifer. The differences in drawdown are the most significant at intermediate times and gradually diminish at later times. The LPZ significantly reduces air flow from the atmosphere to the aquifer. The pore air velocity in the LPZ is very low. The air pressure at the observation point under the LPZ when air begins to enter is significantly lower than the air pressure of the homogeneous aquifer at the same point. After that, the air pressure increases quickly and then increases slowly. The time for the air pressure to reach the atmospheric pressure is significantly longer. This article is protected by copyright. All rights reserved.

Description: Over the last half century, since logging for timber production became intensive, Borneo has lost much of its pristine tropical forests. The long-term consequences of associated decline in precipitation are evident, and might, in turn, cause much more severe deforestation. This article is protected by copyright. All rights reserved.

Description: Accepting the concept of standardization introduced by the standardized precipitation index (SPI), similar methodologies have been developed to construct some other standardized drought indices such as the standardized precipitation evapotranspiration index (SPEI). In this study, the authors provided deep insight into the SPEI and recognized potential deficiencies/limitations in relating to the climatic water balance it used. By coupling another well-known Palmer drought severity index (PDSI), we proposed a new standardized Palmer drought index (SPDI) through a moisture departure probabilistic approach, which allows multiscalar calculation for accurate temporal and spatial comparison of the hydro-meteorological conditions of different locations. Using datasets of monthly precipitation, temperature and soil available water capacity, the moisture deficit/surplus was calculated at multiple temporal scales and a couple of techniques were adopted to adjust corresponding time series to a generalized extreme value distribution out of several candidates. Results of the historical records (1900–2012) for diverse climates by multiple indices showed that the SPDI was highly consistent and correlated with the SPEI and self-calibrated PDSI (SC-PDSI) at most analyzed time scales. Furthermore, a simple experiment of hypothetical temperature and/or precipitation change scenarios also verified the effectiveness of this newly-derived SPDI index in response to climate change impacts. Being more robust and preferable in spatial consistency and comparability as well as combining the simplicity of calculation with sufficient accounting of the physical nature of water supply and demand relating to droughts, the SPDI is promising to serve as a competent reference and alternative for drought assessment and monitoring. This article is protected by copyright. All rights reserved.

Description: Modeled hydrologic processes are represented in a set of numerical equations, the complexity of which can be measured by the total number of variables needed. A single dominant hydrologic process could control the hydrologic response of a watershed, and so the identification of the corresponding dominant variable(s) would aid in identifying a parsimonious model and in collecting more reliable data. By accounting for both model complexity and serial correlation in the variables, a model is used to identify the dominant variables for representing watershed scale streamflow, sediment transport, and phosphorus yields. Long-term water quantity and quality data was used to show that rainfall and non-linear soil water storage were the dominant variables for weekly streamflow, suspended sediment, and particulate phosphorus. Model accuracy did not consistently improve when other statistically significant variables were included. The results suggest that improved model performance may not justify the added model complexity. As such, identification of dominant variables would be the priority for developing parsimonious hydrologic models, especially at watershed scales. This article is protected by copyright. All rights reserved.

Description: We applied graphical methods and multivariate statistics to understand impacts of an unsewered slum catchment on nutrients and hydrochemistry of groundwater in Kampala, Uganda. Data was collected from 56 springs (groundwater), 22 surface water sites and 13 rain samples. Groundwater was acidic and dominated by Na, Cl and NO 3 . These ions were strongly correlated indicating pollution originating from wastewater infiltration from on-site sanitation systems. Results also showed that rain, which was acidic, impacted on groundwater chemistry. Using Q-mode hierarchical cluster analysis, we identified three distinct water quality groups. The first group had springs dominated by Ca-Cl-NO 3 , low values of EC, pH and cations, and relatively high NO 3 values. These springs were shown to have originated from the acidic rains because their chemistry closely corresponded to ion concentrations that would occur from rainfall recharge, which was around 3.3 times concentrated by evaporation. The second group had springs dominated by Na-K-Cl-NO 3 and Ca-Cl-NO 3 , low pH but with higher values of EC, NO 3 and cations. We interpreted these as groundwater affected by both acid rain and infiltration of wastewater from urban areas. The third group had the highest EC values (average of 688 μS/cm), low pH and very high concentrations of NO 3 (average of 2.15 mmol/L) and cations. These springs were exclusively located in slum areas and we interpreted these springs as groundwater affected by infiltration of wastewater from poorly sanitized slums areas. Surface water was slightly reducing and eutrophic due to wastewater effluents, but the contribution of groundwater to nutrients in surface water was minimal because o-PO 4 was absent whereas NO 3 was lost by denitification. Our findings suggest that groundwater chemistry in the catchment is strongly influenced by anthropogenic inputs derived from nitrogen-containing rains and domestic wastewater. This article is protected by copyright. All rights reserved.

Description: We examined the water balance a forested ombrotrophic peatland and adjacent burned peatland in the boreal plain of western Canada over a three-year period. Complete combustion of foliage and fine branches dramatically increased shortwave radiation inputs to the peat surface while halting all tree transpiration at the burned site. End-of-winter snowpack was 7-25% higher at the burned site likely due to decreased ablation from the tree canopy at the unburned site. Shrub regrowth at the burned site was rapid post-fire, and shading by the shrub canopy in the burned site approached that of the unburned site within three years after fire. Site-averaged surface resistance to evaporation was not different between sites, though surface resistance in hollows was lower in the burned site. Water loss at both burned and unburned sites is largely driven by surface evaporative losses. Evaporation at the burned site marginally exceeded the sum of pre-fire transpiration and interception at the unburned site, suggesting that ET during the growing season was 2 0–40  mm greater at the the burned peatland. While the net change in water storage during the growing season was largely unchanged by fire, the lack of low-density surface peat in the burned site appears to have decreased specific yield, leading to greater water table decline at the burned site despite similar net change in storage. This article is protected by copyright. All rights reserved.

Description: The eddy covariance (EC) method was used in a 30-month study to quantify evapotranspiration (ET) and vegetation coefficient (K CW ) for a wetland on a ranch in subtropical south Florida. To evaluate the errors in ET estimates, the EC-based ET (ET C.-EC ) and the Food and Agricultural Organization (FAO) Penman-Monteith (PM) based ET (ET C.-PM ) estimates (with literature crop coefficient K C ) were compared to each other. The ET C.-EC and FAO-PM reference ET were used to develop K CW . Regression models were developed to estimate K CW using climatic and hydrologic variables. Annual and daily ET C.-EC values were 1152 mm and 3.27 mm, respectively. The FAO-PM model underestimated ET by 25% with ET C.-EC being statistically higher than ET C.-PM . The K CW varied from 0.79 (December) to 1.06 (November). The mean K CW for dry (November-April) season (0.95) was much higher than values reported for wetlands in literature, while wet (May-October) season K CW (0.97) was closer to literature values. Higher than expected K CW values during dry season were due to higher temperature, lower humidity and perennial wetland vegetation. Regression analyses showed that factors affecting the K CW were different during the dry (soil moisture, temperature, and relative humidity) and wet (net radiation, inundation, and wind speed) seasons. Separate regression models for the dry and wet seasons were developed. ET and K CW from this study, one of the first for the agricultural wetlands in subtropical environment, will help improve the ET estimates for similar wetlands. This article is protected by copyright. All rights reserved.

Description: Research on runoff processes to date has focused on the differences between the main divisions of runoff partitioning. Indeed, our major advancements in runoff theory have come with new differentiations of different forms of overland flow and subsurface stormflow. These studies of ‘how runoff processes are different’ have resulted in our current summaries of runoff regimes conceptually (e.g. the Variable Source Area (VSA) concept) and codified in our models (e.g. TOPMODEL and its derivatives). This article is protected by copyright. All rights reserved.

Description: Supraglacial channels are an important mechanism for surface water transport over the ablation zone of western Greenland. The first assessment of the spatio-temporal distribution of surface melt channels and their relationship to supraglacial lakes over the Jakobshavn Isbrae region of Western Greenland was analyzed using Landsat ETM + panchromatic images during the 2007 melt season. A total of 1188 melt channels were delineated and show an increase in the number of melt channels throughout the season, reaching a peak on August 9. Water-filled melt channels advanced to maximum elevation of 1647 m on August 9, and attained a minimum average slope of 0.009 on July 8. The ablation zone demonstrates two hydrologic modes, where crevasse and moulin terminating channels dominate at elevations 〈800 m and higher order channel networks 〉800 m. Development of higher order networks is interrupted by flow divergence due to partitioning of melt water into vertical infiltration through moulins and crevasse fields prevalent at lower elevations. Tributary and Connector networks between 800 to 1200 m in elevation are correlated with fewer lake occurrences, relatively lower surface velocities (~50 ma -1 ) and ice flow dominated by internal deformation over basal sliding. High order channels are associated with lake basins that exceed melt water storage capacity. Evolution of channel networks is coupled to changes in melt water production, runoff, and ice dynamics with implication for the englacial and subglacial environments. This article is protected by copyright. All rights reserved.

Description: This paper reviews the use of the Generalised Likelihood Uncertainty Estimation (GLUE) methodology in the 20 years since the paper by Beven and Binley (1992) in Hydrological Processes, which is now one of the most highly cited papers in hydrology. The original conception, the on-going controversy it has generated, the nature of different sources of uncertainty and the meaning of the GLUE prediction uncertainty bounds, are discussed. The hydrological, rather than statistical, arguments about the nature of model and data errors and uncertainties that are the basis for GLUE are emphasised. The application of the IHDM model to the Gwy catchment at Plynlimon presented in the original paper is revisited, using a much larger sample of models, a wider range of likelihood evaluations and new visualisation techniques. It is concluded that there are good reasons to reject this model for that data set. This is a positive result in a research environment in that it requires improved models or data to be made available. In practice, there may be ethical issues of using outputs from models for which there is evidence for model rejection in decision making. Finally some suggestions for what is needed in the next 20 years are provided. This article is protected by copyright. All rights reserved.

Description: A short-term flood inundation prediction model has been formulated based on the combination of the super-tank model, forced with downscaled rainfall from a global numerical weather prediction model, and a one-dimensional hydraulic model. Different statistical methods for downscaled rainfall have been explored, taking into account the availability of historical data. It has been found that the full implementation of a statistical downscaling model considering physically based corrections to the numerical weather prediction model output for rainfall prediction performs better compared with an altitudinal correction method. The integration of the super-tank model into the one-dimensional hydraulic model demonstrates a minimal requirement for the calibration of rainfall-runoff and flood propagation models. Updating the model with antecedent rainfall and regular forecast renewal has enhanced the model's capabilities as a result of the data assimilation processes of the runoff and numerical weather prediction models. The results show that the predicted water levels demonstrate acceptable agreement with those measured by stream gauges and comparable to those reproduced using the actual rainfall. Moreover, the predicted flood inundation depth and extent exhibit reasonably similar tendencies to those observed in the field. However, large uncertainties are observed in the prediction results in lower, flat portions of the river basin where the hydraulic conditions are not properly analysed by the one-dimensional flood propagation model. This article is protected by copyright. All rights reserved.

Description: This study aims to analyse the combined impacts of future discharges and sea levels on erosion-sedimentation potential, and its seasonal changes, in a ~43 km long coastal river reach of South-West Finland. To our knowledge, this kind of combined study has not been performed before. In addition to surveying the present erosion-sedimentation conditions, the daily erosion-sedimentation potential is simulated with a one dimensional hydrodynamic model for the 1971–2000 and 2070–2099 periods by applying four discharge scenarios. Different sea level stages are also employed in the simulations. All scenarios forecast increasing autumn and winter discharges, but diminishing summer discharges. This indicates increasing river channel erosion, particularly during winters and autumns. Although discharge changes have altogether a greater influence on erosion-sedimentation potential, the importance of sea level changes on sedimentation is noticeable in the estuary. The rising sea level scenarios increase the sedimentation potential. In total, by 2070–2099 the erosion potential may increase in most parts of the study area. This article is protected by copyright. All rights reserved.

Description: Northern peatlands are a large source of atmospheric methane (CH 4 ) and both a source and sink of atmospheric carbon dioxide (CO 2 ). The rate and temporal variability in gas exchanges with peat soils is directly related to the spatial distribution of these free-phase gases within the peat column. In this paper we present results from surface and borehole ground penetrating radar (GPR) surveys – constrained with direct soil and gas sampling – that compare the spatial distribution of gas accumulations in two raised bogs: one in Wales (UK), the other in Maine (USA). Although the two peatlands have similar average thickness, physical properties of the peat matrix differ, particularly in terms of peat type and degree of humification. We hypothesize that these variations in physical properties are responsible for the differences in gas distribution between the two peatlands characterized by: 1) gas content up to 10.8 % associated with woody peat and presence of wood layers in Caribou Bog (Maine), and 2) a more homogenous distribution with gas content up to 5.7 % at the surface (i.e. 〈 0.5 m deep) in Cors Fochno (Wales). Our results highlight the variability in biogenic gas accumulation and distribution across peatlands and suggest that the nature of the peat matrix has a key role in defining how biogenic gas accumulates within, and is released to the atmosphere from, peat soils. This article is protected by copyright. All rights reserved.

Description: The confounding effects of step change invalidate the stationarity assumption of commonly used trend analysis methods such as the Mann-Kendall test technique, so previous studies have failed to explain inconsistencies between detected trends and observed large precipitation anomalies. The objectives of this study were to: 1) formulate a trend analysis approach that considers nonstationarity due to step changes; 2) use this approach to detect trends and extreme occurrences of precipitation in a mid-latitude Eurasian steppe watershed in north China; and 3) examine how runoff responds to precipitation trends in the study watershed. Our results indicate that annual precipitation underwent a marginal step jump around 1995. The significant annual downward trend after 1994 was primarily due to a decrease in summer rainfall; other seasons exhibited no significant precipitation trends. At a monthly scale, July rainfall after 1994 exhibited a significant downward trend, whereas precipitation in other months had no trend. The percentage of wet days also underwent a step jump around 1994 following a significant decreasing trend, although the precipitation intensity exhibited neither a step change nor any significant trend. However, both low- and high-frequency precipitation events in the study watershed occurred more often after than before 1994, probably as either a result or an indicator of climate change. In response to these precipitation changes, the study watershed had distinctly different precipitation-runoff relationships for observed annual precipitations of less than 300 mm, between 300 and 400 mm, and greater than 400 mm. This article is protected by copyright. All rights reserved.

Description: Hydrological models are useful tools to analyze present and future conditions of water quantity and quality. The integrated modeling of water and nutrients needs an adequate representation of the different discharge components. In common with many lowlands, groundwater contribution to the discharge in the North German lowlands is a key factor for a reasonable representation of the water balance especially in low flow periods. Several studies revealed that the widely used SWAT model performs poorly for low flow periods. This paper deals with the extension of the groundwater module of the SWAT model to enhance low flow representation. The current two-storage concept of SWAT was further developed to a three-storage-concept. This was realized due to modification of the groundwater module by splitting the active roundwater storage into a fast and a slow contributing aquifer. The results of this study show that the groundwater module with three storages leads to good prediction of the overall discharge especially for the recession limbs and the low flow periods. The improved performance is reflected in the signature measures for the mid segment (PBIAS: −2:4% vs. −15:9%) and the low segment (PBIAS: 14:8% vs. 46:8%) of the flow duration curve. The three-storage groundwater module is more process oriented than the original version due to the introduction of a fast and a slow groundwater flow component. The three-storage version includes a modular approach, since groundwater storages can be activated or deactivated independently for subbasin and HRU level. This article is protected by copyright. All rights reserved.

Description: The ecological condition and biodiversity values of floodplain wetlands are highly dependent on the hydrological connectivity of wetlands to adjacent rivers. This paper describes a method for quantifying connectivity between floodplain wetlands and the main rivers in a wet tropical catchment of northern Australia. We used a 1-D hydrodynamic model to simulate time-varying water depths across the stream network (i.e. rivers, streams and man-made drains). The timing and duration of connectivity of seven wetlands (4 natural and 3 artificial) with the two main rivers in the catchment were then calculated for different hydrological conditions. Location and areal extent of the wetlands and the stream network were identified using high resolution laser altimetry (LiDAR) and these data formed key inputs to the hydrodynamic model. The model was calibrated using measured water depths and discharges across the floodplain. An algorithm was developed to identify contiguous water bodies at daily time steps and this gave the temporal history of connection and disconnection between wetlands and the rivers. Simulation results show that connectivity of individual wetlands to both rivers varies from 26 to 365 days during an average hydrological condition. Location, especially proximity to a main river, and wetland type (natural stream or artificial drain) were identified as key factors influencing these levels of connectivity. Some natural wetlands maintain connection with the river for most or all of the year, whereas the connectivity of some artificial wetlands varies from 26 to 36 days according to their patterns of network connection to adjacent rivers – a result that has important implications for the accessibility of these types of wetland to aquatic biota. Using readily available river gauge data, we also show how connectivity modelling can be used to identify periods when connectivity has fallen below critical thresholds for fish movement. These connectivity patterns within the floodplain network are central to the setting of river flows that will meet environmental requirements for biota that use floodplain wetlands during their life history. This article is protected by copyright. All rights reserved.

Description: Pumping tests are one of the most commonly used in-situ testing techniques for assessing aquifer hydraulic properties. Numerous researches have been conducted to predict the effects of aquifer heterogeneity on the groundwater levels during pumping tests. The objectives of the present work were 1) to predict drawdown conditions and estimate aquifer properties during pumping tests undertaken in radially symmetric heterogeneous aquifers and 2) to identify a method for assessing the transmissivity field along the radial coordinate in radially symmetric and fully heterogeneous transmissivity fields. The first objective was achieved by expanding an existing analytical drawdown formulation which was valid for a radially symmetric confined aquifer with two concentric zones around the pumping well to an N concentric zone confined aquifer having a constant transmissivity value within each zone. The formulation was evaluated for aquifers with three and four concentric zones to assess the effects of the transmissivity field on the drawdown conditions. The specific conditions under which aquifer properties could be identified using traditional methods of analysis were also evaluated. The second objective was achieved by implementing the Inverse Solution Algorithm (ISA) which was developed for petroleum reservoirs to groundwater aquifer settings. The results showed that the drawdown values are influenced by a volumetric integral of a weighting function and the transmissivity field within the cone of depression. The weighting function migrates in tandem with the expanding cone of depression. The ability of the ISA to predict radially symmetric and log-normally distributed transmissivity fields was assessed against analytical and numerical benchmarks. The results of this investigation indicated that the ISA method is a viable technique for evaluating the radial transmissivity variations of heterogeneous aquifer settings. This article is protected by copyright. All rights reserved.

Description: Stormwater management increasingly recognises the need to emulate, to the maximum extent possible, the flow regime of receiving waters in their pre-development state. Hydrological models play a central role in assessing the catchment-scale impacts of alternative stormwater management strategies. However, because of the complexity of physical processes involved in urban hydrology, particularly subsurface flows, the predictive performance of such models is often low. We investigated how the structure of hydrological models influenced the prediction of urbanisation and stormwater management impacts on baseflow. We calibrated three conceptual models of the same reference catchment and compared the modelled flow regime from different stormwater management scenarios, using each of the three model structures. Scenarios were assessed using six metrics, characterising the whole streamflow regime and in particular baseflow. While the three models of the reference catchment represented the observed hydrograph well, the most complex structure, developed using a thorough diagnostic of the catchment behaviour, better captured the change in hydrological regime during dry years. Predictions of baseflow changes due to urbanisation varied significantly according to the model structure. Similarly, the models showed distinct responses to the stormwater management scenarios applied, especially for scenarios involving infiltration of stormwater at source. Our results confirm the importance of predicting the consequences of land use changes with conceptual models that are consistent with the hydrological behaviour of the study catchment. Future work should help quantify the uncertainties due to model structure, and thus provide practical guidance to the use of catchment models for assessing stormwater management strategies. This article is protected by copyright. All rights reserved.

Description: Although fire-induced soil water repellency (SWR) and its effects on soil hydrology and geomorphology have been studied with detail, very few studies have considered the effect of rock fragments resting on the soil surface or partly embedded in soil. In this research, we have studied the effect of rock fragments on the strength and spatial distribution of fire-induced SWR at different fire severities. A fire-affected area was selected for this experiment and classified into different zones according to fire severity (unburned, low, moderate and high) and rock fragment cover (low, 〈20%, and high, 〉60%). During 7 days after fire, SWR and infiltration rates were assessed in the soil surface covered by individual rock fragments and in the midpoint between two adjacent rock fragments (with maximum spacing of 20 cm). SWR increased with fire severity. Rock fragments resting on the soil surface increased the heterogeneity of the spatial distribution of fire-induced SWR. SWR increased significantly with rock fragment cover in bare areas under moderate and high fire severity, but quantitatively important changes were only observed under high fire severity. In areas with a low rock fragment cover, water repellency from soil surfaces covered by rock fragments increased relative to bare soil surfaces, with increasing SWR. In areas with a high rock fragment cover, SWR increased significantly from non-covered to covered soil surfaces only after low-severity burning. Rock fragment cover did not affect infiltration rates, although it decreased significantly in soil surfaces after high-severity burning in areas under low and high rock fragment cover. This article is protected by copyright. All rights reserved.

Description: Extremely low water level events have frequently occurred in the middle-lower Yangtze River (MLYR) in recent years (2006–2011). Most of these drought events coincided with the initial operation of the Three Gorges Dam (TGD). The TGD was therefore the focus of controversy about the causes of the hydrological droughts of the rivers and lakes of the region. We quantified the effects of the TGD's operation on water levels from 2006 to 2011 using a newly developed hydrodynamic model. The operation of the TGD significantly exacerbated the severe hydrological droughts that occurred in late September to November due to water impoundment, but it increased water levels from April to early June in the MLYR due to the drawdown of TGD water levels. Evidence suggests that the recent extremely low water levels were mainly due to the remarkable decline in inflows to the MLYR resulting from precipitation changes and possible human activities. Nevertheless, it should be noted that the effects of the TGD on downstream rivers and lakes will be intensified in the foreseeable future. This article is protected by copyright. All rights reserved.

Description: Water availability is one of the key environmental factors that control ecosystem functions in temperate forests. Changing climate is likely to alter the ecohydrology and other ecosystem processes that affect forest structures and functions. We constructed a multi-year water budget (2004–2010) and quantified environmental controls on an evapotranspiration (ET) in a 70-year-old mixed-oak woodland forest in northwest Ohio, USA. ET was measured using the eddy-covariance (EC) technique along with precipitation (P), soil volumetric water content (VWC), and shallow groundwater table fluctuation. Three biophysical models were constructed and validated to calculate potential ET (PET) for developing predictive monthly ET models. We found that the annual variability in ET was relatively stable and ranged from 578 mm in 2009 to 670 mm in 2010. In contrast, ET/P was more variable and ranged from 0.60 in 2006 to 0.96 in 2010. Mean annual ET/PET_FAO was 0.64 while mean annual PET_FAO/P was 1.15. Annual ET/PET_FAO was relatively stable, varying from 0.60 in 2005 to 0.72 in 2004. Soil water storage and shallow groundwater recharge during the non-growing season were essential in supplying ET during the growing season when ET exceeded P. Spring leaf area index (LAI), summer photosynthetically active radiation (PAR), and autumn and winter air temperatures (T a ) were the most significant controls of monthly ET. Moreover, LAI regulated ET during the whole growing season and higher temperatures increased ET even during dry periods. Our empirical modeling showed that the interaction of LAI and PET explained 〉90% of the variability in measured ET. Altogether we found that increases in T a and shifts in P distribution are likely to impact forest hydrology by altering shallow groundwater fluctuations, soil water storage, and ET and, consequently, alter the ecosystem functions of temperate forests. This article is protected by copyright. All rights reserved.

Description: Subsurface stormflow is thought to occur mainly in humid environments with steep terrains. However in semi-arid areas, preferential flow through macropores can also result in a significant contribution of subsurface stormflow to catchment runoff for varying catchment conditions. Most hydrological models neglect this important subsurface preferential flow. Here we use the process-oriented hydrological model Hillflow-3D, which includes a macropore flow approach, to simulate rainfall-runoff in the semi-arid Parapuños catchment in Spain, where macropore flow was observed in previous research. The model was extended for this study to account for sorptivity under very dry soil conditions. The results of the model simulations with and without macropore flow are compared. Both model versions give reasonable results for average rainfall situations, although the approach with the macropore concept provides slightly better results. The model results for scenarios of extreme rainfall events (〉 13.3 mm per 30 min) however show large differences between the versions with and without macropores. These model results compared to measured rainfall-runoff data show that the model with the macropore concept is better. Our conclusion is that preferential flow is important in controlling surface runoff in case of specific, high intensity rainfall events. Therefore preferential flow processes must be included in hydrological models where we know that preferential flow occurs. Hydrological process models with a less detailed process description may fit observed average events reasonably well but can result in erroneous predictions for more extreme events. This article is protected by copyright. All rights reserved.

Description: The central route of the South-North Water Transfer Project (CTP) is designed to divert approximately 9.5 billion m 3 of water per year from the Han River, a major tributary of the Yangtze River, to the Hai River Basin in the North China. The main purpose of this study is to assess the impact of CTP on groundwater table in the Hai River basin. Our study features a large-scale distributed hydrological model that couples a physically based groundwater module (GWM), which is subbasin-based, with a conceptual surface water module (SWM), which is grid-based. There are several grids in each subbasin and water exchange among grid are considered. Our model couples SWM and GWM, and calculates human water use at the same time. The simulation results indicate that even with the water supply by CTP, the groundwater table will continue to decline in the Hai River basin. However, the CTP water can evidently reduce the decline rate, helping alleviate groundwater over-exploitation in Hai River region. This article is protected by copyright. All rights reserved.

Description: Complexity-reduction modelling can be useful for increasing the understanding of how the climate affects basin soil moisture response upon historical times not covered by detailed hydrological data. For this purpose, here is presented and assessed an empirical regression-based model, the European Soil Moisture Empirical Downscaling (ESMED), in which different climatic variables, easily available on the web, are addressed for simplifying the inherent complexity in the long-time studies. To accommodate this simplification, the Palmer Drought Severity Index, the precipitation, the elevation and the geographical location were used as input data in the ESMED model for predicting annual soil moisture budget. The test area was a large region including central Europe and Mediterranean countries and the spatial resolution was initially set at 50 km. ESMED model calibration was made according with the soil moisture values retrieved from the Terrestrial Water Budget Data archive by selecting randomly 285 grid points (out of 2606). Once parameterized, ESMED model was performed at validation stage both spatially and temporally. The spatial validation was made for the grid points not selected in the calibration stage while the comparison with the soil moisture outputs of the GLDAS-NOAH10 simulations upon the period 1950-2010 was carried out for the temporal validation. Moreover, ESMED results were found to be in good agreement with a root-zone soil moisture product obtained from active and passive microwave sensors from various satellite missions. ESMED model was thus found to be reliable for both the temporal and spatial validation and, hence, it might represent a useful tool to characterize the long-term dynamics of soil moisture-weather interaction. This article is protected by copyright. All rights reserved.

Description: Tundra snow cover is important to monitor as it influences local, regional, and global scale surface water balance, energy fluxes, as well as ecosystem and permafrost dynamics. Observations are already showing a decrease in spring snow cover duration at high latitudes but the impact of changing winter season temperature and precipitation on variables such as snow water equivalent (SWE) is less clear. A multi-year project was initiated in 2004 with the objective to quantify tundra snow cover properties over multiple years at a scale appropriate for comparison with satellite passive microwave remote sensing data and regional climate and hydrological models. Data collected over seven late winter field campaigns (2004 to 2010) show the patterns of snow depth and SWE are strongly influenced by terrain characteristics. Despite the spatial heterogeneity of snow cover, several inter-annual consistencies were identified. A regional average density of 0.293 g/cm 3 was derived and shown to have little difference with individual site densities when deriving SWE from snow depth measurements. The inter-annual patterns of SWE show that despite variability in meteorological forcing, there were many consistent ratios between the SWE on flat tundra and the SWE on lakes, plateaus, slopes. A summary of representative inter-annual snow stratigraphy from different terrain categories is also presented. This article is protected by copyright. All rights reserved.

Description: The main objective of this paper is to provide comparative-quantitative examinations on capabilities of 2DH and pseudo-3D modelling approaches for simulating spatial and temporal variability of the flow and salinity in Lake Urmia, Iran. The water quality in the lake has been an environmentally important subject partly because this shallow hyper saline aquatic ecosystem is considered to be one of the largest natural habitats of a unique multi-cellular organism, Artemia urmiana . This brine shrimp is the major food source for many of the protected and rare shorebirds that visit the lake. Artemia urmiana can grow and survive in certain ranges of salinity and their disappearance could lead to an alteration of existing equilibria. The lake has also experienced considerable man-made changes during the past 3 decades. A newly built crossing embankment almost divided the lake into two northern and southern halves. A relatively small opening of 1.25kmin the new embankment provides water connections between the two halves. As a result, the flow and salinity regimes have been significantly changed. This might have had adverse serious impacts on the lake ecosystem. In the current study the 2DH hydrodynamic model has been found to provide reasonable predictions for the flow regime in the lake, while its salinity predictions have not been consistent with the field observations.Thepseudo-3D model has produced results fairly close to the salinity measurements and its temporal and spatial variations. The pseudo-3Dmodelhas been used for evaluating the embankment effects on the lake hydrodynamics and on the salinity conditions. The effectiveness of introducing a different number or length of openings in the embankment for restoring the pre-embankment conditions has also been examined. These remedy options have been found not to offer substantial improvements to the lake existing ecosystem. This article is protected by copyright. All rights reserved.

Description: Subsurface flow and heat transport near Freienbrink, NE Germany, was simulated in order to study groundwater-surface water exchange between a floodplains aquifer and a section of the lowland River Spree and an adjacent oxbow. Groundwater exfiltration was the dominant process and only fast surface water level rises resulted in temporary infiltration into the aquifer. The main groundwater flow paths are identified based on a 3D groundwater flow model. To estimate mass fluxes across the aquifer-surface water interfaces, a 2D flow and heat transport modelling approach along a transect of 12 piezometers was performed. Results of steady-state and transient water level simulations show an overall high accuracy with a Spearman coefficient ρ = 0.9996 and RMSE = 0.008 m. Based on small groundwater flow velocities of about 10 -7 to 10 -6  ms -1 mean groundwater exfiltration rates of 233 l m -2 d -1 are calculated. Short periods of surface water infiltration into the aquifer do not exceed 10 days and the infiltration rates are in the same range. The heat transport was modelled with slightly less accuracy (ρ = 0.8359 and RMSE = 0.34 °C). In contrast to the predominant groundwater exfiltration, surface water temperatures determine the calculated temperatures in the upper aquifer below both surface water bodies down to 10 m during the whole simulation period. These findings emphasize prevailing of heat conduction over advection in the upper aquifer zones, which seems to be typical for lowland streams with sandy aquifer materials and low hydraulic gradients. Moreover, this study shows the potential of coupled numerical flow and heat transport modelling to understand groundwater-surface water exchange processes in detail. This article is protected by copyright. All rights reserved.

Description: Particles eroded from hillslopes and exported to rivers are recognized to be composite particles of high internal complexity. Their architecture and composition is known to influence their transport behaviour within the water column relative to discrete particles. To-date, hillslope erosion studies consider aggregates to be stable once they are detached from the soil matrix. However, lowland rivers and estuaries studies often suggest that particle structure and dynamics are controlled by flocculation within the water column. In order to improve the understanding of particle dynamics along the continuum from hillslopes to the lowland river environment, soil particle behaviour was tested under controlled laboratory conditions. Seven flume erosion and deposition experiments, designed to simulate a natural erosive event, and five shear cell experiments were performed using three contrasting materials: two of them were poorly developed and as such can not be considered as soils, whilst the third one was a calcareous brown soil. These experiments revealed that soil aggregates were prone to disaggregation within the water column and that flocculation may affect their size distribution during transport. Large differences in effective particle size were found between soil types during the rising limb of the bed shear stress sequence. Indeed, at the maximum applied bed shear stress, the aggregated particles median diameter was found to be three times larger for the well-developed soil than for the two others. Differences were smaller in the falling limb, suggesting that soil aggregates underwent structural changes. However, characterization of particles strength parameters showed that these changes did not fully turn soil aggregates into flocs, but rather into hybrid soil aggregate-floc particles. This article is protected by copyright. All rights reserved.

Description: Submarine groundwater discharges (SGD) were investigated in a marine watershed in south-eastern Korea using water budget analysis and a 222 Rn mass balance model. Multi-layered TOPMODEL added hydrological assumption was used to estimate groundwater components in the water budget analysis. Field observations of soil moisture, rainfall, runoff and groundwater fluctuations were used for calibration and validation of the hydrologic model. Based on observed hydrological data and terrain analyses, parameters for the hydrologic model were delineated and used to describe several hydrologic responses in the watershed. SGD estimations by 222 Rn mass balance method were also performed at Il-Gwang bay in July, 2010, and May, June, July and Nov. 2011. The estimated groundwater through hydrologic modeling and water balance analysis was 1.3x10 6  m 3 /year, which rapidly increased during typhoon season due to heavy rainfall and permeable geologic structure. The estimated groundwater was approximately 3.7-27.1 % of SGD as evaluated by 222 Rn mass balance method ranges 3.44 and 17.45 m 3  m -2  year -1 . Even though SGD is predominantly influenced by tide fluctuation, the head gradient (difference) from hydrologic processes associated with heavy rainfalls can also extra significant influences. Comprehensive understanding of SGD evaluation can be improved through a simultaneous application of both these approaches. This article is protected by copyright. All rights reserved.

Description: Fallout radionuclides 137 Cs and 210 Pb are well established as tracers of surface and sub-surface soil erosion contributing sediment to river systems. However without additional information it has not been possible to distinguish sub-surface soil erosion sources. Here we use the fallout radionuclide 7 Be (half-life 53 days) in combination with 137 Cs and excess 210 Pb to trace the form of erosion contributing sediment in large river catchments in eastern Australia; the Logan River (area 3,700 km 2 ), Bowen River (9,400 km 2 ) and Mitchell River (4,700 km 2 ). We show that the combination of 137 Cs, excess 210 Pb and 7 Be can discriminate horizontally-aligned sub-surface erosion sources (rilled and scalded hillslopes and the floors of incised drainage lines and gully ‘badland’ areas) from vertical erosion sources (channel banks and gully walls). Specifically, sub-surface sources of sediment eroded during high rainfall and high river flow events have been distinguished by the ability of rainfall-derived 7 Be to label horizontal soil surfaces, but not vertical. Our results indicate that in the two northern catchments erosion of horizontal sub-surface soil sources contributed almost as much fine river sediment as vertical channel banks, and several times the contribution of hillslope topsoils. This result improves on source discrimination provided previously and indicates that in some areas erosion of hillslope soils may contribute significantly to sediment yield, but not as topsoil loss. We find that in north-eastern Australia scalded areas on hillslopes and incising drainage lines may be sediment sources of comparable importance to vertical channel banks. Previous studies have used the combination of 137 Cs, excess 210 Pb and 7 Be to estimate soils losses at the hillslope scale. Here we show that with timely and judicious sampling of soil and sediment during and immediately after high flow events 7 Be measurements can augment fallout 137 Cs and 210 Pb to provide important erosion source information over large catchments. This article is protected by copyright. All rights reserved.

Description: The hydrology of tropical dry forests have been poorly characterised when compared to their humid temperate and wet tropical counterparts. Despite accounting for more than 42% of all tropical forests and roughly 19% of the Earth's total forest, tropical dry forest represent less than 1% of the forest hydrology literature. The need for substantial hydrological research in tropical dry forests is extremely important, given that many tropical dry forest regions are currently water stressed due to high population densities and rapid land use change. Furthermore, future climate change scenarios are expected to have significant implications for the hydrological functioning of these catchments and will likely enhance pressures on already limited water resources. This paper provides an overview on the state of hydrological knowledge, particularly runoff generation, of tropical dry forests. We further highlight the research gaps and identify research priorities for tropical dry forests, and issue a call for increased hydrological research efforts in these forests. This article is protected by copyright. All rights reserved.

Description: We use two hydrological models of varying complexity to study the Juncal River Basin in the Central Andes of Chile with the aim to understand the degree of conceptualization and the spatial structure that are needed to model present and future streamflow. We use a conceptual semi-distributed model based on elevation bands (WEAP), frequently used for water management, and a physically oriented, fully-distributed model (TOPKAPI-ETH) developed for research purposes mainly. We evaluate the ability of the two models to reproduce the key hydrological processes in the basin with emphasis on snow accumulation and melt, streamflow and the relationships between internal processes. Both models are capable of reproducing observed runoff and the evolution of MODIS snow cover adequately. In spite of WEAP's simple and conceptual approach for modeling snowmelt, its lack of glacier representation and snow gravitational redistribution as well as a proper routing algorithm, this model can reproduce historical data with similar goodness-of-fit as the more complex TOPKAPI-ETH. We show that the performance of both models can be improved by using measured precipitation gradients of higher temporal resolution. In contrast to the good performance of the conceptual model for the present climate, however, we demonstrate that the simplifications in WEAP lead to error compensation which results in different predictions in simulated melt and runoff for a potentially warmer future climate. TOPKAPI-ETH, using a more physical representation of processes, depends less on calibration and thus is less subject to a compensation of errors through different model components. Our results show that data obtained locally in ad-hoc short-term field campaigns are needed to complement data extrapolated from long-term records for simulating changes in the water cycle of high elevation catchments, but that these data can only be efficiently used by a model applying a spatially distributed physical representation of hydrological processes. This article is protected by copyright. All rights reserved.

Description: Land use/cover (LULC) and climate change are two main factors affecting watershed hydrology. In this paper, individual and combined impacts of LULC and climate change on hydrologic processes were analyzed applying the model Soil and Water Assessment Tool (SWAT) in a coastal Alabama watershed in USA. Temporally and spatially downscaled Global Circulation Model (GCM) outputs predict a slight increase in precipitation in the study area, which is also projected to experience substantial urban growth in the future. Changes in flow frequency and volume in the 2030s (2016–2040) compared to a baseline period (1984–2008) at daily, monthly and annual time scales were explored. A redistribution of daily streamflow is projected when either climate or LULC change was considered. High flows are predicted to increase, while low flows are expected to decrease. Combined change effect results in a more noticeable and uneven distribution of daily streamflow. Monthly average streamflow and surface runoff are projected to increase in spring and winter, but especially in fall. LULC change does not have a significant effect on monthly average streamflow, but the change affects partitioning of streamflow, causing higher surface runoff and lower baseflow. The combined effect leads to a dramatic increase in monthly average streamflow with a stronger increasing trend in surface runoff and decreasing trend in baseflow. This article is protected by copyright. All rights reserved.

Description: The fresh groundwater lenses (FGLs) of small islands can be highly vulnerable to climate change impacts, including sea-level rise (SLR). Many real cases of atoll or sandy islands involve two-layer hydrogeological conceptualizations. In this paper, the influential factors that effect FGL in two-layer small islands subject to SLR are investigated. An analytical solution describing FGLs in circular islands, composed of two geological layers, is developed for the simplified case of steady-state and sharp-interface conditions. An application of the developed model is demonstrated to estimate the FGL thickness of some real-world islands by comparison with existing FGL thickness data. Furthermore, numerical modeling is applied to extend the analysis to consider dispersion effects and to confirm comparable results for both cases. Sensitivity analyses are used to assess the importance of land-surface inundation (LSI) caused by SLR, relative to other parameters (i.e. thickness of aquifer layers, hydraulic conductivity, recharge rate, and land-surface slope) that influence the FGL. Dimensionless parameters are used to generalize the findings. The results demonstrate that LSI has a considerable impact on a FGL influenced by SLR, as expected, although the FGL volume is more sensitive to recharge, aquifer thickness and hydraulic conductivity than SLR impacts, considering typical parameter ranges. The methodology presented in this study provides water resource managers with a rapid-assessment tool for evaluating the likely impacts of SLR and accompanying LSI on FGLs. This article is protected by copyright. All rights reserved.

Description: Mapping groundwater discharge zones at broad spatial scales remains a challenge, particularly in data sparse regions. We applied a regional scale mapping approach based on thermal remote sensing to map discharge zones in a complex watershed with a broad diversity of geological materials, land cover and topographic variation situated within the Prairie Parkland of northern Alberta, Canada. We acquired winter thermal imagery from the USGS Landsat archive to demonstrate the utility of this data source for applications that can complement both scientific and management programs. We showed that the thermally determined potential discharge areas were corroborated with hydrological (spring locations) and chemical (conservative tracers of groundwater) data. This study demonstrates how thermal remote sensing can form part of a comprehensive mapping framework to investigate groundwater resources over broad spatial scales. This article is protected by copyright. All rights reserved.

Description: Obtaining representative meteorological data for watershed-scale hydrological modeling can be difficult and time consuming. Land-based weather stations do not always adequately represent the weather occurring over a watershed, because they can be far from the watershed of interest, have gaps in their data series, or recent data is not available. This study presents a method for using the Climate Forecast System Reanalysis (CFSR) global meteorological data set to obtain historical weather data and demonstrates the application to modeling five watersheds representing different hydroclimate regimes. CFSR data are available globally for each hour since 1979 at a 38 km resolution. Results show that utilizing the CFSR precipitation and temperature data to force a watershed model provide stream discharge simulations that are as good as or better than models forced using traditional weather gaging stations, especially when stations are more than 10 km from the watershed. These results further demonstrate that adding CFSR data to the suite of watershed modeling tools provides new opportunities for meeting the challenges of modeling un-gaged watersheds and advancing real-time hydrological modeling. This article is protected by copyright. All rights reserved.

Description: River water temperature is a very important variable in ecological studies, especially for the management of fisheries and aquatic resources. Temperature can impact on fish distribution, growth, mortality and community dynamics. River evaporation has been identified as an important heat loss and a key process in the thermal regime of rivers. However, its quantification remains a challenge, mainly due to the difficulty of making direct measurements. The objectives of this study were to characterize the evaporative heat flux at different scales (brook vs. river) and to improve the estimation of the evaporative heat flux in a stream temperature model at the hourly timescale. Using a mass-balance approach with floating minipans, we measured river evaporation at an hourly timescale in a medium-sized river (Little Southwest Miramichi) and a small brook (Catamaran Brook) in New Brunswick, Canada. With these direct measurements of evaporation, we developed mass transfer equations to estimate hourly evaporation rates from microclimate conditions measured 2 m above the stream. During the summer 2012, river evaporation was more important for the medium-sized river with a mean daily evaporation rate of 3.0 mm day -1 in the Little Southwest Miramichi River compared to 1.0 mm day -1 in Catamaran Brook. Evaporation was the main heat loss mechanism in the two studied streams and was responsible for 42 % of heat losses in the Little Southwest Miramichi River and 34 % of heat losses in Catamaran Brook during the summer. This article is protected by copyright. All rights reserved.

Description: Although catchment storage is an intrinsic control on the rainfall-runoff response of streams, direct measurement remains a major challenge. Coupled models that integrate long-term hydrometric and isotope tracer data are useful tools that can provide insights into the dynamics of catchment storage and the volumes of water involved. In this study, we use a tracer-aided hydrological model to characterize catchment storage as a dynamic control on system function related to streamflow generation, which also allows direct estimation of the non-stationarity of water ages. We show that in a wet Scottish upland catchment dominated by runoff generation from riparian peats (histosols) with high water storage, non-stationarity in water age distributions are only clearly detectable during more extreme wet and dry periods. This is explained by the frequency and longevity of hydrological connectivity and the associated relative importance of flow paths contributing younger or older waters to the stream. Generally, these saturated riparian soils represent large mixing zones that buffer the time variance of water age and integrate catchment-scale partial mixing processes. Although storage simulations depend on model performance, which is influenced by input variability and the degree of isotopic damping in the stream, a longer-term storage analysis of this model indicates a system which is only sensitive to more extreme hydroclimatic variability. This article is protected by copyright. All rights reserved.

Description:    Natural hazard assessments are always subject to uncertainties due to missing knowledge about the complexity of hazardous processes as well as their natural variability. Decision-makers in the field of natural hazard management need to understand the concept, components, sources, and implications of existing uncertainties in order to reach informed and transparent decisions. Until now, however, only few hazard maps include uncertainty visualizations which would be much needed for an enhanced communication among experts and decision-makers in order to make informed decisions possible. In this paper, an analysis of how uncertainty is currently treated and communicated by Swiss natural hazards experts is presented. The conducted expert survey confirmed that the communication of uncertainty has to be enhanced, possibly with the help of uncertainty visualizations. However, in order to visualize the spatial characteristics of uncertainty, existing uncertainties need to be quantified. This challenge is addressed by the exemplary simulation of a snow avalanche event using a deterministic model and quantified uncertainties with a sensitivity analysis. Suitable visualization methods for the resulting spatial variability of the uncertainties are suggested, and the advantages and disadvantages of their implementation in an interactive cartographic information system are discussed. Content Type Journal Article Pages 1-17 DOI 10.1007/s11069-011-9864-y Authors Melanie Kunz, Institute of Cartography and Geoinformation, ETH Zurich, 8093 Zurich, Switzerland Adrienne Grêt-Regamey, Planning of Landscape and Urban Systems, Institute for Spatial and Landscape Planning, ETH Zurich, 8093 Zurich, Switzerland Lorenz Hurni, Institute of Cartography and Geoinformation, ETH Zurich, 8093 Zurich, Switzerland Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    The present work addresses the problem of lack of coordination between policies and actors with joint competence for risk management, i.e., civil protection, spatial planning, and sectoral planning (e.g., forest policy in the case of forest fire risk). Spatial planning in particular is assigned a minor or no role at all though it might perfectly operate as the coordinating policy platform; the reason is that spatially relevant analysis and policy guidance is an omnipresent component of the risk management cycle. However, disconnected risk relevant policies turning a blind eye to spatial planning might cause several adverse repercussions: Breaks in the response-preparedness-prevention-remediation chain (which should function as a continuum), minimal attention to prevention, risk expansion and growth instead of mitigation, lack of synergies between involved actors as well as duplicated or even diverging measures and funding. The authors bear witness to the above suggestions by examining three cases of European (regional and local) risk management systems faced with failures when confronting natural hazards (floods and forest fires). These three systems are embedded in different types of political-administrative structures, namely those of the city of Dortmund (Germany) facing floods, Eastern Attica region (Greece), and Lazio Region (Italy) facing forest fires. Content Type Journal Article Pages 1-30 DOI 10.1007/s11069-011-9843-3 Authors Kalliopi Sapountzaki, Department of Geography, Harokopio University of Athens, Athens, Greece Sylvia Wanczura, Institute of Fire Service and Rescue Technologies, Fire Department of Dortmund, Dortmund, Germany Gabriella Casertano, Regione Lazio, Dipartimento della Protezione Civile, Rome, Italy Stefan Greiving, Technical University of Dortmund, Dortmund, Germany Gavriil Xanthopoulos, National Agricultural Research Foundation, Institute of Mediterranean Forest Ecosystems and Forest Products Technology, Athens, Greece Floriana F. Ferrara, T6 Ecosystems srl, Rome, Italy Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    A practical, DEM-based practical method is proposed to enhance flood risk management in fluvial areas by quantifying relative risk as a function of vulnerability to inland and evacuation difficulty. Both measures are based mainly on the topography of the region, so the method does not require detailed data on the physical characteristics of the land. First, we use the deterministic 8-node method on a digital elevation map (DEM) to trace storm waterways. Second, we repeat the process on a reversed DEM to trace evacuation routes that avoid the waterways and zones dangerously close to the rivers. Finally, on the basis of such two flow lines of evacuee and storm water, we proposed the protocol to evaluate the flood risk at every point on the map taking into account both the minimum time required for floodwater to arrive and duration of an evacuation from that location. The time that must be allocated for safe evacuation is defined as the potential flood risk of evacuation (PFRE). The method is demonstrated on a fluvial area of the Kaki River in Nagaoka city, Japan. In addition, we illustrated the application of the PFRE map to divide the region into areas of greater or lesser evacuation urgency. Content Type Journal Article Pages 1-12 DOI 10.1007/s11069-011-9852-2 Authors Jinyoung Kim, Department of Urban Engineering, The University of Tokyo, Tokyo, 113-8656 Japan Yuji Kuwahara, Ibaraki University, Ibaraki, Japan Manish Kumar, Tezpur University, Napaam, Assam 784 028, India Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    Earthquake-induced landslides can cause a large number of casualties and great economic loss. Presently, research methods for studying landslides are largely based on a framework of solid mechanics. In this work, visualization software, known as Visual SPH, is developed in Visual Basic 6.0, and can be used to analyze the fluidized movement of landslides based on smoothed particle hydrodynamics (SPH). First, the accuracy of the software is verified through simulating a benchmark problem of a dam-break. Then, fluidized movement analysis of the Tangjiashan landslide is carried out using this software. The simulation derives plots of displacement versus time, reproduces the entire flow process of the Tangjiashan landslide, and determines the run-out, which coincides exactly with the characteristics of flow-like landslides observed in the field. All of these are essential in the design of supporting structures and site selection for reconstruction in earthquake-prone regions. Content Type Journal Article Pages 1-14 DOI 10.1007/s11069-011-9859-8 Authors Yu Huang, Key Laboratory of Geotechnical and Underground Engineering of the Ministry of Education, Tongji University, Shanghai, 200092 China Zili Dai, Department of Geotechnical Engineering, Tongji University, 200092 Shanghai, China Weijie Zhang, Department of Geotechnical Engineering, Tongji University, 200092 Shanghai, China Zhiyi Chen, Key Laboratory of Geotechnical and Underground Engineering of the Ministry of Education, Tongji University, Shanghai, 200092 China Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    Marine overwash from the north a few centuries ago transported hundreds of angular cobbles and boulders tens to hundreds of meters southward from limestone outcrops in the interior of Anegada, 140 km east–northeast of Puerto Rico. We examined two of several cobble and boulder fields as part of an effort to interpret whether the overwash resulted from a tsunami or a storm in a location where both events are known to occur. One of the cobble and boulder field extends 200 m southward from limestone outcrops that are 300 m inland from the island’s north shore. The other field extends 100 m southward from a limestone knoll located 800 m from the nearest shore. In the two fields, we measured the size, orientation, and spatial distribution of a total of 161 clasts and determined their stratigraphic positions with respect to an overwash sand and shell sheet deposit. In both fields, we found the spacing between clasts increased southward and that clast long-axis orientations are consistent with a transport trending north–south. Almost half the clasts are partially buried in a landward thinning and fining overwash sand and none were found embedded in the shelly mud of a pre-overwash marine pond. The two cobble and boulder fields resemble modern tsunami deposits in which dispersed clasts extend inland as a single layer. The fields contrast with coarse clast storm deposits that often form wedge-shaped shore-parallel ridges. These comparisons suggest that the overwash resulted from a tsunami and not from a storm. Content Type Journal Article Pages 1-13 DOI 10.1007/s11069-011-9848-y Authors Steve Watt, U.S. Geological Survey, 400 Natural Bridges Drive, Santa Cruz, CA 95060, USA Mark Buckley, U.S. Geological Survey, 400 Natural Bridges Drive, Santa Cruz, CA 95060, USA Bruce Jaffe, U.S. Geological Survey, 400 Natural Bridges Drive, Santa Cruz, CA 95060, USA Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    Improvements in hurricane forecasts allowing for more timely evacuations from storm-surge zones are credited with reducing lethality of US landfalling hurricanes. The deadly reach of a hurricane, however, is not limited to storm-surge zones. About 80% of direct US hurricane fatalities since 1970 occurred outside of landfall counties, with most of these fatalities caused by inland flooding. We construct a geographic information system database combining the location and cause of fatalities, estimated wind speeds, and rainfall amounts for the entire track of the storm for landfalling US hurricanes between 1970 and 2007. We analyze the determinants of total fatalities and deaths due to freshwater drowning and wind. Inclusion of inland fatalities results in no downward trend in lethality over the period, in contrast to prior research. Local storm conditions significantly affect lethality, as one-inch and one-knot increases in rainfall and wind increase total fatalities by 28 and 4%. Rainfall significantly increases freshwater-drowning deaths and is insignificant for wind deaths, while the opposite relation holds for wind speed. While coastal counties do not exhibit a significantly higher amount of lethality risk versus inland counties for total or wind-driven fatalities, freshwater-drowning fatalities occur most frequently in inland counties along the center of the storm path and its outer county tiers as we have defined them. Content Type Journal Article Pages 1-19 DOI 10.1007/s11069-011-9849-x Authors Jeffrey Czajkowski, Department of Economics, Austin College, Sherman, TX, USA Kevin Simmons, Department of Economics, Austin College, Sherman, TX, USA Daniel Sutter, Department of Economics and Finance, University of Texas - Pan American, Edinburg, TX, USA Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    This study explored demographic variations among victims injured in the 2008 Wenchuan earthquake. The study analyzed records of 2,148 patients who were hospitalized with earthquake-related injuries in 14 hospitals between May 12 and May 21 of 2008. We tested the effects of gender and age on the incidence of fractures and amputations. Although more women than men were hospitalized as earthquake victims, the results indicated that a higher incidence of pelvic fracture in women was the only statistically significant injury with regard to gender. Age had a more complex effect on the rate and pattern of fractures. Children younger than 15 years old were the most likely to sustain skull fractures and amputations. The elderly were the most likely to be associated with limb fractures. Adult men were the most likely to suffer spine fractures and adult women rib fractures. The paper concludes with an interpretation of the results in the context of China’s rural conditions and discusses implications for post-disaster emergency operations and planning. Content Type Journal Article Pages 1-19 DOI 10.1007/s11069-011-9836-2 Authors Ying Cao, College of Architecture and Environment, Sichuan University, No. 24 South Section 1, Yihuan Road, Chengdu, 610065 China Nabil Kamel, School of Geographical Sciences and Urban Planning, Arizona State University, Coor Hall, 975 S. Myrtle Ave., Fifth Floor, P.O. Box 875302, Tempe, AZ 85287-5302, USA Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    Following the catastrophic and devastating Atlantic Hurricane seasons in 2004 and 2005, there has been increased interest in formulating planning directives and policy aimed at minimizing the societal impacts of future storms. Not all populations will evacuate an area forecast to be affected by a hurricane, so emergency managers must plan for these people who remain behind. Such planning includes making food, water, ice, and other provisions available at strategic locations throughout an affected area. Recent research has tackled problems related to humanitarian and relief goods distribution with respect to hurricanes. Experience shows that the torrential rains and heavy winds associated with hurricanes can severely damage transportation network infrastructure rendering it unusable. Scanning the literature on hurricane disaster relief provision, there are no studies that expressly consider the potential damage that may be caused to a transportation network by strong storms. This paper examines the impacts of simulated network failures on hurricane disaster relief planning strategies, using a smaller Florida City as an example. A relief distribution protocol is assumed where goods distribution points are set up in pre-determined locations following the passage of a storm. Simulation results reveal that modest disruptions to the transportation network produce marked changes in the number and spatial configuration of relief facilities. At the same time, the transportation network appears to be robust and is able to support relief service provision even at elevated levels of hypothesized disruption. Content Type Journal Article Pages 1-16 DOI 10.1007/s11069-011-9855-z Authors Mark W. Horner, The Florida State University, Tallahassee, FL, USA Michael J. Widener, University at Buffalo, 105 Wilkeson Quad, Buffalo, NY 14261-0001, USA Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    Although the risk of flooding poses a serious threat to the Dutch public, citizens are not very inclined to engage in self-protective behaviors. Current risk communication tries to enhance these self-protective behaviors among citizens, but is nonetheless not very successful. The level of citizens engaging in self-protective actions remains rather low. Therefore, this research strives to determine the factors that might enhance or lessen the intention to engage in self-protection among citizens. The study was a 2 (flood risk: high vs low) × 2 (efficacy beliefs: high vs low) between subject experiment. It was conducted to test how varying levels of flood risk and efficacy beliefs influence two different self-protective behaviors, namely information seeking and the intention to engage in risk mitigating or preventive behaviors. Furthermore, the relationship between information seeking and the intention to take self-protective actions was discussed. Results showed that high levels of flood risk lead to higher levels of both information seeking and the intention to engage in self-protective behaviors than low levels of flood risk. For efficacy beliefs, the same trend occurred. Also, results showed that information seeking seems to coincide with the intention to take preventive actions and acted as a mediator between the levels of perceived risk and efficacy and the intention to take self-protective actions. Content Type Journal Article Pages 1-16 DOI 10.1007/s11069-011-9845-1 Authors Milou Kievik, Saxion, University for Applied Sciences, PO Box 70.000, 7500 KB Enschede, The Netherlands Jan M. Gutteling, University of Twente, PO Box 217, 7500 AE Enschede, The Netherlands Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    Much progress has been made in the area of tropical cyclone prediction using high-resolution mesoscale models based on community models developed at National Centers for Environmental Predication (NCEP) and National Center for Atmospheric Research (NCAR). While most of these model research and development activities are focused on predicting hurricanes in the Atlantic and Eastern Pacific domains, there has been much interest in using these models for tropical cyclone prediction in the North Indian Ocean region, particularly for Bay of Bengal storms that are known historically causing severe damage to life and property. In this study, the advanced operational hurricane modeling system developed at NCEP, known as the Hurricane Weather Research and Forecast (HWRF) model, is used to simulate two recent Bay of Bengal tropical cyclones—Nargis of November 2007 and Sidr of April 2008. The advanced NCEP operational vortex initialization procedure is adapted for simulating these Bay of Bengal tropical cyclones. Two additional regional models, the NCAR Advanced Research WRF and NCAR/Penn State University Mesoscale Model version 5 (MM5) are also used in simulating these storms. Results from these experiments highlight the superior performance of HWRF model over other models in predicting the Bay of Bengal cyclones. These results also suggest the need for a sophisticated vortex initialization procedure in conjunction with a model designed exclusively for tropical cyclone prediction for operational considerations. Content Type Journal Article Pages 1-19 DOI 10.1007/s11069-011-9839-z Authors D. V. Bhaskar Rao, TLGVRC, Jackson State University, Jackson, MS 39217, USA Vijay Tallapragada, NCEP/EMC, 5200 Auth Road, Camp Springs, MD 20746, USA Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    Landslides interacting with large infrastructures represent a major problem for the economy, society as a whole, and the safety of workers. Continuous monitoring for 23 months using an integrated platform with a ground-based SAR interferometer (GB-InSAR), a weather station, and an automatic camera gave us the opportunity to analyze the response of an unstable slope to the different phases of work. The deformational behavior of both the natural slope and the man-made structures was recorded and interpreted in relation to the working stages and the rainfall conditions during the whole monitoring period. A typical pattern of displacement was identified for shallow landslides, debris produced by the excavation and gabions, metallic walls, and anchored bulkheads. Furthermore, insights into the dynamics and behavior of the slope and the man-made structures that interact with the landslide were obtained. Extreme rainfall is the main trigger of shallow landslides and gabion deformations, while anchored bulkheads are less influenced by rainfalls. Movement of debris that is produced by excavations and temporary metallic barrier deformation are closely related to each other. The herein proposed monitoring platform is very efficient in monitoring unstable slopes that are affected by human activities. Moreover, the recorded patterns of displacement in the slope and the man-made structures can be used as reference data for similar studies and engineering designs. Content Type Journal Article Pages 1-20 DOI 10.1007/s11069-011-9840-6 Authors Francesca Bozzano, CERI, Research Centre on Prevention, Prediction and Control of Geological Risks, P.zza U. Pilozzi 9, 00038 Valmontone, Italy Ivan Cipriani, Dipartimento di Scienze della Terra, “Sapienza” Università di Roma, P.le Aldo Moro 5, 00185 Rome, Italy Paolo Mazzanti, CERI, Research Centre on Prevention, Prediction and Control of Geological Risks, P.zza U. Pilozzi 9, 00038 Valmontone, Italy Alberto Prestininzi, CERI, Research Centre on Prevention, Prediction and Control of Geological Risks, P.zza U. Pilozzi 9, 00038 Valmontone, Italy Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    In the midst of the ever-increasing natural and human-induced disasters, where many of the preparedness and mitigation measures show inefficiencies, there is narrow margin for decision-makers to make mistakes by misallocating budgets, designing infeasible reconstruction plans, and in other terms, making decisions not in line with the public preferences. In particular, public participation in post-disaster measures seems undoubtedly necessary to reduce the possible economic, social, political, and cultural conflicts around the stressful community after a major disaster. This paper aims at evaluating the role of public participation in increasing the reconstruction phase efficiency through a case study of the reconstruction process in Bam, a southeastern Iranian city, after the 2003 earthquake. It is attempted to identify the major motivators of the public participation through a combination of quantitative and qualitative studies. Statistical data are generated through a set of questionnaires being filled by a number of 200 randomly selected survivors. The numerical results were then discussed through the Focus Group technique sessions to determine the main contributors to the public participation. It is later found that the answers are found among the performance of the reconstruction authorities, financial policies, emotional resiliency of the survivors, public information mechanisms, public satisfaction, the pace of reconstruction, and temporary housing policies. Content Type Journal Article Pages 1-16 DOI 10.1007/s11069-011-9842-4 Authors Babak Omidvar, Graduate Faculty of Environment, University of Tehran, P.O. Box 14155-6135, Tehran, Iran Hossein Zafari, Housing Foundation of Islamic Revolution, Tehran, Iran Mehdi Khakpour, Graduate Faculty of Environment, University of Tehran, P.O. Box 14155-6135, Tehran, Iran Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    After the 2004 Sumatra–Andaman tsunamigenic earthquake, waters from the ocean moved upstream along rivers, bays, harbors, and lagoons and inundated many coastal and inland locations in the southern, eastern, and northern parts of Sri Lanka. The tsunami waters were observed to move upwards inland and then recede downwards to the ocean after varying inundation periods in different coastal areas. Subsequent massive tsunami waves came with the wave height varying from 3 to 8 m inland with speed of about 30–40 kmph. The oceanic waves carrying heterogeneous sediments with water deposited them in coastal as well as inland locations about 1 km from the present coastline. Given the chaotic nature of tsunami oceanic waves, pre-tsunami deposits, such as beach sands, debris from coral reefs and buildings, parts of vehicles and ships, and tree trunks are found incorporated in authentic tsunami sediments. Thus, the texture, structure, and composition of sediments deposited by tsunami waters differed from one location to another. Therefore, in identifying paleo-tsunami sediments, care was taken to compare them with diagnostic unmixed uncontaminated recent tsunami sediments having characteristic textures and marine microfossil assemblages, such as foraminifera, radiolarians, and diatoms where preserved in coastal depressions. The radiocarbon ages of the carbonate and the organic fractions of these sediments are stratigraphically inconsistent, indicating mixing of sediments by the tsunami waves. The concentrations of organic carbon and nitrogen and their isotopic signatures confirm marine origin of these sediments. Content Type Journal Article Pages 1-13 DOI 10.1007/s11069-011-9824-6 Authors Kapila Dahanayake, Department of Geology, University of Peradeniya, Peradeniya, Sri Lanka Nayomi Kulasena, Roughton International (UK), Southern Transport Development Project, Baddegama, Sri Lanka G. V. Ravi Prasad, AMS Radiocarbon Laboratory, Institute of Physics, Bhubaneswar, 751005 India Koushik Dutta, AMS Radiocarbon Laboratory, Institute of Physics, Bhubaneswar, 751005 India D. K. Ray, AMS Radiocarbon Laboratory, Institute of Physics, Bhubaneswar, 751005 India Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    Assessment and inventory on soil erosion hazard are essential for the formulation of successful hazard mitigation plans and sustainable development. The objective of this study was to assess and map soil erosion hazard in Lesser Himalaya with a case study. The Dabka watershed constitutes a part of the Kosi Basin in the Lesser Himalaya, India, in district Nainital has been selected for the case illustration. The average rate of erosion hazard is 0.68 mm/year or 224 tons/km 2 /year. Anthropogenic and geo-environmental factors have together significantly accelerated the rate of erosion. This reconnaissance study estimates the erosion rate over the period of 3 years (2006–2008) as 1.21 mm/year (398 tons/km 2 /year) in the barren land having geological background of diamictite, siltstone and shale rocks, 0.92 mm/year (302 tons/km 2 /year) in the agricultural land with lithology of diamictite, slates, siltstone, limestone rocks, while in the forest land, it varies between 0.20 mm/year (66 tons/km 2 /year) under dense forest land having the geology of quartzwacke and quartrenite rocks and 0.40 mm/year (132 tons/km 2 /year) under open forest/shrubs land having geological setup of shale, dolomite and gypsum rocks. Compared to the intensity of erosion in the least disturbed dense forest, the erosion rate is about 5–6 times higher in the most disturbed agricultural land and barren land, respectively. The erosion hazard zones delineated following scalogram modelling approach. Integrated scalogram modelling approach resulted in severe classes of soil erosion hazard in the study area with numerical values of Erosion Hazard Index (EHI) ranging between 01 (very low hazard) and 5 (very high hazard). Content Type Journal Article Pages 1-21 DOI 10.1007/s11069-011-9833-5 Authors Pradeep K. Rawat, Department of Geography, Kumaun University, Nainital, 263002 India P. C. Tiwari, Department of Geography, Kumaun University, Nainital, 263002 India C. C. Pant, Department of Geology, Kumaun University, Nainital, 263002 India A. K. Sharama, Department of Geology, Kumaun University, Nainital, 263002 India P. D. Pant, Department of Geology, Kumaun University, Nainital, 263002 India Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    Soil erosion is a serious environmental problem in Indravati catchment. It carries the highest amount of sediments compared with other catchments in India. This catchment spreading an area of 41,285 km 2 is drained by river Indravati, which is one of the northern tributaries of the river Godavari in its lower reach. In the present study, USLE is used to estimate potential soil erosion from river Indravati catchment. Both magnitude and spatial distribution of potential soil erosion in the catchment is determined. The derived soil loss map from USLE model is classified into six categories ranging from slight to very severe risk depending on the calculated soil erosion amount. The soil erosion map is linked to elevation and slope maps to identify the area for conservation practice in order to reduce the soil loss. From the model output predictions, it is found that average erosion rate predicted is 18.00 tons/ha/year and sediment yield at the out let of the catchment is 22.30 Million tons per annum. The predicted sediment yield verified with the observed data. Content Type Journal Article Pages 1-17 DOI 10.1007/s11069-011-9832-6 Authors Sreenivasulu Vemu, Department of Civil Engineering, JNT University, Kakinada, 533003 AP, India Udaya Bhaskar Pinnamaneni, Department of Civil Engineering, JNT University, Kakinada, 533003 AP, India Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    This paper introduces a new metric for tropical cyclone track shape within the tropical South Pacific (TSP) basin, based on measurements of track sinuosity. A sinuosity index (SI) is developed by applying a simple cube-root transformation to original track sinuosity values. Based on the resulting near-normal SI distribution, an ordinal four-category (quartile) naming system is then proposed for track-type classification. Track sinuosity patterns are also investigated over the last four decades (1970–2008). Analytical findings suggest that cyclone track sinuosity is an important parameter influencing the potential vulnerability of island archipelagoes to cyclone hazard. Principally, sinuously moving cyclones show some tendency for greater longevity and intensity than straighter-tracking storms and make up a larger proportion of systems forming in the western tropical South Pacific than those generated farther east. Although no long-term statistical trend can be established, track sinuosity is highly variable through time, implying that the TSP basin and the islands therein will continue to experience large but irregular inter-annual fluctuations in cyclone track morphology. Content Type Journal Article Pages 1-14 DOI 10.1007/s11069-011-9827-3 Authors James P. Terry, Department of Geography, National University of Singapore, AS2, 1 Arts Link, Kent Ridge, 117570 Singapore Gennady Gienko, Department of Geomatics, School of Engineering, University of Alaska Anchorage, 3211 Providence Drive, Anchorage, AK 99508-4514, USA Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    The level of damage of flood events does not solely depend on exposure to flood waters. Vulnerabilities due to various socio-economic factors such as population at risk, public awareness, and presence of early warning systems, etc. should also be taken into account. Federal and state agencies, watershed management coalitions, insurance companies, need reliable decision support system to evaluate flood risk, to plan and design flood damage assessment and mitigation systems. In current practice, flood damage evaluations are generally carried out based on results obtained from one dimensional (1D) numerical simulations. In some cases, however, 1D simulation is not able to accurately capture the dynamics of the flood events. The present study describes a decision support system, which is based on 2D flood simulation results obtained with CCHE2D-FLOOD. The 2D computational results are complemented with information from various resources, such as census block layer, detailed survey data, and remote sensing images, to estimate loss of life and direct damages (meso or micro scale) to property under uncertainty. Flood damage calculations consider damages to residential, commercial, and industrial buildings in urban areas, and damages to crops in rural areas. The decision support system takes advantage of fast raster layer operations in a GIS platform to generate flood hazard maps based on various user-defined criteria. Monte Carlo method based on an event tree analysis is introduced to account for uncertainties in various parameters. A case study illustrates the uses of the proposed decision support system. The results show that the proposed decision support system allows stake holders to have a better appreciation of the consequences of the flood. It can also be used for planning, design, and evaluation of future flood mitigation measures. Content Type Journal Article Pages 1-19 DOI 10.1007/s11069-011-9822-8 Authors Honghai Qi, National Center for Computational Hydroscience and Engineering (NCCHE), The University of Mississippi, University, Oxford, MS 38677, USA M. S. Altinakar, National Center for Computational Hydroscience and Engineering (NCCHE), The University of Mississippi, University, Oxford, MS 38677, USA Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    In this paper, we present a probabilistic seismic hazard analysis (PSHA) for mainland Spain that takes into account recent new results in seismicity, seismic zoning, and strong ground attenuation not considered in the latest PSHA of the Spanish Building Code. Those new input data have been obtained as a three-step project carried out in order to improve the existing hazard map for mainland Spain. We have produced a new earthquake catalogue for the area, in which the earthquakes are given in moment magnitude through specific deduced relationships for our territory based on intensity data (Mezcua et al. in Seismol Res Lett 75:75–81, 2004 ). In addition, we included a new seismogenetic zoning based on the recent partial zoning studies performed by different authors. Finally, as we have developed a new strong ground motion model for the area García Blanco ( 2009 ), it was considered in the hazard calculation together with other attenuations gathered from different authors using data compatible with our region. With this new data, a logic tree process is defined to quantify the epistemic uncertainty related to those parts of the process. A sensitivity test has been included in order to analyze the different models of ground motion and seismotectonic zonation used in this work. Finally, after applying a weighting scheme, a mean hazard map for PGA, based on rock type condition for 10% exceedance probability in 50 years, is presented, including 15th and 85th percentile hazard maps. The main differences with the present official building code hazard map are analyzed. Content Type Journal Article Pages 1-22 DOI 10.1007/s11069-011-9819-3 Authors Julio Mezcua, Instituto Geográfico Nacional, General Ibáñez de Ibero 3, 28003 Madrid, Spain Juan Rueda, Instituto Geográfico Nacional, General Ibáñez de Ibero 3, 28003 Madrid, Spain Rosa M. García Blanco, Departamento de Ingeniería Topográfica y Cartografía, Universidad Politécnica de Madrid, Campus Sur, Autovía de Valencia km 7.5, 28031 Madrid, Spain Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    Earthquake-induced landslides are responsible worldwide for significant socioeconomic losses and historically have a prominent position in the list of natural hazards affecting the Iran plateau. As a step toward the development of tools for the assessment and the management of this kind of hazard at regional scale, an empirical estimator of coseismic displacements along potential sliding surfaces was obtained through a regression analysis for the Zagros region, a mountainous Iranian region subjected to earthquake-induced landslides. This estimator, based on the Newmark’s model, allows to evaluate the expected permanent displacement (named “Newmark displacement”) induced by seismic shaking of defined energy on potential sliding surface characterized by a given critical acceleration. To produce regression models for Newmark displacement estimators, a data set was constructed for different critical acceleration values on the basis of 108 accelerometric recordings from 80 Iranian earthquakes with moment magnitudes between 3.6 and 7. The empirical estimator has a general form, proposed by Jibson (Eng Geol 91:209–218, 2007 ), relating Newmark displacement to Arias intensity (as parameter representing the energy of the seismic forces) and to critical acceleration (as parameter representing the dynamic shear resistance of the sliding mass). As an example of application, this relation was employed to provide a basic document for earthquake-induced landslide hazard assessment at regional scale, according to a method proposed by Del Gaudio et al. (Bull Seismol Soc Am 93:557–569, 2003 ), applied to the whole Iranian territory, including Zagros region. This method consists in evaluating the shear resistance required to slopes to limit the occurrence of seismically induced failures, on the basis of the Newmark’s model. The obtained results show that the exposure to landslide seismic induction is maximum in the Alborz Mountains region, where critical accelerations up to ~0.1 g are required to limit the probability of seismic triggering of coherent type landslides within 10% in 50 years. Content Type Journal Article Pages 1-15 DOI 10.1007/s11069-011-9829-1 Authors Ali M. Rajabi, Civil Engineering Department, Qom University, PO Box 14115-175, Qom, Iran Mohammad R. Mahdavifar, Geotechnical Engineering Research Center, International Institute of Earthquake Engineering and Seismology (IIEES), PO Box 19395-3913, Tehran, Iran M. Khamehchiyan, Engineering Geology Department, Tarbiat Modares University, PO Box 14115-175, Tehran, Iran V. Del Gaudio, Dipartimento di Geologia e Geofisica, Universita` degli Studi di Bari, Campus, via E.Orabona 4, 70125 Bari, Italy Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    This paper presents the development of the Coastal Community Social Vulnerability Index (CCSVI) in order to quantify the social vulnerability of hurricane-prone areas under various scenarios of climate change. The 2004–2005 Atlantic hurricane seasons is estimated to have caused 150 billion dollars in damages , and in recent years , the annual hurricane damage in the United States is estimated at around 6 billion. Hurricane intensity or/and frequency may change due to the increase in sea surface temperature as a result of climate change. Climate change is also predicted to cause a rise in sea levels, potentially resulting in higher storm surges. The CCSVI combines the intensity of hurricanes and hurricane-induced surge to create a comprehensive index that considers the effects of a changing climate. The main contributing factors of social vulnerability (such as race, age, gender, and socioeconomic status) in hurricane-prone areas are identified through a principal components analysis. The impact of social characteristics on the potential hurricane damage under various scenarios of climate change are evaluated using Miami-Dade County, Florida, as a case study location. This study finds that climate change may have a significant impact on the CCSVI. Content Type Journal Article Pages 1-21 DOI 10.1007/s11069-011-9817-5 Authors Sigridur Bjarnadottir, Department of Civil and Environmental Engineering, Michigan Technological University, Houghton, MI, USA Yue Li, Department of Civil and Environmental Engineering, Michigan Technological University, Houghton, MI, USA Mark G. Stewart, Center for Infrastructure Performance and Reliability, School of Engineering, The University of Newcastle, Callaghan, NSW, Australia Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    In this paper, the applicability of the Newmark method at regional, sub-regional and site scales has been investigated in the Lorca Basin (Murcia). This basin is located in one of the most seismically active regions of Spain. The area is very interesting for studying earthquake-induced slope instabilities as there are well-known cases associated with specific earthquakes. For the regional and sub-regional scales, a geographic information system has been used to develop an implementation of Newmark sliding rigid block method. Soil and topographic amplification effects have been particularly considered. Subsequently, ‘Newmark displacement’ maps for deterministic seismic scenarios have been produced. Some specific studies have also been performed using limit equilibrium methods to estimate the safety factor and the critical acceleration of certain slope instabilities at a site scale. These instabilities were the rock slides related to recent seismic series at the Lorca Basin: 2002 Bullas ( M w  = 5.0) and 2005 La Paca ( M w  = 4.8). Finally, the safety factor, critical acceleration and Newmark displacement values estimated at different scales have been compared to determine which scale is most suitable for the Newmark method. Content Type Journal Article Pages 1-16 DOI 10.1007/s11069-011-9820-x Authors Martín Jesús Rodríguez-Peces, Departamento de Geodinámica, Facultad de Ciencias, Universidad de Granada, C/Fuentenueva s/n, 18002 Granada, Spain José Luis Pérez-García, Departamento de Ingeniería Cartográfica, Geodésica y Fotogrametría, Universidad de Jaén, Campus de las Lagunillas, s/n. Edif. A3, 23071 Jaén, Spain Julián García-Mayordomo, Instituto Geológico y Minero de España (IGME), C/La Calera, 1, 28760 Tres Cantos, Madrid, Spain José Miguel Azañón, Departamento de Geodinámica, Facultad de Ciencias, Universidad de Granada, C/Fuentenueva s/n, 18002 Granada, Spain Juan Miguel Insua-Arévalo, Departamento de Geodinámica, Facultad de Ciencias Geológicas, Universidad Complutense de Madrid, 28040 Madrid, Spain Jorge Delgado-García, Departamento de Ingeniería Cartográfica, Geodésica y Fotogrametría, Universidad de Jaén, Campus de las Lagunillas, s/n. Edif. A3, 23071 Jaén, Spain Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    The present study explored the effect of assimilation of Advanced TIROS Vertical Sounder (ATOVS) temperature and humidity profiles and Spectral sensor microwave imager (SSM/I) total precipitable water (TPW) on the simulation of a monsoon depression which formed over the Arabian Sea during September 2005 using the Weather Research and Forecast model. The three-dimensional variational (3DVAR) data assimilation technique has been employed for the purpose of assimilation of satellite observations. Statistical scores like “equitable threat score,” “bias score,” “forecast impact,” and “improvement parameter” have been used to examine the impact of the above-mentioned satellite observations on the numerical simulation of a monsoon depression. The diagnostics of this study include verification of the vertical structure of depression, in terms of temperature anomaly profiles and relative vorticity profiles with observations/analysis. Additional diagnostics of the study include the analysis of the heat budget and moisture budget. Such budget studies have been performed to provide information on the role of cumulus convection associated with the depression. The results of this study show direct and good evidence of the impact of the assimilation of the satellite observations using 3DVAR on the dynamical and thermodynamical features of a monsoon depression along with the effect of inclusion of satellite observation on the spatial pattern of the simulated precipitation associated with the depression. The “forecast impact” parameter calculated for the wind speed provides good evidence of the positive impact of the assimilation of ATOVS temperature and humidity profiles and SSM/I TPW on the model simulation, with the assimilation of the ATOVS profiles showing better impact in terms of a more positive value of the “forecast impact” parameter. The results of the study also indicate the improvement of the forecast skill in terms of “equitable threat score” and “bias score” due to the assimilation of satellite observation. Content Type Journal Article Pages 1-23 DOI 10.1007/s11069-011-9857-x Authors M. Govindankutty, Department of Physics and Meteorology, Indian Institute of Technology, Kharagpur, India A. Chandrasekar, Department of Earth and Space Sciences, Indian Institute of Space Science and Technology, Thiruvananthapuram, India Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description: The water retention curve (θ(ψ)), which defines the relationship between soil volumetric water content (θ) and matric potential (ψ), is of paramount importance in characterizing the hydraulic behaviour of soils. However, few methods are so far available for estimating θ(ψ) in undisturbed soil samples. We present a new design of TDR-pressure cell (TDR-Cell) for estimating θ(ψ) in undisturbed soil samples. The TDR-Cell consists of a 50-mm-long and 50-mm internal diameter stainless steel cylinder (which constitutes the outer frame of a coaxial line) attached to a porous ceramic disc and closed at the ends with two aluminium lids. A 49-mm-long and 3-mm-diameter stainless steel rod, which runs longitudinally through the centre of the cylinder, constitutes the inner rod of a coaxial TDR probe. The TDR-Cell was used to determine the θ(ψ) curves of a packed sand and seven undisturbed soil samples from three profiles of agricultural soils. These θ(ψ) curves were subsequently compared to those obtained from the corresponding 2-mm sieved soils using the pressure plate method. Measurements of bulk electrical conductivity, σ a , as a function of the water content, σ a (θ), of the undisturbed soil samples were also performed. An excellent correlation (R 2 = 0.988) was found between the θ values measured by TDR on the different undisturbed soils and the corresponding θ obtained from the soil gravimetric water content. A typical bimodal θ(ψ) function was found for most of the undisturbed soil samples. Comparison between the θ(ψ) curves measured with the TDR-Cell and those obtained from the 2-mm sieved soils showed that the pressure plate method overestimates θ at low ψ values. The σ a (θ) relationship was well described by a simple power expression (R 2 〉 0.95), in which the power factor, defined as tortuosity, ranged between 1.18 and 3.75. Copyright © 2011 John Wiley & Sons, Ltd.

Description: Snowfall is an important part of the yearly water balance for the Catskill Mountains in New York State, the location of water supply reservoirs for New York City. Recent studies have shown that the effects of climate change on the hydrology of the Catskills will most likely create (1) a decrease in the proportion of precipitation falling as snow, (2) a shift in the timing of snowmelt that will cause snowmelt-supplemented streamflow events to occur earlier in the fall and winter, and (3) a decrease in the magnitude of traditionally high April streamflow. The shift in timing of snowmelt-influenced streamflow events is measured by the winter-early spring centre of volume (WSCV), defined as the Julian Day on which half the total streamflow volume from January to May occurs. Studies of streamflow, precipitation, and temperature trends in the last 50 years have shown that the WSCV is already earlier by about 5–10 days. This study investigates the use of watershed-scale snowpack and snowmelt algorithms that are incorporated in two existing watershed water quality models, Generalized Watershed Loading Functions-Variable Source Area (GWLF-VSA) and Soil and Water Assessment Tool (SWAT), to capture the potential effects of climate change on the timing and magnitude of streamflow during the late fall, winter, and early spring for the Catskill Mountain region. The GWLF-VSA model reasonably simulated the recent shifts in the winter streamflow timing, with simulations over the previous 50-year period yielding shifts in WSCV of 2–15 days. The SWAT model yielded similar results as the GWLF-VSA simulations. Scenarios of potential climate change 100 years in the future showed a similar shift in direction of timing winter streamflow, but at a larger magnitude than observed to date with WSCV occurring 15–20 days earlier. Copyright © 2011 John Wiley & Sons, Ltd.

Description:    Highways and railroads situated within rugged terrain are often subjected to the hazard of rockfalls. The task of assessing roadside rockmasses for potential hazards typically involves an on-site visual investigation of the rockmass by an engineer or geologist. At that time, numerous parameters associated with discontinuity orientations and spacing, block size (volume) and shape distributions, slope geometry, and ditch profile are either measured or estimated. Measurements are typically tallied according to a formal hazard rating system, and a hazard level is determined for the site. This methodology often involves direct exposure of the evaluating engineer to the hazard and can also create a potentially non-unique record of the assessed slope based on the skill, knowledge and background of the evaluating engineer. Light Detection and Ranging (LiDAR)–based technologies have the capability to produce spatially accurate, high-resolution digital models of physical objects, known as point clouds. Mobile terrestrial LiDAR equipment can collect, at traffic speed, roadside data along highways and rail lines, scanning continual distances of hundreds of kilometres per day. Through the use of mobile terrestrial LiDAR, in conjunction with airborne and static systems for problem areas, rockfall hazard analysis workflows can be modified and optimized to produce minimally biased, repeatable results. Traditional rockfall hazard analysis inputs include two distinct, but related sets of variables related to geological or geometric control. Geologically controlled inputs to hazard rating systems include kinematic stability (joint identification/orientation) and rock block shape and size distributions. Geometrically controlled inputs include outcrop shape and size, road, ditch and outcrop profile, road curvature and vehicle line of sight. Inputs from both categories can be extracted or calculated from LiDAR data, although there are some limitations and special sampling and processing considerations related to structural character of the rockmass, as detailed in this paper. Content Type Journal Article Pages 1-34 DOI 10.1007/s11069-011-9872-y Authors Matthew J. Lato, Norwegian Geotechnical Institute, Oslo, Norway Mark S. Diederichs, Geological Sciences and Geological Engineering, Queen’s University, Kingston, ON, Canada D. Jean Hutchinson, Geological Sciences and Geological Engineering, Queen’s University, Kingston, ON, Canada Rob Harrap, Geological Sciences and Geological Engineering, Queen’s University, Kingston, ON, Canada Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    A rock slide on top of the clayey–silty–sandy–pebbly masses in the Gschliefgraben (Upper Austria province, Lake Traunsee) having occurred in 2006, together with the humid autumn of 2007, triggered a mudslide comprising a volume up to 4 million m³ and moving with a maximum displacement velocity of 5 m/day during the winter of 2007–2008. The possible damage was estimated up to 60 million € due to the possible destruction of houses and of a road to a settlement with intense tourism. The movement front ran ahead in the creek bed. Therefore, it was assumed that water played an important role. Inclinometer measurements showed that a less permeable layer was sliding on a thin, more permeable layer. During the last centuries, mudslides had already pushed farms into the lake, as reported by chronicles. Thus, the inhabitants of 46 houses had to be evacuated for safety reasons. They could return to their homes after displacement velocities had decreased. It was decided (a) to prevent soaking of water into the uppermost, less permeable layer by transversal drainages, (b) to lower the pore water pressures by longitudinal trenches filled with blocky material, (c) to pump water out of the more permeable layer by well drillings upslope of the houses in order to create a stable block below the houses and (d) to remove material thrust over the stable blocks in order to avoid damage to the houses. These mitigation measures costing 11.5 million € led to a deceleration of the process to displacement velocities of some cm/year up to now. The houses and the road were not damaged. Content Type Journal Article Pages 1-13 DOI 10.1007/s11069-011-9865-x Authors R. Poisel, Vienna University of Technology, Vienna, Austria R. Hofmann, Geotechnical Engineering Hofmann, Perchtoldsdorf, Austria A. Preh, Vienna University of Technology, Vienna, Austria T. Sausgruber, Austrian Service for Torrent, Erosion and Avalanche Control, Geological Service, Innsbruck, Austria M. Schiffer, Austrian Service for Torrent, Erosion and Avalanche Control, Project Management Salzkammergut, Bad Ischl, Austria Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    The rainfall events that occurred in the Cubatão do Norte River watershed, Santa Catarina State, Brazil, in 2008, were characterized by both high intensities and amounts and triggered landslides in this watershed. The objective of the present study was to analyze the influence of landslides on the turbidity and the total solid concentration (TS) in this river using turbidity, TS, and river discharge data obtained from March 23, 2008, to June 11, 2010. The comparison between turbidity and discharge patterns implies that the landslide process was not continuous and increased the turbidity intermittently and irregularly. The sediment yield during the landslide occurrence was approximately five times higher than without the landslide, even though the discharges were similar. The turbidity/discharge relationship during the landslide occurrence was markedly different from that before and after the occurrence. The correlation coefficients between turbidity and TS showed that the landslide significantly changed the sediment yield in this watershed. The result indicates that sediment yield estimations at the watershed level should be treated more carefully when landslides occur. Content Type Journal Article Pages 1-10 DOI 10.1007/s11069-011-9818-4 Authors Masato Kobiyama, Departamento de Engenharia Sanitária e Ambiental, Universidade Federal de Santa Catarina, Caixa Postal 476, Florianopolis, SC 88040-900, Brazil Aline de Almeida Mota, Departamento de Engenharia Sanitária e Ambiental, Universidade Federal de Santa Catarina, Caixa Postal 476, Florianopolis, SC 88040-900, Brazil Fernando Grison, Departamento de Engenharia Sanitária e Ambiental, Universidade Federal de Santa Catarina, Caixa Postal 476, Florianopolis, SC 88040-900, Brazil Joana Nery Giglio, Departamento de Engenharia Sanitária e Ambiental, Universidade Federal de Santa Catarina, Caixa Postal 476, Florianopolis, SC 88040-900, Brazil Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    This study describes a critical assessment of the risk control decision model from a methodological perspective and identifies major shortcomings with the employment of enhanced formal evaluation and decision-making methods. This in turn could have major applications for natural disaster risk control. The methodology follows the description of interpretive structural modeling (ISM), which is an interactive learning process in which a set of different and directly related elements is structured to form a comprehensive systemic model. The next step explores the potentials of different mathematical programming approaches in order to improve decision making, i.e., for the development of an economic utility constrained-maximization model that addresses the issue of optimal budget allocation under a trade-off framework. Several aspects of risk and uncertainty are discussed within the context of an economic utility constrained-maximization model with a major focus on the importance of risk and uncertainty in research evaluation, and how the strategy determines the insurance and risk control plans. Content Type Journal Article Pages 1-19 DOI 10.1007/s11069-011-9861-1 Authors Chun-Pin Tseng, Armaments Bureau, Chung Shan Institute of Science and Technology, Taoyuan, Taiwan Cheng-Wu Chen, Institute of Maritime Information and Technology, National Kaohsiung Marine University, Kaohsiung, 80543 Taiwan Yu-Ping Tu, Department of Logistics Management, Shu-Te University, Kaohsiung, 82445 Taiwan, Republic of China Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    District-wide drought climatology over India for the southwest monsoon season (June–September) has been examined using two simple drought indices; Percent of Normal Precipitation (PNP) and Standardized Precipitation Index (SPI). The season drought indices were computed using long times series (1901–2003) of southwest monsoon season rainfall data of 458 districts over the country. Identification of all India (nation-wide) drought incidences using both PNP and SPI yielded nearly similar results. However, the district-wide climatology based on PNP was biased by the aridity of the region. Whereas district-wide drought climatology based on SPI was not biased by aridity. This study shows that SPI is a better drought index than PNP for the district-wide drought monitoring over the country. SPI is also suitable for examining break and active events in the southwest monsoon rainfall over the country. The trend analysis of district-wide season (June–September) SPI series showed significant negative trends over several districts from Chattisgarh, Bihar, Kerala, Jharkhand, Assam and Meghalaya, Uttaranchal, east Madhya Pradesh, Vidarbha etc., Whereas significant positive trends in the SPI series were observed over several districts from west Uttar Pradesh, west Madhya Pradesh, South & north Interior Karnataka, Konkan and Goa, Madhya Maharashtra, Tamil Nadu, East Uttar Pradesh, Punjab, Gujarat etc. Content Type Journal Article Pages 1-17 DOI 10.1007/s11069-011-9867-8 Authors D. S. Pai, India Meteorological Department, Pune, India Latha Sridhar, India Meteorological Department, Pune, India Pulak Guhathakurta, India Meteorological Department, Pune, India H. R. Hatwar, India Meteorological Department, Pune, India Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    The Indus flood in 2010 was one of the greatest river disasters in recent history, which affected more than 14 million people in Pakistan. Although excessive rainfall between July and September 2010 has been cited as the major causative factor for this disaster, the human interventions in the river system over the years made this disaster a catastrophe. Geomorphic analysis suggests that the Indus River has had a very dynamic regime in the past. However, the river has now been constrained by embankments on both sides, and several barrages have been constructed along the river. As a result, the river has been aggrading rapidly during the last few decades due to its exceptionally high sediment load particularly in reaches upstream of the barrages. This in turn has caused significant increase in cross-valley gradient leading to breaches upstream of the barrages and inundation of large areas. Our flow accumulation analysis using SRTM data not only supports this interpretation but also points out that there are several reaches along the Indus River, which are still vulnerable to such breaches and flooding. Even though the Indus flood in 2010 was characterized by exceptionally high discharges, our experience in working on Himalayan rivers and similar recent events in rivers in Nepal and India suggest that such events can occur at relatively low discharges. It is therefore of utmost importance to identify such areas and plan mitigation measures as soon as possible. We emphasize the role of geomorphology in flood analysis and management and urge the river managers to take urgent steps to incorporate the geomorphic understanding of Himalayan rivers in river management plans. Content Type Journal Article Pages 1-12 DOI 10.1007/s11069-011-9869-6 Authors Kumar Gaurav, Department of Civil Engineering, Indian Institute of Technology Kanpur, Kanpur, UP 208016, India R. Sinha, Department of Civil Engineering, Indian Institute of Technology Kanpur, Kanpur, UP 208016, India P. K. Panda, Department of Civil Engineering, Indian Institute of Technology Kanpur, Kanpur, UP 208016, India Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    Tropical cyclone is one of the most devastating weather phenomena all over the world. The Environmental Modeling Center (EMC) of the National Center for Environmental Prediction (NCEP) has developed a sophisticated mesoscale model known as Hurricane Weather Research and Forecasting (HWRF) system for tropical cyclone studies. The state-of-the-art HWRF model (atmospheric component) has been used in simulating most of the features our present study of a very severe tropical cyclone “Mala”, which developed on April 26 over the Bay of Bengal and crossed the Arakan coast of Myanmar on April 29, 2006. The initial and lateral boundary conditions are obtained from Global Forecast System (GFS) analysis and forecast fields of the NCEP, respectively. The performance of the model is evaluated with simulation of cyclone Mala with six different initial conditions at an interval of 12 h each from 00 UTC 25 April 2006 to 12 UTC 27 April 2006. The best result in terms of track and intensity forecast as obtained from different initial conditions is further investigated for large-scale fields and structure of the cyclone. For this purpose, a number of important predicted fields’ viz. central pressure/pressure drop, winds, precipitation, etc. are verified against observations/verification analysis. Also, some of the simulated diagnostic fields such as relative vorticity, pressure vertical velocity, heat fluxes, precipitation rate, and moisture convergences are investigated for understanding of the characteristics of the cyclone in more detail. The vector displacement errors in track forecasts are calculated with the estimated best track provided by the India Meteorological Department (IMD). The results indicate that the model is able to capture most of the features of cyclone Mala with reasonable accuracy. Content Type Journal Article Pages 1-25 DOI 10.1007/s11069-011-9863-z Authors Sujata Pattanayak, Centre for Atmospheric Sciences, Indian Institute of Technology, Delhi, Hauz Khas, New Delhi, 110016 India U. C. Mohanty, Centre for Atmospheric Sciences, Indian Institute of Technology, Delhi, Hauz Khas, New Delhi, 110016 India S. G. Gopalakrishnan, Hurricane Research Division, NOAA, Atlantic Oceanographic and Meteorological Laboratory, 4301 Rickenbacker Causeway, Miami, FL 33149, USA Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description: In semi-arid areas high-intensity rainfall events are often held responsible for the main part of soil erosion. Long-term landscape evolution models usually use average annual rainfall as input, making the evaluation of single events impossible. Event-based soil erosion models are better suited for this purpose, but cannot be used to simulate longer timescales and are usually applied to plots or small catchments. In this study, the openLISEM event-based erosion model was applied to the medium sized ( 50 km 2 ) Prado catchment in SE Spain. Our aim was to (i) test the model's performance for medium sized catchments; (ii) test the ability to simulate four selected typical Mediterranean rainfall events of different magnitude, and (iii) explore the relative contribution of these different storms to soil erosion using scenarios of future climate variability. Results show that due to large differences in the hydrologic response between storms of different magnitudes, each event needed to be calibrated separately. The relation between rainfall event characteristics and the calibration factors might help in determining optimal calibration values if event characteristics are known. Calibration of the model features some drawbacks for large catchments due to spatial variability in K sat values. Scenario calculations show that, although 50% of soil erosion occurs as a result of high frequency, low intensity rainfall events, large magnitude, low frequency events potentially contribute significantly to total soil erosion. The results illustrate the need to incorporate temporal variability in rainfall magnitude-frequency distributions in landscape evolution models. Copyright © 2011 John Wiley & Sons, Ltd.

Description: Sustainable water resources management require scientifically sound information on precipitation, as it plays a key role in hydrological responses in a catchment. In recent years, mesoscale weather models in conjunction with hydrological models have gained great attention as they can provide high resolution downscaled weather variables. Many cumulus parameterization schemes (CPSs) have been developed and incorporated into three-dimensional Pennsylvania State University/National Center for Atmospheric Research (PSU/NCAR) mesoscale model 5 (MM5). This study has performed a comprehensive evaluation of four CPSs (the Anthes–Kuo, Grell, Betts–Miller and Kain–Fritsch93 schemes) to identify how their inclusion influences the mesoscale model's precipitation estimation capabilities. The study has also compared these four CPSs in terms of variability in rainfall estimation at various horizontal and vertical levels. For this purpose, the MM5 was nested down to resolution of 81km for domain 1 (domain span 21 x 81 km) and 3km for 4 (domain span 16 x 3km) respectively with vertical resolutions at 23, 40 and 53 vertical levels. The study was carried out at the Brue catchment in Southwest England using both the ERA-40 reanalysis data and the land based observation data. The performances of four CPs were evaluated in terms of their ability to simulate the amount of cumulative rainfall in four months in 1995 representing the four seasonal months viz January (winter), March (spring), July (summer) and October (autumn). It is observed that the Anthes–Kuo scheme has produced inferior precipitation values during spring and autumn seasons while simulations during winter and summer were consistently good. The Betts–Miller scheme has produced some reasonable results, particularly at the small scaled domain (3 km grid size) during winter and summer. The KF2 scheme was the best scheme for the larger scale (81 km grid size) domain during winter season at both 23 and 53 vertical levels. This scheme tended to underestimate rainfall for other seasons including the small scale domain (3 km grid size) in the mesoscale. The Grell scheme was the best scheme in simulating rainfall rates, and was found to be superior to other three schemes with consistently better results in all four seasons and different domain scales. Copyright © 2011 John Wiley & Sons, Ltd.

Description: Nonlinear dynamics and spatial variability in hydrological systems make difficult the formulation of scaling theories. Therefore, the development of knowledge related to scale effects, scaling techniques, parameterization and linkages of parameters across scales is highly relevant. The main purpose of this work is to analyse the spatial effect of the static storage capacity parameter H u and the saturated hydraulic conductivity parameter k s from microscale (sub-grid level) to mesoscale (grid level) and its implication to the definition of an optimum cell size. These two parameters describe the upper soil water characteristics in the infiltration process conceptualization of the TETIS hydrological model. At microscale, the spatial heterogeneity of H u and k s was obtained generating random parameter fields through probability distribution functions and a spatial dependence model with pre-established correlation lengths. The effective parameters at mesoscale were calculated by solving the inverse problem for each parameter field. Results indicate that the adopted inverse formulation allows transferring the non-linearity of the system from microscale to the mesoscale via non-stationary effective parameters. Their values at each cell and time step are in the range of zero to the mean value of the parameter at microscale. The stochastic simulations showed that the variance of the estimated effective parameters decreases when the ratio between mesoscale cell size and correlation length at microscale increases. For a ratio greater than 1, we found cell sizes having the characteristics of a representative elementary area (REA); in such case, the microscale variability pattern did not affect the system response at mesoscale. Copyright © 2011 John Wiley & Sons, Ltd.

Description:    Seasonal water scarcity in southern China has been an issue of concern for many years. The increased frequency of low precipitation in the growing season of rice created a flurry of discussions in the academic and policy arenas. These events severely disrupted the supply of irrigation water for agriculture in paddy field areas and posed a substantial threat to farmers’ livelihoods. Within a broader context of accessing farmers’ resilience to agricultural drought, this paper focuses on the response mechanisms and adaptive strategies adopted by farming households in three types of areas (Plain, Hill, Mountain) in Dingcheng, Hunan Province. With the increasing drought frequency and the pressure from the demand for livelihood improvement, farmers’ response mechanisms have evolved, expanding from short-term adjustments to long-term adaptations, and switching focus from securing reliable water sources to improving irrigation efficiency and diversifying both on- and off- farm productions. The three types of geographic units have different resilience profiles and have developed diverse patterns of adaptive processes that update the conceptual model of Disaster Resilience of “Loss-Response” of Location. It presents a temporal dimension to the study of resilience, which is largely missing from the current literature and provides insights into how to enhance farmers’ response capacities in the face of agricultural drought in southern China. Content Type Journal Article Pages 1-13 DOI 10.1007/s11069-011-9873-x Authors Yehong Sun, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing, 100101 China Hongjian Zhou, National Disaster Reduction Center of China (NDRCC), No. 6, Guangbaidonglu, Chaoyang District, Beijing, 100022 China Jing’ai Wang, College of Geography and Remote Sensing Science, Beijing Normal University, Beijing, 100875 China Yi Yuan, National Disaster Reduction Center of China (NDRCC), No. 6, Guangbaidonglu, Chaoyang District, Beijing, 100022 China Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    Using newly digitised sea-level data for the ports of Southampton (1935–2005) and Portsmouth (1961–2005) on the south coast of the UK, this study investigates the relationship between the 100 highest sea-level events recorded at the two cities and the incidence of coastal floods in the adjoining Solent region. The main sources of flood data are the daily newspapers The Southern Daily Echo , based in Southampton and The News , based in Portsmouth, supported by a range of local publications and records. The study indicates a strong relationship between the highest measured sea levels and the incidence of coastal floods and highlights the most vulnerable areas to coastal flooding which include parts of Portsmouth, Southampton, Hayling Island, Fareham and Cowes. The most severe flood in the dataset resulted from the storm surge events of 13–17 December 1989 when eight consecutive extreme high waters occurred. The data suggest that while extreme sea-level events are becoming more common, the occurrence of flood events is not increasing. This is attributed to improved flood remediation measures combined with a reduction of storm intensity since the 1980s. However, several recent events of significance were still recorded, particularly 3 November 2005 when Eaststoke on Hayling Island (near Portsmouth) was flooded due to high sea levels combined with energetic swell waves. Content Type Journal Article Pages 1-24 DOI 10.1007/s11069-011-9868-7 Authors Amy C. Ruocco, School of Civil Engineering and the Environment, Tyndall Centre for Climate Change Research, University of Southampton, Highfield Campus, Southampton, SO17 1BJ UK Robert J. Nicholls, School of Civil Engineering and the Environment, Tyndall Centre for Climate Change Research, University of Southampton, Highfield Campus, Southampton, SO17 1BJ UK Ivan D. Haigh, The School of Environment Systems Engineering and UWA Oceans Institute, University of Western Australia, M470, 35 Stirling Highway, Crawley, WA 6009, Australia Matthew P. Wadey, School of Civil Engineering and the Environment, Tyndall Centre for Climate Change Research, University of Southampton, Highfield Campus, Southampton, SO17 1BJ UK Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description: Field survey report on the 11th March 2011 tsunami in Pacific coast of Mexico Content Type Journal Article Pages 1-6 DOI 10.1007/s11069-011-9871-z Authors M. P. Jonathan, Centro Interdisciplinario de Investigaciones y Estudios sobre Medio Ambiente y Desarrollo (CIIEMAD), Instituto Politécnico Nacional (IPN), Calle 30 de Junio de 1520, Barrio la Laguna Ticomán, Del. Gustavo A. Madero, 07340 Mexico, DF, Mexico P. D. Roy, Instituto de Geología, Universidad Nacional Autónoma de México (UNAM), Ciudad Universitaria, Coyoacán, 04510 Mexico, DF, Mexico J. L. Sánchez-Zavala, Instituto de Geología, Universidad Nacional Autónoma de México (UNAM), Ciudad Universitaria, Coyoacán, 04510 Mexico, DF, Mexico S. Srinivasalu, Department of Geology, Anna University, Chennai, 600025 India M. C. Macías-Romo, Instituto de Geología, Universidad Nacional Autónoma de México (UNAM), Ciudad Universitaria, Coyoacán, 04510 Mexico, DF, Mexico C. Lakshumanan, Centre for Remote Sensing, Bharathidasan University, Khajamalai Campus, Thiruchirappalli, 620023 India Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description: :   In North America and Europe, most fatalities due to snow avalanches occur in the backcountry during recreational pursuits. Of these, more than 90% of the fatal avalanches are triggered by the victims themselves. This pattern suggests that the primary cause of avalanche fatalities for human-triggered avalanches is a failure in human perception. For the latter, people thought that the state of stability or instability of the snow cover was different than it actually was. In this paper, the strength and weight of evidence used to make decisions in backcountry travel are discussed from: (a) the perspective of the favored hypothesis to proceed for good recreational enjoyment based on stability evaluation and (b) the null hypothesis based on an assessment of instability. Based on the facts about snow slab avalanche release, it is argued that instability analysis is the best framework for avalanche forecasting, whereas human action is most closely related to the favored hypothesis (stability evaluation). Using scaling laws derived from: (a) fracture mechanics about the size of imperfections causing avalanches and (b) avalanche dimensions, it is suggested that a snow slab could show stability over more than 99% of the total area. From the concepts of Bayesian probability, it is shown that overconfidence about stability can arise when the weight of the likelihood is high and the weight of prior is low. Similarly, underconfidence (excessive conservatism) often results when the weight of the prior is high with little regard for the likelihood, which may be low. Overconfidence about stability is considered to be a prime source of accidents. Content Type Journal Article Pages 1-11 DOI 10.1007/s11069-011-9856-y Authors D. M. McClung, Department of Geography, University of British Columbia, Vancouver, BC V6T 1Z2, Canada Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    A susceptibility map for an area, which is representative in terms of both geologic setting and slope instability phenomena of large sectors of the Sicilian Apennines, was produced using slope units and a multiparametric univariate model. The study area, extending for approximately 90 km 2 , was partitioned into 774 slope units, whose expected landslide occurrence was estimated by averaging seven susceptibility values, determined for the selected controlling factors: lithology, mean slope gradient, stream power index at the foot, mean topographic wetness index and profile curvature, slope unit length, and altitude range. Each of the recognized 490 landslides was represented by its centroid point. On the basis of conditional analysis, the susceptibility function here adopted is the density of landslides, computed for each class. Univariate susceptibility models were prepared for each of the controlling factors, and their predictive performance was estimated by prediction rate curves and effectiveness ratio applied to the susceptibility classes. This procedure allowed us to discriminate between effective and non-effective factors, so that only the former was subsequently combined in a multiparametric model, which was used to produce the final susceptibility map. The validation of this map latter enabled us to verify the reliability and predictive performance of the model. Slope unit altitude range and length, lithology and, subordinately, stream power index at the foot of the slope unit demonstrated to be the main controlling factors of landslides, while mean slope gradient, profile curvature, and topographic wetness index gave unsatisfactory results. Content Type Journal Article Pages 1-11 DOI 10.1007/s11069-011-9846-0 Authors E. Rotigliano, Dipartimento di Scienze della Terra e del Mare, Università degli Studi di Palermo, Palermo, Italy C. Cappadonia, Dipartimento di Scienze della Terra e del Mare, Università degli Studi di Palermo, Palermo, Italy C. Conoscenti, Dipartimento di Scienze della Terra e del Mare, Università degli Studi di Palermo, Palermo, Italy D. Costanzo, Dipartimento di Scienze della Terra e del Mare, Università degli Studi di Palermo, Palermo, Italy V. Agnesi, Dipartimento di Scienze della Terra e del Mare, Università degli Studi di Palermo, Palermo, Italy Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    The purpose of this study is to evaluate and compare the results of applying the statistical index and the logistic regression methods for estimating landslide susceptibility in the Hoa Binh province of Vietnam. In order to do this, first, a landslide inventory map was constructed mainly based on investigated landslide locations from three projects conducted over the last 10 years. In addition, some recent landslide locations were identified from SPOT satellite images, fieldwork, and literature. Secondly, ten influencing factors for landslide occurrence were utilized. The slope gradient map, the slope curvature map, and the slope aspect map were derived from a digital elevation model (DEM) with resolution 20 × 20 m. The DEM was generated from topographic maps at a scale of 1:25,000. The lithology map and the distance to faults map were extracted from Geological and Mineral Resources maps. The soil type and the land use maps were extracted from National Pedology maps and National Land Use Status maps, respectively. Distance to rivers and distance to roads were computed based on river and road networks from topographic maps. In addition, a rainfall map was included in the models. Actual landslide locations were used to verify and to compare the results of landslide susceptibility maps. The accuracy of the results was evaluated by ROC analysis. The area under the curve (AUC) for the statistical index model was 0.946 and for the logistic regression model, 0.950, indicating an almost equal predicting capacity. Content Type Journal Article Pages 1-32 DOI 10.1007/s11069-011-9844-2 Authors Dieu Tien Bui, Department of Mathematical Sciences and Technology, Norwegian University of Life Sciences, P.O. Box 5003IMT, 1432 Ås, Norway Owe Lofman, Department of Mathematical Sciences and Technology, Norwegian University of Life Sciences, P.O. Box 5003IMT, 1432 Ås, Norway Inge Revhaug, Department of Mathematical Sciences and Technology, Norwegian University of Life Sciences, P.O. Box 5003IMT, 1432 Ås, Norway Oystein Dick, Department of Mathematical Sciences and Technology, Norwegian University of Life Sciences, P.O. Box 5003IMT, 1432 Ås, Norway Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    Changing characteristics of hazardous weather-related events have been arousing considerable public interest in recent years. Guangdong is an economically developed province in China and is prone to natural hazards. Using monthly precipitation data covering a period of 1956–2008 from 127 rain gauge stations, the probabilistic behaviors of SPI-based droughts were investigated with copulas functions. Results indicated a higher risk of droughts along the coastal regions and the western Guangdong, particularly the Pearl River Delta (PRD) region. Joint probabilities of droughts with higher intensity and longer duration were found to have relatively even geographical distribution across Guangdong. The northern parts of Guangdong are higher in altitude and have a lower risk of droughts. Identification of regions characterized by droughts of different severity and durations across Guangdong is important for scientific management of water resource and agricultural activities and also the development of social resilience under the influence of climate changes. Content Type Journal Article Pages 1-14 DOI 10.1007/s11069-011-9850-4 Authors Chun-Ling Liu, Pearl River Hydraulic Research Institute, Guangzhou, 510611 China Qiang Zhang, Department of Water Resources and Environment, Sun Yat-Sen University, Guangzhou, 510275 China Vijay P. Singh, Department of Biological and Agricultural Engineering and Department of Civil and Environmental Engineering, Texas A & M University, College Station, TX 77843-2117, USA Ying Cui, Department of Water Resources and Environment, Sun Yat-Sen University, Guangzhou, 510275 China Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    The landslide studies can be categorized as pre- and postdisaster studies. The predisaster studies include spatial prediction of potential landslide zones known as landslide susceptibility zonation (LSZ) mapping to identify the areas/locales susceptible to landslide hazard. The LSZ maps provide an assessment of the safety of existing habitations and infrastructural/functional elements and help plan further developmental activities in the hilly regions. Landslides are one of the natural geohazards that affect at least 15% of land area of India. Different types of landslides occur frequently in geodynamical active domains of the Himalayas. In India, various techniques have been developed and adopted for LSZ mapping of different regions. However, the technique for LSZ mapping is not yet standardized. The present research is an attempt in this direction only. In our earlier work (Kanungo et al. 2006 ), a detailed study on conventional, artificial neural network (ANN)- black box-, fuzzy set-based and combined neural and fuzzy weighting techniques for LSZ mapping in Darjeeling Himalayas has been documented. In this paper, other techniques such as combined neural and certainty factor concept along with combined neural and likelihood ratio techniques have been assessed in comparison with combined neural and fuzzy technique for the preparation of LSZ maps of the same study area in parts of Darjeeling Himalayas. It is observed from the present study that the LSZ map produced using combined neural and fuzzy approach appears to be the most accurate one as in this case only 2.3% of the total area is found to be categorized as very high susceptibility zone and contains 30.1% of the existing landslide area. This approach can serve as one of the key objective approaches for spatial prediction of landslide hazards in hilly terrain. Content Type Journal Article Pages 1-22 DOI 10.1007/s11069-011-9847-z Authors D. P. Kanungo, Geotechnical Engineering Division, CSIR-Central Building Research Institute, Roorkee, India S. Sarkar, Geotechnical Engineering Division, CSIR-Central Building Research Institute, Roorkee, India Shaifaly Sharma, Geotechnical Engineering Division, CSIR-Central Building Research Institute, Roorkee, India Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    Movement of seasonal eddies in the Bay of Bengal (BOB) and its relation with cyclonic heat potential (CHP) and cyclogenesis points have been investigated in this study using 6 years (2002–2007) of global ocean monthly analysis datasets based on the Simple Ocean Data Assimilation (SODA) package (SODA v2.0.4) of Carton et al. ( 2005 ) and Indian Meteorological Department cyclogenesis points. The region dominated by anticyclonic eddies with CHP greater than 70 × 10 7  J/m 2 as well as good correlations (〉0.9) with sea surface height (SSH) and 26°C isothermal depth ( D 26 ) can be a potential region of cyclogenesis. The region dominated by cyclonic eddies with CHP greater than 50 × 10 7  J/m 2 and good correlation (〉0.9) with both SSH and D 26 can serve as a potential region of high-level depression. Potential cyclogenesis regions are the southern BOB (5°N–12°N) for the post-monsoon season and the head of BOB (north of 15°N) during southwest monsoon. Seven potential regions are identified for the eddy formation for different seasons, which are consistent with the cyclogenesis points. The CHP distributions alone are able to explain the cyclone tracks for the pre-monsoon and post-monsoon seasons but not for the monsoon season. Content Type Journal Article Pages 1-19 DOI 10.1007/s11069-011-9858-9 Authors Bishnu Kumar, Centre for Oceans, Rivers, Atmosphere and Land Sciences (CORAL), Indian Institute of Technology Kharagpur, Kharagpur, India Arun Chakraborty, Centre for Oceans, Rivers, Atmosphere and Land Sciences (CORAL), Indian Institute of Technology Kharagpur, Kharagpur, India Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    To overcome the disadvantages of traditional flow analysis methods for liquefied soils that exhibit fluidization and large deformation characteristics, Smoothed particle hydrodynamics (SPH) is adopted in this study to analyze the flow processes of liquefied soils. Bingham model with the use of the Mohr–Coulomb yield criterion, the concepts of equivalent Newtonian viscosity, and the Verlet neighbor list method are introduced into the framework of SPH to build an algorithm for the analysis of flowing liquefied soils. This modeling involves a simulation of physical model test of flowing liquefied soils that can be compared with numerical results. In addition, a shaking table test is selected from the literature for SPH analysis to verify the validation of the SPH method and extend its applications. The SPH simulation can reproduce the flow processes of liquefied soils and constrain estimates of the horizontal displacement, vertical displacement, and velocity of soils after liquefaction. According to the dynamic behaviors of the materials involved, designs can be implemented to improve the seismic safety of structures. Content Type Journal Article Pages 1-14 DOI 10.1007/s11069-011-9851-3 Authors Yu Huang, Key Laboratory of Geotechnical and Underground Engineering of the Ministry of Education, Tongji University, Shanghai, 200092 China Weijie Zhang, Department of Geotechnical Engineering, Tongji University, Shanghai, 200092 China Wuwei Mao, Department of Geotechnical Engineering, Tongji University, Shanghai, 200092 China Chen Jin, Department of Geotechnical Engineering, Tongji University, Shanghai, 200092 China Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    The 27 November 1945 earthquake in the Makran Subduction Zone triggered a destructive tsunami that has left important problems unresolved. According to the available reports, high waves persisted along the Makran coast and at Karachi for several hours after the arrival of the first wave. Long-duration sea-level oscillations were also reported from Port Victoria, Seychelles. On the other hand, only one high wave was reported from Mumbai. Tide-gauge records of the tsunami from Karachi and Mumbai confirm these reports. While the data from Mumbai shows a single high wave, Karachi data shows that high waves persisted for more than 7 h, with maximum wave height occurring 2.8 h after the arrival of the first wave. In this paper, we analyze the cause of these persistent high waves using a numerical model. The simulation reproduces the observed features reasonably well, particularly the persistent high waves at Karachi and the single high wave at Mumbai. It further reveals that the persistent high waves along the Makran coast and at Karachi were the result of trapping of the tsunami-wave energy on the continental shelf off the Makran coast and that these coastally-trapped edge waves were trapped in the along-shore direction within a ∼300-km stretch of the continental shelf. Sensitivity experiments establish that this along-shore trapping of the tsunami energy is due to variations in the shelf width. In addition, the model simulation indicates that the reported long duration of sea-level oscillations at Port Victoria were mainly due to trapping of the tsunami energy over the large shallow region surrounding the Seychelles archipelago. Content Type Journal Article Pages 1-10 DOI 10.1007/s11069-011-9854-0 Authors S. Neetu, National Institute of Oceanography, Dona Paula, Goa, 403004 India I. Suresh, National Institute of Oceanography, Dona Paula, Goa, 403004 India R. Shankar, The Institute of Mathematical Sciences, Chennai, 600113 India B. Nagarajan, Indian Institute of Surveying and Mapping, Survey of India, Hyderabad, 500039 India R. Sharma, Geodetic and Research Branch, Survey of India, Dehra Dun, 248001 India S. S. C. Shenoi, Indian National Centre for Ocean Information Services, Hyderabad, 500055 India A. S. Unnikrishnan, National Institute of Oceanography, Dona Paula, Goa, 403004 India D. Sundar, National Institute of Oceanography, Dona Paula, Goa, 403004 India Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    A very severe cyclonic storm “Aila” hit West Bengal on 26 May 2009. The storm intensified when it encountered with a warm core (SST = 31°C) anti-cyclonic eddy (ACE4) in the north Bay of Bengal. The storm intensity increased by 43% due to this eddy, which is comparable with that (34%) obtained from a best fit line (derived from several numerical experiments over north-west Pacific Ocean). The shallow mixed layer of the large-scale ocean and deep mixed layer inside the eddy appear to be crucial parameters besides translation speed of the storm (Uh), ambient relative humidity and thermal stratification below mixed layer, in the storm intensification. From the eddy size and Uh, the eddy feedback factor is found to be about 0.4 (i.e. 40%), which is close to the above. Since there exists an inverse relationship between Uh and UOHC (upper ocean heat content), slow (fast) moving storms require high (low) UOHC. The warm ACE4 with a high UOHC of 149 kj/cm 2 (300% higher than the climatological value) and deep warm layer (D26 = 126 m) opposes the cooling induced by the storm and helps for the intensification of the storm through the supply of large enthalpy (latent + sensible) flux. Content Type Journal Article Pages 1-11 DOI 10.1007/s11069-011-9837-1 Authors Y. Sadhuram, National Institute of Oceanography, Regional Centre (CSIR), 176, Lawsons Bay colony, Visakhapatnam, 530017 India K. Maneesha, National Institute of Oceanography, Regional Centre (CSIR), 176, Lawsons Bay colony, Visakhapatnam, 530017 India T. V. Ramana Murty, National Institute of Oceanography, Regional Centre (CSIR), 176, Lawsons Bay colony, Visakhapatnam, 530017 India Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    Flooding is one of the major natural hazards in Taiwan, and most of the low-lying areas in Taiwan are flood-prone areas. In order to minimize loss of life and economic losses, a detailed and comprehensive decision-making tool is necessary for both flood control planning and emergency service operations. The objectives of this research were (i) to develop a hierarchical structure through the analytic hierarchy process (AHP) to provide preferred options for flood risk analysis, (ii) to map the relative flood risk using the geographic information system (GIS), and (iii) to integrate these two methodologies and apply them to one urban and one semi-rural area in central Taiwan. Fushin Township and the floodplain of Fazih River (1 km on either side of the channel) in Taichung City were selected for this study. In this paper, the flood risk is defined as the relative flood risk due to broken dikes or the failure of stormwater drainage systems. Seven factors were considered in relation to the failure of stormwater drainage, and five to that of broken dikes. Following well-defined procedures, flood maps were drawn based on the data collected from expert responses to a questionnaire, the field survey, satellite images, and documents from flood management agencies. The relative values of flood risk are presented using a 200-m grid for the two study areas. It is concluded that integration of AHP and GIS in flood risk assessment can provide useful detailed information for flood risk management, and the method can be easily applied to most areas in Taiwan where required data sets are readily available. Content Type Journal Article Pages 1-16 DOI 10.1007/s11069-011-9831-7 Authors Yi-Ru Chen, School of Engineering, Griffith University, 170 Kessels Road, Nathan, Brisbane, QLD 4111, Australia Chao-Hsien Yeh, Department of Water Resources Engineering and Conservation, Feng Chia University, No. 100 Wenhwa Road, Seatwen, Taichung, 40724 Taiwan Bofu Yu, School of Engineering, Griffith University, 170 Kessels Road, Nathan, Brisbane, QLD 4111, Australia Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description: The scale issue is of central concern in hydrological processes to understand the potential upscaling or downscaling methodologies, and to develop models for scaling the dominant processes at different scales and in different environments. In this study, a typical permafrost watershed in the Qinghai-Tibet Plateau was selected. Its hydrological processes were monitored for four years from 2004 to 2008; measuring the effects of freezing and thawing depth of active soil layers on runoff processes. To identify the nature and cause of variation in the runoff response in different size catchments, catchments ranging from 1.07 km 2 to 112 km 2 were identified in the watershed. The results indicated that the variation of runoff coefficients showed a “V” shape with increasing catchment size during the spring and autumn seasons, when the active soil was subjected to thawing or freezing processes. A two-stage method was proposed to create runoff scaling models to indicate the effects of scale on runoff processes. In summer, the scaling transition model followed an exponential function for mean daily discharge, whereas the scaling model for flood flow exhibited a linear function. In autumn, the runoff process transition across multiple scales followed an exponential function with air temperature as the driving factor. These scaling models demonstrate relatively high simulation efficiency and precision, and provide a practical way for upscaling or downscaling runoff processes in a medium-size permafrost watershed. For permafrost catchments of this scale, the results show that the synergistic effect of scale and vegetation cover is an important driving factor in the runoff response. Copyright © 2011 John Wiley & Sons, Ltd.

Description: Regression based regional flood frequency analysis (RFFA) methods are widely adopted in hydrology. This paper compares two regression based RFFA methods using a Bayesian Generalized Least Squares (GLS) modelling framework; the two are quantile regression technique (QRT) and parameter regression technique (PRT). In this study, the QRT focuses on the development of prediction equations for a flood quantile in the range of 2 to 100 years average recurrence intervals (ARI), while the PRT develops prediction equations for the first three moments of the log Pearson Type 3 (LP3) distribution, which are the mean, standard deviation and skew of the logarithms of the annual maximum flows; these regional parameters are then used to fit the LP3 distribution to estimate the desired flood quantiles at a given site. It has been shown that using a method similar to stepwise regression and by employing a number of statistics such as the model error variance, average variance of prediction, Bayesian information criterion and Akaike information criterion, the best set of explanatory variables in the GLS regression can be identified. In this study, a range of statistics and diagnostic plots have been adopted to evaluate the regression models. The method has been applied to 53 catchments in Tasmania, Australia. It has been found that catchment area and design rainfall intensity are the most important explanatory variables in predicting flood quantiles using the QRT. For the PRT, a total of four explanatory variables were adopted for predicting the mean, standard deviation and skew. The developed regression models satisfy the underlying model assumptions quite well; of importance, no outlier sites are detected in the plots of the regression diagnostics of the adopted regression equations. Based on ‘one-at-a-time cross validation’ and a number of evaluation statistics, it has been found that for Tasmania the QRT provides more accurate flood quantile estimates for the higher ARIs while the PRT provides relatively better estimates for the smaller ARIs. The RFFA techniques presented here can easily be adapted to other Australian states and countries to derive more accurate regional flood predictions. Copyright © 2011 John Wiley & Sons, Ltd.

Description: Synthetic data have long been employed in hydrology for model development and testing. The objective of this study was to generate a synthetic dataset of hydrologic response with higher spatial and temporal resolution than could presently be obtained in the field, spanning a longer period than the typical duration of monitoring campaigns in experimental catchments. The synthetic dataset was generated for a rangeland catchment with the Integrated Hydrology Model (InHM), and is presented for future use by the community. The InHM boundary-value problem is based upon the previously reported hypothetical reality of Tarrawarra-like hydrologic response. Whereas the emphasis in developing the hypothetical reality was on parameterizing InHM to reproduce observations from the Tarrawarra catchment, the emphasis in generating the synthetic dataset is on developing an internally valid hydrologic-response dataset that extends well beyond the period of observations at Tarrawarra. The synthetic dataset spans eleven years of continuous forcing and response data (e.g., integrated response, distributed fluxes, state variable dynamics). The dataset should be useful for a wide range of problems including evaluation of simple rainfall runoff modeling techniques, design of measurement networks, development of data-assimilation algorithms, and studies on information theory. The dataset is available at: ftp://pangea.stanford.edu/pub/loague/. Copyright © 2011 John Wiley & Sons, Ltd.

Description:    In the framework of a regional landslide susceptibility study in southern Sicily, a test has been carried out in the Tumarrano river basin (about 80 km 2 ) aimed at characterizing its landslide susceptibility conditions by exporting a “source model”, defined and trained inside a limited (about 20 km 2 ) representative sector (the “source area”). Also, the possibility of exploiting Google Earth ™ software and photo-images databank to produce the landslide archives has been checked. The susceptibility model was defined, according to a multivariate geostatistic approach based on the conditional analysis, using unique condition units (UCUs), which were obtained by combining four selected controlling factors: outcropping lithology, steepness, plan curvature and topographic wetness index. The prediction skill of the exported model, trained with 206 landslides, is compared with the one estimated for the whole studied area, by using a complete landslide archive (703 landslides), to see to what extent the largest time/money costs needed are accounted for. The investigated area stretches in the fore-deep sector of southern Sicily, where clayey rocks, mainly referring to the Numidian Flysch and the Terravecchia Formations, largely crop out. The results of the study confirm both the exploitability of Google Earth ™ to produce landslide archive and possibility to adopt in assessing the landslide susceptibility for large basin, a strategy based on the exportation of models trained in limited representative sectors. Content Type Journal Article Pages 1-12 DOI 10.1007/s11069-011-9870-0 Authors D. Costanzo, Dipartimento di Scienze della Terra e del Mare, Università degli Studi di Palermo, Palermo, Italy C. Cappadonia, Dipartimento di Scienze della Terra e del Mare, Università degli Studi di Palermo, Palermo, Italy C. Conoscenti, Dipartimento di Scienze della Terra e del Mare, Università degli Studi di Palermo, Palermo, Italy E. Rotigliano, Dipartimento di Scienze della Terra e del Mare, Università degli Studi di Palermo, Palermo, Italy Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    The convection and planetary boundary layer (PBL) processes play significant role in the genesis and intensification of tropical cyclones (TCs). Several convection and PBL parameterization schemes incorporate these processes in the numerical weather prediction models. Therefore, a systematic intercomparison of performance of parameterization schemes is essential to customize a model. In this context, six combinations of physical parameterization schemes (2 PBL Schemes, YSU and MYJ, and 3 convection schemes, KF, BM, and GD) of WRF-ARW model are employed to obtain the optimum combination for the prediction of TCs over North Indian Ocean. Five cyclones are studied for sensitivity experiments and the out-coming combination is tested on real-time prediction of TCs during 2008. The tracks are also compared with those provided by the operational centers like NCEP, ECMWF, UKMO, NCMRWF, and IMD. It is found that the combination of YSU PBL scheme with KF convection scheme (YKF) provides a better prediction of intensity, track, and rainfall consistently. The average RMSE of intensity (13 hPa in CSLP and 11 m s −1 in 10-m wind), mean track, and landfall errors is found to be least with YKF combination. The equitable threat score (ETS) of YKF combination is more than 0.2 for the prediction of 24-h accumulated rainfall up to 125 mm. The vertical structural characteristics of cyclone inner core also recommend the YKF combination for Indian seas cyclones. In the real-time prediction of 2008 TCs, the 72-, 48-, and 24-h mean track errors are 172, 129, and 155 km and the mean landfall errors are 125, 73, and 66 km, respectively. Compared with the track of leading operational agencies, the WRF model is competing in 24 h (116 km error) and 72 h (166 km) but superior in 48-h (119 km) track forecast. Content Type Journal Article Pages 1-23 DOI 10.1007/s11069-011-9862-0 Authors Krishna K. Osuri, Indian Institute of Technology, Centre for Atmospheric Sciences, Hauz Khas, New Delhi, 110016 India U. C. Mohanty, Indian Institute of Technology, Centre for Atmospheric Sciences, Hauz Khas, New Delhi, 110016 India A. Routray, National Centre for Medium Range Weather Forecasting, A-50, Institutional Area, Sector-62, Noida, UP, India Makarand A. Kulkarni, Indian Institute of Technology, Centre for Atmospheric Sciences, Hauz Khas, New Delhi, 110016 India M. Mohapatra, India Meteorological Department, Lodi Road, New Delhi, India Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description: This work presents a new design of disc infiltrometer, which, associated with a microflowmeter (MF) and a solenoid valve set, makes it possible to automate the infiltration rate ( Q ) measurements at different soil pressure heads (ψ). The MF consists of a 13.8-cm long and 1.5 mm i.d. pipe, with a pressure transducer connecting the two ends of the MF, inserted in a water-flow pipe that connects the Mariotte tube and the water-supply reservoir of the disc infiltrometer. Water flow is calculated from the head losses in the MF. Changes in ψ in the bubble tower, automatically effected when the infiltration rate reaches steady-state, are controlled by a datalogger connected to four solenoid valves. The new design was tested in laboratory and field conditions, and the results showed that the MF allows the soil water infiltration rates to be correctly estimated for different soil characteristics. The solenoid valve set plus datalogger system satisfactorily monitored the changes in ψ and allowed the measurement time to be optimized. Copyright © 2011 John Wiley & Sons, Ltd.

Description: There are several methods for determining the spatial distribution and magnitude of groundwater inputs to streams. We compared results of conventional methods (dye dilution gauging, acoustic Doppler velocimeter (ADV) differential gauging, and geochemical end-member mixing) to Distributed Temperature Sensing (DTS) using a fiber-optic cable installed along 900 m of Nine Mile Creek in Syracuse, New York, USA during low-flow conditions (discharge = 1.4 m 3 s -1 ). With the exception of differential gauging, all methods identified a focused, contaminated groundwater inflow and produced similar groundwater discharge estimates for that point, with a mean of 66.8 Ls -1 between all methods although the precision of these estimates varied. ADV discharge measurement accuracy was reduced by non-ideal conditions and failed to identify, much less quantify, the modest groundwater input, which was only 5% of total stream flow. These results indicate ambient tracers, such as heat and geochemical mixing, can yield spatially and quantitatively refined estimates of relatively modest groundwater inflow even in large rivers. DTS heat tracing, in particular, provided the finest spatial characterization of groundwater inflow, and may be more universally applicable than geochemical methods, for which a distinct and consistent groundwater end member may be more difficult to identify. Copyright © 2011 John Wiley & Sons, Ltd.

Description: Remote sensing estimates of snow water equivalent (SWE) in mountainous areas are subject to large uncertainties. As a prerequisite for testing passive microwave algorithm estimations of SWE, this study aims to collect snow depth (SD) data and provide an understanding of its complex spatial structure as part of the Canadian International Polar Year observations theme. Snow accumulation, redistribution and ablation are controlled by processes that depend on a variety of topographic factors as well as land surface characteristics, which leads us to modelling SD as a function of proxy variables derived from digital elevation model and Landsat data. Field measurements were performed at 3924 locations compromising 184 sites in 50 transects over two years. These measurements were used to predict SD over the study area using a spatial linear mixed-effects model, a model type capable of handling the hierarchical structure of the field data. The model, built using stepwise variable selection, uses as predictor variables transformed elevation, slope, the logarithm of slope, potential incoming solar radiation and its transform; the normalized difference vegetation index, and a transformed tasseled cap brightness from Landsat imagery. A second, simpler model links SD with density giving SWE. The cross-validated root mean squared error of the SD distribution model was 14 cm around an overall mean of 80 cm over a domain of 250 x 250 km. This instantaneous end-of-season peak-accumulation snow map will enable the validation of satellite remote sensing over a generally inaccessible area. Copyright © 2011 John Wiley & Sons, Ltd.

Description:    Anytime that tornadoes interface with people, the results are often grim. The varied characteristics of tornado windfields and the diverse techniques of home construction complicate the expected outcomes of such interactions. Since it is now known that tornadoes do not produce unfathomable winds, engineers have developed techniques for increasing a home’s tornado resistance; thus, increasing the odds of survival of the inhabitants. Even considering these advances, tornadoes still wreak havoc and cause much causality every year. The story of tornado survival is partly chance, partly knowledge, and partly science. This essay utilizes analogies between a popular motion picture and the state of the art of tornado-resistant construction to illustrate the good, the bad, and the ugly facts of tornado survival. First examined are the positive aspects of tornado probability. The odds of experiencing a tornado are so low that massive expenditures are not necessarily cost effective. Next, various techniques of tornado-resistant construction practices are included for new and existing structures. Introduction of equations for predicting debris penetration are compared with various building techniques to describe the viability of tornado protection. The negative aspect of tornado/structure interfaces describes the fact that despite contemporary construction and prevention exercise, the unpredictability of the tornado windfield causes tremendous damage and loss of life. The final section examines the interface between tornadoes and mobile homes. This proves very ugly in every aspect of the word. Statistics and equations confirm the inherent dangers of such an interaction. Also, described is a novel approach to the design of a portable tornado shelter. Content Type Journal Article Pages 1-21 DOI 10.1007/s11069-011-9875-8 Authors Bobby G. McPeak, Texas Tech University, Lubbock, TX 79409, USA Atila Ertas, Texas Tech University, Lubbock, TX 79409, USA Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X

Description:    Risk assessment of natural hazards is often based on the actual or forecast weather situation. For estimating such climate-related risks, it is important to obtain weather data as frequently as possible. One commonly used climate interpolation routine is DAYMET, which in its current form is not able to update its database for periods of less than a year. In this paper, we report the construction of a new climate database with a standard interface and implement a framework for providing daily updated weather data for online daily weather interpolations across regions. We re-implement the interpolation routines from DAYMET to be compliant with the data handling in the new framework. We determine the optimal number of stations used in two possible interpolation routines, assess the error bounds using an independent validation dataset and compare the results with a previous validation study based on the original DAYMET implementation. Mean absolute errors are 1°C for maximum and minimum temperature, 28 mm for precipitation, 3.2 MJ/m² for solar radiation and 1 hPa for vapour pressure deficit, which is in the range of the original DAYMET routine. Finally, we provide an example application of the methodology and derive a fire danger index for a 1 km grid over Austria. Content Type Journal Article Pages 1-14 DOI 10.1007/s11069-011-9880-y Authors Richard Petritsch, Institute of Silviculture, BOKU University of Natural Resources and Life Sciences Vienna, Peter Jordan Str. 82, 1190 Vienna, Austria Hubert Hasenauer, Institute of Silviculture, BOKU University of Natural Resources and Life Sciences Vienna, Peter Jordan Str. 82, 1190 Vienna, Austria Journal Natural Hazards Online ISSN 1573-0840 Print ISSN 0921-030X