Passive-margin prolonged volcanism,East Australian Plate: outbursts,progressions, plate controls and suggested causes

Prolonged intraplate volcanism along the 4000 km-long East Australian margin for ca 100 Ma raises many genetic questions. Studies of the age-progressive pulses embedded in general basaltic activity have spawned a host of models. Zircon U–Pb dating of inland Queensland central volcanoes gives a stronger database to consider the structure and origin of Australian age-progressive volcanic chains. This assists appraisal of this volcanism in relation to plate motion and plate margin tectonic models. Inland Queensland central volcanoes progressed south-southeast from 34 to 31 Ma (～5.4 cm/yr) until a surge in activity led to irregular southerly progression 31 to 28 Ma. A new inland southeastern Queensland central volcano line (25 to 22 Ma), from Bunya Mountains to North Main Range, followed 3 Ma behind the adjacent coastal progression. The Australian and Tasman Sea age-progressive chains are compared against recent plate motion modelling (Indian Ocean hotspots). The chain lines differ from general vector traces owing to west-facing swells and cessations in activity. Tectonic processes on the eastern plate margin may regulate these irregularities. These include subduction, rapid roll-back and progressive detachment of the Loyalty slab (43 to 15 Ma). West-flowing Pacific-type asthenosphere, related to perturbed mantle convection, may explain the west-facing volcanic surges. Such westward Pacific flow for over 28 Ma is known at the Australian–Antarctic Discordance, southeast of the present Australian plume sites under Bass Strait–West Tasman Sea. Most basaltic activity along eastern Australia marks asthenospheric melt injections into Tasman rift zone mantle and not lithospheric plate speed. The young (post-10 Ma) fields (Queensland, Victoria–South Australia) reflect new plate couplings, which altered mantle convection and stress regimes. These areas receive asthenospheric inputs from deep thermal zones off northeast Queensland and under Bass Strait.     

Recent plumbing system of the Krakatau volcano revealed by teleseismic earthquake distribution

Spatial and temporal analysis of global seismological data 1964–2005 reveals a distinct teleseismic earthquake activity producing a columnar-like formation in the continental wedge between the Krakatau volcano at the surface and the subducting slab of the Indo-Australian plate. These earthquakes occur continuously in time, are in the body-wave (m _b) magnitude range 4.5–5.3 and in the depth range 1–100 km. The Krakatau earthquake cluster is vertical and elongated in the azimuth N30°E, suggesting existence of a deep-rooted fault zone cutting the Sunda Strait in the SSW-NNE direction. Possible continuation of the fault zone in the SW direction was activated by an intensive 2002/2003 aftershock sequence, elongated in the azimuth of N55°E. Beneath the Krakatau earthquake cluster, an aseismic gap exists in the Wadati-Benioff zone of the subducting plate at the depths 100–120 km. We interpret this aseismic gap as a consequence of partial melting inhibiting stress concentration necessary to generate stronger earthquakes, whereas the numerous earthquakes observed in the overlying lithospheric wedge beneath the volcano probably reflect magma ascent in the recent plumbing system of the Krakatau volcano. Focal depth of the deepest events (~100 km) of the Krakatau cluster constrains the location of the primary magma generation to greater depths. The ascending magmatic fluids stress fault segments within the Sunda Strait fault zone and change their friction parameters inducing the observed tectonic earthquakes beneath Krakatau.     

Tidal triggering of earthquakes in the subducting Philippine Sea plate beneath the locked zone of the plate interface in the Tokai region, Japan

