Showing papers on "Peak ground acceleration published in 2014"

Three (nearly) complete inventories of landslides triggered by the May 12, 2008 Wenchuan Mw 7.9 earthquake of China and their spatial distribution statistical analysis

Abstract: The May 12, 2008 Wenchuan earthquake of China (Mw 7.9 or Ms 8.0) triggered hundreds of thousands of landslides. Mapping such a large number of landslides is a major task, considering the large size of the affected area and the availability of pre- and post-earthquake remote sensing images. This paper compares three (nearly) complete landslide inventories that were compiled from visual image interpretation. The three inventories differ in the manner in which the landslides are represented, either as polygons, centroid points, or top points. Landslides in the three inventories use one-to-one correspondence. Each of the three inventories includes a large proportion of the 197,481 landslides triggered by the earthquake. These landslides were delineated as individual solid polygons and points using visual interpretation of high-resolution aerial photographs and satellite images acquired following the earthquake and verified by selected field checking throughout a broad area of approximately 110,000 km2. These landslides cover a total area of approximately 1,160 km2. Based on the inventories of landslide polygons and landslide centroid points, two types of density maps were constructed. Correlations of landslide occurrence with seismic, geologic, and topographic parameters were analyzed using the three landslide inventories. Statistical analysis of their spatial distribution was performed using both the landslide area percentage (LAP), defined as the percentage of the area affected by the landslides and the landslide number density (LND), defined as the number of landslides per square kilometer. There are two types of LNDs: the LND-centroid (based on the centroid point of the landslide) and the LND-top (based on the top point of the landslide). We used the three indexes to determine how the occurrence of the landslides correlates with elevation, slope angle, slope aspect, slope position, slope curvature, lithology, distance from the epicenter, seismic intensity, distance from the Yingxiu-Beichuan surface fault rupture, peak ground acceleration (PGA), and coseismic surface displacements (including horizontal, vertical, and total displacements). Both the LAP and the two types of LND values were observed to have continuous positive or negative correlations with the slope angle, slope curvature, distance from the epicenter and from the Yingxiu-Beichuan surface fault rupture, seismic intensity, and coseismic surface displacement. In addition, the highest values of the LAP and LND values appear at ranges from 1,200 to 3,000 m in elevation. Moreover, the landslides have preferred orientations, dominated by the eastern, southeastern, and southern directions. In addition, the sandstone, siltstone (Z), and granitic rocks experienced more concentrated landslides. No obvious correlations were observed between the LAP and LND values and slope position. Finally, we studied the orders of eight earthquake-triggered landslide impact factor effect on landslide occurrence. The 197,481 landslides triggered by the 2008 Wenchuan earthquake were delineated. Three landslide inventories were constructed: polygon, centroid, and top point inventories. The landslides were spatially analyzed with topographic, lithology, and seismic parameters.

Earthquake induced landslide susceptibility mapping using an integrated ensemble frequency ratio and logistic regression models in West Sumatera Province, Indonesia

Abstract: An 8 Richter Scale (RS) earthquake struck West Sumatra on Wednesday, 30 September 2009, at 17.16 pm which led to huge number of landslides. Hence a comprehensive landslide susceptibility mapping (LSM) should be produced in order to reduce the damages to people and infrastructures. In the international landslide literature, various statistical methods such as frequency ratio (FR) and logistic regression (LR) have been widely used individually for LSM, but they have some weaknesses. FR which is able to perform bivariate statistical analysis (BSA) assesses the influence of classes of each conditioning factor on landslide occurrence. However, the correlation between the factors is mostly neglected. On the other hand, LR is able to analyze the relationship among the factors while it is not capable to evaluate the classes of each landslide conditioning factor. This paper aims to propose an ensemble method of FR and LR in order to overcome their weak points. For LSM, a landslide inventory map with a total of 87 landslide locations was extracted from various sources. Then the landslide inventory was randomly divided into two datasets 70% for training the models and the remaining 30% was used for validation purpose. The landslide conditioning factors consist of: altitude, curvature, river, SPI, rainfall, soil type, soil texture, land use/cover (LULC), peak ground acceleration (PGA), geology, slope, aspect, lineament and topographic wetness index (TWI). Four PGA of 7.5, 8, 8.6 and 9 were acquired and PGA 8 which was related to the 2009 earthquake was used to generate the model. Finally, the produced landslide susceptibility maps were validated using an area under the (ROC) curve method. For the model which was derived by PGA 8, the validation results showed 84% and 78% success and prediction rates respectively. Furthermore, the prediction rates for the models made by PGA 7.2, 8.6 and 9 are 79%, 78% and 81% respectively. The result proved the reasonable efficiency of the proposed method for earthquake induced landslide susceptibility mapping. Also the proposed ensemble method can be used in other hazard studies as it is capable to produce rapid and accurate assessment for disaster management and decision making.

