Geologic hazards associated with seismogenic faulting in southern Siberia and Mongolia: forms and location patterns

The forms and location patterns of geologic hazards induced by earthquakes in southern Siberia, Mongolia, and northern Kazakhstan in1950 through 2008 have been investigated statistically, using a database of coseismic effects created as a GIS MapInfo application, with a handy input box for large data arrays. The database includes 689 cases of macroseismic effects from M_S = 4.1–8.1 events at 398 sites. Statistical analysis of the data has revealed regional relationships between the magnitude of an earthquake and the maximum distance of its environmental effects (soil liquefaction and subsidence, secondary surface rupturing, and slope instability) to the epicenter and to the causative fault. Thus estimated limit distances to the fault for the M_S = 8.1 largest event are 40 km for soil subsidence (sinkholes), 80 km for surface rupture, 100 km for slope instability (landslides etc.), and 130 km for soil liquefaction. These distances are 3.5–5.6 times as short as those to the epicenter, which are 150, 450, 350, and 450 km, respectively. Analysis of geohazard locations relative to nearest faults in southern East Siberia shows the distances to be within 2 km for sinkholes (60% within 1.5 km), 4.5 km for landslides (90% within 1.5 km), 8 km for liquefaction (69% within 1 km), and 35.5 km for surface rupture (86% within 2 km). The frequency of hazardous effects decreases exponentially away from both seismogenic and nearest faults. Cases of soil liquefaction and subsidence are analyzed in more detail in relation to rupture patterns. Equations have been suggested to relate the maximum sizes of secondary structures (sinkholes, dikes, etc.) with the earthquake magnitude and shaking intensity at the site. As a result, a predictive model has been created for locations of geohazard associated with reactivation of seismogenic faults, assuming an arbitrary fault pattern. The obtained results make basis for modeling the distribution of geohazards for the purposes of prediction and estimation of earthquake parameters from secondary deformation.     

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.     

Influence of the 2008 Wenchuan earthquake (Mw 7.9) on the occurrence probability of future earthquakes on neighboring faults

The Longquan–Shan fault and the Huya fault are two major neighboring faults of the Longmen–Shan fault zone where the 12 May 2008 Wenchuan earthquake (M_w 7.9) occurred. To study the influence of the Wenchuan event on these two active faults, we calculate changes of Coulomb stress on the Longquan–Shan fault and the Huya fault caused by the Wenchuan mainshock. Our results indicate that the Coulomb stress in the northern section (Zone A) of the Longquan–Shan fault is increased by 0.07–0.10 bars, that in the middle section (Zone B) by 0.04–0.11 bars, and that in the southern section (Zone C) shows almost no change. For the Huya fault, the Coulomb stress is decreased by 0.01–0.03 bars in the northern section (Zone A), 0.10–0.35 bars in the middle section (Zone B), and nearly 0.5 bars in the southern section (Zone C). The epicenter distribution of small earthquakes (M_L ⩾ 1.5) on the Longquan–Shan fault and the Huya fault after the Wenchuan earthquake is consistent with the distribution of the Coulomb stress change. This implies that the Wenchuan earthquake may have triggered small events on the Longquan–Shan fault, but inhibited those on the Huya fault. We then use the rate/state friction law to calculate the occurrence probability of future earthquakes in the study region for the next decade. They include the distribution of b-values, magnitude of completeness (M_c), the background seismicity rate, a value of Aσ_n and the duration for the transient effect (t_a) in the study region. We also estimate the earthquake occurrence probabilities on the neighboring faults after the Wenchuan earthquake. Our results show that, the occurrence probability of future earthquakes in the Longquan–Shan has a slight increase, being 7% for M ⩾ 5.0 shocks during the next decade, but the earthquake probability in the Huya region is reduced obviously, being 5–20%, 7–26% and 3–9% for M ⩾ 5.0 shocks during the next decade in sections A, B and C of the Huya fault, respectively.     

Low b-value prior to the Indo-Myanmar subduction zone earthquakes and precursory swarm before the May 1995 M 6.3 earthquake

