Two-dimensional spatial analysis of the seismic b-value and the Bouguer gravity anomaly in the southeastern part of the Zagros Fold-and-Thrust Belt,Iran: Tectonic implications

Environmental managers and protection agencies try to assess the magnitudes of earthquakes in regions of seismic activity. For several decades they have used the seismic b-values and Bouguer anomalies for evaluating the crustal character and stress regimes. We have analyzed geostatistically data on both variables to map their spatial distributions in the southeast of the Zagros of Iran. We found a strong correlation between the distribution of the b-value and the Bouguer gravity anomaly in the region. The large Bouguer gravity anomaly values and small b-values all accord with there being a thinner crustal root and a larger concentration of stress in the center. The small to moderate Bouguer gravity anomaly values and intermediate to large b-values accord with the thicker crustal root and the smaller concentration of stress in the northeast. We conclude the southeast of the Zagros, consists of heterogeneous crust, such that accounts for its varied tectonics.     

Crustal shear-wave velocity structure beneath northeast India from teleseismic receiver function analysis

We investigated the seismic shear-wave velocity structure of the crust beneath nine broadband seismological stations of the Shillong–Mikir plateau and its adjoining region using teleseismic P-wave receiver function analysis. The inverted shear wave velocity models show ∼34–38 km thick crust beneath the Shillong Plateau which increases to ∼37–38 km beneath the Brahmaputra valley and ∼46–48 km beneath the Himalayan foredeep region. The gradual increase of crustal thickness from the Shillong Plateau to Himalayan foredeep region is consistent with the underthrusting of Indian Plate beyond the surface collision boundary. A strong azimuthal variation is observed beneath SHL station. The modeling of receiver functions of teleseismic earthquakes arriving the SHL station from NE backazimuth (BAZ) shows a high velocity zone within depth range 2–8 km along with a low velocity zone within ∼8–13 km. In contrast, inversion of receiver functions from SE BAZ shows high velocity zone in the upper crust within depth range ∼10–18 km and low velocity zone within ∼18–36 km. The critical examination of ray piercing points at the depth of Moho shows that the rays from SE BAZ pierce mostly the southeast part of the plateau near Dauki fault zone. This observation suggests the effect of underthrusting Bengal sediments and the underlying oceanic crust in the south of the plateau facilitated by the EW-NE striking Dauki fault dipping 30⁰ toward northwest.     

Long wavelength gravity anomalies over India: Crustal and lithospheric structures and its flexure

Long wavelength gravity anomalies over India were obtained from terrestrial gravity data through two independent methods: (i) wavelength filtering and (ii) removing crustal effects. The gravity fields due to the lithospheric mantle obtained from two methods were quite comparable. The long wavelength gravity anomalies were interpreted in terms of variations in the depth of the lithosphere–asthenosphere boundary (LAB) and the Moho with appropriate densities, that are constrained from seismic results at certain points. Modeling of the long wavelength gravity anomaly along a N–S profile (77°E) suggest that the thickness of the lithosphere for a density contrast of 0.05 g/cm³ with the asthenosphere is maximum of ∼190 km along the Himalayan front that reduces to ∼155 km under the southern part of the Ganga and the Vindhyan basins increasing to ∼175 km south of the Satpura Mobile belt, reducing to ∼155–140 km under the Eastern Dharwar craton (EDC) and from there consistently decreasing south wards to ∼120 km under the southernmost part of India, known as Southern Granulite Terrain (SGT).The crustal model clearly shows three distinct terrains of different bulk densities, and thicknesses, north of the SMB under the Ganga and the Vindhyan basins, and south of it the Eastern Dharwar Craton (EDC) and the Southern Granulite Terrain (SGT) of bulk densities 2.87, 2.90 and 2.96 g/cm³, respectively. It is confirmed from the exposed rock types as the SGT is composed of high bulk density lower crustal rocks and mafic/ultramafic intrusives while the EDC represent typical granite/gneisses rocks and the basement under the Vindhyan and Ganga basins towards the north are composed of Bundelkhand granite massif of the lower density. The crustal thickness along this profile varies from ∼37–38 km under the EDC, increasing to ∼40–45 km under the SGT and ∼40–42 km under the northern part of the Ganga basin with a bulge up to ∼36 km under its southern part. Reduced lithospheric and crustal thicknesses under the Vindhyan and the Ganga basins are attributed to the lithospheric flexure of the Indian plate due to Himalaya. Crustal bulge due to lithospheric flexure is well reflected in isostatic Moho based on flexural model of average effective elastic thickness of ∼40 km. Lithospheric flexure causes high heat flow that is aided by large crustal scale fault system of mobile belts and their extensions northwards in this section, which may be responsible for lower crustal bulk density in the northern part. A low density and high thermal regime in north India north of the SMB compared to south India, however does not conform to the high S-wave velocity in the northern part and thus it is attributed to changes in composition between the northern and the southern parts indicating a reworked lithosphere. Some of the long wavelength gravity anomalies along the east and the west coasts of India are attributed to the intrusives that caused the breakup of India from Antarctica, and Africa, Madagascar and Seychelles along the east and the west coasts of India, respectively.     

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.     

