GIS技术在区域地震滑坡危险性预测中的应用以龙陵地震滑坡为例

地震滑坡是一种有着严重危害的次生地震灾害形式,形成机制复杂,涉及因素众多。运用G IS丰富的空间分析功能,对地震滑坡的影响因素进行研究,并进行潜在地震滑坡区的预测,是地震滑坡研究领域的一种新的发展趋势。本文在对1976年龙陵地震引发的地震滑坡分布特征研究的基础上,结合前人有关中国西南地区地震滑坡特征的研究成果,应用G IS对该区潜在地震滑坡危险区进行了预测。     

Maximum Expected Magnitude Assessment in a Geo Computer Environment: Case Study

Different approaches to seismic hazard assessment are compared. Each of them could be applied more or less successfully for territories which are stable in time and have a high level of seismic activity. A long-term seismic catalogue, not only including historical but also paleoearthquake data, is an essential requirement. But, in practice, such an ideal situation is very rare. Initial data is usually poor and short-term. Seismic hazard assessment could be more complicated for regions which are transient between relatively stable platforms and active mountain massifs. A new step in geoinformation technology for seismic hazard assessment based on a GEO computer environment is presented, its application is illustrated by the real case hazard evaluation for the territory of the Stavropol region, which is situated between the Russian platform and the Great Caucasus. The regional catalogue covers a time period of about 150 years. Though the region under consideration is not large, seismic activity is variable in space, from almost aseismic zones to rather active areas. GEO allows us to incorporate different techniques and all available information in the analysis, including those which are very difficult to formalize. The space distribution of the maximum expected earthquake magnitude is determined as a function of geological and geophysical data. An important feature of GEO is that it makes it possible to control the result of complicated algorithms through some relatively simple physical reasons.     

泥石流危险度野外快速评价方法探讨

文章系统地总结和分析了当前比较常用的泥石流危险度评价的研究方法,指出现行泥石流危险度评价中评价因子的提取往往需要详细野外勘查,工作量大、费用高,不便于对众多泥石流沟的普查;并提出一种在现场踏勘基础上,基于泥石流规模、频率和承灾体特征的快速便捷的新方法——泥石流危险度野外快速评价,并结合其他方法进行了实例对比验证。结果表明：泥石流危险度野外快速评价结果可靠,可作为快速评价泥石流危险性的新方法进行推广应用。     

Analysis of seismic hazard in landslide-prone regions: criteria and example for an area of Daunia (southern Italy)

In relation to the assessment of earthquake-induced landslide hazard, this paper discusses general principles and describes implementation criteria for seismic hazard estimates in landslide-prone regions. These criteria were worked out during the preparation of a hazard map belonging to the official Italian geological cartography and they are proposed as guidelines for future compilation of similar maps. In the presented case study, we used a procedure for the assessment of seismic hazard impact on slope stability adopting Arias intensity Ia as seismic shaking parameter and critical acceleration a _c as parameter representing slope strength to failures induced by seismic shaking. According to this procedure, after a preliminary comparison of estimated historical maximum values of Ia with values proposed in literature as landslide-triggering thresholds, a probabilistic approach, based on the Newmark’s model, is adopted: it allows to estimate the minimum critical acceleration a _c required for a slope to keep under a prefixed value, the probability of failures induced by seismic shakings expected in a given time interval. In this way, one can prepare seismic hazard maps where seismic shaking is expressed in an indirect way through a parameter (the critical acceleration) representing the “strength” that seismic shakings mobilise in slope materials (strength demand) with a prefixed exceedance probability. This approach was applied to an area of Daunia (Apulia—southern Italy) affected by frequent landslide phenomena. The obtained results indicate that shakings with a significant slope destabilisation potential can be expected particularly in the north-western part of the area, which is exposed to the seismic activity of Apennine tectonic structures.     

基于简化 Newmark位移模型的区域地震滑坡危险性快速评估以汶川MS 8.0级地震为例