Some 455 events (m_b ⩾ 4.5) in the Indo-Myanmar subduction zone are compiled using the ISC/EHB/NEIC catalogues (1964–2011) for a systematic study of seismic precursors, b-value and swarm activity. Temporal variation of b-value is studied using the maximum likelihood method beside CUSUM algorithm. The b-values vary from 0.95 to 1.4 for the deeper (depth ⩾60 km) earthquakes, and from 0.85 to 1.3 for the shallower (depth <60 km) earthquakes. A sudden drop in the b-value, from 1.4 to 0.9, prior to the occurrence of larger earthquake(s) at the deeper depth is observed. It is also noted that the CUSUM gradient reversed before the occurrence of larger earthquakes. We further examined the seismicity pattern for the period 1988–1995 within a radius of 150 km around the epicentre (latitude: 24.96°N; longitude: 95.30°E) of a deeper event M 6.3 of May 6, 1995 in this subduction zone. A precursory swarm during January 1989 to July 1992 and quiescence during August 1992 to April 1995 are identified before this large earthquake. These observations are encouraging to monitor seismic precursors for the deeper events in this subduction zone.     

First-order rupture features of the 2011 MW 9.0 Tohoku (Japan) earthquake from surface waves

This study analyzed the rupture directivity of the 2011 Tohoku earthquake by using 100-s Rayleigh-wave travel-times, influenced by the finite source, to derive the fault parameters of the earthquake. The results demonstrated that the earthquake exhibited a slow rupture propagation with a rupture velocity of approximately 1.5–2.0 km/s and asymmetric bilateral faulting. The two rupture directions were N60°E and N127°E, with rupture lengths of approximately 276 km and 231 km, respectively. The rupture toward N60°E had a source duration of approximately 183 s, longer than that toward N127°E (approximately 156 s). Overall, the entire source duration of the earthquake faulting lasted approximately 183 s. Regarding historical seismicity in eastern Japan, the 2011 Tohoku earthquake not only ruptured a locked area in which large earthquakes have rarely occurred, but also ruptured the source regions of several historical earthquakes. With the exception of its slow rupture velocity and generation of a tsunami, the rupture features of the 2011 Tohoku earthquake were inconsistent with those of typical tsunami earthquakes.     

Seismotectonics and large earthquake generation in the Himalayan region

Bounded by the western and eastern syntaxes, the Himalayan region has experienced at least five M ~ 8 earthquakes during a seismically very active phase from 1897 through 1952. However, there has been a paucity of M ~ 8 earthquakes since 1952. Examining of various catalogues and seismograms from the Gottingen Observatory, it is established that this quiescence of M ~ 8 earthquakes is real. While it has not been possible to forecast earthquakes, there has been a success in making a medium term forecast of an M 7.3 earthquake in the adjoining Indo-Burmese arc. Similarly we find that in the central Himalayan region, earthquakes of M > 6.5 have been preceded by seismic swarms and quiescences. In the recent past, based on GPS data, estimates have been made of the accumulated strains and it is postulated that a number of M ~ 8 earthquakes are imminent in the Himalayan region. We examine these estimates and find that while earthquakes of M ~ 8 may occur in the region, however, the available GPS data and their interpretation do not necessarily suggest their size and time of occurrence and whether an earthquake in a particular segment will occur sooner in comparison to that in the neighboring segment. We also comment on the inference of occurrence of M ~ 8 earthquakes based on M8 algorithm for the region. We conclude that while an M ~ 8 earthquake could occur any time anywhere in the Himalayan region, there is no indication as of now as to where and when it would occur. We impress on the need for preparedness to mitigate the pending earthquake disaster in the region.     

The large aftershocks triggered by the 2011 Mw 9.0 Tohoku-Oki earthquake,Japan

The Mw 9.0 Tohoku-Oki earthquake that occurred off the Pacific coast of Japan on March 11, 2011, was followed by thousands of aftershocks, both near the plate interface and in the crust of inland eastern Japan. In this paper, we report on two large, shallow crustal earthquakes that occurred near the Ibaraki-Fukushima prefecture border, where the background seismicity was low prior to the 2011 Tohoku-Oki earthquake. Using densely spaced geodetic observations (GPS and InSAR datasets), we found that two large aftershocks in the Iwaki and Kita-Ibarake regions (hereafter referred to as the Iwaki earthquake and the Kita-Ibarake earthquake) produced 2.1 m and 0.44 m of motion in the line-of-sight (LOS), respectively. The azimuth-offset method was used to obtain the preliminary location of the fault traces. The InSAR-based maximum offset and trace of the faults that produced the Iwaki earthquake are consistent with field observations. The fault location and geometry of these two earthquakes are constrained by a rectangular dislocation model in a multilayered elastic half-space, which indicates that the maximum slips for the two earthquakes are 3.28 m and 0.98 m, respectively. The Coulomb stress changes were calculated for the faults following the 2011 Mw 9.0 Tohoku-Oki earthquake based on the modeled slip along the fault planes. The resulting Coulomb stress changes indicate that the stresses on the faults increased by up to 1.1 MPa and 0.7 MPa in the Iwaki and Kita-Ibarake regions, respectively, suggesting that the Tohoku-Oki earthquake triggered the two aftershocks, supporting the results of seismic tomography.     

