Assessing the potential of VEGETATION sensor data for mapping snow and ice cover: a Normalized Difference Snow and Ice Index

The VEGETATION (VGT) sensor in SPOT 4 has four spectral bands that are equivalent to Landsat Thematic Mapper (TM) bands (blue, red, near-infrared and mid-infrared spectral bands) and provides daily images of the global land surface at a 1-km spatial resolution. We propose a new index for identifying and mapping of snow ice cover, namely the Normalized Difference Snow/Ice Index (NDSII), which uses reflectance values of red and mid-infrared spectral bands of Landsat TM and VGT. For Landsat TM data, NDSII is calculated as NDSIITM =(TM3 -TM5)/(TM3 +TM5); for VGT data, NDSII is calculated as NDSIIVGT =(B2- MIR)/(B2 + MIR). As a case study we used a Landsat TM image that covers the eastern part of the Qilian mountain range in the Qinghai-Xizang (Tibetan) plateau of China. NDSIITM gave similar estimates of the area and spatial distribution of snow/ice cover to the Normalized Difference Snow Index (NDSI=(TM2-TM5)/(TM2+TM5)) which has been proposed by Hall et al. The results indicated that the VGT sensor might have the potential for operational monitoring and mapping of snow/ice cover from regional to global scales, when using NDSIIVGT.

Molecular systematics of pikas (genus Ochotona) inferred from mitochondrial DNA sequences

The phylogenetic relationships among worldwide species of genus Ochotona were investigated by sequencing mitochondrial cytochrome b and ND4 genes. Parsimony and neighbor-joining analyses of the sequence data yielded congruent results that strongly indicated three major clusters: the shrub-steppe group, the northern group, and the mountain group. The subgeneric classification of Ochotona species needs to be revised because each of the two subgenera in the present classification contains species from the mountain group. To solve this taxonomic problem so that each taxon is monophyletic, i.e., represents a natural clade, Ochotona could be divided into three subgenera, one for the shrub-steppe species, a second for the northern species, and a third for the mountain species. The inferred tree suggests that the differentiation of this genus in the Palearctic Region was closely related to the gradual uplifting of the Tibet (Qinghai-Xizang) Plateau, as hypothesized previously, and that vicariance might have played a major role in the differentiation of this genus on the Plateau, On the other hand, the North American species, O. princeps, is most likely a dispersal event, which might have happened during the Pliocene through the opening of the Bering Strait. The phylogenetic relationships within the shrub-steppe group are worth noting in that instead of a monophyletic shrub-dwelling group, shrub dwellers and steppe dwellers are intermingled with each other. Moreover, the sequence divergence within the sister tars of one steppe? dweller and one shrub dweller is very low. These findings support the hypothesis that pikes have entered the steppe environment several times and that morphological similarities within steppe dwellers were due to convergent evolution. (C) 2000 Academic Press.

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

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.

Pg震相走时的各向异性影响及各向异性结构成像

As the first arrival of seismic phase in deep seismic sounding, Pg is the important data for studying the attributes of the sedimentary layers and the shape of crystalline basement because of its high intensity and reliable detection. Conventionally, the sedimentary cover is expressed as isotropic, linear increasing model in the interpretation of Pg event. Actually, the sedimentary medium should be anisotropic as preferred cracks or fractures and thin layers are common features in the upper crust, so the interpretation of Pg event needs to be taken account of seismic velocity anisotropy. Traveltime calculation is the base of data processing and interpretation. Here, we only study the type of elliptical anisotropy for the poor quality and insufficiency of DSS data. In this thesis, we first investigate the meaning of elliptical anisotropy in the study of crustal structure and attribute, then derive Pg event’s traveltime-offset relationship by assuming a linear increasing velocity model with elliptical anisotropy and present the invert scheme from Pg traveltime-offset dataset to seismic velocity and its anisotropy of shallow crustal structure. We compare the Pg traveltime calculated by our analytic formula with numerical calculating method to test the accuracy. To get the lateral variation of elliptical anisotropy along the profiling, a tomography inversion method with the derived formula is presented, where the profile is divided into rectangles. Anisotropic imaging of crustal structure and attribute is efficient method for crust study. The imaging result can help us interprete the seismic data and discover the attribute of the rock to analyze the interaction between layers. Traveltime calculation is the base of image. Base on the ray tracing equations, the paper present a realization of three dimension of layer model with arbitrary anisotropic type and an example of Pg traveltime calculation in arbitrary anisotropic type is presented. The traveltime calculation method is complex and it only adapts to nonlinear inversion. Perturbation method of travel-time calculation in anisotropy is the linearization approach. It establishes the direct relation between seismic parameters and travetime and it is fit for inversion in anisotropic structural imaging. The thesis presents a P-wave imaging method of layer media for TTI. Southeastern China is an important part of the tectonic framework concerning the continental margin of eastern China and is commonly assumed to comprise the Yangtze block and the Cathaysia block, the two major tectonic units in the region. It’s a typical geological and geophysical zone. In this part, we fit the traveltime of Pg phase by the raytracing numerical method. But the method is not suitable here because the inefficiency of numerical method and the method itself. By the analytic method, we fit the Pg and Sg and get the lateral variation of elliptical anisotropy and then discuss its implication. The northeastern margin of Qinghai-Tibetan plateau is typical because it is the joint area of Eurasian plate and Indian plate and many strong earthquakes have occurred there in recent years．We use the Pg data to get elliptical anisotropic variation and discuss the possible meaning.

