Near-surface earthquakes in the Baltic shield, notably Sweden

A recent series of Swedish earthquakes at a focal depth not exceeding 2–3 km, the largest with I_O = V + (MSK scale) andM_L = 3.2 shows that relatively strong seismic activity can occur in the uppermost part of the Baltic Shield. During the last 15 years several near-surface earthquakes have occurred in this region, as indicated by recorded Rg-waves and/or macroseismic data. Many events are located along the coast of central Sweden, suggesting a seismic belt of minor, near-surface activity, which should be considered in connection with the radioactive waste storage in the Swedish bedrock. The appearance of Rg, common in seismic records of explosions and rockbursts, is not a sufficient discriminator between artificial events and earthquakes.     

Recent activity and seismotectonics of the Eastern Pyrenees

Results from a recent earthquake in the Eastern Pyrenees are presented and the seismotectonics of the region is analyzed from the presently available data. On 26 September 1984 an earthquake (M_L = 4.4) took place in the area of the historical destructive earthquake of 1428. Several portable stations installed in the epicentral area to record aftershocks permitted of defining a precise location at 42°19.2′N, 2°10.2′E and 5 km depth. A maximum felt intensity of V (MSK) is obtained from macroseismic data. The epicentral location lies within a block bounded by E-W-trending structures and the focal solution shows right-lateral shearing with a NW-SE pressure axis.The seismicity in the Eastern Pyrenees shows a complex pattern which can be associated with both E-W fractures and NE-SW fault systems. Focal solutions of another two recent earthquakes of M_L ~ 4, with differences in horizontal pressure axis, are also discussed.     

Site response of the Ganges basin inferred from re-evaluated macroseismic observations from the 1897 Shillong, 1905 Kangra,and 1934 Nepal earthquakes

We analyze previously published geodetic data and intensity values for the M _s = 8.1 Shillong (1897), M _s = 7.8 Kangra (1905), and M _s = 8.2 Nepal/Bihar (1934) earthquakes to investigate the rupture zones of these earthquakes as well as the amplification of ground motions throughout the Punjab, Ganges and Brahmaputra valleys. For each earthquake we subtract the observed MSK intensities from a synthetic intensity derived from an inferred planar rupture model of the earthquake, combined with an attenuation function derived from instrumentally recorded earthquakes. The resulting residuals are contoured to identify regions of anomalous intensity caused primarily by local site effects. Observations indicative of liquefaction are treated separately from other indications of shaking severity lest they inflate inferred residual shaking estimates. Despite this precaution we find that intensites are 1–3 units higher near the major rivers, as well as at the edges of the Ganges basin. We find evidence for a post-critical Moho reflection from the 1897 and 1905 earthquakes that raises intensities 1–2 units at distances of the order of 150 km from the rupture zone, and we find that the 1905 earthquake triggered a substantial subsequent earthquake at Dehra Dun, at a distance of approximately 150 km. Four or more M = 8 earthquakes are apparently overdue in the region based on seismic moment summation in the past 500 years. Results from the current study permit anticipated intensities in these future earthquakes to be refined to incorporate site effects derived from dense macroseismic data.     

Implications of north pacific plate tectonics in central Alaska: Focal mechanisms of earthquakes

The focal mechanisms for 86 selected earthquakes (3.0 m_b 5.5) located in central Alaska have been investigated from P-wave first motions; the data were gathered by local seismic networks. The results show a depth-dependent characteristic to the fault-plane solutions. For earthquakes having focal depths shallower than 60–70 km, the focal mechanisms indicate either strike-slip or normal faults, while for earthquakes with foci at intermediate depths the focal mechanisms correspond to thrust faults. The nature of the seismicity indicates the hinge line of the Pacific lithospheric plate under the study area to be striking N17°E from Cook Inlet towards interior Alaska. The comparison of the focal mechanisms with the seismicity shows that the strike-slip and normal faults are the predominant processes of stress release along the shallow section of the plate. The earthquakes with intermediate foci systematically occur along the inclined section of the plate. If the gently dipping nodal planes for these earthquakes are chosen as the fault planes, the focal mechanisms correspond to underthrust motions at the foci. In these, the slip vectors are oriented either to the west or north with the resultant being in the N30°W direction. The tension axes for the underthrust solutions are also found to be parallel to the local dip of the plate, indicating that the subducted plate in interior Alaska is undergoing gravitational sinking.     

