India-Madagascar paleo-fit based on flexural isostasy of their rifted margins

The present study contributes new constraints on, and definitions of, the reconstructed plate margins of India and Madagascar based on flexural isostasy along the Western Continental Margin of India (WCMI) and the Eastern Continental Margin of Madagascar (ECMM). We have estimated the nature of isostasy and crustal geometry along the two margins, and have examined their possible conjugate structure. Here we utilize elastic thickness (Te) and Moho depth data as the primary basis for the correlation of these passive margins. We employ the flexure inversion technique that operates in spatial domain in order to estimate the spatial variation of effective elastic thickness. Gravity inversion and flexure inversion techniques are used to estimate the configuration of the Moho/Crust-Mantle Interface that reveals regional correlations with the elastic thickness variations. These results correlate well with the continental and oceanic segments of the Indian and African plates. The present study has found a linear zone of anomalously low-Te (1-5 km) along the WCMI (similar to 1680 km), which correlates well with the low-Te patterns obtained all along the ECMM. We suggest that the low-Te zones along the WCMI and ECMM represent paleo-rift inception points of lithosphere thermally and mechanically weakened by the combined effects of the Marion hotspot and lithospheric extension due to rifting. We have produced an India-Madagascar paleo-fit representing the initial phase of separation based on the Te estimates of the rifted conjugate margins, which is confirmed by a close-fit correlation of Moho geometry and bathymetry of the shelf margins. The matching of tectonic lineaments, lithologies and geochronological belts between India and Madagascar provide an additional support for the present plate reconstruction. (C) 2014 International Association for Gondwana Research. Published by Elsevier B.V. All rights reserved.

The thermal and mechanical structure of a two-dimensional plume in the earths mantle

On the condition that the distribution of velocity and temperature at the mid-plane of a mantle plume has been obtained (pages 213–218, this issue), the problem of determining the lateral structure of the plume at a given depth is reduced to solving an eigenvalue problem of a set of ordinary differential equations with five unknown functions, with an eigenvalue being related to the thermal thickness of the plume at this depth. The lateral profiles of upward velocity, temperature and viscosity in the plume and the thickness of the plume at various depths are calculated for two sets of Newtonian rheological parameters. The calculations show that the precondition for the existence of the plume, δT/L 1 (L = the height of the plume, δT = lateral distance from the mid-plane), can be satisfied, except for the starting region of the plume or near the base of the lithosphere. At the lateral distance, δT, the upward velocity decreases to 0.1 – 50% of its maximum value at different depths. It is believed that this model may provide an approach for a quantitative description of the detailed structure of a mantle plume.

大洋中脊的形成过程

本文从热传导方程出发,得到了大洋中脊下岩石层温度分布的分析表达式及数值计算结果。结果表明,软流层上涌流动所提供的热源可以使大洋中脊下岩石层逐步融化;岩石层的相对移动速度对大洋中脊岩石层温度场及融化深度影响较大。 

Plate tectonic constraints on flat subduction and paleomagnetic constraints on rifting

Plate tectonics shapes our dynamic planet through the creation and destruction of lithosphere. This work focuses on increasing our understanding of the processes at convergent and divergent boundaries through geologic and geophysical observations at modern plate boundaries. Recent work had shown that the subducting slab in central Mexico is most likely the flattest on Earth, yet there was no consensus about what caused it to originate. The first chapter of this thesis sets out to systematically test all previously proposed mechanisms for slab flattening on the Mexican case. What we have discovered is that there is only one model for which we can find no contradictory evidence. The lack of applicability of the standard mechanisms used to explain flat subduction in the Mexican example led us to question their applications globally. The second chapter expands the search for a cause of flat subduction, in both space and time. We focus on the historical record of flat slabs in South America and look for a correlation between the shallowing and steepening of slab segments with relation to the inferred thickness of the subducting oceanic crust. Using plate reconstructions and the assumption that a crustal anomaly formed on a spreading ridge will produce two conjugate features, we recreate the history of subduction along the South American margin and find that there is no correlation between the subduction of a bathymetric highs and shallow subduction. These studies have proven that a subducting crustal anomaly is neither a sufficient or necessary condition of flat slab subduction. The final chapter in this thesis looks at the divergent plate boundary in the Gulf of California. Through geologic reconnaissance mapping and an intensive paleomagnetic sampling campaign, we try to constrain the location and orientation of a widespread volcanic marker unit, the Tuff of San Felipe. Although the resolution of the applied magnetic susceptibility technique proved inadequate to contain the direction of the pyroclastic flow with high precision, we have been able to detect the tectonic rotation of coherent blocks as well as rotation within blocks.

