Large earthquake doublets and fault plane heterogeneity in the northern Solomon Islands subduction zone

In the Solomon Islands and New Britain subduction zones, the largest earthquakes commonly occur as pairs with small separation in time, space and magnitude. This doublet behavior has been attributed to a pattern of fault plane heterogeneity consisting of closely spaced asperities such that the failure of one asperity triggers slip in adjacent asperities. We analyzed body waves of the January 31, 1974,M _w =7.3, February 1, 1974,M _w =7.4, July 20, 1975 (1437)M _w =7.6 and July 20, 1975 (1945),M _w =7.3 doublet events using an iterative, multiple station inversion technique to determine the spatio-temporal distribution of seismic moment release associated with these events. Although the 1974 doublet has smaller body wave moments than the 1975 events, their source histories are more complicated, lasting over 40 seconds and consisting of several subevents located near the epicentral regions. The second 1975 event is well modeled by a simple point source initiating at a depth of 15 km and rupturing an approximate 20 km region about the epicenter. The source history of the first 1975 event reveals a westerly propagating rupture, extending about 50 km from its hypocenter at a depth of 25 km. The asperities of the 1975 events are of comparable size and do not overlap one another, consistent with the asperity triggering hypothesis. The relatively large source areas and small seismic moments of the 1974 doublet events indicate failure of weaker portions of the fault plane in their epicentral regions. Variations in the roughness of the bathymetry of the subducting plate, accompanying subduction of the Woodlark Rise, may be responsible for changes in the mechanical properties of the plate interface.To understand how variations in fault plane coupling and strength affect the interplate seismicity pattern, we relocated 85 underthrusting earthquakes in the northern Solomon Islands Are since 1964. Relatively few smaller magnitude underthrusting events overlap the Solomon Islands doublet asperity regions, where fault coupling and strength are inferred to be the greatest. However, these asperity regions have been the sites of several previous earthquakes withM _s 7.0. The source regions of the 1974 doublet events, which we infer to be mechanically weak, contain many smaller magnitude events but have not generated any otherM _s 7.0 earthquakes in the historic past. The central portion of the northern Solomon Islands Arc between the two largest doublet events in 1971 (studied in detail bySchwartz et al., 1989a) and 1975 contains the greatest number of smaller magnitude underthrusting earthquakes. The location of this small region sandwiched between two strongly coupled portions of the plate interface suggest that it may be the site of the next large northern Solomon Islands earthquake. However, this region has experienced no known earthquakes withM _s 7.0 and may represent a relatively aseismic portion of the subduction zone.     

Geochemical Constraints on Possible Subduction Components in Lavas of Mayon and Taal Volcanoes, Southern Luzon, Philippines

Mayon is the most active volcano along the east margin of southernLuzon, Philippines. Petrographic and major element data indicatethat Mayon has produced a basaltic to andesitic lava seriesby fractional crystallization and magma mixing. Trace elementdata indicate that the parental basalts came from a heterogeneousmantle source. The unmodified composition of the mantle wedgeis similar to that beneath the Indian Ocean. To this mantlewas added a subduction component consisting of melt from subductedpelagic sediment and aqueous fluid dehydrated from the subductedbasaltic crust. Lavas from the highly active Taal Volcano onthe west margin of southern Luzon are compositionally more variablethan Mayon lavas. Taal lavas also originated from a mantle wedgemetasomatized by aqueous fluid dehydrated from the subductedbasaltic crust and melt plus fluid derived from the subductedterrigenous sediment. More sediment is involved in the generationof Taal lavas. Lead isotopes argue against crustal contamination.Some heterogeneity of the unmodified mantle wedge and differencesin whether the sediment signature is transferred into the lavasource through an aqueous fluid or melt phase are needed toexplain the regional compositional variation of Philippine arclavas. KEY WORDS: Mayon Volcano; Philippines; sediment melt; subduction component; Taal Volcano     

Estimating seismic moments and Lg Q using Lg spectra

A multi-event and multi-station inverse method is presented in the paper to simultaneously estimate the seismic moments (M ₀) and source corner frequencies (f _c) of several Jiashi (Xinjiang, China) earthquakes, as well as the apparent Lg Q models for the paths from Jiashi to eight seismic stations (WMQ, AAK, TLG, MAKZ, KUR, VOS, ZRN and CHK) in Central Asia. The resultant seismic moments correlate well with the M ₀ values obtained by Harvard University using the centroid moment tensor (CMT) inversion and the surface-wave magnitudes as well. After the correction by a typical value of average radiation coefficient for regional SV waves, the M ₀ values from Lg spectral inversion are still close to the corresponding values obtained from CMT inversion. The obtained apparent Q _0Lg values (Lg Q at 1 Hz) are consistent with the tectonic features of corresponding propagation paths. The Q _0Lg values are 351±87, 349±86 and 300±27 for the paths from Jiashi to AAK, TLG and MAKZ, respectively. They are smaller than Q _0Lg values for the paths to KUR, VOS, ZRN and CHK, which are 553±72, 569±58, 550±57 and 603±65, respectively. These results agree with the condition that the paths to AAK, TLG and MAKZ mainly propagate through the mountainous Tianshan area where relatively strong seismic activities and large variations of topography are exhibited, while the paths to KUR, VOS, ZRN and CHK mainly propagate through the stable area of Kazak platform. The Q _0Lg value for the path to WMQ is 462±56. This is also in agreement with the condition that the path to WMQ is basically along the border area between Tianshan Mountain and Tarim Basin, and along this path the variations of topography and crustal thickness are moderate in comparison with that along the path to MAKZ. Foundation item: Foundation of Verification Researches for Army Control Technology (413290102).     

