Complex slab structure and arc magmatism beneath the Japanese Islands

A dense nationwide seismic network recently constructed in Japan has resulted in the production of a large amount of high-quality data that have enabled the high-resolution imaging of deep seismic structures in the Japanese subduction zone. Seismic tomography, precise locations of earthquakes, and focal mechanism research have allowed the identification of the complex structure of subducting slabs beneath Japan, revealing that the subducting Philippine Sea slab underneath southwestern Japan has an undulatory configuration down to a depth of 60–200 km, and is continuous from Kanto to Kyushu without disruption or splitting, even within areas north of the Izu Peninsula. Analysis of the geometry of the Pacific and Philippine Sea slabs identified a broad contact zone beneath the Kanto Plain that causes anomalously deep interplate and intraslab earthquake activity. Seismic tomographic inversions using both teleseismic and local events provide a clear image of the deep aseismic portion of the Philippine Sea slab beneath the Japan Sea north of Chugoku and Kyushu, and beneath the East China Sea west of Kyushu down to a depth of ∼450 km. Seismic tomography also allowed the identification of an inclined sheet-like seismic low-velocity zone in the mantle wedge beneath Tohoku. A recent seismic tomography work further revealed clear images of similar inclined low-velocity zones in the mantle wedge for almost all other areas of Japan. The presence of the inclined low-velocity zones in the mantle wedge across the entirety of Japan suggests that it is a common feature to all subduction zones. These low-velocity zones may correspond to the upwelling flow portion of subduction-induced convection systems. These upwelling flows reach the Moho directly beneath active volcanic areas, suggesting a link between volcanism and upwelling.     

P-wave tomography and origin of the Changbai intraplate volcano in Northeast Asia

We present the first seismic image of the upper mantle beneath the active intraplate Changbai volcano in Northeast Asia determined by teleseismic travel time tomography. The data are measured at a new seismic network consisting of 19 portable stations and 3 permanent stations. Our results show a columnar low-velocity anomaly extending to 400-km depth with a P-wave velocity reduction of up to 3%. High velocity anomalies are visible in the mantle transition zone, and deep-focus earthquakes occur at depths of 500–600 km under the region, suggesting that the subducting Pacific slab is stagnant in the transition zone, as imaged clearly by global tomography. These results suggest that the intraplate Changbai volcano is not a hotspot like Hawaii but a kind of back-arc volcano related to the deep subduction and stagnancy of the Pacific slab under Northeast Asia.     

Amount of Asian lithospheric mantle subducted during the India/Asia collision

Body wave seismic tomography is a successful technique for mapping lithospheric material sinking into the mantle. Focusing on the India/Asia collision zone, we postulate the existence of several Asian continental slabs, based on seismic global tomography. We observe a lower mantle positive anomaly between 1100 and 900 km depths, that we interpret as the signature of a past subduction process of Asian lithosphere, based on the anomaly position relative to positive anomalies related to Indian continental slab. We propose that this anomaly provides evidence for south dipping subduction of North Tibet lithospheric mantle, occurring along 3000 km parallel to the Southern Asian margin, and beginning soon after the 45 Ma break-off that detached the Tethys oceanic slab from the Indian continent. We estimate the maximum length of the slab related to the anomaly to be 400 km. Adding 200 km of presently Asian subducting slab beneath Central Tibet, the amount of Asian lithospheric mantle absorbed by continental subduction during the collision is at most 600 km. Using global seismic tomography to resolve the geometry of Asian continent at the onset of collision, we estimate that the convergence absorbed by Asia during the indentation process is ~ 1300 km. We conclude that Asian continental subduction could accommodate at most 45% of the Asian convergence. The rest of the convergence could have been accommodated by a combination of extrusion and shallow subduction/underthrusting processes. Continental subduction is therefore a major lithospheric process involved in intraplate tectonics of a supercontinent like Eurasia.     

An overview of electrical conductivity structures of the crust and upper mantle beneath the northwestern Pacific, the Japanese Islands, and continental East Asia