Seismic behavior in central Taiwan: Response to stress evolution following the 1999 Mw 7.6 Chi-Chi earthquake

Following the 1999 M_w 7.6 Chi-Chi earthquake, a large amount of seismicity occurred in the Nantou region of central Taiwan. Among the seismic activities, eight M_w ⩾ 5.8 earthquakes took place following the Chi-Chi earthquake, whereas only four earthquakes with comparable magnitudes took place from 1900 to 1998. Since the seismicity rate during the Chi-Chi postseismic period has never returned to the background level, such seismicity activation cannot simply be attributed to modified Omori’s Law decay. In this work, we attempted to associate seismic activities with stress evolution. Based on our work, it appears that the spatial distribution of the consequent seismicity can be associated with increasing coseismic stress. On the contrary, the stress changes imparted by the afterslip; lower crust–upper mantle viscoelastic relaxation; and sequent events resulted in a stress drop in most of the study region. Understanding seismogenic mechanisms in terms of stress evolution would be beneficial to seismic hazard mitigation.     

Precursors and fluid flows in the case of the 1996, ML = 5.2 Saint-Paul-de-Fenouillet earthquake (Pyrenees,France): A complete pre-, co- and post-seismic scenario

Earthquake precursors are now regularly described but often detected only after a major or moderate seismic event. Presence and influence of fluids in the seismogenic processes are often observed at the time of earthquake studies. Even today, the understanding of the physical processes involved in the source region is a real challenge for seismic hazard assessment. Here, the aftershock sequence of the M_L = 5.2, 1996 Saint-Paul-de-Fenouillet (Eastern Pyrenees, France) earthquake is first re-examined with P-wave cross-correlations, resulting in extracting three multiplets and in determining new locations. Multiplets and spatio-temporal distribution analysis of the aftershocks allow for quantifying the hydraulic diffusivity D at a maximum value of 5 m²/s and the permeability K at 10^? 15 m² in the upper Pyrenean crust. Second, a model is established in order to explain the hydrogeochemical transient anomalies, which occurred during the 15 day-period preceding the 1996 earthquake. These anomalies consist on a temporal and spatial sequence of gas emissions in the epicentral area and on chloride and lead concentration variations in a bottled mineral water 25 km north to the main shock epicenter. The proposed model processed in a standard elastic half-space, consists of creep on a low-angle crustal normal-fault, generating volumetric strain field changes over a distance of 25 km from the epicentral area. This model is able to constrain not only the mechanisms and the locations of the geochemical anomalies, but also their timing and probable casual links to the triggering of the impending major event. Also, the active extension proposed here is compatible with seismological observations in the Pyrenees. Thus, the possibility of such creep, which can be considered as a slow-slip event, is discussed in the Pyrenean tectonic and geological context. The model is discussed and compared to previous proposed models on precursor processes of earthquakes, especially concerning the preparation zone concept. Finally, a complete seismic scenario over the period beginning 15 days before the quake and ending 5 days after is proposed and discussed.     

Seismicity of Eastern Algeria: a revised and extended earthquake catalogue

A new revised earthquake catalogue for Eastern Algeria [22°N-38°N–4°E-9.5°E], from 419 to 2008, has been compiled as a starting point for a catalogue dealing with the whole country. The re-evaluation is based on a critical method, and formal criteria are adopted in the parameterization of the earthquake catalogue. Two digital versions are produced: (a) a basic version, presented here, including all the compiled data, the authors of each quoted parameter and the quality rank of the macroseismic information; (b) a second version, obtained after removing doubtful events and parameters, is produced in order to carry out statistical analyses. The newly compiled earthquake catalogue provides a more reliable picture of the seismicity in Eastern Algeria than the previous ones as shown by the comparison of the recent catalogues over the common territory and time span. This catalogue can be easily updated and adapted to the needs of the user.     