大型反倾岩体变形破坏机制及稳定性分析———以金沙江宽谷坝址龙蟠边坡为例

A large number of catastrophic accidents were aroused by the instability and destruction of anti-dip rock masses in the worldwide engineering projects, such as hydropower station, mine, railways and so on. Problems in relation to deformation and failure about anti-dip rock slopes are significant for engineering geology research. This dissertation takes the Longpan slope in the Jinsha River as a case to study the deformation mechanism of large-scale anti-dip rock masses and the slope stability analysis method. The primary conclusions are as follows. The Dale Reach of Jinsha River, from Longpan to the debouchment of Chongjiang tributary, is located in the southeastern margin of the Qinghai-Tibet Plateau. Longpan slope is the right embankment of Dale dam, it is only 26 km to the Shigu and 18 km to Tiger Leaping Gorge. The areal geology tectonic structures here area are complicated and blurry. Base on the information of geophysical exploration (CSAMT and seismology) and engineering geological investigation, the perdue tectonic pattern of Dale Reach is put forward for the first time in this paper. Due to the reverse slip of Longpan fault and normal left-rotation of Baihanchang fault, the old faulted valley came into being. The thick riverbed sediments have layered characters of different components and corresponding causes, which attribute to the sedimentary environments according with the new tectonic movements such as periodic mountain uplifting in middle Pleistocene. Longpan slope consists of anti-dip alternate sandstone and slate stratums, and the deformable volume is 6.5×107m3 approximately. It was taken for an ancient landslide or toppling failure in the past so that Dale dam became a vexed question. Through the latest field surveying, displacement monitoring and rock masses deforming characters analyses, the geological mechanism is actually a deep-seated gravitational bending deformation. And then the discrete element method is used to simulate the deforming evolution process, the conclusion accords very well with the geo-mechanical patterns analyses. In addition strength reduction method based on DEM is introduced to evaluate the factor of safety of anti-dip rock slope, and in accordance with the expansion way of the shear yielding zones, the progressive shear failure mechanism of large-scale anti-dip rock masses is proposed for the first time. As an embankment or a close reservoir bank to the lower dam, the stability of Longpan slope especially whether or not resulting in sliding with high velocity and activating water waves is a key question for engineering design. In fact it is difficult to decide the unified slip surface of anti-dip rock slope for traditional methods. The author takes the shear yielding zones acquired form the discrete element strength reduction calculation as the potential sliding surface and then evaluates the change of excess pore pressure and factor of stability of the slope generated by rapid drawdown of ponded water. At the same time the dynamic response of the slope under seismic loading is simulated through DEM numerical modeling, the following results are obtained. Firstly the effective effect of seismic inertia force is resulting in accumulation of shear stresses. Secondly the discontinuous structures are crucial to wave transmission. Thirdly the ultimate dynamic response of slope system takes place at the initial period of seismic loading. Lastly but essentially the effect of earthquake load to bringing on deformation and failure of rock slope is the coupling effect of shear stresses and excess pore water pressure accumulation. In view of limitations in searching the critical slip surface of rock slope of the existing domestic and international software for limit equilibrium slope stability analyses, this article proposes a new method named GA-Sarma Algorithm for rock slope stability analyses. Just as its name implies, GA-Sarma Algorithm bases on Genetic Algorithm and Sarma method. GA-Sarma Algorithm assumes the morphology of slip surface to be a broken line with traceability to extend along the discontinuous surface structures, and the slice boundaries is consistent with rock mass discontinuities such as rock layers, faults, cracks, and so on. GA-Sarma Algorithm is revolutionary method that is suitable for global optimization of the critical slip surface for rock slopes. The topics and contents including in this dissertation are closely related to the difficulties in practice, the main conclusions have been authorized by the engineering design institute. The research work is very meaningful and useful for the engineering construction of Longpan hydropower station.

班公湖带多不杂超大型富金斑岩铜矿床的成岩成矿年代学、岩石学及高氧化岩浆¬流体-成矿作用

