Global landslide and avalanche hotspots

Allocating resources for natural hazard risk management has high priority in development banks and international agencies working in developing countries. Global hazard and risk maps for landslides and avalanches were developed to identify the most exposed countries. Based on the global datasets of climate, lithology, earthquake activity, and topography, areas with the highest hazard, or “hotspots”, were identified. The applied model was based on classed values of all input data. The model output is a landslide and avalanche hazard index, which is globally scaled into nine levels. The model results were calibrated and validated in selected areas where good data on slide events exist. The results from the landslide and avalanche hazard model together with global population data were then used as input for the risk assessment. Regions with the highest risk can be found in Colombia, Tajikistan, India, and Nepal where the estimated number of people killed per year per 100 km² was found to be greater than one. The model made a reasonable prediction of the landslide hazard in 240 of 249 countries. More and better input data could improve the model further. Future work will focus on selected areas to study the applicability of the model on national and regional scales.     

First-order regional seismotectonic model for South Africa

A first-order seismotectonic model was created for South Africa. This was done using four logical steps: geoscientific data collection, characterisation, assimilation and zonation. Through the definition of subunits of concentrations of earthquake foci and large neotectonic and structural domains, seismotectonic structures, systems and domains were created. Relatively larger controls of seismicity exist between the Great Escarpment and the coast. In the south, this region is characterised by large aeromagnetic anomalies and large EW trending faults. In the west, it is characterised by the NW–SE trending Wegener stress anomaly, radial-trending dykes and earthquake clusters. In the east, it is characterised by a large neotectonic domain where several large historical earthquakes occurred. In the centre of South Africa, several clusters of earthquake activity are found, often related to mining activity. Further north, seismicity is related to both mining activity and neotectonic deformation. This work contributes to the development of a seismotectonic model for South Africa by (1) bringing together, digitally, several data sets in a common GIS platform (geology, geophysics, stress, seismicity, neotectonics, topography, crustal and mantle structure and anisotropy), (2) understanding the significance of data sets for seismotectonic zonation and limitations thereof and (3) obtaining a reasonable regional model for use in seismic hazard assessments.     

Probabilistic seismic hazard assessment in the northeastern part of Algeria

This work involves updating the evaluation of seismic hazard in Northeast Algeria by a probabilistic approach. This reassessment attempts to resolve inconsistencies between seismic zoning in regional building codes and is further motivated by the need to refine the input data that are used to evaluate seismic hazard scenarios. We adopted a seismotectonic model that accounts for differences in interpretations of regional seismicity. We then performed a probabilistic assessment of regional seismic hazard in Northeast Algeria. Based on a homogeneous earthquake catalog and geological and seismotectonic data gathered in the first part of the study, a seismotectonic zoning map was created and seven risk areas were identified. For each area, peak ground acceleration hazard maps were produced. Details of the calculations are provided, including hazard curves at periods of 0.1, 0.2, 0.33, 0.5, 1.0, and 2.0 s and uniform hazard spectra at urban locations in the area, including Sétif, Constantine, Kherrata, Bejaia, and Jijel.     

东南亚地震活动与b 值的时空分布特征

东南亚地区是“21世纪海上丝绸之路”（以下简称“海洋丝路”）的重要组成部分,该区历史上曾发生十余次巨大地震,地震及其次生地质灾害是威胁东南亚地区经济社会发展和国际合作的主要自然灾害。系统梳理该区地震活动的时空分布特征及评估未来灾害风险格局,对于推进“一带一路”倡议实施及区域经济社会可持续发展具有重要意义。文章基于东南亚地区1900年以来M≥5地震的时空分布统计分析和地震b值计算,揭示出该区的地震活动在时间上表现出活跃期与平静期交替变化的特征;空间上表现出明显的聚集效应,成丛性强且主要集中在5个地震统计区内,其中印尼—马来多岛弧盆系地震区和菲律宾群岛地震区的地震活动最为活跃。总体而言,东南亚5个地震区的b值偏低,在0.42~0.91之间。该区内的地震b值也存在时空差异,受大地震事件、俯冲带年龄、活动断裂带和震源深度等众多因素影响,但主控因素在不同区域有所不同。地震b值时空变化特征对区域地震活动预测具有启示作用。上述认识为推进“海洋丝路”工程建设和“一带一路”防灾减灾对策提供了科学支撑。     

