Relocation of aftershocks, focal mechanisms and stress inversion: Implications toward the seismo-tectonics of the causative fault zone of Mw7.6 2001 Bhuj earthquake (India)

The HYPODD relocation of 1172 aftershocks, recorded on 8–17 three-component digital seismographs, delineate a distinct south dipping E–W trending aftershock zone extending up to 35 km depth, which involves a crustal volume of 40 km × 60 km × 35 km. The relocated focal depths delineate the presence of three fault segments and variation in the brittle–ductile transition depths amongst the individual faults as the earthquake foci in the both western and eastern ends are confined up to 28 km depth whilst in the central aftershock zone they are limited up to 35 km depth. The FPFIT focal mechanism solutions of 444 aftershocks (using 8–12 first motions) suggest that the focal mechanisms ranged between pure reverse and pure strike slip except some pure dip slip solutions. Stress inversion performed using the P and T axes of the selected focal mechanisms reveals an N181°E oriented maximum principal stress with a very shallow dip (= 14°). The stress inversions of different depth bins of the P and T axes of selected aftershocks suggest a heterogeneous stress regime at 0–30 km depth range with a dominant consistent N–S orientation of the P-axes over the aftershock zone, which could be attributed to the existence of varied nature and orientation of fractures and faults as revealed by the relocated aftershocks.     

Aftershock activity of Bhuj earthquake of january 26th, 2001

Following a large-sized Bhuj earthquake (M _s = 7.6) of January 26th, 2001, a small aperture 4-station temporary local network was deployed, in the epicentral area, for a period of about three weeks and resulted in the recording of more than 1800 aftershocks (-0.07 ≤M _L <5.0). Preliminary locations of epicenters of 297 aftershocks (2.0 ≤M _L <5.0) have brought out a dense cluster of aftershock activity, the center of which falls 20 km NW of Bhachau. Epicentral locations of after-shocks encompass a surface area of about 50 × 40 km² that seems to indicate the surface projection of the rupture area associated with the earthquake. The distribution of aftershock activity above magnitude 3, shows that aftershocks are nonuniformly distributed and are aligned in the north, northwest and northeast directions. The epicenter of the mainshock falls on the southern edge of the delineated zone of aftershock activity and the maximum clustering of activity occurs in close proximity of the mainshock. Well-constrained focal depths of 122 aftershocks show that 89% of the aftershocks occurred at depths ranging between 6 and 25 km and only 7% and 4% aftershocks occur at depths less than 5 and more than 25 km respectively. The Gutenberg-Richter (GR) relationship, logN = 4.52 - 0.89M_L, is fitted to the aftershock data (1.0<-M _L<5.0) and theb-value of 0.89 has been estimated for the aftershock activity.     

Atlas of macroseismic maps for French earthquakes with their principal characteristics

The SIRENE macroseismic database has been utilized to draw isoseismal maps for the 140 best-documented French earthquakes, characterized by epicentral intensities of at least V (MSK) and located in all parts of the country. A study of focal depths derived from available local intensity data using an intensity versus distance decay law (Sponheuer) shows that the focal depths of most of the events considered do not exceed about 10 km. Their distribution correlates fairly well with regional dynamic geology features. A relationship is then computed between magnitude, intensity and focal distance, based on 73 instrumenta]ly recorded earthquakes (M _L between 3.3 and 6.3) and on 217 mean radius values (from 2 to 380 km) for isoseismals of intensity VIII to III (MSK). This relationship is applied to historical earthquakes contained in the database SIRENE which are characterised by their intensity only. These results are used in the evaluation as well deterministic as probabilistic of the seismic hazard on the national territory.     

Earthquake focal mechanisms of the induced seismicity in 2006 and 2007 below Basel (Switzerland)