Duobuza copper deposit, newly discovered typical gold-rich porphyry copper deposit with superlarge potential, is located in the Tiegelong Mesozoic tectonic -magmatic arc of the southern edge of Qiangtang block and the northern margin of Bangonghu-Nujiang suture. Quartz diorite porphyrite and grandiorite porphyry, occurred in stock, are the main ore-bearing porphyries. As the emplacement of porphyry stock, a wide range of hydrothermal alteration has developed. Within the framework of the ore district, abundant hydrothermal magnetite developed, and the relationship between precipitation of copper and gold and hydrothermal magnetite seems much close. Correspondingly, a series of veinlets and network veinlets occurred in all alteration zones. Therefore, systematic research on such a superlarge high-grade Duobuza gold-rich porphyry copper deposit can fully revealed the metallogenic characteristics of gold-rich porphyry copper deposits in this region, establish metallogenetic model and prospecting criteria, and has important practical significance on the promotion of regional exploration. In addition, this research on it can enrich metallogenic theory of strong oxidation magma-fluid to gold-rich porphyry copper deposit, and will be helpful to understand the metallogenic characteristics in early of subduction of Gangdese arc stages and its entire evolution history of the Qinghai-Tibet Plateau, the temporal and spatial distribution of ore deposits and their geodynamics settings. Northern ore body of Duobuza copper deposit have been controlled with width (north-south) about 100 ~ 400 m, length (east-west) about 1400 m, dip of 200 °, angle of dip 65 °~ 80 °. And controlled resource amount is of 2.7 million tons Cu with grade 0.94% and 13 tons Au with 0.21g/tAu. Overall features of ore body are large scale, higher grade copper, gold-rich. Ore occurred in the body of granodiotite porphyry and quartz diorite porphyrite and its contact zone with wall rock. Through the detailed mapping and field work studies, some typies of alteration are identificated as follows: albitization, biotititation, sericitization, silication, epidotization, chloritization, carbonatization, illitization, kaolinization and so on. The range of alteration is more than 10km2. Wall alteration zone can be divided into potassic alteration, moderate argillization alteration, argillization, illite-hydromuscovite or propylitization from ore-bearing porphyry center outwards, but phyllic alteration has not well developed and only sericite-quartz veins occurred in local area. Moreover, micro-fracture is development in ore district , and correspondingly a series of veinlets are development as follows: biotite vein (EB type), K-feldspar-biotite-chalcopyrite-quartz vein, magnetite-antinolite-K-feldspar vein, quartz-chalcopyrite-magnetite veins (A-type), quartz-magnetite-biotite-K-feldspar vein, chalcopyrite veinlets in potassic alteration zone; (2) chalcopyrite occurring in the center vein–quartz vein (B type), chalcopyrite veinlets, chalcopyrite-gypsum vein in intermediate argillization alteration; (3) chalcopyrite- pyrite-quartz vein, pyrite-quartz vein, chalcopyrite-gypsum veins, quartz-gypsum- molybdenite-chalcopyrite vein in argillization alteration; (4) gypsum veins, quartz-(molybdenite)-chalcopyrite vein, quartz-pyrite vein, gypsum- chalcopyrite vein, potassium feldspar veinlets, Carbonate veins, quartz-magnetite veins in the wall rock. In short, various veins are very abundant within the framework of the ore district. The results of electronic probe microscopy analysis (EMPA) indicate that Albite (Ab 91.5~99.7%) occurred along the rim of plagioclase phenocryst and fracture, and respresents the earliest stages of alteration. K-feldspar (Or 75.1~96.9%) altered plagioclase phenocryst and matrix or formed secondary potassium feldspar veinlets. Secondary biotite occurred mainly in phenocryst, matrix and veinlets, belong to magnesium-rich biotite formed under the conditions of high-oxidation magma- hydrothermal. Chloritization developed in all alteration zones and alterd iron- magnesium minerals such as biotite and hornblende and then formed chlorite veinlets. As the temperature rises, Si in the tetrahedral site of chlorite decreased, and chlorite component evolved from diabantite to ripiolite. The consistent 280℃~360℃ of formation temperature hinted that chlorite formed on the same temperature range in all alteration zones. However, formation temperature range of chlorite from the gypsum-carbonate-chlorite vein was 190℃~220℃, and it may be the product of the latest stage of hydrothermal activity. The closely relationship between biotite and rutile indicate that most of rutiles are precipitated in the process of biotite alteration and recrystallization. In addition, the V2O3 concentration of rutile from ore body in Duobuza gold-rich porphyry copper deposit is >0.4％, indicate that V concentration in rutile has important significance on marking main ore body of porphyry copper deposit. Apatites from Duobuza deposit all are F-rich. And apatite in the wall rock contained low MnO content and relatively high FeO content, which may due to the basaltic composition of the wall rocks. The MnO in apatite from altered porphyry show a strong positive correlation with FeO. In addition, Cl/F ratio of apatite from wall rock was highest, followed by the potassic alteration zone and potassic alteration zone overprinted by moderate argillization alteration was the lowest. SO2 in Apatite are in the scope of 0 to 0.66%, biotite in the apatite has the highest SO2, followed by the potassic alteration zone, potassic alteration zone overprinted by moderate argillization alteration, and the lowest in the surrounding rocks, which may be caused by the decrease of oxygen fugacity of hydrothermal fluid and S exhaust by sulfide precipitation in potassic alteration. Magnetite in the wall rock have higher Cr2O3 and lower Al2O3 features compared with altered porphyry, this may be due to basalt wall rock generally has high Cr content. And magnetites have higher TiO2 content in potassic alteration than moderate argillization alteration overprinted by potassic alteration, argillization and wall rock, suggested that its formation temperature in potassic alteration was the highest among them. The ore minerals mainly are chalcopyrite and bornite, and Au contents of chalcopyrite, bornite, and pyrite are similar with chalcopyrite slightly higher. The Eu* negative anomaly of disseminated chalcopyrite was relatively lower than chalcopyrite in veinlets. Within a drill hole, the Eu* negative anomaly of disseminated chalcopyrite was gradually larger from bottom to top. Magnetite has the same distribution model, with obvious negative Eu* abnormal, and ΣREE in great changes. The gypsum has the highest ΣREE content and the obvious negative anomaly, and biotite obviously has the Eu* abnormal. Based on the petrographic and geochemical characteristics, five series of magmatic rocks can be broadly classified; they are volcanic rocks of the normal island arc, high-Nb basaltic rocks, adakites, altered porphyry and diorite. The Sr, Nd, Hf isotopes and geochemistry of various series of magmatic rock show that they may be the result of mixing between basic magma and various degrees of acid magma coming from lower crust melted by high temperature basic underplating from partial melting of the subduction sediment melt metasomatic mantle wedge. Furthermore S isotope and Pb isotope of the sulfide, ore-bearing porphyries and volcanic rocks indicated ore-forming source is the mantle wedge metasomatied by subduction sediment melt. Oxygen fugacity of magma estimated by Fe2O3/FeO of whole rock and zircon Ce4+/Ce3+ indicated that the oxidation of basalt-andesitic rocks is higher than ore-forming porphyry, and might imply high-oxidation characteristics of underplated basic magma. Its high oxidative mechanism is likely mantle sources metasomatied by subduction sediment magma, including water and Fe3+. And such high oxidation of basaltic magma is conducive to the mantle of sulfides in the effective access to melt. And the An component of dark part within plagioclase phenocryst zoning belong to bytownite (An 74%), and its may be a result of magma composition changes refreshment by basaltic magma injection. SHRIMP zircon U-Pb and LA-ICP-MS zircon U-Pb geochronology study showed that the intrusions and volcanic rocks from Duobuza porphyry copper deposit belong to early Cretaceous magma series (126~105Ma). The magma evolution series are as follows: the earliest diorite and diorite porphyrite → ore-bearing porphyry and barren grandiorite porphyry →basaltic andesite → diorite porphyrite → andesite → basaltic andesite, and magma component shows a evolution trend from intermediate to intermediate-acid to basic. Based on the field evidences, the formation age of high-Nb basalt may be the latest. The Ar-Ar geochronology of altered secondary biotite, K-feldspar and sericite shows that the main mineralization lasting a interval of about 4 Ma, the duration limit of whole magma-hydrothermal evolution of about 6 Ma, and possibly such a long duration limit may result in the formation of Duobuza super-large copper deposit. Moreover, tectonic diagram and trace element geochemistry of volcanic rocks and diorite from Duobuza porphyry copper deposit confirm that it formed in a continental margin arc environment. Zircon U-Pb age of volcanic rocks and porphyry fall in the range of 105~121Ma, and Duobuza porphyry copper deposit locating in the north of the Bangonghu- Nujiang suture zone, suggested that Neo-Tethys ocean still subducted northward at least early Cretaceous, and its closure time should be later than 105 Ma. Three major inclusion types and ten subtypes are distinguished from quartz phenocrysts and various quartz veins. Vapor generally coexisting with brine inclusions, suggest that fluid boiling may be the main ore-forming mechanism. Raman spectrums of fluid inclusions display that the content of vapor and liquid inclusion mainly contain water, and vapor occasionally contain a little CO2. In addition, the component of liquid inclusions mainly include Cl-, SO42-, Na+, K+, a small amount of Ca2+, F-; and Cl- and Na+ show good correlation. Vapor mainly contains water, a small amount of CO2, CH4 and C2H6 and so on. The daughter minerals identified by Laman spectroscopy and SEM include gypsum, chalcopyrite, halite, sylvite, rutile, potassium feldspar, Fe-Mn-chloride and other minerals, and ore-forming fluid belong to a complex hydrothermal system containing H2O-NaCl-KClFeCl2CaCl2. H and O isotopic analysis of quartz phenocryst, vein quartz, magnetite, chlorite and gypsum from all alteration zones show that the ore-forming fluid of Duobuza gold-rich porphyry copper deposit consisted mainly of magmatic water, without addition of meteric water. Duobuza gold-rich porphyry copper deposit formed by the primary magmatic fluid (600-950C), which has high oxidation, ultra-high salinity and metallogenic element-rich, exsolution direct from the magma, and it is representative of the typical orthomagmatic end member of the porphyry continuum. Moreover, the fluid evolution model of Duobuza gold-rich porphyry copper deposit has been established. Furthermore, two key factors for formation of large Au-rich porphyry copper deposit have been summed up, which are ore-forming fluids earlier separated from magma and high oxidation magma-mineralization fluid system.