The first peak ground motion attenuation relationships for North of Vietnam

The first attenuation relationships of peak ground acceleration (PGA) and peak ground velocity (PGV) for northern Vietnam are obtained in this study. Ground motion data are collected by a portable broadband seismic network in northern Vietnam as a part of cooperation between the Institute of Geophysics, Vietnamese Academy of Science and Technology, Vietnam and Institute of Earth Sciences, Academia Sinica, Taiwan. The database comprises a total of 330 amplitude records by 14 broadband stations from 53 shallow earthquakes, which were occurred in and around northern Vietnam in the period between 01/2006 and 12/2009. These earthquakes are of local magnitudes between 1.6 and 4.6, focal depths less than 30 km, and epicentral distances less than 500 km. The new attenuation relationships for PGA and PGV are:

log₁₀(PGA)=-0.987+0.7521M_L-log₁₀(R)-0.00475R,     

The seismicity of Kohistan, Pakistan: Source studies of the Hamran (1972.9.3), Darel (1981.9.12) and Patan (1974.12.28) earthquakes

Long period body waves are examined to show that the Hamran (1972.9.3), Darel (1981.9.12) and Patan (1974.12.28) earthquakes in Kohistan had focal depths of about 8–10 km. All involved high angle reverse faulting (thrusting) and had seismic moments of about 2.2 to 2.7·10²⁵ dyne cm. These shallow depths contrast with the deeper hypocentres found in the Hindu Kush and northeast Karakoram to the north and in Hazara to the south. The Hamran and Patan shocks were assigned depths of 45 km by the ISC, indicating that even well-recorded events in this region may have focal depths in error by 30 km     

A study of seismic reflection imaging using microearthquake sources

We present a method for processing three-component digital recordings of microearthquakes to obtain near-vertical reflection profiles in regions of shallow seismicity. The processing includes magnitude and focal depth normalization and event stacking, where stacking is by small localized groups, with ray theoretical time and distance corrections applied to compensate for varying focal depths. In areas with high seismicity, this procedure allows earthquakes to be treated as “controlled” sources to probe layered structures of the deep crust and upper mantle.The validity of our approach is examined using aftershocks of the Borah Peak, Idaho earthquake (M_s = 7.3). Several thousand events occurred in a NNW-trending zone about 10 km wide, 75 km long, and 15 km deep. A small (~ 10 ×10 km) array of nine University of Wisconsin three-component triggered short-period digital seismographs was installed in the region of aftershock activity. Over a 10-day period, about 1000 useable events were recorded, of which about 120 have been used for this study. Hypocenters have been computed using both P- and S-wave arrivals, the latter being essential for stable solutions of events outside the network.The Borah Peak data have been processed to obtain shear-wave reflection profiles for the central station (Station 8) of the digital station array. The stacked shear wave (transverse) record sections reveal coherent reflections from horizons at mid-crustal to Moho depths. The most prominent reflections are from crustal discontinuities in the depth range 18–28 km. Coherent reflections can be obtained only through stacking, which is necessary to improve the signal to noise ratio. The major sources of data scatter, as manifested by “smearing” of reflections on the stacked records, are crustal heterogeneity and errors in the determination of focal depth and origin time.     

Rupture length and velocity for earthquakes in the Mid-Atlantic Ridge from directivity effect in body and surface waves

