Pn wave velocity and Moho geometry in north eastern India

Pn velocity has been computed across the NE India and Moho geometry constrained, using regional earthquake travel times recorded by a network of 30 seismological stations operated during February-May 1993. Using an appropriate velocity model and the arrival times at the network stations, preliminary hypocentres of 16 regional earthquakes provided by NEIC were also improved. The average Pn wave velocity in NE India has been found to be 8.5 ±0.2 km/s. It varies from 8.3 to 8.5 km/s beneath the Shillong Plateau, Mikhir hills and Assam valley, which is significantly higher than those in other parts of India. The crustal thickness in NE India is also high, varying from 45–49 km under the Shillong plateau and the adjoining region to 55–65 km in the convergence zone. The presence of a thick crust and high Pn velocity suggests that the lithosphere in NE India is colder, as also indicated by the observed deeper level (45-51 km) seismicity of the region.     

青藏高原地壳与上地幔成层速度结构与深部层间物质的运移轨迹

在印度洋板块与欧亚板块碰撞、挤压作用下,促使深部物质重新分异、调整和运移,并导致了地壳的短缩增厚,而且造成了高原的整体隆升和深部壳、幔物质的侧向流展。基于青藏高原腹地和周边地域地壳与上地幔的成层速度结构,特别是其特异层序的展布研究表明,青藏高原地壳巨厚,但岩石圈却相对较薄;地壳中于深20±5km处存在一低速层,层速度为5.7±0.1km/s,厚度为8±2km;上地幔软流圈顶部深度为110±10km;下地壳与上地幔盖层物质以地壳低速层为上滑移面,以岩石圈漂曳的上地幔软流圈顶面为下滑移面,在印度洋板块N-NNE向力源作用下在同步运移,即形成了青藏高原腹地和周边地域特异的大陆地球动力学环境。     

Estimation of crustal discontinuities from reflected seismic waves recorded at Shillong and Mikir Hills Plateau,Northeast India

In this study, an attempt is made to determine seismic velocity structure of the crust and upper mantle beneath the Shillong-Mikir Hills Plateau in northeast India region. The principle of the technique is to relate seismic travel times with crustal thickness above the Conrad and Moho discontinuities. Broadband digital waveforms of the local earthquakes make a precise detection of the seismic phases possible that are reflected at these discontinuities. The results show that the Conrad discontinuity is at 18–20 (±0.5) km beneath the Shillong-Mikir Hills Plateau and the Moho discontinuity is at 30 ± 1.0 km beneath the Shillong Plateau and at 35 ± 1.0 km beneath the Mikir Hills.     

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

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

Depths to density discontinuities beneath the Adamawa Plateau region, Central Africa, from spectral analyses of new and existing gravity data

New gravity data from the Adamawa Uplift region of Cameroon have been integrated with existing gravity data from central and western Africa to examine variations in crustal structure throughout the region. The new data reveal steep northeast-trending gradients in the Bouguer gravity anomalies that coincide with the Sanaga Fault Zone and the Foumban Shear Zone, both part of the Central African Shear Zone lying between the Adamawa Plateau and the Congo Craton. Four major density discontinuities in the lithosphere have been determined within the lithosphere beneath the Adamawa Uplift in central Cameroon using spectral analysis of gravity data: (1) 7–13 km; (2) 19–25 km; (3) 30–37 km; and (4) 75–149 km. The deepest density discontinuities determined at 75–149 km depth range agree with the presence of an anomalous low velocity upper mantle structure at these depths deduced from earlier teleseismic delay time studies and gravity forward modelling. The 30–37 km depths agree with the Moho depth of 33 km obtained from a seismic refraction experiment in the region. The intermediate depth of 20 km obtained within region D may correspond to shallower Moho depth beneath parts of the Benue and Yola Rifts where seismic refraction data indicate a crustal thickness of 23 km. The 19–20 km depths and 8–12 km depths estimated in boxes encompassing the Adamawa Plateau and Cameroon Volcanic Line may may correspond to mid-crustal density contrasts associated with volcanic intrusions, as these depths are less than depths of 25 and 13 km, respectively, in the stable Congo Craton to the south.     

多玛-德庆-达孜断面壳幔密度结构特征