塔里木地块奥陶纪古地磁新结果及其构造意义

The Tarim Block is located between the Tianshan Mountains in the north and the Qinghai-Tibet Plateau in the south and is one of three major Precambrian cratonic blocks of China. Obviously, the Paleozoic paleogeographic position and tectonic evolution for the Tarim Block are very important not only for the study of the formation and evolution of the Altaids, but also for the investigation of the distributions of Paleozoic marine oil and gas in the Tarim Basin. According to the distributions of Paleozoic strata and suface outcrops in the Tarim Block, the Aksu-Keping-Bachu area in the northwestern part of the Tarim Block were selected for Ordovician paleomagnetic studies. A total of 432 drill-core samples form 44 sampling sites were collected and the samples comprise mainly limestones, argillaceous limestones and argillaceous sandstones Based on systematic study of rock magnetism and paleomagnetism, all the samples could be divided into two types: the predominant magnetic minerals of the first type are hematite and subordinate magnetite. For the specimens from this type, characteristic remanent magnetization (ChRM) could generally be isolated by demagnetization temperatures larger than 600℃; we assigned this ChRM as component A; whilst magnetite is the predominant magnetic mineral of the second type; progressive demagnetization yielded another ChRM (component B) with unblocking temperatures of 550-570℃. The component A obtained from the majority of Ordovician specimens has dual polarity and a negative fold test result; we interpreted it as a remagnetization component acquired during the Cenozoic period. The component B can only be isolated from some Middle-Late Ordovician specimens with unique normal polarity, and has a positive fold test result at 95% confidence. The corresponding paleomagnetic pole of this characteristic component is at 40.7°S, 183.3°E with dp/dm = 4.8°/6.9° and is in great difference with the available post-Late Paleozoic paleopoles for the Tarim Block, indicating that the characteristic component B could be primary magnetization acquired in the formation of the rocks. The new Ordovician paleomagnetic result shows that the Tarim Block was located in the low- to intermediate- latitude regions of the Southern Hemisphere during the Middle-Late Ordovician period, and is very likely to situate, together with the South China Block, in the western margin of the Australian-Antarctic continents of East Gondwana. However, it may have experienced a large northward drift and clockwise rotation after the Middle-Late Ordovician period, which resulted in the separation of the Tarim Block from the East Gondwanaland and subsequent crossing of the paleo-equator; by the Late Carboniferous period the Tarim Block may have accreted to the southern margin of the Altaids.