Lithosphere,crust and basement ridges across Ganga and Indus basins and seismicity along the Himalayan front,India and Western Fold Belt,Pakistan

Spectral analysis of the digital data of the Bouguer anomaly of North India including Ganga basin suggest a four layer model with approximate depths of 140, 38, 16 and 7 km. They apparently represent lithosphere–asthenosphere boundary (LAB), Moho, lower crust, and maximum depth to the basement in foredeeps, respectively. The Airy’s root model of Moho from the topographic data and modeling of Bouguer anomaly constrained from the available seismic information suggest changes in the lithospheric and crustal thicknesses from ∼126–134 and ∼32–35 km under the Central Ganga basin to ∼132 and ∼38 km towards the south and 163 and ∼40 km towards the north, respectively. It has clearly brought out the lithospheric flexure and related crustal bulge under the Ganga basin due to the Himalaya. Airy’s root model and modeling along a profile (SE–NW) across the Indus basin and the Western Fold Belt (WFB), (Sibi Syntaxis, Pakistan) also suggest similar crustal bulge related to lithospheric flexure due to the WFB with crustal thickness of 33 km in the central part and 38 and 56 km towards the SE and the NW, respectively. It has also shown the high density lower crust and Bela ophiolite along the Chamman fault. The two flexures interact along the Western Syntaxis and Hazara seismic zone where several large/great earthquakes including 2005 Kashmir earthquake was reported.The residual Bouguer anomaly maps of the Indus and the Ganga basins have delineated several basement ridges whose interaction with the Himalaya and the WFB, respectively have caused seismic activity including some large/great earthquakes. Some significant ridges across the Indus basin are (i) Delhi–Lahore–Sargodha, (ii) Jaisalmer–Sibi Syntaxis which is highly seismogenic. and (iii) Kachchh–Karachi arc–Kirthar thrust leading to Sibi Syntaxis. Most of the basement ridges of the Ganga basin are oriented NE–SW that are as follows (i) Jaisalmer–Ganganagar and Jodhpur–Chandigarh ridges across the Ganga basin intersect Himalaya in the Kangra reentrant where the great Kangra earthquake of 1905 was located. (ii) The Aravalli Delhi Mobile Belt (ADMB) and its margin faults extend to the Western Himalayan front via Delhi where it interacts with the Delhi–Lahore ridge and further north with the Himalayan front causing seismic activity. (iii) The Shahjahanpur and Faizabad ridges strike the Himalayan front in Central Nepal that do not show any enhanced seismicity which may be due to their being parts of the Bundelkhand craton as simple basement highs. (iv) The west and the east Patna faults are parts of transcontinental lineaments, such as Narmada–Son lineament. (v) The Munghyr–Saharsa ridge is fault controlled and interacts with the Himalayan front in the Eastern Nepal where Bihar–Nepal earthquakes of 1934 has been reported. Some of these faults/lineaments of the Indian continent find reflection in seismogenic lineaments of Himalaya like Everest, Arun, Kanchenjunga lineaments. A set of NW–SE oriented gravity highs along the Himalayan front and the Ganga and the Indus basins represents the folding of the basement due to compression as anticlines caused by collision of the Indian and the Asian plates. This study has also delineated several depressions like Saharanpur, Patna, and Purnia depressions.     