Modelagem sedimentar de sistemas deltaicos altamente eficientes: aplicação para o preenchimento sedimentar do terciário superior da bacia de Bengala, porção nordeste da Índia

A bacia de Bengala, localizada a Nordeste da Índia tem uma história evolutiva extraordinária, diretamente controlada bela fragmentação do Gondwana. O início da formação desta bacia é considerada como sendo relacionada ao final do evento da quebra, datado em 126 Ma quando a Índia separou do continente Antártico e da Austrália. Desde então, a placa continental Indiana viajou do pólo sul a uma velocidade muito rápida (16 cm/a) chocando-se com o hemisfério norte e fundindo-se com a Placa Eurasiana. Durante a viagem passou por cima de um hot spot, onde hoje estão localizadas as ilhas Seicheles, resultando em um dos maiores derrames de lava basáltica do mundo, conhecido como Deccan Trap. Na região onde a bacia de Bengala foi formada, não houve aporte significativo de sedimentos siliciclásticos, resultando na deposição de uma espessa plataforma carbonática do Cretáceo tardio ao Eoceno. Após este período, devido a colisão com algumas microplacas e a amalgamação com a Placa Eurasiana, um grande volume sedimentar siliciclástico foi introduzido para a bacia, associado também ao soerguimento da cadeia de montanhas dos Himalaias. Atualmente, a Bacia de Bengala possui mais de 25 km de sedimentos, coletados neste depocentro principal. Nesta dissertação foram aplicados conceitos básicos de sismoestratigrafia na interpretação de algumas linhas regionais. As linhas sísmicas utilizadas foram adquiridas recentemente por programa sísmico especial, o qual permitiu o imageamento sísmico a mais de 35km dentro da litosfera (crosta continental e transicional). O dado permitiu interpretar eventos tectônicos, como a presença dos Seawards Dipping Reflectors (SDR) na crosta transicional, coberto por sedimentos da Bacia de Bengala. Além da interpretação sísmica amarrada a alguns poços de controle, o programa de modelagem sedimentar Beicip Franlab Dionisos, foi utilizado para modelar a história de preenchimento da bacia para um período de 5,2 Ma. O nível relativo do nível do mar e a taxa de aporte sedimentar foram os pontos chaves considerados no modelo. Através da utilização dos dados sísmicos, foi possível reconhecer dez quebras de plataformas principais, as quais foram utilizadas no modelo, amarrados aos seus respectivos tempos geológico, provenientes dos dados dos poços do Plioceno ao Holoceno. O resultado do modelo mostrou que a primeira metade modelada pode ser considerada como um sistema deposicional retrogradacional, com algum picos transgressivos. Este sistema muda drasticamente para um sistema progradacional, o qual atuou até o Holoceno. A seção modelada também mostra que no período considerado o total de volume depositado foi em torno de 2,1 x 106 km3, equivalente a 9,41 x 1014 km3/Ma.

Petrogênese das rochas máficas alcalinas do litoral entre São Sebastião (SP) e Parati (RJ)