中国北方风成堆积的风化作用与环境变迁—U-Th、~(10)Be及元素地球化学的研究

The surface of the Earth is continuously undergoing changes as a result of weathering-erosion, plate tectonics and volcanic processes. Continental weathering-erosion with its complex rock-water interactions is the central process of global biochemical cycling of elements, and affects the long-term ocean atmosphere budget of carbon dioxide both through the consumption of carbonic acid during silicate weathering and through changes in the weathering and burial rates of organic carbon. Rates of the weathering-erosion depend on a variety of factors, in particular rock properties and chemical composition, climate (especially rainfall), structure, and elevation. They are quite variable on a regional scale. Thus, environmental changes in a region could be indicated by the history of weathering-erosion in the region. Recent attention has focused on increased silicate weathering of tectonically uplifted areas in the India-Asia collision zone as a possible cause for falling atmospheric CO_2 levels in the Cenozoic era. The wind blown dust deposits in the Loess Plateau is derived from the arid and semiarid regions in northwestern China, in turn, where the deposits have been derived from the Qinghai-Xizang Plateau and the high mountains around. Therefore, geochemistry of the wind blown loess-paleosol and red clay sequences may provide insight both to paleoenvironmental changes on the Loess Plateau, and to the uplift and weathering-erosion histories of the Qinghai-Xizang Plateau. In this paper, uranium-thorium series nuclides and cosmogenic ~(10)Be have been employed as tracers of weathering intensities and histories of the dust sediments in the Loess Plateau. Major elements, such as Na, Al, Fe etc., are also used to estimate degree of chemical alteration of the dust sediments and to rebuild the history of weathering on the Loess Plateau. First of all, using a low-level HPGe γ-ray detector, we measured U and Th series nuclides in 170 loess and paleosol samples from five sites in the Loess Plateau, going back 2.6 Ma. The results show that ~(238)U activities are disequilibrium with its daughter nuclide ~(230)Th in young loess-paleosol sequence, indicating that weathering was happened both in dust deposition site and in dust source regions. Using concentrations of ~(238)U and ~(232)Th in the samples, we estimated the amounts of ~(238)U leached out of from paleosols due to weathering. Further, based on analyses of ~(230)Th in paleosols deposited in the past ca. 140 ka, we determined when the paleosols weathered in the source regions. We conclude that most of the weathering in the dust-source regions may have occurred during the interglacials before dust deposition.

藏东及滇东南新生代钾质岩浆作用及其深部制约

