Low b-value prior to the Indo-Myanmar subduction zone earthquakes and precursory swarm before the May 1995 M 6.3 earthquake

Some 455 events (m_b ⩾ 4.5) in the Indo-Myanmar subduction zone are compiled using the ISC/EHB/NEIC catalogues (1964–2011) for a systematic study of seismic precursors, b-value and swarm activity. Temporal variation of b-value is studied using the maximum likelihood method beside CUSUM algorithm. The b-values vary from 0.95 to 1.4 for the deeper (depth ⩾60 km) earthquakes, and from 0.85 to 1.3 for the shallower (depth <60 km) earthquakes. A sudden drop in the b-value, from 1.4 to 0.9, prior to the occurrence of larger earthquake(s) at the deeper depth is observed. It is also noted that the CUSUM gradient reversed before the occurrence of larger earthquakes. We further examined the seismicity pattern for the period 1988–1995 within a radius of 150 km around the epicentre (latitude: 24.96°N; longitude: 95.30°E) of a deeper event M 6.3 of May 6, 1995 in this subduction zone. A precursory swarm during January 1989 to July 1992 and quiescence during August 1992 to April 1995 are identified before this large earthquake. These observations are encouraging to monitor seismic precursors for the deeper events in this subduction zone.     

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.     

Evaluation of earthquake potential along the Northern Anatolian Fault Zone in the Marmara Sea using comparisons of GPS strain and seismotectonic parameters

Seismotectonic parameters including the Gutenberg-Richter b-value and multifractal dimensions D₂ and D₁₅ of seismicity patterns (both spatial and temporal) were compared to GPS-derived maximum shear and dilatation strains measured in the Marmara Sea region of western Turkey along the Northern Anatolian Fault Zone (NAFZ). Comparisons of seismotectonic parameters and GPS-derived maximum shear and dilatation strain along the NAFZ in the vicinity of the 1999 M7.4 Izmit earthquake reveal a positive correlation (r = 0.5, p = 0.05) between average dilatation and the Gutenberg-Richter b-value. Significant negative correlation (r = − 0.56, p = 0.03 and r = − 0.56, p = 0.02) was also observed between the spatial fractal dimension D₂ and GPS-derived maximum geodetic and shear strain. This relationship suggests that, as maximum geodetic and shear strains increase, seismicity becomes increasingly clustered.Anomalous interrelationships are observed in the Marmara Sea region prior to the Izmit event along a bend in the NAFZ near the eastern end of the Marmara Sea known as the Northern Boundary Fault (NBF). An asperity is located near the northwest end of the NBF. Along the 50-km length of the NBF, GPS strains become slightly compressive. The correlation between b-value and GPS-derived dilatation suggests that regions in compression have increased probability of larger magnitude rupture. The NBF appears to serve as an impediment to the transfer of strain from east to west along the NAFZ. Recurrence times for large earthquakes along the NBF are larger than in surrounding areas. Temporal clustering of seismicity in the vicinity of the NBF may represent foreshocks of an impending rupture.     

Weak tidal correlation of NW-Bohemia/Vogtland earthquake swarms

We analyze the possible effect of solid Earth tidal stresses upon a vertical strike-slip fault in NW-Bohemia/Vogtland, central Europe, typical by occurrence of swarm earthquakes. The horizontal components of solid Earth tidal stresses were found strongly to prevail and to reach the level of 2 kPa. We examined tidal triggering as influence of tidal stresses to launching the swarm activity in relative absence of other stress disturbances. The onset times of 46 swarms of mostly M_L < 3 earthquakes that occurred in the period 1991–2005 displayed an increased occurrence near the fortnightly maximum of tidal extensive normal stress. The statistical test however did not prove a statistically significant correlation indicating a triggering effect of fault extension due to tidal loading. We also examined tidal effects to the already running seismic activity of the prominent 2000 swarm by comparing the tidal stress distribution in the investigated period with the distribution of tidal stresses in the occurrence times of each earthquake. The results show that these distributions are almost similar, which indicates that individual earthquakes occur independent of tidal stresses. The unclear tidal correlation of the swarm seismicity may be interpreted by small amplitudes and rates of tidal stress changes compared to the amplitudes and rates of coseismic stress perturbations and of pressure bursts of deep generated fluids.     

Seismicity pattern in north Sumatra-Great Nicobar region: In search of precursor for the 26 December 2004 earthquake

