Numerical simulations of CO2 migration during charnockite genesis

Charnockite is considered to be generated either through the dehydration of granitic magma by CO2 purging or by solid-state dehydration through CO2 metasomatism during granulite facies metamorphism. To understand the extent of dehydration, CO2 migration is quantitatively modeled in silicate melt and metasomatic fluid as a function of temperature, H2O wt%, pressure, basal CO2 flux and dynamic viscosity. Numerical simulations show that CO2 advection through porous and permeable high-grade metamorphic rocks can generate dehydrated patches close to the CO2 flow path, as illustrated by the occurrences of ``incipient charnockites.'' CO2 reaction-front velocity constrained by field observations is 0.69 km/m.y., a reasonable value, which matches well with other studies. On the other hand, temperature, rate of cooling, and basal CO2 flux are the critical parameters affecting CO2 diffusion through a silicate melt. CO2 diffusion through silicate melt can only occur at temperature greater than 840 degrees C and during slow cooling (<= 3.7 x 10(-5) degrees C/yr), features that are typical of magma emplacement in the lower crust. Stalling of CO2 fluxing at similar to 840 degrees C explains why some deep-level plutons contain both hydrous and anhydrous (charnockitic) mineral assemblages. CO2 diffusion through silicate melt is virtually insensitive to pressure. Addition of CO2 basal flux facilitates episodic dehydrated melt migration by generating fracture pathways.

Neoarchean continental growth through arc magmatism in the Nilgiri Block, southern India

The formation and growth of continental crust in the Archean have been evaluated through models of subduction-accretion and mantle plume. The Nilgiri Block in southern India exposes exhumed Neoarchean lower crust, uplifted to heights of 2500 m above sea level along the north western margin of the Peninsula. Major lithologies in this block include charnockite with or without garnet, anorthosite-gabbro suite, pyroxenite, amphibolite and hornblende-biotite gneiss (TTG). All these rock types are closely associated as an arc magmatic suite, with diffuse boundaries and coeval nature. The charnockite and hornblende-biotite gneisses (TTG) show SiO2 content varying from 64 to 73 wt.%. The hornblende-biotite gneisses (TTG) are high-Al type with Al2O3 >15 wt.% whereas the charnockites show Al2O3 <15 wt.%. The composition of charnockite is mainly magnesian and calcic to calc-alkaline. The mafic-ultramafic rocks show composition close to that of tholeiitic series. The low values of K(2)o (<3 wt.%), (K/Rb)/K2O (<500), Zr/Ti, and trace element ratios like (La/Yb)n/(Sr/Y), (Y/Nb), (Y + Nb)/Rb, (Y+Ta)/Rb, Yb/Ta indicate a volcanic arc signature for these rocks. The geochemical signature is consistent with arc magmatic rocks generated through oceanic plate subduction. The primitive mantle normalized trace element patterns of these rocks display enrichment in large ion lithophile elements (LILE) and comparable high field strength elements (HFSE) in charnockite and hornblende-biotite gneisses (TTG) consistent with subduction-related origin. Primitive mantle normalized REE pattern displays an enrichment in LREE in the chamockite and hornblende-biotite gneisses (TTG) as compared to a flat pattern for the mafic rocks. The chondrite normalized REE patterns of zircons of all the rock types reveal cores with high HREE formed at ca. 2700 Ma and rims with low HREE formed at 2500-2450 Ma. Log-transformed La/Th-Nb/Th-Sm/Th-Yb/Th discrimination diagram for the mafic and ultramafic rocks from Nilgiri displays a transition from mid-oceanic ridge basalt (MORB) to island arc basalt (IAB) suggesting a MORB source. The U-Pb zircon data from the charnockites, mafic granulites and hornblende-biotite gneisses (TTG) presented in our study show that the magma generation during subduction and accretion events in this block occurred at 2700-2500 Ma. Together with the recent report on Neoarchean supra-subduction zone ophiolite suite at its southern margin, the Nilgiri Block provides one of the best examples for continental growth through vertical stacking and lateral accretion in a subduction environment during the Neoarchean. (c) 2014 Elsevier B.V. All rights reserved.

Neoarchean arc magmatism followed by high-temperature, high-pressure metamorphism in the Nilgiri Block, southern India