The dissertation focuses on the petrology, geochemistry of the volcanic rocks in east Tibet and southeast Yunnan. It lucubrates the Magmatic process, forming mechanism and the possible tectonic settings of the volcanic rocks. The volcanic rocks of Nangqen basin in east Tibet, Qinghai province are mainly Cenozoic intermediate-acid shoshonites. The rocks are LREE enriched and the LREE/HREE = 3～34; (La/Yb)_N = 18.17-53.59, and ΣREE 222～1260μg/g. There are no Eu anomaly, and Nb, Ta, Zr, Hf, Ti are markedly depleted. The isotopic composition is ~(87)Sr/~(86)Sr = 0.70497～0.70614, ~(206)Pb/~(204)Pb = 18.622～18.974, ~(208)Pb/~(204)Pb = 38.431～38.996, ~(207)Pb/~(204)Pb = 15.511～15.613, respectively. K-Ar age of the whole rocks and the single mineral are between 32.0-36.5Ma. Based on the trace elements and isotopic elements, we get the conclusion that the partial melting is one of the dominated forming mechanisms for the volcanic rocks in Naneqen basin. The magma did not experience the crustal contamination en route to the surface; however, the complex mixture took place in the upper mantle before the melt was formed. There are at least two kinds of mixed sources that can be identified. The basalt in southeast Yunnan province is studied. They are distributed in Maguan, Tongguan, and Pingbian County, which is located on the both sides of the Red River belt, and the ultrabasic xenolith are cursory introduced. The volcanic rocks belongs to the alkali series, which can be subdivided into trachybasalt and basanite(Ol normal molecule >5). The volcanic rocks are characteristics by high Ti and low Mg#. According to the magma calculation model, the original rocks of the basalt in southeast Yunnan province are Spinel Lherzolite in Tongguan, Garnet Lherzolite in Pingbian and Maguan, while Togguan undergoes 2-5 percent and percent of partial melting, whereas volcanism in Maguan and Pingbian was so complex to calculate. The fractional crystallization took place during the magma evoltion in southeast Yunnan. The basalt is enriched in LREE with LREE/HREE=9.23-20.19. All of the trace elements display weak Nb, Ta peak, and the depletion of Zr, Hf and Ti in Maguan and pingbian represent the presence of Garnet in the source. The composition of the isotope ratio are ~(87)Sr/~(86)Sr = 0.70333-0.70427, ~(143)Nd/~(144)Nd = 0.512769-0.512940, ~(206)Pb/~(204)Pb = 18.104-18.424, ~(207)Pb/~(204)Pb = 15.483 -15.527; ~(208)Pb/~(204)Pb = 37.938-38.560, respectively, which shows the characteristics of the HIMU type OIB. The volcanic rocks of the southwest Yunnan are derived from the enriched, OIB type mantle sources by synthesizing all the data from trace and isotope elements. It is similar to that of the volcanic rocks in Hawaii, a typical kind of the mixtures of the recycled oceanic crust plume and depleted asthenosphere. To sum up, the volcanic rocks in southeast Yunnan are formed by the intraplate hotpot volcanism.

新第三纪红河盆地特征及红河断裂活动对其形成与发育的控制

Based on geodynamic analysis of sedimentary basins, combined sedimentology with structural geology and other methods, the author studied the Honghe basin located in Yunnan province of Southwestern China. Sandstone slice grain size analysis, combined with field geology and indoors study indicate that a set of inland alluvial fan diposits, fluvial deposites, delta deposits and some lacustrine sediments are in Honghe basin. Studying on shape of the Honghe basin, sedimentary and structural characteristic and distribution of different kinds of conglomerate and its structural significance, we hold the idea that the formation and evolution of Honghe basin are controlled by the activity of Red River faut. Correlation of lithostratic cross section in Honghe basin and studying on activity of Red River fault indicate that Honghe basin was formed in two stages. It is a complex basin constitutes of the first-stage trans-releasing basin and the second-stage trans-downfaulted basin. Due to the uplift of Qinghai-Xizang plateau and deformation of orogeny, the western Yunnan and adjacent area move to SE direction as a result of Tectonic Escape. Right lateral strike slip occurred along Red River fault, trans-releasing basin formed at the bend part of the fault due to stress relexation. As the block escaping, it moves away from the other block of the Red River fault, the upper block move down obliquely and trans-downfaulted basin formed. Combined the age of phytolite and regional structural events, we think the first-stage transreleasing basin was formed in late Miocene, on the other words, the dextral strike slip of Red River fault may began in late Miocene (10-7Ma). The second-stage trans-downfaulted basin may be formed in early stage of Pliocene (about 4.7Ma). Subsequently, the bilateral faults dipping to the inside of the plateau and thrusting outwards occurred in the marginal region of Qinghai-Xizang plateau during its uplifting as a fan-shaped mountain body, this results in the uplift of the strata to the east of Red River fault and supply large quantity of provenance for the Honghe basin. In last Pliocene (about 3Ma), strong uplift of Qinghai-Xizang plateau leads to massive clastic sediment entered Honghe basin and causes its closure. As a kind of trans-tentional basin, trans-releasing basin is different to pull-apart basin. The author compared the Mosha trans-releasing basin with Jinggu pull-apart basin in SW Yunan, China, and described their character correspondingly. Otherwise, the author combined the predecessors' studding with conclusion of own study, discussed the kinematics of Ailaoshan-Red River belt in Cenozoic, and the relationship between the formation of Honghe basin and uplifting of Qinghai-Xizang plateau.

