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.     

Current accelerating seismic excitation along the northern boundary of the Aegean microplate

According to previous observations [Geophys. Res. Lett. 27 (2000) 3957], the generation of large (M≥7.0) earthquakes in the western part of the north Anatolian fault system (Marmara Sea) is followed by strong earthquakes along the Northern Boundary of the Aegean microplate (NAB: northwestermost Anatolia–northern Aegean–central Greece–Ionian islands). Therefore, it can be hypothesized that a seismic excitation along this boundary should be expected after the occurrence of the Izmit 1999 earthquake (M=7.6). We have applied the method of accelerating seismic crustal deformation, which is based on concepts of critical point dynamics in an attempt to locate more precisely those regions along the NAB where seismic excitation is more likely to occur. For this reason, a detailed parametric grid search of the broader NAB area was performed for the identification of accelerating energy release behavior.Three such elliptical critical regions have been identified with centers along this boundary. The first region, (A), is centered in the eastern part of this boundary (40.2°N, 27.2°E: southwest of Marmara), the second region, (B), has a center in the middle part of the boundary (38.8°N, 23.4°E: East Central Greece) and the third region, (C), in the westernmost part of the boundary (38.2°N, 20.9°E: Ionian Islands). The study of the time variation of the cumulative Benioff strain in two of the three identified regions (A and B) revealed that intense accelerating seismicity is observed especially after the occurrence of the 1999 Izmit mainshock. Therefore, it can be suggested that the seismic excitation, at least in these two regions, has been triggered by the Izmit mainshock.Estimations of the magnitudes and origin times of the expected mainshocks in these three critical regions have also been performed, assuming that the accelerating seismicity in these regions will lead to a critical point, that is, to the generation of mainshocks.     

Seismic behaviour of the North Anatolian Fault beneath the Sea of Marmara (NW Turkey): implications for earthquake recurrence times and future seismic hazard

Possible long-term seismic behaviour of the Northern strand of the North Anatolian Fault Zone, between western extreme of the 1999 İzmit rupture and the Aegean Sea, after 400 AD is studied by examining the historical seismicity, the submarine fault mapping and the paleoseismological studies of the recent scientific efforts. The long-term seismic behaviour is discussed through two possible seismicity models devised from M _S ≥ 7.0 historical earthquakes. The estimated return period of years of the fault segments for M1 and M2 seismic models along with their standard deviations are as follows: F4 segment 255 ± 60 and 258 ± 12; F5 segment 258 ± 60 and 258 ± 53; F6 segment 258 ± 60 and 258 ± 53; F7 segment 286 ± 103 and 286 ± 90; F8 segment 286 ± 90 and 286 ± 36. As the latest ruptures on the submarine segments have been reported to be during the 1754–1766 earthquake sequence, and the 1912 mainshock rupture has been evidenced to extend almost all over the western part of the Sea of Marmara, our results imply imminent seismic hazard and, considering the mean recurrence time, a large earthquake to strike the eastern part of the Sea of Marmara in the next two decades.     

A study of microearthquake seismicity and focal mechanisms within the Sea of Marmara (NW Turkey) using ocean bottom seismometers (OBSs)

We have carried out seismological observations within the Sea of Marmara (NW Turkey) in order to investigate the seismicity induced after Gölcük–İzmit (Kocaeli) earthquake (M_w 7.4) of August 17, 1999, using ocean bottom seismometers (OBSs). High-resolution hypocenters and focal mechanisms of microearthquakes have been investigated during this Marmara Sea OBS project involving deployment of 10 OBSs within the Çınarcık (eastern Marmara Sea) and Central-Tekirdağ (western Marmara Sea) basins during April–July 2000. Little was known about microearthquake activity and their source mechanisms in the Marmara Sea. We have detected numerous microearthquakes within the main basins of the Sea of Marmara along the imaged strands of the North Anatolian Fault (NAF). We obtained more than 350 well-constrained hypocenters and nine composite focal mechanisms during 70 days of observation. Microseismicity mainly occurred along the Main Marmara Fault (MMF) in the Marmara Sea. There are a few events along the Southern Shelf. Seismic activity along the Main Marmara Fault is quite high, and focal depth distribution was shallower than 20 km along the western part of this fault, and shallower than 15 km along its eastern part. From high-resolution relative relocation studies of some of the microearthquake clusters, we suggest that the western Main Marmara Fault is subvertical and the eastern Main Marmara Fault dips to south at 45°. Composite focal mechanisms show a strike-slip regime on the western Main Marmara Fault and complex faulting (strike-slip and normal faulting) on the eastern Main Marmara Fault.     

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.     