Crustal structure of the western part of the Southern Granulite Terrain of Indian Peninsular Shield derived from gravity data

The Southern Granulite Terrain (SGT) is composed of high-grade granulite domain occurring to the south of Dharwar Craton (DC). The structural units of SGT show a marked change in the structural trend from the dominant north–south in DC to east–west trend in SGT and primarily consist of different crustal blocks divided by major shear zones. The Bouguer anomaly map prepared based on nearly 3900 gravity observations shows that the anomalies are predominantly negative and vary between −125 mGal and +22 mGal. The trends of the anomalies follow structural grain of the terrain and exhibit considerable variations within the charnockite bodies. Two-dimensional wavelength filtering as well as Zero Free-air based (ZFb) analysis of the Geoid-Corrected Bouguer Anomaly map of the region is found to be very useful in preparing regional gravity anomaly map and inversion of this map gave rise to crustal thicknesses of 37–44 km in the SGT. Crustal density structure along four regional gravity profiles cutting across major shear zones, lineaments, plateaus and other important geological structures bring out the following structural information. The Bavali Shear Zone extending at least up to 10 km depth is manifested as a plane separating two contrasting upper crustal blocks on both sides and the gravity high north of it reveals the presence of a high density mass at the base of the crust below Coorg. The steepness of the Moyar and Bhavani shears on either side of Nilgiri plateau indicates uplift of the plateau due to block faulting with a high density mass at the crustal base. The Bhavani Shear Zone is manifested as a steep southerly dipping plane extending to deeper levels along which alkaline and granite rocks intruded into the top crustal layer. The gravity high over Palghat gap is due to the upwarping of Moho by 1–2 km with the presence of a high density mass at intermediate crustal levels. The gravity low in Periyar plateau is due to the granite emplacement, mid-crustal interface and the thicker crust. The feeble gravity signature across the Achankovil shear characterized by sharp velocity contrast indicates that the shear is not a superficial structure but a crustal scale zone of deformation reaching up to mid-crustal level.     

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.     

Crustal structure and nature of emplacement of the 85°E Ridge in the Mahanadi offshore based on constrained potential field modeling: Implications for intraplate plume emplaced volcanism

An integrated interpretation of multi-channel seismic reflection, gravity and magnetic datasets belonging to northern most part of the 85°E Ridge in the Mahanadi offshore is carried out to study the crustal structure and mode of its emplacement. The basement structure map of the ridge reveals that it is 130–150 km wide and is composed of an eastern high which appears as a continuous, broad and smooth topographyand the western high characterized by several steep isolated highs. The seismic velocities reported for the first time over the ridge indicate several sedimentary sequences ranging in velocities between 1.6 and 4.0 km/s above the acoustic basement top. The salient aspects of the sedimentary velocities are; a low velocity layer (2.6–3.2 km/s) within the Cretaceous sequence in the intervening depressions encompassing the flank region, and a regionally widespread higher velocity layer (3.5–3.8 km/s) belonging to the Eocene–Oligocene section overlying the ridge. A layer having a velocity of 4.2–4.7 km/s probably made of volcanoclastic rocks is observed immediately below the acoustic basement. The sediment isopach maps presented here for three major horizons are used to compute the 3-D sediment gravity effect to obtain a crustal Bouguer anomaly map of the region. Detailed analysis of the gravity and magnetic anomaly maps clearly demonstrates the continuity of ridge up to the Mahanadi coast at Chilka Lake. Seismically constrained gravity and magnetic models indicate that the ridge is composed of volcanic material that was emplaced on continental crust in the shelf-slope areas and over the oceanic crust in the deep offshore areas. The modeled crustal structure below the ridge further indicates volcanic emplacement of the ridge on a relatively younger lithosphere. We propose two alternative models for the emplacement of the ridge.     