We found a characteristic space–time pattern of the tidal triggering effect on earthquake occurrence in the subducting Philippine Sea plate beneath the locked zone of the plate interface in the Tokai region, central Japan, where a large interplate earthquake may be impending. We measured the correlation between the Earth tide and earthquake occurrence using microearthquakes that took place in the Philippine Sea plate for about two decades. For each event, we assigned the tidal phase angle at the origin time by theoretically calculating the tidal shear stress on the fault plane. Based on the distribution of the tidal phase angles, we statistically tested whether they concentrate near some particular angle or not by using Schuster's test. In this test, the result is evaluated by p-value, which represents the significance level to reject the null hypothesis that earthquakes occur randomly irrespective of the tidal phase angle. As a result of analysis, no correlation was found for the data set including all the earthquakes. However, we found a systematic pattern in the temporal variation of the tidal effect; the p-value significantly decreased preceding the occurrence of M ≥ 4.5 earthquakes, and it recovered a high level afterwards. We note that those M ≥ 4.5 earthquakes were considerably larger than the normal background seismicity in the study area. The frequency distribution of tidal phase angles in the pre-event period exhibited a peak at the phase angle where the tidal shear stress is at its maximum to accelerate the fault slip. This indicates that the observed small p-value is a physical consequence of the tidal effect. We also found a distinctive feature in the spatial distribution of p-values. The small p-values appeared just beneath the strongly coupled portion of the plate interface, as inferred from the seismicity rate change in the past few years.     

Stress field in the Ryukyu–Kyushu Wadati–Benioff zone by inversion of earthquake focal mechanisms

The study addresses the space distribution of the stress field in the Kyushu–Ryukyu Wadati–Benioff zone based on homogeneous data of earthquake focal mechanisms and the inverse technique by Gephart and Forsyth [J. Geophys. Res. 89 (1984) 9305]. The used data set consists of 148 Harvard CMT solutions and 22 earthquake focal mechanisms listed in previous studies. The stress field parameters are determined for 0–40, 41–100 and h>100 km depth ranges. The top 100-km layer of the Wadati–Benioff zone (WBZ) is characterized by strike normal maximum compression σ₁ and steeper than the slab minimum compression σ₃, the last indicating for unbalanced slab pull force. The Tokara channel ‘divides’ the subduction into two parts of different stress regime at depth greater than 100 km. To the south of the channel the slab is under slab parallel σ₁ and slab normal σ₃ while its northern part, beneath Kyushu, is under slab parallel extension and slab normal compression. The results of recent studies on the regional velocity structure and geochemistry of the volcanic lava indicate that the most plausible reason for the observed stress field difference below 100 km in the northern and rest part of the arc is the presence of hot low viscosity upper mantle west of Kyushu.The results of this study indicate that the forces involved in the contemporary subduction dynamics in the Ryukyu–Kyushu Wadati–Benioff zone are related to the convergence between the Philippine Sea Plate and the Eurasian plate, the trench suction force, slab pull, the slab anchor force and, in the southern-central part of the arc, mantle resistance.     

The role of slab melting in the petrogenesis of high-Mg andesites: evidence from Simbo Volcano,Solomon Islands