Geodetic and seismological investigation of crustal deformation near Izmir (Western Anatolia)

The Aegean region including western Turkey, mainland Greece, and the Hellenic Arc is the most seismological and geodynamical active domain in the Alpine Himalayan Belt. In this study, we processed 3 years of survey-mode GPS data and present the analysis of a combination of geodetic and seismological data around Izmir, which is the third most populated city in Turkey. The velocities obtained from 15 sites vary between 25 mm/yr and 28 mm/yr relative to the Eurasian plate. The power law exponent of earthquake size distribution (b-value) ranges from 0.8 to 2.8 in the Izmir region between 26.2°E and 27.2°E. The lowest b-value zones are found along Karaburun Fault (b = 0.8) and, between Seferihisar and Tuzla Faults (b = 0.8). A localized stress concentration is expected from numerical models of seismicity along geometrical locked fault patches. Therefore, areas with lowest b-values are considered to be the most likely location for a strong earthquake, a prediction that is confirmed by the 2005 M_w = 5.9 Seferihisar earthquake sequences, with epicentres located to the south of the Karaburun Fault. The north–south extension of the Izmir area is corroborated by extension rates up to 140 nanostrain/yr as obtained from our GPS data. We combined the 3-year GPS velocity field with the published velocity field to determine the strain rate pattern in the area. The spatial distribution of b-value reflects the normal background due to the tectonic framework and is corroborated by the geodetic data. b-Values correlate with strain pattern. This relationship suggests that decrease of b-values signifies accumulating strain.     

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.     

Characteristics of the 1912 co-seismic rupture along the North Anatolian Fault Zone (Turkey): implications for the expected Marmara earthquake

Evidence of right‐lateral offsets associated with the 1912 earthquake (M_w 7.4) along the North Anatolian Fault (Gaziköy–Saros segment) allow us to survey (using DGPS) the co‐seismic and cumulative slip distribution. The damage distribution and surface breaks related with the earthquake show an elongated zone of maximum intensity (X MSK) parallel to the fault rupture on land but this may extend offshore to the north‐east and south‐west. Detailed mapping of the fault using topographic maps and aerial photographs indicates the existence of pull‐apart basins and pressure ridges. At several localities, the average 1912 offset along strike is 3.5–4 m and cumulative slip is 2–6 times that of individual movement. The fault rupture geometry and slip distribution suggest the existence of three subsegments with a combined total length of 110–120 km, a fault length and maximum slip similar to those of the 1999 Izmit earthquake. The amount of slip at the north‐easternmost section and in the coastal region of the Sea of Marmara reaches an average 4 m, thereby implying the offshore extension of the 1912 rupture. The results suggest that the 1912 event generated up to 150 km of surface faulting, which would imply a M_w 7.2–7.4 earthquake and which, added with rupture lengths of the 1999 earthquakes, help to constrain the remaining seismic gap in the Sea of Marmara.     

On the swelling potential of compacted high plasticity clays

We conducted a study of the spatial distributions of seismicity and earthquake hazard parameters for Turkey and the adjacent areas, applying the maximum likelihood method. The procedure allows for the use of either historical or instrumental data, or even a combination of the two. By using this method, we can estimate the earthquake hazard parameters, which include the maximum regional magnitude Mˆ_max, the activity rate of seismic events and the well-known bˆ value, which is the slope of the frequency-magnitude Gutenberg-Richter relationship. These three parameters are determined simultaneously using an iterative scheme. The uncertainty in the determination of the magnitudes was also taken into consideration. The return periods (RP) of earthquakes with a magnitude M ≥ m are also evaluated. The whole examined area is divided into 24 seismic regions based on their seismotectonic regime. The homogeneity of the magnitudes is an essential factor in such studies. In order to achieve homogeneity of the magnitudes, formulas that convert any magnitude to an M_S-surface scale are developed. New completeness cutoffs and their corresponding time intervals are also assessed for each of the 24 seismic regions. Each of the obtained parameters is distributed into its respective seismic region, allowing for an analysis of the localized seismicity parameters and a representation of their regional variation on a map. The earthquake hazard level is also calculated as a function of the form Θ = (Mˆ_max,RP_6.0), and a relative hazard scale (defined as the index K) is defined for each seismic region. The investigated regions are then classified into five groups using these parameters. This classification is useful for theoretical and practical reasons and provides a picture of quantitative seismicity. An attempt is then made to relate these values to the local tectonics.     