Gradual unlocking of plate boundary controlled initiation of the 2014 Iquique earthquake

Abstract: A long foreshock series unlocked the South American plate boundary until eventually initiating the M 8.1 Iquique, Chile, earthquake. Two groups publishing in this issue of Nature analyse the seismic context of the Iquique earthquake that occurred off the coast of northern Chile on 1 April 2014 in a seismic zone that had been quiescent since a significant event in 1877. Gavin Hayes et al. identify areas of remaining or elevated earthquake hazard along the megathrust fault in the region, and conclude that the 2014 Iquique event was not the earthquake that had been anticipated. Given that significant sections of the northern Chile subduction zone have not ruptured in almost 150 years, they suggest that it is likely that future megathrust earthquakes will occur south and potentially north of the 2014 Iquique sequence. Bernd Schurr et al. show that the April 2014 earthquake broke a central fraction of the 'northern Chile seismic gap', the last major segment of the South American plate boundary that had yet to rupture in the past century. From July 2013 up to the April earthquake they identify three seismic clusters along this part of the plate boundary, each lasting a few weeks, with earthquakes of increasing peak magnitudes. They conclude that these seismic clusters and their slip transients reflect a gradual weakening of the central part of the seismic gap that was instrumental in initiating the final failure. On 1 April 2014, Northern Chile was struck by a magnitude 8.1 earthquake following a protracted series of foreshocks. The Integrated Plate Boundary Observatory Chile monitored the entire sequence of events, providing unprecedented resolution of the build-up to the main event and its rupture evolution. Here we show that the Iquique earthquake broke a central fraction of the so-called northern Chile seismic gap, the last major segment of the South American plate boundary that had not ruptured in the past century1,2. Since July 2013 three seismic clusters, each lasting a few weeks, hit this part of the plate boundary with earthquakes of increasing peak magnitudes. Starting with the second cluster, geodetic observations show surface displacements that can be associated with slip on the plate interface. These seismic clusters and their slip transients occupied a part of the plate interface that was transitional between a fully locked and a creeping portion. Leading up to this earthquake, the b value of the foreshocks gradually decreased during the years before the earthquake, reversing its trend a few days before the Iquique earthquake. The mainshock finally nucleated at the northern end of the foreshock area, which skirted a locked patch, and ruptured mainly downdip towards higher locking. Peak slip was attained immediately downdip of the foreshock region and at the margin of the locked patch. We conclude that gradual weakening of the central part of the seismic gap accentuated by the foreshock activity in a zone of intermediate seismic coupling was instrumental in causing final failure, distinguishing the Iquique earthquake from most great earthquakes. Finally, only one-third of the gap was broken and the remaining locked segments now pose a significant, increased seismic hazard with the potential to host an earthquake with a magnitude of >8.5.

Collapse Fragility of Steel Structures Subjected to Earthquake Mainshock-Aftershock Sequences

Abstract: This paper investigates the collapse probability of mainshock-damaged steel buildings in aftershocks, as an essential part of developing a framework to integrate aftershock seismic hazard into performance-based engineering (PBE). Analytical studies were conducted utilizing structural degradation models derived from existing publicly available NEEShub data. During earthquake events, aftershocks have the potential to cause severe damage to buildings and threaten life safety even when only minor damage is present from the mainshock. While aftershocks are normally somewhat smaller in magnitude, their ground motion intensity is not always smaller. Aftershocks may have a higher peak ground acceleration than the mainshock, even longer duration, and significantly different energy content as a result of the change in their location relative to the site. To date, the description of seismic hazard in PBE has not included the probability of aftershocks. In this study, the structural degradation model of a fou...

Pan-European ground-motion prediction equations for the average horizontal component of PGA, PGV, and 5 %-damped PSA at spectral periods up to 3.0 s using the RESORCE dataset