The Nilgiri Block, southern India is an exhumed lower crust formed through arc magmatic processes in the Neoarchean. The main lithologies in this terrane include charnockites, gneisses, volcanic tuff, metasediments, banded iron formation and mafic-ultramafic bodies. Mafic-ultramafic rocks are present towards the northern and central part of the Nilgiri Block. We examine the evolution of these mafic granulites/metagabbros by phase diagram modeling and U-Pb sensitive high resolution ion microprobe (SHRIMP) dating. They consist of a garnet-clinopyroxene-plagioclase-hornblende-ilmenite +/- orthopyroxene +/- rutile assemblage. Garnet and clinopyroxene form major constituents with labradorite and orthopyroxene as the main mineral inclusions. Labradorite, identified using Raman analysis, shows typical peaks at 508 cm(-1), 479 cm(-1), 287 cm(-1) and 177 cm(-1). It is stable along with orthopyroxene towards the low-pressure high-temperature region of the granulite fades (M1 stage). Subsequently, orthopyroxene reacted with plagioclase to form the peak garnet + clinopyroxene + rutile assemblage (M2 stage). The final stage is represented by amphibolite facies-hornblende and plagioclase-rim around the garnet-clinopyroxene assemblage (M3 stage). Phase diagram modeling shows that these mafic granulites followed an anticlockwise P-T-t path during their evolution. The initial high-temperature metamorphism (M1 stage) was at 850-900 degrees C and similar to 9 kbar followed by high-pressure granulite fades metamorphism (M2 stage) at 850-900 degrees C and 14-15 kbar. U-Pb isotope studies of zircons using SHRIMP revealed late Neoarchean to early paleoproterozoic ages of crystallization and metamorphism respectively. The age data shows that these mafic granulites have undergone arc magmatism at ca. 25392 +/- 3 Ma and high-temperature, high-pressure metamorphism at ca. 2458.9 +/- 8.6 Ma. Thus our results suggests a late Neoarchean arc magmatism followed by early paleoproterozoic high-temperature, high-pressure granulite facies metamorphism due to the crustal thickening and suturing of the Nilgiri Block onto the Dharwar Craton. (C) 2015 Elsevier B.V. All rights reserved.

Arquitetura crustal da bacia de Almada no contexto das bacias da margem lesste da América do Sul

A Bacia de Almada, localizada no estado da Bahia, compartilha características similares com as outras bacias da margem leste do Brasil, quando é analisada segundo aspectos como os processos sedimentares e o regime de esforço dominante durante a sua formação. Observa-se uma diferença marcante em relação as outras bacias quando é analisada sob a ótica da composição da crosta transicional, uma vez que não se registra atividade vulcânica durante a fase rifte. A aquisição de um extenso levantamento sísmico 3D, com cabos de 6 km de comprimento e 9.2 segundos de tempo de registro (tempo sísmico duplo), resultaram em imagens sísmicas de boa qualidade das estruturas profundas do rifte. Adicionalmente, estudos de modelagem gravimétrica foram integrados com a análise sísmica para corroborar o modelo geológico. A Bacia de Almada é parte dos sistemas de rifte continentais, desenvolvidos durante o Berriasiano até o Aptiano, que antecederam a quebra do continente do Gondwana, evoluindo posteriormente para uma margem passiva divergente. O processo do rifteamento desenvolveu cinco sub-bacias de orientação NNE-SSO, desde posições terrestres até marinhas profundas, produzindo um arcabouço estrutural complexo. Os perfis da sísmica profunda mostram o afinamento progressivo da crosta continental até espessuras da ordem de 5 km, abaixo da sub-bacia mais oriental, com fatores de estiramento crustal próximo a 7 antes do desenvolvimento de crosta oceânica propriamente dita. As imagens sísmicas de boa qualidade permitem também o reconhecimento de sistemas de falhas lístricas que se iniciam na crosta superior, evoluem atravessando a crosta e conectando as sub-bacias para finalizar em um descolamento horizontal na crosta inferior estratificada. Adicionalmente, a bacia apresenta um perfil assimétrico, compatível com mecanismos de cisalhamento simples. As margens vulcânicas (VM) e não vulcânicas (NVM), são os extremos da análise composicional das margens divergentes continentais. Na Bacia de Almada não se reconhecem os elementos arquiteturais típicos das VM, tais como são as grandes províncias ígneas, caracterizadas por cunhas de refletores que mergulham em direção ao mar e por intenso vulcanismo pré- e sin-rifte nas bacias. Embora a margem divergente do Atlântico Sul seja interpretada tradicionalmente como vulcânica, o segmento do rifte ao sul do Estado da Bahia apresenta características não-vulcânicas, devido à ausência destes elementos arquiteturais e aos resultados obtidos nas perfurações geológicas que eventualmente alcançam a seqüência rifte e embasamento. Regionalmente a margem divergente sul-americana é majoritariamente vulcânica, embora a abundância e a influência do magmatísmo contemporâneo ao rifte seja muito variável. Ao longo da margem continental, desde a Bacia Austral no sul da Argentina, até a Bacia de Pernambuco no nordeste do Brasil, podem ser reconhecidos segmentos de caráter vulcânico forte, médio e não vulcânico. Nos exemplos clássicos de margens não vulcânicas, como a margem da Ibéria, a crosta transicional é altamente afinada podendo apresentar evidências de exumação de manto. Na Bacia de Almada, a crosta transicional apresenta importante estiramento embora não haja evidências concretas de exumação de manto. Os mecanismos responsáveis pela geração e intrusão dos grandes volumes de magma registrados nas margens divergentes são ainda sujeitos a intenso debate. Ao longo da margem divergente sul-americana há evidências da presença dos mecanismos genéticos de estiramento litosférico e impacto de plumas. Alternativamente estes dois mecanismos parecem ter tido um papel importante na evolução tectônica da margem sudeste e sul, diferenciando-as da margem continental onde foi implantada a Bacia de Almada.