To stimulate the reservoir for a “hot dry rock" geothermal project, that was initiated by a private/public consortium in the city of Basel, approximately 11500 m³ of water were injected between December 2nd and 8th, 2006, at high pressures into a 5 km deep well. More than 10500 seismic events were recorded during the injection phase, and minor sporadic seismic activity was still occurring more than two years later. The present article documents the focal mechanisms of the 28 strongest events, with M_L between 1.7 and 3.4, that have been obtained by the Swiss Seismological Service (SED) during and after the stimulation. The analysis is based on data that was recorded by a six-station borehole network, operated by the project developers, as well as by several permanent and temporary surface networks. The hypocenters of the events are located inside the stimulated rock volume at depths between 4 and 5 km within the crystalline basement. Of the 28 faultplane solutions two are normal faulting mechanisms and one is a strike-slip mechanism with a strong normal component. All others are typical strike-slip mechanisms with mostly NS and EW striking nodal planes. As a consequence, the T-axes are all nearly horizontal and oriented in a NE or SW direction (mean azimuth 46 ± 11 degrees) and the P-axes of the strike-slip events point in a NW or SE direction (mean azimuth 138 ± 13 degrees). Overall, the observed focal mechanisms agree with what would be expected from both the stress observations within the well and the stress field derived from the previously known natural seismicity.     

A study of seismic reflection imaging using microearthquake sources

We present a method for processing three-component digital recordings of microearthquakes to obtain near-vertical reflection profiles in regions of shallow seismicity. The processing includes magnitude and focal depth normalization and event stacking, where stacking is by small localized groups, with ray theoretical time and distance corrections applied to compensate for varying focal depths. In areas with high seismicity, this procedure allows earthquakes to be treated as “controlled” sources to probe layered structures of the deep crust and upper mantle.The validity of our approach is examined using aftershocks of the Borah Peak, Idaho earthquake (M_s = 7.3). Several thousand events occurred in a NNW-trending zone about 10 km wide, 75 km long, and 15 km deep. A small (~ 10 ×10 km) array of nine University of Wisconsin three-component triggered short-period digital seismographs was installed in the region of aftershock activity. Over a 10-day period, about 1000 useable events were recorded, of which about 120 have been used for this study. Hypocenters have been computed using both P- and S-wave arrivals, the latter being essential for stable solutions of events outside the network.The Borah Peak data have been processed to obtain shear-wave reflection profiles for the central station (Station 8) of the digital station array. The stacked shear wave (transverse) record sections reveal coherent reflections from horizons at mid-crustal to Moho depths. The most prominent reflections are from crustal discontinuities in the depth range 18–28 km. Coherent reflections can be obtained only through stacking, which is necessary to improve the signal to noise ratio. The major sources of data scatter, as manifested by “smearing” of reflections on the stacked records, are crustal heterogeneity and errors in the determination of focal depth and origin time.     

Physicochemical parameters of crystallization of melts in intermediate suprasubduction chambers (by the example of Tolbachik and Ichinskii Volcanoes,Kamchatka Peninsula)

In study of plagioclases, amphiboles, and melt inclusions, we have determined the physicochemical parameters of crystallization of melts in the intermediate suprasubduction chambers of volcanoes representing different types of subduction magmatism on the Kamchatka Peninsula: the young basaltic systems of Tolbachik Volcano (Klyuchevskaya group) and ancient Ichinskii Volcano (Sredinnyi Ridge) with alternating basaltic and felsic eruptions. For Tolbachik Volcano, we have found that plagioclase lapilli formed from basaltic melts at 1075-1115 °C and low (< 1 kbar) pressures at depths of 2-3 km. Andesite minerals crystallized within a wider range of temperatures and pressures (1220-1020 °C and 3.3-1.6 kbar) in an intermediate chamber at depths of < 10 km. The melts were generated in basaltic magma chambers (detected well by geophysical methods at depths of 18-20 km) with minimum temperatures of ~ 1290 °C. For Ichinskii Volcano, three levels of intermediate chambers are distinguished. Andesites formed at depths of < 23 km at < 1225 °C. Dacitic melts were generated from an intermediate chamber (14 km) at 1135-1045 °C as a result of differentiation of andesitic magmas. Dacites formed in the uppermost horizons (9-3 km) at 1130-1030 °C. Despite the similarity between differentiation processes in the intermediate chambers of the Kamchatka volcanoes, each volcano is characterized by specific magmatism. The lavas of basaltic volcanoes (Tolbachik) and those of andesitic volcanoes (Ichinskii) differ in genesis and differentiation.     