Abstract: This article presents a set of Ground-Motion Prediction Equations (GMPEs) for Europe and the Middle East, derived from the RESORCE strong motion data bank, following a standard regression approach. The parametric GMPEs are derived for the peak ground acceleration, peak ground velocity, and 5 %-damped pseudo-absolute acceleration response spectra computed over 23 periods between 0.02 and 3 s, considering the average horizontal-component ground-motions. The GMPEs are valid for distances less than 300 km, hypocentral depth up to 35 km and over the magnitude range 4–7.6. Two metrics for the source-to-station distance (i.e. Joyner-Boore and hypocentral) are considered. The selected dataset is composed by 2,126 recordings (at a period of 0.1 s) related to 365 earthquakes, that includes strong-motion data from 697 stations.The EC8 soil classification (four classes from A to D) discriminates recording sites and four classes (normal, reverse, strike-slip, and unspecified) describe the style of faulting. A subset which contains only stations with measured Vs30 and earthquakes with specified focal mechanism (1,224 records from 345 stations and 255 earthquakes) is used to test of the accuracy of the median prediction and the variability associated to the broader data set. A random effect regression scheme is applied and bootstrap analyses are performed to estimate the 95 % confidence levels for the parameters. The total standard deviation sigma is decomposed into between-events and within-event components, and the site-to-site component is evaluated as well. The results show that the largest contribution to the total sigma is coming from the within-event component. When analyzing the residual distributions, no significant trends are observed that can be ascribed to the earthquake type (mainshock-aftershock classification) or to the non-linear site effects. The proposed GMPEs have lower median values than global models at short periods and large distances, while are consistent with global models at long periods $$(\hbox {T} > 1)$$ s. Consistency is found with two regional models developed for Turkey and Italy, as the considered dataset is dominated by waveforms recorded in these regions.

A seismological model for earthquakes induced by fluid extraction from a subsurface reservoir

Abstract: A seismological model is developed for earthquakes induced by subsurface reservoir volume changes. The approach is based on the work of Kostrov (1974) and McGarr (1976) linking total strain to the summed seismic moment in an earthquake catalog. We refer to the fraction of the total strain expressed as seismic moment as the strain partitioning function, α. A probability distribution for total seismic moment as a function of time is derived from an evolving earthquake catalog. The moment distribution is taken to be a Pareto Sum Distribution with confidence bounds estimated using approximations given by Zaliapin et al. (2005). In this way available seismic moment is expressed in terms of reservoir volume change and hence compaction in the case of a depleting reservoir. The Pareto Sum Distribution for moment and the Pareto Distribution underpinning the Gutenberg-Richter Law are sampled using Monte Carlo methods to simulate synthetic earthquake catalogs for subsequent estimation of seismic ground motion hazard. We demonstrate the method by applying it to the Groningen gas field. A compaction model for the field calibrated using various geodetic data allows reservoir strain due to gas extraction to be expressed as a function of both spatial position and time since the start of production. Fitting with a generalized logistic function gives an empirical expression for the dependence of α on reservoir compaction. Probability density maps for earthquake event locations can then be calculated from the compaction maps. Predicted seismic moment is shown to be strongly dependent on planned gas production.

The structural role played by masonry infills on RC building performances after the 2011 Lorca, Spain, earthquake

Abstract: On May 11, 2011 an earthquake of magnitude 5.1 ( $$M_{w}$$ ) struck Murcia region causing nine casualties and damage to buildings and infrastructures. Even if the main characteristics of the event would classify it as a moderate earthquake, the maximum Peak Ground Acceleration (PGA) registered (equal to 0.37 g) exceeded significantly local code provisions in terms of hazard at the site. This high PGA was a result of directivity effects in the near source region. An overview of earthquake characteristics and damage observed is provided. Notwithstanding the lack of proper structural design characterizing building stock in the area, most of the losses were caused by non-structural damage. According to in field observations, it emerges that masonry infills provided additional, “not designed”, strength to reinforced concrete (RC) buildings. Observed damage data, collected after the earthquake, are shown and compared to the results of a simplified approach for nonstructural damage assessment of RC infilled structures (FAST vulnerability approach). The latter comparison provided a fair accordance between observed data and analytical results.

Strong ground motion prediction using virtual earthquakes

Abstract: Sedimentary basins increase the damaging effects of earthquakes by trapping and amplifying seismic waves Simulations of seismic wave propagation in sedimentary basins capture this effect; however, there exists no method to validate these results for earthquakes that have not yet occurred We present a new approach for ground motion prediction that uses the ambient seismic field We apply our method to a suite of magnitude 7 scenario earthquakes on the southern San Andreas fault and compare our ground motion predictions with simulations Both methods find strong amplification and coupling of source and structure effects, but they predict substantially different shaking patterns across the Los Angeles Basin The virtual earthquake approach provides a new approach for predicting long-period strong ground motion

Comprehensive observation and modeling of earthquake and temperature‐related seismic velocity changes in northern Chile with passive image interferometry