利用人机交互重震联合反演的方法研究了多玛-德庆-达孜断面的二维壳幔横向密度结构特征。模拟结果显示,剖面下的地壳内部大部分存在低速层,念青唐古拉山两侧的德庆、羊八井附近存在两条深大正断层,切割并抬升了其下的中地壳低密度层,低密度层整体被向上抬升5～10 km,使得念青唐古拉山深部表现为一个地垒构造。念青唐古拉山位于莫霍面由浅变深的缓坡上,向东逐渐变深。软流圈在念青唐古拉山下形状发生变化,表现为“上凸”特征。     

The Central Southern Alps (N. Italy) paleoseismic zone: a comparison between field observations and predictions of fault mechanics

The internal sectors of the Orobic Alps (Northern Italy) are characterised by Alpine age regional shortening showing a transition, through time, from plastic to brittle deformation. Thrust faults cut Alpine ductile folds and are marked by cataclasites and, locally, by pseudotachylytes, suggesting that motion was accommodated by seismic frictional slip. In the Eastern Orobic Alps the thrusting initiated at depths deeper than 10 km (the emplacement depth of the Adamello pluton) and possibly continued at shallower depths. This demonstrates that thrust motion occurred between 10 km depth and the brittle-ductile transition, i.e., at mid-crustal depths. The Orobic Alps exhumed paleoseismic zone shows different geometries along strike. In the central sectors of the Orobic Alps, thrust faults, associated with pseudotachylytes, have average dips around 40° and show no pervasive veining. Much steeper thrusts (dips up to about 85°) occur in the eastern Orobic Alps. In this area, faults are not associated with pervasive veining, i.e., fluid circulation was relatively scarce. This suggests that faulting did not occur with supralithostatic fluid pressure conditions. These reverse faults are severely misoriented (far too steep) for fault reactivation in a sublithostatic fluid pressure regime. We suggest that thrust motion likely started when the faults were less steep and that the faults were progressively rotated up to the present day dips. Domino tilting is probably responsible for this subsequent fault steepening, as suggested by a decrease of the steepness of thrust faults from north to south and by systematic rotations of previous structures consistently with tilting of thrust blocks. When the faults became inclined beyond the fault lock-up angle, no further thrusting was accommodated along them. At later stages regional shortening was accommodated by newly formed lower angle shear planes (dipping around 30–40°), consistently with predictions from fault mechanics.     

Three-dimensional P- and S-wave velocity structures beneath the Japan Islands obtained by high-density seismic stations by seismic tomography

We construct fine-scale 3D P- and S-wave velocity structures of the crust and upper mantle beneath the whole Japan Islands with a unified resolution, where the Pacific (PAC) and Philippine Sea (PHS) plates subduct beneath the Eurasian (EUR) plate. We can detect the low-velocity (low-V) oceanic crust of the PAC and PHS plates at their uppermost part beneath almost all the Japan Islands. The depth limit of the imaged oceanic crust varies with the regions. High-V_P/V_S zones are widely distributed in the lower crust especially beneath the volcanic front, and the high strain rate zones are located at the edge of the extremely high-V_P/V_S zone; however, V_P/V_S at the top of the mantle wedge is not so high. Beneath northern Japan, we can image the high-V subducting PAC plate using the tomographic method without any assumption of velocity discontinuities. We also imaged the heterogeneous structure in the PAC plate, such as the low-V zone considered as the old seamount or the highly seismic zone within the double seismic zone where the seismic fault ruptured by the earthquake connects the upper and lower layer of the double seismic zone. Beneath central Japan, thrust-type small repeating earthquakes occur at the boundary between the EUR and PHS plates and are located at the upper part of the low-V layer that is considered to be the oceanic crust of the PHS plate. In addition to the low-V oceanic crust, the subducting high-V PAC plate is clearly imaged to depths of approximately 250 km and the subducting high-V PHS zone to depths of approximately 180 km is considered to be the PHS plate. Beneath southwestern Japan, the iso-depth lines of the Moho discontinuity in the PHS plate derived by the receiver function method divide the upper low-V layer and lower high-V layer of our model at depths of 30–50 km. Beneath Kyushu, the steeply subducting PHS plate is clearly imaged to depths of approximately 250 km with high velocities. The high-V_P/V_S zone is considered as the lower crust of the EUR plate or the oceanic crust of the PHS plate at depths of 25–35 km and the partially serpentinized mantle wedge of the EUR plate at depths of 30–45 km beneath southwestern Japan. The deep low-frequency nonvolcanic tremors occur at all parts of the high-V_P/V_S zone—within the zone, the seaward side, and the landward side where the PHS plate encounters the mantle wedge of the EUR plate. We prove that we can objectively obtain the fine-scale 3D structure with simple constraints such as only 1D initial velocity model with no velocity discontinuity.     