Lamprófiros e diabásios alcalinos afloram no litoral dos estados de São Paulo e Rio de Janeiro e integram o Enxame de Diques da Serra do Mar (EDSM). Essas rochas ocorrem sob a forma de diques e intrudem o Orógeno Ribeira, de idade Neoproterozóica/Cambro-Ordoviciana, inserindo-se no contexto geodinâmico de abertura do Oceano Atlântico Sul durante o Cretáceo Superior. Essas intrusões são subverticais e orientam-se preferencialmente a NE-SW, seguindo a estruturação das rochas encaixantes. Os lamprófiros são classificados como monchiquitos e camptonitos e exibem, respectivamente, textura hipocristalina e holocristalina. Apresentam também textura panidiomórfica, fenocristais de clinopiroxênio e olivina, imersos em matriz formada essencialmente por esses mesmos minerais, além de biotita, kaersutita e minerais opacos. O camptonito apresenta ainda plagioclásio na matriz. Os diabásios alcalinos são hipocristalinos a holocristalinos, equigranulares a inequigranulares, com fenocristais de olivina e/ou clinopiroxênio e/ou plagioclásio, em uma matriz composta essencialmente por esses minerais. As rochas estudadas caracterizam séries alcalinas miaskíticas, com os lamprófiros sendo tanto sódicos, potássicos e ultrapotássicos e os diabásios alcalinos como predominantemente sódicos. Modelagens petrogenéticas envolvendo possíveis processos evolutivos mostram que é improvável que os lamprófiros sejam cogenéticos por processos evolutivos envolvendo tanto cristalização fracionada, com ou sem assimilação concomitante, quanto hibridização. O mesmo ocorre para os diabásios alcalinos. A discriminação de fontes mantélicas foi feita com base nos teores de elementos traços de amostras representativas de líquidos parentais e indica que esse magmatismo alcalino está relacionado a fontes lherzolíticas com fusão parcial na zona de estabilidade do espinélio, isto é, a poucas profundidades. Os dados litogeoquímicos e isotópicos do sistema Sr-Nd das rochas estudadas sugerem tanto o envolvimento de fontes férteis, associadas ao manto sublitosférico, quanto de fontes enriquecidas, relacionadas ao manto litosférico subcontinental. Modelagens de mistura binária revelam que a petrogênese dos lamprófiros e diabásios alcalinos envolveu uma grande participação de um componente fértil misturado com contribuições menores de um componente enriquecido. Idades TDM (760-557 Ma) obtidas sugerem remobilização do manto litosférico no Neoproterozóico, talvez relacionadas à subducção da Placa São Francisco preteritamente à colisão do Orógeno Ribeira. Altas razões CaO/Al2O3 para os líquidos lamprofíricos menos evoluídos, altos teores de Zr, correlações negativas Zr/Hf e Ti/Eu e associação com carbonatitos indicam condições metassomáticas de alto CO2/H2O. Em escala local, modelos geodinâmicos baseados na astenosfera não isotérmica parecem mais aplicáveis. No entanto, modelos geodinâmicos baseados na astenosfera isotérmica (com o envolvimento de plumas) parecem mais indicados num contexto regional, considerando-se outras províncias alcalinas contemporâneas e correlatas.

Geochemical and isotopic characteristics of volcanic rocks from the northern East China Sea shelf margin and the Okinawa Trough

Volcanic rocks both from the northern East China Sea (NECS) shelf margin and the northern Okinawa Trough are subalkaline less aluminous, and lower in High Field Strength Elements (HFSE). These rocks are higher in Large Ion Lithophile Elements (LILE), thorium and uranium contents, positive lead anomalies, negative Nb-Ta anomalies, and enrichment in Light Rare Earth Elements (LREE). Basalts from the NECS shelf margin are akin to Indian Ocean Mid-Ocean Ridge Basalt (MORB), and rhyolites from the northern Okinawa Trough have the highest Pb-207/Pb-208 and Pb-208/Pb-204 ratios. The NECS shelf margin basalts have lower Sr-87/Sr-86 ratios, epsilon(Nd) and sigma O-18 than the northern Okinawa Trough silicic rocks. According to K-40-Ar-40 isotopic ages of basalts from the NECS shelf margin, rifting of the Okinawa Trough may have been active since at least 3.65-3.86 Ma. The origin of the NECS shelf margin basalt can be explained by the interaction of melt derived from Indian Ocean MORB-like mantle with enriched subcontinental lithosphere. The basalts from both sides of the Okinawa Trough may have a similar origin during the initial rifting of the Okinawa Trough, and the formation of basaltic magmas closely relates to the thinning of continental crust. The source of the formation of the northern Okinawa Trough silicic rocks was different from that of the middle Okinawa Trough, which could have been generated by the interaction of basaltic melt with an enriched crustal component. From the Ryukyu island arc to East China, the Cenozoic basalts have apparently increasing trends of MgO contents and ratios of LREE to Heavy Rare Earth Elements (HREE), suggesting that the trace element variabilities of basalts may have been influenced by the subduction of the Philippine Sea plate, and that the effects of subduction of the Philippine Sea plate on the chemical composition of basaltic melts have had a decreasing effect from the Ryukyu island arc to East China.

Rifting process and formation mechanisms of syn-rift stage prolongation in the deepwater basin, northern South China Sea