Abstract: We report on earthquake and temperature-related velocity changes in high-frequency autocorrelations of ambient noise data from seismic stations of the Integrated Plate Boundary Observatory Chile project in northern Chile. Daily autocorrelation functions are analyzed over a period of 5 years with passive image interferometry. A short-term velocity drop recovering after several days to weeks is observed for the Mw 7.7 Tocopilla earthquake at most stations. At the two stations PB05 and PATCX, we observe a long-term velocity decrease recovering over the course of around 2 years. While station PB05 is located in the rupture area of the Tocopilla earthquake, this is not the case for station PATCX. Station PATCX is situated in an area influenced by salt sediment in the vicinity of Salar Grande and presents a superior sensitivity to ground acceleration and periodic surface-induced changes. Due to this high sensitivity, we observe a velocity response of several regional earthquakes at PATCX, and we can show for the first time a linear relationship between the amplitude of velocity drops and peak ground acceleration for data from a single station. This relationship does not hold true when comparing different stations due to the different sensitivity of the station environments. Furthermore, we observe periodic annual velocity changes at PATCX. Analyzing data at a temporal resolution below 1 day, we are able to identify changes with a period of 24 h, too. The characteristics of the seismic velocity with annual and daily periods indicate an atmospheric origin of the velocity changes that we confirm with a model based on thermally induced stress. This comprehensive model explains the lag time dependence of the temperature-related seismic velocity changes involving the distribution of temperature fluctuations, the relationship between temperature, stress and velocity change, plus autocorrelation sensitivity kernels.

Towards fully data driven ground-motion prediction models for Europe

Abstract: We have used the Artificial Neural Network method (ANN) for the derivation of physically sound, easy-to-handle, predictive ground-motion models from a subset of the Reference database for Seismic ground-motion prediction in Europe (RESORCE). Only shallow earthquakes (depth smaller than 25 km) and recordings corresponding to stations with measured $$V_{s30}$$ properties have been selected. Five input parameters were selected: the moment magnitude $$M_{W}$$ , the Joyner–Boore distance $$R_{JB}$$ , the focal mechanism, the hypocentral depth, and the site proxy $$V_{S30}$$ . A feed-forward ANN type is used, with one 5-neuron hidden layer, and an output layer grouping all the considered ground motion parameters, i.e., peak ground acceleration (PGA), peak ground velocity (PGV) and 5 %-damped pseudo-spectral acceleration (PSA) at 62 periods from 0.01 to 4 s. A procedure similar to the random-effects approach was developed to provide between and within event standard deviations. The total standard deviation ( $$\sigma $$ ) varies between 0.298 and 0.378 (log $$_{10}$$ unit) depending on the period, with between-event and within-event variabilities in the range 0.149–0.190 and 0.258–0.327, respectively. Those values prove comparable to those of conventional GMPEs. Despite the absence of any a priori assumption on the functional dependence, our results exhibit a number of physically sound features: magnitude scaling of the distance dependency, near-fault saturation distance increasing with magnitude, amplification on soft soils and even indications for nonlinear effects in softer soils.

Ground Motion and Seismic Source Aspects of the Canterbury Earthquake Sequence

Abstract: This paper provides an overview of the ground motion and seismic source aspects of the Canterbury earthquake sequence. Common reported attributes among the largest earthquakes in this sequence are ...

Consideration of aging and SSI effects on seismic vulnerability assessment of RC buildings

Abstract: At present, the seismic vulnerability assessment of reinforced concrete (RC) buildings is made considering fixed base conditions; moreover, the mechanical properties of the building remain intact in time. In this study we investigate whether these two fundamental hypotheses are sound as aging and soil-structure interaction (SSI) effects might play a crucial role in the seismic fragility analysis of RC structures. Among the various aging processes, we consider the chloride-induced corrosion based on probabilistic modeling of corrosion initiation time and corrosion rate. Different corrosion aspects are considered in the analysis including the loss of reinforcement cross-sectional area, the degradation of concrete cover and the reduction of steel ultimate deformation. SSI is modeled by applying the direct one-step approach, which accounts simultaneously for inertial and kinematic interactions. Two-dimensional incremental dynamic analysis is performed to assess the seismic performance of the initial uncorroded ( $$\hbox {t}=0$$ years) and corroded ( $$\hbox {t}=50$$ years) RC moment resisting frame structures, having been designed with different seismic code levels. The time-dependent fragility functions are derived in terms of the spectral acceleration at the fundamental mode of the structure $$\hbox {S}_{\mathrm{a}}(\hbox {T}_{1}$$ , 5 %) and the outcropping peak ground acceleration for the immediate occupancy and collapse prevention limit states. Results show an overall increase in seismic vulnerability over time due to corrosion highlighting the important influence of deterioration due to aging effects on the structural behavior. Moreover, the consideration of SSI and site effects may significantly alter the expected structural performance leading to higher vulnerability values.

Seismic vulnerability assessment of school buildings in Tehran city based on AHP and GIS

