An application of regional time and magnitude predictable model for long-term earthquake prediction in the vicinity of October 8, 2005 Kashmir Himalaya earthquake

A regional time and magnitude predictable model has been applied to estimate the recurrence intervals for large earthquakes in the vicinity of 8 October 2005 Kashmir Himalaya earthquake (25°–40°N and 65°–85°E), which includes India, Pakistan, Afghanistan, Hindukush, Pamirs, Mangolia and Tien-Shan. This region has been divided into 17 seismogenic sources on the basis of certain seismotectonics and geomorphological criteria. A complete earthquake catalogue (historical and instrumental) of magnitude Ms ≥ 5.5 during the period 1853–2005 has been used in the analysis. According to this model, the magnitude of preceding earthquake governs the time of occurrence and magnitude of future mainshock in the sequence. The interevent time between successive mainshocks with magnitude equal to or greater than a minimum magnitude threshold were considered and used for long-term earthquake prediction in each of seismogenic sources. The interevent times and magnitudes of mainshocks have been used to determine the following predictive relations: logT _t = 0.05 M _min + 0.09 M _p − 0.01 log M ₀ + 01.14; and M _f = 0.21 M _min − 0.01 M _p + 0.03 log M ₀ + 7.21 where, T _t is the interevent time of successive mainshocks, M _min is minimum magnitude threshold considered, M _p is magnitude of preceding mainshock, M _f is magnitude of following mainshock and M ₀ is the seismic moment released per year in each seismogenic source. It was found that the magnitude of following mainshock (M _f) does not depend on the interevent time (T _t), which indicates the ability to predict the time of occurrence of future mainshock. A negative correlation between magnitude of following mainshock (M _f) and preceding mainshock (M _p) indicates that the larger earthquake is followed by smaller one and vice versa. The above equations have been used for the seismic hazard assessment in the considered region. Based on the model applicability in the studied region and taking into account the occurrence time and magnitude of last mainshock in each seismogenic source, the time-dependent conditional probabilities (PC) for the occurrence of next shallow large mainshocks (Ms ≥ 6.5), during next 20 years as well as the expected magnitudes have been estimated.     

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.     

Magnetic anisotropy of hematite natural crystals: high field experiments

Hematite is one of the most important carriers of remanent magnetization in natural samples. Its strong magnetocrystalline anisotropy makes it difficult to determine one single value for the magnetic anisotropy constant and other magnetic properties. In particular, the anisotropy of hematite within its basal plane is controversial because an assumed triaxial anisotropy compatible with the crystallographic structure has not always been detected. This study presents a comparative analysis of rock magnetic properties, compositional analysis and determination of the magnetic anisotropy constant. Different models with anisotropy constant within the basal plane ranging from 0 to 13 (J m⁻³) are considered in the evaluation of the factors that control the presence or absence of a triaxial anisotropy. A linear relationship between saturation magnetization (M _s ) and coercitivity (B _c ) in measurements at increasing temperatures is observed in samples where anisotropy is either uniaxial or biaxial while those with significant triaxial anisotropies have a power–law relationship between M _s and B _c .     

On the validity of time-predictable model for earthquake generation in north-east India

North-east India is seismically very active and has experienced many widelydistributed shallow, large earthquakes. Earthquake generation model for the region was studied using seismicity data [(1906–1984) prepared by National Geophysical Data Centre (NGDC), Boulder Colorado, USA]. For establishing statistical relations surface wave magnitudes (M _s≥5·5) have been considered. In the region four seismogenic sources have been identified which show the occurrences of atleast three earthquakes of magnitude 5·5≤M _s≤7·5 giving two repeat times. It is observed that the time interval between the two consecutive main shock depends on the preceding main shock magnitude (M _p) and not on the following main shock magnitude (M _f) revealing the validity of time predictable model for the region. Linear relation between logarithm of repeat time (T) and preceding main shock magnitude (M _p) is established in the form of logT=cM _p+a. The values ofc anda are estimated to be 0–36 and 1–23, respectively. The relation may be used for seismic hazard evaluation in the region.     

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.     

Influence of different strain rates on hydro-mechanical behaviour of reconstituted unsaturated soil