We analyse the seismicity pattern including b-value in the north Sumatra-Great Nicobar region from 1976 to 2004. The analysis suggests that there were a number of significant, intermediate and short-term precursors before the magnitude 7.6 earthquake of 2 November 2002. However, they were not found to be so prominent prior to the magnitude 9.0 earthquake of 26 December 2004 though downward migration of activity and a 50-day short-term quiescence was observed before the event. The various precursors identified include post-seismic and intermediate-term quiescence of 13 and 10 years respectively, between the 1976 (magnitude 6.3) and 2002 earthquakes with two years (1990–1991) of increase in background seismicity; renewed seismicity, downward migration of seismic activity and foreshocks in 2002, just before the mainshock. Spatial variation in b-value with time indicates precursory changes in the form of high b-value zone near the epicenter preceding the mainshocks of 2004 and 2002 and temporal rise in b-value in the epicentral area before the 2002 earthquake.     

The Aegean Sea (Greece) earthquake sequence of 25 March 1986: An application of the ν-value method for earthquake prediction

A technique based on the ν-value, which is defined by % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaSGbaeaaca% GGOaGafqiXdqNbaebacaGGPaWaaWbaaSqabeaacaaIYaaaaaGcbaWa% a0aaaeaacqaHepaDdaahaaWcbeqaaiaaikdaaaaaaaaaaaa!3CEB!\[{{(\bar \tau )^2 } \mathord{\left/ {\vphantom {{(\bar \tau )^2 } {\overline {\tau ^2 } }}} \right. \kern-\nulldelimiterspace} {\overline {\tau ^2 } }}\], where τ is the time interval between two adjacent earthquakes and indicates the pattern of time sequences of earthquakes, has been applied to the 25 March 1986 Aegean Sea (Greece) earthquake (M _L = 5.2) in an attempt to discover temporal changes in seismicity. The analysis of several earthquake sequences revealed that low ν-values preceded the occurrence of relatively large earthquakes. The ν-value technique may be used for monitoring the seismicity changes.     

Analysis of Izmit aftershocks 25 days before the November 12th 1999 Düzce earthquake, Turkey

We investigate spatial clustering of 2414 aftershocks along the Izmit M_w = 7.4 August 17, 1999 earthquake rupture zone. 25 days prior to the Düzce earthquake M_w = 7.2 (November 12, 1999), we analyze two spatial clusters, namely Sakarya (SC) and Karadere–Düzce (KDC). We determine the earthquake frequency–magnitude distribution (b-value) for both clusters. We find two high b-value zones in SC and one high b-value zone in KDC which are in agreement with large coseismic surface displacements along the Izmit rupture. The b-values are significantly lower at the eastern end of the Izmit rupture where the Düzce mainshock occurred. These low b-values at depth are correlated with low postseismic slip rate and positive Coloumb stress change along KDC. Since low b-values are hypothesized with high stress levels, we propose that at the depth of the Düzce hypocenter (12.5 km), earthquakes are triggered at higher stresses compared to shallower crustal earthquake. The decrease in b-value from the Karadere segment towards the Düzce Basin supports this low b-value high stress hypothesis at the eastern end of the Izmit rupture. Consequently, we detect three asperity regions which are correlated with high b-value zones along the Izmit rupture. According to aftershock distribution the half of the Düzce fault segment was active before the 12 November 1999 Düzce mainshock. This part is correlated with low b-values which mean high stress concentration in the Düzce Basin. This high density aftershock activity presumably helped to trigger the Düzce event (M_w = 7.2) after the Izmit M_w 7.4 mainshock.     

Verification of our previous definition of preferred earthquake nucleation areas in Kanto-Tokai, Japan

We have proposed that points of future initiation of rupture may be mapped, based on minima in local recurrence times, which are equivalent to local maxima in the probability for main shocks to occur. These minima are often controlled by anomalously low b-values (logN = a − bM). Of the Kanto-Tokai area, approximately 12% showed anomalously short recurrence times and was proposed as asperities, based on seismicity up to 1999. During the period 1999–2003.5, about 75% of the earthquakes with M ≥ 3.5 fell into the asperities, earlier defined (for example 19 out of 23 M ≥ 3.8 events). The probability for this to occur by chance is approximately 2 10^− 14. This supports our idea that the most likely volumes to produce main shocks may be mapped by minima in local recurrence times.     

Seismic potential of Southern Italy