The earthquake of 6 September 2002 and the seismic history of Palermo (Northern Sicily,Italy): Implications for the seismic hazard assessment of the city

On September 6, 2002, a M_L = 5.6 earthquake, occurring some tens of kilometres offshore from the Northern Sicilian coast (Southern Tyrrhenian Sea), slightly damaged the city of Palermo and surroundings (degree 6 in the European Macroseismic Scale 1998). The macroseismic investigation of the shock and a detailed study of effects of the main earthquakes which affected Palermo in the past have been performed in order to evaluate the seismic response of the city. Moreover, the comparison of the recent event, which is instrumentally constrained, with historical earthquakes allows us to infer new insights on the seismogenic sources of the area, that seem located offshore in the Tyrrhenian sea.In the last 500 years, Palermo has never been completely destroyed but has suffered effects estimated between intensities 6 and 8 EMS-98 many times (1693, 1726, 1751, 1823, 1940, 1968, 2002). The damage scenarios of the analysed events have shown that damage distribution is strongly conditioned by soil response in the different parts of the city and by a high building vulnerability, mainly in the historical centre and in the south-eastern zone of the modern city. As a matter of fact, Palermo has always suffered greater effects than those reported for other nearby localities. The hazard assessment obtained using observed site intensities has shown that the probability of occurrence for intensity 8 (the strongest intensity observed in Palermo) exceeds 99% for 550 years, while the estimated mean return period is 152 ± 40 years. These results, in connection with building vulnerability due to the urban expansion before the introduction of seismic code, suggest that the city is exposed to a relatively high seismic risk.This paper has not been submitted elsewhere in identical or similar form, nor will it be during the first 3 months after its submission to Journal of Seismology.     

Thermal consequences of the formation of a slab window beneath the Mid-Cretaceous southwest Japan arc: A 2-D numerical analysis

Abstract   The development of voluminous granitic magmatism and widespread high-grade metamorphism in Mid-Cretaceous southwest Japan have been explained by the subduction of a spreading ridge (Kula–Pacific or Farallon–Izanagi plate boundaries) beneath the Eurasian continent and the formation of a slab window. In the present study, the thermal consequences of the formation of a slab window beneath a continental margin are evaluated through a 2-D numerical simulation. The model results are evaluated by comparison with the Mid-Cretaceous geology of southwest Japan. Of particular interest are the absence of an amphibolite- to granulite-facies metamorphic belt near the Wadati–Benioff plane, and significant melting of the lower crustal-mafic rocks sufficient to form a large amount of granitic magma. Because none of the model results simultaneously satisfied these two geological interpretations, it is suggested that subduction of plate boundaries in Mid-Cretaceous southwest Japan was not associated with the opening of a slab window. According to previous studies, and the results of the present study, two different tectonic scenarios could reasonably explain the geological interpretations for Mid-Cretaceous southwest Japan: (i) The spreading ridge did not subduct beneath the Eurasian continent, but was located off the continental margin, implying the continuous subduction of very young oceanic lithosphere; (ii) ridge subduction beneath the continental margin occurred after active spreading had ceased. Consequently, in both tectonic scenarios, the subduction of plate boundaries at the Mid-Cretaceous southwest Japan was not associated with a slab window, but very young (hot) oceanic lithosphere.     

Joint bulk-sound and shear tomography for Western Pacific subduction zones

Detailed regional body wave tomographic inversion of the Western Pacific region has been performed using P and S travel times from common sources and receivers, with a joint inversion in terms of bulk-sound and shear wave-speed variations in the mantle. This technique allows the separation of the influence of bulk and shear moduli, and hence a more direct comparison with mineral physics information. The study region is parameterized with cells of side 0.5° to 2° and 19 layers to a depth of 1500 km, while the rest of the mantle was parameterized with 5×5° cells with 16 layers between the surface and the core–mantle boundary. A simultaneous inversion is made for regional and global structures to minimize the influence of surrounding structures on the regional image. A nested iterative inversion scheme is employed with local linearization and three-dimensional ray tracing through the successive model updates. The results of the regional tomographic inversion reveal the penetration of a subducted slab below the 660 km discontinuity at the Kurile–Kamchatka trench, while flattening of slabs above this depth is observed in the Japan and Izu–Bonin subduction zones on both the bulk-sound and shear wave-speed images. The penetration of a subducted slab down to a depth of at least 1200 km is seen below the southern part of the Bonin trench, Mariana, Philippine, and Java subduction zones. Fast shear wave-speed perturbations associated with the subducted slabs, down to the 410 km transition zone, are larger than the comparable bulk-sound perturbations for all these subduction zones except the Philippines. The bulk-sound signature for the subducted slab is more pronounced than for shear in the Philippines, Talaud, New Guinea, Solomon, and Tonga subduction zones, where penetration of the slab into the middle mantle is observed. Variation in the amplitude ratio between bulk-sound and shear wave-speed anomalies correlates well with the subduction parameters of the descending slab. Slabs younger than 90 Ma at the trench show bulk-sound dominance in the upper mantle, while older slabs have a stronger shear wave-speed signature. Spreading of the fast shear wave-speed zone between 800 and 1000 km is observed in the areas of deep subducted slab penetration, but has no comparable expression in the bulk-sound images. This high-velocity feature may reflect physical or chemical disequilibria introduced to the lower mantle by subducted slabs.     