Seismicity and seismotectonics of the Gulf of Aqaba region

The historical seismicity of the last ten centuries and the instrumental data that occurred in the Gulf of Aqaba region during the period 1982–2008 are evaluated. It is found that 12 historical earthquakes have occurred with average recurrence periods 70–90 and 333–500 years for M?≥?6.0 and 7.0, respectively. Those with M?≤?6.5 appear to be incomplete and require further investigation. More than 98 % of the instrumental data has occurred in the form of swarms and sequences. The first have released about 32 % of the total energy and are most likely related to subsurface volcanic activities. Their epicentral distribution indicates that all regional faults of the gulf area are active in the present, but with clear concentration within the area bound by latitudes 28.2°–29.8° and longitudes 34.4°–35.2°. Regional strike-slip faults of the northern two basins appear to be as twice active as the normal, or more. An appreciable level of seismic hazard is envisaged as the “a” value is 6.0–6.2 while the “b” value shows a temporal variation, mostly in the range 0.8–1.05. More than 95 % of the seismic energy was released from earthquakes shallower than 22 km. This indicates a brittle upper crust and a ductile lower crust and upper mantle. Tectonic movements at depths?>?22 km appear to be aseismic. The epicentral distribution of the five swarms indicates that the lengths of the causative faults varied in the range 45–70 km. The maximum expected magnitude is Mw?=?6.8–7.2. This implies a seismic slip rate of about 0.54–0.8 Cm/year and some 20–30 % of aseismic tectonic movements. This and the sequence nature of the seismicity of this region result in a noticeable hazard reduction. Combining the seismicity data of the Gulf of Aqaba region with other geophysical, geological, tectonic, and environmental data, clearly indicate that the seismicity of this region is as old as the initiation of the gulf itself. No apparent southward or northward migration of activity is observed.     

Detailed seismicity mapping of the Altai-Sayan zone using large averaging areas

An improved technique is suggested for quantifying seismic activity over averaging areas of an arbitrary size. The example of the Altai-Sayan seismic zone is used to substantiate the choice of a 1° N × 2° E averaging area instead of the traditional one of 40×40 km². Maps compiled with averaging areas of different sizes can be spliced and correlated using a correction coefficient estimated in different models. The new seismicity map of the Altai-Sayan area covers more than 90% of the territory and provides a generalized image of activity being advantageous over the classic maps as it allows better correlation of regional seismicity with the tectonic setting. With larger averaging areas and, correspondingly, a greater amount of data in each area, one can obtain time series of seismic activity to be analyzed using mathematical statistics as a basis for mathematical modeling and simulation of the seismic process.     

Seismic Hazards Assessment for Radioactive Waste Disposal Sites in Regions of Low Seismic Activity

A comparative tectonic quiescence and lack of earthquakes make the stable centres of continents attractive for siting long-term radioactive waste storage facilities. The low rates of deformation in such regions, however, make it difficult to characterize their long-term seismotectonic behaviour, leading to uncertain estimates for the very low probability hazard estimates required by society. In an attempt to overcome the deficiency of both contemporary seismicity and paleoseismic data in central Canada, we have used earthquake histories from regions with similar seismotectonic characteristics from around the world. Substituting space for time, we estimate a long-term rate per 10⁶ km² of 0.004 magnitude ≥6 earthquakes per annum, of which 33–100% might rupture to the surface.     

Seismic hazard in Northern Algeria using spatially smoothed seismicity. Results for peak ground acceleration