Seismicity and source parameters of local earthquakes in Bilaspur region of Himachal Lesser Himalaya

The pattern of local seismicity (110 events) and the source parameters of 26 local events (1.0?≤?Mw?≤?2.5) that occurred during May 2008 to April 2009 in Bilaspur region of Himachal Lesser Himalaya were determined. The digital records available from one station have been used to compute the source parameters and f _max based on the Brune source model (1970) and a high-frequency diminution factor (Boore 1983) above f _max. The epicentral distribution of events within 30 km of local network is broadly divided into three clusters of seismic activity: (1) a cluster located to the south of the Jamthal (JAMT) station and falls to the north of the Main Boundary Thrust (MBT) which seems to reflect the contemporary local seismicity of the segment of the MBT, (2) an elongated zone of local seismicity NE–SW trending, delineated NE of JAMT station that falls in the Lesser Himalaya between the MBT and the Main Central Thrust, and (3) NE–SW trending zone of local seismic activity located at about 10 km east of NHRI station and about 15 km northeast of NERI station and extending over a distance of about 20 km. Majority of events occur at shallow depths up to 20 km, and the maximum number of events occurs in the focal depth range between 10 and 15 km. The entire seismic activity is confined to the crust between 5 and 45 km. The average values of these source parameters range from 3.29?×?10¹⁷ to 3.73?×?10¹⁹?dyne-cm for seismic moment, 0.1 to 9.7 bars for stress drops, and 111.78 to 558.92 m for source radii. The average value of f _max for these events varies from 7 to 18 Hz and seems to be source dependent.     

Aftershocks of the june 20, 1978, Greece earthquake: A multimode faulting sequence

A 10-station portable seismograph network was deployed in northern Greece to study aftershocks of the magnitude (m_b) 6.4 earthquake of June 20, 1978. The main shock occurred (in a graben) about 25 km northeast of the city of Thessaloniki and caused an east-west zone of surface rupturing 14 km long that splayed to 7 km wide at the west end. The hypocenters for 116 aftershocks in the magnitude range from 2.5 to 4.5 were determined. The epicenters for these events cover an area 30 km (east-west) by 18 km (north-south), and focal depths ranges from 4 to 12 km. Most of the aftershocks in the east half of the aftershock zone are north of the surface rupture and north of the graben. Those in the west half are located within the boundaries of the graben. Composite focalmechanism solutions for selected aftershocks indicate reactivation of geologically mapped normal faults in the area. Also, strike-slip and dip-slip faults that splay off the western end of the zone of surface ruptures may have been activated.The epicenters for four large (M 4.8) foreshocks and the main shock were relocated using the method of joint epicenter determination. Collectively, those five epicenters form an arcuate pattern convex southward, that is north of and 5 km distant from the surface rupturing. The 5-km separation, along with a focal depth of 8 km (average aftershock depth) or 16 km (NEIS main-shock depth), implies that the fault plane dips northward 58° or 73°, respectively. A preferred nodal-plane dip of 36° was determined by B.C. Papazachos and his colleagues in 1979 from a focal-mechanism solution for the main shock. If this dip is valid for the causal fault and that fault projects to the zone of surface rupturing, a decrease of dip with depth is required.     

Conductivity Structure along the Gediz Graben,West Anatolia,Turkey: Tectonic Implications

Groups of grabens in west Anatolia have contrasting E-W and NE-SW orientations and are the subject of debate as to their relative ages and relationships. We investigated the E-W-trending Gediz graben and its neighboring NE-SW-trending Gördes, Demirci, and Selendi grabens, which form an important graben system representative of the region. We studied gravity data from one profile and magnetotelluric (MT) data from two profiles, 73 km and 93 km long. The data supports the hypothesis that the Gediz graben was superimposed onto the (older) NE-SW grabens. 2D gravity and MT modelling revealed an undulating graben floor, varying in depth between 500 and 3000-4000 m (gravity-MT); within the graben two apparent basins 3–4 and 1.5-2.5 km deep (gravity-MT) are separated by a subsurface horst. The residual gravity map appears to indicate the continuation of NE-SW grabens from north of Gediz graben to beyond its southern border.

