A New Method for the Estimation of Avalanche Distance Exceeded Probabilities

A crucial point in any methodology for avalanche hazard assessment is the evaluation of avalanche distance exceeded probability, i.e., the annual probability that any assigned location along a given path is reached or exceeded by an avalanche. Typically this problem is faced by estimating the snow volume in the starting zone that is likely to accumulate an average every T years by statistical analysis of snowfall record, and then using this volume as input to an appropriately calibrated avalanche dynamics model to determine the runout distancesfor this design event. This methodology identifies the areas that canbe affected by an avalanche for the considered value of the return period (i.e. the average interval of time for a certain event to repeat itself), ¯T. However, it does not allow us to evaluate the actual avalanche encounter probability for any given point in the runout zone. In the present work this probability is computed by numerical integration of the expression P(x) = ∫₀ ^∞ P^*(V)f(V) dV, where f is the probabilitydensity function (PDF) of the avalanche release volume V, and P^* is the probability of the point x being reached or passed by an avalanche if the release volume is V; this latter probability is calculated by avalanche dynamics simulations. The procedure is implemented using a one-dimensional hydraulic-continuum avalanche dynamic model, calibrated on data from different Italian Alpine ranges, and is applied to a real world hazard mapping problem.     

Local And Global Diffusion Along Resonant Lines in Discrete Quasi-integrable Dynamical Systems

We detect and measure diffusion along resonances in a quasi-integrable symplectic map for different values of the perturbation parameter. As in a previously studied Hamiltonian case (Lega et al., 2003) results agree with the prediction of the Nekhoroshev theorem. Moreover, for values of the perturbation parameter slightly below the critical value of the transition between Nekhoroshev and Chirikov regime we have also found a diffusion of some orbits along macroscopic portions of the phase space. Such a diffusion follows in a spectacular way the peculiar structure of resonant lines.     

人为因素导致的地面沉降

天然气、石油、地下水、地热热水和卤水等流体运移出地下储存地层时,地层耗尽了流体而产生压缩变形,这些变形传递到地表表面就形成了人为的地面沉降。在本文中:(1)列举了世界上主要的地面沉降区域;(2)论述了引发地面沉降的力学机制,这些可以量测的地面沉降发生在地下含水层体系、天然气气田和石油油田之上;(3)描述了目前可以用于地面沉降测量和岩石变形原位测试的技术手段;(4)简要介绍了几个地面沉降预测的数值模型;(5)说明了用于控制地面沉降发展和减轻相关环境影响的几种防治措施。     

The Holocene Secche di Lazzaro phreatomagmatic succession (Stromboli,Italy): evidence of pyroclastic density current origin deduced by facies analysis and AMS flow directions

The edifice of Stromboli volcano gravitationally collapsed several times during its volcanic history (>100 ka–present). The largest Holocene event occurred during the final stage of the Neostromboli activity (∼13–5 ka), and was accompanied by the emplacement of phreatomagmatic and lahar deposits, known as the Secche di Lazzaro succession. A stratigraphic and paleomagnetic study of the Secche di Lazzaro deposits allows the interpretation of the emplacement and the eruptive processes. We identify three main units within the succession that correspond to changing eruption conditions. The lower unit (UA) consists of accretionary lapilli-rich, thinly bedded, parallel- to cross-stratified ash deposits, interpreted to indicate the early stages of the eruption and emplacement of dilute pyroclastic density currents. Upward, the second unit (UB) of the deposit is more massive and the beds thicker, indicating an increase in the sedimentation rate from pyroclastic density currents. The upper unit (UC) caps the succession with thick, immediately post-eruptive lahars, which reworked ash deposited on the volcano’s slope. Flow directions obtained by Anisotropy of Magnetic Susceptibility (AMS) analysis of the basal bed of UA at the type locality suggest a provenance of pyroclastic currents from the sea. This is interpreted to be related to the initial base-surges associated with water–magma interaction that occurred immediately after the lateral collapse, which wrapped around the shoulder of the sector collapse scar. Upward in the stratigraphy (upper beds of UA and UB) paleoflow directions change and show a provenance from the summit vent, probably related to the multiple collapses of a vertical, pulsatory eruptive column.     

Expected loss‐based alarm threshold set for earthquake early warning systems

Earthquake early warning systems (EEWS) seem to have potential as tools for real‐time seismic risk management and mitigation. In fact, although the evacuation of buildings requires warning time not available in many urbanized areas threatened by seismic hazard, they may still be used for the real‐time protection of critical facilities using automatic systems in order to reduce the losses subsequent to a catastrophic event. This is possible due to the real‐time seismology, which consists of methods and procedures for the rapid estimation of earthquake features, as magnitude and location, based on measurements made on the first seconds of the P‐waves. An earthquake engineering application of earthquake early warning (EEW) may be intended as a system able to issue the alarm, if some recorded parameter exceeds a given threshold, to activate risk mitigation actions before the quake strikes at a site of interest. Feasibility analysis and design of such EEWS require the assessment of the expected loss reduction due to the security action and set of the alarm threshold. In this paper a procedure to carry out these tasks in the performance‐based earthquake engineering probabilistic framework is proposed. A merely illustrative example refers to a simple structure assumed to be a classroom. Structural damage and non‐structural collapses are considered; the security action is to shelter occupants below the desks. The cost due to a false alarm is assumed to be related to the interruption of didactic activities. Results show how the comparison of the expected losses, for the alarm‐issuance and non‐issuance cases, allows setting the alarm threshold on a quantitative and consistent basis, and how it may be a tool for the design of engineering applications of EEW. Copyright © 2007 John Wiley & Sons, Ltd.     