The petrogenesis of high-Mg andesites (HMA) in subduction zones involves shallow melting of refractory mantle sources or, alternatively, the interaction of ascending slab-derived melts with mantle peridotite. To unravel the petrogenesis of HMA, we report major, trace element and Sr–Nd–Hf–Pb isotope data for a newly found occurrence of HMA in the New Georgia group, Solomon Islands, SW-Pacific. Volcanism in the Solomon Islands was initiated by subduction of the Pacific plate beneath the Indian–Australian plate until a reversal of subduction polarity occurred ca. 10 Ma ago. Currently, the Indian–Australian plate is subducted northeastwards along the San Cristobál trench, forming the younger and still active southwestern Solomon island arc. However, a fossil slab of Pacific crust is still present beneath the arc. The edifice of the active volcano Simbo is located directly in the San Cristobál trench on top of the subducting Indian–Australian plate. Simbo Island lies on top of a strike-slip fault of the adjacent Woodlark spreading centre that is subducted beneath the Pacific plate. Geochemical and petrological compositions of volcanic rocks from Simbo are in marked contrast to those of volcanic rocks from islands north of the trench (mostly arc basalts). Simbo-type rocks are opx-bearing HMA, displaying 60–62 wt% SiO₂ but rather primitive Mg–Ni–Cr characteristics with 4–6 wt% MgO, up to 65 ppm Ni, up to 264 ppm Cr and Mg# from 67 to 75. The compositions of the Simbo andesites are explained by a binary mixture of silicic and basaltic melts. Relict olivine phenocrysts with Fo_88–90 and reaction-rims of opx also support a mixing model. The basaltic endmember is similar to back-arc basalts from the Woodlark Ridge. A slab melt affinity of the silicic mixing component is indicated by Gd_(N)/Yb_(N) of up to 2.2 that is higher if compared to MORB and other arc basalts from the Solomon Islands. ⁸⁷Sr/⁸⁶Sr, ɛNd and ɛHf values in the analysed rocks range from 0.7035 to 0.7040, +6.4 to +7.9 and +12 to +14.4, respectively. These values reveal the presence of the Indian–Australian mantle domain beneath Simbo (i.e. the Indian–Australian plate) and also beneath all other volcanic islands of the New Georgia group, which are located north of the San Cristobál trench. ²⁰⁶Pb/²⁰⁴Pb, ²⁰⁷Pb/²⁰⁴Pb and ²⁰⁸Pb/²⁰⁴Pb values (18.43–18.52, 15.49–15.55 and 18.13–18.34, respectively) confirm the presence of slab melts from the subducted Pacific plate beneath southern Simbo where the highest Gd_(N)/Yb_(N) ratios are reported. A spatial shift towards an Indian–Australian slab signature is observed when approaching the active San Cristobál trench on northern Simbo, reflecting the decreasing influence of slab melts from the old subducted Pacific plate.     

Collision and stress trajectories in Taiwan: a finite element model

Field analyses of Plio-Quaternary compressional deformations in Taiwan have enabled us to reconstruct the paleostress trajectories resulting from the collision of the Luzon arc (Philippine Sea plate) with the Chinese continental margin (Eurasian plate). The direction of the maximum compressional stress σ₁ shows a fan-shaped pattern that we interpret as resulting from the collision of a rigid body (the Luzon arc) indenting a more deformable material (the thick sediments of the Chinese continental margin). Simple analytical models qualitatively explain the fan-shaped pattern, but the influence of various parameters such as boundary conditions and rheology cannot be quantitatively accounted for by this approach. Consequently, we have used a finite element technique to compute the stresses and strains induced by the push of a rigid body against a two-dimensional, viscous material. The boundary conditions are the velocities based on plate kinematics. A motion in the N300°E direction best explains the stress trajectories observed in central Taiwan. Viscosity contrasts as well as small changes in the shape of the northern edge of the indenter have little influence on the computed stress pattern. The most important parameter is the direction of convergence. Our model quantitatively explains the general pattern of the stress trajectories observed in the collision zone of Taiwan, between the Philippine Sea plate and Eurasia.     

Characteristics of Intraslab Stress Field and Tectonic Implication in the Nankai Trough, Japan

1. IntroductionThe Nankai Trough region (Fig. 1.1) of southwest Japan is one of the most tectonically complex subduction zones in the world. The subduction of the Philippine Sea plate (PH) beneath the Eurasian plate (EU) has caused a series of large and great interplate earthquakes. It is generally accepted that great earthquakes have occurred at intervals of 100-150 years along the Nankai subduction zone since the 684 Hakuho earthquake (Fig. 1.2). However, a large earthquake (M>7.5) has…     

Static and dynamic fault parameters of the Saitama earthquake of July 1, 1968