Gravity field and isostatic state of Ethiopia and adjacent areas

Over 35,000 onshore and offshore gravity stations have been compiled in order to test isostatic models against geologic structures over a part of the Afro–Arabian shield. The area of Ethiopia covers an important part of this system because it contains the major section of the ≈5000 km Afro–Arabian rift and includes the transition between the Arabo-Nubian-Shield (ANS) and the Mozambique Belt (MB).Isostatic residual anomalies have been calculated using both Airy and Vening-Meinesz (flexural rigidity D = 10²² Nm) models. The isostatic residual anomalies outline the major Precambrian belts, the Cenozoic rifts and associated major structures. Positive residual anomalies associated with the main Ethiopian Rift (MER) and Kenyan rift systems could be the expressions of an axial intrusive body and swarms of local faults and fractures. The residual anomalies indicate relative stability in the MER and increased tectonic activity in the areas of the Red Sea, Gulf of Aden and Afar. Near-zero isostatic residuals flank the MER and Kenya rifts and are found within the Danakil Alps and some plateau regions.The small mean isostatic residual anomaly (about 8 mGal) and the isostatic analysis show a slight positive bias indicating under compensation. The undercompensation may imply that there are upper crustal features that are not compensated regionally (probably supported by the rigidity of the lithosphere) and isostatic disequilibrium in the region. Therefore, the high topography of Ethiopia and East African plateau is partly compensated by thicker crust (broad negative isostatic regional anomaly) and partly by dynamic forces.The results of the qualitative interpretation form the basis of continuing three-dimensional gravity modelling and quantitative analysis that also integrates data from eastern Sudan.     

The 2007 Talala,Saurashtra, western India earthquake sequence: Tectonic implications and seismicity triggering

A damaging and widely felt moderate (M_w 5.0) earthquake occurred in the Talala region of Saurashtra, Gujarat (western India) on November 6, 2007. The highly productive sequence comprised about 1300 micro earthquakes (M > 0.5) out of which 325 of M ? 1.5 that occurred during November 6, 2007–January 10, 2008 were precisely located. The spatial aftershock distribution revealed a NE–SW striking fault in accordance with the centroid moment tensor solution, which in turn implies left-lateral motion. The orientation and sense of shear are consistent with similarly orientated geological fault identified in the area from satellite imagery and field investigation.The aftershocks temporal decay, b-value of frequency–magnitude distribution, spatial fractal dimension, D, and slip ratio (ratio of the slip occurred on the primary fault to the total slip) were examined with the purpose to identify the properties of the sequence. The high b-value (1.18 ± 0.01) may be attributed to the paucity of the larger (M ? 4.0) aftershocks and reveals crustal heterogeneity and low stress regime. The high p-value (1.10 ± 0.39), implying fast decay rate of aftershocks, evidences high surface heat flux. A value of the spatial fractal dimension (D) equal to 2.21 ± 0.02 indicates random spatial distribution and source in a two-dimensional plane that is being filled-up by fractures. A slip ratio of 0.42 reveals that more slip occurred on secondary fault systems.The static Coulomb stress changes due to the coseismic slip of the main shock, enhanced off fault aftershock occurrence. The occurrence of a moderate earthquake (M_w 4.3) on October 5, 2008 inside a region of positive Coulomb stress changes supports the postulation on aftershock triggering. When the stress changes were resolved on a cross section including the stronger (M4.8) foreshock plane that is positioned adjacent to the main fault, it became evident that the activity continued there due to stress transfer from the main rupture.     