The hydro-mechanical behaviour of a reconstituted unsaturated soil under different suctions and strain rates was studied through various rate-controlled unsaturated/undrained triaxial tests. The fully saturated reconstituted specimens were desaturated to four different initial suctions (s₀?=?0, 100 kPa, 200 kPa and 300 kPa) and then triaxially sheared (conventional triaxial compression) at three different strain rates in undrained conditions (\(\dot{\varepsilon }_{1} = 0.001\) h^?1, 0.01 h^?1, and 0.1 h^?1). The observed hydro-mechanical behaviour during shearing including the volumetric strain, deviatoric stress, degree of saturation and suction is presented and discussed in this paper. The results indicate that when the strain rate rises at the given initial suctions (or pore water pressures), the maximum deviatoric stress (q_max), critical net stress ratio (M) and critical state suction (s_c) increase but the degree of saturation (S_rc) and volumetric strain at the critical state (ε^c_v ) reduce. The critical effective stress ratio (M′) is not dependent on the strain rate for saturated and unsaturated samples. The critical state lines for unsaturated soils with the constant strain rates are parallel with each other in the e???lnp′ space.

     

The tectono‐metamorphic evolution of the very low‐grade hangingwall constrains two‐stage gneiss dome formation in the Montagne Noire (Southern France)

The Montagne Noire in the southernmost French Massif Central is made of an ENE‐elongated gneiss dome flanked by Palaeozoic sedimentary rocks. The tectonic evolution of the gneiss dome has generated controversy for more than half a century. As a result, a multitude of models have been proposed that invoke various tectonic regimes and exhumation mechanisms. Most of these models are based on data from the gneiss dome itself. Here, new constraints on the dome evolution are provided based on a combination of very low‐grade petrology, K–Ar geochronology, field mapping and structural analysis of the Palaeozoic western Mont Peyroux and Faugères units, which constitute part of the southern hangingwall of the dome. It is shown that southward‐directed Variscan nappe‐thrusting (D₁) and a related medium‐P metamorphism (M₁) are only preserved in the area furthest away from the gneiss dome. The regionally dominant pervasive tectono‐metamorphic event D₂/M₂ largely transposes D₁ structures, comprises a higher metamorphic thermal gradient than M₁ (transition low‐P and medium‐P metamorphic facies series) and affected the rocks between c. 309 and 300 Ma, post‐dating D₁/M₁ by more than 20 Ma. D₂‐related fabrics are refolded by D₃, which in its turn, is followed by dextral‐normal shearing along the basal shear zone of both units at c. 297 Ma. In the western Mont Peyroux and Faugères units, D₂/M₂ is largely synchronous with shearing along the southern dome margin between c. 311 and 303 Ma, facilitating the emplacement of the gneiss dome into the upper crust. D₂/M₂ also overlaps in time with granitic magmatism and migmatization in the Zone Axiale between c. 314 and 306 Ma, and a related low‐P/high‐T metamorphism at c. 308 Ma. The shearing that accompanied the exhumation of the dome therefore was synchronous with a peak in temperature expressed by migmatization and intrusion of melts within the dome, and also with the peak of metamorphism in the hangingwall. Both, the intensity of D₂ fabrics and the M₂ metamorphic grade within the hangingwall, decrease away from the gneiss dome, with grades ranging from the anchizone–epizone boundary to the diagenetic zone. The related zonation of the pre‐D₃ metamorphic field gradients paralleled the dome. These observations indicate that D₂/M₂ is controlled by the exhumation of the Zone Axiale, and suggest a coherent kinematic between the different crustal levels at some time during D₂/M₂. Based on integration of these findings with regional geological constraints, a two‐stage exhumation of the gneiss dome is proposed: during a first stage between c. 316 and 300 Ma dome emplacement into the upper crust was controlled by dextral shear zones arranged in a pull‐apart‐like geometry. The second stage from 300 Ma onwards was characterized by northeast to northward extension, with exhumation accommodated by north‐dipping detachments and hangingwall basin formation along the northeastern dome margin.     

Global regression relations for conversion of surface wave and body wave magnitudes to moment magnitude

