Seismological features of island arc crust as inferred from recent seismic expeditions in Japan

Crustal studies within the Japanese islands have provided important constraints on the physical properties and deformation styles of the island arc crust. The upper crust in the Japanese islands has a significant heterogeneity characterized by large velocity variation (5.5–6.1 km/s) and high seismic attenuation (Qp=100–400 for 5–15 Hz). The lateral velocity change sometimes occurs at major tectonic lines. In many cases of recent refraction/wide-angle reflection profiles, a “middle crust” with a velocity of 6.2–6.5 km/s is found in a depth range of 5–15 km. Most shallow microearthquakes are concentrated in the upper/middle crust. The velocity in the lower crust is estimated to be 6.6–7.0 km/s. The lower crust often involves a highly reflective zone with less seismicity, indicating its ductile rheology. The uppermost mantle is characterized by a low Pn velocity of 7.5–7.9 km/s. Several observations on PmP phase indicate that the Moho is not a sharp boundary with a distinct velocity contrast, but forms a transition zone from the upper mantle to the lower crust. Recent seismic reflection experiments revealed ongoing crustal deformations within the Japanese islands. A clear image of crustal delamination obtained for an arc–arc collision zone in central Hokkaido provides an important key for the evolution process from island arc to more felsic continental crust. In northern Honshu, a major fault system with listric geometry, which was formed by Miocene back arc spreading, was successfully mapped down to 12–15 km.     

Structural characteristics controlling the seismicity crustal structure of southern Japan Trench fore-arc region, revealed by ocean bottom seismographic data

The Japan Trench is a plate convergent zone where the Pacific Plate is subducting below the Japanese islands. Many earthquakes occur associated with plate convergence, and the hypocenter distribution is variable along the Japan Trench. In order to investigate the detailed structure in the southern Japan Trench and to understand the variation of seismicity around the Japan Trench, a wide-angle seismic survey was conducted in the southern Japan Trench fore-arc region in 1998. Ocean bottom seismometers (15) were deployed on two seismic lines: one parallel to the trench axis and one perpendicular. Velocity structures along two seismic lines were determined by velocity modeling of travel time ray-tracing method. Results from the experiment show that the island arc Moho is 18–20 km in depth and consists of four layers: Tertiary and Cretaceous sedimentary rocks, island arc upper and lower crust. The uppermost mantle of the island arc (mantle wedge) extends to 110 km landward of the trench axis. The P-wave velocity of the mantle wedge is laterally heterogeneous: 7.4 km/s at the tip of the mantle wedge and 7.9 km/s below the coastline. An interplate layer is constrained in the subducting oceanic crust. The thickness of the interplate layer is about 1 km for a velocity of 4 km/s. Interplate layer at the plate boundary may cause weak interplate coupling and low seismicity near the trench axis. Low P-wave velocity mantle wedge is also consistent with weak interplate coupling. Thick interplate layer and heterogeneous P-wave velocity of mantle wedge may be associated with the variation of seismic activity.     

Imaging of Vp, Vs, and Poisson’s ratio anomalies beneath Kyushu, southwest Japan: Implications for volcanism and forearc mantle wedge serpentinization

We determine detailed 3-D V_p and V_s structures of the crust and uppermost mantle beneath the Kyushu Island, southwest Japan, using a large number of arrival times from local earthquakes. From the obtained V_p and V_s models, we further calculate Poisson’s ratio images beneath the study area. By using this large data set, we successfully image the 3-D seismic velocity and Poisson’s ratio structures beneath Kyushu down to a depth of 150 km with a more reliable spatial resolution than previous studies. Our results show very clear low V_p and low V_s anomalies in the crust and uppermost mantle beneath the northern volcanoes, such as Abu, Kujyu and Unzen. Low-velocity anomalies are seen in the mantle beneath most other volcanoes. In contrast, there are no significant low-velocity anomalies in the crust or in the upper mantle between Aso and Kirishima. The subducting Philippine Sea slab is imaged generally as a high-velocity anomaly down to a depth of 150 km with some patches of normal to low seismic wave velocities. The Poisson’s ratio is almost normal beneath most volcanoes. The crustal seismicity is distributed in both the high- and low-velocity zones, but most distinctly in the low Poisson’s ratio zone. A high Poisson’s ratio region is found in the forearc crustal wedge above the slab in the junction area with Shikoku and Honshu; this high Poisson’s ratio could be caused by fluid-filled cracks induced by dehydration from the Philippine Sea slab. The Poisson’s ratio is normal to low in the forearc mantle in middle-south Kyushu. This is consistent with the absence of low-frequency tremors, and may indicate that dehydration from the subducting crust is not vigorous in this region.     