Continuation of the San Andreas fault system into the upper mantle: Evidence from spinel peridotite xenoliths in the Coyote Lake basalt, central California

The Coyote Lake basalt, located near the intersection of the Hayward and Calaveras faults in central California, contains spinel peridotite xenoliths from the mantle beneath the San Andreas fault system. Six upper mantle xenoliths were studied in detail by a combination of petrologic techniques. Temperature estimates, obtained from three two-pyroxene geothermometers and the Al-in-orthopyroxene geothermometer, indicate that the xenoliths equilibrated at 970–1100 °C. A thermal model was used to estimate the corresponding depth of equilibration for these xenoliths, resulting in depths between 38 and 43 km. The lattice preferred orientation of olivine measured in five of the xenolith samples show strong point distributions of olivine crystallographic axes suggesting that fabrics formed under high-temperature conditions. Calculated seismic anisotropy values indicate an average shear wave anisotropy of 6%, higher than the anisotropy calculated from xenoliths from other tectonic environments. Using this value, the anisotropic layer responsible for fault-parallel shear wave splitting in central California is less than 100 km thick. The strong fabric preserved in the xenoliths suggests that a mantle shear zone exists below the Calaveras fault to a depth of at least 40 km, and combining xenolith petrofabrics with shear wave splitting studies helps distinguish between different models for deformation at depth beneath the San Andrea fault system.     

Source zones of the catastrophic Simushir earthquakes on November 15, 2006 (M w = 8.3) and January 13, 2007 (M w = 8.1) and the deep crust structure beneath the Middle Kuril segment

The data on catastrophic earthquakes with magnitudes of 8.3 and 8.1 that occurred in the Simushir Island area on November 15, 2006, and January 13, 2007, respectively, were compared with the results of land-sea deep seismic studies by different methods (deep seismic sounding, the correlation method of refracted waves, the earthquake converted-wave method, the common mid-point) in the Central Kuril segment. The structure of the Earth’s crust and the hypocentral zones of these earthquakes were analyzed. It was established that the hypocenter of the main shock of the first earthquake was located at the bend of the seismofocal zone under the island slope of the trench on the outer side of the subsiding lithospheric plate in the rapidly rising granulite-basite (ìbasalticî) crustal layer, which, at depths of 7–15 km, replaced the granulite-gneiss layer. This was accompanied by an increase of the seismic wave velocity from 6.4 to 7.1 km/s. The focus of the second earthquake was located beneath the axis of the deep-sea trench. The aftershocks were concentrated in two bands 60–120 km wide that extend along the trench, as well as in the third zone orthogonal to the island arc. It was shown that the epicenters of the earthquakes are linked with regional faults. The main shock of the first earthquake (November 15, 2006) was interpreted as a thrust fault and the second one (January 13, 2007) was attributed to a normal fault.     

Seismotectonics of the Nepal Himalaya from a local seismic network

The National Seismological Network of Nepal consists of 17 short period seismic stations operated since 1994. It provides an exceptional view of the microseismic activity over nearly one third of the Himalayan arc, including the only segment, between longitudes 78°E and 85°E, that has not produced any M>8 earthquakes over the last century. It shows a belt of seismicity that follows approximately the front of the Higher Himalaya with most of the seismic moment being released at depths between 10 and 20 km. This belt of seismicity is interpreted to reflect interseismic stress accumulation in the upper crust associated with creep in the lower crust beneath the Higher Himalaya. The seismic activity is more intense around 82°E in Far-Western Nepal and around 87°E in Eastern Nepal. Western Nepal, between 82.5 and 85°E, is characterized by a particularly low level of seismic activity. We propose that these lateral variations are related to segmentation of the Main Himalayan Thrust Fault. The major junctions between the different segments would thus lie at about 87°E and 82°E with possibly an intermediate one at about 85°E. These junctions seem to coincide with some of the active normal faults in Southern Tibet. Lateral variation of seismic activity is also found to correlate with lateral variations of geological structures suggesting that segmentation is a long-lived feature. We infer four 250–400 km long segments that could produce earthquakes comparable to the M=8.4 Bihar–Nepal earthquake that struck eastern Nepal in 1934. Assuming the model of the characteristic earthquake, the recurrence interval between two such earthquakes on a given segment is between 130 and 260 years.     