A homogenous earthquake catalog is a basic input for seismic hazard estimation, and other seismicity studies. The preparation of a homogenous earthquake catalog for a seismic region needs regressed relations for conversion of different magnitudes types, e.g. m _b , M _s , to the unified moment magnitude M _w. In case of small data sets for any seismic region, it is not possible to have reliable region specific conversion relations and alternatively appropriate global regression relations for the required magnitude ranges and focal depths can be utilized. In this study, we collected global events magnitude data from ISC, NEIC and GCMT databases for the period 1976 to May, 2007. Data for m_b magnitudes for 3,48,423 events for ISC and 2,38,525 events for NEIC, M _s magnitudes for 81,974 events from ISC and 16,019 events for NEIC along with 27,229 M _w events data from GCMT has been considered. An epicentral plot for M _w events considered in this study is also shown. M _s determinations by ISC and NEIC, have been verified to be equivalent. Orthogonal Standard Regression (OSR) relations have been obtained between M _s and M _w for focal depths (h < 70 km) in the magnitude ranges 3.0 ≤ M _s  ≤ 6.1 and 6.2 ≤ M _s  ≤ 8.4, and for focal depths 70 km ≤ h ≤ 643 km in the magnitude range 3.3 ≤ M _s  ≤ 7.2. Standard and Inverted Standard Regression plots are also shown along with OSR to ascertain the validation of orthogonal regression for M _s magnitudes. The OSR relations have smaller uncertainty compared to SR and ISR relations for M _s conversions. ISR relations between m _b and M _w have been obtained for magnitude ranges 2.9 ≤ m _b  ≤ 6.5, for ISC events and 3.8 ≤ m _b  ≤ 6.5 for NEIC events. The regression relations derived in this study based on global data are useful empirical relations to develop homogenous earthquake catalogs in the absence of regional regression relations, as the events catalog for most seismic regions are heterogeneous in magnitude types.     

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

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

A curious case of agreement between conventional thermobarometry and phase equilibria modelling in granulites: New constraints on P–T estimates in the Antarctica segment of the Musgrave–Albany–Fraser–Wilkes Orogen

The Windmill Islands region in Wilkes Land, east Antarctica, preserves granulite facies metamorphic mineral assemblages that yield seemingly comparable P–T estimates from conventional thermobarometry and mineral equilibria modelling. This is uncommon in granulite facies terranes, where conventional thermobarometry and phase equilibria modelling generally produce conflicting P–T estimates because peak mineral compositions tend to be modified by retrograde diffusion processes. In situ U–Pb monazite geochronology and calculated metamorphic phase diagrams show that the Windmill Islands experienced two phases of high thermal gradient metamorphism during the Mesoproterozoic. The first phase of metamorphism is recorded by monazite ages in two widely separated samples and occurred at c. 1,305 Ma. This event was regional in extent, involved crustally derived magmatism and reached conditions of ~3.2–5 kbar and 690–770°C corresponding to very high thermal gradients of >150°C/kbar. The elevated thermal regime is interpreted to reflect a period of extension or increased extension in a back‐arc setting that existed prior to c. 1,330 Ma. The first metamorphic event was overprinted by granulite facies metamorphism at c. 1,180 Ma that was coeval with the intrusion of charnockite. This event involved peak temperatures of ~840–850°C and pressures of ~4–5 kbar. A phase of granitic magmatism at c. 1,250–1,210 Ma, prior to the intrusion of the charnockite, is interpreted to reflect a phase of compression within an overall back‐arc setting. Existing conventional thermobarometry suggests conditions of ~4 kbar and 750°C for M₁ and 4–7 kbar and 750–900°C for M₂. The apparent similarities between the phase equilibria modelling and existing conventional thermobarometry may suggest either that the terrane cooled relatively quickly, or that the P–T ranges obtained from conventional thermobarometry are sufficiently imprecise that they cover the range of P–T conditions obtained in this study. However, without phase equilibria modelling, the veracity of existing conventional P–T estimates cannot be evaluated. The calculated phase diagrams from this study allow the direct comparison of P–T conditions in the Windmill Islands with phase equilibria models from other regions in the Musgrave–Albany–Fraser–Wilkes Orogen. This shows that the metamorphic evolution of the Wilkes Land region is very similar to that of the eastern Albany–Fraser Orogen and Musgrave Province in Australia, and further demonstrates the remarkable consistency in the timing of metamorphism and the thermal gradients along the ~5,000 km strike length of this system.     

