Fault-plane solution and tectonic implications of the Pattan, Pakistan earthquake of December 28, 1974

A large destructive earthquake occurred on December 28, 1974 on the western bank of the Indus River near the village of Pattan. The earthquake reportedly killed 5,300 persons, injured 17,000 and left 60,000 people homeless. A seismicity map of the region is presented for the period January, 1963, to March, 1974 on a Mercator projection. Two main linear trends are recognized on the epicenter map. The northwest trend, beginning at 32.3°N, 76.6°E terminates at the southwest alignment of epicenters beginning at 36.0°N, 73.5°E and ending at 33.0°N, 71.0°E. The Pattan earthquake occurred near the junction of the two linear trends. A fault-plane solution for this earthquake has been determined from an analysis of teleseismic P-wave first-motion and S-wave polarization data. The strike and dip of the two nodal planes are N65°E, 68°SE and N50°E, 23°NW, respectively. The solution is compatible with and indicates underthrusting of the Indian plate in this region in the NNW direction along a thrust zone striking northeast.     

Earthquake mechanisms and tectonics in the Assam-Burma region

Eleven new focal mechanisms from earthquakes in the Assam-Burma region have been determined using P-wave first-motion directions reported in the Bulletins of the International Seismological Centre (Edinburgh). Out of them, eight mechanisms indicate thrust faulting, two normal faultings and one strike-slip faulting. In the thrust type of mechanism solutions, sense of motion on the shallow dipping of the two nodal planes is consistent with underthrusting beneath the arc-like mountain ranges. Seismic slip vectors strike in almost northerly direction along the eastern Himalayas and in almost easterly direction along the Burmese arc. A predominance of thrust faulting is consistent with geological evidences of thrusting and uplift in the Himalayas and the Assam-Burma region.     

Apparent velocity and azimuths of short period P-Waves recorded at indian WWSSN stations network

Over 200 earthquakes in the distance range 30°–90° and azimuthal range 0°–360°, recorded at Indian WWSSN stations, have been used in the present study. We have treated the four WWSSN recording stations i.e. New Delhi, Poona, Shillong and Kodai-Kanal, as Super Large Aperture Seismic Array (SLASA) Network with New Delhi being its cross-over point. Short period P-wave data as obtained from these stations have been analysed using a least square technique. Slowness and azimuthal anomalies have been computed for all these events. Relative time residuals have also been calculated. A velocity model has been derived on the basis of the slowness and travel-time data. The results do not indicate presence of any triplication in the travel-time curve. Variations in the relative residuals refer to the tectonic features beneath the recording stations. The P-wave velocity increases continuously in the lower mantle region and there is no indication for the presence of any appreciable velocity gradient.     

Focal-mechanism solutions of earthquakes and tectonics of the hindukush region

Seventeen focal-mechanism solutions have been obtained for earthquakes occurring in the Hindukush region using P-wave first motion directions observed from short as well as long period records. These solutions have indicated a thrust type of faulting. Some of the solutions show small components of strike slip motions. The trend of nodal planes in these solutions was found to vary between northeast and southeast directions. The dip of the compressional axes rarely exceeds 25°. Orientation of tensional axes was found to be almost vertical in all cases. These findings together with the spatial distribution of earthquakes in the Hindukush region suggested that earthquakes were caused by down-dip extension within a sinking slab.     

On the seismicity and tectonics associated with the Sinkiang—Tibetan region

A plausible seismo-tectonic boundary of the Sinkiang—Tibetan region is defined on the basis of the trend of higher magnitude earthquakes (M7.0) and energy released by them for the period 1905–1965. In order to study the nature of forces at the northwestern and eastern sides of the region focal mechanisms for eleven shocks have been determined using P-wave first-motion directions reported in the Bulletin of the International Seismological Centre (Edinburgh). Of these, seven mechanisms show thrust faulting, three strike-slip and one normal faulting. The sense of motions of underthrusting blocks in thrust-faulting mechanisms for the two sides are directed towards the Sinkiang—Tibetan region. The slip vectors of strike-slip faulting are also in agreement with the direction of movement of thrust faulting. Thus, the seismicity, energy released, slip vectors and the orientation of T-axes reflect that the northwestern and eastern sides of the Sinkiang—Tibetan region are the plausible seismo-tectonic boundary and the major earthquakes and higher crustal thickness are the results of the movements of surrounding plates towards the region.     