Silicic arc volcanism in Central Luzon, the Philippines: Characterization of its space, time and geochemical relationship

Abstract   The silicic volcanic rocks in Central Luzon show a temporal and spatial relationship with its geochemistry. Volcanic centers dated to approximately 5 Ma are silicic in geochemical composition whereas those between <5–1 Ma expose basaltic to andesitic rocks. Volcanic centers dated <1 Ma are characterized by a wide range of geochemistry encompassing basaltic through andesitic to dacitic signatures. Aside from changes in geochemistry through time, the areas (i.e. fore-arc to back-arc region) where the volcanic centers are formed also vary. The shift in the location of the volcanic centers in Central Luzon is attributed to changes in the dip of subduction of the South China Sea crust along the Manila Trench. Flat subduction resulted from the subduction of the Scarborough Seamount Chain, an oceanic bathymetric high along the Manila Trench west of northern Luzon. However, collision of Luzon with Taiwan in the north and Palawan in the south resulted in steepening of the subduction angle. The silicic volcanic centers in the forearc (Ce/Yb = 20–140) and back-arc (Ce/Yb = 20–60) regions are generally characterized by higher Ce/Yb compared to the basaltic-andesitic volcanic rocks in the main volcanic arc (Ce/Yb = 20) and back-arc (Ce/Yb = 20–30) regions. This across-arc geochemical variation highlights the contributions from the slab, mantle and crust coupled with the effects of geochemical processes that include partial melting, fractionation, magma mixing and mantle–melt interaction.     

Insights into operation of the subduction factory from the oxygen isotopic values of the southern Izu–Bonin–Mariana Arc

Abstract Oxygen is the most abundant element in the earth, and isotopic analysis of this element in island arc lavas potentially provides sensitive constraints on the proportion of oxygen recycled from subducted material, relative to that extracted from the mantle. Here we report on 225 new oxygen isotopic analyses of whole‐rock and glass samples, and clinopyroxene separates, from lavas collected from the southernmost 1500 km of the Izu–Bonin–Mariana (IBM) convergent margin. Whole‐rock samples clustered around a mean of 6.11 ± 0.47‰, whereas Mariana Trough glasses and mafic melts, calculated to be in equilibrium with mafic phenocrysts, clustered narrowly around a mean of 5.7‰. These data demonstrate that unequivocal identification of magmatic oxygen requires analysis of fresh glass or mafic minerals, and that the source of southern IBM Arc melts is entirely, or almost entirely, in equilibrium with normal mantle oxygen. If the elemental enrichments characteristic of the subduction component originate in subducted materials, these oxygen isotopic data are most consistent with the interaction of a small amount of sediment melt (<4%; mostly less than 1%) with mantle peridotite to yield the hybrid mantle that melts to form IBM Arc magmas.     

The non‐stationary cross‐correlation coefficients in the seismic analysis of MDOF structural systems

A method concerning the evaluation, in a very compact form, of the non‐stationary modal cross‐correlation coefficients of MDOF structural systems subjected to seismic excitations is presented. It is available both in the case when the excitation is considered as a white‐noise process and when it is considered as a filtered process. The evaluation of these coefficients is required when a transient seismic analysis is performed by the use of the modal response spectrum approach. This is necessary when the strong‐motion phase of the earthquake is significantly short with respect to the fundamental period of the structure. Copyright © 2003 John Wiley & Sons, Ltd.     

施甸5级震群前中短-短临的地震活动参数动态演化

通过2001年施甸5级震群前地震学各项参数异常的动态追踪，研究了中强地震孕育、发生的中短—短临的演化过程。结果表明：在云南地区出现b值、小震调制比、波速比等地震活动性参数异常以及中强地震的有规律迁移的中期异常背景上，出现了孕震区的频度、b值、应变加速等群体中短异常，以及云南省内的3、4级地震连续增频、3级日频次异常和显的3、4级地震平静的短临异常。特别是3、4级地震异常平静，标志了明显的临震信息，并初步探讨了中强地震连续发生中如何识别单个地震的震兆问题。