Paleoproterozoic crustal evolution of the Tarim Craton: Constrained by zircon U–Pb and Hf isotopes of meta-igneous rocks from Korla and Dunhuang

The Tarim Craton is one of three large cratons in China. Presently, there is only scant information concerning its crustal evolutionary history because most of the existing geochronological studies have lacked a combined isotopic analysis, especially an in situ Lu–Hf isotope analysis of zircon. In this study, Precambrian basement rocks from the Kuluketage and Dunhuang Blocks in the northeastern portion of the Tarim Craton have been analyzed for combined in situ laser ablation ICP-(MC)-MS zircon U–Pb and Lu–Hf isotopic analyses, as well as whole rock elements, to constrain their protoliths, forming ages and magma sources. Two magmatic events from the Kuluketage Block at ∼2.4 Ga and ∼1.85 Ga are revealed, and three stages of magmatic events are detected in the Dunhuang Block, i.e., ∼2.0 Ga, ∼1.85 Ga and ∼1.75 Ga. The ∼1.85 Ga magmatic rocks from both areas were derived from an isotopically similar crustal source under the same tectonic settings, suggesting that the Kuluketage and Dunhuang Blocks are part of the uniform Precambrian basement of the Tarim Craton. Zircon Hf model ages of the ∼2.4 Ga magmatism indicate that the crust of the Tarim Craton may have been formed as early as the Paleoarchean period. The ∼2.0 Ga mafic rock from the Dunhuang Block was formed in an active continental margin setting, representing an important crustal growth event of the Tarim Craton in the mid-Paleoproterozoic that coincides with the global episode of crust formation during the assembly of the Columbia supercontinent. The ∼1.85 Ga event in the Kuluketage and Dunhuang Blocks primarily involved the reworking of the old crust and most likely related to the collisional event associated with the assembly of the Columbia supercontinent, while the ∼1.75 Ga magmatism in the Dunhuang Block resulted from a mixture of the reworked Archean crust with juvenile magmas and was most likely related to a post-collisional episode.     

Relationship of gravity anomalies and tectonics in Assam, India

Gravity data from Assam compiled on Bouguer, Hayford and Airy isostatic anomaly maps have been interpreted in terms of tectonics of the area. The gravity anomalies suggest that the Dauki fault is very deep-seated. A gravity high of about 60 mGal near Haflong is interpreted as being the expression of an intrusive body with a density contrast of about + 0.15 g/cm₃ with respect to the surroundings. From isostatic considerations, approximate crustal thicknesses over the Shillong Plateau, the Upper Assam valley and the Surma valley are estimated to be 40, 29 and 22 km respectively, suggesting a sharp change in crustal thickness from the Shillong Plateau to the Surma valley across the Dauki fault.     

Maximum magnitudes in aftershock sequences in Taiwan

In this work, Båth’s Law, the b-value in Gutenberg–Richter Law (G–R Law) in the form of the 1/β relationship, and both the a- and b-values in the G–R Law were introduced in order to estimate maximum aftershock magnitudes of earthquake sequences in the Taiwan region. The averaged difference of magnitude between the mainshock and the maximum aftershock is 1.20, and is consistent with Båth’s Law, however, with a large uncertainty. The large uncertainty implies that the difference may result from a variable controlled by other factors, such as the aftershocks number of an earthquake sequence and magnitude threshold for mainshock. With 1/β, since 86% of the earthquake sequences with a M ⩾ 6.0 mainshock follow this relationship, the upper bound of the maximum magnitude can be estimated for an earthquake sequence with a large mainshock. The a- and b-values in the G–R Law was also considered by evaluating maximum aftershock magnitudes. As there are low residuals between the model and the observations, the results suggest that the G–R Law is a good index for maximum aftershock magnitude determinations. In order to evaluate the temporal decays of maximum aftershock magnitudes, modified Omori’s Law was introduced. Using the approaches mentioned above, the maximum magnitudes and the temporal evolution of an earthquake sequence could be modeled. Among them, the model of the G–R Law has the best fit with observations for most of earthquake sequences. It shows its feasibility. The results of this work may benefit seismic hazards mitigation in the form of rapid re-evaluations for short-term seismic hazards immediately following devastating earthquakes.     