Analysis of the characteristics of self-similarity of seismicity and the active fault network of Eurasia

An analysis of the possible relationship between fractal dimensions of the active fault network, spatial distribution of earthquake epicenters, and parameter b in the Gutenberg-Richter law is presented. The quantitative characteristics of self-similarity of the seismic process and the active fault network of seismically active areas of Eurasia are obtained. This self-similarity manifests itself over a range of at least two orders of spatial scales and magnitudes. The obtained estimations of the fractal dimensions of the fault network D _f and epicenter field D _e are close for all the areas analyzed. It is established that the average value connecting values D and b for all the investigated areas is slightly higher than the theoretical value (2.0) and varies within the range of 1.7–2.4.     

Quantitative finite strain and kinematic flow analyses along the Zagros transpression zone, Iran

The Ghouri area in southwest Iran exposes a cross section through the Zagros orogenic belt. The area provides an opportunity to investigate quantitative finite strain (Rs), kinematic vorticity number (W_k), proportions of pure shear and simple shear components, sense of shear indicators, steeply plunging lineations, and other moderate to steeply plunging stretching lineations in a transpressional zone. Based on a classical strain analysis of deformed microfossils with oblate strain ellipsoid shape, the Zagros orogenic belt is classified as a pure-shear dominated zone of transpression, but asymmetry of shear-sense indicators suggests that a significant component of simple shear was involved along the deformation zone boundaries. The long axes of the microfossils and stretched pebbles of a deformed conglomerate were used to indicate the stretching direction in this zone. The stretching lineations have a steep to moderate plunge but a constant strain magnitude. Characteristics of dextral inclined transpressional kinematics in the Zagros continental collision zone were quantified and indicate an estimated k-value < 1, an angle between the maximum horizontal axis of the instantaneous strain ellipsoid and the zone boundary (θ = 32°), asymmetrical dextral shear-sense indicators, and an angle of relative plate motion (α = 25°).     

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.     

Maximum Likelihood Estimation of Seismic Hazard for Sweden

The maximum magnitude, the activity rate, and the Gutenberg-Richterb parameter as earthquake hazard parameters, have been evaluated for Sweden. The maximum likelihood method permits the combination of historical and instrumental data. The catalog used consists of 1100 earthquakes in the time interval 1375–1989. The extreme part of the catalog contains only the strongest historical earthquakes, whereas the complete part is divided into several subcatalogs, each assumed complete above a specified threshold magnitude. The uncertainty in magnitude determination was taken into account. For southern Sweden, the calculations giveb-values of 1.04 (0.05) for the whole area south of 60° N and 0.98 (0.06) for a subregion of enhanced seismicity in the Lake Vänern area. For the whole area north of 60° N, theb-value is 1.35 (0.06) and for the seismicity zone along the Gulf of Bothnia 1.26 (0.06). The number of annually expected earthquakes with magnitude equal to or larger than 2.4 [M_L(UPP) or M_M(UPP)] is 1.8 for the whole southern Sweden, 1.3 for the Lake Vänern region, 3.7 for northern Sweden, and 2.4 for the region along the Gulf of Bothnia. The maximum expected regional magnitude is calculated to 4.9 (0.5) for a time span of 615 years for southern Sweden and the Lake Vänern subregion, and 4.3 (0.5) for a time span of 331 years for northern Sweden and the Gulf of Bothnia subregion. However, several historical earthquakes with magnitude above 5 in nearby areas of Norway indicate that the seismic potential may be higher.     

Geochronological problems related to polymetamorphism in the Limpopo Complex, South Africa