Seismicity and velocity structures along the south-Alpine thrust front of the Venetian Alps (NE-Italy)

In this paper we show the seismicity and velocity structure of a segment of the Alpine retro-belt front along the continental collision margin of the Venetian Alps (NE Italy). Our goal is to gain insight on the buried structures and deep fault geometry in a “silent” area, i.e., an area with poor instrumental seismicity but high potential for future earthquakes, as indicated by historical earthquakes (1695 Me = 6.7 Asolo and 1936 Ms = 5.8 Bosco del Cansiglio). Local earthquakes recorded by a dense temporary seismic network are used to compute 3-D Vp and Vp/Vs tomographic images, yielding well resolved images of the upper crust underneath the south-Alpine front. We show the presence of two main distinct high Vp S-verging thrust units, the innermost coincides with the piedmont hill and the outermost is buried under a thick pile of sediments in the Po plain.Background seismicity and Vp/Vs anomalies, interpreted as cracked fluid-filled volumes, suggest that the NE portion of the outermost blind thrust and its oblique/lateral ramps may be a zone of high fluid pressure prone to future earthquakes.Three-dimensional focal mechanisms show compressive and transpressive solutions, in agreement with the tectonic setting, stress field maps and geodetic observations. The bulk of the microseismicity is clustered in two different areas, both in correspondence of inherited lateral ramps of the thrust system. Tomographic images highlight the influence of the paleogeographic setting in the tectonic style and seismic activity of the region.     

Seismicity of the Baikal rift system from regional network observations

In the paper we report the state-of-the-art of seismicity study in the Baikal rift system and the general results obtained. At present, the regional earthquake catalog for fifty years of the permanent instrumental observations consists of over 185,000 events. The spatial distribution of the epicenters, which either gather along well-delineated belts or in discrete swarms is considered in detail for different areas of the rift system. At the same time, the hypocenters are poorly constrained making it difficult to identify the fault geometry. Clustered events like aftershock sequences or earthquake swarms are typical patterns in the region; moreover, aftershocks of M ⩾ 4.7 earthquakes make up a quarter of the whole catalog. The maximum magnitude of earthquakes recorded instrumentally is M_LH7.6 for a strike-slip event in the NE part of the Baikal rift system and M_LH6.8 for a normal fault earthquake in the central part of the rift system (Lake Baikal basin). Predominant movement type is normal faulting on NE striking faults with a left lateral strike-slip component on W–E planes. In conclusion, some shortcomings of the seismic network and data processing are pointed out.     

Location of the pilot borehole for investigations of reservoir triggered seismicity at Koyna,India

Artificial water reservoir triggered earthquakes are now known to have occurred at over 120 sites globally. The part played by the reservoirs in triggering is not exactly known due to lack of near field observations of triggered earthquakes. Koyna, located near the west coast of India, where triggered earthquakes have been occurring since 1962 provides an excellent site for near field observations of the target M ≥ 2 earthquakes. A 6 borehole seismic network has been deployed recently in the Koyna region at depths of 981–1522 m to improve the hypocenter locations. During May–December 2015, a total of 1039 earthquakes of M_L ≥ 0.5 were located using the borehole seismic network. The region is also monitored through a dense network of 23 surface broad-band stations. Our analysis indicates a significant improvement in the estimation of absolute locations of earthquakes with errors of the order of ± 300 m, combining both the networks. Based on seismicity, and logistics, a block of 2 × 2 km² area has been chosen for drilling the first pilot borehole of ~ 3 km depth, where M ≥ 2 earthquakes have been occurring frequently since 2005.     

Possible roles of the Zipingpu Reservoir in triggering the 2008 Wenchuan earthquake

To better understand the role the Zipingpu Reservoir may have played in triggering the 2008 M_s8 (M_w7.9) Wenchuan earthquake in China, this study evaluates changes of Coulomb failure stress (ΔCFS) and assesses their role in local seismicity and their potential impact on the Wenchuan earthquake. In addition, key aspects associated with reservoir-triggered earthquake (RTS), including mechanisms of stress triggering and permeability of fault zones, is briefly reviewed. ΔCFS was calculated at the faults involved in the Wenchuan earthquake due to the combined effects of gravitational loading and pore-pressure diffusion from the impoundment history of the reservoir. ΔCFS on the major source fault is larger than 0.1 MPa in the upper 10 km below the reservoir and reached a few tens of kPa at the focal depth. Such levels of ΔCFS are large enough to modulate the secular stress buildup of a few kPa/yr in the Longmen-shan thrust zone. Based on detailed analysis of numerical results and local seismicity, the author suggests that it is not proper to rule out the possibility of the Wenchuan earthquake being a RTS only based on very limited knowledge from a few cases of historical RTS so far.     