This article reviews the electrical conductivity structures of the oceanic upper mantle, subduction zones, and the mantle transition zone beneath the northwestern Pacific, the Japanese Islands, and continental East Asia, which have particularly large potential of water circulation in the global upper mantle. The oceanic upper mantle consists of an electrically resistive lid and a conductive layer underlying the lid. The depth of the top of the conductive layer is related to lithospheric cooling in the older mantle, whereas it is attributable to the difference in water distribution beneath the vicinity of the seafloor spreading-axis. The location of a lower crustal conductor in a subduction zone changes according to the subduction type. The difference can be explained by the characteristic dehydration from the subducting slab in each subduction zone and by advection from the backarc spreading. The latest one-dimensional electrical conductivity model of the mantle transition zone beneath the Pacific Ocean predicts values of 0.1–1.0 S/m. These values support a considerably dry oceanic mantle transition zone. However, one-dimensional electrical profiles may not be representative of the mantle transition zone there, since there exists a three-dimensional structure caused by the stagnant slab. Three-dimensional electromagnetic modeling should be made in future studies.     

The making and breaking of supercontinents: Some speculations based on superplumes, super downwelling and the role of tectosphere

The mechanisms of formation and disruption of supercontinents have been topics of debate. Based on the Y-shaped topology, we identify two major types of subduction zones on the globe: the Circum-Pacific subduction zone and the Tethyan subduction zone. We propose that the process of formation of supercontinents is controlled by super downwelling that develops through double-sided subduction zones as seen in the present day western Pacific, and also as endorsed by both geologic history and P-wave whole mantle tomography. The super-downwelling swallows all material like a black hole in the outer space, pulling together continents into a tight assembly. The fate of supercontinents is dictated by superplumes (super-upwelling) which break apart the continental assemblies. We evaluate the configuration of major supercontinents through Earth history and propose the tectonic framework leading to the future supercontinent Amasia 250 million years from present, with the present day Western Pacific region as its frontier. We propose that the tectosphere which functions as the buoyant keel of continental crust plays a crucial role in the supercontinental cycle, including continental fragmentation, dispersion and amalgamation. The continental crust is generally very thin, only about one tenth of the thickness of the tectosphere. If the rigidity and buoyancy is derived from the tectosphere, with the granitic upper crust playing only a negligible role, then supercontinent cycle may reflect the dispersion and amalgamation of the tectosphere. Therefore, supercontinent cycle may correspond to super-tectosphere cycle.     

Mantle dynamics and mantle melting beneath Niuafo’ou Island and the northern Lau back-arc basin

A suite of young volcanic basaltic lavas erupted on the intra-plate island of Niuafo’ou and at active rifts and spreading centres (the King’s Triple Junction and the Northeastern Lau Spreading Centre) in the northern Lau Basin is used to examine the pattern of mantle flow and the dynamics of melting beneath this complex back-arc system. All lavas contain variable amounts of a subduction related component inherited from the Tonga subduction zone to the east. All lavas have higher ⁸⁷Sr/⁸⁶Sr, lower ¹⁴³Nd/¹⁴⁴Nd and more radiogenic Pb isotope compositions than basalts erupted at the Central Lau Spreading Centre in the central Lau Basin, and are interpreted as variable mixtures of subduction-modified, depleted upper mantle, and mantle residues derived from melting beneath the Samoan Islands which has leaked through a tear in the subducting Pacific Plate beneath the Vitiaz Lineament at the northern edge of the Lau Basin. Our data can be used to map out the present-day distribution of Samoan mantle in this region, and show that it influences the compositions of lavas erupted as far as 400 km from the Samoan Islands. The distribution of Samoan-influenced lavas implies south- and southwest-wards mantle flow rates of >4 cm/year. U-series disequilibria in historic Niuafo’ou lavas have average (²³⁰Th/²³⁸U) = 1.13, (²³¹Pa/²³⁵U) = 2.17, (²²⁶Ra/²³⁰Th) = 2.11, and together with major and trace element data require ∼5% partial melting of mantle at between 2 and 3 GPa, with a residual porosity of 0.002 and an upwelling rate of 1 cm year⁻¹. We suggest that intraplate magmatism in the northern Lau Basin results from decompression melting during southward flow of mantle from beneath old (110–120 Ma), relatively thick Pacific oceanic lithosphere to beneath young (<5 Ma), thinner oceanic lithosphere beneath the northern Lau Basin.     