To improve estimates of the long-term average seismic potential of the slowly straining South Central Mediterranean plate boundary zone, we integrate constraints on tectonic style and deformation rates from geodetic and geologic data with the traditional constraints from seismicity catalogs. We express seismic potential (long-term average earthquake recurrence rates as a function of magnitude) in the form of truncated Gutenberg–Richter distributions for seven seismotectonic source zones. Seismic coupling seems to be large or even complete in most zones. An exception is the southern Tyrrhenian thrust zone, where most of the African–European convergence is accommodated. Here aseismic deformation is estimated to range from at least 25% along the western part to almost 100% aseismic slip around the Aeolian Islands. Even so, seismic potential of this zone has previously been significantly underestimated, due to the low levels of recorded past seismicity. By contrast, the series of 19 M6–7 earthquakes that hit Calabria in the 18th and 19th century released tectonic strain rates accumulated over time spans up to several times the catalog duration, and seismic potential is revised downward. The southern Tyrrhenian thrust zone and the extensional Calabrian faults, as well as the northeastern Sicilian transtensional zone between them (which includes the Messina Straits, where a destructive M7 event occurred in 1908), all have a similar seismic potential with minimum recurrence times of M ≥ 6.5 of 150–220 years. This potential is lower than that of the Southern Apennines (M ≥ 6.5 recurring every 60 to 140 years), but higher than that of southeastern Sicily (minimum M ≥ 6.5 recurrence times of 400 years). The high seismicity levels recorded in southeastern Sicily indicate some clustering and are most compatible with a tectonic scenario where the Ionian deforms internally, and motions at the Calabrian Trench are small. The estimated seismic potential for the Calabrian Trench and Central and Western Sicily are the lowest (minimum M ≥ 6.5 recurrence times of 550–800 years). Most zones are probably capable of generating earthquakes up to magnitudes 7–7.5, with the exception of Central and Western Sicily where maximum events sizes most likely do not exceed 7.     

Plate subduction,and generation of earthquakes and magmas in Japan as inferred from seismic observations: An overview

A dense nationwide seismic network recently constructed in Japan has been yielding large volumes of high-quality data that have made it possible to investigate the seismic structure in the Japanese subduction zone with unprecedented resolution. In this article, recent studies on the subduction of the Philippine Sea and Pacific plates beneath the Japanese Islands and the mechanism of earthquake and magma generation associated with plate subduction are reviewed. Seismic tomographic studies have shown that the Philippine Sea plate subducting beneath southwest Japan is continuous throughout the entire region, from Kanto to Kyushu, without disruption or splitting even beneath the Izu Peninsula as suggested in the past. The contact of the Philippine Sea plate with the Pacific plate subducting below has been found to cause anomalously deep interplate and intraslab earthquake activity in Kanto. Detailed waveform inversion studies have revealed that the asperity model is applicable to interplate earthquakes. Analyses of dense seismic and GPS network data have confirmed the existence of episodic slow slip accompanied in many instances by low-frequency tremors/earthquakes on the plate interface, which are inferred to play an important role in stress loading at asperities. High-resolution studies of the spatial variation of intraslab seismicity and the seismic velocity structure of the slab crust strongly support the dehydration embrittlement hypothesis for the generation of intraslab earthquakes. Seismic tomography studies have shown that water released by dehydration of the slab and secondary convection in the mantle wedge, mechanically induced by slab subduction, are responsible for magma generation in the Japanese islands. Water of slab origin is also inferred to be responsible for large anelastic local deformation of the arc crust leading to inland crustal earthquakes that return the arc crust to a state of spatially uniform deformation.     

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.     

Seismic activity changes progressing simultaneously with slow-slip in the Tokai area

Anomalous movements were detected simultaneously in both the seismic and the GPS observations in the Tokai area, the central part of the Japanese islands from the late 1990s to 2000. The anomalies are of great concern since the pending risk of a large megathrust earthquake in this area has been predicted for more than 20 years. The GPS data revealed that a slow-slip on the plate interface had commenced beneath Lake Hamana, the center of which is positioned around the edge of the assumed focal zone. On the other hand, the seismic data indicated that a delicate but clear quiescence appeared over a wide area that spreads into the main focal zone. Analyses of the seismicity changes in space and time confirmed that the contrast in the seismicity rate is distinct inside the focal zone. While the integrated seismicity indicated quiescence, some locations were distinguished as activated zones, possibly indicating the appearance of asperities. The rise of the seismicity rate in a quasi-stationary manner suggests an increase in the stress rate at that location. The following hypothesis is proposed based on the simultaneously detected evidences. The slow-slip progressing beneath Lake Hamana will induce a stress shift that invades the interior of the main locked zone, which will increase the contrast of the seismicity rate, possibly reflecting inhomogeneity in the locking strength. Even in this stage, the activated zones still maintain a locked state to prevent overall breakage. Investigations of the b-value changes and of tidal dependence in seismicity that reveal the stress-concentrated state also support the hypothesis. If this is the case, the observed change in seismicity would indicate the process of stress redistribution in the locking state, which represents the preparatory process toward final breakage. Tracking such seismicity changes would yield valid information for predictions of the next Tokai earthquake.     