The dimensions and rupture velocities of four earthquakes, two in the Mid-Atlantic Ridge and two in Iceland with strike–slip mechanisms and magnitudes (M_w) between 6.2 and 6.8 were studied using the directivity effects of Rayleigh and body waves. For Rayleigh waves we used the directivity function for different pairs of stations and for body waves the waveforms of P and SH waves corresponding to a simple extended line source. We have found that three have very shallow depths about 3 km and one 8 km, fault lengths between 12 km and 21 km, and a low rupture velocity of about 1.5 km/s to 2.0 km/s which supports the idea of the presence of slow earthquakes in transform faults.     

Aftershock activity of Bhuj earthquake of january 26th, 2001

Following a large-sized Bhuj earthquake (M _s = 7.6) of January 26th, 2001, a small aperture 4-station temporary local network was deployed, in the epicentral area, for a period of about three weeks and resulted in the recording of more than 1800 aftershocks (-0.07 ≤M _L <5.0). Preliminary locations of epicenters of 297 aftershocks (2.0 ≤M _L <5.0) have brought out a dense cluster of aftershock activity, the center of which falls 20 km NW of Bhachau. Epicentral locations of after-shocks encompass a surface area of about 50 × 40 km² that seems to indicate the surface projection of the rupture area associated with the earthquake. The distribution of aftershock activity above magnitude 3, shows that aftershocks are nonuniformly distributed and are aligned in the north, northwest and northeast directions. The epicenter of the mainshock falls on the southern edge of the delineated zone of aftershock activity and the maximum clustering of activity occurs in close proximity of the mainshock. Well-constrained focal depths of 122 aftershocks show that 89% of the aftershocks occurred at depths ranging between 6 and 25 km and only 7% and 4% aftershocks occur at depths less than 5 and more than 25 km respectively. The Gutenberg-Richter (GR) relationship, logN = 4.52 - 0.89M_L, is fitted to the aftershock data (1.0<-M _L<5.0) and theb-value of 0.89 has been estimated for the aftershock activity.     

Active faulting in the western Pyrénées (France): Paleoseismic evidence for late Holocene ruptures

We have identified a 50-km-long active fault scarp, called herewith the Lourdes Fault, between the city of Lourdes and Arette village in the French Pyrénées. This region was affected by large and moderate earthquakes in 1660 (Io = VIII–IX, MSK 64,), in 1750 (Io = VIII, MSK 64) and in 1967 (M_d = 5.3, Io = VIII, MSK 64). Most earthquakes in this area are shallow and the few available focal mechanism solutions do not indicate a consistent pattern of active deformation. Field investigations in active tectonics indicate an East–West trending and up to 50-m-high fault scarp, in average, made of 3 contiguous linear fault sub-segments. To the north, the fault controls Quaternary basins and shows uplifted and tilted alluvial terraces. Deviated and abandoned stream channels of the southern block are likely due to the successive uplift of the northern block of the fault. Paleoseismic investigations coupled with geomorphic studies, georadar prospecting and trenching along the fault scarp illustrate the cumulative fault movements during the late Holocene. Trenches exhibit shear contacts with flexural slip faulting and thrust ruptures showing deformed alluvial units in buried channels. ¹⁴C dating of alluvial and colluvial units indicates a consistent age bracket from two different trenches and shows that the most recent fault movements occurred between 4221 BC and 2918 BC. Fault parameters and paleoseismic results imply that the Lourdes Fault and related sub-segments may produce a M_W 6.5 to 7.1 earthquake. Fault parameters imply that the Lourdes Fault segment corresponds to a major seismic source in the western Pyrénées that may generate earthquakes possibly larger than the 1660 historical event.     

Characteristics of the 1912 co-seismic rupture along the North Anatolian Fault Zone (Turkey): implications for the expected Marmara earthquake