Seismic anisotropy tomography: New insight into subduction dynamics

Body-wave and surface-wave tomography, receiver-function imaging, and shear-wave splitting measurements have shown that seismic anisotropy and heterogeneity coexist in all parts of subduction zones, providing important constraints on the mantle flow and subduction dynamics. P-wave anisotropy tomography is a new but powerful tool for mapping three-dimensional variations of azimuthal and radial seismic anisotropy in the crust and mantle. P-wave azimuthal-anisotropy tomography has been applied widely to the Circum-Pacific subduction zones, Mainland China and North America, whereas P-wave radial-anisotropy tomography was applied to only a few areas including Northeast Japan, Southwest Japan and North China Craton. These studies have revealed complex anisotropy in the crust and mantle lithosphere associated with the surface geology and tectonics, anisotropy reflecting subduction-driven corner flow in the mantle wedge, frozen-in fossil anisotropy in the subducting slabs formed at the mid-ocean ridge, as well as olivine fabric transitions due to changes in water content, stress and temperature. Shear-wave splitting tomography methods have been also proposed, but their applications are still limited and preliminary. There is a discrepancy between the surface-wave and body-wave tomographic models in radial anisotropy of the mantle wedge beneath Japan, which is a puzzle but an intriguing topic for future studies.     

Seismic anisotropy in the wedge above the Philippine Sea slab beneath Kanto and southwest Japan derived from shear wave splitting

We conduct shear wave splitting measurements on waveform data from the Hi-net and the broadband F-net seismic stations in Kanto and SW Japan generated by shallow and intermediate-depth earthquakes occurring in the subducting Philippine Sea and Pacific slabs. We obtain 1115 shear wave splitting parameter pairs. The results are divided into those from the shallow (depth < 50 km) and the deep (depth > 50 km) events. The deep events beneath Kanto are further divided into PHS1 and PHS2 (upper and lower planes of the double seismic zone in the Philippine Sea slab, respectively), PAC1 and PAC2 (western and eastern Pacific slab, respectively) events. The results from the shallow events represent the crustal anisotropy, and their fast directions are more or less aligned in the σ_Hmax directions, implying that the anisotropy is produced by the alignment of the vertical cracks in the crust induced by the compressive stresses. In Kanto, Kii Peninsula and Kyushu regions, the results from the deep events suggest a contribution from the mantle wedge anisotropy. Events from all groups beneath Kanto show NW, NE and EW fast directions. This complex pattern seems to be produced by the corner flows induced by both the WNW PAC plate subduction and the oblique NNW PHS slab subduction with the associated olivine lattice-preferred orientations (LPOs), and the anisotropy frozen in the PHS slab. The deep events beneath Kii Peninsula show NE and NW fast directions and may be produced by the corner flow produced by the NNW PHS slab subduction with the associated olivine LPOs. The NE directions might also be produced by the segregated melts in the thin layers parallel to the PHS slab subduction. The deep events beneath N Kyushu show NNW fast directions, which may result from the southeastward flow in the upper mantle inferred from the stresses in the upper plate. Results from the deep events beneath middle-south Kyushu show dominantly E–W fast directions, in both the fore- and back-arcs. They may be produced by the corner flow of the westward PHS slab subduction with the olivine LPOs. Because the source regions with multiple fast directions are not resolved in this study, further detailed analyses of shear wave splitting are necessary for a better understanding of the stress state, the induced mantle flow, and the melt-segregation processes.     

An overview of electrical conductivity structures of the crust and upper mantle beneath the northwestern Pacific,the Japanese Islands,and continental East Asia