Heating during thrust faulting in the rocky mountains: friction or fiction?

Measurement of the levels of organic maturation below and within shear zones of thrust faults in the Rocky Mountains has revealed no general thermal metamorphism that can be attributed to faulting, with the exception of very localized areas. With few exceptions the vitrinite reflectance values obtained are in the range expected if the maximum level of organic maturation was produced as a result of increasing temperature during progressive burial accompanying sedimentation. Only at Marias Pass evidence has been found to suggest additional maturation as a result of post-orogenic burial below the thrust sheets. Anomalously high vitrinite reflectances obtained from the Lewis thrust, McConnell thrust, Coleman thrust and two unnamed thrusts, are restricted to very narrow films immediatly adjacent to, or within the shear zone which, considering any reasonable thermal conductivity, indicates elevated temperatures were very short lived. The anomalously high vitrinite reflectances within these films, when compared to laboratory heated coals, suggest temperatures in the order of 350°-650°C were locally generated during faulting. Such high temperatures are considered to have been generated during stick-slip faulting at macro-asperities or at ramps on the fault plane where local, and possibly transient, high frictional stresses existed. The absence of evidence for extensive frictional heating supports previous arguments that stable-sliding and/or low frictional stress must exist during thrusting.     

Seismotectonics of the northern Aegean area

Data concerning the focal mechanism and the spatial distribution of earthquakes have been used to investigate the active tectonics of the northern Aegean and the surrounding area.A thrust region, which includes the northernmost part of the Aegean and at least part of the Marmara Sea, has been defined. An amphitheatrical Benioff zone dipping towards the thrust region from south, east and probably from west, at a mean angle of about 30°, has been detected.The thrust region is surrounded by a region of normal faulting. An eastward progression of the seismic activity in this normal faulting region between 1954 and 1971 has been observed.A correspondence between the earthquake occurrence in the thrust and normal faulting regions has also been observed. Each large shock produced by tensional mechanism in the region of normal faulting is preceded or followed by one or more shocks of compressional mechanism in the thrust region.The focal mechanism, the distribution of the earthquake foci with intermediate focal depth, as well as some magnetic and gravimetric observations can be interpreted by assuming that dense oceanic crust sinks in the northern part of this area and that the adjacent lithosphere moves by segmentation to fill the void with the consequence of producing tensile stresses associated with normal faulting. Such a mechanism of lithospheric interaction suggests that accretion probably takes place in this area.     

Three-dimensional P-wave velocity image under the Carpathian arc

An inversion of P-wave travel time residuals from selected earthquakes in the distance range 30°–98° to two seismic station networks was used to model P-wave velocity anomalies down to 250 km depth. In the first inversion experiment a region between 43.5°–47.5°N and 21°–29°E was modelled, using 35 seismic stations, while in the second one a region between 44°–47°N and 25°–29°E was modelled, using 19 seismic stations. The 4-layer block model of the first inversion offers 19% reduction in residual variance, while the 5-layer block model of the second one offers 26% reduction, the rest being explained by noise and smaller scale heterogeneities. The obtained velocity anomalies correlate remarkably well with the gravity anomalies and with the tectonic model for the Vrancea region of Fuchs et al. (1979).     

Crustal structure of the Eastern Alps along the TRANSALP profile from wide-angle seismic tomography