Crustal and upper mantle seismic structure of the Australian Plate, South Island, New Zealand

Seismic reflection and refraction data were collected west of New Zealand's South Island parallel to the Pacific–Australian Plate boundary. The obliquely convergent plate boundary is marked at the surface by the Alpine Fault, which juxtaposes continental crust of each plate. The data are used to study the crustal and uppermost mantle structure and provide a link between other seismic transects which cross the plate boundary. Arrival times of wide-angle reflected and refracted events from 13 recording stations are used to construct a 380-km long crustal velocity model. The model shows that, beneath a 2–4-km thick sedimentary veneer, the crust consists of two layers. The upper layer velocities increase from 5.4–5.9 km/s at the top of the layer to 6.3 km/s at the base of the layer. The base of the layer is mainly about 20 km deep but deepens to 25 km at its southern end. The lower layer velocities range from 6.3 to 7.1 km/s, and are commonly around 6.5 km/s at the top of the layer and 6.7 km/s at the base. Beneath the lower layer, the model has velocities of 8.2–8.5 km/s, typical of mantle material. The Mohorovicic discontinuity (Moho) therefore lies at the base of the second layer. It is at a depth of around 30 km but shallows over the south–central third of the profile to about 26 km, possibly associated with a southwest dipping detachment fault. The high, variable sub-Moho velocities of 8.2 km/s to 8.5 km/s are inferred to result from strong upper mantle anisotropy. Multichannel seismic reflection data cover about 220 km of the southern part of the modelled section. Beneath the well-layered Oligocene to recent sedimentary section, the crustal section is broadly divided into two zones, which correspond to the two layers of the velocity model. The upper layer (down to about 7–9 s two-way travel time) has few reflections. The lower layer (down to about 11 s two-way time) contains many strong, subparallel reflections. The base of this reflective zone is the Moho. Bi-vergent dipping reflective zones within this lower crustal layer are interpreted as interwedging structures common in areas of crustal shortening. These structures and the strong northeast dipping reflections beneath the Moho towards the north end of the (MCS) line are interpreted to be caused by Paleozoic north-dipping subduction and terrane collision at the margin of Gondwana. Deeper mantle reflections with variable dip are observed on the wide-angle gathers. Travel-time modelling of these events by ray-tracing through the established velocity model indicates depths of 50–110 km for these events. They show little coherence in dip and may be caused side-swipe from the adjacent crustal root under the Southern Alps or from the upper mantle density anomalies inferred from teleseismic data under the crustal root.     

http://www.sciencedirect.com/science/article/pii/S1674987112000849

New methods are presented for processing and interpretation of shallow marine differential magnetic data,including constructing maps of offshore total magnetic anomalies with an extremely high resolution of up to 1-2 nT,mapping weak anomalies of 5-10 nT caused by mineralization effects at the contacts of hydrocarbons with host rocks,estimating depths to upper and lower boundaries of anomalous magnetic sources,and estimating thickness of magnetic layers and boundaries of tectonic blocks. Horizontal dimensions of tectonic blocks in the so-called "seismic gap" region in the central Kuril Arc vary from 10 to 100 km,with typical dimensions of 25-30 km.The area of the "seismic gap" is a zone of intense tectonic activity and recent volcanism.Deep sources causing magnetic anomalies in the area are similar to the "magnetic belt" near Hokkaido. In the southern and central parts of Barents Sea,tectonic blocks with widths of 30-100 km,and upper and lower boundaries of magnetic layers ranging from depths of 10 to 5 km and 18 to 30 km are calculated.Models of the magnetic layer underlying the Mezen Basin in an inland part of the White Sea-Barents Sea paleorift indicate depths to the lower boundary of the layer of 12-30 km.Weak local magnetic anomalies of 2-5 nT in the northern and central Caspian Sea were identified using the new methods,and drilling confirms that the anomalies are related to concentrations of hydrocarbon.Two layers causing magnetic anomalies are identified in the northern Caspian Sea from magnetic anomaly spectra.The upper layer lies immediately beneath the sea bottom and the lower layer occurs at depths between 30-40 m and 150-200 m.     

