Geological and morphological study of the Jiufengershan landslide triggered by the Chi-Chi Taiwan earthquake

The Jiufengershan rock and soil avalanche is one of the largest landslides triggered by the Chi-Chi earthquake Taiwan 1999. The landslide destabilized the western limb of the Taanshan syncline along a weak stratigraphic layer. It involved a flatiron remnant, which was almost entirely mobilized during the earthquake. The avalanche was slowed down by NS trending ridges located downstream along the Jiutsaihu creek. The landslide affected a 60 m thick and 1.5 km long sedimentary pile composed of shales and sandstones, which dip 22°SE toward a transverse valley. The triggering mechanism and the sliding process were analyzed by means of geological and morphological data from aerial photographs and observed in the field. A high-resolution airborne Light Detection and Ranging (LiDAR) image taken 2.5 years after the landslide allows the identification of morphological structures along the sliding surface and the landslide accumulation. The sliding surface shows several deformation structures such as fault scarps and folds. These structures are interpreted in terms of basal shear stresses created during the avalanche. Three major joint sets were identified at the sliding surface. The isopach map of the landslide was calculated from the comparison between elevation models before and after the earthquake. The coseismic volume of mobilized material and landslide deposit data are 42 × 10⁶ m³ and 50 × 10⁶ m³, respectively. The geometry of the landslide accumulation in the field has an irregular star shape. The morphology of the deposit area shows a sequence of smooth reliefs and depressions that contrast with the neighboring ridges.     

Tsunami as agents of geomorphic change in mid-ocean reef islands

Low-lying atoll islands appear highly vulnerable to the effects of climate change and extreme natural events. Potentially disastrous effects of future sea-level rise have been inferred in many studies, and the actual impacts of tropical storms on island destruction and formation have been well documented. In contrast, the role of tsunami in the geomorphic development of atoll islands has not been investigated. The Sumatran earthquake of 26 December 2004 generated a tsunami that reached the Maldives 2500 km away, with waves up to 2.5 m high. Observations on the geomorphic changes resulting from the tsunami are detailed here, based on pre-and post-tsunami profile measurements of island, beach and reef topography, and GPS surveys of the planform shape of islands and beaches of 11 uninhabited islands in South Maalhosmadulu atoll, Maldives. Erosional and depositional impacts were observed on all islands and these have been quantified. In general the changes were of a minor nature with a maximum reduction in island area of 9% and average of 3.75%. Rather, the tsunami accentuated predictable seasonal oscillations in shoreline change, including localised erosion reflected in fresh scarps and seepage gullies. Depositional features in the form of sand sheets and sand lobes emplaced on the vegetated island surfaces provide clear evidence that the tsunami waves washed over parts of all the islands. Both erosional scarps and overwash deposits were concentrated at the tsunami-exposed eastern sides of the islands. Impacts on leeward shores were primarily accretionary, in the form of spit and cuspate foreland extension. Whereas the nature and magnitude of intra-and inter-island impacts was variable, an east to west decline in aggregate effects was noted. Detailed consideration of the morphodynamic interaction between the tsunami waves and island morphology, show that this cross-atoll gradient resulted not just from the reduction in tsunami energy as it passed through the atoll, but also from variations in elevation of the encircling island ridge, and the quantity and distribution of sediment in the antecedent beach. A conceptual model identifying the sequence of changes to individual islands supports the observational data and the pattern of geomorphic changes resulting from the tsunami. This model leads to consideration of the longer-term impacts of the tsunami on the future stability of islands. Four scenarios are presented, each of which has a different island-beach sediment budget, and different relaxation time to achieve dynamic equilibrium.     

Remote sensing-based assessment of tsunami vulnerability and risk in Alexandria, Egypt