The objective of the TRANSALP project is an investigation of the Eastern Alps with regard to their deep structure and dynamic evolution. The core of the project is a 340-km-long seismic profile at 12°E between Munich and Venice. This paper deals with the P-wave velocity distribution as derived from active source travel time tomography. Our database consists of Vibroseis and explosion seismic travel times recorded at up to 100 seismological stations distributed in a 30-km-wide corridor along the profile. In order to derive a velocity and reflector model, we simultaneously inverted refractions and reflections using a derivative of a damped least squares approach for local earthquake tomography. 8000 travel time picks from dense Vibroseis recordings provide the basis for high resolution in the upper crust. Explosion seismic wide-angle reflection travel times constrain both deeper crustal velocities and structure of the crust–mantle boundary with low resolution. In the resulting model, the Adriatic crust shows significantly higher P-wave velocities than the European crust. The European Moho is dipping south at an angle of 7°. The Adriatic Moho dips north with a gentle inclination at shallower depths. This geometry suggests S-directed subduction. Azimuthal variations of the first-break velocities as well as observations of shear wave splitting reveal strong anisotropy in the Tauern Window. We explain this finding by foliations and laminations generated by lateral extrusion. Based on the P-wave model we also localized almost 100 local earthquakes recorded during the 2-month acquisition campaign in 1999. Seismicity patterns in the North seem related to the Inn valley shear zone, and to thrusting of Austroalpine units over European basement. The alignment of deep seismicity in the Trento-Vicenza region with the top of the Adriatic lower crust corroborates the suggestion of a deep thrust fault in the Southern Alps.     

Focal mechanisms of three earthquakes in Finland and their relation to surface faults

Focal mechanisms for three recent earthquakes in Finland are determined using P-wave polarities together with SV/P and SH/P phase amplitude ratios. The events occurred on May 11, 2000 in Toivakka, Central Finland (M_L=2.4), on September 15, 2000 in Kuusamo, northeastern Finland (M_L=3.5), and on May 2, 2001 in Kolari, western Finnish Lapland (M_L=2.9).In order to obtain reliable estimates of the source parameters, one-dimensional crust and upper mantle velocity models are derived for the epicenter areas from deep-seismic sounding results. The starting models are modified by one-dimensional ray tracing using the earthquake observations. The events are relocated by employing P- and S-phase arrival times from the nearest seismic stations and the final velocity models. Synthetic waveforms, calculated with the reflectivity method, are used to further constrain and verify the source and structural parameters.The Toivakka earthquake indicates thrust- or reverse-faulting mechanism at a depth of 5 km. After comparison with aeromagnetic and topographic data we suggest the eastward dipping nodal plane (358°/42°) was the fault plane. The best-fitting fault plane solution of the Kolari earthquake suggests pure thrust-faulting at a depth of 5 km. The nodal plane striking 035°/30° correlates well with surface observations of the postglacial, possibly listric fault systems in the source area. The Kuusamo earthquake (focal depth 14 km) has a normal-faulting mechanism with the nodal planes trending 133°/47° or 284°/47°. Preference is given to the SW-dipping nodal plane, as it seems to coincide with topographic and magnetic lineament directions that have been active after the last ice age.The three earthquakes have occurred in old Precambrian faults and shear zones, which have been reactivated. The reactivated faults are favourably oriented in the local stress field.     

High-angle reverse faulting associated with the 1946 Nankaido earthquake

Analysis of vertical crustal deformation data in the southwestern part of Shikoku, southwest Japan, suggests that the Nankaido earthquake of 1946 (M_w = 8.1), which is a principal interplate thrust earthquake, was accompanied by subsidiary faulting on a splay fault adjacent to the coast of Shikoku. Discarding crustal movement resulting from the main thrusting of the Nankaido earthquake, local leveling data are explained by slip on a simple rectangular thrust fault located just offshore of Shikoku. Although it is difficult to constrain the fault location, a possible result is a high-angle thrust dipping landward at an angle of about 70°, with a dislocation of about 1.5 m, and source dimensions of 30 × 13 km along strike and dip. respectively. This result indicates that the fault may be one of the steeply dipping subsidiary faults branching from the main low-angle thrust, as was the case in the Alaska earthquake of 1964. Although several lines of evidence suggest that this faulting occurred as slow aseismic slip, its discrimination from the main seismic event is extremely difficult. This kind of high-angle thrusting just offshore of the coast would play an important role for the formation of the marine terraces during the late Quaternary period.     