The MT model revealed three main zones of varying thickness within the crust. The britde upper crust comprises two zones: sedimentary fill (apparent resistivity 15-50 ohm.m) and Menderes massif basement (200 ohm.m). The third zone is highly conductive lower crust (10 ohm.m), identified by our MT modeling at an average depth of 10 km. This conductive layer was considered in conjunction with two other regional features, high heat flow values and shallow earthquake focal depths. A heat flow map shows a very high average value of 108 mWm^?2 for west Anatolia and 120-300 mWm^?2 for the Gediz graben area specifically, compared with the world average of 80 mWm^?2. Seismological records showing shallow earthquake focal depths together with the high conductivity zone were taken to indicate a partially melted, viscoelastic lower crust.     

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.     

The October 20, 2011 Mw 5.1 Talala earthquake in the stable continental region of India

The Talala (Sasangir) area in the Saurashtra region of Gujarat, western India, is experiencing tremors since 2001. The swarm type of earthquake activity in 2001, 2004, and every year from 2007 onward has occurred after the monsoon and lasted 2?C3?months each time. In 2007 some 200 shocks (largest M_w 5.0) and in 2011 about 400 shocks down to M1 are well recorded with 1?C2?km location error. The focal depths are about 2?C10?km and shocks are accompanied by blast-like subterranean sounds. The epicenter (21.09?N 70.45E, focal depth: 5?km from location program, 3?km from MTS) of the October 20, 2011 mainshock occurred about 12-km WNW of Talala town or 8-km SSW of the 2007?M _w 5.0 earthquake epicenter. The epicentral trends deciphered from local earthquake data indicate two ENE trends (Narmada trend) for about 50?km length and a conjugate 15-km-long NNW trend (Aravali trend). The focal mechanisms by moment-tensor analysis of full wave forms of two 2007 events of M_w 4.8 and 5.0 and the 2011 event of M_w 5.1 indicate rupture along either of the two trends. The ENE trends follow a gravity low between the gravity highs of Girnar mounts. Seismic reflections also indicate a fault in the area named Girnar Fault. Most of Saurashtra region including the Talala area is covered by Deccan Trap Basalt forming plateaus and conical ridges. There is no major fault within Saurashtra Peninsula though it is believed to have major faults along the boundaries that are non-seismic. The intensity of the October 20, 2011 Talala earthquake is estimated to be 6.5 in MM scale while isoseismals of 6, 5, and 4 and felt distance give M_w 5.1 based on Johnston??s 1994 empirical regressions. The source parameters of the 2011 Talala earthquake are estimated using data from 14 broadband seismograph stations. Estimated seismic moment, moment magnitude, stress drop, corner frequency, and source radius are found to be 10^16.6 N-m, 5.1, 1.6?MPa, 1.3?Hz, and 2,300?m, respectively. The b and p values are obtained to be low, being 0.67 and 0.71, respectively. PGA of 35?cm/sec² is noted and the decay rate of acceleration has been estimated from strong motion data recorded at 5 stations with epicentral distances ranging from 32 to 200?km.     

Fluid and enthalpy production during regional metamorphism