Based on the latest seismic and geological data, tectonic subsidence of three seismic lines in the deepwater area of Pearl River Mouth Basin (PRMB), the northern South China Sea (SCS), is calculated. The result shows that the rifting process of study area is different from the typical passive continental margin basin. Although the seafloor spreading of SCS initiated at 32 Ma, the tectonic subsidence rate does not decrease but increases instead, and then decreases at about 23 Ma, which indicates that the rifting continued after the onset of seafloor spreading until about 23 Ma. The formation thickness exhibits the same phenomenon, that is the syn-rift stage prolonged and the post-rift thermal subsidence delayed. The formation mechanisms are supposed to be three: (1) the lithospheric rigidity of the northern SCS is weak and its ductility is relatively strong, which delayed the strain relaxation resulting from the seafloor spreading; (2) the differential layered independent extension of the lithosphere may be one reason for the delay of post-rift stage; and (3) the southward transition of SCS spreading ridge during 24 to 21 Ma and the corresponding acceleration of seafloor spreading rate then triggered the initiation of large-scale thermal subsidence in the study area at about 23 Ma.

青岛晚白垩世岩石圈结构特征：劈石口晚白垩世基性脉岩及其捕虏体提供的证据

The occurrence of Late Cretaceous mafic dykes and their entrained peridotite and granulite xenoliths as well as clinopyroxene xenocrysts in the Qingdao region provide us a precious opportunity to unveil the nature and characteristics of the Late Mesozoic lithospheric mantle and lower crust beneath the Jiaodong region, and the change of the magma sources. These studies are of important and significant for understanding the lithospheric evolution in the eastern North China Craton. There were two periods of magma activities in Late Mesozoic in Qingdao Laoshan region, one was around 107Ma in the Early Cretaceous and the other around 86Ma in the Late Cretaceous according to the whole rock K-Ar age determination. The Early Cretaceous mafic dykes and the Late Cretaceous mafic dyke (i.e. Pishikou mafic dike) have completely different geochemical characteristics. The Early Cretaceous mafic dykes are enriched in LILE, strongly depleted in HFSE (Nb, Ta, Zr, Hf) and characterized by the highly radiogenic Sr and Nd isotopic compositions. These geochemical features indicate that the Early Cretaceous mafic dykes were derived from an enriched lithospheric mantle. In contrast, the Late Cretaceous mafic dyke is enriched in LILE, without HFSE depletion (Nb, Ta, Zr, Hf) and has less radiogenic Nd and Sr isotopic compositions. These geochemical features indicate that the Late Cretaceous mafic dyke was derived from the asthenosphere modified by subducted pelagic sediment contamination. The intrusive age of the Late Cretaceous mafic dyke provides further information for the termination of the lithosphere thinning for the eastern North China Crtaon. Pishikou Late Cretaceous mafic dyke contains abundant peridotitic xenoliths, granulite xenoliths and clinopyroxene xenocrysts. The peridotitic xenoliths can be divided into two types: high Mg# peridotites and low Mg# peridotites, according to their textural and mineral features. The high-Mg# peridotites have high Fo (up to 92.2) olivines and high Cr# (up to 55) spinels. The clinopyroxenes in the high# peridotites are rich in Cr2O3 and poor in Al2O3. The low-Mg# peridotites are typified by low Mg# (Fo <90) in olivines and low Cr# (Cr# <0.14) in spinels. The clinopyroxenes in the low-Mg# peridotites are rich in Al2O3 and Na2O and poor in Cr2O3. These two type peridotites have similar equilibrated temperatures of 950C-1100C. The Clinopyroxenes in the high-Mg# peridotites generally have high and variable REE contents (REE = 5.6-84 ppm) and LREE-enriched chondrite-normalized patterns ((La/Yb)N>1). In contrast, the clinopyroxenes in the low-Mg# peridotites have low REE contents (REE = 12 ppm) and LREE-depleted patterns ((La/Yb)N<1). The textural, mineral and elemental features of the low-Mg# peridotites are similar to those of the low-Mg peridotites from the Junan, representing the newly-accreted lithospheric mantle. However, the mineralogical and petrological features of the high-Mg# peridotites are similar to those of the high-Mg# peridotites from the Junan region, representing samples from the old refractory lithospheric mantle that was strongly and multiply affected by melts of different origins Late Cretaceous mafic dike in the Qingdao region also contains two types of granulite xenoliths according to the mineral constituents: the pyroxene-rich granulites and the plagioclase-rich granulites. Equilibrated temperatures calculated from the cpx-opx geothermometers are in a range of 861C - 910C for the pyroxene-rich granulites and of 847C - 890C for the plagioclase-rich granulites. The equilibrated pressure for the plagioclase-rich granulites is in a range of 9.9-11.7 kbar. Combined with the results of the peridotitic xenoliths, a 40C temperature gap exists between the peridotite and the granulite. The petrological Moho was 33~36 km at depths, broadly consistent with the seismic Moho estimated from the geophysical data. This indicates that there was no obvious crust-mantle transition zone in the Qingdao region in the Late Mesozoic. Pishikou Late Cretaceous mafic dyke entrained lots of clinopyroxene xenocrysts which are characterized by the chemical zoning. According to the zoning features, two types of clinopyroxene xenoliths can be classified, the normal zoning and the revise zoning. The normally-zoned clinopyroxene xenocrysts have LREE-depleted REE patterns in the cores. In contrast, the revisely-zoned clinopyroxenes have LREE-enriched REE patterns in the cores. According to the rim and core compositions of xenocrysts, all the rims are balanced with the host magma. Meanwhile, the origins of the cores were complicated, in which the normally-zoned clinopyroxenes were derived form the lithospheric mantle and some of the reversely-zoned clinopyroxnes were originated from the lower crust. Other revisely-zoned clinopyroxenes had experienced complex geological evolution and need to be further investigated. According to the above results, a simplified lithospheric profile has been established beneath the Qingdao region and a constraint on the nature and characteristics of the lithospheric mantle and lower crust has been made.