P–T evolution and episodic zircon growth in barroisite eclogites of the Lanterman Range,northern Victoria Land,Antarctica

Prograde P–T–t paths of eclogites are often ambiguous owing to high variance of mineral assemblages, large uncertainty in isotopic age determinations and/or variable degree of retrograde equilibration. We investigated these issues using the barroisite eclogites from the Lanterman Range, northern Victoria Land, Antarctica, which are relatively uncommon but free of retrogression. These eclogites revealed three stages of prograde metamorphism, defining two distinctive P–T trajectories, M_1–2 and M₃. Inclusion minerals in garnet porphyroblasts suggest that initial prograde assemblages (M₁) consist of garnet+omphacite+barroisite/Mg‐pargasite+epidote+phengite+paragonite+rutile/titanite+quartz, and subsequent M₂ assemblages of garnet+omphacite+barroisite+phengite+rutile±quartz. The inclusion‐rich inner part of garnet porphyroblasts preserves a bell‐shaped Mn profile of the M₁, whereas the inclusion‐poor outer part (M₂) is typified by the outward decrease in Ca/Mg and X_Fe (=Fe²⁺/(Fe²⁺+Mg)) values. A pseudosection modelling employing fractionated bulk‐rock composition suggests that the eclogites have initially evolved from ~15 to 20 kbar and 520–570°C (M₁) to ~22–25 kbar and 630–650°C (M₂). The latter is in accordance with P–T conditions estimated from two independent geothermobarometers: the garnet–clinopyroxene–phengite (~25 ± 3 kbar and 660 ± 100°C) and Zr‐in‐rutile (~650–700°C at 22–27 kbar). The second segment (M_3A–B) of prograde P–T path is recorded in the grossular‐rich overgrowth rim of garnet. Apart from disequilibrium growth of the M_3A garnet, ubiquitous overgrowth of the M_3B garnet permits us to estimate the P–T conditions at ~26 ± 3 kbar and 720 ± 80°C. The cathodoluminescence (CL) imaging of zircon grains separated from a barroisite eclogite revealed three distinct zones with bright rim, dark mantle and moderately dark core. Eclogitic phases such as garnet, omphacite, epidote and rutile are present as fine‐grained inclusions in the mantle and rim of zircon, in contrast to their absence in the core. The sensitive high‐resolution ion microprobe U–Pb dating on metamorphic mantle domains and neoblasts yielded a weighted mean ²⁰⁶Pb/²³⁸U age of 515 ± 4 Ma (tσ), representing the time of the M₂ stage. On the other hand, overgrowth rims as well as bright‐CL neoblasts of zircon were dated at 498 ± 11 Ma (tσ), corresponding to the M₃. Average burial rates estimated from the M₂ and M₃ ages are too low (<2 mm/year) for cold subduction regime (~5–10°C/km), suggesting that an exhumation stage intervened between two prograde segments of P–T path. Thus, the P–T–t evolution of barroisite eclogites is typified by two discrete episodes with an c. 15 Ma gap during the middle Cambrian subduction of the Antarctic Ross Orogeny.     

OH in zoned amphiboles of eclogite from the western Tianshan,NW-China

