Evaluation of approaches to calculate debris-flow parameters for hazard assessment

Many different runout prediction methods can be applied to estimate the mobility of future debris flows during hazard assessment. The present article reviews the empirical, analytical, simple flow routing and numerical techniques. All these techniques were applied to back-calculate a debris flow, which occurred in 1982 at La Guingueta catchment, in the Eastern Pyrenees. A sensitivity analysis of input parameters was carried out, while special attention was paid to the influence of rheological parameters. We used the Voellmy fluid rheology for our analytical and numerical modelling, since this flow resistance law coincided best with field observations. The simulation results indicated that the “basal” friction coefficients rather affect the runout distance, while the “turbulence” terms mainly influence flow velocity. A comparison of the velocity computed on the fan showed that the analytical model calculated values similar to the numerical ones. The values of our rheological parameters calibrated at La Guingueta agree with data back-calculated for other debris flows. Empirical relationships represent another method to estimate total runout distance. The results confirmed that they contain an important uncertainty and they are strictly valid only for the conditions, which were the basis for their development. With regards to the simple flow routing algorithm, this methods could satisfactorily simulate the total area affected by the 1982 debris flow, but it was not able to directly calculate total runout distance and velocity. Finally, a suggestion on how different runout prediction methods can be applied to generate debris-flow hazard maps is presented. Taking into account the definition of hazard and intensity, the best choice would be to divide the resulting hazard maps into two types: “final hazard maps” and “preliminary hazard maps”. Only the use of numerical models provided final hazard maps, because they could incorporate different event magnitudes and they supplied output-values for intensity calculation. In contrast, empirical relationships and flow routing algorithms, or a combination of both, could be applied to create preliminary hazard maps. The present study only focussed on runout prediction methods. Other necessary tasks to complete the hazard assessment can be looked up in the “Guidelines for landslide susceptibility, hazard and risk zoning” included in this Special Issue.     

Seismic Hazard Analysis — Quo vadis?

The paper is dedicated to the review of methods of seismic hazard analysis currently in use, analyzing the strengths and weaknesses of different approaches. The review is performed from the perspective of a user of the results of seismic hazard analysis for different applications such as the design of critical and general (non-critical) civil infrastructures, technical and financial risk analysis. A set of criteria is developed for and applied to an objective assessment of the capabilities of different analysis methods. It is demonstrated that traditional probabilistic seismic hazard analysis (PSHA) methods have significant deficiencies, thus limiting their practical applications. These deficiencies have their roots in the use of inadequate probabilistic models and insufficient understanding of modern concepts of risk analysis, as have been revealed in some recent large scale studies. These deficiencies result in the lack of ability of a correct treatment of dependencies between physical parameters and finally, in an incorrect treatment of uncertainties. As a consequence, results of PSHA studies have been found to be unrealistic in comparison with empirical information from the real world. The attempt to compensate these problems by a systematic use of expert elicitation has, so far, not resulted in any improvement of the situation. It is also shown that scenario-earthquakes developed by disaggregation from the results of a traditional PSHA may not be conservative with respect to energy conservation and should not be used for the design of critical infrastructures without validation. Because the assessment of technical as well as of financial risks associated with potential damages of earthquakes need a risk analysis, current method is based on a probabilistic approach with its unsolved deficiencies.

Traditional deterministic or scenario-based seismic hazard analysis methods provide a reliable and in general robust design basis for applications such as the design of critical infrastructures, especially with systematic sensitivity analyses based on validated phenomenological models. Deterministic seismic hazard analysis incorporates uncertainties in the safety factors. These factors are derived from experience as well as from expert judgment. Deterministic methods associated with high safety factors may lead to too conservative results, especially if applied for generally short-lived civil structures. Scenarios used in deterministic seismic hazard analysis have a clear physical basis. They are related to seismic sources discovered by geological, geomorphologic, geodetic and seismological investigations or derived from historical references. Scenario-based methods can be expanded for risk analysis applications with an extended data analysis providing the frequency of seismic events. Such an extension provides a better informed risk model that is suitable for risk-informed decision making.     

An operational multiscale system for hazards prediction, mapping, and response

By definition, a crisis is a situation that requires assistance to be managed. Hence, response to a crisis involves the merging of local and non-local emergency response personnel. In this situation, it is critical that each participant: (1) know the roles and responsibilities of each of the other participants; (2) know the capabilities of each of the participants; and (3) have a common basis for action. For many types of natural disasters, this entails having a common operational picture of the unfolding events, including detailed information on the weather, both current and forecasted, that may impact on either the emergency itself or on response activities. The Consequences Assessment Tool Set (CATS) is a comprehensive package of hazard prediction models and casualty and damage assessment tools that provides a linkage between a modeled or observed effect and the attendant consequences for populations, infrastructure, and resources, and, hence, provides the common operational picture for emergency response. The Operational Multiscale Environment model with Grid Adaptivity (OMEGA) is an atmospheric simulation system that links the latest methods in computational fluid dynamics and high-resolution gridding technologies with numerical weather prediction to provide specific weather analysis and forecast capability that can be merged into the geographic information system framework of CATS. This paper documents the problem of emergency response as an end-to-end system and presents the integrated CATS–OMEGA system as a prototype of such a system that has been used successfully in a number of different situations.     