青藏高原及邻区面波频散反演--岩石圈-软流层体系结构研究

The Qinghai-Tibet Plateau lies in the place of the continent-continent collision between Indian and Eurasian plates. Because of their interaction the shallow and deep structures are very complicated. The force system forming the tectonic patterns and driving tectonic movements is effected together by the deep part of the lithosphere and the asthenosphere. It is important to study the 3-D velocity structures, the spheres and layers structures, material properties and states of the lithosphere and the asthenosphere for getting knowledge of their formation and evolution, dynamic process, layers coupling and exchange of material and energy. Based on the Rayleigh wave dispersion theory, we study the 3-D velocity structures, the depths of interfaces and thicknesses of different layers, including the crust, the lithosphere and the asthenosphere, the lithosphere-asthenosphere system in the Qinghai-Tibet Plateau and its adjacent areas. The following tasks include: (1)The digital seismic records of 221 seismic events have been collected, whose magnitudes are larger than 5.0 over the Qinghai-Tibet Plateau and its adjacent areas. These records come from 31 digital seismic stations of GSN , CDSN、NCDSN and part of Indian stations. After making instrument response calibration and filtering, group velocities of fundamental mode of Rayleigh waves are measured using the frequency-time analysis (FTAN) to get the observed dispersions. Furthermore, we strike cluster average for those similar ray paths. Finally, 819 dispersion curves (8-150s) are ready for dispersion inversion. (2)From these dispersion curves, pure dispersion data in 2°×2° cells of the areas (18°N-42°N, 70°E-106°E) are calculated by using function expansion method, proposed by Yanovskaya. The average initial model has been constructed by taking account of global AK135 model along with geodetic, geological, geophysical, receiving function and wide-angle reflection data. Then, initial S-wave velocity structures of the crust and upper mantle in the research areas have been obtained by using linear inversion (SVD) method. (3)Taking the results of the linear inversion as the initial model, we simultaneously invert the S wave velocities and thicknesses by using non-linear inversion (improved Simulated Annealing algorithm). Moreover, during the temperature dropping the variable-scale models are used. Comparing with the linear results, the spheres and layers by the non-linear inversion can be recognized better from the velocity value and offset. (4)The Moho discontinuity and top interface of the asthenosphere are recognized from the velocity value and offset of the layers. The thicknesses of the crust, lithosphere and asthenosphere are gained. These thicknesses are helpful to studying the structural differentia between the Qinghai-Tibet Plateau and its adjacent areas and among geologic units of the plateau. The results of the inversion will provide deep geophysical evidences for studying deep dynamical mechanism and exploring metal mineral resource and oil and gas resources. The following conclusions are reached by the distributions of the S wave velocities and thicknesses of the crust, lithosphere and asthenosphere, combining with previous researches. (1)The crust is very thick in the Qinghai-Tibet Plateau, varying from 60 km to 80 km. The lithospheric thickness in the Qinghai-Tibet Plateau is thinner (130-160 km) than its adjacent areas. Its asthenosphere is relatively thicker, varies from 150 km to 230 km, and the thickest area lies in the western Qiangtang. India located in south of Main Boundary thrust has a thinner crust (32-38 km), a thicker lithosphere of about 190 km and a rather thin asthenosphere of only 60 km. Sichuan and Tarim basins have the crust thickness less than 50km. Their lithospheres are thicker than the Qinghai-Tibet Plateau, and their asthenospheres are thinner. (2)The S-wave velocity variation pattern in the lithosphere-asthenosphere system has band-belted distribution along east-westward. These variations correlate with geology structures sketched by sutures and major faults. These sutures include Main Boundary thrust (MBT), Yarlung-Zangbo River suture (YZS), Bangong Lake-Nujiang suture (BNS), Jinshajiang suture (JSJS), Kunlun edge suture (KL). In the velocity maps of the upper and middle crust, these sutures can be sketched. In velocity maps of 250-300 km depth, MBT, BNS and JSJS can be sketched. In maps of the crustal thickness, the lithospheric thickness and the asthenospheric thickness, these sutures can be still sketched. In particular, MBT can be obviously resolved in these velocity maps and thickness maps. (3)Since the collision between India and Eurasian plate, the “loss” of surface material arising from crustal shortening is caused not only by crustal thickening but also by lateral extrusion material. The source of lateral extrusion lies in the Qiangtang block. These materials extrude along the JSJS and BNS with both rotation and dispersion in Daguaiwan. Finally, it extends toward southeast direction. (4)There is the crust-mantle transition zone of no distinct velocity jump in the lithosphere beneath the Qiangtang Terrane. It has thinner lithosphere and developed thicker asthenosphere. It implies that the crust-mantle transition zone of partial melting is connected with the developed asthenosphere. The underplating of asthenosphere may thin the lithosphere. This buoyancy might be the main mechanism and deep dynamics of the uplift of the Qinghai-Tibet hinterland. At the same time, the transport of hot material with low velocity intrudes into the upper mantle and the lower crust along cracks and faults forming the crust-mantle transition zone.