以汶川MS8.0级地震重灾区的11县市为例,初步提出了基于简化Newmark位移模型的地震滑坡危险性应急快速评估方法。利用汶川地震即时地震动参数、工程地质岩性经验分组及地形坡度数据,借助ArcGIS空间数据建模工具编制了地震滑坡危险性快速评估流程模块。计算了区域浅表层饱和岩土体斜坡的静态安全系数Fs、临界加速度ac,并借此分析了地震滑坡易发性。利用经验式获得了汶川地震Arias强度和区域滑坡位移DN分布,实现了汶川地震重灾区地震滑坡危险性的快速评估,为应急救灾决策提供了参考。通过对比评估结果和震后滑坡调查成果,可知数十处灾难性滑坡绝大部分位于-高危险区的龙门山主中央断裂带两侧约20km地带中,显示了评估方法的可靠性; 同时,分析指出了空间数据精度及更新不足导致局部评估结果欠佳的局限性,并提出了改进建议。     

The SISMA prototype system: integrating Geophysical Modeling and Earth Observation for time-dependent seismic hazard assessment

An innovative approach to seismic hazard assessment is illustrated that, based on the available knowledge of the physical properties of the Earth structure and of seismic sources, on geodetic observations, as well as on the geophysical forward modeling, allows for a time-dependent definition of the seismic input. According to the proposed approach, a fully formalized system integrating Earth Observation data and new advanced methods in seismological and geophysical data analysis is currently under development in the framework of the Pilot Project SISMA, funded by the Italian Space Agency. The synergic use of geodetic Earth Observation data (EO) and Geophysical Forward Modeling deformation maps at the national scale complements the space- and time-dependent information provided by real-time monitoring of seismic flow (performed by means of the earthquake prediction algorithms CN and M8S) and permits the identification and routine updating of alerted areas. At the local spatial scale (tens of km) of the seismogenic nodes identified by pattern-recognition analysis, both GNSS (Global Navigation Satellite System) and SAR (Synthetic Aperture Radar) techniques, coupled with expressly developed models for interseismic phase, allow us to retrieve the deformation style and stress evolution within the seismogenic areas. The displacement fields obtained from EO data provide the input for the geophysical modeling, which eventually permits to indicate whether a specific fault is in a “critical state.” The scenarios of expected ground motion (shakemaps) associated with the alerted areas are then defined by means of full waveforms modeling, based on the possibility to compute synthetic seismograms by the modal summation technique (neo-deterministic hazard assessment). In this way, a set of deterministic scenarios of ground motion, which refer to the time interval when a strong event is likely to occur within the alerted area, can be defined both at national and at local scale. The considered integrated approach opens new routes in understanding the dynamics of fault zones as well as in modeling the expected ground motion. The SISMA system, in fact, provides tools for establishing warning criteria based on deterministic and rigorous forward geophysical models and hence allows for a well-controlled real-time prospective testing and validation of the proposed methodology over the Italian territory. The proposed approach complements the traditional probabilistic approach for seismic hazard estimates, since it supplies routinely updated information useful in assigning priorities for timely mitigation actions and hence it is particularly relevant to Civil Defense purposes.     

The Role of Data Processing and Uncertainty Management in Seismic Hazard Evaluations: Insights from Estimates in the Garfagnana–Lunigiana Area (Northern Italy)

Significantly different estimates of seismic hazard may result for the same site as aneffect of different methodological choices underlying the adopted procedures. In orderto explore this aspect, two approaches devoted to probabilistic seismic hazard assessment are considered for the evaluation of hazard in a seismic area in Northern Italy. In particular, results of a standard procedure are compared with those obtained by an innovative approach. Fundamental features of this last methodology are the extensive use of intensity data relative to seismic effects observed at the site of interest during past earthquakes and the basic role attributed to the parameterisation of uncertainty which affects the considered pieces ofinformation. The analysis indicates that the new approach supplies results significantlydifferent from those obtained from standard methodology and that these differences strongly depend on strategies adopted for data processing and for the management of uncertainties which affect input parameters.     

Estimation of tsunami hazard from volcanic activity — Suggested methodology with Augustine volcano,Alaska as an example