地震波场模拟及叠前深度偏移的李群方法

The processes of seismic wave propagation in phase space and one way wave extrapolation in frequency-space domain, if without dissipation, are essentially transformation under the action of one parameter Lie groups. Consequently, the numerical calculation methods of the propagation ought to be Lie group transformation too, which is known as Lie group method. After a fruitful study on the fast methods in matrix inversion, some of the Lie group methods in seismic numerical modeling and depth migration are presented here. Firstly the Lie group description and method of seismic wave propagation in phase space is proposed, which is, in other words, symplectic group description and method for seismic wave propagation, since symplectic group is a Lie subgroup and symplectic method is a special Lie group method. Under the frame of Hamiltonian, the propagation of seismic wave is a symplectic group transformation with one parameter and consequently, the numerical calculation methods of the propagation ought to be symplectic method. After discrete the wave field in time and phase space, many explicit, implicit and leap-frog symplectic schemes are deduced for numerical modeling. Compared to symplectic schemes, Finite difference (FD) method is an approximate of symplectic method. Consequently, explicit, implicit and leap-frog symplectic schemes and FD method are applied in the same conditions to get a wave field in constant velocity model, a synthetic model and Marmousi model. The result illustrates the potential power of the symplectic methods. As an application, symplectic method is employed to give synthetic seismic record of Qinghai foothills model. Another application is the development of Ray+symplectic reverse-time migration method. To make a reasonable balance between the computational efficiency and accuracy, we combine the multi-valued wave field & Green function algorithm with symplectic reverse time migration and thus develop a new ray+wave equation prestack depth migration method. Marmousi model data and Qinghai foothills model data are processed here. The result shows that our method is a better alternative to ray migration for complex structure imaging. Similarly, the extrapolation of one way wave in frequency-space domain is a Lie group transformation with one parameter Z and consequently, the numerical calculation methods of the extrapolation ought to be Lie group methods. After discrete the wave field in depth and space, the Lie group transformation has the form of matrix exponential and each approximation of it gives a Lie group algorithm. Though Pade symmetrical series approximation of matrix exponential gives a extrapolation method which is traditionally regarded as implicit FD migration, it benefits the theoretic and applying study of seismic imaging for it represent the depth extrapolation and migration method in a entirely different way. While, the technique of coordinates of second kind for the approximation of the matrix exponential begins a new way to develop migration operator. The inversion of matrix plays a vital role in the numerical migration method given by Pade symmetrical series approximation. The matrix has a Toepelitz structure with a helical boundary condition and is easy to inverse with LU decomposition. A efficient LU decomposition method is spectral factorization. That is, after the minimum phase correlative function of each array of matrix had be given by a spectral factorization method, all of the functions are arranged in a position according to its former location to get a lower triangular matrix. The major merit of LU decomposition with spectral factorization (SF Decomposition) is its efficiency in dealing with a large number of matrixes. After the setup of a table of the spectral factorization results of each array of matrix, the SF decomposition can give the lower triangular matrix by reading the table. However, the relationship among arrays is ignored in this method, which brings errors in decomposition method. Especially for numerical calculation in complex model, the errors is fatal. Direct elimination method can give the exact LU decomposition But even it is simplified in our case, the large number of decomposition cost unendurable computer time. A hybrid method is proposed here, which combines spectral factorization with direct elimination. Its decomposition errors is 10 times little than that of spectral factorization, and its decomposition speed is quite faster than that of direct elimination, especially in dealing with a large number of matrix. With the hybrid method, the 3D implicit migration can be expected to apply on real seismic data. Finally, the impulse response of 3D implicit migration operator is presented.

各向异性介质中的地震面波和S波分裂

Surface wave propagation in the anisotropic media and S-wave splitting in China mainland are focused in this M.S. dissertation. We firstly introduced Anderson parameters in the research of surface wave propagation in the anisotropic media were deduced, respectively. By applying the given initial model to the forward calculation of Love wave, we compared dispersion curves of Love wave in the anisotropic media with the one in the isotropic media. the results show that, although the two kind of results are similar with each other, the effect of anisotropy can not be neglected. Furthermore, the variation of anisotropy factors will result in the variation of dispersion curves, especially for high-mode one. The method of grid dispersion inversion was then described for further tectonic inversion. We also deduced inversion equation on the condition that the layered media is anisotropic, and calculated the phase-velocity partial derivatives with respect to the model parameters, P- and S-wave velocities, density, anisotropic parameters for Rayleigh wave and Love wave. Having analyzed the results of phase-velocity partial derivatives, we concluded that the derivatives within each period decreased with the depth increasing, the phase-velocity of surface wave is sensitive to the S-wave velocities and anisotropic factors and is not sensitive to the densities of layers. Dispersion data of Love wave from the events occurred during the period from 1991 to 1998 around the Qinghai and Tibet Plateau, which magnitudes are more than 5.5, have been used in the grid dispersion inversion. Those data have been preprocessed and analyzed in the F-T domain. Then the results of 1°*1° grid dispersion inversion, the pure path dispersion data, in the area of Qianghai and Tibet Plateau were obtained. As an example, dispersion data have been input for the tectonic inversion in the anisotropic media, and the results of anisotropic factors under the region of Qianghai and Tibet Plateau were initially discussed. As for the other part of this dissertation. We first introduced the phenomena of S-wave splitting and the methods for calculation the splitting parameters. Then, We applied Butterworth band-pass filter to S-wave data recorded at 8 stations in China mainland, and analyzed S-wave splitting at different frequency bands. The results show the delay time and the fast polarization directions of S-wave splitting depend upon the frequency bands. There is an absence of S-wave splitting at the station of Wulumuqi (WMQ) for the band of 0.1-0.2Hz. With the frequency band broaden, the delay time of S-wave splitting decreases at the stations of Beijing (BJI), Enshi (ENH), Kunming (KMI) and Mudanjiang (MDJ); the fast polarization direction at Enshi (ENH) changes from westward to eastward, and eastward to westward at Hailaer (HIA). The variations of delay time with bands at Lanzhou (LZH) and qiongzhong (QIZ) are similar, and there is a coherent trend of fast polarization directions at BJI, KMI and MDJ respectively. Initial interpretations to the results of frequency band-dependence of S-wave splitting were also presented.