Seismic hazard in terms of peak ground acceleration (PGA) has been evaluated in northern Algeria using spatially smoothed seismicity data. We present here a preliminary seismic zoning in northern Algeria as derived from the obtained results.Initially, we have compiled an earthquake catalog of the region taking data from several agencies. Afterwards, we have delimited seismic areas where the b and m_max parameters are different. Finally, by applying the methodology proposed by Frankel [Seismol. Res. Lett. 66 (1995) 8], and using four complete and Poissonian seismicity models, we are able to compute the seismic hazard maps in terms of PGA with 39.3% and 10% probability of exceedance in 50 years.A significant result of this work is the observation of mean PGA values of the order of 0.20 and 0.45 g, for return periods of 100 and 475 years, respectively, in the central area of the Tell Atlas.     

On the swelling potential of compacted high plasticity clays

We conducted a study of the spatial distributions of seismicity and earthquake hazard parameters for Turkey and the adjacent areas, applying the maximum likelihood method. The procedure allows for the use of either historical or instrumental data, or even a combination of the two. By using this method, we can estimate the earthquake hazard parameters, which include the maximum regional magnitude Mˆ_max, the activity rate of seismic events and the well-known bˆ value, which is the slope of the frequency-magnitude Gutenberg-Richter relationship. These three parameters are determined simultaneously using an iterative scheme. The uncertainty in the determination of the magnitudes was also taken into consideration. The return periods (RP) of earthquakes with a magnitude M ≥ m are also evaluated. The whole examined area is divided into 24 seismic regions based on their seismotectonic regime. The homogeneity of the magnitudes is an essential factor in such studies. In order to achieve homogeneity of the magnitudes, formulas that convert any magnitude to an M_S-surface scale are developed. New completeness cutoffs and their corresponding time intervals are also assessed for each of the 24 seismic regions. Each of the obtained parameters is distributed into its respective seismic region, allowing for an analysis of the localized seismicity parameters and a representation of their regional variation on a map. The earthquake hazard level is also calculated as a function of the form Θ = (Mˆ_max,RP_6.0), and a relative hazard scale (defined as the index K) is defined for each seismic region. The investigated regions are then classified into five groups using these parameters. This classification is useful for theoretical and practical reasons and provides a picture of quantitative seismicity. An attempt is then made to relate these values to the local tectonics.     

A new probabilistic seismic hazard study of Spain

In this paper, we present a probabilistic seismic hazard analysis (PSHA) for mainland Spain that takes into account recent new results in seismicity, seismic zoning, and strong ground attenuation not considered in the latest PSHA of the Spanish Building Code. Those new input data have been obtained as a three-step project carried out in order to improve the existing hazard map for mainland Spain. We have produced a new earthquake catalogue for the area, in which the earthquakes are given in moment magnitude through specific deduced relationships for our territory based on intensity data (Mezcua et al. in Seismol Res Lett 75:75–81, 2004). In addition, we included a new seismogenetic zoning based on the recent partial zoning studies performed by different authors. Finally, as we have developed a new strong ground motion model for the area García Blanco (2009), it was considered in the hazard calculation together with other attenuations gathered from different authors using data compatible with our region. With this new data, a logic tree process is defined to quantify the epistemic uncertainty related to those parts of the process. A sensitivity test has been included in order to analyze the different models of ground motion and seismotectonic zonation used in this work. Finally, after applying a weighting scheme, a mean hazard map for PGA, based on rock type condition for 10% exceedance probability in 50 years, is presented, including 15th and 85th percentile hazard maps. The main differences with the present official building code hazard map are analyzed.     

Multi-scale imaging of a slow active fault zone: contribution for improved seismic hazard assessment in the Swiss Alpine foreland