Lava flow hazard at Nyiragongo Volcano,DRC

Mt. Nyiragongo is one of the most dangerous volcanoes in the world for the risk associated with the propagation of lava flows. In 2002 several vents opened along a huge system of fractures, pouring out lava which reached and destroyed a considerable part of Goma, a town of about 500,000 inhabitants on the shore of Lake Kivu. In a companion paper (Favalli et al. in Bull Volcanol, this issue, 2008) we employed numerical simulations of probable lava flow paths to evaluate the lava flow hazard on the flanks of the volcano, including the neighbouring towns of Goma (DRC) and Gisenyi (Rwanda). In this paper we use numerical simulations to investigate the possibility of significantly reducing the lava flow hazard in the city through the construction of protective barriers. These barriers are added to the DEM of the area as additional morphological elements, and their effect is evaluated by repeating numerical simulations with and without the presence of barriers. A parametric study on barrier location, size, shape and orientation led to the identification of barriers which maximize protection while minimizing their impact. This study shows that the highest hazard area corresponding to eastern Goma, which was largely destroyed by lava flows in 2002, cannot be effectively protected from future lava flows towards Lake Kivu and should be abandoned. On the contrary, the rest of the town can be sheltered from lava flows by means of two barriers that deviate or contain the lava within the East Goma sector. A proposal for the future development of the town is formulated, whereby “new” Goma is completely safe from the arrival of lava flows originating from vents outside its boundaries. The proposal minimizes the risk of further destruction in town due to future lava flows.     

鄂尔多斯盆地西南部三叠纪末抬升剥蚀事件及热年代学记录

三叠纪末期大型鄂尔多斯盆地遭受了中生代成盆以来首次较大规模的抬升剥蚀,显著改造了中晚三叠世延长期盆地面貌,并控制了侏罗纪早期沉积格局和油藏分布,对盆地演化及矿产资源分布产生了重要影响。本文利用地质及大量钻井资料揭示了该期构造事件对盆地的剥蚀改造特征,盆地及周邻地区磷灰石裂变径迹年代学记录并约束了此次构造抬升的时限与过程;综合周邻区域构造研究成果,探讨了其发育的动力学背景。结果表明,鄂尔多斯盆地三叠纪末期的剥蚀具西南强、东北弱的特点,西南部大范围内延长组地层残缺不全,剥蚀量最大可达1000余米;前侏罗纪沉积古地貌总体呈西南高、北东低的特点;其抬升时间始于205～190Ma,西南部稍早于盆地腹部,抬升速率大于1℃/Myr,可持续至中侏罗世(约160Ma)。该期抬升剥蚀事件范围可涉及至盆地西南缘更广阔的区域,与同期秦岭造山带内出现的快速抬升冷却事件具有较好的时空耦合关系,是对秦岭造山带区域构造环境转变的响应和纪录。该研究丰富和发展了三叠纪末期构造事件在华北克拉通的影响,对该区油气、煤炭资源的进一步勘探和评价提供了新的思路,具有一定的现实意义。     

Robust numerical implementation of a 3D rate-dependent model for reservoir geomechanical simulations

A 3D elasto-plastic rate-dependent model for rock mechanics is formulated and implemented into a Finite Element (FE) numerical code. The model is based on the approach proposed by Vermeer and Neher (A soft soil model that accounts for creep. In: Proceedings of the International Symposium “Beyond 2000 in Computational Geotechnics,” pages 249-261, 1999). An original strain-driven algorithm with an Inexact Newton iterative scheme is used to compute the state variables for a given strain increment.The model is validated against laboratory measurements, checked on a simplified test case, and used to simulate land subsidence due to groundwater and hydrocarbon production. The numerical results prove computationally effective and robust, thus allowing for the use of the model on real complex geological settings.     

Parameterisation of the Planetary Boundary Layer for Diagnostic Wind Models

The planetary boundary-layer (PBL) parameterization is a key issue for the definition of initial wind flow fields in diagnostic models. However, PBL theories usually treat separately stable, neutral, and convective stability conditions, so that their implementation in diagnostic wind models is not straightforward. In the present paper, an attempt is made to adopt a comprehensive PBL parameterisation, covering stable/neutral and unstable atmospheric conditions, which appears suitable to diagnostic models. This parameterisation is implemented into our diagnostic mass-consistent code. A validation of the consistency between the implemented PBL parameterisations has been checked through an analysis of the sensitivity of the vertical wind profiles to atmospheric stability.