利用台阵面波资料研究华北克拉通北缘燕山褶皱带及其邻区的地壳与上地幔结构

The continent of eastern China, especially the North China Craton (NCC), has endured intensive tectonic renovation during Mesozoic and Cenozoic, with the presence of widespread magmatism, high heat flow and development of large sedimentary basins and mountain ranges. The cratonic lithosphere of the region has been destroyed remarkably, which is characterized by not only a significant reduction in thickness but also complex modifications in physical and chemical properties of the lithosphere. As for the tectonic regime controlling the evolution of the NCC, various models have been put forward, including the impingement of mantle plumes (“mushroom cloud” model), the collision of south China block and north China block, the subduction of the Pacific plate, etc. Lithosphere delamination and thermal erosion were proposed as the two end-member mechanisms of the lithospheric thinning. However, given the paucity of the data, deep structural evidence is currently still scarce for distinguishing and testifying these models. To better understand the deep structure of the NCC, from 2000 to the present, temporary seismic array observations have been conducted in the NCC by the Seismological Laboratory of the Institute of the Geology and Geophysics, Chinese Academy of Sciences under the North China Interior Structure Project (NCISP). Many arrays extend from the North China Craton and the off-craton regions, and traverse a lot of main tectonic boundaries. A total of more than 300 broadband seismic stations have been deployed along several profiles that traversed the major tectonic units within the craton’s interior, at the boundary areas and in the neighboring off-craton regions. These stations recorded abundant high-quality data, which provides an unprecedented opportunity for us to unravel the deep structural features of the NCC using seismological methods. Among all the seismological methods, the surface wave method appears to be an efficient and widely adopted technique in studying the crustal and upper mantle structures. In particular, it can provide the absolute values of S-wave velocity that are difficult to obtain with other methods. Benefiting from the deployment of dense seismic arrays, progresses have been made in improving the spatial resolution of surface wave imaging, which makes it possible to resolve the fine-scale velocity structures of the crust and upper mantle based on surface wave analysis. Meanwhile, the differences in the S-wave velocities derived from Rayleigh and Love wave data can provide information on the radial anisotropy beneath the seismic arrays. In this thesis, using the NCISP-III broadband data and based on phase velocity dispersion analysis and inversion of fundamental mode Rayleigh and Love waves, I investigated the lateral variations in the S-wave velocity structure of the crust and uppermost mantle beneath the Yanshan Belt and adjacent regions at the northeastern boundary of the NCC. Based on the constructed structural images, I discussed possible deep processes of the craton destruction in the study region.

地震层析成像与天山地区的壳幔深部结构研究