Abstract: . The objective of the current study is to evaluate the seismic vulnerability of school buildings in Tehran city based on the analytic hierarchy process (AHP) and geographical information system (GIS). To this end, the peak ground acceleration, slope, and soil liquefaction layers were utilized for developing a geotechnical map. Also, the construction materials of structures, age of construction, the quality, and the seismic resonance coefficient layers were defined as major factors affecting the structural vulnerability of school buildings. Then, the AHP method was applied to assess the priority rank and weight of criteria (layers) and alternatives (classes) of each criterion via pairwise comparison in all levels. Finally, the geotechnical and structural spatial layers were overlaid to develop the seismic vulnerability map of school buildings in Tehran. The results indicated that only in 72 (about 3%) out of 2125 school buildings of the study area will the destruction rate be very high and therefore their reconstruction should seriously be considered.

Ground Motion Prediction Equations in the San Jacinto Fault Zone: Significant Effects of Rupture Directivity and Fault Zone Amplification

Abstract: We present a new set of Ground Motion Prediction Equations (GMPEs) for horizontal Peak Ground Acceleration, Peak Ground Velocity, and 5 % damped pseudo-spectral acceleration (PSA), developed for the San Jacinto Fault Zone (SJFZ) area. Besides using these equations to quantify seismic shaking in the area, the results allow us to examine the physics and local properties controlling the observed ground motions. The analyzed dataset includes ~30,000 observations from ~800 events spanning a magnitude range of 1.5 < M < 6.0 and recorded by up to 140 stations at epicentral distances ranging from essentially zero to 150 km. The local GMPE is developed for the SJFZ by applying classical regression techniques with predictive variables that include first distance and magnitude, and then site characteristics, rupture directivity, and fault zone amplification. The significance of these effects is determined by measuring the uncertainty-reduction of the GMPE due to each factor. The results show that, in contrast to many regional studies, traditional site characteristic has a relatively minor effect on peak amplitudes in our study area. However, rupture directivity is a significant factor controlling the amplitudes of ground motion even for small events. The dense seismic network and newly developed directivity tool enable us to extract efficiently directivity effects with statistical significance, using the ground-motion dataset during the regression analysis process. The obtained rupture directivities are consistent with the main focal mechanism orientations and surface trace orientations, known from other studies, and predictions for bimaterial ruptures in the trifurcation area of the SJFZ. Fault zone amplification is a second important factor, showing strong impact on the peak ground motion values, with increasing role for the lower frequency range (<10 Hz) examined in the 5 % damped PSA values. We also observe signatures of large amplitude-variances, which indicate additional source-related control on the distribution of amplitudes (besides rupture directivity) for aftershocks close in time and location to the M L 5.1 earthquake of March 2013. Using the full set of records we present the most complete set of GMPEs for the SJFZ area, including a higher-amplitude prediction for regions in the direction of rupture.

Understanding the Magnitude Dependence of PGA and PGV in NGA‐West 2 Data

Abstract: The Next Generation Attenuation‐West 2 (NGA‐West 2) 2014 ground‐motion prediction equations (GMPEs) model ground motions as a function of magnitude and distance, using empirically derived coefficients (e.g., Bozorgnia et al. , 2014); as such, these GMPEs do not clearly employ earthquake source parameters beyond moment magnitude ( M ) and focal mechanism. To better understand the magnitude‐dependent trends in the GMPEs, we build a comprehensive earthquake source‐based model to explain the magnitude dependence of peak ground acceleration and peak ground velocity in the NGA‐West 2 ground‐motion databases and GMPEs. Our model employs existing models (Hanks and McGuire, 1981; Boore, 1983, 1986; Anderson and Hough, 1984) that incorporate a point‐source Brune model, including a constant stress drop and the high‐frequency attenuation parameter κ , random vibration theory, and a finite‐fault assumption at the large magnitudes to describe the data from magnitudes 3 to 8. We partition this range into four different magnitude regions, each of which has different functional dependences on M . Use of the four magnitude partitions separately allows greater understanding of what happens in any one subrange, as well as the limiting conditions between the subranges. This model provides a remarkably good fit to the NGA data for magnitudes from 3< M <8 at close rupture distances ( R rup≤20 km). We explore the trade‐offs between Δ σ and κ in ground‐motion models and data, which play an important role in understanding small‐magnitude data, for which the corner frequency is masked by the attenuation of high frequencies. That this simple, source‐based model matches the NGA‐West 2 GMPEs and data so well suggests that considerable simplicity underlies the parametrically complex NGA GMPEs. Online Material: Figures providing detail on the V S 30 distribution in the subset of the Next Generation Attenuation‐West 2 (NGA‐West 2) data used, the quarter‐wavelength amplifications used in the model, the statistical test for the large magnitude portion of the model, and the magnitude dependence.

Seismic Activities and Earthquake Potential in the Tibetan Plateau