Models for regional metamorphism have been constructed to determine the thermal effects of reaction enthalpy and the amount of fluid generated by dehydration metamorphism. The model continental crust contains an average of 2.9 wt % water and dehydrates by a series of reactions between temperatures of 300 and 750° C. Large scale metamorphism is induced by instantaneous collision belt thickening events which double the crustal thickness to 70 km. After a 20 Ma time lag, erosion due to isostatic rebound restores the crust to its original thickness in 100 Ma. At crustal depths greater than 10 km, where most metamorphism takes place, fluid pressure is unlikely to deviate significantly from lithostatic pressure. This implies that lower crustal porosity can only be maintained if rock pores are filled by fluid. Therefore, porosities are primarily dependent on the rate of metamorphic fluid production or consumption and the crustal permeability. In the models, permeability is taken as a function of porosity; this permits estimation of both fluid fluxes and porosities during metamorphism. Metamorphic activity, as measured by net reaction enthalpy, can be categorized as endothermic or exothermic depending on whether prograde dehydration or retrograde hydration reactions predominate. The endothermic stage begins almost immediately after thickening, peaks at about 20 Ma, and ends after 40 to 55 Ma. During this period the maximum and average heat consumption by reactions are on the order 11.2·10^–14 W/cm³ and 5.9·10^–14 W/ cm³, respectively. The maximum rates of prograde isograd advance decrease from 2.4·10^–8 cm/s, for low grade reactions at 7 Ma, to 7·10^–10 cm/s, for the highest grade reaction between 45 and 58 Ma. Endothermic cooling reduces the temperature variation in the metamorphic models by less than 7% (40 K); in comparison, the retrograde exothermic heating effect is negligible. Dehydration reactions are generally poor thermal buffers, but under certain conditions reactions may control temperature over depth and time intervals on the order of 1 km and 3 Ma. The model metamorphic events reduce the hydrate water content of the crust to values between 1.0 and 0.4 wt % and produce anhydrous lower crustal granulites up to 15 km in thickness. In the first 60 Ma of metamorphism, steady state fluid fluxes in the rocks overlying prograde reaction fronts are on the order of 5·10^–11 g/cm²-s. These fluid fluxes can be accommodated by low porosities (<0.6%) and are thus essentially determined by the rate of devolitalization. The quantity of fluid which passes through the metamorphic column varies from 25000 g/cm², within 10 km of the base of the crust, to amounts as large as 240000 g/cm², in rocks initially at a depth of 30 km. Measured petrologic volumetric fluid-rock ratios generated by this fluid could be as high as 500 in a 1 m thick horizontal layer, but would decrease in inverse proportion of the thickness of the rock layer. Fluid advection causes local heating at rates of about 5.9·10^–14 W/cm³ during prograde metamorphism and does not result in significant heating. The amount of silica which can be transported by the fluids is very sensitive to both the absolute temperature and the change in the geothermal gradient with depth. However, even under optimal conditions, the amount of silica precipitated by metamorphic fluids is small (<0.1 vol %) and inadequate to explain the quartz veining observed in nature. These results are based on equilibrium models for fluid and heat transport that exclude the possibility of convective fluid recirculation. Such a model is likely to apply at depths greater than 10 km; therefore, it is concluded that large scale heat and silica transport by fluids is not extensive in the lower crust, despite large time-integrated fluid fluxes.     

Seismogenesis of clustered seismicity beneath the Kangra–Chamba sector of northwest Himalaya: Constraints from 3D local earthquake tomography

To investigate subsurface structure and seismogenic layers, 3D velocity inversion was carried out in the source zone of 1905 Kangra earthquake (M8.0) in the northwestern Himalaya. P-wave and S-wave phase data of 159 earthquakes recorded by a network of 21 stations were used for this purpose. Inverted velocity tomograms up to a depth range of 18 km show significant variations of 14% in Vp and Vs and 6% in the Vp/Vs across the major tectonic zones in the region. Synthesis of seismicity pattern, velocity structure, distinctive focal mechanisms coupled with nature of stress distribution allows mapping of three different source regions that control regional seismotectonics. Accumulating strains are partly consumed by sliding of Chamba Nappe to the southwest through reverse-fault movements along Chamba/Panjal/Main Boundary Thrusts. This coupled with normal-fault type displacements along Chenab Normal Fault in the north account for low magnitude widespread seismicity in upper 8–10 km of the crust. At intermediate depths from 8 to 15 km, adjusting to residual compressive stresses, the detachment or lower end of the MBT slips to produce thrust dominated seismicity. Nucleation of secondary stresses in local NE–SW oriented structure interacts in complex manner with regional stresses to generate normal type earthquakes below the plane of detachment and therefore three seismic regimes at different depths produce intense seismicity in a block of 30 × 30 km² centered NE to the epicenter of Kangra earthquake.     

On the origin of saline fluids in the KTB (Continental Deep Drilling Project of Germany)