Evidence of right‐lateral offsets associated with the 1912 earthquake (M_w 7.4) along the North Anatolian Fault (Gaziköy–Saros segment) allow us to survey (using DGPS) the co‐seismic and cumulative slip distribution. The damage distribution and surface breaks related with the earthquake show an elongated zone of maximum intensity (X MSK) parallel to the fault rupture on land but this may extend offshore to the north‐east and south‐west. Detailed mapping of the fault using topographic maps and aerial photographs indicates the existence of pull‐apart basins and pressure ridges. At several localities, the average 1912 offset along strike is 3.5–4 m and cumulative slip is 2–6 times that of individual movement. The fault rupture geometry and slip distribution suggest the existence of three subsegments with a combined total length of 110–120 km, a fault length and maximum slip similar to those of the 1999 Izmit earthquake. The amount of slip at the north‐easternmost section and in the coastal region of the Sea of Marmara reaches an average 4 m, thereby implying the offshore extension of the 1912 rupture. The results suggest that the 1912 event generated up to 150 km of surface faulting, which would imply a M_w 7.2–7.4 earthquake and which, added with rupture lengths of the 1999 earthquakes, help to constrain the remaining seismic gap in the Sea of Marmara.     

Tomographic image of crustal structures across the Chelungpu fault: Is the seismogenic layer structure- or depth-dependent?

In 1999, a large earthquake (M_w = 7.6) occurred along the Chelungpu fault in the fold-and-thrust belt of western Taiwan. To shed more light on the subsurface structures and the seismogenic layers, three-dimensional velocity structures were inverted by using the travel times of both P- and S-waves from 2391 aftershocks recorded by the Central Weather Bureau during the 15 months that followed. From tomography, a typical image of the large-scale thrusting structures in the upper crust across the Chelungpu fault was obtained. In general, high velocities beneath the Western Foothills and Central Ranges are separated from low velocities beneath the Coastal Plain by an east-dipping boundary that is roughly consistent with the Chelungpu fault on the surface. The contrast in velocity on either side of the Chelungpu fault is indicative of about a 7- to 9-km vertical offset in the upper crust. The relocated hypocenter for the Chi-Chi earthquake shifts by 2.2 km toward the northwest, and its focal depth decreases by 0.7 km. A plot of focal depths versus rock velocities where the aftershocks occurred shows earthquakes are more inclined to occur in rock with a velocity of around 5.6 km/s. This strongly suggests the seismogenic layer in the fold-and-thrust belt of Taiwan is more structure-dependent than depth-dependent.     

Seismicity and seismotectonic implications in the southern Baltic Sea area

An updated list of earthquakes and earthquake parameters (location, homogenized magnitude, macroseismic data) for the southern Baltic Sea area reveals activity north of the sea, whereas there are very few epicentres in the sea itself and in the region south of it. This is the first study to combine seismological data for the whole region to cover also the sea. Macroseismic data for the 1930 earthquake were reinvestigated leading to an intensity of V–VI (MM or MSK scale), a radius of perceptibility of 135 km and an unusually big focal depth of about 40 km. It is difficult to correlate individual earthquakes with specific faults, but some seismotectonic relations are suggested, e.g. for the Tornquist zone, the predominant structure of the region. Only few reliable focal-mechanism solutions exist. Possible seismogenic processes (ridge push, isostasy, etc) are discussed.     

Andean earthquakes triggered by the 2010 Maule,Chile (Mw 8.8) earthquake: Comparisons of geodetic,seismic and geologic constraints

The Maule, Chile, (M_w 8.8) earthquake on 27 February 2010 triggered deformation events over a broad area, allowing investigation of stress redistribution within the upper crust following a mega-thrust subduction event. We explore the role that the Maule earthquake may have played in triggering shallow earthquakes in northwestern Argentina and Chile. We investigate observed ground deformation associated with the M_w 6.2 (GCMT) Salta (1450 km from the Maule hypocenter, 9 h after the Maule earthquake), M_w 5.8 Catamarca (1400 km; nine days), M_w 5.1 Mendoza (350 km; between one to five days) earthquakes, as well as eight additional earthquakes without an observed geodetic signal. We use seismic and Interferometric Synthetic Aperture Radar (InSAR) observations to characterize earthquake location, magnitude and focal mechanism, and characterize how the non-stationary, spatially correlated noise present in the geodetic imagery affects the accuracy of our parameter estimates. The focal mechanisms for the far-field Salta and Catamarca earthquakes are broadly consistent with regional late Cenozoic fault kinematics. We infer that dynamic stresses due to the passage of seismic waves associated with the Maule earthquake likely brought the Salta and Catamarca regions closer to failure but that the involved faults may have already been at a relatively advanced stage of their seismic cycle. The near-field Mendoza earthquake geometry is consistent with triggering related to positive static Coulomb stress changes due to the Maule earthquake but is also aligned with the South America-Nazca shortening direction. None of the earthquakes considered in this study require that the Maule earthquake reactivated faults in a sense that is inconsistent with their long-term behavior.     