The Tien Shan is the most prominent intracontinental mountain belt on the earth. The active crustal deformation and earthquake activities provide an excellent place to study the continental geodynamics of intracontinental mountain belt. The studies of deep structures in crust and upper mantle are significantly meaningful for understanding the geological evolution and geodynamics of global intracontinental mountain belts. This dissertation focuses on the deep structures and geodynamics in the crust and upper mantle in the Tien Shan mountain belt. With the arrival time data from permanent and temporal seismic stations located in the western and central Tien Shan, using seismic travel time tomographic method, we inversed the P-wave velocity and Vp/Vs structures in the crust and uppermost mantle, the Pn and Sn velocities and Pn anisotropic structures in the uppermost mantle, and the P-wave velocity structures in the crust and mantle deep to 690km depth beneath the Tien Shan. The tomographic results suggest that the deep structures and geodynamics have significant impacts not only on the deformations and earthquake activities in the crust, but also on the mountain building, collision, and dynamics of the whole Tien Shan mountain belt. With the strongly collision and deformations in the crust, the 3-D P-wave velocity and Vp/Vs ratio structures are highly complex. The Pn and Sn velocities in the uppermost mantle beneath the Tien Shan, specially beneath the central Tien Shan, are significantly lower than the seismic wavespeed beneath geological stable regions. We infer that the hot upper mantle from the small-scale convection could elevate the temperature in the lower crust and uppermost mantle, and partially melt the materials in the lower crust. The observations of low P-wave and S-wave velocities, high Vp/Vs ratios near the Moho and the absences of earthquake activities in the lower crust are consistent with this inference. Based on teleseismic tomography images of the upper mantle beneath the Tien Shan, we infer that the lithosphere beneath the Tarim basin has subducted under the Tien Shan to depths as great as 500 km. The lithosphere beneath the Kazakh shield may have subducted to similar depths in the opposite direction, but the limited resolution of this data set makes this inference less certain. These images support the plate boundary model of converge for the Tien Shan, as the lithospheres to the north and south of the range both appear to behave as plates.

拉萨地块中新生代火山作用-从大洋到大陆俯冲转换的构造演化纪录

Late Mesozoic－Cenozoic volcanic rocks are well exposed in Lhasa Terrane, southern Tibet. This research attempts to apply 40Ar/39Ar geochronology, major, trace element and Sr-Nd-O isotopic geochemistry data to constrain the spatio-temporal variations, the composition of source, geodynamic setting. The results indicate that Lhasa Terrane mainly went through three tectonic-magmatic cycle: (1) Phase of Oceanic subduction (140-80Ma). Along with the subducting beneath the Eurasian Plate of Neo-Tethys slab, the oceanic sediment and/or the subducting slab released fluids/melts to metasomatize the subcontinental lithospheric mantle, and induced the mantle wedge partially melt and produced the calc-alkaline continental arc volcanic rocks; (2) Phase of continental-continental collision. Following the subducting of the Neo-Tethys slab, the Indian Plate collided with the Eurasian Plate dragged by the dense Neo-Tethys oceanic lithosphere. The oceanic lithosphere detached from continental lithosphere during roll-back and break-off and the asthenosphere upwelled. The resulting conducted thermal perturbation leads to the melting of the overriding mantle lithosphere and produced the syn-collisional magmatism: the Linzizong Formation and dykes; (3) Following by the detachment of the Tethys oceanic lithosphere, the Indian Lithosphere subducted northward by the drive from the expanding of Indian Ocean. The dense Indian continental lithospheric mantle (±the thickened lower crust) break off, disturb the asthenosphere, and lead to the melting of the overriding mantle lithosphere, which has been metasomatized by the melts/fluids from the subducting oceanic/continental lithosphere and the asthenosphere, and produced the rift-related ultrapotassic rocks.

中国华北和美国加州地区地震层析成像研究