Tsunamis can cause catastrophic loss of life, destruction of property, engineered structures and coastal infrastructure, and they can lead to major economic losses. Even though tsunamis are relatively rare in the Mediterranean Sea, their potential danger to cities along the Mediterranean coast cannot be neglected. In order to create awareness among the potentially affected people it is important to know the risk and vulnerability of the population and infrastructure related to a possible tsunami impact. In this work a hazard, vulnerability and risk analysis for buildings in two districts of Alexandria was carried out. Relevant input parameters were derived mainly from remote sensing and field data and were analyzed with a geographical information system (GIS). Based on historical records of past tsunamis, two inundation scenarios of 5 m and 9 m were defined and modeled applying a bath-type model. The resulting tsunami building risk zone maps showed that 12% of the buildings in El Gomrok district are at high or very high risk for the 5 m scenario, while the risk for El Montazah area is low. For the 9 m scenario, on the other hand, the majority of the buildings in both districts, 56% of El Gomrok, and 60% of El Montazah, are in the high or very high risk zone. An analysis of the building use indicated that the majority of these buildings are residential and commercial types, highlighting that the potential consequences of a tsunami could be severe. Due to the scarcity of historical data no frequency could be associated with the two selected scenarios. While both are credible we consider the 5 m scenario as possible but unlikely and the 9 m scenario as unlikely.     

The role of strong earthquakes and tsunami in the Late Holocene evolution of the Fortore River coastal plain (Apulia, Italy): A synthesis

Morphological analysis of the Fortore River coastal plain and the Lesina Lake coastal barrier integrated with radiocarbon age data indicates that the evolution of the coastal landscape has been strongly affected by a number of strong earthquakes and related tsunamis which occurred during the last 3000 years. The first seismic event struck this coastal area in the V century BC. It produced strong erosion of the Fortore River coastal plain and significant emersion of Punta delle Pietre Nere, as well as the large tsunami responsible for the development of the Sant'Andrea washover fan. The second event occurred in 493 AD; it induced severe erosion of the Fortore River coastal plain and triggered the large tsunami that hit the Lesina Lake coastal barrier, producing the Foce Cauto washover fan. Then later in 1627, an earthquake was responsible for the further coseismic uplift of Punta delle Pietre Nere, the subsidence of Lesina village area and the development of a tsunami which produced two washover fans.Morphological analysis points out that seismic events strong enough to control the morphological evolution of local coastal landscapes show a statistical return period of about 1000 years. These major events produced important coseismic vertical movements and large tsunamis. However, the correct identification of the tectonic structure responsible for the generation of these strong earthquakes is still an unsolved problem.     

汶川5·12地震次生泥石流沟应急判识方法与指标

5·12汶川大地震诱发了崩塌、滑坡、泥石流等次生灾害,崩塌、滑坡堆积物给泥石流的形成提供了大量松散固体物质,将导致灾区部分山洪沟转化为泥石流沟,为此,给出了一种泥石流沟的判识方法和指标.调查发现,汶川灾区的地形地貌和降雨条件满足泥石流的暴发条件,提出用流域单位面积的松散固体物质方量来判识泥石流沟;调查西部山区的50条泥石流沟,提出以0.1 m3/m2的松散固体物质量作为泥石流沟的判别指标,以2m3/m2的松散固体物质量作为粘性泥石流沟的判别指标.     

Uncertainty due to DEM error in landslide susceptibility mapping

Terrain attributes such as slope gradient and slope shape, computed from a gridded digital elevation model (DEM), are important input data for landslide susceptibility mapping. Errors in DEM can cause uncertainty in terrain attributes and thus influence landslide susceptibility mapping. Monte Carlo simulations have been used in this article to compare uncertainties due to DEM error in two representative landslide susceptibility mapping approaches: a recently developed expert knowledge and fuzzy logic-based approach to landslide susceptibility mapping (efLandslides), and a logistic regression approach that is representative of multivariate statistical approaches to landslide susceptibility mapping. The study area is located in the middle and upper reaches of the Yangtze River, China, and includes two adjacent areas with similar environmental conditions – one for efLandslides model development (approximately 250 km²) and the other for model extrapolation (approximately 4600 km²). Sequential Gaussian simulation was used to simulate DEM error fields at 25-m resolution with different magnitudes and spatial autocorrelation levels. Nine sets of simulations were generated. Each set included 100 realizations derived from a DEM error field specified by possible combinations of three standard deviation values (1, 7.5, and 15 m) for error magnitude and three range values (0, 60, and 120 m) for spatial autocorrelation. The overall uncertainties of both efLandslides and the logistic regression approach attributable to each model-simulated DEM error were evaluated based on a map of standard deviations of landslide susceptibility realizations. The uncertainty assessment showed that the overall uncertainty in efLandslides was less sensitive to DEM error than that in the logistic regression approach and that the overall uncertainties in both efLandslides and the logistic regression approach for the model-extrapolation area were generally lower than in the model-development area used in this study. Boxplots were produced by associating an independent validation set of 205 observed landslides in the model-extrapolation area with the resulting landslide susceptibility realizations. These boxplots showed that for all simulations, efLandslides produced more reasonable results than logistic regression.     