This article reviews the electrical conductivity structures of the oceanic upper mantle, subduction zones, and the mantle transition zone beneath the northwestern Pacific, the Japanese Islands, and continental East Asia, which have particularly large potential of water circulation in the global upper mantle. The oceanic upper mantle consists of an electrically resistive lid and a conductive layer underlying the lid. The depth of the top of the conductive layer is related to lithospheric cooling in the older mantle, whereas it is attributable to the difference in water distribution beneath the vicinity of the seafloor spreading-axis. The location of a lower crustal conductor in a subduction zone changes according to the subduction type. The difference can be explained by the characteristic dehydration from the subducting slab in each subduction zone and by advection from the backarc spreading. The latest one-dimensional electrical conductivity model of the mantle transition zone beneath the Pacific Ocean predicts values of 0.1–1.0 S/m. These values support a considerably dry oceanic mantle transition zone. However, one-dimensional electrical profiles may not be representative of the mantle transition zone there, since there exists a three-dimensional structure caused by the stagnant slab. Three-dimensional electromagnetic modeling should be made in future studies.     

Configuration of the subducting Philippine Sea slab in the eastern part of southwestern Japan with seismic array and Hi-net data

It is important to know the shape of a subducting slab in order to understand the mechanisms of inter-plate earthquakes and the process of subduction. Seismicity data and converted phases have been used to detect plate boundaries. The configuration of the Philippine Sea slab has been obtained at the western part of southwestern Japan. At the eastern part of southwestern Japan, however, the configuration of the Philippine Sea slab has not yet been confirmed. A spatially high-density seismic network makes it possible to detect the boundaries of the Philippine Sea slab. We used a spatially high-density temporal seismic array in the area. The configuration of the Philippine Sea plate is obtained at the eastern part of southwestern Japan using the temporal seismic array and permanent seismic network data and comparing the seismic structure obtained from the results of refraction surveys. The configuration of the Philippine Sea plate obtained by this study does not bend sharply compared to previous models obtained from receiver function analyses. We delineated the upper boundary of the slab to a depth of about 45 km. The weak image of the boundary, which corresponds to the upper mantle reflector beneath the source area of the 2000 Western Tottori earthquake, was detected using the spatially dense array.     

武夷山成矿带及邻区上地幔速度结构及其对燕山期岩浆-成矿活动的启示

武夷山成矿带燕山期岩浆-成矿活动的深部动力学机制一直是学者们研究的热点.已有的研究结果表明,武夷山成矿带及邻区的上地幔存在着显著的低速异常,可能与地表的岩浆-成矿活动存在着密切的联系.本研究利用分布在华夏地块98个固定地震台站以及59个流动地震台站所记录到的278个远震事件,采用远震层析成像方法构建了武夷山成矿带及邻区...     

长江中下游成矿带及邻区地壳结构——MASH成矿过程的P波接收函数成像证据？

