Do static stress changes of a moderate‐magnitude earthquake significantly modify the regional seismic hazard? Hints from the L'Aquila 2009 normal‐faulting earthquake (Mw 6.3, central Italy)

We investigated the Coulomb stress changes in the active faults surrounding a moderate‐magnitude normal‐faulting earthquake (2009 L'Aquila, Mw 6.3) and the associated variations in the expected ground motion on regional probabilistic seismic hazard maps. We show that the static stress variations can locally increase the seismic hazard by modifying the expected mean recurrence time on neighbouring faults by up to ~290 years, with associated variations in the probability of occurrence of the maximum expected earthquake of up to ~2%. Our findings suggest that the increase in seismic hazard on neighbouring faults following moderate‐magnitude earthquakes is probably not sufficient to necessitate systematic upgrades of regional probabilistic seismic hazard maps, but must be considered to better address and schedule strategies for local‐scale mitigation of seismic risk.     

A seismic landslide hazard analysis with topographic effect,a case study in the 99 Peaks region,Central Taiwan

It has been known that ground motion amplitude will be amplified at mountaintops; however, such topographic effects are not included in conventional landslide hazard models. In this study, a modified procedure that considers the topographic effects is proposed to analyze the seismic landslide hazard. The topographic effect is estimated by back analysis. First, a 3D dynamic numerical model with irregular topography is constructed. The theoretical topographic amplification factors are derived from the dynamic numerical model. The ground motion record is regarded as the reference motion in the plane area. By combining the topographic amplification factors with the reference motions, the amplified acceleration time history and amplified seismic intensity parameters are obtained. Newmark’s displacement model is chosen to perform the seismic landslide hazard analysis. By combining the regression equation and the seismic parameter of peak ground acceleration and Arias intensity, the Newmark’s displacement distribution is generated. Subsequently, the calculated Newmark’s displacement maps are transformed to the hazard maps. The landslide hazard maps of the 99 Peaks region, Central Taiwan are evaluated. The actual landslide inventory maps triggered by the 21 September 1999, Chi-Chi earthquake are compared with the calculated hazard maps. Relative to the conventional procedure, the results show that the proposed procedures, which include the topographic effect can obtain a better result for seismic landslide hazard analysis. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Liquefaction hazard mapping at the town of Edessa,Northern Greece

The town of Edessa is located on Northern Greece at a region that is characterized as low seismicity zone due to the fact that few moderate events of M < 6 occurred during the last century. According to the Greek Seismic Code, the expected acceleration having a 10% probability of being exceeded in 50 years is equal to 0.16g. However, an amplification of ground motion is likely to occur due the local geology that is consisted of Holocene fluvio-torrential deposits. The basic aim of this paper is to evaluate the site amplification due to geological conditions and to assess the liquefaction hazard. In order to achieve this, 1-D site response analyses were performed. The data that were employed for the construction of the numerical models have been collected from borings with standard penetrations tests (SPT) that were drilled for construction purposes. Afterward, the liquefaction potential of the subsoil layers was evaluated taking into consideration two seismic scenarios. The first scenario was based on the seismic parameters, earthquake magnitude and PGA, assigned by the Greek Seismic Code. On the second seismic model, we employed the values of acceleration, resulted from the 1-D analyses and the earthquake magnitude as it was defined by the Greek Seismic Code. In order to compile the liquefaction hazard maps, we initially estimated the liquefaction potential index (LPI) of the soil columns using the parameters provided by SPT, for both seismic loadings, and afterward we correlated these values with the proposed classification of the severity of liquefaction-induced deformations. In addition, having computed the value of probability based on the LPI, liquefaction manifestations probability maps were compiled for both scenarios. The result of this study was that liquefaction-induced ground disruptions are likely to occur at the center of the city, among the branches of Voda River, only when the amplified values of acceleration are taken into account to the computation of liquefaction potential.     

Seismic hazard maps of Mexico, the Caribbean, and Central and South America

The growth of megacities in seismically active regions around the world often includes the construction of seismically unsafe buildings and infrastructures due to an insufficient knowledge of existing seismic hazard and/or economic constraints. Minimization of the loss of life, property damage, and social and economic disruption due to earthquakes depends on reliable estimates of seismic hazard. We have produced a suite of seismic hazard estimates for Mexico, the Caribbean, and Central and South America. One of the preliminary maps in this suite served as the basis for the Caribbean and Central and South America portion of the Global Seismic Hazard Map (GSHM) published in 1999, which depicted peak ground acceleration (pga) with a 10% chance of exceedance in 50 years for rock sites. Herein we present maps depicting pga and 0.2 and 1.0 s spectral accelerations (SA) with 50%, 10%, and 2% chances of exceedance in 50 years for rock sites. The seismicity catalog used in the generation of these maps adds 3 more years of data to those used to calculate the GSH Map. Different attenuation functions (consistent with those used to calculate the U.S. and Canadian maps) were used as well. These nine maps are designed to assist in global risk mitigation by providing a general seismic hazard framework and serving as a resource for any national or regional agency to help focus further detailed studies required for regional/local needs. The largest seismic hazard values in Mexico, the Caribbean, and Central and South America generally occur in areas that have been, or are likely to be, the sites of the largest plate boundary earthquakes. High hazard values occur in areas where shallow-to-intermediate seismicity occurs frequently.     

Strong ground-motion prediction and uncertainties estimation for Delhi,India

An endeavour has been made in the present study to estimate the uncertainties in strong ground-motion estimation for Delhi region. The strong ground motion (SGM) has been estimated using the logic-tree approach with equal weights chosen for the branches of the logic tree because a scientific preference has not yet been developed for the area. Coefficient of variation (COV) maps have been generated along with the mean SGM which aid in visualizing the effect of our lack of knowledge on the final results and highlight those areas where improved scientific understanding can have an impact on future hazard maps. Monte Carlo simulation has been sued to consider the effect of the variation in the seismic hazard parameters. The spectral ground generated for various return periods suggests higher values in the north-west region with a decreasing trend towards the south-east part of Delhi. The hazard gradient is highest at spectral acceleration (Sa) 0.1 s and lowest at spectral acceleration (Sa) 1.0 s for all return periods. The highest COV values (~0.60) were observed in case of PGA while lowest COV values (~0.15) were observed for spectral acceleration (Sa) 0.3 s at all return periods. One of the conspicuous observations is that the trends of COV maps have been found to be governed by the boundaries of the seismogenic sources. While COV values are governed mostly by Line Sources in smaller return and time periods, in higher return periods trend of the COV maps were found to be governed by the boundaries of the Areal Sources.     

Spatial mapping of earthquake hazard parameters in the Hindukush–Pamir Himalaya and adjacent regions: Implication for future seismic hazard

The study deals spatial mapping of earthquake hazard parameters like annual and 100-years mode along with their 90% probability of not being exceeded (NBE) in the Hindukush–Pamir Himalaya and adjoining regions. For this purpose, we applied a straightforward and most robust method known as Gumbel’s third asymptotic distribution of extreme values (GIII). A homogeneous and complete earthquake catalogue during the period 1900–2010 with magnitude M_W ⩾ 4.0 is utilized to estimate these earthquake hazard parameters. An equal grid point mesh, of 1° longitude X 1° latitude, is chosen to produce detailed earthquake hazard maps. This performance allows analysis of the localized seismicity parameters and representation of their regional variations as contour maps. The estimated result of annual mode with 90% probability of NBE is expected to exceed the values of M_W 6.0 in the Sulaiman–Kirthar ranges of Pakistan and northwestern part of the Nepal and surroundings in the examined region. The 100-years mode with 90% probability of NBE is expected to exceed the value of M_W 8.0 in the Hindukush–Pamir Himalaya with Caucasus mountain belt, the Sulaiman–Kirthar ranges of Pakistan, northwestern part of the Nepal and surroundings, the Kangra–Himanchal Pradesh and Kashmir of India. The estimated high values of earthquake hazard parameters are mostly correlated with the main tectonic regimes of the examined region. The spatial variations of earthquake hazard parameters reveal that the examined region exhibits more complexity and has high crustal heterogeneity. The spatial maps provide a brief atlas of the earthquake hazard in the region.     

Seismic-event-based methodology to obtain earthquake-induced translational landslide regional hazard maps

A seismic-event-based methodology to generate earthquake-induced translational landslide maps using Newmark method is proposed. The steps are: (1) to construct a GIS-based geotechnical database; (2) to identify those areas that are susceptible to the occurrence of translational landslides based on available geological information; (3) to compute a static safety factor; (4) to compute the critical acceleration that defines the threshold acceleration required to cause a displacement; (5) to characterize the seismic hazard as a set of stochastic events, collectively exhaustive and mutually exclusive, that fully describes the hazard spatial distribution and annual frequency of occurrence (in accordance with the earthquake location, depth and magnitude) with the appropriate ground-motion prediction equations; (6) to compute the Newmark displacement; and finally, (7) to carry out a probabilistic translational landslide hazard analysis to estimate an exceedance rate of a given displacement. This methodology is applied to Mexico, and maps for return periods of 150 and 500 years are presented. Results shown in maps are estimations of where translational landslides may occur and should be useful to carry out local studies to elaborate recommendations of site specific hazard reduction plans as well as to calculate insurance rates. In addition, these results are useful to identify civil protection actions, risk management at regional and local level, and land use planning, as well as for promoting more detailed vulnerability and risk studies at different scales.     

Site-specific earthquake response study for hazard assessment in Kolkata city,India

The assessment of local site effects on seismic ground motions is of great importance in earthquake engineering practice. Several destructive earthquakes in the past have demonstrated that the amplification of ground motion and associated damage to structures due to local site conditions is a significant consideration in earthquake hazard analysis. A recent paper published in this journal highlights the hazard posed by earthquakes in the megacity of Kolkata in India due to its seismic and geological settings. The seismic hazard assessment study speculates that the deep alluvial deposit in the city may increase the seismic hazard probably due to the amplification of the seismic energies. This paper focuses on the seismic response studies of the various soil strata (i.e. for local subsurface conditions) obtained from various construction sites in the city for predicted earthquake. It is very well recognized that site response studies (a part of seismic microhazard zonation for urban areas) are the first step towards performance-based foundation design or seismic risk analysis and mitigation strategy. One of the problems for carrying out site-specific study in Kolkata is the lack of recorded strong motion data in the city. Hence, this paper outlines a methodology to carry out site-specific study, where no strong motion data or seismic data are available. The methodology uses wavelet-based spectrum compatibility approach to generate synthetic earthquake motions and equivalent linear method for seismic site response analysis. The Mega City of Kolkata has been considered to explain the methodology. Seismic hazard zonation map by the Bureau of Indian Standards classifies the City of Kolkata as moderate seismic zone (Zone III) with a zone factor 0.16. On the other hand, GSHAP(Global Seismic Hazard Assessment Program) map which is based on 10% probability of exceedance in 50 years specifies a maximum peak ground acceleration (PGA) of 1.6 m/s² (0.163 g) for this region. In the present study, the seismic response has been carried out based on GSHAP. The results of the analysis indicate the amplification of ground motion in the range of 4.46–4.82 with the fundamental period ranging from 0.81 to 1.17 s. Furthermore, the maximum spectral accelerations vary in the range of 0.78–0.95 g.     

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.     

Landslide risk analysis and its application in regional planning: an example from the highlands of the Outer Western Carpathians, Czech Republic

Using detailed field mapping, an analysis of landslide risk has been undertaken in the flysch highlands of the Outer Western Carpathians. The standardized Czech methodology of expert derived susceptibility zonation widely used for land development planning purposes and deterministic modeling of shallow landslides was used to separately assess the susceptibility of different landslide types. The two susceptibility zonation maps were used to define landslide hazard using information about landslide reactivation and the return periods of precipitation that triggered the respective landslide types. A risk matrix was then used to qualitatively analyze the landslide risk to selected assets. The monetary value of these assets, according to actual market prices, was calculated and analyzed with respect to the risk classification. Since the study area is an important residential and recreational area, the practical application of the derived results was checked through a series of interviews conducted with personnel of the local government planning and construction office. This demonstrated a willingness to apply the landslide hazard maps as well as restraints of its successful application. The main one is the absence of legally binding regulations to enforce the spatial planers to use this information.     

Assessing flood hazard using flood marks and analytic hierarchy process approach: a case study for the 2013 flood event in Quang Nam,Vietnam

The production of flood hazard assessment maps is an important component of flood risk assessment. This study analyses flood hazard using flood mark data. The chosen case study is the 2013 flood event in Quang Nam, Vietnam. The impacts of this event included 17 deaths, 230 injuries, 91,739 flooded properties, 11,530 ha of submerged and damaged agricultural land, 85,080 animals killed and widespread damage to roads, canals, dykes and embankments. The flood mark data include flood depth and flood duration. Analytic hierarchy process method is used to assess the criteria and sub-criteria of the flood hazard. The weights of criteria and sub-criteria are generated based on the judgements of decision-makers using this method. This assessment is combined into a single map using weighted linear combination, integrated with GIS to produce a flood hazard map. Previous research has usually not considered flood duration in flood hazard assessment maps. This factor has a rather strong influence on the livelihood of local communities in Quang Nam, with most agricultural land within the floodplain. A more comprehensive flood hazard assessment mapping process, with the additional consideration of flood duration, can make a significant contribution to flood risk management activities in Vietnam.     

Applying the HAZUS-MH software tool to assess seismic risk in downtown Ottawa,Canada

The aim of this paper is to present earthquake loss estimations for a portion of downtown Ottawa, Canada, using the HAZUS-MH (Hazards United States Multi-Hazard) software tool. The assessment is performed for a scenario earthquake of moment magnitude 6.5, at an epicentral distance of 15 km, occurring during business hours. A level 2 HAZUS-MH analysis was performed where the building inventory, microzonation studies, and site-specific ground motion hazard maps (2% exceedence probability in 50 years) were all improved based on local information. All collected data were assembled into a set of standard geodatabases that are compatible with the HAZUS-MH software using a GIS-specific procedure. The results indicate that the greatest losses are expected in unreinforced masonry buildings and commercial buildings. Sensitivity studies show that soil classes, the vulnerability of schools, and the spatial scale of loss estimations are also important factors to take into account.     

公路工程建设用地地质灾害危险性评估工作方法简介

文章详细地介绍了公路工程建设用地地质灾害危险性评估工作一些常见的基本工作方法,包括野外图件的选用,地质灾害危险性评估区范围的确定,各类地质灾害的测绘、描述、记录以及评估报告编写应注意的几个问题。     

Seismic hazard analysis of India using areal sources

In view of the major advancement made in understanding the seismicity and seismotectonics of the Indian region in recent times, an updated probabilistic seismic hazard map of India covering 6–38°N and 68–98°E is prepared. This paper presents the results of probabilistic seismic hazard analysis of India done using regional seismic source zones and four well recognized attenuation relations considering varied tectonic provinces in the region. The study area was divided into small grids of size 0.1° × 0.1°. Peak Horizontal Acceleration (PHA) and spectral accelerations for periods 0.1 s and 1 s have been estimated and contour maps showing the spatial variation of the same are presented in the paper. The present study shows that the seismic hazard is moderate in peninsular shield, but the hazard in most parts of North and Northeast India is high.     

Tsunami vulnerability assessment in urban areas using numerical model and GIS

Natural disasters can neither be predicted nor prevented. Urban areas with a high population density coupled with the construction of man-made structures are subjected to greater levels of risk to life and property in the event of natural hazards. One of the major and densely populated urban areas in the east coast of India is the city of Chennai (Madras), which was severely affected by the 2004 Tsunami, and mitigation efforts were severely dampened due to the non-availability of data on the vulnerability on the Chennai coast to tsunami hazard. Chennai is prone to coastal hazards and hence has hazard maps on its earth-quake prone areas, cyclone prone areas and flood prone areas but no information on areas vulnerable to tsunamis. Hence, mapping has to be done of the areas where the tsunami of December 2004 had directly hit and flooded the coastal areas in Chennai in order to develop tsunami vulnerability map for coastal Chennai. The objective of this study is to develop a GIS-based tsunami vulnerability map for Chennai by using a numerical model of tsunami propagation together with documented observations and field measurements of the evidence left behind by the tsunami in December 2004. World-renowned and the second-longest tourist beach in the world “Marina” present in this region witnessed maximum death toll due to its flat topography, resulting in an inundation of about 300 m landward with high flow velocity of the order of 2 m/s.     

Probabilistic seismic hazard of Pakistan, Azad-Jammu and Kashmir

The seismic hazard study for Pakistan and Azad Jammu and Kashmir has been conducted by using probabilistic approach in terms of peak ground acceleration (PGA) in m/s² and also seismic hazard response spectra for different cities. A new version of Ambraseys et al. (Bull Earthq Eng 3:1–53, 2005) ground acceleration model is used, and parameterization is based on most recent updated earthquake catalogs that consisted of 14,000 events. The threshold magnitude was fixed at M _w 4.8, but seismic zones like northern Pakistan–Tajikistan, Hindukush and northern Afghanistan–Tajikistan border had M _w 5.2. The average normalized ‘a’ and ‘b’ values for all zones are 6.15 and 0.95, respectively. Seismicity of study area was modeled, and ground motion was computed for eight frequencies (0.025, 0.1, 0.2, 0.5, 1.0, 1.5, 2.0, 2.5 s) for different annual exceedance rates of 0.02, 0.01, 0.005, 0.002 and 0.001 (return periods 50, 100, 200, 500 and 1,000 years) for stiff rocks at the gridding of 0.1° × 0.1°. Seismic hazard maps based on computed PGA for 0.02, 0.01 and 0.002 annual exceedance are prepared. These maps indicate the earthquake hazard of Pakistan and surrounding areas in the form of acceleration contour lines, which are in agreement with geological and seismotectonic characteristics of the study area. The maximum seismic hazard values are found at Muzaffarabad, Gilgit and Quetta areas.     

Probability-weighted hazard maps for comparing different flood risk management strategies: a case study

The study proposes an original methodology for producing probability-weighted hazard maps based on an ensemble of numerical simulations. These maps enable one to compare different strategies for flood risk management. The methodology was applied over a 270-km² flood-prone area close to the left levee system of a 28-km reach of the river Reno (Northern Central Italy). This reach is characterised by the presence of a weir that allows controlled flooding of a large flood-prone area during major events. The proposed probability-weighted hazard maps can be used to evaluate how a structural measure such as the mentioned weir alters the spatial variability of flood hazard in the study area. This article shows an application by constructing two different flood hazard maps: a first one which neglects the presence of the weir using a regular levee system instead, and a second one that reflects the actual geometry with the weir. Flood hazard maps were generated by combining the results of several inundation scenarios, simulated by coupling 1D- and 2D-hydrodynamic models.     

Probabilistic seismic hazard analysis for Bangalore

This article presents the results of probabilistic seismic hazard analysis (PSHA) for Bangalore, South India. Analyses have been carried out considering the seismotectonic parameters of the region covering a radius of 350 km keeping Bangalore as the center. Seismic hazard parameter ‘b’ has been evaluated considering the available earthquake data using (1) Gutenberg–Richter (G–R) relationship and (2) Kijko and Sellevoll (1989, 1992) method utilizing extreme and complete catalogs. The ‘b’ parameter was estimated to be 0.62 to 0.98 from G–R relation and 0.87 ± 0.03 from Kijko and Sellevoll method. The results obtained are a little higher than the ‘b’ values published earlier for southern India. Further, probabilistic seismic hazard analysis for Bangalore region has been carried out considering six seismogenic sources. From the analysis, mean annual rate of exceedance and cumulative probability hazard curve for peak ground acceleration (PGA) and spectral acceleration (Sa) have been generated. The quantified hazard values in terms of the rock level peak ground acceleration (PGA) are mapped for 10% probability of exceedance in 50 years on a grid size of 0.5 km × 0.5 km. In addition, Uniform Hazard Response Spectrum (UHRS) at rock level is also developed for the 5% damping corresponding to 10% probability of exceedance in 50 years. The peak ground acceleration (PGA) value of 0.121 g obtained from the present investigation is slightly lower (but comparable) than the PGA values obtained from the deterministic seismic hazard analysis (DSHA) for the same area. However, the PGA value obtained in the current investigation is higher than PGA values reported in the global seismic hazard assessment program (GSHAP) maps of Bhatia et al. (1999) for the shield area.     

Landslide Activity Maps for Landslide Hazard Evaluation: Three Case Studies from Southern Italy

This paper focuses on the application of landslideactivity maps for evaluating the mass movementhazard in selected areas of the Southern Apenninesof Italy: Bisaccia, Calitri, and Buoninventre. Theavailability of multi-year aerial photo coveragehelped to assess the morphological changes whichoccurred in the last 40 years. This information,integrated with historical data on slope instabilityand field checks, were used to produce landslideactivity maps. These maps represent a short-cut inthe assessment of mass movement hazard, because theyfocus on the effects of slope instability ratherthan on the causative conditions and processes;however, if kept simple and prepared at large scale,they may help the local administrators and land-useplanners to reduce the socio-economic costs oflandslides. Furthermore, the comparative study oflandslide activity represents a relativelyinexpensive and quick method for evaluating theperformance of the engineering control efforts.The quantification of landslide activity in terms ofareal frequency can represent an additional step,useful to determine the relative landslide hazard(zonation in more or less hazardous areas). Forexample, the estimates of areal frequency of activelandsliding for the last 40 years demonstrated thegreat influence of the 1980 Irpinia earthquake(M_s = 6.9) on the stability of slopes situatedclose to its epicenter (within a radius of about20 km).