Abstract: The Tibetan Plateau is a region with intensive tectonic deformation and high level earthquake activity. On the basis of historical records, there have been 18 earthquakes with magnitude over 8 and more than 100 events with magnitude between 7 and 7.9 in the plateau. All of these events occured either along plate boundary of the Himalaya or along boundaries among the intra-plate fault blocks as well as sub-blocks. There have been three phases of earthquake clustering in the Tibetan Plateau, 1920 to 1937, 1947 to 1976, and 1995 to present, since instrumental seismic records began in 1900. Each phase of the cluster is characterized by events of magnitude over 8 and a number of events of magnitude over 7 to form an earthquake series. The three phases are the Haiyuan-Gulang earthquake series in 1920s to 1930s, the Zayu-Damxung earthquake series in 1950s to 1970s, and the Kunlun-Wenchuan series since the beginning of this century. Each earthquake series also has its major active region of occurrences. For example, the major active region of occurrence is associated with the newest phase of earthquake series is the Bayan Har fault block. The episodes of seismic activity in the Tibetan Plateau coincide with the high level activities of global mega-earthquake (Mw ≥ 8) activity. The Kunlun-Wenchuan earthquake series concurs with the global earthquake activity since 2001 suggesting geodynamic connections between them. Based on careful comparative studies, we suggest the current phase of earthquake activity continues today, giving that the global mega-earthquakes (Mw8∼9) occur frequently along plate boundaries and major intra-plate earthquakes (Mw7∼8) occur along boundaries of crustal fault block. Through studies of relationships between global and regional earthquake activities and of temporal and spatial distribution of the earthquakes in the Tibetan Plateau, this paper summarized general tendencies of current earthquake activity, assesses recent earthquake potential, and proposes suggestions for intensive earthquake monitoring.

Landslides triggered by the 12 January 2010 Port-au-Prince, Haiti, M w = 7.0 earthquake: visual interpretation, inventory compiling, and spatial distribution statistical analysis

Abstract: . The 12 January 2010 Port-au-Prince, Haiti, earthquake (Mw= 7.0) triggered tens of thousands of landslides. The purpose of this study is to investigate the correlations of the occurrence of landslides and the thicknesses of their erosion with topographic, geologic, and seismic parameters. A total of 30 828 landslides triggered by the earthquake covered a total area of 15.736 km2, distributed in an area more than 3000 km2, and the volume of landslide accumulation materials is estimated to be about 29 700 000 m3. These landslides are of various types, mostly belonging to shallow disrupted landslides and rock falls, but also include coherent deep-seated landslides and rock slides. These landslides were delineated using pre- and post-earthquake high-resolution satellite images. Spatial distribution maps and contour maps of landslide number density, landslide area percentage, and landslide erosion thickness were constructed in order to analyze the spatial distribution patterns of co-seismic landslides. Statistics of size distribution and morphometric parameters of co-seismic landslides were carried out and were compared with other earthquake events in the world. Four proxies of co-seismic landslide abundance, including landslides centroid number density (LCND), landslide top number density (LTND), landslide area percentage (LAP), and landslide erosion thickness (LET) were used to correlate co-seismic landslides with various environmental parameters. These parameters include elevation, slope angle, slope aspect, slope curvature, topographic position, distance from drainages, lithology, distance from the epicenter, distance from the Enriquillo–Plantain Garden fault, distance along the fault, and peak ground acceleration (PGA). A comparison of these impact parameters on co-seismic landslides shows that slope angle is the strongest impact parameter on co-seismic landslide occurrence. Our co-seismic landslide inventory is much more detailed than other inventories in several previous publications. Therefore, we carried out comparisons of inventories of landslides triggered by the Haiti earthquake with other published results and proposed possible reasons for any differences. We suggest that the empirical functions between earthquake magnitude and co-seismic landslides need to be updated on the basis of the abundant and more complete co-seismic landslide inventories recently available.

Statistical analysis of landslides caused by the Mw 6.9 Yushu, China, earthquake of April 14, 2010