Quantitative sinkhole hazard assessment. A case study from the Ebro Valley evaporite alluvial karst (NE Spain)

Quantitative sinkhole hazard assessments in karst areas allow calculation of the potential sinkhole risk and the performance of cost-benefit analyses. These estimations are of practical interest for planning, engineering, and insurance purposes. The sinkhole hazard assessments should include two components: the probability of occurrence of sinkholes (sinkholes/km² year) and the severity of the sinkholes, which mainly refers to the subsidence mechanisms (progressive passive bending or catastrophic collapse) and the size of the sinkholes at the time of formation; a critical engineering design parameter. This requires the compilation of an exhaustive database on recent sinkholes, including information on the: (1) location, (2) chronology (precise date or age range), (3) size, and (4) subsidence mechanisms and rate. This work presents a hazard assessment from an alluvial evaporite karst area (0.81 km²) located in the periphery of the city of Zaragoza (Ebro River valley, NE Spain). Five sinkholes and four locations with features attributable to karstic subsidence where identified in an initial investigation phase providing a preliminary probability of occurrence of 0.14 sinkholes/km² year (11.34% in annual probability). A trenching program conducted in a subsequent investigation phase allowed us to rule out the four probable sinkholes, reducing the probability of occurrence to 0.079 sinkholes/km² year (6.4% in annual probability). The information on the severity indicates that collapse sinkholes 10–15 m in diameter may occur in the area. A detailed study of the deposits and deformational structures exposed by trenching in one of the sinkholes allowed us to infer a modern collapse sinkhole approximately 12 m in diameter and with a vertical throw of 8 m. This collapse structure is superimposed on a subsidence sinkhole around 80 m across that records at least 1.7 m of synsedimentary subsidence. Trenching, in combination with dating techniques, is proposed as a useful methodology to elucidate the origin of depressions with uncertain diagnosis and to gather practical information with predictive utility about particular sinkholes in alluvial karst settings: precise location, subsidence mechanisms and magnitude, and timing and rate of the subsidence episodes.     

Impact of lightning on organic matter-rich soils: influence of soil grain size and organic matter content on underground fires

A detailed study was carried out on a piece of land that had been struck by lightning during the violent rainstorm that raged over the Island of S?o Miguel (Azores Archipelago) in late October 2006. Temperature and gas measurements (CO₂, CO, H₂S and CH₄) were performed in four study trenches, dug in an area of ∼3 m², where an underground fire had been initiated by the impact with a lightning stroke, followed by the emission of a column of gases and smoke. The soil under study was originally a well-pedogenized about 80 cm thick bed, made of volcanic clayey to silty tephra fallouts and contained 5.5–9.7% of organic matter. The underground fire was monitored for one week and revealed a peak release of 404 ppm CO and 3.4% CO₂ originating from a horizon located about 45 cm under the soil surface. Measurements of temperature, performed one week after the impact, indicated a maximum value of 326°C inside the soil, while 516.5°C were measured on the surface of a lava block interred about 20 cm under the surface. Subsequently, a stratigraphic and sedimentologic study proved the role of the grain-size of the soil and of the organic matter content of the different horizons of the impact area, in determining the ratio between anoxic/oxidised combustion conditions and in the progress of the process itself. It was also noticed that combustion was not total all over in the soil bed and that the process had slightly migrated toward SW during the observation period. The combustion process went on for about ten days, in spite of several other violent rainstorms, until it was artificially extinguished through the excavations made to obtain study trenches. This particular circumstance evidenced the potential natural hazard represented by this kind of atmospheric event, especially in a land where the volcanic nature of the soil may easily mislead inexperienced observers and, consequently, delay proper action.     

Methodology for preliminary assessment of Natech risk in urban areas

Concern for natural hazard-triggered technological disasters (Natech disasters) in densely populated and industrialized areas is growing. Residents living in urban areas subject to high natural hazard risk are often unaware of the potential for secondary disasters such as hazardous materials releases from neighboring industrial facilities, chemical storage warehouses or other establishments housing hazardous materials. Lessons from previous disasters, such as the Natech disaster during the Kocaeli earthquake in Turkey in 1999 call for the need to manage low frequency/high consequence events, particularly in today’s densely populated areas. However, there is little guidance available on how local governments and communities can assess Natech risk. To add to the problem, local governments often do not have the human or economic resources or expertise to carry out detailed risk assessments. In this article, we propose a methodology for preliminary assessment of Natech risk in urban areas. The proposed methodology is intended for use by local government officials in consultation with the public. The methodology considers possible interactions between the various systems in the urban environment: the physical infrastructure (e.g., industrial plants, lifeline systems, critical facilities), the community (e.g., population exposed), the natural environment (e.g., delicate ecosystems, river basins), and the risk and emergency management systems (e.g., structural and nonstructural measures). Factors related to vulnerability and hazard are analyzed and qualitative measures are recommended. Data from hazardous materials releases during the Kocaeli, Turkey earthquake of August 17, 1999 are used as a case study to demonstrate the applicability of the methodology. Limitations of the proposed methodology are discussed as well as future research needs.