Understanding the tectonic behaviour of the Shillong Plateau,India using remote sensing data

Shillong Plateau in India is tectonically and geologically interesting entity in the subducted front of Indian Plate below Burmese Plate to the southeast and Tibetan Plate to the north and associated with thrusts and shears along the plate boundaries. Horse-tail geometry in the foothills of the Arunachal Himalaya, east of Jia Bhareli river, associated with south-convex foothill ranges in the eastern Himalaya and exactly similar structural geometry in the eastern part of Shillong Plateau in Meghalaya seems to develop due to resistance received by the plateau in its eastward journey. Wide separation of Karbi Anglong Plateau and Shillong Plateau to the southeast as compared to northwestern part defines the shape of Kopili graben. Low seismic activity in southeastern part of Shillong Plateau might be related to stress released field generated by its clockwise rotation. Satellite derived images and digital elevation data from Landsat ETM+ and SRTM data shows that the central part of Shillong Plateau possesses young topography with strong structural fabrics along with relatively high topography aligning NE-SW following Kolkota-Pabna-Mymansingh High and if extended passes through western part of Arunachal Pradesh in eastern Himalayas. This alignment has been observed in Precambrian gneissic complex west of the Proterozoic intracratonic Shillong Basin. The epicentral plot for the period 1918 to 2009 shows their high concentration within the Shillong Plateau aligning along this trend. The active geodynamics of Shillong Plateau is reflected in its seismic activity pattern in relation with the structural fabrics, northward migration of the Brahmaputra in the north front of the Plateau and by shrinking pattern of Chandubi Lake in the Kulsi river catchment, a north-flowing tributary of the Brahmaputra in the north-central part of the plateau.     

Seismic refraction study of the continental edge off the eastern united states

Three long, strike-parallel, seismic-refraction profiles were made on the continental shelf edge, slope and upper rise off New Jersey during 1975. The shelf edge line lies along the axis of the East Coast Magnetic Anomaly (ECMA), while the continental rise line lies 80 km seaward of the shelf edge. Below the unconsolidated sediments (1.7–3.6 km/sec), high-velocity sedimentary rocks (4.2–6.2 km/sec) were found at depths of 2.6–8.2 km and are inferred to be cemented carbonates. Although multichannel seismic-reflection profiles and magnetic depth-to-source data predicted the top of oceanic basement at 6–8 km beneath the shelf edge and 10–11 km beneath the rise, no refracted events occurred as first arrivals from either oceanic basement (layer 2, approximately 5.5 km/ sec) or the upper oceanic crust (layer 3A, approximately 6.8 km/sec). Second arrivals from 10.5 km depth beneath the shelf edge are interpreted as events from a 5.9 km/sec refractor within igneous basement. Other refracted events from either layers 2 or 3A could not be resolved within the complex second arrivals. A well-defined crustal layer with a compressional velocity of 7.1–7.2 km/sec, which can be interpreted as oceanic layer 3B, occurred at 15.8 km depth beneath the shelf and 12.9 km beneath the upper rise. A well-reversed mantle velocity of 8.3 km/sec was measured at 18–22 km depth beneath the upper continental rise. Comparison with other deep-crustal profiles along the continental edge of the Atlantic margin off the United States, specifically in the inner magnetically quiet zone, indicates that the compressional wave velocities and layer depths determined on the U.S.G.S. profiles are very similar to those of nearby profiles. This suggests that the layers are continuous and that the interpretation of the oceanic layer 3B under the shelf edge east of New Jersey implies progradation of the shelf outward over the oceanic crust in that area. This agrees with magnetic anomaly evidence which shows the East Coast Magnetic Anomaly landward of the shelf edge off New Jersey and with previous seismic reflection data which reveal extensive outbuilding of the shelf edge during the Jurassic and Lower Cretaceous, probably by carbonate bank-margin accretion.     

Mohorovičić discontinuity beneath Fennoscandia and adjacent parts of the Norwegian Sea and the North Sea