利用长江中下游成矿带多学科深部探测剖面于2009年11月至2011年3月间采集的天然地震数据,通过天然地震接收函数成像等分析研究,得到了研究区地壳和上地幔结构的清晰图像。接收函数成像结果显示研究区内Moho面深度存在着明显的起伏变化,在长江中下游成矿带(指剖面穿过的长江中下游成矿带宁芜矿集区,下同)下方存在着"幔隆构造"。在剖面东南端(即扬子克拉通北缘),Moho面相对稳定,深度约为30km;在茅山和江南断裂附近,Moho面存在上下起伏现象;在剖面中部或宁芜矿集区下方,Moho面存在明显隆起,深度只有28km;在郯庐断裂带下方,Moho面明显加深,深度达到36km;进一步向北到华北地台南缘,Moho面深度逐渐恢复到了32km左右的平均深度水平。其次,我们在接收函数成像结果中发现,长江中下游成矿带与其周边下地壳结构存在着明显的差异,成矿带的下地壳具有显著的地震波方位各向异性。扬子克拉通北缘的下地壳呈高速的近水平状结构,地震波各向异性特征不明显;与此相比,长江中下游成矿带的下地壳虽然也呈近水平状结构特征,但是,对于沿成矿带走向方向传播的地震波,其下地壳具有高速特征,而对于垂直于成矿带走向方向上传播的地震波,其下地壳却又表现为低速特征,这意味着成矿带的下地壳存在着平行于成矿带走向(即近北东—南西)方向的地震波各向异性,我们解释其是下地壳熔融并沿成矿带走向水平流动导致矿物晶体定向排列的结果。最后,在郯庐断裂以西的华北地台南缘观测到一条从上地壳延伸到中下地壳的南南东向倾斜的转换震相,我们推测它可能是合肥盆地内地壳伸展构造的反映。此外,我们发现接收函数成像结果中观测到的"幔隆构造"与远震P波层析成像结果在成矿带下方150km深度上显示的上地幔低速异常(江国明等,另文发表)存在着良好的对应关系,我们解释它们是软流圈物质上涌的遗迹。综合天然地震接收函数成像、远震P波层析成像和前人关于岩浆岩等方面的研究成果,我们认为长江中下游成矿带现今的下地壳可能是中生代发生成矿作用的多级岩浆房系统的一部分,成矿带的形成可能是类似MASH过程的产物。首先,软流圈物质上涌导致了长江中下游成矿带及其周边拉张环境的形成,在其上部地壳中形成了一系列伸展构造;然后,软流圈物质通过底侵进入长江中下游成矿带的原下地壳并与原下地壳物质发生同化作用,形成类埃达克质岩浆;接着,类埃达克质岩浆沿着伸展、拆离构造上升到地壳浅部形成不同层次的岩浆房和侵入岩体,并与围岩作用形成矿床。     

中国天然地震探测深部构造的研究现状与发展

在青藏高原等地区为研究下地壳上地幔的构造特征,最近10a投入了较大工作量的天然地震探测。它可以获得较大深度上的地震速度图像、上地幔的各向异性,划分出低速体低速层,为研究地球动力学特征提供了新的资料。随着仪器的发展完善,这种方法在我国深部研究中成为重要而且经济的手段。地震层析、接收函数的研究、以及剪切波分裂的分析都取得了很大的成绩。在这些方面还有许多工作要深入研究,继续深入发展具有重要意义。     

Fossil flat-slab subduction beneath the Illinois basin, USA

The Illinois basin is one of several well-studied intracratonic sedimentary basins within the North American craton whose formational mechanisms and subcrustal structure are not well understood. We study the S-velocity structure of the upper mantle beneath the Illinois basin and its surrounding area through seismic tomography. We utilize continental scale waveform data of seismic S and surface waves, enhanced by regional earthquakes located near the Illinois basin. Our 3D tomographic model, IL05, confirms the existence of a slow S-velocity structure in the uppermost mantle beneath the Illinois basin region. This anomalously slow region exists from the base of the crust to depths of  90 km, and is slower than the North American cratonic average by about 200 m/s. This anomalous uppermost mantle beneath the Illinois basin is underlain by a faster lithosphere, typical of the surrounding craton, to depths of  200 km. Excluding the formation of the Reelfoot Rift, this area of North American has been stable for over 1.0 Gy. Thus, we do not expect thermal anomalies from before that time to persist into present day S-velocity anomalies and we consider a delamination origin as an explanation of Illinois basin subsidence unlikely. We cannot rule out that the slow mid-lithosphere beneath the Illinois basin is caused by an uppermost mantle enriched by a deep, but weak plume. We attribute the slow mid-lithosphere to the presence of either oceanic, hydrous crust, or, a relatively cool mantle wedge with preserved hydrous minerals in the Illinois basin's uppermost mantle, related to a fossilized flat subduction zone.     

Seismic structure and origin of active intraplate volcanoes in Northeast Asia