Pervasive wetland flooding in the glacial drift prairie of North Dakota (USA)

This article describes a unique flood hazard, produced by the dramatic expansion of wetlands in Nelson County, located within the North American Prairie Pothole Region of North Dakota, USA. There has been an unprecedented increase in the number, average size, and permanence of prairie wetlands, and a significant increase in the size of a closed lake (Stump Lake) due to a decade-long wet spell that began in 1993 following a prolonged drying trend. Base-line land cover information from the 1992 USGS National Land Cover Characterization dataset, and a Landsat TM scene acquired 9 July 2001 are used to assess the growth of the closed lake and wetland pond surface areas, and to analyze the type and area of various land cover classes inundated between 1992 and 2001. The open water profile in Nelson County changed from one marked by relatively comparable coverage of closed lake and wetland pond areas in 1992, to one in which wetland open water accounted for the vast majority of total open water in 2001. The bulk of the wetland pond area expansion occurred by displacing existing wetland vegetation and agricultural cropland. Producers responded to the flood hazard by filing Federal Crop Insurance Corporation (FCIC) claims and enrolling cropland in the Conservation Reserve Program (CRP), a federal land retirement program. Land taken out of agricultural production has had an enormous impact upon the agricultural sector that forms the economic base of the rural economy. In 2001 the land taken out of production due to CRP enrollment and preventive planting claims represented nearly 42% of Nelson County’s 205.2 K ha base agricultural land. The patterns obtained from this detailed study of Nelson County are likely to be the representative of the more publicized flood disaster occurring within the Devils Lake Basin of North Dakota.     

乌鲁木齐河流域地下水水质监测网设计

文章运用地下水易污性编图及污染源分布图法进行了乌鲁木齐河流域地下水水质监测网设计。共设计了130监测孔,现有46个监测孔,另需要84个新的监测孔。按监测类型分为面源监测点22个,点源监测点87个,重点水源地与泉水监测点21个。按监测运行分长期监测点55个,流域普查监测点75个。普查监测点监测频率为1次/5年,长期监测点监测频率为1次/年。首期有针对性地在污染严重的柴窝堡新化厂排污区、乌鲁木齐河谷老排污区、米泉污灌区、米泉工业污染区、老龙河污染区取了25个污染水样测试分析,结果显示地下水已经受到严重污染。     

加强技术创新服务地质调查——庆祝探矿工艺研究所成立30周年

为加强地质探矿技术力量,改变地质装备落后现状,经国务院批准,1978年8月8日成立了探矿工艺研究所。30年来,经过新老两代工艺所人的艰苦奋斗、刻苦钻研、勇于创新,在科研攻关、成果转化、基地建设和管理等方面都取得了可喜的成绩,为国家的地质科学事业做出了贡献,目前已成为一个“以勘查技术为基础,以承担地质调查任务和科研项目、提供勘查与监测技术方法和技术服务为手段,以服务国家经济建设、社会发展和地质工作为目标,以探矿工艺技术和地质灾害监测防治技术研发为主业”的地质调查科研机构。回顾了工艺所艰苦创业、改革进取的发展历程,对科技体制改革、专业结构调整、人事分配制度改革、地质调查和科研工作进行了总结,展示了丰硕的科研成果和强劲的科技实力;指出了工艺所在“十一五”期间地质调查和科研工作的重点领域和优先发展方向。     

杭州泗岭铝质A型花岗岩的发现及其构造意义

初步研究表明,以往被认为是典型的S型花岗岩或改造-重熔型花岗岩的泗岭岩体应为铝质A型花岗岩。该岩体高硅、碱、钾而明显贫钙、镁;高(K2O Na2O)/CaO值(平均16.79)和AKI指数(平均0.92);FeO*/MgO比值大(平均13.95),高于M型、S型和I型花岗岩,而与世界A型花岗岩平均值(13.4)相近;岩石为弱碱性、准铝-弱过铝质,富含稀土元素、HFSE元素,高Rb、F和Nb,反映其具A型花岗岩的成分特征。泗岭岩体氧同位素δ18O为 8‰～ 9‰,反映岩浆主要起源于下地壳。该岩体侵位于拉张的构造环境,是晚白垩世早期一次重要的构造-岩浆热事件,同时也反映了区域性孝丰—三门湾大断裂与学川—湖州大断裂在晚白垩世早期处于拉伸阶段。