The source mechanism of the Saitama earthquake (36.07°N,139.40°E, M_s = 5.4) of July 1, 1968, is studied on the basis of P-wave first motion, aftershock, long-period surface-wave data and low-magnification long-period seismograms recorded in the nearfield. A precise location of the aftershocks is made using P and S—P time data obtained by a micro-earthquake observatory network. The synthetic near-field seismograms based on the Haskell model are directly compared with the observed near-field seismograms for wave form and amplitude to determine the dynamic fault parameters. The results obtained are as follows: source geometry, reverse dip slip with considerable right-lateral strike-slip component; dip direction, N6°E; dip angle 30°; fault dimension, 10 × 6 km²; rupture velocity, 3.4 km/sec in the direction S30°E; average dislocation, 92 cm; average dislocation velocity, 92 cm/sec; seismic moment, 1.9 · 10²⁵ dyn-cm; stress drop, 100 bar. The effective stress is about the same as the stress drop. For major earthquakes in the Japanese Islands, the dislocation velocity, .D, is found to be proportional to the stress drop, σ. This relation can be expressed by .D - (β/μ)σ, where β is the shear velocity and μ is the rigidity. This result has an importance in engineering seismology because the stress drop scales the seismic motion in the vicinity of an earthquake fault.     

Three-dimensional P- and S-wave velocity structures beneath the Japan Islands obtained by high-density seismic stations by seismic tomography

We construct fine-scale 3D P- and S-wave velocity structures of the crust and upper mantle beneath the whole Japan Islands with a unified resolution, where the Pacific (PAC) and Philippine Sea (PHS) plates subduct beneath the Eurasian (EUR) plate. We can detect the low-velocity (low-V) oceanic crust of the PAC and PHS plates at their uppermost part beneath almost all the Japan Islands. The depth limit of the imaged oceanic crust varies with the regions. High-V_P/V_S zones are widely distributed in the lower crust especially beneath the volcanic front, and the high strain rate zones are located at the edge of the extremely high-V_P/V_S zone; however, V_P/V_S at the top of the mantle wedge is not so high. Beneath northern Japan, we can image the high-V subducting PAC plate using the tomographic method without any assumption of velocity discontinuities. We also imaged the heterogeneous structure in the PAC plate, such as the low-V zone considered as the old seamount or the highly seismic zone within the double seismic zone where the seismic fault ruptured by the earthquake connects the upper and lower layer of the double seismic zone. Beneath central Japan, thrust-type small repeating earthquakes occur at the boundary between the EUR and PHS plates and are located at the upper part of the low-V layer that is considered to be the oceanic crust of the PHS plate. In addition to the low-V oceanic crust, the subducting high-V PAC plate is clearly imaged to depths of approximately 250 km and the subducting high-V PHS zone to depths of approximately 180 km is considered to be the PHS plate. Beneath southwestern Japan, the iso-depth lines of the Moho discontinuity in the PHS plate derived by the receiver function method divide the upper low-V layer and lower high-V layer of our model at depths of 30–50 km. Beneath Kyushu, the steeply subducting PHS plate is clearly imaged to depths of approximately 250 km with high velocities. The high-V_P/V_S zone is considered as the lower crust of the EUR plate or the oceanic crust of the PHS plate at depths of 25–35 km and the partially serpentinized mantle wedge of the EUR plate at depths of 30–45 km beneath southwestern Japan. The deep low-frequency nonvolcanic tremors occur at all parts of the high-V_P/V_S zone—within the zone, the seaward side, and the landward side where the PHS plate encounters the mantle wedge of the EUR plate. We prove that we can objectively obtain the fine-scale 3D structure with simple constraints such as only 1D initial velocity model with no velocity discontinuity.     

Integrated 3D density modelling and segmentation of the Dead Sea Transform