Abstract: The April 14, 2010 Yushu, China, earthquake (Mw 6.9) triggered a great number of landslides. At least 2,036 co-seismic landslides, with a total coverage area of 1.194 km2, were delineated by visual interpretation of aerial photographs and satellite images taken following the earthquake, and verified by field inspection. Based on the mapping results, a statistical analysis of the spatial distribution of these landslides is performed using the landslide area percentage (LAP), defined as the percentage of the area affected by the landslides, and landslide number density (LND), defined as the number of landslides per square kilometer. The purpose is to clarify how the landslides correlate the control factors, which are the elevation, slope angle, slope aspect, slope position, distance from drainages, lithology, distance from the surface rupture, and peak ground acceleration (PGA). The results show that both LAP and LND have strongly positive correlations with slope angle and negative correlations with distance from the surface rupture and distance from drainages. The highest LAP and LPD values are in places of elevations from 3,800 to 4,000 m. The slopes producing landslides are mostly facing toward NE, E, and SE. The geological units of Q4 al-pl, N, and T3 kn 1 have the highest concentrations of co-seismic landslides. No apparent correlations are present between LAP and LND values and PGA. On both sides of the surface rupture, the landslide distributions are almost similar except a few exceptions, likely associated with the nature of the strike-slip seismogenic fault for this event. The bivariate statistical analysis shows that, in descending order, the earthquake-triggered landslide impact factors are distance from surface rupture > slope angle > distance from drainages > lithology > PGA. Besides, as the detailed co-seismic landslides inventories related to strike-slip earthquakes are still few compared with that of thrusting-fault earthquakes, this case study would shed new light on the subject. For instance, the landslide spatial distribution on both sides of the strike-slip seismogenic fault is rather different from that of thrusting-fault earthquakes. It reminds us to take different strategies of measures for prevention and mitigation of landslides induced by earthquakes with different mechanisms.

Application of Horizontal‐to‐Vertical Spectral Ratios of Earthquake Ground Motions to Identify Subsurface Structures at and around the K‐NET Site in Tohoku, Japan

Abstract: We propose an optimal way to use horizontal‐to‐vertical spectral ratios (HVRs) for subsurface structure exploration, based on the diffuse field concept (Kawase et al. , 2011; Sanchez‐Sesma et al. , 2011). This approach is applicable to both earthquake and microtremor ground motions. We show here analyses of the observed ground‐motion data at and around a K‐NET station in Miyagi Prefecture, Japan, where very large peak horizontal ground acceleration was observed during the earthquake of 11 March 2011 off the Pacific coast of Tohoku, Japan. We compare HVRs of the strong motions for the mainshock and the largest peak acceleration aftershock with those averaged over tens of weak motions to observe soil nonlinearity effects on the HVRs. Then, we determine detailed velocity profiles from the HVRs at the K‐NET Tsukidate station and the temporary aftershock observation sites. We find that HVRs can be explained quite well by the identified velocity profiles at all the target sites. The observed peak at 9 Hz for the averaged weak‐motion data originates in the topmost layers lying over the engineering bedrock.

Earthquake-induced ground failures in Italy from a reviewed database