Three-dimensional P-wave velocity structure beneath the Changbai and other intraplate volcanic areas in Northeast Asia is determined by inverting 1378 high-quality P-wave arrival times from 186 teleseismic events recorded by 61 broadband seismic stations. Low-velocity (low-V) anomalies are revealed beneath the Changbai, Longgan, Xianjindao volcanoes. High-velocity (high-V) anomalies are found in the mantle transition zone, where deep-focus earthquakes under Hunchun occur at depths of 500–600 km. The high-V anomaly reflects the deep subduction of the Pacific slab under NE Asia which may have contributed to the formation of the Changbai, Longgang, Xianjindao and Jingpohu intraplate volcanoes. A low-V anomaly is also revealed in the mantle transition zone, which may have a close relationship with the occurrence of deep earthquakes under the Hunchun area. Our results support the Big Mantle Wedge (BMW) model by Zhao et al. [Zhao, D., Lei, J., Tang, Y., 2004. Origin of the Changbai volcano in northeast China: evidence from seismic tomography, Chin. Sci. Bull. 49, 1401–1408; Zhao, D., Maruyama, S., Omori, S., 2007. Mantle dynamics of western Pacific and East Asia: insight from seismic tomography and mineral physics. Gondwana Res. 11, 120–131.] who proposed that the intraplate volcanoes in NE Asia are caused by the back-arc magmatism associated with the deep dehydration process of the subducting slab and convective circulation process in the BMW above the stagnant Pacific slab.     

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.     

The structure of the crust and upper mantle to depths of 320 km beneath the Kaapvaal craton, from P wave arrivals generated by regional earthquakes and mining-induced tremors

Travel times from earthquakes recorded at two seismic networks were used to derive an average P wavespeed model for the crust and upper mantle to depths of 320 km below southern Africa. The simplest model (BPI1) has a Moho depth of 34 km, and an uppermost mantle wavespeed of 8.04 km/s, below which the seismic wavespeeds have low positive gradients. Wavespeed gradients decrease slightly around 150 km depth to give a ‘knee’ in the wavespeed-depth model, and the wavespeed reaches 8.72 km/s at a depth of 320 km. Between the Moho and depths of 270 km, the seismic wavespeeds lie above those of reference model IASP91 of Kennett [Research School of Earth Sciences, Australian National University, Canberra, Australia (1991)] and below the southern African model of Zhao et al. [Journal of Geophysical Research 104 (1999) 4783]. At depths near 300 km all three models have similar wavespeeds. The mantle P wavespeeds for southern Africa of Qiu et al. [Geophysical Journal International 127 (1996) 563] lie close to BPI1 at depths between 40 and 140 km, but become lower at greater depths. The seismic wavespeeds in the upper mantle of model BPI1 agree satisfactorily with those estimated from peridotite xenoliths in kimberlites from within the Kaapvaal craton.The crustal thickness of 34 km of model BPI1 is systematically lower than the average thickness of 41 km computed over the same region from receiver functions. This discrepancy can be partly explained by an alternative model (BPI2) in which there is a crust–mantle transition zone between depths of 35 and 47 km, below which seismic wavespeed increases to 8.23 km/s. A low-wavespeed layer is then required at depths between 65 and 125 km.     

South African seismicity, April 1997 to April 1999, and regional variations in the crust and uppermost mantle of the Kaapvaal craton

Events induced by deep gold-mining activity on the edge of the Witwatersrand basin dominate the seismicity of South Africa. The deployment of 54 broad-band seismic stations at 84 separate locations across southern Africa between April 1997 and April 1999 (Kaapvaal network) enabled the seismicity of South Africa to be better defined over a 2-year period. Seismic events located by the South African national network, and by localized seismic networks deployed in mines or across gold-mining areas, were used to evaluate earthquake location procedures and to show that the Kaapvaal network locates mining-induced tremors with an average error of 1.56±0.10 km compared with 9.50±0.36 km for the South African network. Travel times of seismic events from the mines recorded at the Kaapvaal network indicate regional variations in the thickness of the crust but no clearly resolved variations in seismic wavespeeds in the uppermost mantle. Greater average crustal thicknesses (48–50 km compared with 41–43 km) are observed in the northern parts of the Kaapvaal craton that were affected by the Bushveld magmatism at 2.05 Ga. Estimates of average crustal thickness for the southern part of the Kaapvaal craton from receiver functions (38 km) agree well with those from refracted arrivals from mining-induced earthquakes if the crustal thicknesses below the sources are assumed to be 40–43 km. In contrast, the average crustal thickness inferred from refracted arrivals for the northern part of the Kaapvaal craton is larger by about 7 km (51 km) than that inferred from receiver functions (44 km), suggesting a thick mafic lower crust of variable seismic properties due to variations in composition and metamorphic grade. Pn wavespeeds are high (8.3–8.4 km/s), indicating the presence of highly depleted magnesium-rich peridotite throughout the uppermost mantle of the craton. Seismic Pg and Sg phases indicate that the upper crust around the Witwatersrand basin is comparatively uniform in composition when averaged over several kilometres.     