Seismicity of the Bering Glacier region and its relation to tectonic and glacial processes

We relocated over 1000 earthquakes of magnitude > 0.1 occurring between 1973 and 2001 in the Bering Glacier region of southern Alaska. We used first-motion data from these events to determine focal mechanisms and directly invert for stress orientations. Our results indicate that much of the seismicity in the region is occurring within the North American plate in a zone where an inferred structure, which lies beneath Bering Glacier, intersects the Chugach-St. Elias fault system. Stress-field analysis indicates that the events in the Bering Glacier surge reservoir region are likely occurring on northeast-trending thrust faults, consistent with previous modeling that suggested thrust faulting would be enhanced in regions of ice draw down. We also observe a stress field compatible with either high-angle normal or reverse faulting in regions located northwest of the Bering Glacier. This may indicate localized complexities in interactions between the Bering Glacier structure and the Chugach-St. Elias fault system.     

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.     

Thermal maturity of a fold–thrust belt based on vitrinite reflectance analysis in the Western Foothills complex, western Taiwan

Burial depth, cumulative displacement, and peak temperature of frictional heat of a fault system are estimated by thermal analysis in the fold–thrust belt of the Western Foothills complex, western Taiwan based on the vitrinite reflectance technique. The regional thermal structure across the complex reveals that the rocks were exposed to maximum temperatures ranging from 100 °C to 180 °C, which corresponds to a burial depth of 3.7–6.7 km. A large thermal difference of 90 °C were observed at the Shuilikeng fault which make the eastern boundary of the fold–thrust belt where it is in contact with metamorphic rock of Hsuehshan Range. The large thermal difference corresponds to cumulative displacements on the Shuilikeng fault estimated to be in the range of 5.2–6.9 km. However, thermal differences in across the Shuangtung and Chelungpu faults cannot be determined apparently due to small vertical offsets. The large displacement observed across the Shuilikeng fault is absent at the other faults which are interpreted to be younger faults within the piggyback thrust system. Localized high temperatures adjacent to fault zones were observed in core samples penetrating the Chelungpu fault. Three major fracture zones were observed at core lengths of 225 m, 330 m, and 405 m and the two lower zones which comprise dark gray narrow shear zones. A value of vitrinite reflectance of 1.8%, higher than the background value of 0.8%, is limited at a narrow shear zone of 1 cm thickness at the fracture zone at 330 m. The estimated peak temperature in the range of 550–680 °C in the shear zone is far higher than the background temperature of 130 °C, and it is interpreted as due to frictional heating during seismic faulting.     

An analysis of the gravity field over north Singhbhum, Bihar, India

About seven hundred gravity stations were established 2–3 miles apart over the Precambrian terrain of Singhbhum that lies between latitude 22° 15′ to 23°°15′N and longitude 85° to 87°E. Bouguer anomalies ranging from +4 to −62 mGal are found in the area. The observed Bouguer anomaly map was analyzed into regional and residual components. The residual anomaly map shows an excellent correlation with geology. The Singhbhum granite batholith is associated with several gravity lows. The residual anomaly map outlines nine plutonic granitic masses within the Singhbhum batholith. Negative residuals are also observed over some intrusive granites outside the batholith. Residual gravity highs are noted over the Dalma hill as well as over the Dhanjori lava complex on the eastern part of the Singhbhum batholith.Two-dimensional models suggestive of subsurface configuration of several major geologic units in the area are presented. These indicate that some of the plutonic granites within the Singhbhum batholith are of relatively large dimensions. The basin containing the Iron Ore Group of rocks to the west of the batholith, as well as the basin containing Singhbhum Group of rocks outside the Copper Belt thrust, may have sedimentary thicknesses of the order of 6–7 km. The Dalma lavas attain their maximum thickness of about 2.5 km in the form of a syncline, underneath which the Singhbhum Group of rocks is also found to be the thickest. The Copper Belt thrust, a major Precambrian fracture around the Singhbhum batholith, is moderately north-dipping near the surface but possibly attains a steeper slope at depth. The thrust appears to be quite deep seated. A threedimensional computer-based model for the Dhanjori lava—gabbro complex on the eastern part of the Singhbhum batholith has been deduced. Maximum thickness of these basic rocks is found to exist underneath a thin cap of granophyre. The geological implication of these results is discussed.Variation in the regional anomalies seems to be attributable to a mass deficiency under the Singhbhum batholith. The batholith may extend subsurfacially towards the north across the Copper Belt thrust. The northern tip of the batholith probably became dissected along the line of intersection of the two orogenic trends in the area and subsided. Over this subsided part, the Singhbhum Group of rocks was deposited at a later stage. Gravity data suggest a fairly large amount of subsidence in the area.     