Crustal seismicity in central Chile

Both the genesis and rates of activity of shallow intraplate seismic activity in central Chile are poorly understood, mainly because of the lack of association of seismicity with recognizable fault features at the surface and a poor record of seismic activity. The goal of this work is to detail the characteristics of seismicity that takes place in the western flank of the Andes in central Chile. This region, located less than 100 km from Santiago, has been the site of earthquakes with magnitudes up to 6.9, including several 5+ magnitude shocks in recent years. Because most of the events lie outside the Central Chile Seismic Network, at distances up to 60 km to the east, it is essential to have adequate knowledge of the velocity structure in the Andean region to produce the highest possible quality of epicentral locations. For this, a N–S refraction line, using mining blasts of the Disputada de Las Condes open pit mine, has been acquired. These blasts were detected and recorded as far as 180 km south of the mine. Interpretation of the travel times indicates an upper crustal model consisting of three layers: 2.2-, 6.7-, and 6.1-km thick, overlying a half space; their associated P wave velocities are 4.75–5.0 (gradient), 5.8–6.0 (gradient), 6.2, and 6.6 km/s, respectively.Hypocentral relocation of earthquakes in 1986–2001, using the newly developed velocity model, reveals several regions of concentrated seismicity. One clearly delineates the fault zone and extensions of the strike-slip earthquake that took place in September 1987 at the source of the Cachapoal River. Other regions of activity are near the San José volcano, the source of the Maipo River, and two previously recognized lineaments that correspond to the southern extension of the Pocuro fault and Olivares River. A temporary array of seismographs, installed in the high Maipo River (1996) and San José volcano (1997) regions, established the hypocentral location of events with errors of less than 1 km. These events are clustered along no particular lineament approximately 25 km away from the San José and Maipo volcanoes. Recurrence intervals, based on a frequency magnitude relationship for lower magnitude events, indicate that earthquakes with magnitudes of 4.7 and 7 have a repeat time of 1 and 1200 years, respectively. Focal mechanisms of the two largest events indicate horizontal maximum and minimum compressive stresses with σ₁ varying from a NW–SE orientation in the north to E–W at the southern extreme.     

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.     

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     

Seismological structure and implications of collision between the Ontong Java Plateau and Solomon Island Arc from ocean bottom seismometer–airgun data

A seismic refraction–reflection experiment using ocean bottom seismometers and a tuned airgun array was conducted around the Solomon Island Arc to investigate the fate of an oceanic plateau adjacent to a subduction zone. Here, the Ontong Java Plateau is converging from north with the Solomon Island Arc as part of the Pacific Plate. According to our two-dimensional P-wave velocity structure modeling, the thickness of the Ontong Java Plateau is about 33 km including a thick (15 km) high-velocity layer (7.2 km/s). The thick crust of the Ontong Java Plateau still persists below the Malaita Accreted Province. We interpreted that the shallow part of the Ontong Java Plateau is accreted in front of the Solomon Island Arc as the Malaita Accreted Province and the North Solomon Trench are not a subduction zone but a deformation front of accreted materials. The subduction of the India–Australia Plate from the south at the San Cristobal Trench is confirmed to a depth of about 20 km below sea level. Seismicity around our survey area shows shallow (about 50 km) hypocenters from the San Cristobal Trench and deep (about 200 km) hypocenters from the other side of the Solomon Island Arc. No earthquakes occurred around the North Solomon Trench. The deep seismicity and our velocity model suggest that the lower part of the Ontong Java Plateau is subducting. After the oceanic plateau closes in on the arc, the upper part of the oceanic plateau is accreted with the arc and the lower part is subducted below the arc. The estimation of crustal bulk composition from the velocity model indicates that the upper portion and the total of the Solomon Island Arc are SiO₂ 58% and 53%, respectively, which is almost same as that of the Izu–Bonin Arc. This means that the Solomon Island Arc can be a contributor to growing continental crust. The bulk composition of the Ontong Java Plateau is SiO₂ 49–50%, which is meaningfully lower than those of continents. The accreted province in front of the arc is growing with the convergence of the two plates, and this accretion of the upper part of the oceanic plateau may be another process of crustal growth, although the proportion of such contribution is not clear.     