Abstract: . A database (Italian acronym CEDIT) of earthquake-induced ground failures in Italy is presented, and the related content is analysed. The catalogue collects data regarding landslides, liquefaction, ground cracks, surface faulting and ground changes triggered by earthquakes of Mercalli epicentral intensity 8 or greater that occurred in the last millennium in Italy. As of January 2013, the CEDIT database has been available online for public use ( http://www.ceri.uniroma1.it/cn/gis.jsp ) and is presently hosted by the website of the Research Centre for Geological Risks (CERI) of the Sapienza University of Rome. Summary statistics of the database content indicate that 14% of the Italian municipalities have experienced at least one earthquake-induced ground failure and that landslides are the most common ground effects (approximately 45%), followed by ground cracks (32%) and liquefaction (18%). The relationships between ground effects and earthquake parameters such as seismic source energy (earthquake magnitude and epicentral intensity), local conditions (site intensity) and source-to-site distances are also analysed. The analysis indicates that liquefaction, surface faulting and ground changes are much more dependent on the earthquake source energy (i.e. magnitude) than landslides and ground cracks. In contrast, the latter effects are triggered at lower site intensities and greater epicentral distances than the other environmental effects.

Fully probabilistic seismic displacement analysis of spatially distributed slopes using spatially correlated vector intensity measures

Abstract: SUMMARY Earthquake-induced slope displacement is an important parameter for safety evaluation and earthquake design of slope systems. Traditional probabilistic seismic hazard analysis usually focuses on evaluating slope displacement at a particular location, and it is not suitable for spatially distributed slopes over a large region. This study proposes a computationally efficient framework for fully probabilistic seismic displacement analysis of spatially distributed slope systems using spatially correlated vector intensity measures (IMs). First, a spatial cross-correlation model for three key ground motion IMs, that is, peak ground acceleration (PGA), Arias intensity, and peak ground velocity, is developed using 2686 ground motion recordings from 11 recent earthquakes. To reduce the computational cost, Monte Carlo simulation and data reduction techniques are utilized to generate spatially correlated random fields for the vector IMs. The slope displacement hazards over the region are further quantified using empirical predictive equations. Finally, an illustrative example is presented to highlight the importance of the spatial correlation and the advantage of using spatially correlated vector IMs in seismic hazard analysis of spatially distributed slopes. Copyright © 2013 John Wiley & Sons, Ltd.

Probabilistic assessment of the seismic performance of earth slopes

Abstract: Sliding block displacements are used to evaluate the potential for seismic slope instability. Deterministic approaches are typically used to predict the expected level of sliding block displacement, although they do not rigorously account for uncertainties in the expected ground shaking, dynamic response, or displacement prediction. As a result, there is no concept of the actual hazard associated with the displacement computed by the deterministic approach. This paper summarizes and extends recent developments related to the probabilistic assessment of sliding block displacements. The probabilistic approach generates a hazard curve for displacement in which the annual rate of exceedance for a range of displacement levels is computed. The probabilistic approach is formulated both in terms of scalar hazard analysis (i.e., using one ground motion parameter, peak ground acceleration) and vector hazard analysis (i.e., using two ground motion parameters, peak ground acceleration and peak ground velocity), and applied both to rigid and flexible sliding block conditions. Generally, the vector probabilistic approach predicts displacements that are 2–3 times smaller than the scalar probabilistic approach, revealing the value of characterizing frequency content via peak ground velocity. Comparisons between the deterministic and probabilistic approaches, in either a scalar or vector context, indicate that the deterministic approach can severely underestimate displacements relative to the probabilistic approach because it ignores the aleatory variabilities in the dynamic and sliding responses of the sliding mass. This under-prediction is most significant for longer period sliding masses. Modifications to the deterministic approach are proposed that provide displacements that are more consistent with the probabilistic approach.

Neo-deterministic seismic hazard assessment in North Africa

Abstract: North Africa is one of the most earthquake-prone areas of the Mediterranean. Many devastating earthquakes, some of them tsunami-triggering, inflicted heavy loss of life and considerable economic damage to the region. In order to mitigate the destructive impact of the earthquakes, the regional seismic hazard in North Africa is assessed using the neo-deterministic, multi-scenario methodology (NDSHA) based on the computation of synthetic seismograms, using the modal summation technique, at a regular grid of 0.2 × 0.2°. This is the first study aimed at producing NDSHA maps of North Africa including five countries: Morocco, Algeria, Tunisia, Libya, and Egypt. The key input data for the NDSHA algorithm are earthquake sources, seismotectonic zonation, and structural models. In the preparation of the input data, it has been really important to go beyond the national borders and to adopt a coherent strategy all over the area. Thanks to the collaborative efforts of the teams involved, it has been possible to properly merge the earthquake catalogues available for each country to define with homogeneous criteria the seismogenic zones, the characteristic focal mechanism associated with each of them, and the structural models used to model wave propagation from the sources to the sites. As a result, reliable seismic hazard maps are produced in terms of maximum displacement (D max), maximum velocity (V max), and design ground acceleration.

Earthquake scenario in West Bengal with emphasis on seismic hazard microzonation of the city of Kolkata, India

Abstract: . Seismic microzonation is a process of estimating site-specific effects due to an earthquake on urban centers for its disaster mitigation and management. The state of West Bengal, located in the western foreland of the Assam–Arakan Orogenic Belt, the Himalayan foothills and Surma Valley, has been struck by several devastating earthquakes in the past, indicating the need for a seismotectonic review of the province, especially in light of probable seismic threat to its capital city of Kolkata, which is a major industrial and commercial hub in the eastern and northeastern region of India. A synoptic probabilistic seismic hazard model of Kolkata is initially generated at engineering bedrock (Vs30 ~ 760 m s−1) considering 33 polygonal seismogenic sources at two hypocentral depth ranges, 0–25 and 25–70 km; 158 tectonic sources; appropriate seismicity modeling; 14 ground motion prediction equations for three seismotectonic provinces, viz. the east-central Himalaya, the Bengal Basin and Northeast India selected through suitability testing; and appropriate weighting in a logic tree framework. Site classification of Kolkata performed following in-depth geophysical and geotechnical investigations places the city in D1, D2, D3 and E classes. Probabilistic seismic hazard assessment at a surface-consistent level – i.e., the local seismic hazard related to site amplification performed by propagating the bedrock ground motion with 10% probability of exceedance in 50 years through a 1-D sediment column using an equivalent linear analysis – predicts a peak ground acceleration (PGA) range from 0.176 to 0.253 g in the city. A deterministic liquefaction scenario in terms of spatial distribution of liquefaction potential index corresponding to surface PGA distribution places 50% of the city in the possible liquefiable zone. A multicriteria seismic hazard microzonation framework is proposed for judicious integration of multiple themes, namely PGA at the surface, liquefaction potential index, NEHRP soil site class, sediment class, geomorphology and ground water table in a fuzzy protocol in the geographical information system by adopting an analytical hierarchal process. The resulting high-resolution surface consistent hazard, liquefaction and microzonation maps are expected to play vital roles in earthquake-related disaster mitigation and management of the city of Kolkata.