High frequency wave propagation in the earth: theory and observation

The wave-theoretical analysis of acoustic and elastic waves refracted by a spherical boundary across which both velocity and density increase abruptly and thence either increase or decrease continuously with depth is formulated in terms of the general problem of waves generated at a steady point source and scattered by a radially heterogeneous spherical body. A displacement potential representation is used for the elastic problem that results in high frequency decoupling of P-SV motion in a spherically symmetric, radially heterogeneous medium. Through the application of an earth-flattening transformation on the radial solution and the Watson transform on the sum over eigenfunctions, the solution to the spherical problem for high frequencies is expressed as a Weyl integral for the corresponding half-space problem in which the effect of boundary curvature maps into an effective positive velocity gradient. The results of both analytical and numerical evaluation of this integral can be summarized as follows for body waves in the crust and upper mantle:

1) In the special case of a critical velocity gradient (a gradient equal and opposite to the effective curvature gradient), the critically refracted wave reduces to the classical head wave for flat, homogeneous layers.

2) For gradients more negative than critical, the amplitude of the critically refracted wave decays more rapidly with distance than the classical head wave.

3) For positive, null, and gradients less negative than critical, the amplitude of the critically refracted wave decays less rapidly with distance than the classical head wave, and at sufficiently large distances, the refracted wave can be adequately described in terms of ray-theoretical diving waves. At intermediate distances from the critical point, the spectral amplitude of the refracted wave is scalloped due to multiple diving wave interference.

These theoretical results applied to published amplitude data for P-waves refracted by the major crustal and upper mantle horizons (the Pg, P*, and Pn travel-time branches) suggest that the 'granitic' upper crust, the 'basaltic' lower crust, and the mantle lid all have negative or near-critical velocity gradients in the tectonically active western United States. On the other hand, the corresponding horizons in the stable eastern United States appear to have null or slightly positive velocity gradients. The distribution of negative and positive velocity gradients correlates closely with high heat flow in tectonic regions and normal heat flow in stable regions. The velocity gradients inferred from the amplitude data are generally consistent with those inferred from ultrasonic measurements of the effects of temperature and pressure on crustal and mantle rocks and probable geothermal gradients. A notable exception is the strong positive velocity gradient in the mantle lid beneath the eastern United States (2 x 10^-3 sec^-1), which appears to require a compositional gradient to counter the effect of even a small geothermal gradient.