A preliminary contour map showing the Mohorovičić discontinuity (Moho) beneath Fennoscandia, adjacent parts of the Norwegian Sea and the North Sea has been compiled on the basis of published information from deep seismic soundings.The Moho contour map shows a 10 km thick crust beneath the investigated basin-region of the Norwegian Sea. It seems that the Vøring Plateau has at least in part a continental structure even if the Moho-depth is only 15 km. A shallow Moho (28–30 km) all along the Norwegian coast is a well established feature. A good correlation between the surface elevation of the mountain range running through Norway and parts of Sweden and the depth of the Moho is also well established. The Gulf of Bothnia is a region of a great Mono-depression.     

Moho depth variation beneath southwestern Japan revealed from the velocity structure based on receiver function inversion

The Philippine Sea plate is subducting under the Eurasian plate beneath the Chugoku-Shikoku region, southwestern Japan. We have constructed depth contours for the continental and oceanic Mohos derived from the velocity structure based on receiver function inversion. Receiver functions were calculated using teleseismic waveforms recorded by the high-density seismograph network in southwestern Japan. In order to determine crustal velocity structure, we first improved the linearized time-domain receiver function inversion method. The continental Moho is relatively shallow ( 30 km) at the coastline of the Sea of Japan and at the Seto Inland Sea, and becomes deeper–greater than 40 km–around 35°N and 133.8°E. Near the Seto Inland Sea, a low-velocity layer of thickness 10 km lies under the continental Moho. This low-velocity layer corresponds to the subducting oceanic crust of the Philippine Sea plate. The oceanic Moho continues to descend from south to northwest and exhibits complicated ridge and valley features. The oceanic Moho runs around 25 km beneath the Pacific coast and 45 km beneath the Seto Inland Sea, and it extends to at least to 34.5°N. The depth variation of the Moho discontinuities is in good qualitative agreement with the concept of isostasy. From the configurations of both the continental and oceanic Mohos, we demonstrate that the continental lower crust and the subducting oceanic crust overlap beneath the southern and central part of Shikoku and that a mantle wedge may exist beneath the western and eastern part of Shikoku. The southern edge of the overlapping region coincides with the downdip limit of the slip area of a megathrust earthquake.     

Seismogenesis of the uninterrupted occurrence of the aftershock activity in the 2001 Bhuj earthquake zone, Gujarat, India, during 2001�C2010

Seismogenesis of aftershocks occurring in the Kachchh seismic zone for more than last 10?years is investigated through modeling of fractal dimensions, b-value, seismic velocities, stress inversion, and Coulomb failure stresses, using aftershock data of the 2001 Bhuj earthquake. Three-dimensional mapping of b-values, fractal dimensions, and seismic velocities clearly delineate an area of high b-, D-, and Vp/Vs ratio values at 15?C35?km depth below the main rupture zone (MRZ) of the 2001 mainshock, which is attributed to higher material heterogeneities in the vicinity of the MRZ or deep fluid enrichment due to the release of aqueous fluid/volatile CO₂ from the eclogitisation of the olivine-rich lower crustal rocks. We notice that several aftershocks are occurred near the contacts between high (mafic brittle rocks) and low velocity regions while many of the aftershocks including the 2001 Bhuj mainshock are occurred in the zones of low velocity (low dVp, low dVs and large Vp/Vs) in the 15?C35?km depth range, which are inferred to be the fractured rock matrixes filled with aqueous fluid or volatiles containing CO₂. Further support for this model comes from the presence of hydrous eclogitic layer at sub-lithospheric depths (34?C42?km). The depth-wise stress inversions using the P- and T-axes data of the focal mechanisms reveal an increase in heterogeneity (i.e., misfit) with an almost N?CS ??₁ orientation up to 30?km depth. Then, the misfit decreases to a minimum value in the 30?C40?km depth range, where a 60^o rotation in the ??₁ orientation is also noticed that can be explained in terms of the fluid enrichment in that particular layer. The modeling of Coulomb failure stress changes (??CFS) considering three tectonic faults [i.e., NWF, GF, and Allah bund fault (ABF)] and the slip distribution of the 2001 mainshock on NWF could successfully explain the occurrences of moderate size events (during 2006?C2008) in terms of increase in positive ??CFS on GF and ABF. In a nutshell, we propose that the fluid-filled mafic intrusives are acting as stress accentuators below the Kachchh seismic zone, which generate crustal earthquakes while the uninterrupted occurrence of aftershocks is triggered by stress transfer and aqueous fluid or volatile CO₂ flow mechanisms. Further, our results on the 3-D crustal seismic velocity structure, focal mechanisms, and b-value mapping will form key inputs for understanding wave propagation and earthquake hazard-related risk associated with the Kachchh basin.     