A general approach for the estimation of tsunami height and hazard in the vicinity of active volcanoes has been developed. An empirical relationship has been developed to estimate the height of the tsunami generated for an eruption of a given size. This relationship can be used to estimate the tsunami hazard based on the frequency of eruptive activity of a particular volcano. This technique is then applied to the estimation of tsunami hazard from the eruption of the Augustine volcano in Alaska. Modification of this approach to account for a less than satisfactory data base and differing volcanic characteristics is also discussed with the case of the Augustine volcano as an example. This approach can be used elsewhere with only slight modifications and, for the first time, provides a technique to estimate tsunami hazard from volcanic activity, similar to a well-established approach for the estimation of tsunami hazard from earthquake activity.     

Seismic hazard in the Sannio-Matese area of Italy

The assessment of seismic hazard at five selected sites in the Sannio-Matese region is based on the computer program SRAMSC. Owing to the extensive historical data base for the output parameter, the MSK intensiy is chosen. The seismicity model is made up of five narrow area seismic sources. Circular or elliptical macroseismic fields are assigned to individual sources. A generalized Kövesligethy equation is used for this purpose as the attenuation relationship. The study reveals similar and a rather high hazard at the sites at Benevento, Boiano, and Melfi, which are located in the zone of highest seismic activity. At the Pomigliano and Lucera sites, the assessed hazard is much lower.     

Local seismic hazard assessment in areas of weak to moderate seismicity—case study from Eastern Germany

Generally the seismic hazard of an area of interest is considered independent of time. However, its seismic risk or vulnerability, respectively, increases with the population and developing state of economy of the area. Therefore, many areas of moderate seismic hazard gain increasing importance with respect to seismic hazard and risk analysis. However, these areas mostly have a weak earthquake database, i.e., they are characterised by relative low seismicity and uncertain information concerning historical earthquakes. In a case study for Eastern Thuringia (Germany), acting as example for similar places in the world, seismic hazard is estimated using the probabilistic approach. Because of the lack of earthquakes occurring in the recent past, mainly historical earthquakes have to be used. But for these the actual earthquake sources or active faults, needed for the analysis, are imprecisely known. Therefore, the earthquake locations are represented by areal sources, a common practice. The definition of these sources is performed carefully, because their geometrical shape and size (apart from the earthquake occurrence model) influence the results significantly. Using analysis tools such as density maps of earthquake epicentres, seismic strain and energy release support this. Oversizing of areal sources leads to underestimation of seismic hazard and should therefore be avoided. Large location errors of historical earthquakes on the other hand are represented by several alternative areal sources with final superimposition of the different results. In a very similar way information known from macroseismic observations interpreted as source rather than as site effects are taken into account in order to achieve a seismic hazard assessment as realistic as possible. In very local cases the meaning of source effects exceeds those of site effects very likely. The influence of attenuation parameter variations on the result of estimated local seismic hazard is relatively low. Generally, the results obtained by the seismic hazard assessment coincide well with macroseismic observations from the thoroughly investigated largest earthquake in the region.     

Seismic hazard studies for Gaziantep city in South Anatolia of Turkey

Seismic hazard studies were conducted for Gaziantep city in the South Anatolia of Turkey. For this purpose, a new attenuation relationship was developed using the data of Zaré and Bard and accelerations were predicted employing this new equation. Deterministic approach, total probability theorem and GIS methodology were all together utilized for the seismic assessments. Seismic hazard maps with 0.25° grid intervals considering the site conditions were produced by the GIS technique. The results indicated that the acceleration values by the GIS hazard modelings were matched with the ones from the deterministic approach, however, they were underestimated comparing with the total probability theorem. In addition, the GIS based seismic hazard maps showed that the current seismic map of Turkey fairly yields conservative acceleration values for the Gaziantep region. Therefore, the constructed GIS hazard models are offered as a base map for a further modification of the current seismic hazard map.     

The comparison of the NDSHA, PSHA seismic hazard maps and real seismicity for the Italian territory

Rigorous and objective testing of seismic hazard assessments against the real seismic activity must become the necessary precondition for any responsible seismic risk estimation. Because seismic hazard maps seek to predict the shaking that would actually occur, the reference hazard maps for the Italian seismic code, obtained by probabilistic seismic hazard assessment (PSHA), and the alternative ground shaking maps based on the neo-deterministic approach (NDSHA), are cross-compared and tested against the real seismicity for the territory of Italy. The comparison between predicted intensities and those reported for past earthquakes shows that models generally provide rather conservative estimates, except for PGA with 10 % probability of being exceeded in 50 years, which underestimates the largest earthquakes. In terms of efficiency in predicting ground shaking, measured accounting for the rate of underestimated events and for the territorial extent of areas characterized by high seismic hazard, the NDSHA maps appear to outscore the PSHA ones.     

