Seismotectonics and large earthquake generation in the Himalayan region

Bounded by the western and eastern syntaxes, the Himalayan region has experienced at least five M ~ 8 earthquakes during a seismically very active phase from 1897 through 1952. However, there has been a paucity of M ~ 8 earthquakes since 1952. Examining of various catalogues and seismograms from the Gottingen Observatory, it is established that this quiescence of M ~ 8 earthquakes is real. While it has not been possible to forecast earthquakes, there has been a success in making a medium term forecast of an M 7.3 earthquake in the adjoining Indo-Burmese arc. Similarly we find that in the central Himalayan region, earthquakes of M > 6.5 have been preceded by seismic swarms and quiescences. In the recent past, based on GPS data, estimates have been made of the accumulated strains and it is postulated that a number of M ~ 8 earthquakes are imminent in the Himalayan region. We examine these estimates and find that while earthquakes of M ~ 8 may occur in the region, however, the available GPS data and their interpretation do not necessarily suggest their size and time of occurrence and whether an earthquake in a particular segment will occur sooner in comparison to that in the neighboring segment. We also comment on the inference of occurrence of M ~ 8 earthquakes based on M8 algorithm for the region. We conclude that while an M ~ 8 earthquake could occur any time anywhere in the Himalayan region, there is no indication as of now as to where and when it would occur. We impress on the need for preparedness to mitigate the pending earthquake disaster in the region.     

Location of the pilot borehole for investigations of reservoir triggered seismicity at Koyna,India

Artificial water reservoir triggered earthquakes are now known to have occurred at over 120 sites globally. The part played by the reservoirs in triggering is not exactly known due to lack of near field observations of triggered earthquakes. Koyna, located near the west coast of India, where triggered earthquakes have been occurring since 1962 provides an excellent site for near field observations of the target M ≥ 2 earthquakes. A 6 borehole seismic network has been deployed recently in the Koyna region at depths of 981–1522 m to improve the hypocenter locations. During May–December 2015, a total of 1039 earthquakes of M_L ≥ 0.5 were located using the borehole seismic network. The region is also monitored through a dense network of 23 surface broad-band stations. Our analysis indicates a significant improvement in the estimation of absolute locations of earthquakes with errors of the order of ± 300 m, combining both the networks. Based on seismicity, and logistics, a block of 2 × 2 km² area has been chosen for drilling the first pilot borehole of ~ 3 km depth, where M ≥ 2 earthquakes have been occurring frequently since 2005.     

Effect of local soil conditions on dynamic ground response in the southern coast of Izmir Bay,Turkey

The city of Izmir, located at the western end of Turkey, has experienced many strong earthquakes throughout its history. The southern coast of Izmir Bay, one of the most densely populated areas of Izmir, is located on deep alluvial sediments. It is important to determine the effect of local soil conditions on dynamic ground response in the study area, where thick loose water-saturated alluvial sediments exist. A database including geotechnical and geological information on the study area is constructed. Majority of the site is classified as D and E according to NEHRP provisions. Dynamic site response analyses are performed with EERA by utilizing the field and laboratory test results and earthquake time histories of moderate-scale earthquakes such as 1977 Izmir (M_L = 5.3), 2003 Urla (M_d = 5.6), and 2005 Uzunkuyu-Urla (M_L = 5.9), which occurred in and nearby Izmir. In addition, a scenario ground motion generated by the Izmir Fault with a magnitude of 6.5, having an average distance of 10 km to the study area, is also considered. The output data obtained from the dynamic site response analyses are evaluated, and maps displaying variation in dynamic parameters on ground surface are prepared for the southern coast of Izmir Bay, Turkey. Consequently, the dynamic analyses performed with the soil models constituted for the study area verified the damage occurred in a close distance event of 1977 Izmir earthquake. The scenario earthquake resulted in peak ground accelerations more than 0.6 g at the eastern and western ends of the study area. However, long distance events resulted in spectral amplifications by up to 5 times. With this study, it is emphasized that local soil conditions should be evaluated individually in the area of interest. Generation of a site-specific design spectrum is recommended for the areas located on deep alluvial sediments.     