The Central Southern Alps (N. Italy) paleoseismic zone: a comparison between field observations and predictions of fault mechanics

The internal sectors of the Orobic Alps (Northern Italy) are characterised by Alpine age regional shortening showing a transition, through time, from plastic to brittle deformation. Thrust faults cut Alpine ductile folds and are marked by cataclasites and, locally, by pseudotachylytes, suggesting that motion was accommodated by seismic frictional slip. In the Eastern Orobic Alps the thrusting initiated at depths deeper than 10 km (the emplacement depth of the Adamello pluton) and possibly continued at shallower depths. This demonstrates that thrust motion occurred between 10 km depth and the brittle-ductile transition, i.e., at mid-crustal depths. The Orobic Alps exhumed paleoseismic zone shows different geometries along strike. In the central sectors of the Orobic Alps, thrust faults, associated with pseudotachylytes, have average dips around 40° and show no pervasive veining. Much steeper thrusts (dips up to about 85°) occur in the eastern Orobic Alps. In this area, faults are not associated with pervasive veining, i.e., fluid circulation was relatively scarce. This suggests that faulting did not occur with supralithostatic fluid pressure conditions. These reverse faults are severely misoriented (far too steep) for fault reactivation in a sublithostatic fluid pressure regime. We suggest that thrust motion likely started when the faults were less steep and that the faults were progressively rotated up to the present day dips. Domino tilting is probably responsible for this subsequent fault steepening, as suggested by a decrease of the steepness of thrust faults from north to south and by systematic rotations of previous structures consistently with tilting of thrust blocks. When the faults became inclined beyond the fault lock-up angle, no further thrusting was accommodated along them. At later stages regional shortening was accommodated by newly formed lower angle shear planes (dipping around 30–40°), consistently with predictions from fault mechanics.     

Tikhonov's regularization for deconvolution in the empirical Green function method and vertical directivity effect

Application of Tikhonov's technique, using input errors for the parameter of regularization estimation, enhances the accuracy and stability of the reconstruction of a source time function (STF) by the empirical Green function (EGF) method that gives us an opportunity to use simultaneously for analysis body and surface waves data, and to estimate the horizontal and vertical directivity effects. Knowledge of the last is particularly useful for the choice of an active nodal plane of earthquakes with the dip slip fault orientation that allows us to classify these earthquakes to the interplate or intraplate types and thereby to reach the better understanding of tectonic processes in the region of interest.By way of illustration, an attempt to estimate average parameters of faulting in a first approximation is made herein for two Russian Far East large events with opposite types of focal mechanism orientation, strike slip and dip slip. The former is not a matter of interest in the context of vertical directivity effect but enables us to test the method.The directivity analysis of pulse durations and inverse amplitudes of the relative source time functions (RSTFs) restored at eight globally distributed stations IRIS indicates that the destruction in the source of the Neftegorsk earthquake (05/27/1995 M_W=7.1) propagated roughly horizontally in the direction 8±11° during 19.2±0.4 s along the rupture extending 35.5±4.9 km. The calculated slip distribution along the rupture coincides within the error with the results of field geological measurements on the causal surface fault that proves that the Neftegorsk earthquake source is well described by the model of the linear unilateral fault and gives a good assessment of the method applied.The average parameters of faulting in the Kamchatka earthquake (03/08/1999 M_W=6.9) have been determined from data of 13 station IRIS. It was shown that the destruction in its source propagated downward at an angle of about 60° with horizon, in the direction about S156° E, during 13.4±0.2 s, along the rupture totaling 25.5±2.3 km in length. Therefore, the nodal plane, steeply dipped to the SE, was active and this event can be regarded as an intraplate type. Two asperities can be selected; the first with the maximum slip 3.3 m located at a distance of about 7 km from the onset of rupture, and the second with the maximum slip about 0.9 m centered at approximately 19 km from that.     