Source parameters of intermediate-depth Vrancea (Romania) earthquakes from empirical Green's functions modeling

Several source parameters (source dimensions, slip, particle velocity, static and dynamic stress drop) are determined for the moderate-size October 27th, 2004 (M_W = 5.8), and the large August 30th, 1986 (M_W = 7.1) and March 4th, 1977 (M_W = 7.4) Vrancea (Romania) intermediate-depth earthquakes. For this purpose, the empirical Green's functions method of Irikura [e.g. Irikura, K. (1983). Semi-Empirical Estimation of Strong Ground Motions during Large Earthquakes. Bull. Dis. Prev. Res. Inst., Kyoto Univ., 33, Part 2, No. 298, 63–104., Irikura, K. (1986). Prediction of strong acceleration motions using empirical Green's function, in Proceedings of the 7th Japan earthquake engineering symposium, 151–156., Irikura, K. (1999). Techniques for the simulation of strong ground motion and deterministic seismic hazard analysis, in Proceedings of the advanced study course seismotectonic and microzonation techniques in earthquake engineering: integrated training in earthquake risk reduction practices, Kefallinia, 453–554.] is used to generate synthetic time series from recordings of smaller events (with 4 ≤ M_W ≤ 5) in order to estimate several parameters characterizing the so-called strong motion generation area, which is defined as an extended area with homogeneous slip and rise time and, for crustal earthquakes, corresponds to an asperity of about 100 bar stress release [Miyake, H., T. Iwata and K. Irikura (2003). Source characterization for broadband ground-motion simulation: Kinematic heterogeneous source model and strong motion generation area. Bull. Seism. Soc. Am., 93, 2531–2545.] The parameters are obtained by acceleration envelope and displacement waveform inversion for the 2004 and 1986 events and MSK intensity pattern inversion for the 1977 event using a genetic algorithm. The strong motion recordings of the analyzed Vrancea earthquakes as well as the MSK intensity pattern of the 1977 earthquake can be well reproduced using relatively small strong motion generation areas, which corresponds to small asperities with high stress drops (300–1200 bar) and high particle velocities (3–5 m/s). These results imply a very efficient high-frequency radiation, which has to be taken into account for strong ground motion prediction, and indicate that the intermediate-depth Vrancea earthquakes are inherently different from crustal events.     

Variations of b‐values at the western edge of the Ryukyu Subduction Zone,north‐east Taiwan

Using earthquakes relocated in north‐east Taiwan, we estimated b‐value distribution along a cross‐section located near the Ryukyu slab edge, and four b‐value anomalous areas are evidenced: (1) a high b‐value body lying on top of a low Vp, low Vs and high Vp/Vs sausage‐like body was considered as a region of enhanced partial melt or water supply above which seismicity occurs; (2) beneath the Ilan Plain, an anomalous area characterized by b‐values slightly higher than 1.1 might give evidence to the magma conduits to the Kueishantao Island; (3) above the Ryukyu Wadati‐Benioff zone, at depths ranging from 90 to 110 km, a high b‐value anomaly might correspond to the depth where dehydration occurs in the subducting oceanic plate; and (4) a low b‐value area located within the Ryukyu slab, at depths ranging from 70 to 90 km, might be linked to the compressive mechanisms shown by focal mechanisms and the bending of the subducting plate.     