Present activity and seismogenic potential of a low-angle normal fault system (Città di Castello,Italy): Constraints from surface geology,seismic reflection data and seismicity

We present new constraints on an active low-angle normal fault system in the Città di Castello–Sansepolcro basin (CSB) of the northern Apennines of Italy. New field data from the geological survey of the Carta Geologica d'Italia (CARG project) define the surface geometry of the normal fault system and lead to an interpretation of the CROP 03 deep-crust seismic reflection profile (Castiglion Fiorentino–Urbania segment), with particular attention paid to the geometry of the Plio-Quaternary extensional structures. Surface and sub-surface geological data are integrated with instrumental and historical seismicity in order to define the seismotectonics of the area.Low-angle east-dipping reflectors are the seismic expression of the well-known Altotiberina Fault (AF), a regional extensional detachment on which both east- and west-dipping high-angle faults, bounding the CSB, sole out. The AF breakaway zone is located ～ 10 km west of the CSB. Within the extensional allochthon, synthetic east-dipping planes prevail. Displacement along the AF is ～ 4.5 km, which agrees with the cumulative offset due to its synthetic splays. The evolution of the CSB has mainly been controlled by the east-dipping fault system, at least since Early Pleistocene time; this system is still active and responsible for the seismicity of the area. A low level of seismic activity was recorded instrumentally within the CSB, but several damaging earthquakes have occurred in historical times. The instrumental seismicity and the intensity data points of the largest historical earthquakes (5 events with maximum MCS intensity of IX to IX–X) allow us to propose two main seismogenic structures: the Monte Santa Maria Tiberina (Mmax = 5.9) and Città di Castello (Mmax up to 6.5) normal faults. Both are synthetic splays of the AF detachment, dipping to the NE at moderate (45–50°) to low (25–30°) angles and cutting the upper crust up to the surface. This study suggests that low-angle normal faults (at least with dips of 25–30°) may be seismogenic.     

The AD 365 Crete earthquake and possible seismic clustering during the fourth to sixth centuries AD in the Eastern Mediterranean: a review of historical and archaeological data

Historical and archaeological data are used to test geological claims that, in the fourth to sixth centuries AD, the Eastern Mediterranean experienced an unusual clustering of destructive earthquakes (the ‘Early Byzantine Tectonic Paroxsym’). A review of historical accounts of a notable earthquake at this time, that of 21 July AD 365, indicates that this event destroyed nearly all the towns in Crete and was followed by a tsunami which devastated the Nile Delta. The AD 365 event was also probably responsible for reported or observed destruction in ancient towns of west Cyprus and Libya. This earthquake is most likely to be identified with a Hellenic Arc subduction-zone event of ‘great’ (M>8) magnitude, as testified by up to 9 m of uplift in western Crete dated by previous geological studies to around this time. Historical and archaeological data also support the hypothesis that the fourth to sixth centuries AD was a period of abnormally high seismicity in the Eastern Mediterranean. The high seismicity rates of this period may reflect a reactivation of all plate boundaries in the region (Dead Sea Transform, East Anatolian Fault, North Anatolian Fault, Hellenic Arc, Cyprus Arc Fault).     

Source characteristics of the Yutian earthquake in 2008 from inversion of the co-seismic deformation field mapped by InSAR