Melting phase relations of the Udachnaya-East Group-I kimberlite at 3.0–6.5 GPa: Experimental evidence for alkali-carbonatite composition of primary kimberlite melts and implications for mantle plumes

Experiments on the origin of the Udachnaya-East kimberlite (UEK) have been performed using a Kawai-type multianvil apparatus at 3–6.5 GPa and 900–1500 °C. The studied composition represents exceptionally fresh Group-I kimberlite containing (wt.%): SiO₂ = 25.9, TiO₂ = 1.8, Al₂O₃ = 2.8, FeO = 9.0, MgO = 30.1, CaO = 12.7, Na₂O = 3.4, K₂O = 1.3, P₂O₅ = 1.0, Cl = 0.9, CO₂ = 9.9, and H₂O = 0.5. The super-solidus assemblage consists of melt, olivine (Ol), Ca-rich (26.0–30.2 wt.% CaO) garnet (Gt), Al-spinel (Sp), perovskite (Pv), a CaCO₃ phase (calcite or aragonite), and apatite. The low pressure assemblage (3–4 GPa) also includes clinopyroxene. The apparent solidus was established between 900 and 1000 °C at 6.5 GPa. At 6.5 GPa and 900 °C Na–Ca carbonate with molar ratio of (Na + K)/Ca ≈ 0.44 was observed. The UEK did not achieve complete melting even at 1500 °C and 6.5 GPa, due to excess xenogenic Ol in the starting material. In the studied P–T range, the melt has a Ca-carbonatite composition (Ca# = molar Ca/(Ca + Mg) ratio = 0.62–0.84) with high alkali and Cl contents (7.3–11.4 wt.% Na₂O, 2.8–6.7 wt.% K₂O, 1.6–3.4 wt.% Cl). The K, Na and Cl contents and Ca# decrease with temperature. It is argued that the primary kimberlite melt at depths > 200 km was an essentially carbonatitic (< 5 wt.% SiO₂), but evolved toward a carbonate–silicate composition (up to 15–20 wt.% SiO₂) during ascent. The absence of orthopyroxene among the run products indicates that xenogenic orthopyroxene was preferentially dissolved into the kimberlite melt. The obtained subliquidus phase assemblage (Ol + Sp + Pv + Ca-rich Gt) at P–T conditions of the UEK source region, i.e. where melt was in the last equilibrium with source rock before magma ascent, differs from the Opx-bearing peridotitic mineral assemblage of the UEK source region. This difference can be ascribed to the loss of substantial amounts of CO₂ from the kimberlite magma at shallow depths, as indicated by both petrological and experimental data. Our study implies that alkali-carbonatite melt would be a liquid phase within mantle plumes generated at the core–mantle boundary or shallower levels of the mantle, enhancing the ascent velocity of the plumes. We conclude that the long-term activity of a rising hot mantle plume and associated carbonatite melt (i.e. kimberlite melt) causes thermo-mechanical erosion of the subcontinental lithosphere mantle (SCLM) roots and creates hot and deformed metasomatic regions in the lower parts of the SCLM, which corresponds to depths constrained by P–T estimates of sheared Gt-peridotite xenoliths. The sheared Gt-peridotites undoubtedly represent samples of these regions.     

Geologic hazards associated with seismogenic faulting in southern Siberia and Mongolia: forms and location patterns

The forms and location patterns of geologic hazards induced by earthquakes in southern Siberia, Mongolia, and northern Kazakhstan in1950 through 2008 have been investigated statistically, using a database of coseismic effects created as a GIS MapInfo application, with a handy input box for large data arrays. The database includes 689 cases of macroseismic effects from M_S = 4.1–8.1 events at 398 sites. Statistical analysis of the data has revealed regional relationships between the magnitude of an earthquake and the maximum distance of its environmental effects (soil liquefaction and subsidence, secondary surface rupturing, and slope instability) to the epicenter and to the causative fault. Thus estimated limit distances to the fault for the M_S = 8.1 largest event are 40 km for soil subsidence (sinkholes), 80 km for surface rupture, 100 km for slope instability (landslides etc.), and 130 km for soil liquefaction. These distances are 3.5–5.6 times as short as those to the epicenter, which are 150, 450, 350, and 450 km, respectively. Analysis of geohazard locations relative to nearest faults in southern East Siberia shows the distances to be within 2 km for sinkholes (60% within 1.5 km), 4.5 km for landslides (90% within 1.5 km), 8 km for liquefaction (69% within 1 km), and 35.5 km for surface rupture (86% within 2 km). The frequency of hazardous effects decreases exponentially away from both seismogenic and nearest faults. Cases of soil liquefaction and subsidence are analyzed in more detail in relation to rupture patterns. Equations have been suggested to relate the maximum sizes of secondary structures (sinkholes, dikes, etc.) with the earthquake magnitude and shaking intensity at the site. As a result, a predictive model has been created for locations of geohazard associated with reactivation of seismogenic faults, assuming an arbitrary fault pattern. The obtained results make basis for modeling the distribution of geohazards for the purposes of prediction and estimation of earthquake parameters from secondary deformation.     

Shear wave attenuation characteristics over the Central India Tectonic Zone and its surroundings

The study area lies between latitude 18–26°N and longitude 73–83°E, and mainly covers the Central India Tectonic Zone (CITZ). The frequency-dependent shear wave quality factor (Q_s) has been estimated over the CITZ and its surroundings using Double Spectral Ratio (DSR) method. We have considered 25 local earthquakes with magnitude (M_L) varies from 3.0 to 4.7 recorded at 11 stations running under national seismic network. The Fast Fourier Transformed (FFT) spectra were computed from the recorded waveform having time-window from onset of S-phase to 1.0 s and for a frequency-band of 0.1–10 Hz. Three different shear wave velocities (i.e., 3.87, 3.39 and 3.96 km/s) were obtained over the study area based on a pair of earthquakes recorded at a pair of stations. The low Q_s values of 51–96 at 1 Hz (i.e., Q_s = 51f^0.49; Q_s = 90f^0.488 and Q_s = 96f^0.53) were found in the area covering the Son–Narmada–Tapti (SONATA) lineament, CITZ, eastern part of the Satpura fold belt, Vindhyan and Gondwana basins, Godavari and Mahanadi grabens, and southern part of Gangetic plain. Intermediate Q_s values of the order of 204–277 (i.e., Q_s = 204f^0.56 and Q_s = 277f^0.55) were noted in the cartonic areas, namely, Bundelkhand, Dharwar-Bhandara and Bastar. While the higher Q_s values of 391–628 at 1 Hz (i.e., Q_s = 391f^0.49, Q_s = 409f^0.48, Q_s = 417f^0.48, Q_s = 500f^0.66, Q_s = 585f^0.65 and Q_s = 628f^0.69) were found in the eastern part of the SONATA, CITZ, and the northeastern part of the Satpura fold belt. The low Q_s values might be attributing to the more heterogeneous SONATA rift system. Low Q_s values further may presumably be associated with lower-level of seismicity and apparently account for higher tectonic stress accumulation over long duration. The long-term accumulated stress is generally released through occasional triggering of moderate magnitude earthquakes in the SONATA zone. Surrounding the SONATA region, the higher Q_s values possibly accounts for a more homogeneous subsurface structure along the SONATA zone.     

Influence of the 2008 Wenchuan earthquake (Mw 7.9) on the occurrence probability of future earthquakes on neighboring faults

The Longquan–Shan fault and the Huya fault are two major neighboring faults of the Longmen–Shan fault zone where the 12 May 2008 Wenchuan earthquake (M_w 7.9) occurred. To study the influence of the Wenchuan event on these two active faults, we calculate changes of Coulomb stress on the Longquan–Shan fault and the Huya fault caused by the Wenchuan mainshock. Our results indicate that the Coulomb stress in the northern section (Zone A) of the Longquan–Shan fault is increased by 0.07–0.10 bars, that in the middle section (Zone B) by 0.04–0.11 bars, and that in the southern section (Zone C) shows almost no change. For the Huya fault, the Coulomb stress is decreased by 0.01–0.03 bars in the northern section (Zone A), 0.10–0.35 bars in the middle section (Zone B), and nearly 0.5 bars in the southern section (Zone C). The epicenter distribution of small earthquakes (M_L ⩾ 1.5) on the Longquan–Shan fault and the Huya fault after the Wenchuan earthquake is consistent with the distribution of the Coulomb stress change. This implies that the Wenchuan earthquake may have triggered small events on the Longquan–Shan fault, but inhibited those on the Huya fault. We then use the rate/state friction law to calculate the occurrence probability of future earthquakes in the study region for the next decade. They include the distribution of b-values, magnitude of completeness (M_c), the background seismicity rate, a value of Aσ_n and the duration for the transient effect (t_a) in the study region. We also estimate the earthquake occurrence probabilities on the neighboring faults after the Wenchuan earthquake. Our results show that, the occurrence probability of future earthquakes in the Longquan–Shan has a slight increase, being 7% for M ⩾ 5.0 shocks during the next decade, but the earthquake probability in the Huya region is reduced obviously, being 5–20%, 7–26% and 3–9% for M ⩾ 5.0 shocks during the next decade in sections A, B and C of the Huya fault, respectively.     

Evidence for Holocene activity of the Yilan-Yitong fault,northeastern section of the Tan-Lu fault zone in Northeast China

The Tan-Lu fault zone (TLFZ) is the largest of the major faults in eastern China. Many strong earthquakes have occurred on its section in North China, but no quake greater than M ⩾ 6 has been documented in history at its northeastern section, the Yilan-Yitong fault (YYF) in Northeast China. It is usually considered that this fault has been inactive since late Quaternary and incapable of generating moderate-sized quakes. This conclusion is, however, questioned by our recent work based on high-resolution satellite image interpretation and field investigation. We found a 70-km-long surface scarp near Fangzheng county in Heilongjiang province (HLJP) and a 20-km-long scarp near Shulan county in Jilin province (JLP), and both are associated with the YYF. The trenches across these two scarps reveal a ¹⁴C displacement date of 1730 ± 40 years BP at Fangzheng and of 4410 ± 30 years BP at Shulan. The dextral offsets of the Songhua River and Second Songhua River and nearly horizontal fault striations indicate that the new activity of the YYF has been dominated by dextral strike slipping with a normal component. These new data suggest that, at least for partial sections, the YYF has been active since the Holocene, implying a potential seismic hazard. However, current quake-protection standards in this region are very low due to the previous view that the YYF fault has not been active since the late Quaternary. If an M ⩾ 7 quake takes place on this fault, it will be a devastating event. Therefore, it is necessary to conduct a detailed study on the whole YYF and to reassess its future seismic risk.     

Composition of primary kimberlite melt in a garnet lherzolite mantle source: constraints from melting phase relations in anhydrous Udachnaya-East kimberlite with variable CO2 content at 6.5 GPa

The critical issue in the study of kimberlites, known as principal host rocks of diamonds, is the reconstruction of their primary melt composition, which is poorly constrained due to contamination by xenogenic materials, significant loss of volatiles during eruption, and post-magmatic alteration. It is generally accepted that the last equilibration of primary kimberlite melt with surrounding mantle (garnet lherzolite) occurred beneath cratons at 5–7 GPa (150–230 km depths). However, the subliquidus mineral assemblages obtained in kimberlite melting experiments at mantle pressures differ from lherzolite_, probably owing to unaccounted loss of CO₂. Here we present experiments at 6.5 GPa and 1200–1600 °C on unaltered kimberlite with an addition of 2–22 mol% CO₂ over its natural abundance in the rock (13 mol%), but keeping proportions of other components identical to those in an exceptionally fresh anhydrous kimberlite from Udachnaya-East pipe in Siberia. We found that the partial melt achieves equilibrium with garnet lherzolite at 1500 °C and 19–23 mol% CO₂ in the system. Under these conditions this melt contains (mol%): SiO₂ = 9, FeO = 6–7, MgO = 23–26, CaO = 16, Na₂O = 4, K₂O = 1, and CO₂ = 30–35. We propose, therefore, the alkali-rich carbonatitic composition of primary kimberlite melt and loss of 34–45 mol% (34–46 wt%) CO₂ during ascent of the kimberlite magma to the surface.     

Is earthquake activity along the French Atlantic margin favoured by local rheological contrasts?

The seismological study of recent seismic crises near Oleron Island confirms the coexistence of an extensional deformation and a transtensive regime in the Atlantic margin of France, which is different from the general western European stress field corresponding to a strike-slip regime. We argue that the switch of the principal stress axes σ₁/σ₂ in a NW–SE vertical plane is linked with the existence of crustal heterogeneities. Events of magnitude larger than 5 sometimes occur along the Atlantic margin of France, such as the 7 September 1972 (M_L = 5.2) earthquake near Oleron island and the 30 September 2002 (M_L = 5.7) Hennebont event in Brittany. To test the mechanism of local strain localization, we model the deformation of the hypocentral area of the Hennebont earthquake using a 3D thermo-mechanical finite element code. We conclude that the occurrence of moderate earthquakes located in limited parts of the Hercynian shear zones (as the often reactivated swarms near Oleron) could be due to local reactivation of pre-existing faults. These sporadic seismic ruptures are favoured by stress concentration due to rheological heterogeneities.     

First-order rupture features of the 2011 MW 9.0 Tohoku (Japan) earthquake from surface waves

This study analyzed the rupture directivity of the 2011 Tohoku earthquake by using 100-s Rayleigh-wave travel-times, influenced by the finite source, to derive the fault parameters of the earthquake. The results demonstrated that the earthquake exhibited a slow rupture propagation with a rupture velocity of approximately 1.5–2.0 km/s and asymmetric bilateral faulting. The two rupture directions were N60°E and N127°E, with rupture lengths of approximately 276 km and 231 km, respectively. The rupture toward N60°E had a source duration of approximately 183 s, longer than that toward N127°E (approximately 156 s). Overall, the entire source duration of the earthquake faulting lasted approximately 183 s. Regarding historical seismicity in eastern Japan, the 2011 Tohoku earthquake not only ruptured a locked area in which large earthquakes have rarely occurred, but also ruptured the source regions of several historical earthquakes. With the exception of its slow rupture velocity and generation of a tsunami, the rupture features of the 2011 Tohoku earthquake were inconsistent with those of typical tsunami earthquakes.     

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

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

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.     

The large aftershocks triggered by the 2011 Mw 9.0 Tohoku-Oki earthquake,Japan

The Mw 9.0 Tohoku-Oki earthquake that occurred off the Pacific coast of Japan on March 11, 2011, was followed by thousands of aftershocks, both near the plate interface and in the crust of inland eastern Japan. In this paper, we report on two large, shallow crustal earthquakes that occurred near the Ibaraki-Fukushima prefecture border, where the background seismicity was low prior to the 2011 Tohoku-Oki earthquake. Using densely spaced geodetic observations (GPS and InSAR datasets), we found that two large aftershocks in the Iwaki and Kita-Ibarake regions (hereafter referred to as the Iwaki earthquake and the Kita-Ibarake earthquake) produced 2.1 m and 0.44 m of motion in the line-of-sight (LOS), respectively. The azimuth-offset method was used to obtain the preliminary location of the fault traces. The InSAR-based maximum offset and trace of the faults that produced the Iwaki earthquake are consistent with field observations. The fault location and geometry of these two earthquakes are constrained by a rectangular dislocation model in a multilayered elastic half-space, which indicates that the maximum slips for the two earthquakes are 3.28 m and 0.98 m, respectively. The Coulomb stress changes were calculated for the faults following the 2011 Mw 9.0 Tohoku-Oki earthquake based on the modeled slip along the fault planes. The resulting Coulomb stress changes indicate that the stresses on the faults increased by up to 1.1 MPa and 0.7 MPa in the Iwaki and Kita-Ibarake regions, respectively, suggesting that the Tohoku-Oki earthquake triggered the two aftershocks, supporting the results of seismic tomography.     

Metamorphic P–T path of mafic granulites from Eastern Hebei: Implications for the Neoarchean tectonics of the Eastern Block,North China Craton

This study documents the metamorphic evolution of mafic granulites from the Eastern Hebei Complex in the Eastern Block of the North China Craton. Mafic granulites from Eastern Hebei occur as boudins or enclaves within Neoarchean high-grade TTG gneisses. Petrographic observations reveal three characteristic metamorphic mineral assemblages in the mafic granulites: the pre-peak hornblende + plagioclase + ilmenite + quartz + sphene assemblage (M₁) existing as mineral inclusions within coarse-grained peak assemblage (M₂) represented by garnet + clinopyroxene + orthopyroxene + plagioclase + hornblende + ilmenite + quartz, and post-peak assemblage (M₃) marked by garnet + quartz ± ilmenite symplectites surrounding the peak pyroxene and plagioclase. Based on pseudosection modeling calculated in the NCFMASHTO model system using the program THERMOCALC, P–T conditions of the pre-peak (M₁), peak (M₂) and post-peak (M₃) assemblages are constrained at 600–715 °C/6.0 kbar or below, 860–900 °C/9.6–10.3 kbar, and 790–810 °C/9.6–10.4 kbar, respectively. These P–T estimates, combined with their mineral compositions and reaction relations, define an anticlockwise P–T path incorporating isobaric cooling subsequent to the peak medium-pressure granulite-facies metamorphism for the mafic granulites from Eastern Hebei. Such an anticlockwise P–T path suggests that the end-Neoarchean metamorphism of the Eastern Hebei Complex correlated closely with underplating and intrusion of voluminous mantle-derived magmas. In conjunction with other geological considerations, a mantle-plume model is favored to interpret the Neoarchean tectonothermal evolution of the Eastern Hebei Complex and other metamorphic complexes in the Eastern Block. The prograde amphibolite-facies metamorphism (M₁) was initiated due to the upwelling of the relatively cooler mantle plume head, followed by the peak medium-pressure granulite-facies metamorphism (M₂) as triggered by the uprising hotter plume “tail”, and finally when plume activity ceased, the heated metamorphic crust experienced nearly isobaric cooling (M₃).     

Petrogenesis of the 2115 Ma Haicheng mafic sills from the Eastern North China Craton: Implications for an intra-continental rifting

The Rhyacian (2300–2050 Ma) is a special era of the Paleoproterozoic represented by large layered intrusions in many cratons. It is well known that there are widespread igneous events at ~ 2100 Ma in the Eastern North China Craton; however, their tectonic environments are under debate: whether they were related to an intra-continental rifting or an arc/back-arc setting along a continental margin. These ~ 2100 Ma igneous events comprise several mafic dykes/sills, with some coeval A-type granites and volcanic events in several rifts; among them, the Haicheng mafic sills in the Liaohe rift are unique as their host rock, the Liaohe Group, bears the world's largest magnesium deposit. Most of the mafic sills are E-W-elongated at present coordinates. Exclusive of superimposition caused by deformation, the widths of the individuals are tens to hundreds of meters and the lengths are hundreds to thousands of meters. They have metamorphosed to an assemblage of plagioclase and hornblende, with minor quartz and accessory chlorite, epidote, apatite, ilmenite, and magnetite. However, relic gabbro and ophitic textures with mainly plagioclase and clinopyroxene are well-preserved. SIMS Pb–Pb dating on baddeleyites from one ~ 1000 m thick sill near Xialiulinzi village yields an average ²⁰⁷Pb/²⁰⁶Pb age of 2115 ± 3 Ma (n = 15, MSWD = 2.3), representing the timing of crystallization. SIMS U–Pb dating on zircon yields a similar forming age. They are tholeiitic in composition (MgO: 4.36–8.88 wt.%; SiO₂: 45.76–53.39 wt.%), enriched in light rare earth elements ((La/Yb)_N = 1.72–4.37) and large ion lithophile elements (i.e., Cs, Rb, Sr, and K) but depleted in high field strength elements (i.e., Nb, Ta, and Ti). These features were unlikely caused by crustal contamination during their emplacement, as there are little variations in Nb/La and Th/Nb. The rocks have experienced significant plagioclase-plus clinopyroxene-dominating fractional crystallization. Their enriched Sr–Nd isotope characteristics (⁸⁷Sr/⁸⁶Sr_t = 0.703 ~ 0.705, εNd_t = − 1.9 ~ 0.6) and trace element patterns indicate that their source(s) could be the ancient subcontinental lithospheric mantle; and this source is similar to those coeval sills from other parts of the craton. Their arc-like trace element features could be inherited from their source regions formed via a subduction process at the late Archean rather than at the middle-late Paleoproterozoic. These sill swarms, throughout the craton, might have developed in an integrated intra-continental rift system at ~ 2100 Ma.     

The 10 June 2012 Fethiye Mw 6.0 aftershock sequence and its relation to the 24–25 April 1957 Ms 6.9–7.1 earthquakes in SW Anatolia,Turkey

The 10 June 2012 M_w 6.0 aftershock sequence in southwestern Anatolia is examined. Centroid moment tensors for 23 earthquakes with moment magnitudes (M_w) between 3.7 and 6.0 are determined by applying a waveform inversion method. The mainshock is a shallow focus strike-slip with reverse component event at a depth of 30 km. The seismic moment (M_o) of the mainshock is estimated as 1.28 × 10¹⁸ Nm and rupture duration of the Fethiye mainshock is 38 s. The focal mechanisms of the aftershocks are mainly strike-slip faulting with a reverse component. The geometry of the focal mechanisms reveals a strike-slip faulting regime with NE–SW trending direction of T-axis in the entire activated region. A stress tensor inversion of focal mechanism data is performed to obtain a more accurate picture of the Fethiye earthquake stress field. The stress tensor inversion results indicate a predominant strike-slip stress regime with a NW–SE oriented maximum horizontal compressive stress (S_H). According to variance of the stress tensor inversion, to first order, the Fethiye earthquake area is characterized by a homogeneous interplate stress field. The Coulomb stress change associated with the mainshock and the largest aftershock are also investigated to evaluate any significant enhancement of stresses along the Gulf of Fethiye and surrounding region. Positive lobes with stress more than 0.4 bars are obtained, indicating that these values are large enough to increase the Coulomb stress failure towards NNW–SSE and E–W directions.     

Ultramafic xenoliths from Lake Nyos area,Cameroon volcanic line,West-central Africa: Petrography,mineral chemistry,equilibration conditions and metasomatic features

Spinel-bearing mantle xenoliths have been recovered in the pyroclastic breccia surrounding the Lake Nyos maar. These include spinel lherzolites, spinel harzburgites and olivine websterites. They exhibit coarse granular or protogranular to weakly porphyroclastic textures, and show variations in mineral chemistry, modal compositions and equilibrium temperature. The xenoliths consist of four mineral phases typical of upper mantle origin: olivine (Fo₈₉–Fo_91.5, NiO = 0.29–0.38 wt%, CaO = 0.02–0.17 wt%), enstatite (Mg# = 90–92, Cr₂O₃ = 0.35 ± 0.04 wt%), Cr-diopside (Mg# = 92–98, Cr₂O₃ = 0.7–1.65 wt%, TiO₂ = 0.26–0.6 wt%) and spinel (high Mg# of 70–80, low TiO₂ ≤ 0.4 wt%). Spinels are aluminous (Cr# = 9.7–11) in most lherzolites, and become increasingly chromiferous from websterites (Cr#_Sp = 15.3–19.8) to harzburgites (Cr#_Sp = 19–33.6). The lherzolites are composed of olivine (48–58%), orthopyroxene (22–30%), and clinopyroxene (8–15%). The harzburgites modes are olivine (60–81%), orthopyroxene (11–29%), and clinopyroxene (<5%). The websterites are mainly composed of pyroxene (～62%) with variable amounts of olivine (23–31%). Temperatures of mineral equilibration in the xenoliths have been estimated from the two-pyroxene thermometer of Wells (1977) and range between 850 and 1050 °C, corresponding to about 10–30 kbar at a depth mantle of 30 km at least. These P–T conditions show significant variations between different petrographical types, the maximum conditions being recorded in two spinel lherzolites (NY-05 and NY-23) that have atypical chemical compositions and textures suggesting that they were initially formed in an environment close to the garnet stability field, then re-equilibrated within the spinel stability field prior to their incorporation in the host magma. With the exception of minerals from these two lherzolite nodules, all the minerals exhibit depletion of light REE, a typical feature of abyssal peridotites implying that some xenoliths from the Cameroon volcanic line were probably sampled in a part of the sub-continental mantle that is chemically similar to sub-oceanic mantle. The variations observed in the mineral chemistry and modal compositions of xenoliths suggest that the spinel harzbugite nodules which represent residues of a significant degree of partial melting of lherzolitic mantle were affected by infiltration of alkali-enriched metasomatizing melts (or fluids) within the uppermost mantle to produce pargasitic amphiboles prior to their sampling by the host lava. The features of this metasomatism event occur in the rocks of all three petrographical facies xenoliths from Lake Nyos.     

Paleomagnetic analysis of the Marwar Supergroup,Rajasthan, India and proposed interbasinal correlations

The Marwar Supergroup refers to a 1000–2000 m thick marine and coastal sequence that covers a vast area of Rajasthan in NW–India. The Marwar Basin uncomformably overlies the ∼750–770 Ma rocks of the Malani Igneous Suite and is therefore considered Late Neoproterozoic to Early Cambrian in age. Upper Vindhyan basinal sediments (Bhander and Rewa Groups), exposed in the east and separated by the Aravalli–Delhi Fold Belt, have long been assumed to coeval with the Marwar Supergroup. Recent studies based on detrital zircon populations of the Marwar and Upper Vindhyan sequences show some similarity in the older populations, but the Vindhyan sequence shows no zircons younger than 1000 Ma whereas samples taken from the Marwar Basin show distinctly younger zircons. This observation led to speculation that the Upper Vindhyan and Marwar sequences did not develop coevally.While there are alternative explanations for why the two basins may differ in their detrital zircon populations, paleomagnetic studies may provide independent evidence for differences/similarities between the assumed coeval basins. We have collected samples in the Marwar Basin and present the paleomagnetic results. Previous paleomagnetic studies of Marwar basinal sediments were misinterpreted as being indistinguishable from the Upper Vindhyan sequence. The vast majority of our samples show directional characteristics similar to the previously published studies. We interpret these results to be a recent overprint. A small subset of hematite-bearing rocks from the Jodhpur Formation (basal Marwar) exhibit directional data (Dec = 89° Inc = −1° α95 = 9°) that are distinct from the Upper Vindhyan pole and may offer additional support for temporally distinct episodes of sedimentation in these proximal regions. A VGP based upon our directional data is reported at 1°S 344°E (dp = 5°, dm = 9°). We conclude that the Marwar Supergroup developed near the close of the Ediacaran Period and is part of a larger group of sedimentary basins that include the Huqf Supergroup (Oman), the Salt-Range (Pakistan), the Krol–Tal belt (Himalayas) and perhaps the Molo Supergroup (Madagascar).