Petrology of some volcanic rock series of the Aeolian arc,Southern Tyrrhenian Sea: Calc-alkaline and shoshonitic associations

In the Aeolian island arc two different magmatological associations, calc-alkaline andesite series and shoshonites, occur in close vicinity. Although both associations erupted simultaneously during the last glaciation, there is a general tendency for the calc-alkaline rocks to be older. Shoshonitic activity is still going on.Calc-alkaline lavas include high-Al basalts, andesites and dacites, with the general characteristics of the island-arc type andesite series. Large cation trace elements (Rb, Ba, Sr) however are distinctly enriched.Shoshonite series include trachybasalts and latites, with which potassium-rich rhyolites can be associated. Leucite tephrites and potassic trachytes form a different evolution trend of the shoshonitic association.Petrology relates both associations of the Aeolian Islands to the island arc dynamics which is presently characterized by deep-focus earthquakes in the depth range of 200–350 km. The present-day gap in seismic activity from 50–200 km coincides with the present-day lack of calc-alkaline volcanic activity and is explained by the model of a detached slab which continues to sink into the mantle.     

Seismotectonics of the Calabrian arc

Relocation of intermediate and deep earthquakes of Tyrrhenian Sea area through joint hypocenter determination for the period 1962–1979, has allowed a more detailed definition of the geometry of this peculiar Benioff zone. Earthquakes dip along a quasi-vertical plane to 250 km depth; there is a 50° dip in the 250–340 km depth range, and a low dip angle to 480 km depth. The structure sketched from the hypocenters is almost continuous, but most energy has been released in the 230–340 km depth interval. An evaluation of fault plane solutions of intermediate earthquakes in this area indicates predominance of down-dip compressions in the central part of the slab. At the border, strike-slip motion occurs independent of depth. Some earthquakes that occurred at intermediate depth (less than 100 km) along the Ionian margin of Calabria show predominance of reverse faulting, with the P-axis oriented SE-NW. However, shallow earthquakes in the Calabria-Sicily region indicate a more complex motion, with predominance of normal faulting. A possible interpretation of these features according to the available geological history, which involves subduction of continental lithosphere, is discussed.     

Triggered seismicity in the Andean arc region via static stress variation by the MW = 8.8, February 27, 2010, Maule earthquake

We localized crustal earthquakes in the Andean arc, between 35°S and 36°S, from December 2009 to May 2010. This research shows a seismicity increase, in a narrow longitudinal area, of more than nine times after the great M_w 8.8 Maule earthquake.The localized seismicity defines an area of ∼80 km long and ∼18 km wide and NNW to NNE trend. The M_d magnitudes varied from 0.7 to 3.1 except for two earthquakes with M_w of 3.9 and 4.5, located in the northern end of the area. The focal mechanisms for these two last events were normal/strike-slip and strike-slip respectively.During 2011, a network of 13 temporary stations was installed in the trasarc region in Malargüe, Argentina. Sixty earthquakes were localized in the study region during an 8 month period.We explored how changes in Coulomb conditions associated with the mega-thrust earthquake triggered subsequent upper-plate events in the arc region. We assumed the major proposed structures as receiver faults and used previously published earthquake source parameters and slip distribution for the Maule quake. The largest contribution to static stress change, up to 5 bars, derives from unclamping resulting consistent with co-seismic dilatational deformation inferred from GPS observations in the region and subsidence in nearby volcanoes caused by magma migration.Three different Quaternary tectonic settings–extensional, strike-slip and compressional-have been proposed for the arc region at these latitudes. We found that the unclamping produced by the Maule quake could temporarily change the local regime to normal/strike-slip, or at least it would favor the activation of Quaternary NNE to N-trending dextral strike-slip faults with dextral transtensional movement.     