Crustal velocity inhomogeneities along the Hirapur–Mandla profile, central India and its tectonic implications

Wide-angle seismic and gravity data across the Narmada-Son lineament (NSL) in central India are analyzed to determine crustal structure, velocity inhomogeneities and hence constrain the tectonics of the lineament. We present the 2-D crustal velocity structure from deep wide-angle reflection data by using a ray-trace inverse approach. The main result of the study is the delineation of fault-bounded horst raised to a subsurface depth (1.5 km) and the Moho upwarp beneath the NSL. The crust below the basement consists of three layers with velocities of 6.45–6.7, 6.2–6.5 and 6.7–6.95 km/s and interface depths of about 5.5–8.7, 14–17 and 18–23 km along the profile. The low-velocity (6.2–6.5 km/s) layer goes up to a depth of 5 km and becomes the thickest part (13 km), while the overlying high-velocity (6.45–6.7 km/s) layer becomes the thinnest (3 km) and upper boundary lies at a depth of 1.5 km beneath the NSL. The overall uncertainties of various velocity and boundary nodes are of the order of ±0.12 km/s and ±1.40 km, respectively. The up-lifted crustal block and the up-warping Moho beneath the NSL indicate that the north and south faults bounding the NSL are deeply penetrated through which mafic materials from upper mantle have been intruded into the upper crust. Gravity modeling was also undertaken to assess the seismically derived crustal features and to fill the seismic data gap. The lateral and vertical heterogeneous nature of the structure and velocity inhomogeneities in the crust cause instability to the crustal blocks and played an important role in reactivation of the Narmada south fault during the 1997 Jabalpur earthquake.     

Mantle dynamics of Western Pacific and East Asia: Insight from seismic tomography and mineral physics

Recent results of high-resolution seismic tomography and mineral physics experiments are used to study mantle dynamics of Western Pacific and East Asia. The most important processes in subduction zones are the shallow and deep slab dehydration and the convective circulation (corner flow) processes in the mantle wedge. The combination of the two processes may have caused the back-arc spreading in the Lau basin, affected the morphology of the subducting Philippine Sea slab and its seismicity under southwest Japan, and contributed to the formation of the continental rift system and intraplate volcanism in Northeast Asia, which are clearly visible in our tomographic images. Slow anomalies are also found in the mantle under the subducting Pacific slab, which may represent (a) small mantle plumes, (b) upwellings associated with the slab collapsing down to the lower mantle, or (c) sub-slab dehydration associated with deep earthquakes caused by the reactivation of large faults preserved in the slab. Combining tomographic images and earthquake hypocenters with phase diagrams in the systems of peridotite + water, we proposed a petrologic model for arc volcanism. Arc magmas are caused by the dehydration reactions of hydrated slab peridotite that supply water-rich fluids to the mantle wedge and cause partial melting of the convecting mantle wedge. A large amount of fluids can be released from hydrated MORB at depths shallower than 55 km, which move upwards to hydrate the wedge corner under the fore-arc, and never drag down to the deeper mantle along the slab surface. Slab dehydration reactions at 120 km depth are the antigorite-related 5 reactions which supply water-rich fluids for forming the volcanic front. Phase A and Mg-surssasite breakdown reactions at 200 and 300 km depths below 700 °C cause the second and third arcs, respectively. Moreover, the dehydration reactions of super-hydrous phase B, phases D and E at 500–660 km depths cause the fluid transportation to the mantle boundary layer (MBL) (410–660 km depth). The stagnant slabs extend from Japan to Beijing, China for over 1000 km long, indicating that the arc–trench system covers the entire region from the Japan trench to East Asia. We propose a big mantle wedge (BMW) model herein, where hydrous plumes originating from 410 km depth cause a series of intra-continental hot regions. Fluids derived from MBL accumulated by the double-sided subduction zones, rather than the India–Asia collision and the subsequent indentation into Asia, are the major cause for the active tectonics and mantle dynamics in this broad region.