The Elusive AD 1826 Tsunami, South Westland, New Zealand

In AD 1826 sealers reported earthquake and tsunami activity in Fiordland, although, contemporary or near‐contemporary accounts of tsunami inundation at the time are elusive. A detailed analysis of recent sediments from Okarito Lagoon builds on contextual evidence provided by earlier research concerning past tsunami inundation. Sedimentological, geochemical, micropalaeontological and geochronological data are used to determine palaeoenvironments before, during and after what was most probably tsunami inundation in AD 1826. The most compelling chronological control is provided fry a young cohort of trees growing on a raised shoreline bench stranded fry a drop in the lagoon water level following tsunami inundation.     

Tsunamigenic incisions produced by the December 2004 earthquake along the coasts of Thailand, Indonesia and Sri Lanka

Field observations and satellite images indicate that tsunami waves exhibit specific patterns during flooding and recession forming characteristic incisions in the coastal landscape. To study these incisions we analyze high resolution remote sensing images of the coastline of Indonesia, Thailand and Sri Lanka impacted by the tsunami of December 26th, 2004. The analysis sheds light on the different mechanisms by which currents scour incisions during the flooding and receding phases of a tsunami. During flooding the high velocity flow indents the levees of existing tidal channels and bays, leaving short flood scours. The receding water then dissects the coastline with equally spaced return channels widening toward the coast.     

Local and regional components of western Aegean deformation extracted from 100 years of geodetic displacement measurements

Our objectives are as follows. First, we wish to develop a methodology to recover the long-term component of deformation from any set of distributed, time-averaged geodetic strain measurements that were subject to seismic disturbance, given a catalogue of local seismicity that occurred during the measurement period. Second, using seismic and geodetic data sets that span approximately 100 years, we apply this technique in the western Aegean to assess the role of local seismicity in regional deformation. The methodology is developed using a model for crustal deformation constructed from a long-term, smooth regional strain field combined with instantaneous, local perturbations from upper-crustal earthquakes approximated by static elastic dislocations. By inverting geodetic displacements for the smooth field while simultaneously floating influential but uncertain earthquake source parameters, an estimate of the regional component of deformation that is approximately independent of the seismicity can be made. In the western Aegean we find that the horizontal component of regional deformation can be described with minor inaccuracy by a quadratic relative displacement field. The principal horizontal extensional axes calculated from the regionally smooth displacement field agree in orientation with the T-axes of earthquakes in the region. These observations indicate that the instantaneous elastic strain of the 10 km thick seismogenic layer is driven by a stress field that is smooth on the scale of the geodetic network as a whole, 200-300 km.     

Analysing seismic-source mechanisms by linear-programming methods

Summary. Linear-programming methods are powerful and efficient tools for objectively analysing seismic focal mechanisms and are applicable to a wide range of problems, including tsunami warning and nuclear explosion identification. The source mechanism is represented as a point in the six-dimensional space of moment-tensor components. Each observed polarity provides an inequality constraint, linear with respect to the moment tensor components, that restricts the solution to a half-space bounded by a hyperplane passing through the origin. The intersection of these half-spaces is the convex set of all acceptable solutions. Using linear programming, a solution consistent with the polarity constraints can be obtained that maximizes or minimizes any desired linear function of the moment tensor components; the dilatation, the thrust-like nature, and the strike-slip-like nature of an event are examples of such functions. The present method can easily be extended to fit observed seismic-wave amplitudes (either signed or absolute) subject to polarity constraints, and to assess the range of mechanisms consistent with a set of measured amplitudes.     