3-D seismic structure of the Kachchh,Gujarat, and its implications for the earthquake hazard mitigation

Several pieces of studies on the January 26, 2001, Bhuj earthquake (Mw 7.6) revealed that the mainshock was triggered on the hidden unmapped fault in the western part of Indian stable continental region that caused a huge loss in the entire Kachchh rift basin of Gujarat, India. Occurrences of infrequent earthquakes of Mw 7.6 due to existence of hidden and unmapped faults on the surface have become one of the key issues for geoscientific research, which need to be addressed for evolving plausible earthquake hazard mitigation model. In this study, we have carried out a detailed autopsy of the 2001 Bhuj earthquake source zone by applying three-dimensional (3-D) local earthquake tomography (LET) method to a completely new data set consisting of 576 local earthquakes recorded between November 2006 and April 2009 by a seismic network consisting of 22 numbers of three-component broadband digital seismograph stations. In the present study, a total of 7560 arrival times of P-wave (3820) and S-wave (3740) recorded at least 4 seismograph stations were inverted to assimilate 3-D P-wave velocity (Vp), S-wave velocity (Vs), and Poisson’s ratio (σ) structures beneath the 2001 Bhuj earthquake source zone for reliable interpretation of the imaged anomalies and its bearing on earthquake hazard of the region. The source zone is located near the triple junction formed by juxtapositions of three Indian, Arabian, and Iranian tectonic plates that might have facilitated the process of brittle failure at a depth of 25 km beneath the KRB, Gujarat, which caused a gigantic loss to both property and persons of the region. There may be several hidden seismogenic faults around the epicentral zone of the 2001 Bhuj earthquake in the area, which are detectable using 3-D tomography to minimize earthquake hazard for a region. We infer that the use of detailed 3-D seismic tomography may offer potential information on hidden and unmapped faults beneath the plate interior to unravel the genesis of such big damaging earthquakes. This study may help in evolving a comprehensive earthquake risk mitigation model for regions of analogous geotectonic settings, elsewhere in the world.     

The 2003 earthquake swarm in Anjalankoski, south-eastern Finland

During May 2003 a swarm of 16 earthquakes (M_L = 0.6–2.1) occurred at Anjalankoski, south-eastern Finland. The activity lasted for three weeks, but additional two events were observed at the same location in October 2004. A comparison of the waveforms indicated that the source mechanisms and the hypocentres of the events were nearly identical.A relative earthquake location method was applied to better define the geometry of the cluster and to identify the fault plane associated with the earthquakes. The relocated earthquakes aligned along an ENE–WSW trending zone, with a lateral extent of about 1.0 km by 0.8 km. The relative location and the waveform-modelling of depth sensitive surface wave (Rg) and S-to-P converted body wave (sP) phases indicated that the events were unusually shallow, most likely occurring within the first 2 km of the surface. The revised historical earthquake data confirm that shallow swarm-type seismicity is characteristic to the area.The focal mechanism obtained as a composite solution of the five strongest events corresponds to dip-slip motion along a nearly vertical fault plane (strike = 250°, dip = 80°, rake = 90°). The dip and strike of this nodal plane as well as the relocated hypocentres coincide with an internal intrusion boundary of the Vyborg rapakivi batholith.The events occur under a compressive local stress field, which is explained by large gravitational potential energy differences and ridge-push forces. Pore-pressure changes caused by intrusion of ground water and/or radon gas into the fracture zones are suggested to govern the swarm-type earthquake activity.     