Seismic structure of the crust and uppermost mantle of South America and surrounding oceanic basins

We present a new set of contour maps of the seismic structure of South America and the surrounding ocean basins. These maps include new data, helping to constrain crustal thickness, whole-crustal average P-wave and S-wave velocity, and the seismic velocity of the uppermost mantle (P_n and S_n). We find that: (1) The weighted average thickness of the crust under South America is 38.17 km (standard deviation, s.d. ±8.7 km), which is ∼1 km thinner than the global average of 39.2 km (s.d. ±8.5 km) for continental crust. (2) Histograms of whole-crustal P-wave velocities for the South American crust are bi-modal, with the lower peak occurring for crust that appears to be missing a high-velocity (6.9–7.3 km/s) lower crustal layer. (3) The average P-wave velocity of the crystalline crust (P_cc) is 6.47 km/s (s.d. ±0.25 km/s). This is essentially identical to the global average of 6.45 km/s. (4) The average P_n velocity beneath South America is 8.00 km/s (s.d. ±0.23 km/s), slightly lower than the global average of 8.07 km/s. (5) A region across northern Chile and northeast Argentina has anomalously low P- and S-wave velocities in the crust. Geographically, this corresponds to the shallowly-subducted portion of the Nazca plate (the Pampean flat slab first described by Isacks et al., 1968), which is also a region of crustal extension. (6) The thick crust of the Brazilian craton appears to extend into Venezuela and Colombia. (7) The crust in the Amazon basin and along the western edge of the Brazilian craton may be thinned by extension. (8) The average crustal P-wave velocity under the eastern Pacific seafloor is higher than under the western Atlantic seafloor, most likely due to the thicker sediment layer on the older Atlantic seafloor.     

Crustal character and thickness over the Dangerous Grounds and beneath the Northwest Borneo Trough

New deep reflection seismic, bathymetry, gravity and magnetic data have been acquired in a marine geophysical survey of the southern South China Sea, including the Dangerous Grounds, Northwest Borneo Trough and the Central Luconia Platform. The seismic and bathymetry data map the topography of shallow density interfaces, allowing the application of gravity modeling to delineate the thickness and composition of the deeper crustal layers. Many of the strongest gravity anomalies across the area are accounted for by the basement topography mapped in the seismic data, with substantial basement relief associated with major rift development. The total crustal thickness is however quite constant, with variations only between 25 and 30 km across the Central Luconia Platform and Dangerous Grounds. The Northwest Borneo Trough is underlain by thinned crust (25–20 km total crustal thickness) consistent with the substantial water depths. There is no evidence of any crustal suture associated with the trough, nor any evidence of relict oceanic crust beneath the trough. The crustal thinning also does not extend along the complete length of the trough, with crustal thicknesses of 25 km and more modeled on the most easterly lines to cross the trough. Modeled magnetic field variations are also consistent with the study area being underlain by continental crust, with the magnetic field variations well explained by irregular magnetisations consistent with inhomogeneous continental crust, terminating at the basement unconformity as mapped from the seismic data.     

Neotectonics at the Arabian plate margins