Seismic hazard assessment of slow active fault zones is challenging as usually only a few decades of sparse instrumental seismic monitoring is available to characterize seismic activity. Tectonic features linked to the observed seismicity can be mapped by seismic imaging techniques and/or geomorphological and structural evidences. In this study, we investigate a seismic lineament located in the Swiss Alpine foreland, which was discussed in previous work as being related to crustal structures carrying in size the potential of a magnitude M 6 earthquake. New, low-magnitude (?2.0 ≤ M_L ≤ 2.5) earthquake data are used to image the spatial and temporal distribution of seismogenic features in the target area. Quantitative and qualitative analyses are applied to the waveform dataset to better constrain earthquakes distribution and source processes. Potential tectonic features responsible for the observed seismicity are modelled based on new reinterpretations of oil industry seismic profiles and recent field data in the study area. The earthquake and tectonic datasets are then integrated in a 3D model. Spatially, the seismicity correlates over 10–15 km with a N–S oriented sub-vertical fault zone imaged in seismic profiles in the Mesozoic cover units above a major decollement on top of the mechanically more rigid basement and seen in outcrops of Tertiary series east of the city of Fribourg. Observed earthquakes cluster at shallow depth (<4 km) in the sedimentary cover. Given the spatial extend of the observed seismicity, we infer the potential of a moderate size earthquake to be generated on the lineament. However, since the existence of along strike structures in the basement cannot be excluded, a maximum M 6 earthquake cannot be ruled out. Thus, the Fribourg Lineament constitutes a non-negligible source of seismic hazard in the Swiss Alpine foreland.     

The instrumental seismicity of the Jordan Dead Sea transform

The instrumental seismicity that occurred in the Jordan Dead Sea transform region during the period 1900–2014 is compiled from all available sources. Some 492 phosphate mining explosions (M ≤ 3.9) are recognized and filtered from the data. Excluding these, it is found that 4448 earthquakes have occurred with magnitudes M ≥ 3.0. Only 572, 18 and 2 of these had magnitudes M ≥ 4, 5, and 6 in respective order. Average recurrence periods for the 5 and 6 magnitudes are 6.3 and 57 years. Much of these have occurred in sequences and swarms. The epicentral distribution of the compiled instrumental seismicity data shows very good correlation with the general tectonics of the study region. All tectonic elements are active in the present with a noticeable hazard. The regional strike-slip faults of the transform proper remain the major sources of this hazard. They account for not less than 99% of the seismic energy released from all instrumental data. The calculated a-parameter of the whole transform is 6.6. It varies for all its strike-slip faults mostly in the range 6.0–6.6. The b-value of the whole transform and some of its major segments is 1.0. Others show b-variations in the range of 1.1–1.3. Such a- and b-values imply recurrence periods of 38 years and 395 years for the 6 and 7 magnitude earthquakes. Such values, their variations and the seismic moment calculations clearly indicate an appreciable level of seismic hazard associated with all segments. This hazard appears to be highest for Al Ghab segment, followed by Beqa’a and Wadi Araba segments, respectively. The other three segments appear to be of lower hazard. The seismicity of this region is very shallow. More than 99% of the seismic energy has been released from the brittle granitic upper crust whose thickness is about 21 km and its Poisson’s ratio is 0.25. More than 93.6% of the energy was released from its upper 10 km. Very little energy is released from the underlying ductile basaltic crust whose Poisson’s ratio is 0.29. The calculated seismic slip rate along the Whole Jordan Dead Sea transform is 0.54 cm/year if the fault depth is assumed 10 km. It increases to 0.77 and 1.07 cm/year if the fault depth is reduced to 7 and 5 km, respectively. These slip rates are comparable with the long-term geologically deduced rate of 1 cm/year.     

Giant sandstone pipes record basin-scale liquefaction of buried dune sands in the Middle Jurassic of SE Utah

Giant sandstone pipes crop out over a 20 000 km² area of SE Utah. They range from metres to decametres in diameter and crosscut more than 100 m of continental sediments. New and published observations, encompassing all known pipe localities in SE Utah, suggest that the pipes formed by upward movement of groundwater, sand and sediment blocks in response to liquefaction of buried dune sands. Pipe formation was most likely triggered by seismicity during a phase of Cordilleran shortening and regional tectonic reorganization in the late Middle Jurassic. The depth of liquefaction recorded by the pipes far exceeds depths usually considered in earthquake engineering.     

Uncertainties in the estimation of <Emphasis Type="Italic">M</Emphasis><Subscript>max</Subscript>