New seismic-refraction data were recorded along a 800 km profile extending due south from the Canadian border across the Columbia Plateau into eastern Oregon. The source for the seismic waves was a series of 20 high-energy chemical explosions detonated by the Canadian government in Greenbush Lake, British Columbia. The first arrivals recorded along this profile are on the Pn travel-time branch. In northern Washington and central Oregon their travel time is described by T = Δ/8.0 + 7.7 sec, but in the Columbia Plateau the Pn arrivals are as much as 0.9 sec early with respect to this line. An interpretation of these Pn arrivals together with later crustal arrivals suggest that the crust under the Columbia Plateau is thinner by about 10 km and has a higher average P-wave velocity than the 35-km-thick, 62-km/sec crust under the granitic-metamorphic terrain of northern Washington. A tentative interpretation of later arrivals recorded beyond 500 km from the shots suggests that a thin 8.4-km/sec horizon may be present in the upper mantle beneath the Columbia Plateau and that this horizon may form the lid to a pronounced low-velocity zone extending to a depth of about 140 km.

Caracterização do maciço Santa Clara no município de Cujubim (RO) com base em litogeoquímica, geocronologia e estudos isotópicos

A Suíte Intrusiva Santa Clara está inserida na Província Estanífera de Rondônia, na porção SW do Cráton Amazônico. Essa suíte intrusiva é composta pelos maciços Santa Clara, Oriente Velho, Oriente Novo, Manteiga-Sul, Manteiga-Norte, Jararaca, Carmelo, Primavera e das Antas. Os litotipos que perfazem a Suíte Santa Clara ocorrem hospedados nas rochas do Complexo Jamari, uma associação polideformada composta por gnaisses ortoderivados e paraderivados. Características observadas em campo e em análises petrográficas permitiram subdividir o Maciço Santa Clara em cinco fácies distintas: fácies porfirítica, fácies isotrópica, fácies fina, fácies piterlítica e fácies viborgítica. Os litotipos observados correspondem a hornblenda-biotita granitos e biotita granitos intermediários a ácidos, com composições médias semelhantes àquelas verificadas para sienogranitos e monzogranitos. Geoquimicamente, três magmas podem ser identificados. O magma menos evoluído corresponde às rochas das fácies porfirítica e equigranular, e o mais evoluído compreende as fácies de granulometria fina e piterlítica. A fácies viborgítica representa o terceiro líquido magmático, e aparentemente é diferente de todas as outras fácies em termos de aspectos de campo e geoquímica. A análise litogeoquímica indica que estes granitoides são subalcalinos, bastante empobrecidos em MgO e exibem caráter metaluminoso a fracamente peraluminoso. Os padrões de elementos-traços evidenciam que tais granitóides possuem alto conteúdo em elementos incompatíveis (Rb, Zr, Y, Ta, Ce) e ETR, com exceção do Eu. Além disso, também exibem leve enriquecimento em LILE, forte depleção em elementos como Sr e Ti, e leve empobrecimento de Ba, indicando que o fracionamento de minerais como plagioclásio e titanita foi importante na evolução do líquido magmático analisado. A anomalia negativa de Nb indica envolvimento de material crustal nos processos magmáticos que geraram estes granitoides. Os litotipos analisados possuem características típicas de granitos tipo-A ferroan, e as razões FeOt/MgO entre 4,27 e 26,22 sugerem tratar-se de uma série de granitos félsicos fracionados. Os padrões de ETR observados para os litotipos analisados exibem um considerável enriquecimento em ETRL, e anomalia negativa de Eu, sugerindo fracionamento de feldspato durante o processo de diferenciação do líquido magmático. Diagramas discriminantes de ambientes tectônicos sugerem que os litotipos do Maciço Intrusivo Santa Clara são típicos de ambiente intraplaca, do tipo-A2, isto é, associados a ambientes pós-colisionais/pós-orogênicos. As características isotópicas observadas para os granitoides do Maciço Santa Clara sugerem que os mesmos foram gerados a partir da fusão parcial de uma crosta inferior pré-existente. As idades U-Pb entre 1,07 e 1,06 Ga são compatíveis com um magmatismo ocorrido nos estágios finais da colagem do supercontinente Rodínia (1,2-1,0 Ga) e estágios finais do Ciclo Orogênico Sunsás-Aguapeí (1320-1100 Ma). Sugere-se ainda que na verdade o Maciço Santa Clara seja formado por uma coalescência das três intrusões graníticas que são representadas pelos três magmas anteriormente descritos.