Highly saline fluids were encountered during the German Continental Deep Drilling Project (KTB) from depths ranging between 2 and 3 km to about 9 km. The most reliable data were obtained from samples extracted during a long-term pumping test in the 4000-m deep KTB pilot hole. Some 460 m³ Ca–Na–Cl brines with about 68 g l⁻¹ total dissolved solids (TDS) and some 270 m³ associated gases, mainly N₂ and CH₄ were pumped to the surface from the main fracture system situated near the bottom of the pilot hole. Geochemical and isotopic data support the hydraulic tests which suggest the presence of an open and large fluid reservoir at depth. The pumped fluids from this main fracture system were released from a deep reservoir situated at more than 5500 m depth which is hydraulically connected with the 9101 m deep KTB main hole, drilled some 250 m to the northeast of the pilot hole.While Ca and Sr contents of the extracted brines may be the result of water–rock interaction, Cl is most likely of external origin. The Cl is hypothesized to derive from geotectonic processes rather than to descending infiltration of paleo-seawater (evaporitic brines). The sampled fluids have probably migrated from a deeper reservoir to their present position since the Cretaceous–Tertiary period due to tectonic activity. However, several isotopic studies have identified an admixture of descending paleowaters down to more than 4000 m depth. The high ³⁶Cl/Cl ratio of the fluids sampled during the long-term pumping test point to a host rock highly enriched in U–Th, unlike the sampled KTB country rocks. The fluid reservoir is believed to be in contact with the Falkenberg granite massif situated about 2 km to the E of the KTB holes, capable of supplying sufficient neutron flux for considerable subsurface production of ³⁶Cl. The Na–Cl–(K-, SO₄) precursor fluids of the Ca–Na–Cl brines were produced in the course of extensive tectonic processes since the Late Caledonian within the Bohemian Massif.     

Physical properties of ultrahigh-pressure metamorphic rocks from the Sulu terrain, eastern central China: implications for the seismic structure at the Donghai (CCSD) drilling site

About 30 samples representing major lithologies of Sulu ultrahigh-pressure (UHP) metamorphic rocks were collected from surface exposures and exploration wells, and compressional (Vp) and shear wave (Vs) velocities and their directional dependence (anisotropy) were determined over a range of constant confining pressures up to 600 MPa and temperatures ranging from 20 to 600 °C. Samples range in composition from acidic to ultramafic. P- and S-wave velocities measured at 600 MPa vary from 5.08 to 8.64 km/s and 2.34 to 4.93 km/s, respectively. Densities are in the range from 2.60 to 3.68 g/cm³. To make a direct tie between seismic measurements (refraction and reflection) and subsurface lithologies, the experimental velocity data (corresponding to shallow depths) were used to calculate velocity profiles for the different lithologies and profiles of reflection coefficients at possible lithologic interfaces across the projected 5000-m Chinese Continental Scientific Drilling Program (CCSD) crustal segment. Comparison of calculated in situ velocities with respective intrinsic velocities suggests that the in situ velocities at shallow depths are lowered by an increased abundance of open microcracks. The strongly reflective zone beneath the Donghai drill site can be explained by the impedance contrasts between the different lithologies. Contacts between eclogite/peridotite and felsic rocks (gt-gneiss, granitic gneiss), in particular, may give rise to strong seismic reflections. In addition, shear-induced (lattice preferred orientation (LPO)-related) seismic anisotropy can increase reflectivity. For the explanation of the high velocity bodies (>6.4 km/s) around 1000 m and below 3200-m depth, large proportions of eclogite/peridotite (about 40 and 30 vol.%, respectively) are needed.     

Observations of the electrical asthenosphere beneath Scandinavia