A 3D interpretation of the newly compiled Bouguer anomaly in the area of the “Dead Sea Rift” is presented. A high-resolution 3D model constrained with the seismic results reveals the crustal thickness and density distribution beneath the Arava/Araba Valley (AV), the region between the Dead Sea and the Gulf of Aqaba/Elat. The Bouguer anomalies along the axial portion of the AV, as deduced from the modelling results, are mainly caused by deep-seated sedimentary basins (D > 10 km). An inferred zone of intrusion coincides with the maximum gravity anomaly on the eastern flank of the AV. The intrusion is displaced at different sectors along the NNW–SSE direction. The zone of maximum crustal thinning (depth 30 km) is attained in the western sector at the Mediterranean. The southeastern plateau, on the other hand, shows by far the largest crustal thickness of the region (38–42 km). Linked to the left lateral movement of approx. 105 km at the boundary between the African and Arabian plate, and constrained with recent seismic data, a small asymmetric topography of the Moho beneath the Dead Sea Transform (DST) was modelled. The thickness and density of the crust suggest that the AV is underlain by continental crust. The deep basins, the relatively large intrusion and the asymmetric topography of the Moho lead to the conclusion that a small-scale asthenospheric upwelling could be responsible for the thinning of the crust and subsequent creation of the Dead Sea basin during the left lateral movement. A clear segmentation along the strike of the DST was obtained by curvature analysis: the northern part in the neighbourhood of the Dead Sea is characterised by high curvature of the residual gravity field. Flexural rigidity calculations result in very low values of effective elastic lithospheric thickness (t _e < 5 km). This points to decoupling of crust in the Dead Sea area. In the central, AV the curvature is less pronounced and t _e increases to approximately 10 km. Curvature is high again in the southernmost part near the Aqaba region. Solutions of Euler deconvolution were visualised together with modelled density bodies and fit very well into the density model structures. An erratum to this article can be found at     

Seismic velocity discontinuities in the crust and uppermost mantle beneath the Kanto district, central Japan, identified from receiver function imaging and repeating earthquake activity

We constructed vertical cross-sections of depth-converted receiver function images to estimate the seismic velocity structure of the crust and uppermost mantle beneath the Kanto district, central Japan. Repeating earthquake data for the plate boundary were also used to estimate geometries of the subducting Philippine Sea plate and the subducting Pacific plate. As a result, we present images of some major seismic discontinuities. The upper boundary of the Pacific plate dips to the northwest in northern Kanto and to the west–southwest in southern Kanto with some undulations. On the other hand, the upper boundary of the Philippine Sea plate as a whole dips to the northwest. However, it is concave to the northeast in the southern Boso peninsula. We suggest that the low-velocity mantle wedge may be indicated on the top of both subducting plates. Plate thickness gradually decreases to the northeast. The northeastern end of the Philippine Sea plate is interpreted to be at depths of 45–90 km. The Moho discontinuity in the overriding plate is deeper than 25 km in the northern Kanto. It contacts the subducting Philippine Sea plate in the southwestern part near 35.8°N.     

Repeating earthquake activities associated with the Philippine Sea plate subduction in the Kanto district, central Japan: A new plate configuration revealed by interplate aseismic slips

We detect repeating earthquakes associated with the Philippine Sea plate subduction to reveal the plate configuration. In the Kanto district, we find 140 repeating earthquake groups with 428 events by waveform similarity analysis. Most repeating earthquakes in the eastern part of the Kanto district occur with a regular time interval. They have thrust-type focal mechanisms and are distributed near the upper surface of the Philippine Sea plate. These observations indicate that the repeating earthquakes there occur as a repetition of ruptures on the isolated patches distributed on the plate boundary owing to the concentration of stress caused by aseismic slips in the surrounding areas. This shows that the distributions of repeating earthquakes suggest the aseismic slips in the surrounding areas of small patches. We determine spatial distributions of repeating earthquakes in the eastern part of the Kanto district and find that they correspond to the upper boundary of the Philippine Sea plate, that is, the upper boundary of the oceanic crust layer of the Philippine Sea plate. The plate geometry around Choshi is newly constrained by repeating earthquake data and a rather flat geometry in the eastern part of the Kanto district is revealed. The obtained geometry suggests uplift of the Philippine Sea plate due to the collision with the Pacific plate beneath Choshi.Repeating earthquakes in the western part of the Kanto district have extremely shorter recurrence times, and their focal mechanisms are not of the thrust types. These repeating earthquakes are classified as “burst type” activity and likely to occur on the preexistent fault planes which are distributed around the “collision zone” between the Philippine Sea plate and the inland plate. The variation among the repeating earthquake activities in the Kanto district indicates that regular repetition of repeating earthquakes is possible only on the plate boundary with a smooth and simple geometry.     