Comparação entre o Complexo Juiz de fora e a Unidade Granulítica Ponte de Zinco: geocronologia U-Pb em zircão (LA-ICPMS), geoquímica isotópica e composição das fontes geradoras

A comparação entre o Complexo Juiz de Fora e a Unidade Granulítica Ponte de Zinco (Mangaratiba - RJ) revelou que existem diferenças significativas entre essas unidades. Na Unidade Granulítica Ponte de Zinco, são encontrados ortogranulitos de composição granítica e granodiorítica, que representam o embasamento da Unidade, denominada aqui de Ortogranulitos Ribeirão das Lajes. Dois outros litotipos ortoderivados também foram encontrados: (i) um ortognaisse com granada (MAN-JEF-03a), sendo que a granada ocorre de forma subordinada; (ii) ortognaisse leucocrático (MAN-JEF-04), com características miloníticas. Além dessas, rochas metassedimentares também afloram na Unidade Granulítica Ponte de Zinco. Foi interpretado que um ortogranulito (MAN-JEF-01a) de alto-K, com composição monzogranítica, cristalizou em ca. 2653 37 Ma (U-Pb em zircão por LA-ICPMS), afetado por um evento no Neoproterozóico, que gerou os minerais máficos hidratados, observados na análise petrográfica e como mostram as borda de sobrecrescimento em zircão. Sua idade modelo de Nd de 2,7 Ga e seu εNd positivo de +2,1, apontam para uma gênese mantélica, tendo assimilado rochas crustais, pois são encontrados grãos de zircão herdados de aproximadamente 2996 17 Ma e 3343 3.8 Ma. Os dados de litogeoquímica e sua razão 87Sr/86 Sr(t) (0,70529), são compatíveis com uma geração em um arco continental. O ortognaisses com granada do ponto MAN-JEF-03 possui composição granodiorítica. A idade de cristalização interpretada pela análise geocronológica U-Pb em zircão (LA-ICPMS), foi ca. 2117 15 Ma. Esse litotipo foi metamorfizado no Neoproterozóico, sendo a idade obtida pelo intercepto inferior de 631 40 Ma. Seus dados isotópicos apontam para uma rocha juvenil gerada a partir do manto (TDM ≈ de 2,1 Ga e εNd = +3,4). Sua alta razão 87Sr/86 Sr(t) (0,710 ) juntamente com os grãos de zircão herdados (2,6 Ga) e a presença de enclaves, indicam assimilação de rochas crustais. O ortognaisse leucocrático (MAN-JEF-04) classificado como alto-K, possui composição monzogranítica, idade 2132 9,4 Ma U-Pb em zircão (LA-ICPMS). Um único grão relíquiar de ortopiroxênio encontrado em lâmina, , indica que a rocha já foi submetida a metamorfismo de fácies granulito, porém esse evento não deixou registro nos grãos analisados. O retrometamorfismo pode ter ocorrido em dois momentos, 647 11 Ma e 595 38 Ma, calculados através da concordia age, em sobrecrescimentos homogêneos e, pelo intercepto inferior, respectivamente. Sua baixa razão 87Sr/86 Sr(t) (≈ 0,703) associada com εNd positivo (+2,3) e sua idade modelo de aproximadamente 2,1 Ga, revelam que a rocha foi formada por um material mantélico juvenil. Já as análises geocronológicas em U-Pb em zircão (LA-ICPMS) na região de Juiz de Fora (MG), revelaram a existência de dois litotipos Arqueanos: um ortogranulito granodiorítico (MB-JEF-01b), de baixo-K com idade de 2849 11 Ma e com herança de 2975 10 Ma. Seu εNd positivo (+5,9) aponta para uma gênese a partir do manto depletado, já sua alta razão 87Sr/86Sr(t) (≈0,709) indica contaminação de Rb de fontes externas, talvez causada pela assimilação da crosta, como revelam os zircões herdados e/ou fluidos retrometamorficos. Outro litotipo é uma rocha gabróica do tipo E-MORB, cuja idade foi calculada em 2691 14 Ma, com retrometamorfismo ocorrido no intervalo de 604 67 Ma, obtida pelo intercepto inferior. Seu εNd igual a +3,4 e sua razão 87Sr/86 Sr(t) (≈0,701) mostram extração a partir do manto depletado. Novos dados isotópicos do CJF na região de Três Rios (RJ) e Juiz de Fora (MG), sugerem que os ortogranulitos calcioalcalinos podem representar grupos distintos. Rochas com εNd positivos são consequentemente associadas ao manto depletado, porém rochas com εNd negativos devem ter sido geradas por fusão crustal, que podem ser fusão de crosta inferior, devido a razão 87Sr/86 Sr(t) (0,70514) encontrada no amostra MB-JEF-02a (ortogranulito de alto-K) ou tendo a crosta contribuição nas gênese dessas rochas. Os ortogranulitos básicos possuem εNd positivos com baixas razões 87Sr/86 Sr(t) , o que indica extração a partir do manto depletado, porém sua razões La/YbN e La/NbN maiores que 1, revelam alguma contribuição de uma fonte enriquecida, assim também mostram suas razões Pb/Pb, que são maiores do que as razões calculadas para evolução de Pb na Terra. Essas interpretações ainda podem ser estendidas para um ortoanfibolito da série alcalina, encontrado na região de Três Rios (RJ).