Chemically-zoned amphibole porphyroblast grains in an eclogite (sample ws24-7) from the western Tianshan (NW-China) have been analyzed by electron microprobe (EMP), micro Fourier-transform infrared (micro-FTIR) and micro-Raman spectroscopy in the OH-stretching region. The EMP data reveal zoned amphibole compositions clustering around two predominant compositions: a glaucophane end-member ( ^B Na₂ ^C M²⁺ ₃ M³⁺ ₂ ^T Si₈(OH)₂) in the cores, whereas the mantle to rim of the samples has an intermediate amphibole composition ( ^A _0.5 ^B Ca_1.5Na_0.5 ^C M ²⁺ _4.5 M _0.5^{3+ T} Si_7.5Al_0.5(OH)₂) (A = Na and/or K; M ²⁺ = Mg and Fe²⁺; M ³⁺ = Fe³⁺ and/or Al) between winchite (and ferro-winchite) and katophorite (and Mg-katophorite). Furthermore, we observed complicated FTIR and Raman spectra with OH-stretching absorption bands varying systematically from core to rim. The FTIR/Raman spectra of the core amphibole show three lower-frequency components (at 3,633, 3,649–3,651 and 3,660–3,663 cm⁻¹) which can be attributed to a local O(3)-H dipole surrounded by ^{M(1) M(3)}Mg₃, ^{M(1) M(3)}Mg₂Fe²⁺ and ^{M(1) M(3)} Fe²⁺ ₃, respectively, an empty A site and ^T Si₈ environments. On the other hand, bands at higher frequencies (3,672–3,673, 3,691–3,697 and 3,708 cm⁻¹) are observable in the rims of the amphiboles, and they indicate the presence of an occupied A site. The FTIR and Raman data from the OH-stretching region allow us to calculate the site occupancy of the A, M(1)–M(3), T sites with confidence when combined with EPM data. By contrast M(2)- and M(4) site occupancies are more difficult to evaluate. We use these samples to highlight on the opportunities and limitations of FTIR OH-stretching spectroscopy applied to natural high pressure amphibole phases. The much more detailed cation site occupancy of the zoned amphibole from the western Tianshan have been obtained by comparing data from micro-chemical and FTIR and/or Raman in the OH-stretching data. We find the following characteristic substitutions Si(T-site) (Mg, Fe)[M(1)–M(3)-site] → Al(T-site) Al[M(1)–M(3)-site] (tschermakite), Ca(M4-site)□ (A-site) → Na(M4-site) Na + K(A-site) (richterite), and Ca(M4-site) (Mg, Fe) [M(1)–M(3)-site] → Na(M4-site) Al[M(1)–M(3)-site] (glaucophane) from the configurations observed during metamorphism.     

Long-term prediction of intermediate depth earthquakes in the southern Aegean region based on a time-predictable model

Repeat times of strong intermediate depth (60 km h 180 km) earthquakes have been determined by the use of instrumental and historical data for six seismogenic sources in the Benioff zone of the southern Aegean area. For four of these sources, at least two interevent times (three mainshocks) are available for each source. By using the repeat times for these four sources, the following relation has been determined: logT _t = 0.20M _min + 0.19M _p +a, whereT _t is the repeat time (in years),M _min the surface wave magnitude of the smallest earthquake considered,M _p the magnitude of the preceding mainshock and a parameter which varies from source to source. A multilinear correlation coefficient equal to 0.91 was determined for this relation, which indicates that the time predictable model holds to a satisfactory degree for the strong mainshocks of intermediate focal depth in the southern Aegean.By assuming that the ratioT/T _t, whereT is the observed andT _t the calculated repeat time, follows a lognormal distribution, the conditional probabilities for the occurrence of strong (M _s 6.5) and very strong (M _s 7.5) earthquakes during the period 1991–2001 in these four seismogenic sources have been calculated. These probabilities are very high (P > 0.9) for the strong and high (P > 0.5) for the very strong intermediate depth earthquakes which occur in the three sources of the shallower (h < 100 km) part of the Benioff zone where coupling occurs between the front parts of the Mediterranean lithosphere (downgoing) and the Aegean lithosphere.     

3-D seismic structure of the source area of the 1993 Latur, India, earthquake and its implications for rupture nucleations

The Latur earthquake (M_w 6.1) of 29 September 1993 is a rare stable continental region (SCR) earthquake that occurred on a previously unknown blind fault. In this study, we determined detailed three-dimensional (3-D) P- and S-wave velocity (V_p, V_s) and Poisson's ratio (σ) structures by inverting the first P- and S-wave high-quality arrival time data from 142 aftershocks that were recorded by a network of temporary seismic stations. The source zone of the Latur earthquake shows strong lateral heterogeneities in V_p, V_s and σ structures, extending in a volume of about 90 × 90 × 15 km³. The mainshock occurred within, but near the boundary, of a low-V_p, high-V_s and low-σ zone. This suggests that the structural asperities at the mainshock hypocenter are associated with a partially fluid-saturated fractured rock in a previously unknown source zone with intersecting fault surfaces. This might have triggered the 1993 Latur mainshock and its aftershock sequence. Our results are in good agreement with other geophysical studies that suggest high conductivity and high concentration of radiogenic helium gas beneath the source zone of the Latur earthquake. Our study provides an additional evidence for the presence of fluid related anomaly at the hidden source zone of the Latur earthquake in the SCR and helps us understand the genesis of damaging earthquakes in the SCR of the world.     