Real time earthquake forecasting in Italy

We have applied an earthquake clustering epidemic model to real time data at the Italian Earthquake Data Center operated by the Istituto Nazionale di Geofisica e Vulcanologia (INGV) for short-term forecasting of moderate and large earthquakes in Italy. In this epidemic-type model every earthquake is regarded, at the same time, as being triggered by previous events and triggering following earthquakes. The model uses earthquake data only, with no explicit use of tectonic, geologic, or geodetic information. The forecasts are displayed as time-dependent maps showing both the expected rate density of M_l ≥ 4.0 earthquakes and the probability of ground shaking exceeding Modified Mercalli Intensity VI (PGA ≥ 0.01 g) in an area of 100 × 100 km² around the zone of maximum expected rate density in the following 24 h. For testing purposes, the overall probability of occurrence of an M_l ≥ 4.5 earthquake in the same area of 100 × 100 km² is also estimated. The whole procedure is tested in real time, for internal use only, at the INGV Earthquake Data Center.Forecast verification procedures have been carried out in forward-retrospective way on the 2006–2007 INGV data set, making use of statistical tools as the Relative Operating Characteristics (ROC) diagrams. These procedures show that the clustering epidemic model performs up to several hundred times better than a simple random forecasting hypothesis. The seismic hazard modeling approach so developed, after a suitable period of testing and refinement, is expected to provide a useful contribution to real time earthquake hazard assessment, even with a possible practical application for decision making and public information.     

Aftershock sequences in the mid-ocean ridge environment: an analysis using hydroacoustic data

Hydroacoustic data from autonomous arrays and the U.S. Navy's Sound Surveillance System (SOSUS) provide an opportunity to examine the temporal and spatial properties of seismicity along portions of the slow-spreading Mid-Atlantic Ridge (MAR), intermediate-spreading Juan de Fuca Ridge (JdFR) and fast-spreading East Pacific Rise (EPR). Aftershock and foreshock events are selected from the hydroacoustic earthquake catalog using single-link cluster (SLC) analysis, with a combined space–time metric. In the regions examined, hydroacoustic data improve the completeness level of the earthquake catalog by 1.5–2.0 orders of magnitude, allowing the decay constant, p, of the modified Omori law (MOL) to be determined for individual sequences. A non-parametric goodness-of-fit test indicates six of the seven sequences examined are described well by a MOL model. The p-values obtained for individual ridge and transform sequences using hydroacoustic data are larger than that previously estimated from the analysis of a stacked sequence generated from teleseismic data. For three sequences along the Siqueiros, Discovery and western Blanco Transforms, p-values are estimated to be 0.94–1.29. The spatial distribution of aftershocks suggests that the mainshock rupture is constrained by intra-transform spreading centers at these locations. An aftershock sequence following a 7.1M_s thrust event near the northern edge of the Easter Microplate exhibits p=1.02±0.11. Within the sequence, aftershocks are located to the north of a large topographic ridge, which may represent the surface expression of the shallow-dipping fault that ruptured during the mainshock. Two aftershock sequences near 24°25′N and 16°35′N on the MAR exhibit higher p-values, 1.74±0.23 and 2.37±1.65, although the latter estimate is not well constrained because of the small number of aftershocks. Larger p-values along the ridge crest might reflect a hotter thermal regime in this setting. Additional monitoring, however, will be needed to determine if p-value differences between the ridge and transform sequences are robust. A 1999 sequence on the Endeavour segment of the JdFR, which has been correlated with changes in the hydrothermal system, is described poorly by the MOL model. The failure of the MOL model, the anomalously large number of earthquakes within the sequence and absence of a clearly dominant mainshock are inconsistent with aftershock activity and the simple tectonic origin that has been proposed previously for this sequence.     

Properties of the aftershocks sequences of the 2000 and 2002 earthquakes in Aksehir-Afyon graben,west-central Anatolia,Turkey