Moho depth variation beneath southwestern Japan revealed from the velocity structure based on receiver function inversion

The Philippine Sea plate is subducting under the Eurasian plate beneath the Chugoku-Shikoku region, southwestern Japan. We have constructed depth contours for the continental and oceanic Mohos derived from the velocity structure based on receiver function inversion. Receiver functions were calculated using teleseismic waveforms recorded by the high-density seismograph network in southwestern Japan. In order to determine crustal velocity structure, we first improved the linearized time-domain receiver function inversion method. The continental Moho is relatively shallow ( 30 km) at the coastline of the Sea of Japan and at the Seto Inland Sea, and becomes deeper–greater than 40 km–around 35°N and 133.8°E. Near the Seto Inland Sea, a low-velocity layer of thickness 10 km lies under the continental Moho. This low-velocity layer corresponds to the subducting oceanic crust of the Philippine Sea plate. The oceanic Moho continues to descend from south to northwest and exhibits complicated ridge and valley features. The oceanic Moho runs around 25 km beneath the Pacific coast and 45 km beneath the Seto Inland Sea, and it extends to at least to 34.5°N. The depth variation of the Moho discontinuities is in good qualitative agreement with the concept of isostasy. From the configurations of both the continental and oceanic Mohos, we demonstrate that the continental lower crust and the subducting oceanic crust overlap beneath the southern and central part of Shikoku and that a mantle wedge may exist beneath the western and eastern part of Shikoku. The southern edge of the overlapping region coincides with the downdip limit of the slip area of a megathrust earthquake.     

Elastic and inelastic triggering of earthquakes in the North Anatolian Fault zone

Deformation models used to explain the triggering mechanism often assume pure elastic behaviour for the crust and upper mantle. In reality however, the mantle and possibly the lower crust behave viscoelastically, particularly over longer time scales. Consequently, the stress field of an earthquake is in general time-dependent. In addition, if the elastic stress increase were enough to trigger a later earthquake, this triggered event should occur instantaneously and not many years after the triggering event. Hence, it is appropriate to include inelastic behaviour when analysing stress transfer and earthquake interaction.In this work, we analyse a sequence of 10 magnitude M_s > 6.5 events along the North Anatolian Fault between 1939 and 1999 to study the evolution of the regional Coulomb stress field. We investigate the triggering of these events by stress transfer, taking viscoelastic relaxation into account. We evaluate the contribution of elastic stress changes, of post-seismic viscoelastic relaxation in the lower crust and mantle, and of steady tectonic loading to the total Coulomb stress field. We analyse the evolution of stress in the region under study, as well as on the rupture surfaces of the considered events and their epicentres. We study the state of the Coulomb stress field before the 1999 İzmit and Düzce earthquakes, as well as in the Marmara Sea region.In general, the Coulomb stress failure criterion offers a plausible explanation for the location of these events. However, we show that using a purely elastic model disregards an important part of the actual stress increase/decrease. In several cases, post-seismic relaxation effects are important and greater in magnitude than the stress changes due to steady tectonic loading. Consequently, viscoelastic relaxation should be considered in any study dealing with Coulomb stress changes.According to our study, and assuming that an important part of the rupture surface must be stressed for an earthquake to occur, the most likely value for the viscosity of the lower crust or mantle in this region is 5 · 10¹⁷–10¹⁸ Pa · s. Our results cannot rule out the possibility of other time-dependent processes involved in the triggering of the 1999 Düzce event. However, the stress increase due to viscoelastic relaxation brought 22% of the 1999 Düzce rupture area over the threshold value of Δσ_c ≥ 0.01 MPa (0.1 bar), and took the whole surface closer to failure by an average of 0.2 MPa. Finally, we argue that the Marmara Sea region is currently being loaded with positive Coulomb stresses at a much faster rate than would arise exclusively from steady tectonic loading on the North Anatolian Fault.     