Deep Fractionation of Clinopyroxene in the East Pacific Rise 13 degrees N: Evidence from High MgO MORB and Melt Inclusions

Mid-ocean ridge basalts (MORBs) from East Pacific Rise (EPR) 13 degrees N are analysed for major and trace elements, both of which show a continuous evolving trend. Positive MgO-Al2O3 and negative MgO-Sc relationships manifest the cotectic crystallization of plagioclase and olivine, which exist with the presence of plagioclase and olivine phenocrysts and the absence of clinopyroxene phenocrysts. However, the fractionation of clinopyroxene is proven by the positive correlation of MgO and CaO. Thus, MORB samples are believed to show a "clinopyroxene paradox". The highest magnesium-bearing MORB sample E13-3B (MgO = 9.52%) is modelled for isobaric crystallization with COMAGMAT at different pressures. Observed CaO/Al2O3 ratios can be derived from E13-3B only by fractional crystallization at pressure > 4 +/- 1 kbar, which necessitates clinopyroxene crystallization and is not consistent with cotectic crystallization of olivine plus plagioclase in the magma chamber (at pressure similar to 1 kbar). The initial compositions of the melt inclusions, which could represent potential parental magmas, are reconstructed by correcting for post-entrapment crystallization (PEC). The simulated crystallization of initial melt inclusions also produce observed CaO/Al2O3 ratios only at > 4 +/- 1 kbar, in which clinopyroxene takes part in crystallization. It is suggested that MORB magmas have experienced clinopyroxene fractionation in the lower crust, in and below the Moho transition zone. The MORB magmas have experienced transition from clinopyroxene+plagioclase+olivine crystallization at > 4 +/- 1 kbar to mainly olivine+plagioclase crystallization at < 1 kbar, which contributes to the explanation of the "clinopyroxene paradox".

Analysis on deep crustal structure along the onshore-offshore seismic profile across the Binhai (littoral) fault zone in northeastern South China Sea

The onshore-offshore deep seismic experiment was carried out for the first time and filled the blankness of the seismic surveys in the transition area between South China and northeastern South China Sea. The seismic data were analyzed and processed. The different seismic phases were identified and their travel time arrivals were modeled by ray-tracing to study the P-wave velocity crustal structure of this area. The crustal structure of this area is the continental crust. The crust thickness is gradually decreasing southward along the on-shore-offshore seismic line. The low-velocity layer (5.5 similar to 5.9 km (.) s(-1)) exists generally in the middle crust (about 10.0 similar to 18.0km)with about 2.5 similar to 4.0 km thickness, which is also thinning seaward. No obvious high-velocity layer appears in the lower crust. The Binhai (littoral) fault zone is a low velocity zone, which is located about 35km southeast to the Nan'ao station and corresponding to the gradient belt of gravity & magnetism anomalies. The depth of the fault zone is close to the Moho discontinuity. The littoral fault zone is a boundary between the normal continental crust of South China and the thinned continental crust of the sea area.

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

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.