Continental rifting as a function of lithosphere mantle strength

The role of the uppermost mantle strength in the pattern of lithosphere rifting is investigated using a thermo-mechanical finite-element code. In the lithosphere, the mantle/crust strength ratio (S_M/S_C) that decreases with increasing Moho temperature T_M allows two strength regimes to be defined: mantle dominated (S_M > S_C) and crust dominated (S_M < S_C). The transition between the two regimes corresponds to the disappearance of a high strength uppermost mantle for T_M > 700 °C. 2D numerical simulations for different values of S_M/S_C show how the uppermost mantle strength controls the style of continental rifting. A high strength mantle leads to strain localisation at lithosphere scale, with two main patterns of narrow rifting: “coupled crust–mantle” at the lowest T_M values and “deep crustal décollement” for increasing T_M values, typical of some continental rifts and non-volcanic passive margins. The absence of a high strength mantle leads to distributed deformations and wide rifting in the upper crust. These numerical results are compared and discussed in relation with series of classical rift examples.     

Thermal expansion and high-temperature P21/c–C2/c phase transition in clinopyroxene-type LiFeGe2O6 and comparison to NaFe(Si,Ge)2O6

A synthetic clinopyroxene with composition LiFe³⁺Ge₂O₆, monoclinic s.g. P2₁/c, a = 9.8792(7), b = 8.8095(5), c = 5.3754(3) Å, β = 108.844(6)°, V = 442.75(16) ų, has been studied by in situ low- and high-temperature single-crystal X-ray diffraction. The variation of lattice parameters and the intensity of the b-type reflections (h + k = 2n + 1, only present in the P-symmetry) with increasing temperature showed a displacive phase transition from space group P2₁/c to C2/c at a transition temperature T _tr = 789 K, first order in character, with a sudden volume increase of 1.6% and a decrease of β by 1° at the transition. This spontaneous dilatation is reversible, shows a limited hysteresis of ±10°C, and corresponds to the vanishing of the b-type reflections, thus indicating a symmetry increase to space group C2/c. Below T _tr an expansion is observed for all the cell parameters, while the β angle remained almost constant; at T > T _tr the thermal volume expansion is due to dilatation of the structure in the $(\bar{1}\,0\,1) A synthetic clinopyroxene with composition LiFe³⁺Ge₂O₆, monoclinic s.g. P2₁/c, a = 9.8792(7), b = 8.8095(5), c = 5.3754(3) ?, β = 108.844(6)°, V = 442.75(16) ?³, has been studied by in situ low- and high-temperature single-crystal X-ray diffraction. The variation of lattice parameters and the intensity of the b-type reflections (h + k = 2n + 1, only present in the P-symmetry) with increasing temperature showed a displacive phase transition from space group P2₁/c to C2/c at a transition temperature T _tr = 789 K, first order in character, with a sudden volume increase of 1.6% and a decrease of β by 1° at the transition. This spontaneous dilatation is reversible, shows a limited hysteresis of ±10°C, and corresponds to the vanishing of the b-type reflections, thus indicating a symmetry increase to space group C2/c. Below T _tr an expansion is observed for all the cell parameters, while the β angle remained almost constant; at T > T _tr the thermal volume expansion is due to dilatation of the structure in the ([`1] 0 1)(\bar{1}\,0\,1) plane, mostly along [0 1 0], and pure shear in the (0 1 0) plane due to the decrease of β. From comparison with silicate analogues, the germanate clinopyroxenes are more expansible, while the P2₁/c expands more than the C2/c phase. The evolution of Q ² (calculated as the normalized intensity of b-type reflections) with T in the framework of the Landau theory has been done using a standard expression for a first order phase transition. We observe a jump of Q ₀² = 0.538(2) at T _tr, with T _c of 481(7) K, b/a = −2,290 K, and c/a = 3,192 K, and thus far from being tri-critical point. A closely related composition (LiFe³⁺Si₂O₆) shows an equivalent phase transition at 228 K, which is very close to the tri-critical point and 561 K cooler. This result indicates that a change in the composition of tetrahedral sites can have dramatic effects on the P2₁/c ↔ C2/c displacive phase transition in clinopyroxenes. The major changes observed in the evolution of the crystal structure with T are observed in the M2 polyhedron, with a volume decrease by ca. 13.3%, compared to ca. 1.3% observed in the M1 polyhedron. The tetrahedra behave as rigid units with neither a significant change of volume at T > T _tr (<1‰), nor a change of tilting of the basal plane. No change in coordination is observed at T > T _tr in the M2 polyhedron, which remains sixfold coordinated although a strong deformation of this polyhedron is observed. This deformation is related to a strong change by 51.4° at T _tr of the kinking angle (O3–O3–O3 angle) of the B-chain of tetrahedra, which switches from O-rotated to S-rotated [from 143.3(5)° to 194.7(6)°]. The A-chain is S-rotated at T < T _tr [206.8(5)° at 703 K] and extends by 12° at the transition.     

Relative stress variations as determined by b-values from earthquakes in circum-pacific subduction zones

Assuming a relation of “b” to stress state, the possibility of globe-wide stress variation and transmission was investigated. The NOAA earthquake data file served to determine the temporal change in “b” of log N = a − bM from 1963 to 1975.Periods of six to eight years are observed in the b-values (stress pattern) for most circum-Pacific areas (South America, Tonga, Kermadec, New Hebrides, Kamchatka and Eastern Aleutians).In the Kurils, fore- and aftershock sequences of large earthquakes seem to mask any characteristic global pattern that might exist. These sequences exhibit low b-values (high stress) through the time of foreshocks and early stages of aftershocks, followed by rapid increase in b-values (decrease in stress).Use of a worldwide earthquake data file clearly yields less resolution of the temporal “b” variation than the use of local network studies published by other authors.Incidental to the study, 1124 earthquakes of the NOAA data file yield the M_s − m_b relations: M_s = 1.16m_b − 0.835 for 4.5 m_b 6 and: log₁₀M_s = 0.1432m_b − 0.0629 form_b > 6 with correlative coefficients of 0.994 and 0.992 respectively.     

On the Ratios between Elastic Modulus and Uniaxial Compressive Strength of Heterogeneous Carbonate Rocks

The ratios M _R = E/σ _c for 11 heterogeneous carbonate (dolomites, limestones and chalks) rock formations collected from different regions of Israel were examined. Sixty-eight uniaxial compressive tests were conducted on weak-to-strong (5 MPa < σ _c < 100 MPa) and very strong (σ _c > 100 MPa) rock samples exhibiting wide ranges of elastic modulus (E = 6100–82300 MPa), uniaxial compressive strength (σ _c = 14–273.9 MPa), Poisson's ratio (ν = 0.13–0.49), and dry bulk density (ρ = 1.7–2.7 g/cm³). The observed range of M _R = 60.9–1011.4 and mean value of M _R = 380.5 are compared with the results obtained by Deere (Rock mechanics in engineering practice, Wiley, London, pp 1–20, 1968) for limestones and dolomites, and the statistical analysis of M _R distribution is performed. Mutual relations between E, σ _c, ρ, M _R for all studied rocks, and separately for concrete rock formations are revealed. Linear multiple correlations between E on the one hand and σ _c and ρ on the other for Nekorot and Bina limestone and Aminadav dolomite are obtained. It is established that the elastic modulus and M _R in very strong carbonate samples are more correlated with ρ−σ _c combination and ε _a max, respectively, than in weak to strong samples. The relation between M _R and maximum axial strain (ε _a max) for all studied rock samples (weak-to-strong and very strong) is discussed.     

Thermoelasticity and high-<Emphasis Type="Italic">T</Emphasis> behaviour of anthophyllite

The thermoelastic behaviour of anthophyllite has been determined for a natural crystal with crystal-chemical formula ^ANa_0.01 ^B(Mg_1.30Mn_0.57Ca_0.09Na_0.04) ^C(Mg_4.95Fe_0.02Al_0.03) ^T(Si_8.00)O₂₂ ^W(OH)₂ using single-crystal X-ray diffraction to 973 K. The best model for fitting the thermal expansion data is that of Berman (J Petrol 29:445–522, 1988) in which the coefficient of volume thermal expansion varies linearly with T as α _V,T  = a ₁ + 2a ₂ (T − T ₀): α₂₉₈ = a ₁ = 3.40(6) × 10⁻⁵ K⁻¹, a ₂ = 5.1(1.0) × 10⁻⁹ K⁻². The corresponding axial thermal expansion coefficients for this linear model are: α _a ,₂₉₈ = 1.21(2) × 10⁻⁵ K⁻¹, a _2,a  = 5.2(4) × 10⁻⁹ K⁻²; α _b ,₂₉₈ = 9.2(1) × 10⁻⁶ K⁻¹, a _2,b  = 7(2) × 10⁻¹⁰ K⁻². α _c ,₂₉₈ = 1.26(3) × 10⁻⁵ K⁻¹, a _2,c  = 1.3(6) × 10⁻⁹ K⁻². The thermoelastic behaviour of anthophyllite differs from that of most monoclinic (C2/m) amphiboles: (a) the ε ₁ − ε ₂ plane of the unit-strain ellipsoid, which is normal to b in anthophyllite but usually at a high angle to c in monoclinic amphiboles; (b) the strain components are ε ₁ ≫ ε ₂ > ε ₃ in anthophyllite, but ε ₁ ~ ε ₂ ≫ ε ₃ in monoclinic amphiboles. The strain behaviour of anthophyllite is similar to that of synthetic C2/m ^ANa ^B(LiMg) ^CMg₅ ^TSi₈ O₂₂ ^W(OH)₂, suggesting that high contents of small cations at the B-site may be primarily responsible for the much higher thermal expansion ⊥(100). Refined values for site-scattering at M4 decrease from 31.64 epfu at 298 K to 30.81 epfu at 973 K, which couples with similar increases of those of M1 and M2 sites. These changes in site scattering are interpreted in terms of Mn ↔ Mg exchange involving M1,2 ↔ M4, which was first detected at 673 K.     

Petrology and thermobarometry of mafic granulites and migmatites from the Chafalote Metamorphic Suite: New insights into the Neoproterozoic P–T evolution of the Uruguayan—Sul-Rio-Grandense shield

This paper reports a study of the metamorphic evolution of pelitic, semi-pelitic migmatites and mafic granulites of the Chafalote Metamorphic Suite (CMS), Uruguay, which represents the southernmost exposures of high-grade metamorphic rocks in the Dom Feliciano Belt, Uruguain—Sul-Rio-Grandense shield, South America. This belt is one of the Brasiliano orogens that crop out along the Brazilian and Uruguayan Atlantic margin, and the CMS is one of several disconnected segments of supracrustal rock in a dominantly granitic terrain. Petrological evidence from CMS mafic granulites and semi-pelitic migmatites indicates four distinct metamorphic assemblages. The early prograde assemblage (M₁) is preserved only as inclusions in porphyroblasts of the peak-metamorphic (M₂) assemblage. Peak-metamorphism was followed by near-isothermal decompression (M₃), which resulted in symplectites and coronitic textures in the mafic granulites and compositional zoning of Ca in garnet (decreasing rimwards) and plagioclase (increasing rimwards) in the semi-pelitic migmatites. The retrograde metamorphic assemblage (M₄) is represented by hydration reaction textures replacing minerals of the M₂ and M₃ assemblages. Average P–T calculations using the program THERMOCALC and conventional thermobarometric methods yield peak-metamorphic (M₂) P–T conditions of 7–10 kbar and 830–950 °C, near-decompressional (M₃) P–T conditions of 4.8–5.5 kbar and 788–830 °C and M₄ retrograde P–T conditions of 3–6 kbar and 600–750 °C. The calculated P–T path for the CMS rocks is ‘clockwise’ and incorporates a near-isothermal decompression segment followed by minor cooling, consistent with a history of crustal thickening followed by extensional collapse at ca. 650–600 Ma. The metamorphism recorded by rocks of this crustal segment may be correlated with 650 Ma metamorphism in the Coastal Terrane of the Kaoko Belt in Namibia, being the first unequivocal match between South America and Africa provided by crystalline rocks south of the Congo Craton.