Stochastic seismic slope stability assessment using polynomial chaos expansions combined with relevance vector machine

This paper presents probabilistic assessment of seismically-induced slope displacements considering uncertainties of seismic ground motions and soil properties.A stochastic ground motion model representing both the temporal and spectral non-stationarity of earthquake shakings and a three-dimensional rotational failure mechanism are integrated to assess Newmark-type slope displacements.A new probabilistic approach that incorporates machine learning in metamodeling technique is proposed,by combining relevance vector machine with polynomial chaos expansions(RVM-PCE).Compared with other PCE methods,the proposed RVM-PCE is shown to be more effective in estimating failure probabilities.The sensitivity and relative influence of each random input parameter to the slope displacements are discussed.Finally,the fragility curves for slope displacements are established for sitespecific soil conditions and earthquake hazard levels.The results indicate that the slope displacement is more sensitive to the intensities and strong shaking durations of seismic ground motions than the frequency contents,and a critical Arias intensity that leads to the maximum annual failure probabilities can be identified by the proposed approach.     

Seismic Hazard Assessment in the Tihamat Asir Region, Southwestern Saudi Arabia

The occurrence of earthquakes, faulting of Pleistocene sediments, uplifting of Pleistocene coral reefs, recent incised wadis and lava effusions in addition to hot springs all clearly indicate that southeastern Saudi Arabia is tectonically active. This paper reviews the tectonic features of southwestern Saudi Arabia and provides new approaches and maps for the interpretation of old and recent earthquake data for improved assessment of the regional tectonics. A regionalized variable approach is used to develop earthquake groundmotion hazard maps for the region, based on geostatistical methods using the kriging technique. This hazard must be considered in any design/construction of engineering structures in the region. The application of kriging for estimating the ground shaking in the study region succeeded clearly in accomplishing its ultimate aim where the generated groundmotions are well correlated with the instrumental magnitude and historical intensity of any earthquake occurrence in the study region.     

Assessment of seismic hazard and liquefaction potential of Gujarat based on probabilistic approaches

Gujarat is one of the fastest-growing states of India with high industrial activities coming up in major cities of the state. It is indispensable to analyse seismic hazard as the region is considered to be most seismically active in stable continental region of India. The Bhuj earthquake of 2001 has caused extensive damage in terms of causality and economic loss. In the present study, the seismic hazard of Gujarat evaluated using a probabilistic approach with the use of logic tree framework that minimizes the uncertainties in hazard assessment. The peak horizontal acceleration (PHA) and spectral acceleration (Sa) values were evaluated for 10 and 2?% probability of exceedance in 50?years. Two important geotechnical effects of earthquakes, site amplification and liquefaction, are also evaluated, considering site characterization based on site classes. The liquefaction return period for the entire state of Gujarat is evaluated using a performance-based approach. The maps of PHA and PGA values prepared in this study are very useful for seismic hazard mitigation of the region in future.     

Seismic Hazard and Risk Assessment for Tulbagh,South Africa: Part I – Assessment of Seismic Hazard