青藏高原横波速度与各向异性结构及其壳幔形变印迹

As the most spectacular and youngest case of continental collision on the Earth, to investigate the crust and mantle of Tibetan plateau, and then to reveal its characters of structure and deformation, are most important to understand its deformation mechanism and deep process. A great number of surface wave data were initially collected from events occurred between 1980 and 2002, which were recorded by 13 broadband digital stations in Eurasia and India. Up to 1,525 source-station Rayleigh waveforms and 1,464 Love wave trains were analysed to obtain group velocity dispersions, accompanying with the detail and quantitative assessment of the fitness of the classic Ray Theory, errors from focal and measurements. Assuming the model region covered by a mesh of 2ox2o-sized grid-cells, we have used the damped least-squares approach and the SVD to carry out tomographic inversion, SV- and SH-wave velocity images of the crust and upper mantle beneath the Tibetan Plateau and surroundings are obtained, and then the radial anisotropy is computed from the Love-Rayleigh discrepancy. The main results demonstrate that follows, a) The Moho beneath the Tibetan Plateau presents an undulating shape that lies between 65 and 74 km, and a clear correlation between the elevations of the plateau and the Moho topography suggests that at least a great part of the highly raised plateau is isostatically compensated. b) The lithospheric root presents a depth that can be substantiated at ~140 km (Qiangtang Block) and exceptionally at ~180 km (Lhasa Block), and exhibits laterally varying fast velocity between 4.6 and 4.7 km/s, even ~4.8 km/s under northern Lhasa Block and Qiangtang Block, which may be correlated with the presence of a shield-like upper mantle beneath the Tibetan Plateau and therefore looked as one of the geophysical tests confirming the underthrusting of India, whose leading edge might have exceeded the Bangong-Nujiang Suture, even the Jinsha Suture. c) The asthenosphere is depicted by a low velocity channel at depths between 140 and 220 km with negative velocity gradient and velocities as low as 4.2 km/s; d) Areas in which transverse radial anisotropy is in excess of ~4% and 6% on the average anisotropy are found in the crust and upper mantle underlying most of the Plateau, and up to 8% in some places. The strength, spatial configuration and sign of radial anisotropy seem to indicate the existence of a regime of horizontal compressive forces in the frame of the convergent orogen at the same time that laterally varying lithospheric rheology and a differential movement as regards the compressive driving forces. e) Slow-velocity anomalies of 12% or more in southern Tibet and the eastern edge of the Plateau support the idea of a mechanically weak middle-to-lower crust and the existence of crustal flow in Tibet.

西天山造山带古生代构造演化与地壳增生—来自花岗岩和蛇绿岩的证据