The stress field near the sites of the Meckering (1968) and Calingiri (1970) earthquakes, Western Australia

A comparison of three methods of stress determination using overcoring tests, earthquake focal-mechanism solutions, and fault-displacement analysis, was carried out near the epicentres of the 1968 Meckering and the 1970 Calingiri earthquakes, Western Australia. Shallow overcoring measurements (3–10 m) were made at seven sites in competent granite, along a 200-km north—south traverse through the location of the Meckering earthquake. The in situ measurements indicate a high regional compressive stress, acting about N77°E. This compares with the direction of the P-axes from the Meckering and Calingiri earthquakes of N91°E and N102°E, respectively, and east-west orientation of major principal stress deduced from observations of surface faulting. The highest maximum principal stress, 23 MPa, was measured at the site farthest north from the Meckering epicentre (~ 90 km), and the lowest, 4 MPa, was measured near the epicentre. The magnitude of the stress increases with distance away from the epicentral area at about 0.2 kPa/m. The results were corrected for drilling water temperature, thermal rock stresses, and suction pore pressure. The total correction was less than 2 MPa at most sites.Estimates of shear stress release from the fault-scarp displacements for Meckering and Calingiri were of similar magnitude (~10 MPa) to the difference between observed shear stress at sites close to and sites remote from the Meckering earthquake.As a cause of the earthquake it is proposed that the fault plane is progressively weakened by alteration or weathering, and that a small long-term fluctuation of head of ground water was the trigger.The high stress at shallow depths, the growing evidence of east—west major principal stress in most of the Australian continent, and the existence of horizontal stress higher than gravity lithostatic stress may have a common explanation, in terms of the driving forces of plate tectonics. The Australian results do not agree with the north—south orientation suggested by crude plate-drift theories in which the major stress is in the direction of the drift, and do not appear to fit with orientations from other continents, in the Richardson et al. (1976) interacting plate model which assumes a velocity-dependent plate driving force. An analysis of driving forces at diverging plate boundaries was carried out and suggests no dominant force-velocity relationship, thus forces independent of plate velocity may control the stress orientation in the plate.     

The MW = 6.3 2003 Lefkada earthquake (Greece) and induced stress transfer changes

A large earthquake of magnitude M_W = 6.3 occurred on 14 August 2003 NW of the Lefkada Island, which is situated at the Ionian Sea (western Greece). The source parameters of this event are determined using body-wave modeling. The focal depth was found equal to 9 km, the constrained focal mechanism revealed dextral strike–slip motion (φ = 15°, Δ = 80° and λ = 170°), the duration of the source time function was 8 s and the seismic moment 2.9 × 10²⁵ dyn cm. The earthquake occurred close to the northern end of the Kefallinia transform fault, where the 1994 moderate event and its aftershock sequence were also located. The epicentral distribution of the 2003 aftershock sequence revealed the existence of two clusters. The first one is located close to the epicentral area of the mainshock, while the second southern, close to the northwestern coast of the Kefallinia Island. A gap of seismicity is observed between the two clusters. The length of the activated zone is approximately 60 km. The analysis of data revealed that the northern cluster is directly related to the mainshock, while the southern one was triggered by stress transfer caused by the main event.     

Stress fields and fault reactivation angles of the 2000 western Tottori aftershocks and the 2001 northern Hyogo swarm in southwest Japan