Opening of the Red Sea is accompanied by convergence between the Arabian plate and Eurasia. Regional topography and structure favour gravity glide as the main driving force of plate translation. At the leading edge of the plate, the Zagros Mountains undergo coseismic serial folding which is equivalent to Holocene shortening by ∼20 mm/year and which has led to major episodes of coastal uplift of which the last was ∼1700 years BP. At the Jordan Rift transform, which bounds the Arabian plate on the west, a recurrence interval of ∼1600 years is reported for events of M_L≥5.5. The palaeomagnetic record for the last 3.2 Ma indicates an average spreading rate for the Red Sea of ∼20 mm/year; there is some evidence that hydrothermal activity in the Red Sea is pulsatory, with a period of ∼2000 year, and that it reflects discontinuous spreading. The Holocene neotectonic records of the Zagros, the Jordan Rift and the Red Sea are the product of complex plate interactions and of the accumulation and release of strain in the crust along the plate margins. But they also reflect elastic strain energy storage and release within the Arabian plate, whence parallels in the period of major deformation episodes in the three deforming zones and the apparent discrepancy between the seismic moment predicted by plate kinematics and that recorded in the Zagros. Any associated intraplate deformation, if detected geodetically, would thus help the assessment of seismic hazard.     

Integrated potential-field and seismic constraints on the structure of the Archean metasedimentary English River belt,Western Superior craton,Canada

The English River Subprovince is a prominent belt of metasedimentary rocks in the Archean Western Superior Province. The structure of its western half was investigated by using techniques of enhancement and automatic interpretation of magnetic data, and integration of magnetic-derived information with seismic and gravity data. The results indicate that a suite of exposed felsic plutons that intruded the belt at ca. 2698 Ma extends under most of the metasedimentary rocks that are exposed at the surface. The thickness of the metasedimentary rocks is interpreted to be less than 1 km in areas where it is underlain by the members of this intrusive suite. In other areas, the metasedimentary rocks attain thicknesses of 3–4 km and appear to be underlain by rocks similar to the gneissic rocks that are exposed in the adjacent metaplutonic Winnipeg River Subprovince. The integration of enhanced magnetic data with gravity data indicates that the large gravity anomaly that extends along the English River belt correlates well spatially and morphologically with the extensive suite of felsic intrusions that underlies the belt, suggesting that the crustal component of the gravity anomaly is related to this suite of intrusions. We interpret the source of the gravity anomaly as a dense unit comprising anhydrous mineral assemblages that formed within these felsic intrusions in response to low-pressure, high-temperature metamorphism that affected the belt at ca. 2691 Ma. On the basis of geochronological, geological and geophysical constraints, we propose that this metamorphic episode is linked to the continuation of magmatism at depth after the emplacement of the ca. 2698 Ma felsic plutons, being ultimately related to the advection of mantle heat into the crust during a period of regional extension.     

Crustal structure across southern Mexico inferred from gravity data

We present a gravity model of the crustal structure in southern Mexico based on interpretation of a detailed marine gravity profile perpendicularly across the Middle America Trench offshore from Acapulco, and a regional gravity transect extending into continental Mexico across the Sierra Madre del Sur, the central sector of the Trans-Mexican Volcanic Belt, the Sierra Madre Oriental, the Coastal Plain, and into the Gulf of Mexico. The elastic thickness of the Cocos lithospheric plate was found to be 30 km. In agreement with a previous seismic refraction study, no major differences in crustal structure were observed on both sides of the O’Gorman Fracture Zone. The gravity high seaward of the trench is interpreted as due to the incipient flexure and crustal thinning. The gravity low at the axis of the trench is explained by the increase in water depth and the existence of low-density accreted or continental-derived sediments (2.25 and 2.40 g/cm³). A gravity high of 50 mGal extending about 100 km landward is interpreted as caused by local shoaling of the Moho. The crust attains a thickness of 42 km under the Trans-Mexican Volcanic Belt but thins beneath the Coastal Plain and the continental slope of the Gulf of Mexico. Gravity highs around the Sierra de Tamaulipas are interpreted in terms of relief of the lower–upper crustal interface, implying a shallow basement.     

An attempt to detect temporal variations of crustal structure in the source area of the 2006 Wen-An earthquake in North China