In this part of our study the probabilistic seismic hazard analysis (PSHA) for Tulbagh was performed. The applied procedure is parametric and consists essentially of two steps. The first step is applicable to the area in the vicinity of Tulbagh and requires an estimation of the area-specific parameters, which, in this case, is the mean seismic activity rate, , the Gutenberg-Richter parameter, b, and the maximum regional magnitude, m_max. The second step is applicable to the Tulbagh site, and consists of parameters of distribution of amplitude of the selected ground motion parameter. The current application of the procedure provides an assessment of the PSHA in terms of peak ground acceleration (PGA) and spectral acceleration (SA). The procedure permits the combination of both historical and instrumental data. The historical part of the catalogue only contains the strongest events, whereas the complete part can be divided into several subcatalogues, each assumed complete above a specified threshold of magnitude. In the analysis, the uncertainty in the determination of the earthquake was taken into account by incorporation of the concept of `apparent magnitude'. The PSHA technique has been developed specifically for the estimation of seismic hazard at individual sites without the subjective judgement involved in the definition of seismic source zones, when the specific active faults have not been mapped or identified, and where the causes of seismicity are not well understood. The results of the hazard assessment are expressed as probabilities that specified values of PGA will be exceeded during the chosen time intervals, and similarly for the spectral accelerations. A worst case scenario sketches the possibility of a maximum PGA of 0.30g. The results of the hazard assessment can be used as input to a seismic risk assessment.     

Critical parameter investigation for earthquake hazard assessment in the Sannio-Matese area of Southern Italy

Probabilistic methods have been applied for the assessment of seismic hazard in a selected region of Southern Italy (Sannio-Matese). This method is mainly suitable for engineering and planning purposes and was first introduced in 1968 by Cornell and efficiently codified into a FORTRAN computer program by McGuire. Special attention is paid in this paper to the specific input parameters, i.e. completeness of data catalogues, time variability of seismic activity, different forms of frequency distributions, and regional attenuation characteristics.     

Seismic Hazard Estimation in Northern Algeria

In the present study, the seismic hazard in northern Algeria is estimated using both physical strain energy release and Gumbel's extreme values approaches. For six of the most industrial and populated cities in Algeria, seismic hazard is assessed and examined in greater detail. Gumbel's extreme values approach has been used to estimate seismic hazard in terms of magnitude and P.G.A at each point of an equispaced grid all over the north of Algeria. An average attenuation relationship for PGA has been provided using known relations which have been established in regions with similar attenuation characteristics.The results are presented mainly in the form of graphs and contour maps of magnitudes (respectively PGA) with a 60% probability of not being exceeded in the next 100 and 200 years. Globally, they give main features of northern Algeria in terms of zoning (as well as in terms of magnitude and in terms of PGA). They corroborate the ones obtained through other works, especially in the basin areas (Mitidja, Cheliff, Soumam and Constantine Basin).     

Seismic microzonation of Bangalore,India

In the present study, an attempt has been made to evaluate the seismic hazard considering local site effects by carrying out detailed geotechnical and geophysical site characterization in Bangalore, India to develop microzonation maps. An area of 220 km², encompassing Bangalore Mahanagara Palike (BMP) has been chosen as the study area. Seismic hazard analysis and microzonation of Bangalore are addressed in three parts: in the first part, estimation of seismic hazard is done using seismotectonic and geological information. Second part deals with site characterization using geotechnical and shallow geophysical techniques. In the last part, local site effects are assessed by carrying out one-dimensional (1-D) ground response analysis (using the program SHAKE2000) using both standard penetration test (SPT) data and shear wave velocity data from multichannel analysis of surface wave (MASW) survey. Further, field experiments using microtremor studies have also been carried out for evaluation of predominant frequency of the soil columns. The same has been assessed using 1-D ground response analysis and compared with microtremor results. Further, the Seed and Idriss simplified approach has been adopted to evaluate the soil liquefaction susceptibility and liquefaction resistance assessment. Microzonation maps have been prepared with a scale of 1:20,000. The detailed methodology, along with experimental details, collated data, results and maps are presented in this paper.     

地震滑坡危险性概念和基于力学模型的评估方法探讨

在我国大陆地区运用基于力学模型的Newmark位移分析方法开展地震滑坡危险性定量评估,尚处在起步阶段.为了进一步明确地震滑坡危险性概念和改进推广基于力学模型的评估方法,首先阐明了狭义的地震滑坡危险性预测评估与震后反演评估的关系; 同时为了应对地震应急、震后重建及潜在地震条件下的不同评估需求,初步提出了广义的地震滑坡危险性评估框架.随后申述了基于Newmark位移分析的地震滑坡危险性评估方法的理论基础、方法分类及最新进展,并以汶川地震滑坡危险性快速评估为例,剖析了目前影响评估有效性的不确定性及空间数据质量等问题,指出了基于力学模型的地震滑坡危险性评估方法的改进方向.建议开展潜在地震及其诱发滑坡危险性的耦合评估,建立适用于我国大陆地区地震滑坡位移分析的经验模型,以便为国家层面的地震滑坡危险性区划服务.