The integration of new and published geochronologic data with structural, magmatic/anatectic and pressure–temperature (P–T) process information allow the recognition of high-grade polymetamorphic granulites and associated high-grade shear zones in the Central Zone (CZ) of the Limpopo high-grade terrain in South Africa. Together, these two important features reflect a major high-grade D₃/M₃ event at ~ 2.02 Ga that overprinted the > 2.63 Ga high-grade Neoarchaean D₂/M₂ event, characterized by SW-plunging sheath folds. These major D₂/M₂ folds developed before ~ 2.63 Ga based on U–Pb zircon age data for precursors to leucocratic anatectic gneisses that cut the high-grade gneissic fabric. The D₃/M₃ shear event is accurately dated by U–Pb monazite (2017.1 ± 2.8 Ma) and PbSL garnet (2023 ± 11 Ma) age data obtained from syntectonic anatectic material, and from sheared metapelitic gneisses that were completely reworked during the high-grade shear event. The shear event was preceded by isobaric heating (P = ~ 6 kbar and T = ~ 670–780 °C), which resulted in the widespread formation of polymetamorphic granulites. Many efforts to date high-grade gneisses from the CZ using PbSL garnet dating resulted in a large spread of ages (~ 2.0–2.6 Ga) that reflect the polymetamorphic nature of these complexly deformed high-grade rocks.     

The 1998 Hida Mountain, Central Honshu, Japan, earthquake swarm: Double-difference event relocation, frequency–magnitude distribution and Coulomb stress changes

By using the double-difference relocation technique, we have determined the fine structure of seismicity during the 1998 Hida Mountain earthquake swarm. The distribution of seismic activity defines two main directions (N–S and E–W) that probably correspond to the regional stress pattern. The detailed structure of seismicity reveals intense spatio-temporal clustering and earthquake lineations. Each cluster of events contains a mainshock and subsequent aftershock activity that decays according to the Omori law. The seismicity and the b-value temporal and spatial patterns reflect the evolution of the static stress changes during the earthquake swarm. About 80% of the swarm's best-relocated events occur in regions of increased ΔCFF. The smaller value of b found in the northern part of the swarm region and a larger b-value observed to the south, for the same period of time, could be well explained by the static stress changes caused by the larger events of the sequence. We argue that the state of stress in the crust is the main factor that controls the variation of b-value.     

Spatial variation of aftershock activity across the rupture zone of the 17 August 1999 Izmit earthquake, Turkey

We examined the spatial variation in the aftershock activity from the 17 August 1999 Izmit, Turkey earthquake. We found that this aftershock sequence is non-uniform both in space and time, aspects that need to be taken into account in any further statistical analysis. Other aspects of this aftershock sequence are similar to other aftershock sequences, namely low b-values and a high degree of spatial variation. We have detected three zones of relatively high b-values, two of which coincide with asperities revealed by previous slip inversion studies. The third zone with an anomalous b-value is located beyond the fault rupture and indicates a weakened fractured zone in the Yalova-Tuzla area. This b-value analysis provided no evidence for any significant difference that may exist between the two sides of the mainshock fault plane.     

Seismic source characteristics in Kachchh and Saurashtra regions of Western India: b-value and fractal dimension mapping of aftershock sequences

Seismic source characteristics in the Kachchh rift basin and Saurashtra horst tectonic blocks in the stable continental region (SCR) of western peninsular India are studied using the earthquake catalog data for the period 2006–2011 recorded by a 52-station broadband seismic network known as Gujarat State Network (GSNet) running by Institute of Seismological Research (ISR), Gujarat. These data are mainly the aftershock sequences of three mainshocks, the 2001 Bhuj earthquake (M _w 7.7) in the Kachchh rift basin, and the 2007 and 2011 Talala earthquakes (M _w ≥ 5.0) in the Saurashtra horst. Two important seismological parameters, the frequency–magnitude relation (b-value) and the fractal correlation dimension (D _c) of the hypocenters, are estimated. The b-value and the D _c maps indicate a difference in seismic characteristics of these two tectonic regions. The average b-value in Kachchh region is 1.2 ± 0.05 and that in the Saurashtra region 0.7 ± 0.04. The average D _c in Kachchh is 2.64 ± 0.01 and in Saurashtra 2.46 ± 0.01. The hypocenters in Kachchh rift basin cluster at a depth range 20–35 km and that in Saurashtra at 5–10 km. The b-value and D _c cross sections image the seismogenic structures that shed new light on seismotectonics of these two tectonic regions. The mainshock sources at depth are identified as lower b-value or stressed zones at the fault end. Crustal heterogeneities are well reflected in the maps as well as in the cross sections. We also find a positive correlation between b- and D _c-values in both the tectonic regions.

     

A comparison of seismic and structural measurements of scaling exponents during tensile subcritical crack growth