Landslide susceptibility mapping using geological data, a DEM from ASTER images and an Artificial Neural Network (ANN)

An efficient and accurate method of generating landslide susceptibility maps is very important to mitigate the loss of properties and lives caused by this type of geological hazard. This study focuses on the development of an accurate and efficient method of data integration, processing and generation of a landslide susceptibility map using an ANN and data from ASTER images. The method contains two major phases. The first phase is the data integration and analysis, and the second is the Artificial Neural Network training and mapping. The data integration and analysis phase involve GIS based statistical analysis relating landslide occurrence to geological and DEM (digital elevation model) derived geomorphological parameters. The parameters include slope, aspect, elevation, geology, density of geological boundaries and distance to the boundaries. This phase determines the geological and geomorphological factors that are significantly correlated with landslide occurrence. The second phase further relates the landslide susceptibility index to the important geological and geomorphological parameters identified in the first phase through ANN training. The trained ANN is then used to generate a landslide susceptibility map. Landslide data from the 2004 Niigata earthquake and a DEM derived from ASTER images were used. The area provided enough landslide data to check the efficiency and accuracy of the developed method. Based on the initial results of the experiment, the developed method is more than 90% accurate in determining the probability of landslide occurrence in a particular area.     

Neural dynamic modelling on earthquake magnitude series

Earthquake magnitude prediction is of vital importance for human safety. The earthquake is a very complicated and non-linear dynamic process. It cannot be described adequately by any deterministic models. In this paper a neural dynamic modelling for earthquake magnitude prediction is reported. Historical records of earthquake magnitude series are used to construct the optimal non-linear dynamic model, and the consequent outcome of the earthquake behaviour is then predicted by this model. In turn, the latest recorded data set can be fed back to improve the accuracy of the neural dynamic model. The modelling of experiments of three earthquake magnitude series in China and Japan and their extrapolated predictions are included in this paper. The values predicted by extrapolation are in good agreement with the historical data.     

A numerical method to determine physical source parameters from free oscillation data

A method to determine physical source parameter using free oscillation data is presented. It is assumed that the geometry of the source is known, e.g. from P -wave data. The source is assumed to propagate in the horizontal direction, while unknown parameters to be determined are the azimuth and velocity of propagation, the distance over which the seismic source propagated and the source intensity as a function of propagation distance.
The method consists in the systematic search for the set of source parameters rendering phase corrections which maximize the spectral peak amplitudes within the excitation criterion scheme.
If there is no precursive motion, the average dislocation time function can be determined from the spectrum of the seismic moment and the space source intensity. The source intensity as a function of instantaneous source location is found independently of the P -wave origin time and source dislocation time function. The method does not require to correct the data for attenuation.     

Fault slip in the 1997 Manyi, Tibet earthquake from linear elastic modelling of InSAR displacements

The M _w 7.6 1997 Manyi earthquake occurred in an area of central northern Tibet where sparse vegetation coverage and a lack of human habitation provide excellent conditions for Interferometric Synthetic Aperture Radar (InSAR) studies. We use coseismic pairs of radar images acquired by the ESA ERS-2 satellite to construct interferograms of the surface displacement field due to the earthquake. The location and extent of the coseismic fault rupture are mapped using a combination of optical satellite imagery, high-resolution digital topography, interferometric correlation and azimuth offset measurements; in so doing, we are able to relate prominent geomorphic features in the fault zone to bends in the fault.
Using elastic dislocation models consistent with this mapped fault trace, we then test a range of fault geometries and slip conditions to find the combination which best explains the InSAR displacements. Our favoured model contains a reversal in fault dip, approximately halfway along its length, occurring at the location of a restraining bend. Slip on this model fault is heterogeneous, with two areas of peak slip of 7 m or greater, and components of dip-slip displacement which vary significantly along-strike. The success of this model in fitting the data implies that an observed asymmetry in the coseismic interferograms can be explained in terms of the local fault geometry, rather than by using non-linear elastic rheologies as suggested by earlier authors.     