In recent years seismic tomography has become a powerful tool for studying the three-dimensional crust and mantle structure. In this study, we collected a large number of regional and teleseismic travel-time data and used seismic tomography method to study the relationship between earthquake occurrence and crustal heterogeneity for the 1992 Landers earthquake, heterogeneity and evolution of lithosphere under North China Craton and Southern California, and deep structure and origin of the Changbai intraplate volcano in Northeast China. Our results show a correlation between the seismic rupture zone and crustal heterogeneity. The distribution of the Landers aftershocks is cluster-like and separated or terminated in areas where low-velocity anomalies exist．Most of the large earthquakes with magnitudes >4.0 occurred in or around areas with high P-wave velocity．The possibility is that high-velocity areas are brittle and strong parts which can sustain seismogenic stress，and so can generate earthquakes. Our tomographic images show a very heterogeneous structure in the crust and upper mantle beneath Southern California. Three major anomalies in the upper mantle are revealed clearly beneath the southern Sierra Nevada, Transverse Ranges and Salton Trough. We consider that the high-velocity anomaly beneath the Transverse Ranges was formed through asymmetrical two-side convergence of subcrustal lithosphere and sinking to asthenosphere. Formation of the dense crust root and “drip structure” caused the high-velocity anomaly under the southern Sierra Nevada. The Salton Trough low is the response to the lithospheric extension when the Pacific plate was rifted away from the North American Plate. The tomograpic images beneath the North China Craton show that there exist different lithospheric structures under the different blocks. Complex, prominent low-velocity and high-velocity anomalies are imaged beneath the North China Basin, Trans-North China Orogen (TNCO), and Ordos Block which correspond to rifted, orogenic and cratonic lithospheres, respectively. The thickness of the three-type lithospheres is about 70, 90 and >250 km, respectively. Our results suggest that lithospheric thinning under the eastern part of North China Craton is due to long-term replacement and chemical and thermal erosion of the ancient lithosphere by the hot asthenosphere. The remains of ancient lithosphere exist either in the present upper mantle or have sunk into the mantle transition zone. Our tomographic result of the Changbai volcanic area suggests that the origin of the Changbai volcano is related to the deep dehydration of the subducted Pacific slab and corner flow in the big mantle wedge (BMW) above the stagnant Pacific slab.

喜马拉雅西构造结区岩石圈演化三维有限元动力学数值模拟

The past three decades have seen numerous attempts to numerically model stress and strain patterns in the lithosphere of the Earth on both global and regional scales. This efforts have been indispensable in identifying the features we need to include in our endeavour to develop better models of our planet’s lithosphere and they have also raised our awareness for the many unresolved issue in the deep geodynamical issues that need to be addressed in the future. Nonetheless, in most models, the lithosphere is treated as a single layer with depth-averaged properties, and as the same distribution in the stress and strain fields, and as deforming under plane strain. All these above make a great hander for its reality and degree of recognition. As the beginning in this paper, some principal numerical models and results on the evolution of Tibetan plateau are reviewed and analyzed. Then, the geological and geophysical expedition on the Western Himalayan Syntaxis is briefly reviewed. Furthermore, we analysis the feature in deep geophysical field studies in this area and adjacent regions. Because, for most continents, stress models driven by plate boundary forces have successfully reproduced the main characteristics of the stress and strain field, we present a set of three-dimensional models of lithosphere system for a simplified geometry of the Western Himalayan Syntaxis area and its adjacent regions, where we try to match the first-order characteristics of the stress and strain fields of lithosphere since 10 Ma, and deformation and geodynamical evolution process in former 2Ma. Of course, the kinematic boundary conditions of the stress models driven by plate boundary forces were applied. The rheology plays a significant role in the lithospheric tectonics, which lead to different rheological parameters were used in different works although the have the same constitutive equations in models. So, in this paper we do not aim to produce all characteristics of the Western Himalayan Syntaxis areas’ stress and strain fields by the choices of various parameters, but rather the dynamic response between various rheological parameters and stress and strain fields. We have chosen to concentrate on the importance of rheology and lateral strength variations for lithospheric stress and strain patterns and use our findings to build a model of the Western Himalayan Syntaxis areas. In doing so, we want to go beyond purely elastic models or purely viscoelastic models. Compared the results of the crust viscosity in the Western Himalayan Syntaxis areas, we believed that, when various viscoelastic models are adopted, the selection of the coefficient of viscosity in the Western Syntaxis area has important influence on the its uplifts and evolutions. A wider uplift ranges and gently elevation was observed at the same time when a lower viscosity was used in our models, and vice versa. Data of stress magnitudes are not available, but it is clear that the stress levels must be at or below the failure threshold of rock under compression. Under these criteria, the calculation results show that the viscosity in the Western Syntaxis area should be smaller than 1023Pa.s When elastic model is adopted in relatively rigid Tarim basin, obvious changes are induced to the stress and strain fields of the whole Western Syntaxis area. We found that rigid block of lithosphere reduced stress levels within its interior and that, at the edges of such regions, stress orientation can change. Furthermore there is no evidence that such rigid regions act as stress barriers in that they shield areas in opposite sides of the structure from the influence of one another. In our models, the upper crustal material of the Western Syntaxis area does not turns to move westward. Whereas, because of the stress and strain fields have been decoupling at the interior of the lithosphere, we can get the results that the deep material must not move westward.