Body-waveforms and source parameters of some moderate-sized earthquakes near North Island, New Zealand

P-wave first motion and synthetic seismogram analysis of P- and SH-waveforms recorded at teleseismic distances on the WWSSN are used to estimate source parameters of seven of the largest earthquakes (6.1 ≤ m_b ≤ 6.3) that occurred in the vicinity of North Island, New Zealand since 1965. The source parameters of three other (m_b ≥ 6.1) events determined outside of this study are included and considered in the final analysis. Four of the earthquakes occurred at shallow depths (< 20 km), of which three were located within and to the north of North Island. Two of the shallow events show strike-slip and normal focal mechanisms with T-axes oriented in a manner consistent with their location in an area of known back-arc extension. One of the shallow events occurred in northern South Island and shows a reverse-type mechanism indicating horizontal contraction of the crust in an easterly azimuth. Six events occurred at intermediate depths (h = 39 to 195 km) of which five exhibit thrust mechanisms with T-axes consistently oriented near vertical. In the light of previously published plate tectonic models, the near vertical orientation of T-axes of the intermediate-depth events may be used to infer that the southern Kermadec plate boundary immediately north of North Island is not strongly coupled, and hence, not likely capable of producing great earthquakes. A similar inference cannot be made for the section of the Hikurangi Margin adjacent to North Island since the intermediate-depth events considered in this study lie to the north of this segment of the plate boundary.     

On the P-wave velocity and plate-tectonics implications for the Tyrrhenian deep-earthquake zone

P-wave velocities in the Tyrrhenian mantle have been determined for the 230–480 km depth range. Analysis of P-wave travel times for a set of Tyrrhenian deep earthquakes gives a velocity-distribution law which shows different behaviours in the 230–300 km and 300–480 km depth intervals. For the first interval the velocity gradient is 0.64 · 10⁻² sec⁻¹ and for the second one it is 0.59 · 10⁻² sec⁻¹. At a depth of 300 km the velocity decreases rapidly from 8.75 to 8.43 km/sec.The results have been analyzed in the framework of a Tyrrhenian structural model characterized by a lithospheric slab dipping 55–60° in the WNW direction.It is also pointed out that the analysis of some geodynamic features of the slabs of Pacific island arcs carried out by Oliver et al. (1973) and Sleep (1973) can be applied to the Tyrrhenian mantle geodynamic features.     

Anatolian surface wave evaluated at GEOFON Station ISP Isparta, Turkey

We have studied seismic surface waves of 255 shallow regional earthquakes recently recorded at GEOFON station ISP (Isparta, Turkey) and have selected these 52 recordings with high signal-to-noise ratio for further analysis. An attempt was made by the simultaneous use of the Rayleigh and Love surface wave data to interpret the planar crust and uppermost mantle velocity structure beneath the Anatolian plate using a differential least-square inversion technique. The shear-wave velocities near the surface show a gradational change from approximately 2.2 to 3.6 km s^− 1 in the depth range 0–10 km. The mid-crustal depth range indicating a weakly developed low velocity zone has shear-wave velocities around 3.55 km s^− 1. The Moho discontinuity characterizing the crust–mantle velocity transition appears somewhat gradual between the depth range  25–45 km. The surface waves approaching from the northern Anatolia are estimated to travel a crustal thickness of  33 km whilst those from the southwestern Anatolia and part of east Mediterranean Sea indicate a thicker crust at  37 km. The eastern Anatolia events traveled even thicker crust at  41 km. A low sub-Moho velocity is estimated at  4.27 km s^− 1, although consistent with other similar studies in the region. The current velocities are considerably slower than indicated by the Preliminary Reference Earth Model (PREM) in almost all depth ranges.     

Source zones of the catastrophic Simushir earthquakes on November 15, 2006 (M w = 8.3) and January 13, 2007 (M w = 8.1) and the deep crust structure beneath the Middle Kuril segment