Determination of the fault plane using a single near-field seismic station with a finite-dimension source model

We explore the possibility of determining the actual fault plane of an earthquake from the inversion of near-source displacement seismograms of one station when a finite-dimension source is used instead of a point source model and when the complete displacement is taken into account, including near-field waves. Tests on synthetic seismograms and real data recorded at local distances show that this is possible even with a single, three-component station. A single accelerogram available for the Erzincan, Turkey, 1992 March 13, M _s = 6.8 earthquake is inverted and the solution found is compatible with other seismological studies and with the mechanism expected for the North Anatolian Fault.     

Applying the relative hypocentre location approach: where was the 1980 November 23 Irpinia earthquake?

The purpose of this work is to evaluate under what conditions it is feasible and with what accuracy it is possible to locate the nucleation point of a large earthquake, given the availability of aftershocks located with high precision by the deployment of a local network. We experiment with several approaches and apply them lo the location of the epicentre of the 1980 November 23 Irpinia earthquake ( M ^w= 6.9).
First we use local P ^g phases selected to optimize the azimuthal coverage, obtaining a well-constrained location with a small statistical error, which typically underestimates the true hypocentre uncertainty.
We then exploit the relative location technique, obtaining stable, almost coincident solutions under three conditions: (1) using multiple independent master events to derive an average epicentre; (2) fitting simultaneously the larger data set for all available master events, using a forward approach; (3) conducting an a priori evaluation of the statistics of station and master events to separate model uncertainties and improve the statistical accuracy of the relative locations. Moreover, only by introducing station statistics can we achieve the desired accuracy of ≅ 1 km in constraining the rupture nucleation point of this large earthquake, and we show that the application of the relative location technique to uncleaned, unweighted data for a single master event provides only a crude epicentre with a confidence ellipse deceivingly smaller than the true hypocentre uncertainty.
The revised epicentre for the 1980 November 23 Irpinia earthquake (48.803 °N-15.302°E) validates the class of multidisciplinary reconstructions of the source process such as the model of Valensise et al. (1989), based on the hypocentre of Westaway & Jackson (1987), and is shifted by almost 13 km to the NW of the epicentre recently proposed by Westaway (1992).     

Human Response to Extreme Events: a review of three post-tsunami disaster case studies

Investigating survivors' behaviour prior to, during and after a disaster provides emergency management agencies with greater understanding of the complexities which influence public response. This information can then be used to develop better community-based disaster risk-reduction strategies. In this paper, we review three post-tsunami disaster case studies: the Indian Ocean tsunami (IOT) on 26 December 2004, the Java tsunami on 17 July 2006 and the South Pacific tsunami on 29 September 2009. The 2004 IOT and 2006 Java tsunami surveys involved delayed-response post-disaster research using video interviewing. The 2009 South Pacific tsunami entailed rapid-response post-disaster research using questionnaire interviews. We highlight the major outcomes of each case study and, based on these, make recommendations for improving tsunami education programs in Australia. These include educating the public about tsunami risk, natural warning signs of tsunamis and regionally specific behavioural response. To help facilitate improvements to future post-disaster research, discussion on survey-related issues from each case study is provided.     

基于信息量与逻辑回归模型的次生滑坡灾害敏感性评价——以汶川县北部为例

次生滑坡灾害的影响是震后较长时间里人们持续关注的焦点,对其开展敏感性评价具有重要意义。选取5.12地震的重灾区汶川县北部作为研究区,利用遥感与地理信息技术提取地震滑坡信息,在全面分析滑坡与高程、坡度、坡向、岩性、断裂带、地震烈度以及水系等7个影响因子相关特性的基础上,采用信息量法与逻辑回归模型进行灾害敏感性评价,将研究区划分为极轻度、轻度、中度、高度和极高危险5个级别,并对不同模型的适用性开展分析和对比。结果表明,逻辑回归模型在描述区域滑坡灾害危险度总体特征方面稍具优势。