In the present paper, the parameters affecting the uncertainties on the estimation of M _max have been investigated by exploring different methodologies being used in the analysis of seismicity catalogue and estimation of seismicity parameters. A critical issue to be addressed before any scientific analysis is to assess the quality, consistency, and homogeneity of the data. The empirical relationships between different magnitude scales have been used for conversions for homogenization of seismicity catalogues to be used for further seismic hazard assessment studies. An endeavour has been made to quantify the uncertainties due to magnitude conversions and the seismic hazard parameters are then estimated using different methods to consider the epistemic uncertainty in the process. The study area chosen is around Delhi. The b value and the magnitude of completeness for the four seismogenic sources considered around Delhi varied more than 40% using the three catalogues compiled based on different magnitude conversion relationships. The effect of the uncertainties has been then shown on the estimation of M _max and the probabilities of occurrence of different magnitudes. It has been emphasized to consider the uncertainties and their quantification to carry out seismic hazard assessment and in turn the seismic microzonation.     

Geochemical hazard evaluation of sulphide-rich iron mines: The Rio Marina district (Elba Island,Italy)

Rio Marina mining district (Elba Island) is characterised by hematite + pyrite ore association and was exploited for iron till 1981, leaving waste rock dumps of several millions m³. The effect of open pit mining activity in this site is to produce acid mine drainage (AMD) processes leading to environmental pollution, testified by all the sampled waters (Giove stream, drainage channels, superficial pools and settling basin) which have pH values ranging from 2.08 to 3.35 and heavy metal concentrations that reach 903.16 mg/l for Fe, 45.02 mg/l for Mn, 10.08 mg/l for Zn and 1.75 mg/l for Cu. In the present work a space and time related approach to geochemical hazard evaluation was applied. The geochemical hazard is mainly related to high heavy metal concentration, acid mine drainage processes development and topographic setting. As all these parameters are related in space, hazard evaluation was performed by geostatistical methods. Fifty-four earth material samples (residual soils, waste rocks or debris materials) were collected in a central aligned 100 m mesh square grid. These were analysed for major elements by XRF, for Cu, Pb, Zn by ICP-AES and for AMD potential following the AMIRA procedure. The concentration of heavy metals was compared with Italian law limits. The overlap of Cu, Pb and Zn content maps show that at least one of these heavy metals exceed law limits in all the area. The AMD test results show that more than 50% of samples have a positive NAPP (Net Acid Producing Potential) that could reach 258.9 kg H₂SO₄/t. According to the obtained data, three main geochemical hazard classes were established and their distribution in the mining area was assessed. About 51% of the mining area surface belongs to the major hazard class, where AMD process occurs, about 49% belongs to an intermediate hazard class, where AMD process could occur only if certain conditions are met. Finally, the persistence of the AMD process in the Rio Marina area was evaluated on the basis of yearly rainfall, mining waters pH and NAPP values. A complete leaching of the first 0.25 m of the earth materials can retain the current environmental conditions for several centuries.     

Improved earthquake locations in the Koyna-Warna seismic zone

In estimating the likelihood of an earthquake hazard for a seismically active region, information on the geometry of the potential source is important in quantifying the seismic hazard. The damage from an earthquake varies spatially and is governed by the fault geometry and lithology. As earthquake damage is amplified by guided seismic waves along fault zones, it is important to delineate the disposition of the fault zones by precisely determined hypocentral parameters. We used the double difference (DD) algorithm to relocate earthquakes in the Koyna-Warna seismic zone (KWSZ) region, with the P- and S-wave catalog data from relative arrival time pairs constituting the input. A significant improvement in the hypocentral estimates was achieved, with the epicentral errors <30 m and focal depth errors <75 m i.e. errors have been significantly reduced by an order of magnitude from the parameters determined by HYPO71. The earthquake activity defines three different fault segments. The seismogenic volume is shallower in the south by 3 km, with seismicity in the north extending to a depth of 11 km while in the south the deepest seismicity observed is at a depth of 8 km. By resolving the structure of seismicity in greater detail, we address the salient issues related to the seismotectonics of this region.     