The data on catastrophic earthquakes with magnitudes of 8.3 and 8.1 that occurred in the Simushir Island area on November 15, 2006, and January 13, 2007, respectively, were compared with the results of land-sea deep seismic studies by different methods (deep seismic sounding, the correlation method of refracted waves, the earthquake converted-wave method, the common mid-point) in the Central Kuril segment. The structure of the Earth’s crust and the hypocentral zones of these earthquakes were analyzed. It was established that the hypocenter of the main shock of the first earthquake was located at the bend of the seismofocal zone under the island slope of the trench on the outer side of the subsiding lithospheric plate in the rapidly rising granulite-basite (ìbasalticî) crustal layer, which, at depths of 7–15 km, replaced the granulite-gneiss layer. This was accompanied by an increase of the seismic wave velocity from 6.4 to 7.1 km/s. The focus of the second earthquake was located beneath the axis of the deep-sea trench. The aftershocks were concentrated in two bands 60–120 km wide that extend along the trench, as well as in the third zone orthogonal to the island arc. It was shown that the epicenters of the earthquakes are linked with regional faults. The main shock of the first earthquake (November 15, 2006) was interpreted as a thrust fault and the second one (January 13, 2007) was attributed to a normal fault.     

Empirical relationships between aftershock area dimensions and magnitude for earthquakes in the Mediterranean Sea region

Fault dimension estimates derived from the aftershock area extent of 36 shallow depth (≤ 31 km) earthquakes that occurred in the Mediterranean Sea region have been used in order to establish empirical relationships between length, width, area and surface-wave/moment magnitude. This dataset consists of events whose aftershock sequence was recorded by a dense local or regional network and the reported location errors did not exceed on average 3–5 km. Surface-wave magnitudes for these events were obtained from the NEIC database and/or published reports, while moment magnitudes as well as focal mechanisms were available from the Harvard/USGS catalogues. Contrary to the results of some previously published studies we found no evidence in our dataset that faulting type may have an effect on the fault dimension estimates and therefore we derived relationships for the whole of the dataset. Comparisons, by means of statistical F-tests, of our relationships with other previously published regional and global relationships were performed in order to check possible similarities or differences. Most such comparisons showed relatively low significance levels (< 95%), since the differences in source dimension estimates were large mainly for magnitudes lower than 6.5, becoming smaller with increasing magnitude. Some degree of similarity, however, could be observed between our fault length relationship and the one derived from aftershock area lengths of events in Greece, while a difference was found between our regional and global fault length relationships. A calculation of the ratio defined as the fault length, derived from our relationships, to the length estimated from regional empirical relationships involving surface ruptures showed that it can take a maximum value of about 7 for small magnitudes while it approaches unity at Ms 7.2. When calculating the same ratio using instead global empirical relationships we see the maximum value not exceeding 1.8, while unity is reached at Mw 7.8, indicating the existence of a strong regional variation in the fault lengths of earthquakes occurring in the Mediterranean Sea region. Also, a relationship between the logarithms of the rupture area and seismic moment is established and it is inferred that there is some variation of stress drop as a function of seismic moment. In particular, it is observed that for magnitudes lower than 6.6 the stress drop fluctuates around 10 bar, while for larger magnitudes the stress drop reaches a value as high as 60 bar.     

Large earthquakes in the New Guinea-Solomon Islands area, 1873–1972