Stress field in the Northern Rhine area, Central Europe, from earthquake fault plane solutions

Fault plane solutions (FPS) from 110 earthquakes in the northern Rhine area with local magnitudes, ranging from 1.0 to 6.1, and occurring between 1976 and 2002 are determined. FPS are retrieved from P-wave first motions using a grid search approach allowing a detailed exploration of the parameter space. The influence of the 1D velocity model on take-off angles and resulting FPS is examined. All events were relocated with a recently developed minimum 1D model of the velocity structure [J. Geophys. Res. (2003)]. Rose diagrams of the orientation of P, T and B axes show a clear preference of trends of P and T axes at N292°E and N27°E, respectively. The majority of B axes trend in northerly directions. Plunges of P and T axes are mostly around 45° while most B axes are subhorizontal. The main direction of the maximum horizontal stress directly inferred from the fault plane solutions is N118°E.To calculate the orientations of the principal stress axes and the shape of the stress tensor, the inversion method of Gephard and Forsyth [J. Geophys. Res. 89 (1984) 9305] was applied to the whole data set and to several subsets of data. The subsets were formed by grouping events from various geological and tectonic areas and by grouping events into different depth ranges. The subset areas include the Lower Rhine Embayment, the Rhenish Massif, the middle Rhine area, the Neuwied Basin and the area known as the Stavelot–Venn Massif. Inversion of the entire data set shows some ambiguity between a strike-slip and extensional stress regime, with a vertical axis for the medium principal stress and a trend of N305°E and N35°E for the σ₁ and σ₃ axis, respectively, as the best fitting tensor. Earthquakes from the Lower Rhine Embayment and, to some degree, from the middle Rhine area indicate an extensional stress regime. In the Lower Rhine Embayment, plunge and trend of the σ₁ axis are 76° and N162°E and for the σ₃ axis 7° and N42°E. The best fitting solution for the area of the Stavelot–Venn Massif is a strike-slip regime with subhorizontal σ₁ and σ₃ axes with a trend of N316°E and N225°E, respectively. Stress orientations found here agree overall with the results from earlier studies based on smaller data sets. The directions of the maximum and minimum horizontal stresses inverted from focal mechanisms agree well with the stress field predicted by the European Stress Map. This confirms earlier interpretations that the stress field of the Rhine Graben system is controlled by plate driving forces acting on the plate boundaries. However, amplitudes of the stresses change on a local scale and with depth. Estimates of the absolute magnitude of principal stresses favor a normal faulting regime in the shallow crust (above 12-km depth) and a strike-slip regime in the lower crust.     

Nature of accretion related to Paleo-Tethys subduction recorded in northern Thailand: Constraints from mélange kinematics and illite crystallinity

We reconstructed the accretion process related to Paleo-Tethys subduction recorded in northern Thailand, based on mélange and thrust structures, and metamorphic temperatures derived from illite crystallinity data. Mélange formation was characterized by hydrofracturing and cataclastic deformation, with mud injection under semi-lithified conditions followed by shear deformation and pressure solution. Illite crystallinity data suggest metamorphic temperatures below 250 °C during mélange formation. The combined structural and metamorphic data indicate that during mélange formation, the accretionary complex related to Paleo-Tethys subduction developed at shallow levels within an accretionary prism. Asymmetric shear fabrics in mélange indicate top-to-south shear. After correction for rotation associated with collision between the Indian and Eurasian continents, the trend of the Paleo-Tethys subduction zone is estimated to have been N80 °E. We conclude that the Paleo-Tethys was subducted northward beneath the Indochina Block from the Permian to Triassic.