The b-value of the Gutenberg–Richter’s frequency–magnitude relation and the p-value of the modified Omori law, which describes the decay rate of aftershock activity, were investigated for more than 500 aftershocks in the Aksehir-Afyon graben (AAG) following the 15 December 2000 Sultandagi–Aksehir and the 3 February 2002 Çay–Eber and Çobanlar earthquakes. We used the Kandilli Observatory’s catalog, which contains records of aftershocks with magnitudes ≥2.5. For the Çobanlar earthquake, the estimated b-values for three aftershock sequences are in the range 0.34 ≤  b ≤ 2.85, with the exception of the one that occurred during the first hour (4.77), while the obtained p-values are in the range 0.44 ≤ p ≤ 1.77. The aftershocks of the Sultandagi earthquake have a high p-value, indicating fast decay of the aftershock activity. A regular increase of b can be observed, with b < 1.0 after 0.208 days for the Çay–Eber earthquake. A systematic and similar increase and decrease pattern exists for the b- and p-values of the Çobanlar earthquakes during the first 5 days.     

Usage of an early warning and information system Web-site for real-time seismicity in Iceland

Iceland has been subjected to destructive earthquakes and volcanic eruptions throughout history. Such events are often preceded by changes in earthquake activity over varying timescales. Although most seismicity is confined to micro-earthquakes, large earthquakes have occurred within populated regions. Following the most recent hazardous earthquakes in 2000, the Icelandic Meteorological Office (IMO) developed an early warning and information system (EWIS) Web-site for viewing near-real-time seismicity in Iceland. Here we assess Web-site usage data in relation to earthquake activity, as recorded by the South Iceland Lowland (SIL) seismic network. Between March 2005 and May 2006 the SIL seismic network recorded 12,583 earthquakes. During this period, the EWIS Web-site logged a daily median of 91 visits. The largest onshore event (M _L 4.2) struck 20 km from Reykjavík on 06 March 2006 and was followed by an immediate, upsurge in usage resulting in a total of 1,173 unique visits to the Web-site. The greatest cluster of large (≥M _L 3) events occurred 300 km offshore from Reykjavík in May 2005. Within this swarm, 9 earthquakes ≥M _L 3 were detected on 11 May 2005, resulting in the release of a media bulletin by IMO. During the swarm, and following the media bulletin, the EWIS Web-site logged 1,234 unique visits gradually throughout the day. In summary, the data reveal a spatial and temporal relationship between Web-site usage and earthquake activity. The EWIS Web-site is accessed immediately after the occurrence of a local earthquake, whereas distant, unfelt earthquakes generate gradual interest prompted by media bulletins and, possibly, other contributing factors. We conclude that the Internet is a useful tool for displaying seismic information in near-real-time, which has the capacity to help increase public awareness of natural hazards.     

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.     

Verification test of the mean event size model for moderate earthquakes in the Kanto region, central Japan

In the Kanto region, there have been reports of decreasing b-values prior to earthquakes (M ≥ 5.5). The change of b-values is defined as the difference between the long- and short-term average earthquake magnitudes. A hazard function for moderately large earthquakes has also been proposed. This model was based on earthquakes that occurred between 1982 and 1999, and its effectiveness is measured retrospectively. In order to evaluate the effectiveness of the model through the use of independent data, the verification test is based on earthquakes observed after January 2000. Through the end of 2004, the only event where there was a decrease in mean event size was observed on June 3, 2000. This decrease resulted in a log-likelihood for the proposed model 1.3 units larger than that of the Poisson model, supporting the validity of the proposed model. Without accumulating further examples, we attempted to improve the verification test by expanding the study volume and by lowering the target cutoff magnitude in order to overcome the small sample size. When two other targets from the expanded volume were added, the difference in the log-likelihood (ΔlnL) increased to 3.6. In this case, the information rate per event was about 1.2, larger than that of the model period. Lowering the cutoff magnitude increased ΔlnL to 2.5. These extended tests led to higher confidence in the model with the larger ΔlnL value than did the primary test. From the viewpoint of the per-event information rate, each case involving targets with a magnitude 5.5 and larger resulted in better performance than in the model period.     

Statistical occurrence analysis and spatio-temporal distribution of earthquakes in the Apennines (Italy)

We present two examples of statistical analysis of seismicity conducted by integrating geological, geophysical and seismological data with the aim to characterize the active stress field and to define the spatio-temporal distribution of large earthquakes. Moreover, our data will help to improve the knowledge of the “seismogenic behavior” of the areas and to provide useful information for seismic hazard evaluation.The earthquakes are described by two non-parametric statistical procedures integrating also tectonic-physical parameters to study the spatio-temporal variability.The results show that the areas are characterized by: 1) a stress regime with mainly extensional kinematics; 2) tectonic structures mainly oriented with the active stress field (S_hmin = N44° ± 18° in the southern Apennines and S_hmin = N50° ± 17° in the central Apennines); 3) cluster distribution of seismicity and 4) a high probability of earthquake occurrence (M > 5.5) in the next 10 years.