On 21 March 2008, an Ms7.3 earthquake occurred at Yutian County, Xinjiang Uygur Autonomous Region, which is in the same year as 2008 Mw 7.9 Wenchuan earthquake. These two earthquakes both took place in the Bayar Har block, while Yutian earthquake is located in the west edge and Wenchuan earthquake is in the east. The research on source characteristics of Yutian earthquake can serve to better understand Wenchuan earthquake mechanism. We attempt to reveal the features of the causative fault of Yutian shock and its co-seismic deformation field by a sensitivity-based iterative fitting (SBIF) method. Our work is based on analysis and interpretation to high-resolution satellite (Quickbird) images as well as D-InSAR data from the satellite Envisat ASAR, in conjunction with the analysis of seismicity, focal mechanism solutions and active tectonics in this region. The result shows that the 22 km long, nearly NS trending surface rupture zone by this event lies on a range-front alluvial platform in the Qira County. It is characterized by distinct linear traces and a simple structure with 1–3 m-wide individual seams and maximum 6.5 m width of a collapse fracture. Along the rupture zone are seen many secondary fractures and fault-bounded blocks by collapse, exhibiting remarkable extension. The co-seismic deformation affected a big range 100 km × 40 km. D-InSAR analysis indicates that the interferometric deformation field is dominated by extensional faulting with a small strike-slip component. Along the causative fault, the western wall fell down and the eastern wall, that is the active unit, rose up, both with westerly vergence. The maximum subsidence displacement is ～2.6 m in the LOS, and the maximum uplift is 1.2 m. The maximum relative vertical dislocation reaches 4.1 m, which is 10 km distant from the starting rupture point to south. The 42 km-long seismogenic fault in the subsurface extends in NS direction as an arc, and it dipping angle changes from 70° near the surface to 52° at depth ～10 km. The slip on the fault plane is concentrated in the depth range 0–8 km, forming a belt of length 30 km along strike on the fault plane. There are three areas of concentrating slip, in which the largest slip is 10.5 m located at the area 10 km distant from the initial point of the rupture.     

Detailed fault structure of the Tarutung Pull-Apart Basin in Sumatra,Indonesia, derived from local earthquake data

The Tarutung Basin is located at a right step-over in the northern central segment of the dextral strike-slip Sumatran Fault System (SFS). Details of the fault structure along the Tarutung Basin are derived from the relocations of seismicity as well as from focal mechanism and structural geology. The seismicity distribution derived by a 3D inversion for hypocenter relocation is clustered according to a fault-like seismicity distribution. The seismicity is relocated with a double-difference technique (HYPODD) involving the waveform cross-correlations. We used 46,904 and 3191 arrival differences obtained from catalogue data and cross-correlation analysis, respectively. Focal mechanisms of events were analyzed by applying a grid search method (HASH code). Although there is no significant shift of the hypocenters (10.8 m in average) and centroids (167 m in average), the application of the double difference relocation sharpens the earthquake distribution. The earthquake lineation reflects the fault system, the extensional duplex fault system, and the negative flower structure within the Tarutung Basin. The focal mechanisms of events at the edge of the basin are dominantly of strike-slip type representing the dextral strike-slip Sumatran Fault System. The almost north–south striking normal fault events along extensional zones beneath the basin correlate with the maximum principal stress direction which is the direction of the Indo-Australian plate motion. The extensional zones form an en-echelon pattern indicated by the presence of strike-slip faults striking NE–SW to NW–SE events. The detailed characteristics of the fault system derived from the seismological study are also corroborated by structural geology at the surface.     

Late Holocene relative sea-level changes and the earthquake deformation cycle around upper Cook Inlet,Alaska

Multiple peat-silt couplets preserved in tidal marsh sediment sequences suggest that numerous great plate boundary earthquakes caused the coast around Cook Inlet, Alaska, to subside over the past 3500 years. Field and laboratory analyses of the two youngest couplets record the well-documented earthquake of AD 1964 and the penultimate one, approximately 850 cal yr BP. Diatom assemblages from a range of modern day estuarine environments from tidal flat through salt marsh to acidic bog produce quantitative diatom transfer function models for elevation reconstructions based on fossil samples. Only nine out of 124 fossil assemblages analysed, including previously published data for the AD 1964 earthquake, have a poor modern analogue. Calibration of fossil samples indicate co-seismic subsidence of 1.50±0.32 m for AD 1964, similar to measurements taken after the earthquake, and 1.45±0.34 m for the ∼850 cal yr BP earthquake. Elevation standard errors for individual fossil samples range from ∼0.08 m in peat layers to ∼0.35 m in silt units. Lack of a chronology within fossil silt units prevents identification of changes in the rate of recovery and land uplift between the post-seismic and inter-seismic periods. However, preservation of multiple peat-silt couplets indicates no net emergence over multiple earthquake cycles. Glacio-isostatic movements from Little Ice Age glacier advance and retreat explains a ∼0.15 m relative sea-level oscillation recorded within the peat layer subsequently submerged as a result of the AD 1964 earthquake. Before both this and the ∼850 cal yr BP earthquake, diatom assemblages suggest pre-seismic relative sea-level rise of ∼0.12±0.13 m, representing possible precursors to great earthquakes.