The recent recognition that long period (i.e., of the order of hours) electromagnetic induction studies could play a major role in the detection of the asthenosphere has led to much interest amongst the geophysical and geological communities of the geomagnetic response functions derived for differing tectonic environments. Experiments carried out on the ocean bottom have met with considerable success in delineating the “electrical asthenosphere”, i.e., a local maximum in electrical conductivity (minimum in electrical resistivity) in the upper mantle.In this paper, observations of the time-varying magnetic field recorded in three regions of Scandinavia, northern Sweden (Kiruna—KIR), northern Finland/northeastern Norway (Kevo—KEV) and southern Finland (Sauvamaki—SAU), are analysed in order to obtain estimates of the inductive response function, C(ω), for each region. The estimated response functions are compared with one from the centre of the East European Platform (EEP), and it is shown that the induced eddy currents, at periods of the order of 10³–10⁴ s, in the three regions flow much closer to the surface than under the platform centre. Specifically, at a period of ~3000 s, these currents are flowing at depths of the order of: KEV—120 km; KIR—180 km; SAU—210 km; EEP—280 km; implying that the transition to a conducting zone, of σ -0.2 S/m, occurs at around these depths. One-dimensional inversion of and shows that there must exist a good conducting zone, of σ = 0.1–1.0 S/m, under each of the two regions, of 40 km minimum thickness, at depths of: KEV 105–115 km; KIR 160–185 km. This is to be contrasted with EEP, where the ρ-d profile displays a monotonically decreasing resistivity with depth, reaching σ~0.1 S/m at > 300 km.Finally, a possible temperature range for the asthenosphere, consistent with the deduced conducvitity, is discussed. It is shown that, at present, there is insufficient knowledge of the conditions (water content, melt fraction, etc.) likely to prevail in the asthenosphere to narrow down the probable range of 900°–1500°C.     

Tomographic image of crustal structures across the Chelungpu fault: Is the seismogenic layer structure- or depth-dependent?

In 1999, a large earthquake (M_w = 7.6) occurred along the Chelungpu fault in the fold-and-thrust belt of western Taiwan. To shed more light on the subsurface structures and the seismogenic layers, three-dimensional velocity structures were inverted by using the travel times of both P- and S-waves from 2391 aftershocks recorded by the Central Weather Bureau during the 15 months that followed. From tomography, a typical image of the large-scale thrusting structures in the upper crust across the Chelungpu fault was obtained. In general, high velocities beneath the Western Foothills and Central Ranges are separated from low velocities beneath the Coastal Plain by an east-dipping boundary that is roughly consistent with the Chelungpu fault on the surface. The contrast in velocity on either side of the Chelungpu fault is indicative of about a 7- to 9-km vertical offset in the upper crust. The relocated hypocenter for the Chi-Chi earthquake shifts by 2.2 km toward the northwest, and its focal depth decreases by 0.7 km. A plot of focal depths versus rock velocities where the aftershocks occurred shows earthquakes are more inclined to occur in rock with a velocity of around 5.6 km/s. This strongly suggests the seismogenic layer in the fold-and-thrust belt of Taiwan is more structure-dependent than depth-dependent.     

A new tectonic model for Abu-Dabbab seismogenic zone (Eastern Desert,Egypt): evidence from field-structural,EMR and seismic data

Abu-Dabbab area is the most active seismic zone in the central Eastern Desert of Egypt, where seismic activities are daily recorded. The reported earthquakes are microearthquakes of local magnitudes (ML < 2.0). A spatial distribution of these microearthquakes shows that the earthquakes of the area follow an ENE–WSW trending pattern, which is nearly perpendicular to the Red Sea Rift. Focal mechanisms of different fault styles were recognized with dominant normal faulting (with a strike-slip component) events characterized by focal depths greater than 7 km and reverse ones of shallower focal depths. Several lines of evidence indicating that the brittle-ductile transition zone underlies the Abu-Dabbab area occurs at a relatively shallow depth (10–12 km) and it is acting as a low-angle normal shear zone (LANF). Field-structural, EMR and seismic data (this study) reveal that the maximum compressive stress (σ1) in the area is perturbed from the regional NW–SE direction to ENE–WSW orientation. This stress rotation is evidently akin to the reactivation of the crustal scale Najd Fault System (NFS), where such reactivation is attributed to the ongoing activity/opening of the Red Sea. Our tectonic model proposes that the continuous activity on the brittle-ductile transition zone including the LANF led to stress localization, which triggering a brittle deformation in the upper crustal-levels and associated shallow dipping thrusts. Such bimodal tectonic model suggests that the deep earthquakes are owing to the tectonic movement on the LANF (transtension), whereas the shallow earthquakes are related to a brittle deformation inside the fault blocks of the upper crust (transpression). Deformation creep along this zone didn’t permit continuous accumulation of strain and hence reduce the possible occurrence of large earthquakes.     