We estimated the stress fields of the aftershocks of the 2000 western Tottori earthquake (M_w 6.6) and the northern Hyogo swarm (max M_w 5.2) by a stress tensor inversion of moment tensor solutions reported from the National Research Institute for Earth Science and Disaster Prevention (Japan). The maximum principal stress direction of the western Tottori sequence was estimated as N107°E with a strike–slip regime. In the northern Hyogo swarm, the orientations of the principal stress directions could not be well constrained by the observed data, but after examining the detailed characteristics of the solution, we obtained a most probable solution of N113°E for the σ₁ direction. These solutions are consistent with the maximum horizontal directions roughly estimated from the strike directions of large earthquakes occurring geographically between these two seismic activities. We measured the angle between each fault–slip direction and maximum principal stress direction to investigate the frictional properties of earthquakes. The distribution of the angles was forward modeled to estimate the coefficient of friction and the stress ratio, assuming uniformly distributed fault orientations. For the western Tottori sequence, a homogeneous stress field with a coefficient of friction less than 0.4 was estimated. A high stress level was also suggested because very little change occurred in the stress field during the mainshock. For the northern Hyogo sequence, the coefficient of friction was estimated to be between 0.5 and 1.0.     

Relations between fracture patterns, seismicity and plate motions in the Lesser Antilles

Measurements of orientations of joints that have been made on the Leeward and Windward Islands (Lesser Antilles) from St. Martin to St. Vincent are reported. The preferred orientations of the joints have been determined by a statistical analysis. Assuming the joints to be shearing phenomena, the principal directions of the stress field that caused them, can be calculated. The latter are in very good agreement with the principal stress directions obtained by a statistical analysis of fault plane solutions of earthquakes in the area; the inferred maximum compression direction being NW-SE (ca. N 130°E). This corresponds to a model which assumes that the American plate is subducted in that direction beneath the Caribbean plate.     

Prediction of the present-day stress field in the Australian continental crust using 3D geomechanical–numerical models

The Australian continent has an enigmatic present-day stress pattern with considerable regional variability in maximum horizontal stress (S_Hmax) orientations. Previous attempts to estimate the Australian S_Hmax orientation with geomechanical–numerical models indicate that plate boundary forces provide the major controls on the contemporary stress orientations. However, these models do not satisfactorily predict the observed stress orientation in major basins throughout eastern Australia, where the knowledge of the present-day crustal stresses is of vital importance for development and management of different types of geo-reservoirs. In addition, a new comprehensive stress-data compilation in Australia, which contains 2150 data records and is the key dataset for model calibration, provides motivation to construct a new geomechanical–numerical model for Australia. Herein, we present a 3D geomechanical–numerical model that predicts both the S_Hmax orientation and the relative stress magnitudes throughout the Australian continent. Our best-fit model, with mean absolute deviation of 15°, is in good agreement with observed S_Hmax orientations and the stress regime in most areas, and shows a much better fit in areas where the stress pattern was unable to be predicted by previous published attempts. Interestingly, the best-fit model requires a significant push from the western boundary of Australian continental model, which is possible supporting evidence for the east–west-oriented mantle drag postulated by state-of-the-art global convection models, or may be generated by the excess of gravitational potential energy from Tibetan Plateau, transferred through the Indo-Australian Plate. Hence, our modelling results provide a good first-order prediction of the stress field for areas where no stress information is currently available and can be used to derive initial and boundary conditions for local and reservoir-scale 3D geomechanical models across Australia.     

The earthquake in the Swabian Jura of 16 November 1911 and present concepts of seismotectonics

Beginning with the Swabian Jura earthquake in 1911 the seismic activity in Central Europe is concentrated to this area. A comparison with other events of the same epicentral region shows that the largest earthquake in Germany has the character of a left-lateral horizontal strike slip striking N — NNE. The focal parameters can be assumed within the following intervalls: Seismic moment M_o = 1…8·10¹⁷ Nm; focal area F_o = 18…53 km²; average dislocation d_o = 20…53 cm and stress drop Δp_o = 13…19 bar.