The observed fractal nature of both fault length distributions and earthquake magnitude-frequency distributions suggests that there may be a relationship between the structure of active fault systems and the resulting seismicity. In previous theoretical work, a positive correlation between the exponent D from the fracture length distribution, and the seismic or acoustic emission (AE) b-value has been inferred from a simple dislocation model of the seismic source. Here, we present the first experimental evidence for a correlation between D and b from a series of tensile fracture mechanics tests on crystalline rock, carried out in different environmental conditions, both air-dry and water-saturated, and at ambient temperature and pressure. The microseismic acoustic emissions were monitored during subcritical crack growth under controlled conditions of constant stress intensity, K_I, and quantitative analyses of the resulting fracture patterns were carried out on the same specimens. It is found that AE b-values, ranging from 1.0 to 2.3, correlate negatively with the normalized stress intensity K_I/KIC, where K_IC is the fracture toughness of the specimen. The microcrack length distribution exponent D, ranges from 1.0 to 1.7. Fluid presence has a first-order influence on both the AE and structure produced in these experiments. For experiments at low stress intensity or high fluid content, the activation of the stress corrosion mechanism for K_I < K_IC leads to a greater relative proportion both of small cracks and of low amplitude acoustic emissions, reflected in higher values of D and b. The exponent D is found to correlate positively with the AE b-value.     

Paleomagnetism of the middle Cretaceous Iritono granite in the Abukuma region, northeast Japan

We have studied the paleomagnetism of the middle Cretaceous Iritono granite of the Abukuma massif in northeast Japan together with ⁴⁰Ar–³⁹Ar dating. Paleomagnetic samples were collected from ten sites of the Iritono granite (102 Ma ⁴⁰Ar–³⁹Ar age) and two sites of its associated gabbroic dikes. The samples were carefully subjected to alternating field and thermal demagnetizations and to rock magnetic analyses. Most of natural remanent magnetizations show mixtures of two components: (1) H component, high coercivity (B_c > 50–90 mT) or high blocking temperature (T_b > 350–560 °C) component and (2) L component, relatively low B_c or low T_b component. H component was obtained from all the 12 sites to give a mean direction of shallow inclination and northwesterly declination (I = 29.9°, D = 311.0°, α₉₅ = 2.7°, N = 12). This direction is different from the geocentric axial dipole field at the present latitude (I = 56.5°) and the typical direction of the Cenozoic remagnetization in northeast Japan. Since rock magnetic properties indicate that the H component of the Iritono granite is carried mainly by magnetite inclusions in plagioclase, this component probably retains a primary one. Thus the shallow inclination indicates that the Abukuma massif was located at a low latitude (16.1 ± 1.6°N) about 100 Ma and then drifted northward by about 20° in latitude. The northwesterly deflection is attributed mostly to the counterclockwise rotation of northeast Japan due to Miocene opening of the Japan Sea. According to this model, the low-pressure and high-temperature (low-P/high-T) metamorphism of the Abukuma massif, which has been well known as a typical location, would have not occurred in the present location. On the other hand, the L component is carried mainly by pyrrhotite and its mean direction shows a moderate inclination and a northwesterly declination (I = 42.8°, D = 311.5°, α₉₅ = 3.3°, N = 9). Since this direction is intermediate between the H component and early Cenozoic remagnetization in northeast Japan, some thermal event would have occurred at lower temperature than pyrrhotite Curie point ( 320 °C) during the middle Cretaceous to early Cenozoic time to have resulted in partial remagnetization.     

Comment on “Structural controls on a carbon dioxide-driven mud volcano field in the Northern Apennines (Pieve Santo Stefano, Italy): Relations with pre-existing steep discontinuities and seismicity”

Bonini (2009, Structural controls on a carbon dioxide-driven mud volcano field in the Northern Apennines (Pieve Santo Stefano, Italy): relations with pre-existing steep discontinuities and seismicity. Journal of Structural Geology 31, 44–54) presents a 2D mechanical analysis to infer the failure conditions responsible for the seismicity distribution during an M_w = 4.6 seismic sequence nucleating during 2001 in the Northern Apennines. In my view the mechanical analysis presented in this paper has some weakness or is not well constrained, in particular: 1) the assumption of a dip angle of 50°, is not consistent with the activated structures; 2) the P_f = σ₃ condition, difficult to be attained along a cohesionless fault dipping at 50°; 3) the isotropic stress state, i.e. σ₂ = σ₃, that is not consistent with the active or recent stress field in the area.