青藏高原东缘新生代富钾碱性岩的年龄与岩石地球化学特征

The Cenozoic high-Potassium alkaline magma rocks are widely distributed in the eastern margin of Qinghai-Xizang (Tibet) plateau, such as Mangkang, Deqin, Jianchuan, Haidong in the western Yunnan Province, Yanyuan, Muli in the western Sichuan province, including syenite, minette and trachyte and lamprophyres, which formaed a large alkaline-rich magmatic rock zone along the Jinshajiang-Ailaoshan fault. Here we present detailed Geochemical and Geochronology study for the Cenozoic high-K alkline magma rocks from Mangkang, Deqin, Jianchuan, Haidong and Muli, Yanyuan, Yanbian, Yaoan (western Yangtze craton). Rocks of our target area include syenite, minette and trachyte. Geochronologic dating and compilation of existing age data suggest magmatism occurred at 32-38Ma and distribute along the Jinsha suture, from Jinping, Daping, Zhenyuan to Haidong, Jianchuan, Deqin, Mangkang, Zaduo to Qiangtang Bandaohu, Zhentouya. They all share same geochemical features such as rich in potassium and alkaline, rich in LREE, slightly negative Eu anomalies, negative anomalies in Nb, Ta, Ti,. These indicate they come from EMU mantle. Geochronologic dating on western Yangtze craton minette surgest minette in Haidong, Yanyuan, Yanbian, Yaoan, Deqin ocurred at 30-37Ma, well minette in Muli ocurred at IMa. Need more detailed geology study to find differece of these two episode magmatism.

碳、氧稳定同位素分析在古环境研究中的应用

A number of proxy records of paleoenvironment using stable isotopes could show the history of past environmental changes. These archives include peat and lake sediments, loess-paleosot sequence, fossil mammals and stalagmite, and so on. The stable isotopic composition of carbonate and organic matter and frequency magnetic susceptibility from Tianshuigou and Yuanlei loess-palesol sequence can be used to give estimates of the paleoenvironmental history of Dali, and even of the whole Chinese Loess Plateau during the last 250ka. Features of the High Temperature and Large Precipitation Event in the Tibet Plateau and its adjacent area during 40~30kaBP had been studied by Professor Shi Y. In this dissertation, its impact on Chinese Loess Plateau has been discussed again. Carbon and oxygen isotopic ratios, magnetic susceptibility and frequency magnetic susceptibility in Tianshuigou and Yuanlei profiles show that the Event in this area is not so stronger as the Tibet Plateau. The carbon isotopic composition of organic matter in Tianshuigou, Yuanlei, dingcun and Jingcun loess-palesol sequences are indicative of major changes in the paleovagetation between terrace and plain of the Chinese Loess Plateau. Water is one of the most important factors adjusting the relative biomass of C4 plant in terrestrial ecosystems. Stable carbon isotope ratio of vertebrate tooth enamel is used increasingly to reconstruct environmental and ecological information modern and ancient ecosystems. The SI3C value of tooth enamel bioapatites can distinguish between browsers and grazers. Data from typical grassland of Inter Mongolia, the Alpine meadow of Qinghai-Tibet Plateau and the Yaluzangbu Great Canyon indicate that diets of mammals could record the relative biomass of C4 plant only in the C4 dominated ecosystem. In a C3 dominated ecosystem, diet of mammals would include more C3 plants than vegetation. According to Professor Cerling, proxy records from North and South America, Africa and Pakistan show that at the end of the Miocene (between 8Ma to 6 Ma) there was a global expansion of CA biomass, probably when atmospheric CO2 levels declined. Thus, "C4 world" and "CO2 starvation" are put forward. In this dissertation, carbon isotopes of fossil tooth such as Equus sanmeniensis and Hipparion chiai from Linxia, China reveal that there is a C3 dominated ecosystem in the late Miocene. Diets of ancient mammals in Linxia are not evidence of global expansion of C4 biomass.