Magnetotelluric characterization through the Ambargasta-Sumampa Range: The connection between the northern and southern trace of the Río de La Plata Craton – Pampean Terrane tectonic boundary

The South American Platform was part of the Western Gondwana, a collage of plates of different ages assembled in late Neoproterozoic to Cambrian times. The Transbrasiliano Lineament, a continental shear belt that transversely intersects this platform from NE to SW, has its southern expression in the tectonic boundary between the Río de La Plata Craton and the Pampean Terrane. Magnetotelluric long-period data in a W–E profile (29°30′ S) that crosses the Ambargasta-Sumampa Range and the Chaco-Pampean Plain were obtained to connect information of this mostly inferred tectonic boundary. A 2-D inversion model shows the Chacoparanense basin, Río Dulce lineament, Ambargasta-Sumampa Range and Salina de Ambargasta in the upper crust. At mid-to-lower crust and 40 km to the east of the Ambargasta-Sumampa Range, a discontinuity (500–2000 Ω m) of 20-km-wide separates two highly resistive blocks, the Río de La Plata Craton (6000–20,000 Ω m) in the east, and the Pampean Terrane (5000–20,000 Ω m) in the west. This discontinuity represents the tectonic boundary between both cratons and could be explained by the presence of graphite. The geometry of the Pampean Terrane suggests an east-dipping paleo-subduction. Our results are consistent with gravimetric and seismicity data of the study area. A more conductive feature beneath the range and the tectonic boundary was associated with the NE–SW dextral transpressive system evidenced by the mylonitic belts exposed in the Eastern Pampean Ranges. This belt represents a conjugate of the mega-shear Transbrasiliano Lineament and could be explained by fluid–rock interaction by shearing during hundreds of years. The eastern border of the Ambargasta-Sumampa Range extends the trace of the Transbrasiliano Lineament. The electrical Moho depth (40 km to the west and 35 km to the east) was identified by a high electrical contrast between the crust and upper mantle. The upper mantle shows a resistive structure beneath the Río de La Plata Craton that could have been originated by stationary delamination by the presence of hydrated lithosphere.     

Upper mantle P-wave tomography across the Longmenshan fault belt from passive-source seismic observations along Aba-Longquanshan profile

A passive seismic experiment across the Longmenshan (LMS) fault belt had been conducted between August 2006 and July 2007 for the understanding of geodynamic process between the Eastern Tibet and Sichuan basin. We herein collected 3677 first P arrival times with high precision from seismograms of 288 teleseismic events so as to reconstruct the upper mantle velocity structure. Our results show that the depth of the Lithosphere–asthenosphere boundary (LAB) changes from 70 km beneath Eastern Tibet to about 110 km beneath Longquanshan, Sichuan Basin, which is consistent with the receiver function imaging results. The very thin mantle part of the lithosphere beneath Eastern Tibet may suggest the lithosphere delamination due to strong interaction between the Tibetan eastward escaping flow and the rigid resisting Sichuan basin, which can be further supported by the existences of two high-velocity anomalies beneath LAB in our imaging result. We also find there are two related low-velocity anomalies beneath the LMS fault belt, which may indicate magmatic upwelling from lithosphere delamination and account for the origin of tremendous energy needed by the devastating Wenchuan earthquake.