The Central Asian Orogen Belt (CAOB), which is different from the subductional orogen and the collisional orogen, is known as the most important site of crustal growth in the Phanerozoic, and it has been a ‘hot spot’ for studying the orogenic belts. The Chinese West Tianshan Orogen is occupying the west-southern part of the CAOB and is of great importances to understand the orogenic processes and the continental growth in the Central Asia. The West Tianshan Orogen had undergone complex tectonic evolutional processes in Paleozoic times and large volumes granitic rocks have recorded important information about these processes. Litter is known about Phanerozoic continental growth in the Western Tianshan area so far, compared with the other areas of the CAOB, such as eastern Junggar, western Junggar, Altai and Alakol. The aim of this dissertation is to set up the chronology frame of granitoids in western Tianshan, provide new evidence for the tectonic evolution and discuss the Paleozoic continental growth in this area, on the basis of the studies on the isotopic chronology, major element, trace element and Nd-Sr isotopic geochemistry of granitoids and the isotopic chronology and geochemistry of the ophiolites in this area, especially the Kule Lake ophiolites. 25 precise SHRIMP U-Pb zircon and LA-ICPMS U-Pb zircon ages have been obtained in this dissertation. The granitic rocks in western Tianshan had been formed during two periods: the granitic gneiss with an age of 896Ma, possibly representing the forming age of the Precambrian basement; the granitic rocks with ages varying from 479Ma to 247Ma, recording the Paleozoic orogenic process of western Tianshan. The granitoids in western Tianshan are composed of intermediate-basic rocks, intermediate rocks, intermediate-acid rocks and acid rocks, mainly intermediate-acid rocks and acid rocks. They are mostly granite, granodiorite, quartz syenite and monzodiorite. Different types of granitic rocks are exposed in different tectonic units. The granitoids on the northern margin of the Yili Plate mainly formed in late Paleozoic (413Ma ~ 281Ma), those with ages varying from 413Ma to 297Ma show continental arc affinities and the magnesian calc-alkalic metaluminous diorite of 281Ma display the geochemical characteristics similar to those of granites formed during the post-orogenic period. The granitiods on the southern margin of the Yili Plate include the adakite diorite of 470Ma which was formd by partial melting of thickened lower crust, the post-collisional alkali-feldspar granite of 430Ma, the volcanic arc granite of 348Ma and the Triassic post-collisional granite. The granitoids in the Central Tianshan Plate formed in 479Ma ~ 247Ma, mainly in 433Ma ~ 321Ma. The granitic rocks with ages of 479Ma ~ 321Ma are magnesian calc-alkalic to alkalic rocks with continental arc affinities. A few post-collisional granitoids of 276Ma ~ 247Ma may have inherited the geochemical characteristics of pre-existing arc magma. The granitic rocks in Southern Tianshan (northern margin of the Tarim plate) formed two stages, 420Ma ~ 411Ma and ca. 285Ma. The magnesian calcic to alkalic granites of 420Ma ~ 411Ma may formed during the extension process of the continental margin. The granite of 285Ma includes mostly ferroan calc-alkalic to alkali-calcic rocks with high SiO2 and high alkaline contents, and obviously negative anomaly of Eu, Ba, Sr, P, Ti, similar to the geochemical characteristics of the A-type granite which is formed during post-collisional extension. The Kule Lake ophiolite in southern Tianshan shows the affinity of N-MORB. A SHRIMP zircon U-Pb age of 425±8Ma has obtained for gabbros. Some zircons have given another group of 206Pb/238U age 918Ma, which may indicate the information of the pre-exist old basement rock. The small oceanic basin represented by Kule Lake ophiolite probably developed on the split northern margin of Tarim block. A model for Paleozoic tectonic evolution of the West Tianshan Orogen has been proposed here on the basis of the new results obtained in this dissertation and the previous published data. In Early Cambrian, the Terskey Ocean occurred along the North Nalati fault (NNF), and it separated the Yili plate from the Central Tianshan plate which was probably connected with the Tarim plate. The Terskey Ocean probably subducted towards south under the Central Tianshan plate and towards north under the Yili plate simultaneously. In the early stage of Late Ordovician, the Terskey Ocean had been closed, and the Yili and Central Tianshan plates collided. Meanwhile, extension happened within the joint Central Tianshan and Tarim plates gradually and the Paleo-South Tianshan Ocean had been formed. In Early Silurian, the Paleo-South Tianshan Ocean began to subduct beneath the composite Yili-Central Tianshan plate, which was intruded by volcanic arc granitoids. In Middle Silurian, the Paleo-South Tianshan Ocean, which had reached a certain width, was subducting strongly. And this subduction may have produced voluminous granitoids in the Central Tianshan plate. In the latest stage of Carboniferous, the Paleo-South Tianshan ocean closed, and the Yili-Central Tianshan plate and Tarim plate collided. In Late Cambrian, Paleo-Junggar Ocean occurred to north of the Yili plate; and started to subduct towards south under the Yili plate in Ordovician. This subduction may have produced a magma arc on the northern margin of the Yili plate. In Late Carboniferous, the Paleo-Junggar Ocean had been closed. The Yili-Central and Junggar plates amalgamated together. The West Tianhan Orogen may undergo a post-collisional collapse since Permian. And the magmatic activities may continue to early Triassic. The initial 87Sr/86Sr ration of the granitic rocks in the western Tianshan Mountains varies from 0.703226 to 0.716343, and Nd（t）from -6.50 to 2.03. The characteristics of Sr-Nd isotope indicate that the source of granitic material is not a sole source, which may be produced by mantle-crust magma mixing. In Paleozoic time, lateral growth of the continental crust along active continental margins was dominant, whereas the vertical growth of continental crust resulted from post- collisional mantle derived magmas was not obvious.

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

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.