ALBERT

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Description: A survey of strong high-frequency sea level oscillations along the US East Coast between 2006 and 2011 Natural Hazards and Earth System Science, 13, 473-482, 2013 Author(s): S. Pasquet, I. Vilibić, and J. Šepić A systematic survey of high-frequency sea level oscillations (

Description: A fuzzy decision making system for building damage map creation using high resolution satellite imagery Natural Hazards and Earth System Science, 13, 455-472, 2013 Author(s): H. Rastiveis, F. Samadzadegan, and P. Reinartz Recent studies have shown high resolution satellite imagery to be a powerful data source for post-earthquake damage assessment of buildings. Manual interpretation of these images, while being a reliable method for finding damaged buildings, is a subjective and time-consuming endeavor, rendering it unviable at times of emergency. The present research, proposes a new state-of-the-art method for automatic damage assessment of buildings using high resolution satellite imagery. In this method, at the first step a set of pre-processing algorithms are performed on the images. Then, extracting a candidate building from both pre- and post-event images, the intact roof part after an earthquake is found. Afterwards, by considering the shape and other structural properties of this roof part with its pre-event condition in a fuzzy inference system, the rate of damage for each candidate building is estimated. The results obtained from evaluation of this algorithm using QuickBird images of the December 2003 Bam, Iran, earthquake prove the ability of this method for post-earthquake damage assessment of buildings.

Description: 19 May 2011 Kütahya – Simav earthquake and evaluation of existing sample RC buildings according to the TEC-2007 criteria Natural Hazards and Earth System Science, 13, 505-522, 2013 Author(s): M. H. Arslan, M. Olgun, M. A. Köroğlu, I. H. Erkan, A. Köken, and O. Tan This study examines the damage caused to reinforced concrete structures by the 2011 earthquake that occurred in Simav, Turkey. The study briefly reports on post-earthquake field observations, tectonic characteristics of the earthquake area, geotechnical characteristics of the field, and seismic characteristics of the earthquake. The main part of the study comprises a field study, material experiments, and performance analyses of two reinforced concrete buildings that survived the earthquake with medium level damage. The building performance was calculated and assessed according to the Turkish Earthquake Code requirements for existing building stock, and recommendations were made based on the findings.

Description: Tsunami damping by mangrove forest: a laboratory study using parameterized trees Natural Hazards and Earth System Science, 13, 483-503, 2013 Author(s): A. Strusińska-Correia, S. Husrin, and H. Oumeraci Tsunami attenuation by coastal vegetation was examined under laboratory conditions for mature mangroves Rhizophora sp. The developed novel tree parameterization concept, accounting for both bio-mechanical and structural tree properties, allowed to substitute the complex tree structure by a simplified tree model of identical hydraulic resistance. The most representative parameterized mangrove model was selected among the tested models with different frontal area and root density, based on hydraulic test results. The selected parameterized tree models were arranged in a forest model of different width and further tested systematically under varying incident tsunami conditions (solitary waves and tsunami bores). The damping performance of the forest models under these two flow regimes was compared in terms of wave height and force envelopes, wave transmission coefficient as well as drag and inertia coefficients. Unlike the previous studies, the results indicate a significant contribution of the foreshore topography to solitary wave energy reduction through wave breaking in comparison to that attributed to the forest itself. A similar rate of tsunami transmission (ca. 20%) was achieved for both flow conditions (solitary waves and tsunami bores) and the widest forest (75 m in prototype) investigated. Drag coefficient C D attributed to the solitary waves tends to be constant ( C D = 1.5) over the investigated range of the Reynolds number.

Description: Comparison between different approaches to modeling shallow landslide susceptibility: a case history in Oltrepo Pavese, Northern Italy Natural Hazards and Earth System Science, 13, 559-573, 2013 Author(s): D. Zizioli, C. Meisina, R. Valentino, and L. Montrasio On the 27 and 28 April 2009, the area of Oltrepo Pavese in northern Italy was affected by a very intense rainfall event that caused a great number of shallow landslides. These instabilities occurred on slopes covered by vineyards or recently formed woodlands and caused damage to many roads and one human loss. Based on aerial photographs taken immediately after the event and field surveys, more than 1600 landslides were detected. After acquiring topographical data, geotechnical properties of the soils and land use, susceptibility analysis on a territorial scale was carried out. In particular, different physically based models were applied to two contiguous sites with the same geological context but different typologies and sizes of shallow landslides. This paper presents the comparison between the ex-post results obtained from the different approaches. On the basis of the observed landslide localizations, the accuracy of the different models was evaluated, and the significant results are highlighted.

Description: Bayesian hierarchical modelling of North Atlantic windiness Natural Hazards and Earth System Science, 13, 545-557, 2013 Author(s): E. Vanem and O. N. Breivik Extreme weather conditions represent serious natural hazards to ship operations and may be the direct cause or contributing factor to maritime accidents. Such severe environmental conditions can be taken into account in ship design and operational windows can be defined that limits hazardous operations to less extreme conditions. Nevertheless, possible changes in the statistics of extreme weather conditions, possibly due to anthropogenic climate change, represent an additional hazard to ship operations that is less straightforward to account for in a consistent way. Obviously, there are large uncertainties as to how future climate change will affect the extreme weather conditions at sea and there is a need for stochastic models that can describe the variability in both space and time at various scales of the environmental conditions. Previously, Bayesian hierarchical space-time models have been developed to describe the variability and complex dependence structures of significant wave height in space and time. These models were found to perform reasonably well and provided some interesting results, in particular, pertaining to long-term trends in the wave climate. In this paper, a similar framework is applied to oceanic windiness and the spatial and temporal variability of the 10-m wind speed over an area in the North Atlantic ocean is investigated. When the results from the model for North Atlantic windiness is compared to the results for significant wave height over the same area, it is interesting to observe that whereas an increasing trend in significant wave height was identified, no statistically significant long-term trend was estimated in windiness. This may indicate that the increase in significant wave height is not due to an increase in locally generated wind waves, but rather to increased swell. This observation is also consistent with studies that have suggested a poleward shift of the main storm tracks.

Description: Real-time flood forecasting coupling different postprocessing techniques of precipitation forecast ensembles with a distributed hydrological model. The case study of may 2008 flood in western Piemonte, Italy Natural Hazards and Earth System Science, 13, 211-220, 2013 Author(s): D. Cane, S. Ghigo, D. Rabuffetti, and M. Milelli In this work, we compare the performance of an hydrological model when driven by probabilistic rain forecast derived from two different post-processing techniques. The region of interest is Piemonte, northwestern Italy, a complex orography area close to the Mediterranean Sea where the forecast are often a challenge for weather models. The May 2008 flood is here used as a case study, and the very dense weather station network allows us for a very good description of the event and initialization of the hydrological model. The ensemble probabilistic forecasts of the rainfall fields are obtained with the Bayesian model averaging, with the classical poor man ensemble approach and with a new technique, the Multimodel SuperEnsemble Dressing. In this case study, the meteo-hydrological chain initialized with the Multimodel SuperEnsemble Dressing is able to provide more valuable discharge ranges with respect to the one initialized with Bayesian model averaging multi-model.

Description: Technical Note: Use of remote sensing for landslide studies in Europe Natural Hazards and Earth System Science, 13, 299-309, 2013 Author(s): V. Tofani, S. Segoni, A. Agostini, F. Catani, and N. Casagli Within the framework of FP7, an EU-funded SafeLand project, a questionnaire was prepared to collect information about the use of remote sensing for landslide study and to evaluate its actual application in landslide detection, mapping and monitoring. The questionnaire was designed using a Google form and was disseminated among end-users and researchers involved in landslide studies in Europe. In total, 49 answers from 17 different European countries were collected. The outcomes showed that landslide detection and mapping is mainly performed with aerial photos, often associated with optical and radar imagery. Concerning landslide monitoring, satellite radars prevail over the other types of data. Remote sensing is mainly used for detection/mapping and monitoring of slides, flows and lateral spreads with a preferably large scale of analysis (1:5000–1:25 000). All the compilers integrate remote sensing data with other thematic data, mainly geological maps, landslide inventory maps and DTMs and derived maps. According to the research and working experience of the compilers, remote sensing is generally considered to have a medium effectiveness/reliability for landslide studies. The results of the questionnaire can contribute to an overall sketch of the use of remote sensing in current landslide studies and show that remote sensing can be considered a powerful and well-established instrument for landslide mapping, monitoring and hazard analysis.

Description: Seismic risk assessment of architectural heritages in Gyeongju considering local site effects Natural Hazards and Earth System Science, 13, 251-262, 2013 Author(s): H.-J. Park, D.-S. Kim, and D.-M. Kim A seismic risk assessment is conducted for cultural heritage sites in Gyeongju, the capital of Korea's ancient Silla Kingdom. Gyeongju, home to UNESCO World Heritage sites, contains remarkable artifacts of Korean Buddhist art. An extensive geotechnical survey including a series of in situ tests is presented, providing pertinent soil profiles for site response analyses on thirty cultural heritage sites. After the shear wave velocity profiles and dynamic material properties were obtained, site response analyses were carried out at each historical site and the amplification characteristics, site period, and response spectrum of the site were determined for the earthquake levels of 2400 yr and 1000 yr return periods based on the Korean seismic hazard map. Response spectrum and corresponding site coefficients obtained from site response analyses considering geologic conditions differ significantly from the current Korean seismic code. This study confirms the importance of site-specific ground response analyses considering local geological conditions. Results are given in the form of the spatial distribution of bedrock depth, site period, and site amplification coefficients, which are particularly valuable in the context of a seismic vulnerability study. This study presents the potential amplification of hazard maps and provides primary data on the seismic risk assessment of each cultural heritage.

Description: Simulating future precipitation extremes in a complex Alpine catchment Natural Hazards and Earth System Science, 13, 263-277, 2013 Author(s): C. Dobler, G. Bürger, and J. Stötter The objectives of the present investigation are (i) to study the effects of climate change on precipitation extremes and (ii) to assess the uncertainty in the climate projections. The investigation is performed on the Lech catchment, located in the Northern Limestone Alps. In order to estimate the uncertainty in the climate projections, two statistical downscaling models as well as a number of global and regional climate models were considered. The downscaling models applied are the Expanded Downscaling (XDS) technique and the Long Ashton Research Station Weather Generator (LARS-WG). The XDS model, which is driven by analyzed or simulated large-scale synoptic fields, has been calibrated using ECMWF-interim reanalysis data and local station data. LARS-WG is controlled through stochastic parameters representing local precipitation variability, which are calibrated from station data only. Changes in precipitation mean and variability as simulated by climate models were then used to perturb the parameters of LARS-WG in order to generate climate change scenarios. In our study we use climate simulations based on the A1B emission scenario. The results show that both downscaling models perform well in reproducing observed precipitation extremes. In general, the results demonstrate that the projections are highly variable. The choice of both the GCM and the downscaling method are found to be essential sources of uncertainty. For spring and autumn, a slight tendency toward an increase in the intensity of future precipitation extremes is obtained, as a number of simulations show statistically significant increases in the intensity of 90th and 99th percentiles of precipitation on wet days as well as the 5- and 20-yr return values.

Description: Fuzzy approach to analysis of flood risk based on variable fuzzy sets and improved information diffusion methods Natural Hazards and Earth System Science, 13, 239-249, 2013 Author(s): Q. Li The predictive analysis of natural disasters and their consequences is challenging because of uncertainties and incomplete data. The present article studies the use of variable fuzzy sets (VFS) and improved information diffusion method (IIDM) to construct a composite method. The proposed method aims to integrate multiple factors and quantification of uncertainties within a consistent system for catastrophic risk assessment. The fuzzy methodology is proposed in the area of flood disaster risk assessment to improve probability estimation. The purpose of the current study is to establish a fuzzy model to evaluate flood risk with incomplete data sets. The results of the example indicate that the methodology is effective and practical; thus, it has the potential to forecast the flood risk in flood risk management.

Description: Review of unprecedented ULF electromagnetic anomalous emissions possibly related to the Wenchuan M S = 8.0 earthquake, on 12 May 2008 Natural Hazards and Earth System Science, 13, 279-286, 2013 Author(s): M. Li, J. Lu, M. Parrot, H. Tan, Y. Chang, X. Zhang, and Y. Wang This work presents ground based ultra low frequency (ULF) electromagnetic field measurements in the frequency range 0.1–10 Hz from January 2007 to December 2008. In this time period a strong earthquake series hits the Wenchuan region with a main shock of magnitude M S = 8.0 on 12 May 2008. The Hebei ULF electromagnetic observation network includes eight observation stations in north China and the observation system named E-EM is employed to record the electric potential difference between two electrodes with an analog automatic real-time continuous pen recorder. First, weak electric signals appeared on 11 October 2007 at Ningjin station, most of which are with relative long periods ~0.4–3 s and unequal amplitudes ~0.5–20 mm. Then, similar signals appeared at Gaobeidian station at the end of October. Abnormal behavior with various time intervals appeared randomly and not every day. At the beginning of April 2008, one and a half months before the Wenchuan M S = 8.0 earthquake, the anomalies were gradually subject to an intensive increase mainly in Gaobeidian SN direction and Ningjin EW direction. The abnormal behavior appeared almost every day and the amplitudes of electric signals, with short periods of ~0.1–0.3 s, enhanced to ~3–30 mm. Qingxian station started to record marginally high frequency signals in SN and EW components in the middle of April. On 9 May, 3 days before the main shock, the amplitude of high frequency information increased sharply at the same time in two components at Gaobeidian station and the maximum amplitude was up to 70 mm, i.e. 1.3 mV m −1 for the electric field. This situation did not stop until 17 May, 5 days after the main event. However, this kind of climax phenomena did not happen at Ningjin station and Qingxian station. Then weak anomalous information lasted about four months again, and strong signals appeared again for a short time before several powerful aftershocks. It is the first time that an abnormity with so large an amplitude and so long a duration time in the observation history of this network though several strong earthquakes were recorded. Furthermore, no obvious interferences have been found during this period. So this event is possibly related to this shock although all these three stations are more than 1300 km away from the Wenchuan earthquake epicenter.

Description: Combining criteria for delineating lahar- and flash-flood-prone hazard and risk zones for the city of Arequipa, Peru Natural Hazards and Earth System Science, 13, 339-360, 2013 Author(s): J.-C. Thouret, G. Enjolras, K. Martelli, O. Santoni, J. A. Luque, M. Nagata, A. Arguedas, and L. Macedo Arequipa, the second largest city in Peru, is exposed to many natural hazards, most notably earthquakes, volcanic eruptions, landslides, lahars (volcanic debris flows), and flash floods. Of these, lahars and flash floods, triggered by occasional torrential rainfall, pose the most frequently occurring hazards that can affect the city and its environs, in particular the areas containing low-income neighbourhoods. This paper presents and discusses criteria for delineating areas prone to flash flood and lahar hazards, which are localized along the usually dry (except for the rainy season) ravines and channels of the Río Chili and its tributaries that dissect the city. Our risk-evaluation study is based mostly on field surveys and mapping, but we also took into account quality and structural integrity of buildings, available socio-economic data, and information gained from interviews with risk-managers officials. In our evaluation of the vulnerability of various parts of the city, in addition to geological and physical parameters, we also took into account selected socio-economic parameters, such as the educational and poverty level of the population, unemployment figures, and population density. In addition, we utilized a criterion of the "isolation factor", based on distances to access emergency resources (hospitals, shelters or safety areas, and water) in each city block. By combining the hazard, vulnerability and exposure criteria, we produced detailed risk-zone maps at the city-block scale, covering the whole city of Arequipa and adjacent suburbs. Not surprisingly, these maps show that the areas at high risk coincide with blocks or districts with populations at low socio-economic levels. Inhabitants at greatest risk are the poor recent immigrants from rural areas who live in unauthorized settlements in the outskirts of the city in the upper parts of the valleys. Such settlements are highly exposed to natural hazards and have little access to vital resources. Our study provides good rationale for the risk zoning of the city, which in turn may be used as an educational tool for better understanding the potential effects of natural hazards and the exposure of the population residing in and around Arequipa. We hope that our work and the risk-zonation maps will provide the impetus and basis for risk-management authorities of the Municipality and the regional government of Arequipa to enforce existing regulations in building in hazardous zones and to adopt an effective long-term strategy to reduce risks from lahar, flash flood, and other natural hazards.

Description: Analyzing resilience of urban networks: a preliminary step towards more flood resilient cities Natural Hazards and Earth System Science, 13, 221-230, 2013 Author(s): S. Lhomme, D. Serre, Y. Diab, and R. Laganier In Europe, river floods have been increasing in frequency and severity. These circumstances require the management of flood risk by integrating new concepts like urban resilience. Nevertheless, urban resilience seems to have no accurate meanings. That is why researchers are primarily concerned with defining resilience. Nevertheless, focus on research object seems to be more important than focus on conceptual debate (Resilience of what? Rather than what is resilience?). Thus the methodology designed here is focused on urban considerations. In fact, a system approach emphasizes technical networks' importance concerning urban resilience. Principles and assumptions applied in this research finally lead to the analysis of how urban networks are able to face natural hazards. In this context, a Web-GIS has been developed for analyzing resistance capacity, absorption capacity and recovery capacity of different technical networks. A first application has been carried out on a French agglomeration in order to analyze road network absorption capacity. This application is very specific but, thanks to this example, it is already possible to highlight the methodology's usefulness.

Description: Effect of heterogeneities on evaluating earthquake triggering of volcanic eruptions Natural Hazards and Earth System Science, 13, 231-237, 2013 Author(s): J. Takekawa, H. Mikada, and T. Goto Recent researches have indicated coupling between volcanic eruptions and earthquakes. Some of them calculated static stress transfer in subsurface induced by the occurrences of earthquakes. Most of their analyses ignored the spatial heterogeneity in subsurface, or only took into account the rigidity layering in the crust. On the other hand, a smaller scale heterogeneity of around hundreds of meters has been suggested by geophysical investigations. It is difficult to reflect that kind of heterogeneity in analysis models because accurate distributions of fluctuation are not well understood in many cases. Thus, the effect of the ignorance of the smaller scale heterogeneity on evaluating the earthquake triggering of volcanic eruptions is also not well understood. In the present study, we investigate the influence of the assumption of homogeneity on evaluating earthquake triggering of volcanic eruptions using finite element simulations. The crust is treated as a stochastic media with different heterogeneous parameters (correlation length and magnitude of velocity perturbation) in our simulations. We adopt exponential and von Karman functions as spatial auto-correlation functions (ACF). In all our simulation results, the ignorance of the smaller scale heterogeneity leads to underestimation of the failure pressure around a chamber wall, which relates to dyke initiation. The magnitude of the velocity perturbation has a larger effect on the tensile failure at the chamber wall than the difference of the ACF and the correlation length. The maximum effect on the failure pressure in all our simulations is about twice larger than that in the homogeneous case. This indicates that the estimation of the earthquake triggering due to static stress transfer should take account of the heterogeneity of around hundreds of meters.

Description: Modelling long-term fire occurrence factors in Spain by accounting for local variations with geographically weighted regression Natural Hazards and Earth System Science, 13, 311-327, 2013 Author(s): J. Martínez-Fernández, E. Chuvieco, and N. Koutsias Humans are responsible for most forest fires in Europe, but anthropogenic factors behind these events are still poorly understood. We tried to identify the driving factors of human-caused fire occurrence in Spain by applying two different statistical approaches. Firstly, assuming stationary processes for the whole country, we created models based on multiple linear regression and binary logistic regression to find factors associated with fire density and fire presence, respectively. Secondly, we used geographically weighted regression (GWR) to better understand and explore the local and regional variations of those factors behind human-caused fire occurrence. The number of human-caused fires occurring within a 25-yr period (1983–2007) was computed for each of the 7638 Spanish mainland municipalities, creating a binary variable (fire/no fire) to develop logistic models, and a continuous variable (fire density) to build standard linear regression models. A total of 383 657 fires were registered in the study dataset. The binary logistic model, which estimates the probability of having/not having a fire, successfully classified 76.4% of the total observations, while the ordinary least squares (OLS) regression model explained 53% of the variation of the fire density patterns (adjusted R 2 = 0.53). Both approaches confirmed, in addition to forest and climatic variables, the importance of variables related with agrarian activities, land abandonment, rural population exodus and developmental processes as underlying factors of fire occurrence. For the GWR approach, the explanatory power of the GW linear model for fire density using an adaptive bandwidth increased from 53% to 67%, while for the GW logistic model the correctly classified observations improved only slightly, from 76.4% to 78.4%, but significantly according to the corrected Akaike Information Criterion (AIC c ), from 3451.19 to 3321.19. The results from GWR indicated a significant spatial variation in the local parameter estimates for all the variables and an important reduction of the autocorrelation in the residuals of the GW linear model. Despite the fitting improvement of local models, GW regression, more than an alternative to "global" or traditional regression modelling, seems to be a valuable complement to explore the non-stationary relationships between the response variable and the explanatory variables. The synergy of global and local modelling provides insights into fire management and policy and helps further our understanding of the fire problem over large areas while at the same time recognizing its local character.

Description: Landslide susceptibility assessment by using a neuro-fuzzy model: a case study in the Rupestrian heritage rich area of Matera Natural Hazards and Earth System Science, 13, 395-407, 2013 Author(s): F. Sdao, D. S. Lioi, S. Pascale, D. Caniani, and I. M. Mancini The complete assessment of landslide susceptibility needs uniformly distributed detailed information on the territory. This information, which is related to the temporal occurrence of landslide phenomena and their causes, is often fragmented and heterogeneous. The present study evaluates the landslide susceptibility map of the Natural Archaeological Park of Matera (Southern Italy) (Sassi and area Rupestrian Churches sites). The assessment of the degree of "spatial hazard" or "susceptibility" was carried out by the spatial prediction regardless of the return time of the events. The evaluation model for the susceptibility presented in this paper is very focused on the use of innovative techniques of artificial intelligence such as Neural Network, Fuzzy Logic and Neuro-fuzzy Network. The method described in this paper is a novel technique based on a neuro-fuzzy system. It is able to train data like neural network and it is able to shape and control uncertain and complex systems like a fuzzy system. This methodology allows us to derive susceptibility maps of the study area. These data are obtained from thematic maps representing the parameters responsible for the instability of the slopes. The parameters used in the analysis are: plan curvature, elevation (DEM), angle and aspect of the slope, lithology, fracture density, kinematic hazard index of planar and wedge sliding and toppling. Moreover, this method is characterized by the network training which uses a training matrix, consisting of input and output training data , which determine the landslide susceptibility. The neuro-fuzzy method was integrated to a sensitivity analysis in order to overcome the uncertainty linked to the used membership functions. The method was compared to the landslide inventory map and was validated by applying three methods: a ROC (Receiver Operating Characteristic) analysis, a confusion matrix and a SCAI method. The developed neuro-fuzzy method showed a good performance in the determination of the landslide susceptibility map.

Description: A digital social network for rapid collection of earthquake disaster information Natural Hazards and Earth System Science, 13, 385-394, 2013 Author(s): J. H. Xu, G. Z. Nie, and X. Xu Acquiring disaster information quickly after an earthquake is crucial for disaster and emergency rescue management. This study examines a digital social network – an earthquake disaster information reporting network – for rapid collection of earthquake disaster information. Based on the network, the disaster information rapid collection method is expounded in this paper. The structure and components of the reporting network are introduced. Then the work principles of the reporting network are discussed, in which the rapid collection of disaster information is realised by using Global System for Mobile Communications (GSM) messages to report the disaster information and Geographic information system (GIS) to analyse and extract useful disaster information. This study introduces some key technologies for the work principles, including the methods of mass sending and receiving of SMS for disaster management, the reporting network grouping management method, brief disaster information codes, and the GIS modelling of the reporting network. Finally, a city earthquake disaster information quick reporting system is developed and with the support of this system the reporting network obtained good results in a real earthquake and earthquake drills. This method is a semi-real time disaster information collection method which extends current SMS based method and meets the need of small and some moderate earthquakes.

Description: Dynamic decision making for dam-break emergency management – Part 2: Application to Tangjiashan landslide dam failure Natural Hazards and Earth System Science, 13, 439-454, 2013 Author(s): M. Peng and L. M. Zhang Tangjiashan landslide dam, which was triggered by the M s = 8.0 Wenchuan earthquake in 2008 in China, threatened 1.2 million people downstream of the dam. All people in Beichuan Town 3.5 km downstream of the dam and 197 thousand people in Mianyang City 85 km downstream of the dam were evacuated 10 days before the breaching of the dam. Making such an important decision under uncertainty was difficult. This paper applied a dynamic decision-making framework for dam-break emergency management (DYDEM) to help rational decision in the emergency management of the Tangjiashan landslide dam. Three stages are identified with different levels of hydrological, geological and social-economic information along the timeline of the landslide dam failure event. The probability of dam failure is taken as a time series. The dam breaching parameters are predicted with a set of empirical models in stage 1 when no soil property information is known, and a physical model in stages 2 and 3 when knowledge of soil properties has been obtained. The flood routing downstream of the dam in these three stages is analyzed to evaluate the population at risk (PAR). The flood consequences, including evacuation costs, flood damage and monetized loss of life, are evaluated as functions of warning time using a human risk analysis model based on Bayesian networks. Finally, dynamic decision analysis is conducted to find the optimal time to evacuate the population at risk with minimum total loss in each of these three stages.

Description: Reasons for large fluctuation of radon and CO 2 levels in a dead-end passage of a karst cave (Postojna Cave, Slovenia) Natural Hazards and Earth System Science, 13, 287-297, 2013 Author(s): A. Gregorič, J. Vaupotič, and F. Gabrovšek Measurements of radon concentration were performed at three geomorphologically different locations in Postojna Cave, Slovenia. In the part of the cave open to visitors, annual average radon activity concentrations of 3255 ± 1190 Bq m −3 and 2315 ± 1019 Bq m −3 were found at the lowest point (LP) and in the Lepe jame (Beautiful Caves, BC), respectively. A much higher average of 25 020 ± 12 653 Bq m −3 was characteristic of the dead-end passage Pisani rov (Gaily Coloured Corridor, GC), in which CO 2 concentration also reached very high values of 4689 ± 294 ppm in summer. Seasonal variations of radon and CO 2 levels in the cave are governed by convective airflow, controlled mainly by the temperature difference between the cave and the outside atmosphere. The following additional sources of radon and CO 2 were considered: (i) flux of geogas from the Earth's crust through fractured rocks (radon and CO 2 source), (ii) clay sediments inside the passage (radon source) and (iii) the soil layer above the cave (radon and CO 2 source).

Description: Damages to unreinforced masonry buildings by the Van earthquakes of 23 October and 9 November 2011 Natural Hazards and Earth System Science, 13, 329-337, 2013 Author(s): Y. S. Tama, A. Solak, N. Çetinkaya, G. Şen, S. Yılmaz, and H. Kaplan Van, a city in Eastern Anatolian Turkey, was hit by two earthquakes with magnitudes of M w = 7.2 and M w = 5.6 in October and November 2011. Both earthquakes caused extensive damage to many buildings. Unreinforced masonry buildings, especially in rural areas, suffered from those earthquakes extensively as in many other cases observed in Turkey during other previous earthquakes. This paper presents a site survey of damaged masonry buildings. Reasons for the wide spread damages are discussed in the paper.

Description: Assessment of the warning system against floods on a rural area: the case of the lower Siret River (Romania) Natural Hazards and Earth System Science, 13, 409-416, 2013 Author(s): F. Salit, L. Zaharia, and G. Beltrando The development of non-structural measures such as an early warning system, across the Europe, in flood risk management, requires a better understanding of the public involved and of the territory threatened. This paper aims to conduct an assessment of early warning and information to people with an analysis of the population's behaviour, presented in a form of an event tree. The objective is to understand the strengths and weaknesses of the warning system during a deadly flood in the lower Siret River (Romania) in 2005 and to demonstrate that each warning system has to be adapted to the territory in which it is effective. The behavioural model aims to determine to what extent the warning system can be improved but also to suggest ways to adapt risk education to the study area.

Description: Climate-induced landslide reactivation at the edge of the Most Basin (Czech Republic) – progress towards better landslide prediction Natural Hazards and Earth System Science, 13, 361-374, 2013 Author(s): J. Burda, F. Hartvich, J. Valenta, V. Smítka, and J. Rybář The catastrophic landslide at Eisenberg in North Bohemia was reactivated during January 2011. This study integrates a range of geoscientific evidence in order to constrain the spatial and temporal development of this reactivation. It has investigated long-term geodetic measurements to assess the morphological development of the site over the last two decades. There is evidence to suggest that, over this period, the site had been subjected to progressive deformation caused by the collapse of an old mine gallery. However, climatic data show that the reactivation itself was triggered by a dramatic rise in the water table induced by rapid snowmelt during a period of winter warming. Furthermore, geomorphological mapping has been used to characterise the morphology of the reactivated landslide and geophysical profiling has been used to analyse its internal structure. The results show that fissures are continuing to develop above the reactivated landslide scarp while highly saturated stiff-fissured claystones provide an incipient slide plane. The application of laser scanning has shown minimal evidence for ongoing landslide activity. It is, however, clear that future landslide events will occur here due to the favourable lithological, structural, and geotechnical conditions. Finally, we propose that future landslide activity at the site may be predicted by the height of water table as this defines theoretical pore pressure at the depth of the shear plane.

Description: Owen Ridge deep-water submarine landslides: implications for tsunami hazard along the Oman coast Natural Hazards and Earth System Science, 13, 417-424, 2013 Author(s): M. Rodriguez, N. Chamot-Rooke, H. Hébert, M. Fournier, and P. Huchon The recent discovery of voluminous submarine landslides along the Owen Ridge may represent a source of tsunami hazard for the nearby Oman coast. We assess the severity of this potential hazard by performing numerical simulations of tsunami generation and propagation from the biggest landslide (40 km 3 in volume) observed along the Owen Ridge. A finite-difference model, assimilating the landslide to a visco-plastic flow, simulates tsunami generation. Computation results show that Salalah city (190 000 inhabitants) is impacted by 2.5 m-high tsunami waves one hour after sediment failure. Higher wave elevation values (4 m) are reached in the low populated Sawqara Bay over 80 min after slide initiation. Although large submarine failures along remote oceanic ridges are infrequent, this study reveals an underestimated source of tsunami hazard in the Arabian Sea.

Description: Temporal and spatial variations in ionospheric electron density profiles over South Africa during strong magnetic storms Natural Hazards and Earth System Science, 13, 375-384, 2013 Author(s): Y. B. Yao, P. Chen, S. Zhang, and J. J. Chen Observations from the South African TrigNet global navigation satellite system (GNSS) and vertical total electron content (VTEC) data from the Jason-1 satellite were used to analyze the variations in ionospheric electron density profiles over South Africa before and after the severe geomagnetic storms on 15 May 2005. Computerized ionospheric tomography (CIT) was used to inverse the 3-D structure of ionospheric electron density and its response to the magnetic storms. Inversion results showed that electron density significantly increased at 10:00 UT, 15 May compared with that at the same period on 14 May. Positive ionospheric storms were observed in the inversion region during the magnetic storms. Jason-1 data show that the VTEC observed on descending orbits on 15 May significantly increased, whereas that on ascending orbits only minimally changed. This finding is identical to the CIT result.

Description: Dynamic decision making for dam-break emergency management – Part 1: Theoretical framework Natural Hazards and Earth System Science, 13, 425-437, 2013 Author(s): M. Peng and L. M. Zhang An evacuation decision for dam breaks is a very serious issue. A late decision may lead to loss of lives and properties, but a very early evacuation will incur unnecessary expenses. This paper presents a risk-based framework of dynamic decision making for dam-break emergency management (DYDEM). The dam-break emergency management in both time scale and space scale is introduced first to define the dynamic decision problem. The probability of dam failure is taken as a stochastic process and estimated using a time-series analysis method. The flood consequences are taken as functions of warning time and evaluated with a human risk analysis model (HURAM) based on Bayesian networks. A decision criterion is suggested to decide whether to evacuate the population at risk (PAR) or to delay the decision. The optimum time for evacuating the PAR is obtained by minimizing the expected total loss, which integrates the time-related probabilities and flood consequences. When a delayed decision is chosen, the decision making can be updated with available new information. A specific dam-break case study is presented in a companion paper to illustrate the application of this framework to complex dam-breaching problems.

Description: Scale model for the confluent area of debris flow and main river: a case study of the Wenjia Gully Natural Hazards and Earth System Science, 13, 3083-3093, 2013 Author(s): J. Zhang, Z. X. Guo, S. Y. Cao, and V. P. Singh When debris flow discharges into the main river, the deposition of debris raises the river bed, occupies the path of water conveyance and damages or even destroys buildings, resulting in considerable economic loss and possibly fatalities. Mathematical models are normally employed to compute debris flow. However, most of these models employ empirical formulae and coefficients and their results are seldom reliable. On the other hand, scale model tests associated with debris flow have seldom been conducted due to the lack of corresponding similarity laws and the difficulty of achieving the grain diameter scale. Focusing on pseudo-one-phase flow, this paper discusses the laws of similarity for the confluence of debris flow and main river and conducts a case study of the debris flow that occurred on 13 August 2010, in the Wenjia Gully, China. After satisfying the roughness scale, the kinematic viscosity coefficient scale, and the momentum ratio scale, it was found that the deposition terrain in the model test is consistent with the one in the prototype.

Description: Shallow landslide's stochastic risk modelling based on the precipitation event of August 2005 in Switzerland: results and implications Natural Hazards and Earth System Science, 13, 3169-3184, 2013 Author(s): P. Nicolet, L. Foresti, O. Caspar, and M. Jaboyedoff Due to their relatively unpredictable characteristics, shallow landslides represent a risk for human infrastructures. Multiple shallow landslides can be triggered by widespread intense precipitation events. The event of August 2005 in Switzerland is used in order to propose a risk model to predict the expected number of landslides based on the precipitation amounts and lithological units. The spatial distribution of rainfall is characterized by merging data coming from operational weather radars and a dense network of rain gauges with an artificial neural network. Lithologies are grouped into four main units, with similar characteristics. Then, from a landslide inventory containing more than 5000 landslides, a probabilistic relation linking the precipitation amount and the lithology to the number of landslides in a 1 km 2 cell, is derived. In a next step, this relation is used to randomly redistribute the landslides using Monte Carlo simulations. The probability for a landslide to reach a building is assessed using stochastic geometry and the damage cost is assessed from the estimated mean damage cost using an exponential distribution to account for the variability. Although the model reproduces well the number of landslides, the number of affected buildings is underestimated. This seems to result from the human influence on landslide occurrence. Such a model might be useful to characterize the risk resulting from shallow landslides and its variability.

Description: Tephra hazard assessment at Mt. Etna (Italy) Natural Hazards and Earth System Science, 13, 3221-3233, 2013 Author(s): S. Scollo, M. Coltelli, C. Bonadonna, and P. Del Carlo In this paper we present a probabilistic hazard assessment for tephra fallout at Mt. Etna (Italy) associated with both short- and long-lived eruptions. Eruptive scenarios and eruption source parameters were defined based on the geological record, while an advection–diffusion–sedimentation model was used to capture the variation in wind speed and direction with time after calibration with the field data. Two different types of eruptions were considered in our analysis: eruptions associated with strong short-lived plumes and eruptions associated with weak long-lived plumes. Our probabilistic approach was based on one eruption scenario for both types and on an eruption range scenario for eruptions producing weak long-lived plumes. Due to the prevailing wind direction, the eastern flanks are the most affected by tephra deposition, with the 122 BC Plinian and 2002–2003 eruptions showing the highest impact both on infrastructures and agriculture.

Description: Sediment transport on the inner shelf off Khao Lak (Andaman Sea, Thailand) during the 2004 Indian Ocean tsunami and former storm events: evidence from foraminiferal transfer functions Natural Hazards and Earth System Science, 13, 3113-3128, 2013 Author(s): Y. Milker, M. Wilken, J. Schumann, D. Sakuna, P. Feldens, K. Schwarzer, and G. Schmiedl We have investigated the benthic foraminiferal fauna from sediment event layers associated with the 2004 Indian Ocean tsunami and former storms that have been retrieved in short sediment cores from offshore environments of the Andaman Sea, off Khao Lak, western Thailand. Species composition and test preservation of the benthic foraminiferal faunas exhibit pronounced changes across the studied sections and provide information on the depositional history of the tsunami layer, particularly on the source water depth of the displaced foraminiferal tests. In order to obtain accurate bathymetric information on sediment provenance, we have mapped the distribution of modern faunas in non-tsunamigenic surface sediments and created a calibration data set for the development of a transfer function. Our quantitative reconstructions revealed that the resuspension of sediment particles by the tsunami wave was restricted to a maximum water depth of approximately 20 m. Similar values were obtained for former storm events, thus impeding an easy distinction of different high-energy events.

Description: Mapping tsunami impacts on land cover and related ecosystem service supply in Phang Nga, Thailand Natural Hazards and Earth System Science, 13, 3095-3111, 2013 Author(s): G. Kaiser, B. Burkhard, H. Römer, S. Sangkaew, R. Graterol, T. Haitook, H. Sterr, and D. Sakuna-Schwartz The 2004 Indian Ocean tsunami caused damages to coastal ecosystems and thus affected the livelihoods of the coastal communities who depend on services provided by these ecosystems. The paper presents a case study on evaluating and mapping the spatial and temporal impacts of the tsunami on land use and land cover (LULC) and related ecosystem service supply in the Phang Nga province, Thailand. The method includes local stakeholder interviews, field investigations, remote-sensing techniques, and GIS. Results provide an ecosystem services matrix with capacity scores for 18 LULC classes and 17 ecosystem functions and services as well as pre-/post-tsunami and recovery maps indicating changes in the ecosystem service supply capacities in the study area. Local stakeholder interviews revealed that mangroves, casuarina forest, mixed beach forest, coral reefs, tidal inlets, as well as wetlands (peat swamp forest) have the highest capacity to supply ecosystem services, while e.g. plantations have a lower capacity. The remote-sensing based damage and recovery analysis showed a loss of the ecosystem service supply capacities in almost all LULC classes for most of the services due to the tsunami. A fast recovery of LULC and related ecosystem service supply capacities within one year could be observed for e.g. beaches, while mangroves or casuarina forest needed several years to recover. Applying multi-temporal mapping the spatial variations of recovery could be visualised. While some patches of coastal forest were fully recovered after 3 yr, other patches were still affected and thus had a reduced capacity to supply ecosystem services. The ecosystem services maps can be used to quantify ecological values and their spatial distribution in the framework of a tsunami risk assessment. Beyond that they are considered to be a useful tool for spatial analysis in coastal risk management in Phang Nga.

Description: Calibration of a real-time tsunami detection algorithm for sites with no instrumental tsunami records: application to coastal tide-gauge stations in eastern Sicily, Italy Natural Hazards and Earth System Science, 13, 3129-3144, 2013 Author(s): L. Bressan, F. Zaniboni, and S. Tinti Coastal tide gauges play a very important role in a tsunami warning system, since sea-level data are needed for a correct evaluation of the tsunami threat, and the tsunami arrival has to be recognized as early as possible. Real-time tsunami detection algorithms serve this purpose. For an efficient detection, they have to be calibrated and adapted to the specific local characteristics of the site where they are installed, which is easily done when the station has recorded a sufficiently large number of tsunamis. In this case the recorded database can be used to select the best set of parameters enhancing the discrimination power of the algorithm and minimizing the detection time. This chance is however rare, since most of the coastal tide-gauge stations, either historical or of new installation, have recorded only a few tsunamis in their lifetimes, if any. In this case calibration must be carried out by using synthetic tsunami signals, which poses the problem of how to generate them and how to use them. This paper investigates this issue and proposes a calibration approach by using as an example a specific case, which is the calibration of a real-time detection algorithm called TEDA (Tsunami Early Detection Algorithm) for two stations (namely Tremestieri and Catania) in eastern Sicily, Italy, which were recently installed in the frame of the Italian project TSUNET, aiming at improving the tsunami monitoring capacity in a region that is one of the most hazardous tsunami areas of Italy and of the Mediterranean.

Description: Assessment of flash floods taking into account climate change scenarios in the Llobregat River basin Natural Hazards and Earth System Science, 13, 3145-3156, 2013 Author(s): M. Velasco, P. A. Versini, A. Cabello, and A. Barrera-Escoda Global change may imply important changes in the future occurrence and intensity of extreme events. Climate scenarios characterizing these plausible changes were previously obtained for the Llobregat River basin (NE Spain). This paper presents the implementation of these scenarios in the HBV (Hydrologiska Byråns Vattenbalansavdelning) hydrological model. Then, the expected changes in terms of flash flood occurrence and intensity are assessed for two different sub-basins: the Alt Llobregat and the Anoia (Llobregat River basin). The assessment of future flash floods has been done in terms of the intensity and occurrence of extreme events, using a peak over threshold (POT) analysis. For these two sub-basins, most of the simulated scenarios present an increase of the intensity of the peak discharge values. On the other hand, the future occurrence follows different trends in the two sub-basins: an increase is observed in Alt Llobregat but a decrease occurs in Anoia. Despite the uncertainties that appear in the whole process, the results obtained can shed some light on how future flash floods events may occur.

Description: Sea cliff instability susceptibility at regional scale: a statistically based assessment in the southern Algarve, Portugal Natural Hazards and Earth System Science, 13, 3185-3203, 2013 Author(s): F. M. S. F. Marques, R. Matildes, and P. Redweik Sea cliff evolution is dominated by the occurrence of slope mass movements of different types and sizes, which are a considerable source of natural hazard, making their assessment a relevant issue in terms of human loss prevention and land use regulations. To address the assessment of the spatial component of sea cliff hazards, i.e. the susceptibility, a statistically based study was made to assess the capacity of a set of conditioning factors to express the occurrence of sea cliff failures affecting areas located along their top. The study was based on the application of the bivariate information value and multivariate logistic regression statistical methods, using a set of predisposing factors for cliff failures, mainly related to geology (lithology, bedding dip, faults) and geomorphology (maximum and mean slope, height, aspect, plan curvature, toe protection), which were correlated with a photogrammetry-based inventory of cliff failures that occurred in a 60 yr period (1947–2007). The susceptibility models were validated against the inventory data using standard success rate and ROC curves, and provided encouraging results, indicating that the proposed approaches are effective for susceptibility assessment. The results obtained also stress the need for improvement of the predisposing factors to be used in this type of study and the need for detailed and systematic cliff failure inventories.

Description: Application and reliability of techniques for landslide site investigation, monitoring and early warning – outcomes from a questionnaire study Natural Hazards and Earth System Science, 13, 3157-3168, 2013 Author(s): I. Baroň and R. Supper The presented questionnaire study summarizes an evaluation of approaches, techniques and parameters of slope-instability investigation and monitoring of their occurrence, reliability and the applicability of the monitoring techniques for early warning. The study is based on information collected from 86 monitored landslides in 14 European and Asian countries. Based on the responses, lidar ALS (airborne laser scanners), geophysical logging, aerial photographs, resistivity surveying, GB InSAR (ground-based synthetic aperture radar interferometer) and the refraction seismic were considered the most reliable methods for investigation of structure and character of landslides. Especially lidar ALS and geophysical logging were ranked high despite their application at relatively few landslides. Precipitation amount, pore-water pressure and displacement monitored by wire extensometers, dGPS and total stations, followed by air temperature and EM-emissions monitoring and displacement monitored by the TM 71 crack gauge were considered the most promising parameters for early warning.

Description: Novel method for hurricane trajectory prediction based on data mining Natural Hazards and Earth System Science, 13, 3211-3220, 2013 Author(s): X. Dong and D. C. Pi This paper describes a novel method for hurricane trajectory prediction based on data mining (HTPDM) according to the hurricane's motion characteristics. Firstly, all frequent trajectories in the historical hurricane trajectory database are mined by using association analysis technology and their corresponding association rules are generated as motion patterns. Then, the current hurricane trajectories are matched with the motion patterns for predicting. If no association rule is found for matching, a predicted result according to the hurricane current movement trend would be returned. All experiments are conducted with the Atlantic weather Hurricane/Tropical Data from 1900 to 2008. The experimental results show that if the matching failure part is contained, the prediction accuracy is 57.5%. Whereas, the valve would be to 65% provided all matches are successful.

Description: Potential and limitations of risk scenario tools in volcanic areas through an example at Mount Cameroon Natural Hazards and Earth System Science, 13, 2409-2424, 2013 Author(s): P. Gehl, C. Quinet, G. Le Cozannet, E. Kouokam, and P. Thierry This paper presents an integrated approach to conduct a scenario-based volcanic risk assessment on a variety of exposed assets, such as residential buildings, cultivated areas, network infrastructures or individual strategic buildings. The focus is put on the simulation of scenarios, based on deterministic adverse event input, which are applied to the case study of an effusive eruption on the Mount Cameroon volcano, resulting in the damage estimation of the assets located in the area. The work is based on the recent advances in the field of seismic risk. A software for systemic risk scenario analysis developed within the FP7 project SYNER-G has been adapted to address the issue of volcanic risk. Most significant improvements include the addition of vulnerability models adapted to each kind of exposed element and the possibility to quantify the successive potential damages inflicted by a sequence of adverse events (e.g. lava flows, tephra fall, etc.). The use of an object-oriented architecture gives the opportunity to model and compute the physical damage of very disparate types of infrastructures under the same framework. Finally, while the risk scenario approach is limited to the assessment of the physical impact of adverse events, a specific focus on strategic infrastructures and a dialogue with stakeholders helps in evaluating the potential wider indirect consequences of an eruption.

Description: Simulation systems for tsunami wave propagation forecasting within the French tsunami warning center Natural Hazards and Earth System Science, 13, 2465-2482, 2013 Author(s): A. Gailler, H. Hébert, A. Loevenbruck, and B. Hernandez A model-based tsunami prediction system has been developed as part of the French Tsunami Warning Center (operational since 1 July 2012). It involves a precomputed unit source functions database (i.e., a number of tsunami model runs that are calculated ahead of time and stored). For the Mediterranean basin, the faults of the unit functions are placed adjacent to each other, following the discretization of the main seismogenic faults. An automated composite scenarios calculation tool is implemented to allow the simulation of any tsunami propagation scenario (i.e., of any seismic moment). Uncertainty on the magnitude of the detected event and inaccuracy of the epicenter location are taken into account in the composite scenarios calculation. Together with this forecasting system, another operational tool based on real time computing is implemented as part of the French Tsunami Warning Center. This second tsunami simulation tool takes advantage of multiprocessor approaches and more realistic seismological parameters, once the focal mechanism is established. Three examples of historical earthquakes are presented, providing warning refinement compared to the rough tsunami risk map given by the model-based decision matrix.

Description: Daily precipitation records over mainland Spain and the Balearic Islands Natural Hazards and Earth System Science, 13, 2483-2491, 2013 Author(s): C. Ramis, V. Homar, A. Amengual, R. Romero, and S. Alonso Understanding the spatial distribution of extreme precipitations is of major interest in order to improve our knowledge of the climate of a region and its relationship with society. These analyses inevitably require the use of directly observed values to account for the actual extreme amounts rather than analyzed gridded values. A study of daily rainfall extremes observed over mainland Spain and the Balearic Islands is performed by using records from 8135 rain gauge stations from the Spanish Weather Agency (AEMET). Results show that the heaviest daily precipitations have been observed mainly on the coastal Mediterranean zone from Gibraltar to the Pyrenees. In particular, a record value of 817 mm was recorded in the Valencia region in 1987. The current map of daily records in Spain, which updates the pioneering work of the Spanish meteorologist Font, shows similar distribution of extreme events but with notably higher amounts. Generalized extreme values distributions fit the Mediterranean and Atlantic rain gauge measurements and shows the different characteristics of the extreme daily precipitations in both regions. We identify the most extreme events (above 500 mm per day) and provide a brief description of a typical meteorological situation in which these damaging events occur. An analysis of the low-level circulation patterns producing such extremes – by means of simple indices such as NAO, WeMOi and IBEI – confirms the relevance of local flows in the generation of either Mediterranean or Atlantic episodes. WeMOi, and even more IBEI, are good discriminants of the region affected by the record precipitation event.

Description: Seismicity at the northeast edge of the Mexican Volcanic Belt (MVB) and activation of an undocumented fault: the Peñamiller earthquake sequence of 2010–2011, Querétaro, Mexico Natural Hazards and Earth System Science, 13, 2521-2531, 2013 Author(s): A. Clemente-Chavez, A. Figueroa-Soto, F. R. Zúñiga, M. Arroyo, M. Montiel, and O. Chavez The town of Peñamiller in the state of Querétaro, Mexico, is located at the northeast border of the seismogenic zone known as the Mexican Volcanic Belt (MVB), which transects the central part of Mexico with an east–west orientation. In the vicinity of this town, a sequence of small earthquakes occurred during the end of 2010 and beginning of 2011. Seismicity in the continental regimen of central Mexico is not too frequent; however, it is known that there are precedents of large earthquakes ( M w magnitude greater than 6.0) occurring in this zone. Three large earthquakes have occurred in the past 100 yr: the 19 November 1912 ( M S = 7.0), the 3 January 1920 ( M S = 6.4), and the 29 June 1935 ( M S = 6.9) earthquakes. Prior to the instrumental period, the earthquake of 11 February 1875, which took place near the city of Guadalajara, caused widespread damage. The purpose of this article is to contribute to the available seismic information of this region. This will help advance our understanding of the tectonic situation of the central Mexico MVB region. Twenty-four shallow earthquakes of the Peñamiller seismic sequence of 2011 were recorded by a temporary accelerograph network installed by the Universidad Autónoma de Querétaro (UAQ). The data were analyzed in order to determine the source locations and to estimate the source parameters. The study was carried out through an inversion process and by spectral analysis. The results show that the largest earthquake occurred on 8 February 2011 at 19:53:48.6 UTC, had a moment magnitude M w = 3.5, and was located at latitude 21.039° and longitude −99.752°, at a depth of 5.6 km. This location is less than 7 km away in a south-east direction from downtown Peñamiller. The focal mechanisms are mostly normal faults with small lateral components. These focal mechanisms are consistent with the extensional regimen of the southern extension of the Basin and Range (BR) province. The source area of the largest event was estimated to have a radius of 0.5 km, which corresponds to a normal fault with azimuth of 174° and an almost pure dip slip. Peak ground acceleration (PGA) was close to 100 cm s −2 in the horizontal direction. Shallow earthquakes induced by crustal faulting present a potential seismic risk and hazard within the MVB, considering the population growth. Thus, the necessity to enrich seismic information in this zone is very important since the risk at most urban sites in the region might even be greater than that posed by subduction earthquakes.

Description: Anomalies of total column CO and O 3 associated with great earthquakes in recent years Natural Hazards and Earth System Science, 13, 2513-2519, 2013 Author(s): Y. Cui, J. Du, D. Zhang, and Y. Sun Variations of total column CO and O 3 in the atmosphere over the epicenter areas of 35 great earthquakes that occurred throughout the world in recent years were studied based on the hyper-spectrum data from Atmospheric Infrared Sounder (AIRS). It was found that anomalous increases of CO and/or O 3 concentrations occurred over the epicenter areas of 12 earthquakes among the 35 studied ones. However, increases in both CO and O 3 concentrations were found for 6 earthquakes. The O 3 anomalies appeared in the month when the earthquake occurred and lasted for a few months, whereas CO anomalies occurred irregularly. The duration of CO and O 3 anomalies related to the earthquakes ranged from 1 to 6 months. The anomalies of CO concentration related to the earthquake can be mainly attributed to gas emission from the lithosphere and photochemical reaction, while the anomalous increases in O 3 concentration can be mainly due to the transport of O 3 -enriched air and photochemical reaction. However, more work needs to be done in order to understand the mechanism of the CO and O 3 anomalies further.

Description: Review Article: Economic evaluation of flood damage to agriculture – review and analysis of existing methods Natural Hazards and Earth System Science, 13, 2493-2512, 2013 Author(s): P. Brémond, F. Grelot, and A.-L. Agenais In Europe, economic evaluation of flood management projects is increasingly used to help decision making. At the same time, the management of flood risk is shifting towards new concepts such as giving more room to water by restoring floodplains. Agricultural areas are particularly targeted by projects following those concepts since they are frequently located in floodplain areas and since the potential damage to such areas is expected to be lower than to cities or industries for example. Additional or avoided damage to agriculture may have a major influence on decisions concerning these projects and the economic evaluation of flood damage to agriculture is thus an issue that needs to be tackled. The question of flood damage to agriculture can be addressed in different ways. This paper reviews and analyzes existing studies which have developed or used damage functions for agriculture in the framework of an economic appraisal of flood management projects. A conceptual framework of damage categories is proposed for the agricultural sector. The damage categories were used to structure the review. Then, a total of 42 studies are described, with a detailed review of 26 of them, based on the following criteria: types of damage considered, the influential flood parameters chosen, and monetized damage indicators used. The main recommendations resulting from this review are that even if existing methods have already focused on damage to crops, still some improvement is needed for crop damage functions. There is also a need to develop damage functions for other agricultural damage categories, including farm buildings and their contents. Finally, to cover all possible agricultural damage, and in particular loss of activity, a farm scale approach needs to be used.

Description: Are great Cascadia earthquakes recorded in the sedimentary records from small forearc lakes? Natural Hazards and Earth System Science, 13, 2441-2463, 2013 Author(s): A. E. Morey, C. Goldfinger, C. E. Briles, D. G. Gavin, D. Colombaroli, and J. E. Kusler Here we investigate sedimentary records from four small inland lakes located in the southern Cascadia forearc region for evidence of earthquakes. Three of these lakes are in the Klamath Mountains near the Oregon–California border, and one is in the central Oregon Coast range. The sedimentary sequences recovered from these lakes are composed of normal lake sediment interbedded with disturbance event layers. The thickest of these layers are graded, and appear to be turbidites or linked debrites (turbidites with a basal debris-flow deposit), suggesting rapid deposition. Variations in particle size and organic content of these layers are reflected in the density and magnetic susceptibility data. The frequency and timing of these events, based on radiocarbon ages from detrital organics, is similar to the offshore seismogenic turbidite record from trench and slope basin cores along the Cascadia margin. Stratigraphic correlation of these anomalous deposits based on radiocarbon ages, down-core density, and magnetic susceptibility data between lake and offshore records suggest synchronous triggering. The areal extent and multiple depositional environments over which these events appear to correlate suggest that these deposits were most likely caused by shaking during great Cascadia earthquakes.

Description: Advanced interpretation of land subsidence by validating multi-interferometric SAR data: the case study of the Anthemountas basin (Northern Greece) Natural Hazards and Earth System Science, 13, 2425-2440, 2013 Author(s): F. Raspini, C. Loupasakis, D. Rozos, and S. Moretti The potential of repeat-pass space borne SAR (Synthetic Aperture Radar) interferometry has been exploited to investigate spatial patterns of land subsidence in the Anthemountas basin, in the northern part of Greece. The PSI (Persistent Scatterer Interferometry) approach, based on the processing of long series of SAR acquisitions, has been applied to forty-two images acquired in 1995–2001 by ERS1/2 satellites. Interferometric results have been analysed at a basin scale as support for land motion mapping and at a local scale for the characterisation of ground motion events affecting the village of Perea in the Thermaikos municipality and the "Macedonia" international airport. PSI results revealed a moderate subsidence phenomenon along the wider coastal zone of Anthemountas basin corresponding to intense groundwater extraction. Highest values, exceeding −20 mm yr −1 , were measured in the airport area where the thickest sequence of compressible Quaternary sediments occurs. Intense subsidence has been detected also in the Perea village (maximum deformation of −10 to −15 mm yr −1 ), where a series of fractures, causing damages to both buildings and infrastructure, occurred in 2005–2006.

Description: A shallow-flow model for the propagation of tsunamis over complex geometries and mobile beds Natural Hazards and Earth System Science, 13, 2533-2542, 2013 Author(s): D. A. S. Conde, M. A. V. Baptista, C. Sousa Oliveira, and R. M. L. Ferreira A distinguishable feature of overland tsunami propagation is the incorporation of solids within the flow column, either sediment from the natural environment or remains from built infrastructure. This article describes a 2DH (two-dimensional horizontal) mathematical model particularly suited for tsunami propagation over complex and dynamic geometries, such as river and estuarine mobile beds. The discretization scheme is based on a finite-volume method using a flux-splitting technique featuring a reviewed Roe–Riemann solver, with appropriate source-term formulations to ensure full conservativeness. The model is validated with laboratory data and paleo-tsunami evidence. As a forecasting application, it is applied to a tsunami scenario in the Tagus estuary, an effort justified by the numerous catastrophic tsunamis that are known to have struck this location over the past two millennia. The obtained results show that, despite the significant differences in Lisbon's layout and morphology, a 1755-like tsunami would still inflict a devastating impact on this major city.

Description: Book Review: The Role of Ecosystems in Disaster Risk Reduction Natural Hazards and Earth System Science, 13, 2553-2554, 2013 Author(s): P. Tarolli

Description: Airborne geophysical survey of the catastrophic landslide at Stože, Log pod Mangrtom, as a test of an innovative approach for landslide mapping in steep alpine terrains Natural Hazards and Earth System Science, 13, 2543-2550, 2013 Author(s): I. Baroň, R. Supper, E. Winkler, K. Motschka, A. Ahl, M. Čarman, and Š. Kumelj Airborne geophysics is a promising method for investigating landslides. Here we present a case study of multisensor airborne geophysical survey at the catastrophic landslide Stože near Log pod Mangrtom in Slovenia, which was conducted in the framework of the European FP7th Project "SafeLand". Based on the survey itself and achieved results, we discuss applicability, limits, and benefits and costs of the method for investigating landslides in steep alpine terrains. Despite of several operational constraints, the airborne electromagnetic survey of the area well presented the lithological pattern and water saturation. The high resistivity regions mostly indicated drained slope scree and landslide mass, drained and loosened material of the moraine deposit in the tension zone of the landslide with present cracks and cavities. The minima of the resistivity pattern were attributed to the outcrop of marls rich in clay, to water-saturated moraine deposit above impermeable marls in the tension zone, and to water-saturated porous alluvial gravel and landslide scree along the Koritnica River. The magnetic survey proved to be inapplicable for such a small and rough area. The Potassium and Thorium maps, on the other hand, both well identified the regions of tension inside the landslide zone, outcrops of marls and dolomite, clay-rich colluvium, weathered zones along a regional tectonic fault, and alluvial deposits and deposits of debris flows, and the minima of the 137 Cs clearly revealed the zones of material removal due to recent mass movements.

Description: Book Review: Natural Hazards in the Asia–Pacific Region: Recent Advances and Emerging Concepts Natural Hazards and Earth System Science, 13, 2551-2552, 2013 Author(s): P. Tarolli

Description: Effects of relative density and accumulated shear strain on post-liquefaction residual deformation Natural Hazards and Earth System Science, 13, 2567-2577, 2013 Author(s): J. Kim, M. Kazama, and Y. Kwon The damage caused by liquefaction, which occurs following an earthquake, is usually because of settlement and lateral spreading. Generally, the evaluation of liquefaction has been centered on settlement, that is, residual volumetric strain. However, in actual soil, residual shear and residual volumetric deformations occur simultaneously after an earthquake. Therefore, the simultaneous evaluation of the two phenomena and the clarification of their relationship are likely to evaluate post-liquefaction soil behaviors more accurately. Hence, a quantitative evaluation of post-liquefaction damage will also be possible. In this study, the effects of relative density and accumulated shear strain on post-liquefaction residual deformations were reviewed through a series of lateral constrained-control hollow cylindrical torsion tests under undrained conditions. In order to identify the relationship between residual shear and residual volumetric strains, this study proposed a new test method that integrates monotonic loading after cyclic loading, and K 0 -drain after cyclic loading – in other words, the combination of cyclic loading, monotonic loading, and the K 0 drain. In addition, a control that maintained the lateral constrained condition across all the processes of consolidation, cyclic loading, monotonic loading, and drainage was used to reproduce the anisotropy of in situ ground. This lateral constrain control was performed by controlling the axial strain, based on the assumption that under undrained conditions, axial and lateral strains occur simultaneously, and unless axial strain occurs, lateral strain does not occur. The test results confirmed that the recovery of effective stresses, which occur during monotonic loading and drainage after cyclic loading, respectively, result from mutually different structural restoration characteristics. In addition, in the ranges of 40–60% relative density and 50–100% accumulated shear strain, relative density was found to have greater effects than the number of cycles (accumulated shear strain).

Description: Challenges to social capacity building in flood-affected areas of southern Poland Natural Hazards and Earth System Science, 13, 2555-2566, 2013 Author(s): J. Działek, W. Biernacki, and A. Bokwa Various aspects of beliefs, behaviour and expectations of at-risk populations were analysed in four case study localities in southern Poland that were affected by flooding in 1997 and 2001. They represent localities of different sizes and are characterised by different paths of historical development. Two of them are deep-rooted communities with dense, strong family and neighbourhood ties, while the other two experienced an almost total replacement of their population due to decisions taken after World War II and still suffer from less developed social networks. Historical events also resulted in the disruption of local memories of flooding and transmission of knowledge about natural hazards. A questionnaire survey was conducted in late autumn 2006, followed by structured telephone interviews and focus group interviews in spring 2008. The results of the survey and interviews were analysed with reference to the social capacity framework and its five dimensions: knowledge, motivational, network, economic and governance capacities. Network capacities, that is resources of bonding and bridging social capital, were considered a key notion when analysing and interpreting the results. The differences in the local resources and abilities available in each of the localities to prepare a response to natural hazards were revealed. Consequently, challenges faced in the process of building and strengthening social capacity were identified as well as ways to address these challenges. It was concluded that there are general trends and tendencies that need to be considered in risk management strategies, however the different starting points of each case study community calls for different means and approaches, as well as producing somewhat different expected outcomes.

Description: Precursory pattern of tidal triggering of earthquakes in six regions of China: the possible relation to the crustal heterogeneity Natural Hazards and Earth System Science, 13, 2605-2618, 2013 Author(s): Q. Li and G.-M. Xu We found the possible correlation between the precursory pattern of tidal triggering of earthquakes and the crustal heterogeneities, which is of particular importance to the researchers in earthquake prediction and earthquake hazard prevention. We investigated the connection between the tidal variations and earthquake occurrence in the Liyang, Wunansha, Cangshan, Wenan, Luquan and Yaoan regions of China. Most of the regions show a higher correlation with tidal triggering in several years preceding the large or destructive earthquakes compared to other times, indicating that the tidal triggering may inherently relate to the nucleation of the destructive earthquakes during this time. In addition, the analysis results indicate that the Liyang, Cangshan and Luquan regions, with stronger heterogeneity, show statistically significant effects of tidal triggering preceding the large or destructive earthquakes, while the Wunansha, Wenan and Yaoan regions, with relatively weak heterogeneity, show statistically insignificant effects of it, signifying that the precursory pattern of tidal triggering of earthquakes in these six regions is possibly related to the heterogeneities of the crustal rocks. The above results suggest that when people try to find the potential earthquake hazardous areas or make middle–long-term earthquake forecasting by means of precursory pattern of the tidal triggering, the crustal heterogeneity in these areas has to be taken into consideration for the purpose of increasing the prediction efficiency. If they do not consider the influence of crustal heterogeneity on the tidal triggering of earthquakes, the prediction efficiency might greatly decrease.

Description: A preliminary evaluation of surface latent heat flux as an earthquake precursor Natural Hazards and Earth System Science, 13, 2639-2647, 2013 Author(s): W. Zhang, J. Zhao, W. Wang, H. Ren, L. Chen, and G. Yan The relationship between variations in surface latent heat flux (SLHF) and marine earthquakes has been a popular subject of recent seismological studies. So far, there are two key problems: how to identify the abnormal SLHF variations from complicated background signals, and how to ensure that the anomaly results from an earthquake. In this paper, we proposed four adjustable parameters for identification, classified the relationship and analyzed SLHF changes several months before six marine earthquakes by employing daily SLHF data. Additionally, we also quantitatively evaluate the long-term relationship between earthquakes and SLHF anomalies for the six study areas over a 20 yr period preceding each earthquake. The results suggest the following: (1) before the South Sandwich Islands, Papua, Samoa and Haiti earthquakes, the SLHF variations above their individual background levels have relatively low amplitudes and are difficult to be considered as precursory anomalies; (2) after removing the clustering effect, most of the anomalies prior to these six earthquakes are not temporally related to any earthquake in each study area in time sequence; (3) for each case, apart from Haiti, more than half of the studied earthquakes, which were moderate and even devastating earthquakes (larger than M w = 5.3), had no precursory variations in SLHF; and (4) the correlation between SLHF and seismic activity depends largely on data accuracy and parameter settings. Before any application of SLHF data on earthquake prediction, we suggest that anomaly-identifying standards should be established based on long-term regional analysis to eliminate subjectivity. Furthermore, other factors that may result in SLHF variations should also be carefully considered.

Description: The price of safety: costs for mitigating and coping with Alpine hazards Natural Hazards and Earth System Science, 13, 2619-2637, 2013 Author(s): C. Pfurtscheller and A. H. Thieken Due to limited public budgets and the need to economize, the analysis of costs of hazard mitigation and emergency management of natural hazards becomes increasingly important for public natural hazard and risk management. In recent years there has been a growing body of literature on the estimation of losses which supported to help to determine benefits of measures in terms of prevented losses. On the contrary, the costs of mitigation are hardly addressed. This paper thus aims to shed some light on expenses for mitigation and emergency services. For this, we analysed the annual costs of mitigation efforts in four regions/countries of the Alpine Arc: Bavaria (Germany), Tyrol (Austria), South Tyrol (Italy) and Switzerland. On the basis of PPP values (purchasing power parities), annual expenses on public safety ranged from EUR 44 per capita in the Free State of Bavaria to EUR 216 in the Autonomous Province of South Tyrol. To analyse the (variable) costs for emergency services in case of an event, we used detailed data from the 2005 floods in the Federal State of Tyrol (Austria) as well as aggregated data from the 2002 floods in Germany. The analysis revealed that multi-hazards, the occurrence and intermixture of different natural hazard processes, contribute to increasing emergency costs. Based on these findings, research gaps and recommendations for costing Alpine natural hazards are discussed.

Description: Experiences from site-specific landslide early warning systems Natural Hazards and Earth System Science, 13, 2659-2673, 2013 Author(s): C. Michoud, S. Bazin, L. H. Blikra, M.-H. Derron, and M. Jaboyedoff Landslide early warning systems (EWSs) have to be implemented in areas with large risk for populations or infrastructures when classical structural remediation measures cannot be set up. This paper aims to gather experiences of existing landslide EWSs, with a special focus on practical requirements (e.g., alarm threshold values have to take into account the smallest detectable signal levels of deployed sensors before being established) and specific issues when dealing with system implementations. Within the framework of the SafeLand European project, a questionnaire was sent to about one-hundred institutions in charge of landslide management. Finally, we interpreted answers from experts belonging to 14 operational units related to 23 monitored landslides. Although no standard requirements exist for designing and operating EWSs, this review highlights some key elements, such as the importance of pre-investigation work, the redundancy and robustness of monitoring systems, the establishment of different scenarios adapted to gradual increasing of alert levels, and the necessity of confidence and trust between local populations and scientists. Moreover, it also confirms the need to improve our capabilities for failure forecasting, monitoring techniques and integration of water processes into landslide conceptual models.

Description: Collisions of two breathers at the surface of deep water Natural Hazards and Earth System Science, 13, 3205-3210, 2013 Author(s): A. I. Dyachenko, D. I. Kachulin, and V. E. Zakharov We present results of numerical experiments on long-term evolution and collisions of breathers (which correspond to envelope solitons in the NLSE approximation) at the surface of deep ideal fluid. The collisions happen to be nonelastic. In the numerical experiment it can be observed only after many acts of interactions. This supports the hypothesis of "deep water nonintegrability". The experiments were performed in the framework of the new and refined version of the Zakharov equation free of nonessential terms in the quartic Hamiltonian. Simplification is possible due to exact cancellation of nonelastic four-wave interaction.

Description: Preface: Forecast and projection in climate scenario of Mediterranean intense events: uncertainties and propagation on environment (the MEDUP project) Natural Hazards and Earth System Science, 13, 3043-3047, 2013 Author(s): V. Ducrocq, P. Drobinski, D. Lambert, G. Molinié, and C. Llasat

Description: The role of risk perception in making flood risk management more effective Natural Hazards and Earth System Science, 13, 3013-3030, 2013 Author(s): M. Buchecker, G. Salvini, G. Di Baldassarre, E. Semenzin, E. Maidl, and A. Marcomini Over the last few decades, Europe has suffered from a number of severe flood events and, as a result, there has been a growing interest in probing alternative approaches to managing flood risk via prevention measures. A literature review reveals that, although in the last decades risk evaluation has been recognized as key element of risk management, and risk assessment methodologies (including risk analysis and evaluation) have been improved by including social, economic, cultural, historical and political conditions, the theoretical schemes are not yet applied in practice. One main reason for this shortcoming is that risk perception literature is mainly of universal and theoretical nature and cannot provide the necessary details to implement a comprehensive risk evaluation. This paper therefore aims to explore a procedure that allows the inclusion of stakeholders' perceptions of prevention measures in risk assessment. It proposes to adopt methods of risk communication (both one-way and two-way communication) in risk assessment with the final aim of making flood risk management more effective. The proposed procedure not only focuses on the effect of discursive risk communication on risk perception, and on achieving a shared assessment of the prevention alternatives, but also considers the effects of the communication process on perceived uncertainties, accepted risk levels, and trust in the managing institutions. The effectiveness of this combined procedure has been studied and illustrated using the example of the participatory flood prevention assessment process on the Sihl River in Zurich, Switzerland. The main findings of the case study suggest that the proposed procedure performed well, but that it needs some adaptations for it to be applicable in different contexts and to allow a (semi-) quantitative estimation of risk perception to be used as an indicator of adaptive capacity.

Description: QVAST: a new Quantum GIS plugin for estimating volcanic susceptibility Natural Hazards and Earth System Science, 13, 3031-3042, 2013 Author(s): S. Bartolini, A. Cappello, J. Martí, and C. Del Negro One of the most important tasks of modern volcanology is the construction of hazard maps simulating different eruptive scenarios that can be used in risk-based decision making in land-use planning and emergency management. The first step in the quantitative assessment of volcanic hazards is the development of susceptibility maps (i.e., the spatial probability of a future vent opening given the past eruptive activity of a volcano). This challenging issue is generally tackled using probabilistic methods that use the calculation of a kernel function at each data location to estimate probability density functions (PDFs). The smoothness and the modeling ability of the kernel function are controlled by the smoothing parameter, also known as the bandwidth. Here we present a new tool, QVAST, part of the open-source geographic information system Quantum GIS, which is designed to create user-friendly quantitative assessments of volcanic susceptibility. QVAST allows the selection of an appropriate method for evaluating the bandwidth for the kernel function on the basis of the input parameters and the shapefile geometry, and can also evaluate the PDF with the Gaussian kernel. When different input data sets are available for the area, the total susceptibility map is obtained by assigning different weights to each of the PDFs, which are then combined via a weighted summation and modeled in a non-homogeneous Poisson process. The potential of QVAST, developed in a free and user-friendly environment, is here shown through its application in the volcanic fields of Lanzarote (Canary Islands) and La Garrotxa (NE Spain).

Description: A joint probability approach using a 1-D hydrodynamic model for estimating high water level frequencies in the Lower Rhine Delta Natural Hazards and Earth System Science, 13, 1841-1852, 2013 Author(s): H. Zhong, P.-J. van Overloop, and P. H. A. J. M. van Gelder The Lower Rhine Delta, a transitional area between the River Rhine and Meuse and the North Sea, is at risk of flooding induced by infrequent events of a storm surge or upstream flooding, or by more infrequent events of a combination of both. A joint probability analysis of the astronomical tide, the wind induced storm surge, the Rhine flow and the Meuse flow at the boundaries is established in order to produce the joint probability distribution of potential flood events. Three individual joint probability distributions are established corresponding to three potential flooding causes: storm surges and normal Rhine discharges, normal sea levels and high Rhine discharges, and storm surges and high Rhine discharges. For each category, its corresponding joint probability distribution is applied, in order to stochastically simulate a large number of scenarios. These scenarios can be used as inputs to a deterministic 1-D hydrodynamic model in order to estimate the high water level frequency curves at the transitional locations. The results present the exceedance probability of the present design water level for the economically important cities of Rotterdam and Dordrecht. The calculated exceedance probability is evaluated and compared to the governmental norm. Moreover, the impact of climate change on the high water level frequency curves is quantified for the year 2050 in order to assist in decisions regarding the adaptation of the operational water management system and the flood defense system.

Description: Exploring local risk managers' use of flood hazard maps for risk communication purposes in Baden-Württemberg Natural Hazards and Earth System Science, 13, 1857-1872, 2013 Author(s): S. Kjellgren In response to the EU Floods Directive (2007/60/EC), flood hazard maps are currently produced all over Europe, reflecting a wider shift in focus from "flood protection" to "risk management", for which not only public authorities but also populations at risk are seen as responsible. By providing a visual image of the foreseen consequences of flooding, flood hazard maps can enhance people's knowledge about flood risk, making them more capable of an adequate response. Current literature, however, questions the maps' awareness raising capacity, arguing that their content and design are rarely adjusted to laypeople's needs. This paper wants to complement this perspective with a focus on risk communication by studying how these tools are disseminated and marketed to the public in the first place. Judging from communication theory, simply making hazard maps publicly available is unlikely to lead to attitudinal or behavioral effects, since this typically requires two-way communication and material or symbolic incentives. Consequently, it is relevant to investigate whether and how local risk managers, who are well positioned to interact with the local population, make use of flood hazard maps for risk communication purposes. A qualitative case study of this issue in the German state of Baden-Württemberg suggests that many municipalities lack a clear strategy for using this new information tool for hazard and risk communication. Four barriers in this regard are identified: perceived disinterest/sufficient awareness on behalf of the population at risk; unwillingness to cause worry or distress; lack of skills and resources; and insufficient support. These barriers are important to address – in research as well as in practice – since it is only if flood hazard maps are used to enhance local knowledge resources that they can be expected to contribute to social capacity building.

Description: Preface "2nd International Conference on Ecohydrology and Climate Change" Natural Hazards and Earth System Science, 13, 1853-1856, 2013 Author(s): C. Andrade, J. A. Santos, R. Trigo, and J. Corte-Real

Description: Coulomb static stress changes before and after the 23 October 2011 Van, eastern Turkey, earthquake ( M W = 7.1): implications for the earthquake hazard mitigation Natural Hazards and Earth System Science, 13, 1889-1902, 2013 Author(s): M. Utkucu, H. Durmuş, H. Yalçın, E. Budakoğlu, and E. Işık Coulomb stress changes before and after the 23 October 2011 Van, eastern Turkey, earthquake have been analysed using available data related to the background and the aftershock seismicity and the source faults. The coseismic stress changes of the background seismicity had slightly promoted stress over the rupture plane of the 2011 Van earthquake, while it yielded a stress shadow over the Gürpı nar Fault which has been argued to have produced the 7 April 1646 Van earthquake. The stress shadow over the Gürp\i nar fault has become more pronounced following the occurrence of the 2011 Van earthquake, meaning that the repetition of the 1646 Van earthquake has been further suppressed. Spatial distribution and source mechanisms of the 2011 Van earthquake's aftershocks have been utilised to define four clusters with regard to their relative location to the mainshock rupture. In addition, the aftershock sequence covers a much broader area toward the northeast. Correlations between the observed spatial patterns of the aftershocks and the coseismic Coulomb stress changes caused by the mainshock are determined by calculating the stress changes over both optimally oriented and specified fault planes. It is shown here that there is an apparent correlation between the mainshock stress changes and the observed spatial pattern of the aftershock occurrence, demonstrating the usefulness of the stress maps in constraining the likely locations of the upcoming aftershocks and mitigating earthquake hazard.

Description: A comparison of classical and intelligent methods to detect potential thermal anomalies before the 11 August 2012 Varzeghan, Iran, earthquake ( M w = 6.4) Natural Hazards and Earth System Science, 13, 1077-1083, 2013 Author(s): M. Akhoondzadeh In this paper, a number of classical and intelligent methods, including interquartile, autoregressive integrated moving average (ARIMA), artificial neural network (ANN) and support vector machine (SVM), have been proposed to quantify potential thermal anomalies around the time of the 11 August 2012 Varzeghan, Iran, earthquake ( M w = 6.4). The duration of the data set, which is comprised of Aqua-MODIS land surface temperature (LST) night-time snapshot images, is 62 days. In order to quantify variations of LST data obtained from satellite images, the air temperature (AT) data derived from the meteorological station close to the earthquake epicenter has been taken into account. For the models examined here, results indicate the following: (i) ARIMA models, which are the most widely used in the time series community for short-term forecasting, are quickly and easily implemented, and can efficiently act through linear solutions. (ii) A multilayer perceptron (MLP) feed-forward neural network can be a suitable non-parametric method to detect the anomalous changes of a non-linear time series such as variations of LST. (iii) Since SVMs are often used due to their many advantages for classification and regression tasks, it can be shown that, if the difference between the predicted value using the SVM method and the observed value exceeds the pre-defined threshold value, then the observed value could be regarded as an anomaly. (iv) ANN and SVM methods could be powerful tools in modeling complex phenomena such as earthquake precursor time series where we may not know what the underlying data generating process is. There is good agreement in the results obtained from the different methods for quantifying potential anomalies in a given LST time series. This paper indicates that the detection of the potential thermal anomalies derive credibility from the overall efficiencies and potentialities of the four integrated methods.

Description: Potential flood volume of Himalayan glacial lakes Natural Hazards and Earth System Science, 13, 1827-1839, 2013 Author(s): K. Fujita, A. Sakai, S. Takenaka, T. Nuimura, A. B. Surazakov, T. Sawagaki, and T. Yamanokuchi Glacial lakes are potentially dangerous sources of glacial lake outburst floods (GLOFs), and represent a serious natural hazard in Himalayan countries. Despite the development of various indices aimed at determining the outburst probability, an objective evaluation of the thousands of Himalayan glacial lakes has yet to be completed. In this study we propose a single index, based on the depression angle from the lakeshore, which allows the lakes to be assessed using remotely sensed digital elevation models (DEMs). We test our approach on five lakes in Nepal, Bhutan, and Tibet using images taken by the declassified Hexagon KH-9 satellite before these lakes experienced an outburst flood. All five lakes had a steep lakefront area (SLA), on which a depression angle was steeper than our proposed threshold of 10° before the GLOF event, but the SLA was no longer evident after the events. We further calculated the potential flood volume (PFV); i.e., the maximum volume of floodwater that could be released if the lake surface was lowered sufficiently to eradicate the SLA. This approach guarantees repeatability to assess the possibility of GLOF hazards because it requires no particular expertise to carry out, though the PFV does not quantify the GLOF risk. We calculated PFVs for more than 2000 Himalayan glacial lakes using visible band images and DEMs of ASTER data. The PFV distribution follows a power-law function. We found that 794 lakes did not have an SLA, and consequently had a PFV of zero, while we also identified 49 lakes with PFVs of over 10 million m 3 , which is a comparable volume to that of recorded major GLOFs. This PFV approach allows us to preliminarily identify and prioritize those Himalayan glacial lakes that require further detailed investigation on GLOF hazards and risk.

Description: Influence of flood risk characteristics on flood insurance demand: a comparison between Germany and the Netherlands Natural Hazards and Earth System Science, 13, 1691-1705, 2013 Author(s): I. Seifert, W. J. W. Botzen, H. Kreibich, and J. C. J. H. Aerts The existence of sufficient demand for insurance coverage against infrequent losses is important for the adequate function of insurance markets for natural disaster risks. This study investigates how characteristics of flood risk influence household flood insurance demand based on household surveys undertaken in Germany and the Netherlands. Our analyses confirm the hypothesis that willingness to pay (WTP) for insurance against medium-probability medium-impact flood risk in Germany is higher than WTP for insurance against low-probability high-impact flood risk in the Netherlands. These differences in WTP can be related to differences in flood experience, individual risk perceptions, and the charity hazard. In both countries there is a need to stimulate flood insurance demand if a relevant role of private insurance in flood loss compensation is regarded as desirable, for example, by making flood insurance compulsory or by designing information campaigns.

Description: Synoptic climatology of winter intense precipitation events along the Mediterranean coasts Natural Hazards and Earth System Science, 13, 1707-1722, 2013 Author(s): M. Reale and P. Lionello The link between winter (December-January-February) precipitation events at 15 Mediterranean coastal locations and synoptic features (cyclones and Northern Hemisphere teleconnection patterns) is analyzed. A list of precipitation events has been produced; q percentile thresholds (Th q ) and corresponding frequency N q (for q equal to 25, 50, 90 and 98) have been considered. A negative trend has been detected in total precipitation and N 50 at many locations, while no significant trend in N 25 , N 90 and N 98 has been found. The negative phase of the North Atlantic Oscillation (NAO) and the East Atlantic/West Russia pattern (EAWR) compete for exerting the largest influence on the frequency of the 25th, 50th and 90th percentiles, with EAWR and NAO exerting their largest influence in the central and western Mediterranean areas, respectively. All percentiles show a similar behavior except for the 98th percentile, which shows no convincing link to any teleconnection pattern. The cyclone tracks that are associated with precipitation events have been selected using the ERA-40 reanalysis data, and a strong link between intense precipitation and cyclones is shown for all stations. In general, the probability of detecting a cyclone within a distance of 20° from each station increases with the intensity of the precipitation event and decreases with the duration of a dry period. The origin and track of cyclones producing intense precipitation differ among different areas. When precipitation occurs in the northwestern Mediterranean, cyclones are generally either of Atlantic origin or secondary cyclones associated with the passage of major cyclones north of the Mediterranean Basin, while they are mostly generated inside the region itself for events at the eastern Mediterranean coast. An important fraction of intense events in the southern areas is produced by cyclones that are generated over northern Africa. The analysis of sea level pressure and geopotential height at 500 hPa highlights the important role of cyclone depth, circulation strength, surrounding synoptic condition, and of slow speed of the cyclone center for producing intense precipitation events.

Description: A relative vulnerability estimation of flood disaster using data envelopment analysis in the Dongting Lake region of Hunan Natural Hazards and Earth System Science, 13, 1723-1734, 2013 Author(s): C.-H. Li, N. Li, L.-C. Wu, and A.-J. Hu The vulnerability to flood disaster is addressed by a number of studies. It is of great importance to analyze the vulnerability of different regions and various periods to enable the government to make policies for distributing relief funds and help the regions to improve their capabilities against disasters, yet a recognized paradigm for such studies seems missing. Vulnerability is defined and evaluated through either physical or economic–ecological perspectives depending on the field of the researcher concerned. The vulnerability, however, is the core of both systems as it entails systematic descriptions of flood severities or disaster management units. The research mentioned often has a development perspective, and in this article we decompose the overall flood system into several factors: disaster driver, disaster environment, disaster bearer, and disaster intensity, and take the interaction mechanism among all factors as an indispensable function. The conditions of flood disaster components are demonstrated with disaster driver risk level, disaster environment stability level and disaster bearer sensitivity, respectively. The flood system vulnerability is expressed as vulnerability = f (risk, stability, sensitivity). Based on the theory, data envelopment analysis method (DEA) is used to detail the relative vulnerability's spatiotemporal variation of a flood disaster system and its components in the Dongting Lake region. The study finds that although a flood disaster system's relative vulnerability is closely associated with its components' conditions, the flood system and its components have a different vulnerability level. The overall vulnerability is not the aggregation of its components' vulnerability. On a spatial scale, zones central and adjacent to Dongting Lake and/or river zones are characterized with very high vulnerability. Zones with low and very low vulnerability are mainly distributed in the periphery of the Dongting Lake region. On a temporal scale, the occurrence of a vibrating flood vulnerability trend is observed. A different picture is displayed with the disaster driver risk level, disaster environment stability level and disaster bearer sensitivity level. The flood relative vulnerability estimation method based on DEA is characteristic of good comparability, which takes the relative efficiency of disaster system input–output into account, and portrays a very diverse but consistent picture with varying time steps. Therefore, among different spatial and time domains, we could compare the disaster situations with what was reflected by the same disaster. Additionally, the method overcomes the subjectivity of a comprehensive flood index caused by using an a priori weighting system, which exists in disaster vulnerability estimation of current disasters.

Description: Complementary methods to plan pedestrian evacuation of the French Riviera's beaches in case of tsunami threat: graph- and multi-agent-based modelling Natural Hazards and Earth System Science, 13, 1735-1743, 2013 Author(s): A. Sahal, F. Leone, and M. Péroche Small amplitude tsunamis have impacted the French Mediterranean shore (French Riviera) in the past centuries. Some caused casualties; others only generated economic losses. While the North Atlantic and Mediterranean tsunami warning system is being tested and is almost operational, no awareness and preparedness measure is being implemented at a local scale. Evacuation is to be considered along the French Riviera, but no plan exists within communities. We show that various approaches can provide local stakeholders with evacuation capacities assessments to develop adapted evacuation plans through the case study of the Cannes–Antibes region. The complementarity between large- and small-scale approaches is demonstrated with the use of macro-simulators (graph-based) and micro-simulators (multi-agent-based) to select shelter points and choose evacuation routes for pedestrians located on the beach. The first one allows automatically selecting shelter points and measuring and mapping their accessibility. The second one shows potential congestion issues during pedestrian evacuations, and provides leads for the improvement of urban environment. Temporal accessibility to shelters is compared to potential local and distal tsunami travel times, showing a 40 min deficit for an adequate crisis management in the first scenario, and a 30 min surplus for the second one.

Description: Forcing factors of cloud-to-ground lightning over Iberia: regional-scale assessments Natural Hazards and Earth System Science, 13, 1745-1758, 2013 Author(s): J. A. Santos, M. A. Reis, F. De Pablo, L. Rivas-Soriano, and S. M. Leite Cloud-to-ground lightning in a sector covering the Iberian Peninsula, the Balearic Islands and nearby seas (36–44° N, 10° W–5° E) is analysed in the period from 2003 to 2009 (7 yr). Two Iberian lightning detection networks, composed of 18 sensors over Portugal and Spain, are combined for the first time in the present study. The selected characteristics are cloud-to-ground flashes (CGFs), first stroke peak current, polarity and multiplicity (number of strokes in a given flash). This study examines the temporal (on hourly, monthly and seasonal timescales) and spatial variability of CGFs. The influence of five forcing factors on lightning (elevation, lifted index, convective available potential energy and daily minimum and maximum near-surface air temperatures) over the Iberian sector is also assessed. For regional-scale assessments, six subsectors with different climatic conditions were analysed separately. Despite important regional differences, the strongest lightning activity occurs from late spring to early autumn, and mostly in the afternoon. Furthermore, CGFs are mainly located over high-elevation areas in late spring to summer, while they tend to occur over the sea in autumn. The results suggest that (1) orographically forced thunderstorms over mountainous areas, mostly from May to September, (2) tropospheric buoyancy forcing over western-central and northern regions in summer and over the Mediterranean regions in autumn, and (3) near-surface thermal contrasts from October to February largely control the location of lightning in Iberia. There is no evidence of different forcings by polarity. A clear correspondence between summertime precipitation patterns and CGFs is also found.

Description: Rogue waves, rogue events and extreme wave kinematics in spatio-temporal fields of simulated sea states Natural Hazards and Earth System Science, 13, 1759-1771, 2013 Author(s): A. Sergeeva and A. Slunyaev An approach to the extensive study of rogue wave occurrence in numerical simulations is presented. As a result of numerical simulations of the unidirectional wave evolution, spatio-temporal fields of wave data of the size 20 min × 10 km are obtained with high resolution in time and space and are used for statistical analysis with the focus on extreme waves. Having the exhaustive information on the wave evolution enables us to capture the detailed picture of individual rogue waves; to detect intermittent rogue wave events, which last for a significantly longer time, and hence, to depict the portrait of a rogue wave. Due to the benefit of having full-wave data, the question of relation between extreme wave kinematics and extremely high waves is discussed in the statistical sense.

Description: The UBO-TSUFD tsunami inundation model: validation and application to a tsunami case study focused on the city of Catania, Italy Natural Hazards and Earth System Science, 13, 1795-1816, 2013 Author(s): S. Tinti and R. Tonini Nowadays numerical models are a powerful tool in tsunami research since they can be used (i) to reconstruct modern and historical events, (ii) to cast new light on tsunami sources by inverting tsunami data and observations, (iii) to build scenarios in the frame of tsunami mitigation plans, and (iv) to produce forecasts of tsunami impact and inundation in systems of early warning. In parallel with the general recognition of the importance of numerical tsunami simulations, the demand has grown for reliable tsunami codes, validated through tests agreed upon by the tsunami community. This paper presents the tsunami code UBO-TSUFD that has been developed at the University of Bologna, Italy, and that solves the non-linear shallow water (NSW) equations in a Cartesian frame, with inclusion of bottom friction and exclusion of the Coriolis force, by means of a leapfrog (LF) finite-difference scheme on a staggered grid and that accounts for moving boundaries to compute sea inundation and withdrawal at the coast. Results of UBO-TSUFD applied to four classical benchmark problems are shown: two benchmarks are based on analytical solutions, one on a plane wave propagating on a flat channel with a constant slope beach; and one on a laboratory experiment. The code is proven to perform very satisfactorily since it reproduces quite well the benchmark theoretical and experimental data. Further, the code is applied to a realistic tsunami case: a scenario of a tsunami threatening the coasts of eastern Sicily, Italy, is defined and discussed based on the historical tsunami of 11 January 1693, i.e. one of the most severe events in the Italian history.

Description: Performance of coastal sea-defense infrastructure at El Jadida (Morocco) against tsunami threat: lessons learned from the Japanese 11 March 2011 tsunami Natural Hazards and Earth System Science, 13, 1779-1794, 2013 Author(s): R. Omira, M. A. Baptista, F. Leone, L. Matias, S. Mellas, B. Zourarah, J. M. Miranda, F. Carrilho, and J.-P. Cherel This paper seeks to investigate the effectiveness of sea-defense structures in preventing/reducing the tsunami overtopping as well as evaluating the resulting tsunami impact at El Jadida, Morocco. Different tsunami wave conditions are generated by considering various earthquake scenarios of magnitudes ranging from M w = 8.0 to M w = 8.6. These scenarios represent the main active earthquake faults in the SW Iberia margin and are consistent with two past events that generated tsunamis along the Atlantic coast of Morocco. The behaviour of incident tsunami waves when interacting with coastal infrastructures is analysed on the basis of numerical simulations of near-shore tsunami waves' propagation. Tsunami impact at the affected site is assessed through computing inundation and current velocity using a high-resolution digital terrain model that incorporates bathymetric, topographic and coastal structures data. Results, in terms of near-shore tsunami propagation snapshots, waves' interaction with coastal barriers, and spatial distributions of flow depths and speeds, are presented and discussed in light of what was observed during the 2011 Tohoku-oki tsunami. Predicted results show different levels of impact that different tsunami wave conditions could generate in the region. Existing coastal barriers around the El Jadida harbour succeeded in reflecting relatively small waves generated by some scenarios, but failed in preventing the overtopping caused by waves from others. Considering the scenario highly impacting the El Jadida coast, significant inundations are computed at the sandy beach and unprotected areas. The modelled dramatic tsunami impact in the region shows the need for additional tsunami standards not only for sea-defense structures but also for the coastal dwellings and houses to provide potential in-place evacuation.

Description: Implications of spatial distribution of rockfall reconstructed by dendrogeomorphological methods Natural Hazards and Earth System Science, 13, 1817-1826, 2013 Author(s): K. Šilhán, T. Pánek, and J. Hradecký Rockfall is a dangerous geomorphological process. The prediction of potentially threatened areas requires thorough reconstruction of spatial rockfall activity. Dendrogeomorphic methods allow precise determination of both temporal and spatial occurrences of rockfall without the necessity of long-term monitoring. At the case-study site of Taraktash, located among southern slopes of the Crimean Mountains, 114 Crimean pine trees ( Pinus nigra ssp. pallasiana ) were sampled on a talus slope located under a 150 m high rockwall. Based on their age, the trees were divided into two distinct groups (young and old trees). Considerable disturbance in the age structure of the trees on the talus was probably caused by a series of strong earthquakes. Major differences were identified in the ability of young and old trees to record a rockfall event. We found that in the first decades of their growth, the ability of the studied P. nigra to record rockfall events gradually increased. The trees showed the highest sensitivity at the age of 80 to 90 yr; after that age their sensitivity gradually decreases. Two indicators were selected for the spatial reconstruction of rockfall events (the number of rockfall events per tree and recurrence interval). The highest activity was identified on the talus using selected indicators.

Description: Brief Communication: Likelihood of societal preparedness for global change: trend detection Natural Hazards and Earth System Science, 13, 1773-1778, 2013 Author(s): R. M. Vogel, A. Rosner, and P. H. Kirshen Anthropogenic influences on earth system processes are now pervasive, resulting in trends in river discharge, pollution levels, ocean levels, precipitation, temperature, wind, landslides, bird and plant populations and a myriad of other important natural hazards relating to earth system state variables. Thousands of trend detection studies have been published which report the statistical significance of observed trends. Unfortunately, such studies only concentrate on the null hypothesis of "no trend". Little or no attention is given to the power of such statistical trend tests, which would quantify the likelihood that we might ignore a trend if it really existed. The probability of missing the trend, if it exists, known as the type II error, informs us about the likelihood of whether or not society is prepared to accommodate and respond to such trends. We describe how the power or probability of detecting a trend if it exists, depends critically on our ability to develop improved multivariate deterministic and statistical methods for predicting future trends in earth system processes. Several other research and policy implications for improving our understanding of trend detection and our societal response to those trends are discussed.

Description: Assessing gridded observations for daily precipitation extremes in the Alps with a focus on northwest Italy Natural Hazards and Earth System Science, 13, 1457-1468, 2013 Author(s): M. Turco, A. L. Zollo, C. Ronchi, C. De Luigi, and P. Mercogliano In this study we compare three gridded observed datasets of daily precipitation (EOBS, MAP and NWIOI) over the Great Alpine Region (GAR) and a subregion in northwest Italy (NWI) in order to better understand the past variability of daily climate extremes and to set up a basis for developing regional climate scenarios. The grids are first compared with respect to their temporal similarity by calculating the correlation and relative mean absolute error of the time series. They are then compared with respect to their spatial similarity to the climatology of the ETCCDI indices (characterizing total precipitation, dry and wet spells and extremes with short return periods). The results indicate first that most EOBS gridpoint series in northeastern Italy have to be shifted back by 1 day to have maximum overlap of the measurement period and, second, that both the temporal and spatial similarities of most indices are higher between the NWIOI and MAP than between MAP or the NWIOI and EOBS. These results suggest that, although there is generally good temporal agreement between the three datasets, EOBS should be treated with caution, especially for extreme indices over the GAR region, and it does not provide reliable climatology over the NWI region. The high agreement between MAP and NWIOI, on the other hand, builds confidence in using these datasets. Users should consider carefully the limitations of the gridded observations available: the uncertainties of the observed datasets cannot be neglected in the overall uncertainties cascade that characterizes climate change studies.

Description: The application of FORMOSAT-2 high-temporal- and high-spatial resolution imagery for monitoring open straw burning and carbon emission detection Natural Hazards and Earth System Science, 13, 575-582, 2013 Author(s): C.-C. Liu, P.-Y. Tseng, and C.-Y. Chen Rice is produced in more than 95 countries worldwide and is a staple food for over half of the world's population. Rice is also a major food crop of Taiwan. There are numerous rice crops planted on the western plains of Taiwan, and, after the harvest season, the left-over straw is often burned on-site. The air pollutants from the burning emissions include CO 2 , CO, CH 4 and other suspended particles, most of these being the greenhouse gases which cause global climate change. In this study FORMOSAT-2 satellite images and ground-truth data from 2008 and 2009 are used to conduct supervised classification and calculate the extent of the straw burning areas. It was found that 10% of the paddies in the study area were burned after harvest during this 2-yr period. On this pro rata basis, we calculated the overall carbon emissions from the burning of the straw. The findings showed that these few farmers produced up to 34 000 tons of carbon emissions in 2008, and 40 000 tons in 2009. The study results indicate that remotely sensed images can be used to efficiently evaluate the important characteristics for carbon emission detection. It also provides quantitative results that are relevant to tracking sources of transport pollution, postharvest burning, and Asian dust in Taiwan.