Monitoring massive fracture growth at 2-km depths using surface tiltmeter arrays

Tilt due to massive hydraulic fractures induced in sedimentary rocks at depths of up to 2.2 km have been recorded by surface tiltmeters. Injection of fluid volumes up to 4 · 10⁵ liters and masses of propping agent up to 5 · 10⁵ kg is designed to produce fractures approximately 1 km long, 50–100 m high and about 1 cm wide. The surface tilt data adequately fit a dislocation model of a tensional fault in a half-space. Theoretical and observational results indicate that maximum tilt occurs at a distance off the strike of the fracture equivalent to 0.4 of the depth to the fracture. Azimuth and extent of the fracture deduced from the geometry of the tilt field agree with other kinds of geophysical measurements. Detailed correlation of the tilt signatures with pumping parameters (pressure, rate, volume, mass) have provided details on asymmetry in geometry and growth rate. Whereas amplitude variations in tilt vary inversely with the square of the depth, changes in flow rate or pressure gradient can produce a cubic change in width. These studies offer a large-scale experimental approach to the study of problems involving fracturing, mass transport, and dilatancy processes.     

Mantle-derived mafic-ultramafic xenoliths and the nature of Indian sub-continental lithosphere

Mantle derived xenoliths in India are known to occur in the Proterozoic ultrapotassic rocks like kimberlites from Dharwar and Bastar craton and Mesozoic alkali igneous rocks like lamrophyres, nephelinites and basanites. The xenoliths in kimberlites are represented by garnet harzburgites, lherzolites, wehrlite, olivine clinopyroxenites and kyaniteeclogite varieties. The PT conditions estimated for xenoliths from the Dharwar craton suggest that the lithosphere was at least 185 km thick during the Mid-Proterozoic period. The ultrabasic and eclogite xenoliths have been derived from depths of 100–180 km and 75–150 km respectively. The Kalyandurg and Brahmanpalle clusters have sampled the typical Archaean subcontinental lithospheric mantle (SCLM) with a low geotherm (35 mW/m²) and harzburgitic to lherzolitic rocks with median X_mg ^olivine > 0.93. The base of the depleted lithosphere at 185–195 km depth is marked by a 10–15 km layer of strongly metasomatised peridotites (X_mg ^olivine > ∼0.88). The Anampalle and Wajrakarur clusters 60 km to the NW show a distinctly different SCLM; it has a higher geotherm (37.5 to 40 mW/m²) and contains few subcalcic harzburgites, and has a median X_mg ^olivine = 0.925. In contrast, the kimberlites of the Uravakonda and WK-7 clusters sampled quite fertile (median X_mg ^olivine ∼0.915) SCLM with an elevated geotherm (> 40 mW/m²). The lamrophyres, basanites and melanephelinites associated with the Deccan Volcanic Province entrain both ultramafic and mafic xenoliths. The ultramafic group is represented by (i) spinel lherzolites, harzburgites, and (ii) pyroxenites. Single pyroxene granulite and two pyroxene granulites constitutes the mafic group. Temperature estimates for the West Coast xenoliths indicate equilibration temperatures of 500–900°C while the pressure estimates vary between 6–11 kbar corresponding to depths of 20–35 km. This elevated geotherm implies that the region is characterized by abnormally high heat flow, which is also supported by the presence of linear array of hot springs along the West Coast. Spinel peridotite xenoliths entrained in the basanites and melanephelinites from the Kutch show low equilibrium temperatures (884–972°C). The estimated pressures obtained on the basis of the absence of both plagioclase and garnet in the xenoliths and by referring the temperatures to the West Coast geotherm is ∼ 15 kbar (40–45 km depth). The minimum heat flow of 60 to 70 mW/m² has been computed for the Kutch xenolith (Bhujia hill), which is closely comparable to the oceanic geotherm. Xenolith studies from the West Coast and Kutch indicate that the SCLM beneath is strongly metasomatised although the style of metasomatism is different from that below the Dharwar Craton.