Lei et al. (2008) revealed a low-velocity (low-V) anomaly in the lower crust under the source area of the 4 July 2006 Wen-An earthquake (M 5.1) in North China. In this work we tried to investigate the temporal variations of the crustal structure by using a number of P and PmP (Moho reflected) wave arrival times recorded by 107 digital seismic stations from earthquakes that occurred separately in 2002, 2003, 2004, and 2005–2006 to determine P-wave velocity structures in and around the source area of the Wen-An earthquake in different periods. Our results show that tomographic images inferred from the data sets in different years are all dominated by a low-V anomaly in the lower crust under the Wen-An source area. However, there exist some differences in the P-wave velocity image in the Wen-An source area inferred from the P and P + PmP data sets, suggesting that the PmP data have improved the tomographic images in the middle and lower crust, but the results from the 2005–2006 P and P + PmP data sets all show a relatively lager increase of the low-V anomaly under the Wen-An source area in the amplitude and extent as compared with those from the 2003 and 2004 P and P + PmP data sets. Incorporating the previous results, if this low-V anomaly may indicate the existence of fluids, then our results suggest that the occurrence of the Wen-An earthquake is not only related to the long-term influence of fluids that decrease the effective normal stress on the fault plane, but also closely associated with the drastic increase of such influence. However, this study is just an experimental work and the results are still preliminary because the resolution scale of the present tomographic model is much larger than the Wen-An source area and our extensive tests show that different samplings of seismic rays from different data sets have affected the details of the tomographic images, suggesting that the present sparse data coverage can hardly detect reliably any temporal variations of the velocity anomalies in the Wen-An source area. In future studies it is necessary to improve the resolution of crustal tomography to the size of the rupture zone and utilize identical seismic ray paths from the same pairs of sources and receivers in order to detect any temporal variations of the crust structure in the source area of a large earthquake.     

Bearing of plate geometry and rheology on shallow-focus mega-thrust seismicity with special reference to 26 December 2004 Sumatra event

We propose a model pertaining to the generation of 26th December 2004 off Sumatra mega-event in the backdrop of other similar type earthquakes along subduction zones around the world. Reconstructions of Benioff trajectories through the hypocenters of historical earthquakes including six mega-earthquakes indicate (i) confinement of hypocenters right within the descending lithosphere, and (ii) natural coincidence of foci of the mega-events around the zones of plate flexing. These observations are discussed in detail with special emphasis on the Sumatra margin considering the role of rheological anomaly across the cross-section of the descending lithosphere; yield strength envelope and residual stress accumulation through time. The intraplate origin of shallow mega-thrust earthquakes allowed us to advocate the ‘zone of flexing’ along the profiles of the subducting plates as nodal area for stress concentration. We propose here that at elevated confining pressure and temperature, loading of unidirectional cyclic stress on time-average bending stress enhanced the material yield strength (i.e., strain-hardening), and leads the semi-brittle portion of the lithosphere into near-brittle condition through rheological transformation. Under subsequent rise in neutral surface and increase in compressive stress field, non-coaxial deformation triggered shear failure on 26th December 2004 preferably at the rheological interface between strain-hardened near-brittle layer and deformed ductile layer within the sub-oceanic mantle.A two-stage fracture mechanism viz. a slow (～1.1 km/s) bilateral initiation in an essentially strain-hardened near-brittle domain and a follow-up very rapid progression (3.3 km/s) in the brittle, crustal domain was mainly involved in the generation of 2004 off Sumatra mega-event. Estimation shows an amount of 3.38 × 10²² to 4.50 × 10²² N m seismic moment (M_o) and 8.95–9.03 moment magnitude (M_w) for the southern part of the 1300 km extended rupture i.e. between the North Andaman to the north and the Sumatra at its south. The study necessitates the reassessment of other shallow-focus mega-thrust earthquakes along the subduction margins around the globe.