Indicative Structures of Paleoseismicity in the Acequion River Valley,San Juan Province,Central-Western Argentina

Evidence for paleo-seismicity has been discovered in the Acequión river valley, in West central Argentina. Two Holocene rock avalanches have been observed; the most recent of these dammed a lake, whose sediments contain liquefaction structures. At least five paleo-earthquakes affected this region during the late Quaternary, as deduced from the succession of their secondary effects. The magnitude and the probable tectonic source of these paleo-events are discussed. The observed liquefaction features associated with slumps, joints, fractures, and faults, should be generated by M>5 earthquakes related to the nearby quaternary Cerro Salinas fault, which belongs to the Eastern Precordillera fault system. These data extend the regional seismicity record to the Holocene and highlight the high seismic hazard in this part of Argentina.     

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.     

Seismogenesis of clustered seismicity beneath the Kangra–Chamba sector of northwest Himalaya: Constraints from 3D local earthquake tomography

To investigate subsurface structure and seismogenic layers, 3D velocity inversion was carried out in the source zone of 1905 Kangra earthquake (M8.0) in the northwestern Himalaya. P-wave and S-wave phase data of 159 earthquakes recorded by a network of 21 stations were used for this purpose. Inverted velocity tomograms up to a depth range of 18 km show significant variations of 14% in Vp and Vs and 6% in the Vp/Vs across the major tectonic zones in the region. Synthesis of seismicity pattern, velocity structure, distinctive focal mechanisms coupled with nature of stress distribution allows mapping of three different source regions that control regional seismotectonics. Accumulating strains are partly consumed by sliding of Chamba Nappe to the southwest through reverse-fault movements along Chamba/Panjal/Main Boundary Thrusts. This coupled with normal-fault type displacements along Chenab Normal Fault in the north account for low magnitude widespread seismicity in upper 8–10 km of the crust. At intermediate depths from 8 to 15 km, adjusting to residual compressive stresses, the detachment or lower end of the MBT slips to produce thrust dominated seismicity. Nucleation of secondary stresses in local NE–SW oriented structure interacts in complex manner with regional stresses to generate normal type earthquakes below the plane of detachment and therefore three seismic regimes at different depths produce intense seismicity in a block of 30 × 30 km² centered NE to the epicenter of Kangra earthquake.     

Recent seismicity and crustal stress field in the Lucanian Apennines and surrounding areas (Southern Italy): Seismotectonic implications

We analyzed the instrumental seismicity of Southern Italy in the area including the Lucanian Apennines and Bradano foredeep, making use of the most recent seismological data base available so far. P- and S-wave arrival times, recorded by the Italian National Seismic Network (RSNC) operated by the Istituto Nazionale di Geofisica e Vulcanologia (INGV), were re-picked along with those of the SAPTEX temporary array deployed in the region in the period 2001–2004. For some events located in the upper Val d'Agri, we also used data from the Eni-Agip oil company seismic network. We examined the seismicity occurred during the period between 2001 and 2006, considering 514 events with magnitudes M ≥ 2.0. We computed the V_P/V_S ratio obtaining a value of 1.83 and we carried out an analysis for the one-dimensional (1D) velocity model that approximates the seismic structure of the study area. Earthquakes were relocated and, for well- recorded events, we also computed 108 fault plane solutions. Finally, using 58 solutions, the most constrained, we computed regional stress field in the study area.Earthquake distribution shows three main seismic regions: the westernmost (Lucanian Apennines) characterized by high background seismicity, mostly with shallow hypocenters, the easternmost below the Bradano foredeep and the Murge with deeper and more scattered seismicity, and finally the more isolated and sparse seismicity localized in the Sila Range and in the offshore area along the northeastern Calabrian coast. Focal mechanisms computed in this work are in large part normal and strike-slip solutions and their tensional axes (T-axes) have a generalized NE–SW orientation. The denser station coverage allowed us to improve hypocenters determination compared to those obtained by using only RSNC data, for a better characterization of the crustal and subcrustal seismicity in the study area.