A catalogue of 1873–1972 earthquakes with M > 6.9 for the New Guinea—Solomon Islands region (130–165° E) is compiled. There are 152 events listed. Duda's (1965) results for 1900–1968 are improved for the Papua New Guinea area (141–156° E) because of the availability of historical data for that area.Although there is evidence of rapid Holocene uplift in the main seismic zones, there is little historical evidence for visible uplift or subsidence resulting directly from modern major earthquakes. Coastal subsidences commonly reported as a result of earthquakes are of smaller extent and appear to be due to settlement. However, the occurrence of tsunamigenic earthquakes does suggest that surface deformations do take place off-shore.Using Davies and Brune's (1971) method, regional fault slip rates over 5° -segments of the shallow seismic zone are determined from the seismicity catalogue. The slip rate for the island of New Guinea (Gutenberg and Richter's Region 16) is found to be at least 4.4 cm/y which is almost double the very anomalously low rate of 2.3 cm/y found by Davies and Brune (1971). If allowance is made for shear movement without seismicity and for the approximately ratio of dip-slip versus strike-slip faulting indicated by fault plane solutions, the agreement with Le Pichon's (1970) approach value of 10.7 cm/y for the Pacific—India (Australia) plates is reasonable. The fault slip rate in the area between east New Britain and Bougainville at the Pacific—Bismarck—Solomon triple junction is extremely high (20.6 cm/y at least). The smallest slip rate (1.5 cm/y) is found for westernmost New Guinea (130–135° E).Temporal cumulative summation of moments curves show a periodicity of approximately 25 years in the seismic activity at the triple junction (150–155° E). Elsewhere the rate of seismic activity is aperiodic.     

Double seismic zone beneath the middle of Hokkaido, Japan, in the southwestern side of the Kurile Arc

The vertical section of microearthquakes, determined accurately by using the Hokkaido University network, shows two dipping zones (the double seismic zone) 25–30 km apart in the depth range of 80–150 km beneath the middle of Hokkaido in the southwestern side of the Kurile arc. Hypocentral distribution of large earthquakes (m_b > 4) based on the ISC (International Seismological Centre) bulletin also shows the double seismic zone beneath the same region. The hypocentral distribution indicates that the frequency of events occurring in the lower zone is four times greater than that in the upper zone. The difference in seismic activity between the two zones beneath Hokkaido is in contrast with the region beneath northeastern Honshu in the northeastern Japan arc.Composite focal mechanisms of microearthquakes and individual mechanisms of large events mainly characterize the down-dip extension for the lower zone as is observed beneath northeastern Honshu. For the upper zone, however, the stress field is rather complex and not necessarily similar to that beneath northeastern Honshu. This may be considered to indicate the influence of slab contortion or transformation in the Hokkaido corner between the Kurile and the northeastern Japan arcs.     

New constraints on the current stress field and seismic velocity structure of the eastern Yilgarn Craton from mechanisms of local earthquakes

The Yilgarn Craton has hosted some of the largest earthquakes within the Australian continent in the last 100 years. Earthquakes have mainly been studied in the western part of the craton, and are thought to result from the reactivation of Precambrian structures in an E–W compressive regional stress field imposed by plate-scale processes. Here we present moment tensor solutions for three recent moderate-sized earthquakes around the town of Kalgoorlie that are inconsistent with E–W compression, but instead suggest E–W extension in the eastern Yilgarn Craton. Waveforms of earthquakes at Boulder (M_W = 4.0, 20 April 2010), Kalgoorlie (M_W = 4.3, 26 February 2014) and Coolgardie (M_W = 3.9, 31 October 2014) were inverted for moment tensors. All three earthquakes were shallow (centroid depth ≤4 km) normal-faulting events that occurred along roughly N–S-striking planes, either with a steep westward or a relatively shallow eastward dip. The robustness of the retrieved mechanisms has been thoroughly tested, employing different earth models, assuming different locations for the earthquakes and using different period bands for the inversion. The fit of synthetic long-period waveforms to the observations was in all cases substantially improved by assuming a two-layered crust with high S wavespeeds (about 3.9–4 km/s) overlying substantially slower material. Since there is independent evidence from active source profiles for a P velocity increase between the upper and lower crust, a large difference in v_p/v_s ratio between upper and lower crust is the only way to explain both lines of evidence. This vertical contrast could represent a dominance of felsic material in the upper crust, and substantially more mafic material in the lower crust. Taken together, our results also appear to imply that the regional stress field is E–W extensive in the Kalgoorlie area, and possibly for the entire Kalgoorlie Terrane. This is contrary to current assumptions from continent-scale stress modelling. That the orientations of rupture planes roughly align with the regional structural grain could indicate that Archean structures are reactivated in response to the current stress field.