Preliminary analysis of the first digital recordings obtained in the Portuguese seismographic network — Attenuation studies

During experiments with digital stations in the period 1985–1987, twenty five earthquakes with magnitudes m _b in the range 2.9 to 4.8 and epicentres located within the area 36°–42.3° (N) and 4.5°–13.6° (W) were recorded at Montemor (MOE) and Montachique (MTH). The three-component recordings were obtained by Geotech S13 instruments with 1 second period. A preliminary analysis of the recordings consisted in the determination of amplitudes and spectral contents of P and S waves, and led to the following observations: (1) The attenuation of waves is expressed by the equation V = exp(C ₂).R ^C ₁. exp(C ₃.M), where V stands for acceleration, velocity or displacement; M-magnitude; R-focal distance; C ₁, C ₂ and C ₃ are constants to be obtained by least square fitting. The application of this equation led to C ₁ of the order 1.7 for displacement, 1.8 for velocity and 2.0 for acceleration, with an average mean square error 0.8. (2) The ratios L/T (longitudinal/transversal amplitudes), for velocity and displacement, showed a tendency to reduce with increasing focal distance, being 2 for short distances (<50 km) and 1 for long distances (400 km). (3) The ratios S/P (S-wave/P-wave amplitudes), although with a large dispersion, showed a slight tendency for increasing with focal distance. (4) The predominant frequencies also showed a slight tendency to decrease with increasing focal distance and with magnitude. (5) The dependence of C ₁ with frequency (3 to 12 Hz) is well behaved from 0.95 to 1.75 (for the velocity trace).     

Source parameters of earthquakes in the reservoir-triggered seismic (RTS) zone of Koyna–Warna, Western India

New empirical relations are derived for source parameters of the Koyna–Warna reservoir-triggered seismic zone in Western India using spectral analysis of 38 local earthquakes in the magnitude range M _L 3.5–5.2. The data come from a seismic network operated by the CSIR-National Geophysical Research Institute, India, during March 2005 to April 2012 in this region. The source parameters viz. seismic moment, source radius, corner frequency and stress drop for the various events lie in the range of 10¹³–10¹⁶ Nm, 0.1–0.4 km, 2.9–9.4 Hz and 3–26 MPa, respectively. Linear relationships are obtained among the seismic moment (M ₀), local magnitude (M _L), moment magnitude (M _w), corner frequency (fc) and stress drop (?σ). The stress drops in the Koyna–Warna region are found to increase with magnitude as well as focal depths of earthquakes. Interestingly, accurate depths derived from moment tensor inversion of earthquake waveforms show a strong correlation with the stress drops, seemingly characteristic of the Koyna–Warna region.     

Tsunami magnitudes in Taiwan, the Philippines and Indonesia

The regional characteristics of tsunami magnitudes in the SE Asia region are discussed in relation to earthquake magnitudes during the period from 1960 to 1994. Tsunami magnitudes on the Imamura-Iida scale are investigated by the author's method (Hatori 1979, 1986) using the data of inundation heights near the source area and tide-gauge records observed in Japan. The magnitude values of the Taiwan tsunamis showed relatively to be small. On the contrary, the magnitudes of tsunamis in the vicinities of the Philippines and Indonesia exceed more than 1–2 grade (tsunami heights: 2–5 times) compared to earthquakes with similar size on the circum-Pacific zone. The relation between tsunami magnitude, m, and earthquake magnitude, M _s, is expressed as m = 2.66 M _s– 17.5 for these regions. For example, the magnitudes for the 1976 Mindanao tsunami (M _s= 7.8, 3702 deaths) and the 1992 Flores tsunami (M _s= 7.5, 1713 deaths) were determined to be m = 3 and m = 2.5, respectively. The focal depth of tsunamigenic earthquakes is shallower thand< 36 km, and the detectively of tsunamis is small for deep earthquakes being d > 40 km. For future tsunamis, it is indispensable to take precautions against shallow earthquakes having the magnitudes M _s> 6.5.