Effects of vegetation control on ecosystem water use efficiency within and among four grassland ecosystems in China

Through 2-3-year (2003-2005) continuous eddy covariance measurements of carbon dioxide and water vapor fluxes, we examined the seasonal, inter-annual, and inter-ecosystem variations in the ecosystem-level water use efficiency (WUE, defined as the ratio of gross primary production, GPP, to evapotranspiration, ET) at four Chinese grassland ecosystems in the Qinghai-Tibet Plateau and North China. Representing the most prevalent grassland types in China, the four ecosystems are an alpine swamp meadow ecosystem, an alpine shrub-meadow ecosystem, an alpine meadow-steppe ecosystem, and a temperate steppe ecosystem, which illustrate a water availability gradient and thus provide us an opportunity to quantify environmental and biological controls on ecosystem WUE at different spatiotemporal scales. Seasonally, WUE tracked closely with GPP at the four ecosystems, being low at the beginning and the end of the growing seasons and high during the active periods of plant growth. Such consistent correspondence between WUE and GPP suggested that photosynthetic processes were the dominant regulator of the seasonal variations in WUE. Further investigation indicated that the regulations were mainly due to the effect of leaf area index (LAI) on carbon assimilation and on the ratio of transpiration to ET (T/ET). Besides, except for the swamp meadow, LAI also controlled the year-to-year and site-to-site variations in WUE in the same way, resulting in the years or sites with high productivity being accompanied by high WUE. The general good correlation between LAI and ecosystem WUE indicates that it may be possible to predict grassland ecosystem WUE simply with LAI. Our results also imply that climate change-induced shifts in vegetation structure, and consequently LAI may have a significant impact on the relationship between ecosystem carbon and water cycles in grasslands.

Summer diet of two sympatric species of raptors Upland Buzzard (Buteo hemilasius) and Eurasian Eagle Owl (Bubo bubo) in alpine meadow: problem of coexistence

Summer diets of two sympatric raptors Upland Buzzards (Buteo hemilasius Temminck et Schlegel) and Eurasian Eagle Owls (Bubo bubo L. subsp. Hemachalana Hume) were studied in an alpine meadow (3250 m a.s.l.) on Qinghai-Tibet Plateau, China. Root voles Microtus oeconomus Pallas, plateau pikas Ochotona curzoniae Hodgson, Gansu pikas O. cansus Lyon and plateau zokors Myospalax baileyi Thomas were the main diet components of Upland Buzzards as identified through the pellets analysis with the frequency of 57, 20, 19 and 4%, respectively. The four rodent species also were the main diet components of Eurasian Eagle Owls basing on the pellets and prey leftovers analysis with the frequency of 53, 26, 13 and 5%, respectively. The food niche breadth indexes of Upland Buzzards and Eurasian Eagle Owls were 1.60 and 1.77 respectively (higher value of the index means the food niche of the raptor is broader), and the diet overlap index of the two raptors was larger (C-ue = 0.90) (the index range from 0 - no overlap - to I - complete overlap). It means that the diets of Upland Buzzards and Eurasian Eagle Owls were similar (Two Related Samples Test, Z = -0.752, P = 0.452). The classical resource partitioning theory can not explain the coexistence of Upland Buzzards and Eurasian Eagle Owls in alpine meadows of Qinghai-Tibet Plateau. However, differences in body size, predation mode and activity rhythm between Upland Buzzards and Eurasian Eagle Owls may explain the coexistence of these two sympatric raptors.

Response of chlorophyll fluorescence to dynamic light in three alpine species differing in plant architecture

A study was carried out to examine the effect of dynamic photosynthetically active photon flux density (PPFD) on photoinhibition and energy use in three herbaceous species, prostrate Saussurea superba, erect-leaved S. katochaete, and half-erect-leaved Gentiana straminea, from the Qinghai-Tibet Plateau. Chlorophyll fluorescence response was measured under each of three sets of high-low PPFD combinations: 1700-0, 1400-300, and 1200-500 mu mol m(-2) s(-1), illuminating in four dynamic frequencies: 1, 5, 15, and 60 cycles per 2 h. The total light exposure time was 2h and the integrated PPFD was the same in all treatments. The highest frequency of PPFD fluctuation resulted in the lowest photochemical activity, the highest level of non-photochemical quenching, and the greatest decrease of F-v/F-m (maximal photochemical efficiency of PSII). The 5 and 15 cycles per 2h treatments resulted in higher photochemical activity than the 1 cycle per 2h treatment. The 1700-0 PPFD combination led to the lowest photochemical activity and more serious photoinhibition in all species. S. superba usually exhibited the highest photochemical activity and CO2 uptake rate, the lowest reduction of F-v/F-m,F- and the smallest fraction of energy in thermal dissipation. With similar fractions of thermal dissipation, S. katochaete had relatively less photoinhibition than G. straminea owing to effective F-o quenching. The results suggest that high frequency of fluctuating PPFD generally results in photoinhibition, which is more serious under periods of irradiation with high light intensity. (c) 2005 Elsevier B.V. All rights reserved.

Effect of burrowing activity of root vole (Microtus oeconomus Pallas) on plant diversity in alpine meadow

Herbivory and burrowing activity of mammals may influence the species composition and diversity of plant communities. The effect of corridors and holes systems constructed by root vole (Microtus oeconomus Pallas) on the plant species diversity was studied in the habitat of high - mountain meadow (3250 in a.s.l in Qinghai-Tibet Plateau, China). By using grid method, these disturbances were studied on 16 plots (100 cm x 100 cm) distributed in 4 transects in studied area, in August 2000 and 2001. The disturbance intensity index, D, was calculated as the percent of the ground surface disturbed by voles in the study area. Plant species were identified and counted on the same plots. In total 46 plant species were identified - 39% of this number was considered as sensitive to the vole disturbances as their occurrence and/or abundance decreased along the disturbance intensity. Generally, a significantly negative correlation (r = - 0.911 P < 0.01) between vole aboveground disturbances and plant species diversity (H') was found. The results suggest that root vole ground disturbances, especially in the form of actively utilized holes and corridors have significantly negative influence on plant species diversity in high-mountain grassland habitat.

Local warming about 1.3 degrees C in alpine meadow has no effect on root vole (Microtus oeconomus L.) population during winter

The influence of air and soil warming on root vole (Microtus oeconomus L.) population was studied in winter period in top open chambers (OTC) (0.8-1.8 m(2) warmed by conical fiberglass material and situated in alpine meadow (3250 m) at Qinghai-Tibet Plateau, China. The OTCs were distributed on an area of 30 x 30 m of experimental warming site; another site of the same area was a control one. The root vole population was investigated on two pairs of sites in "low-grazing" and "high-grazing" (by sheep) parts of the meadow; mark-recapture method was used. The winter-season averaged air and soil temperature inside of the chambers were 1.3 degrees C higher than the temperature outside the chambers. The warming in the chambers had no statistically significant effect on root vole numbers, on average body mass of individual, and on average body mass of males and females. In conclusion, as small as 1.3 degrees C warming of soil and air introduced locally and on small (several m(2)) scale, in the alpine meadow habitat in winter period, has possibly no effect on root vole numbers and biomass.

Phylogenetic analyses of Saussurea sect. Pseudoeriocoryne (Asteraceae : Cardueae) based on chloroplast DNA trnL-F sequences

The genus Saussurea is distributed mainly in the temperate and subarctic regions of Eurasia and consists of about 300 species classified into six subgenera and 20 sections. Sect. Pseudoeriocoryne in the subgenus Eriocoryne comprises four species, and is delimited mainly by acaulescence and an inflorescence with congested capitula surrounded by a rosette of leaves. All of these species are endemic to the and Qinghai-Tibet Plateau. Sequences from the chloroplast DNA trnL-F region were obtained for the four species in this section and 26 other species from four subgenera of Saussurea to resolve phylogenetic relationships among these species and to determine whether the shared characters that define sect. Pseudoeriocoryne are synapomorphic or were acquired by convergent evolution. The resulting phylogenies indicated that Saussurea sect. Pseudoeriocoryne as traditionally defined does not constitute a monophyletic group and that each of its species belongs to separate clades. Furthermore, none of these species showed a close relationship with the other species of subgenus Eriocoryne. Our results further indicated that none of the investigated subgenera are monophyletic, and that species from different subgenera clustered together. All these conclusions are provisional and their confirmation would require stronger phylogenetic support. Two possible explanations are suggested for low sequence divergence, poor resolution of internal clades and clustering of species with the rather distinct morphology of Saussurea detected in the present study. The first is rapid radiation and diversification triggered by fast habitat fragmentation due to the recent lifting of the Qinghai-Tibet Plateau and the Quaternary climate oscillations. This could have led to rapid morphological divergence while sequences diverged very little, and also caused the convergent acquisition of similar characteristics in unrelated lineages due to similar selection pressures. The second possible explanation is that both introgressive hybridization and reticulate evolution might have caused the transferring of cpDNA sequences between morphologically dissimilar species, thus leading to homogenization of sequences between lineages. (C) 2004 Elsevier Ltd. All rights reserved.

Photosynthesis of Saussurea superba and Gentiana straminea is not reduced after long-term enhancement of UV-B radiation

Experiments were conducted in an alpine Kobresia humilis meadow near Haibei Alpine Meadow Ecosystem Research Station (37degrees29'-37degrees45'N, 101degrees12'-101degrees33'E; altitude 3200 m). Effects of enhanced ultraviolet-B (UV-B) radiation on photosynthesis of the alpine plants of Saussurea superba and Gentiana straminea were investigated. Both species were exposed to a UV-B-BE density at 15.80 kJ m(-2) per day, simulating nearly 14% ozone (O-3) reduction during the plant growing season. Neither photosynthetic CO2 uptake rate nor photosynthetic O-2 evolution rate were decreased after a long period of enhanced UV-B radiation treatment. On the contrary, there was a tendency to increase of both parameters in both species. The photosynthetic pigments were also increased, when expressed on a leaf area basis. UV-B absorbing compounds, detected by the absorbance values at 300 mm, had a tendency to increase in both species after enhanced UV-B radiation. After long-term exposure of plants to enhanced UV-B radiation, leaf morphology was also affected. Leaf thickness in both S. superba and G. straminea were increased significantly (P < 0.001). This supports our hypothesis that the increase of leaf thickness in both species after long-term exposure of enhanced UV-B radiation could compensate for the photodestruction of photosynthetic pigments when light passes through the leaf. Therefore, photosynthesis is not reduced in either species when expressed on leaf area basis. (C) 2003 Elsevier B.V. All rights reserved.

Photosynthetic depression in relation to plant architecture in two alpine herbaceous species

The Qinghai-Tibet Plateau is characterized by extremely high radiation, which may induce down-regulation of photosynthesis in plants living in this alpine ecosystem. To clarify whether photoinhibition occurs in the alpine environment and to discern its underlying mechanisms, we examined photosynthetic gas exchange and fluorescence emission in response to the changes in photosynthetic photon flux density (PPFD) and leaf temperature under natural regimes for two herbaceous species: prostrate Saussurea superba and erect-leaved Saussurea katochaete from altitude 3250 m on the Qinghai-Tibet Plateau. S. superba intercepted a higher maximum PPFD and experienced much higher leaf temperature than the erect-leaved S. katochaete. S. superba exhibited a much higher light saturation point for photosynthesis than S. katochaete. Under controlled conditions, the former species had higher CO2 uptake rates and neither species showed obvious photosynthetic down-regulation at high PPFD. Under natural environmental conditions, however, apparent photoinhibition, indicated by reduced electron transport rate (ETR), was evident at high PPFD for both species. After a night frost, the photochemistry of S. katochaete was depressed markedly in the early morning and recovered by mid-day. After a frost-free night, it was high in the morning and low at noon due to high radiation. S. superba did not respond to the night frost in terms of daytime photochemical pattern. In both species, photochemical depression was aggravated by high leaf temperature and the erect species was more sensitive to high temperature. This study suggests that the high radiation on the Qinghai-Tibet Plateau is likely to induce rapidly reversible photoinhibition, which is related closely to plant architecture. Photochemistry in the prostrate species seems able to tolerate higher PPFD, without obvious suppression, than the erect species. (C) 2003 Elsevier Science B.V. All rights reserved.

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

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.

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

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.

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

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.

地幔对流理论及其在解释地球物理观测的应用

With the improving of mantle convection theory, the developing of computing method and increasing of the measurement data, we can numerically simulate more clearly about the effects on some geophysical observed phenomenons such as the global heat flow and global lithospheric stress field in the Earth's surface caused by mantle convection, which is the primary mechanism for the transport of heat from the Earth's deep interior to its surface and the underlying force mechanism of dynamics in the Earth.Chapter 1 reviews the historical background and present research state of mantle convection theory.In Chapter 2, the basic conception of thermal convection and the basic theory about mantle flow.The effects on generation and distribution of global lithospheric stres s field induced by mantle flow are the subject of Chapter 3. Mantle convection causes normal stress and tangential stresses at the bottom of the lithosphere, and then the sublithospheric stress field induces the lithospheric deformation as sixrface force and results in the stress field within the lithosphere. The simulation shows that the agreement between predictions and observations is good in most regions. Most of subduction zones and continental collisions are under compressive. While ocean ridges, such as the east Pacific ridge, the Atlantic ridge and the east African rift valley, are under tensile. And most of the hotspots preferentially occur in regions where calculated stress is tensile. The calculated directions of the most compressive principal horizontal stress are largely in accord with that of the observation except for some regions such as the NW-Pacifie subduction zone and Qinghai-Tibet Plateau, in which the directions of the most compressive principal horizontal stress are different. It shows that the mantel flow plays an important role in causing or affecting the large-scale stress field within the lithosphere.The global heat flow simulation based on a kinematic model of mantle convection is given in Chapter 4. Mantle convection velocities are calculated based on the internal loading theory at first, the velocity field is used as the input to solve the thermal problem. Results show that calculated depth derivatives of the near surface temperature are closely correlated to the observed surface heat flow pattern. Higher heat flow values around midocean ridge systems can be reproduced very well. The predicted average temperature as a function of function of depth reveals that there are two thermal boundary layers, one is close to the surface and another is close to the core-mantle boundary, the rest of the mantle is nearly isothermal. Although, in most of the mantle, advection dominates the heat transfer, the conductive heat transfer is still locally important in the boundary layers and plays an important role for the surface heat flow pattern. The existence of surface plates is responsible for the long wavelength surface heat flow pattern.In Chapter 5, the effects on present-day crustal movement in the China Mainland resulted from the mantle convection are introduced. Using a dynamic method, we present a quantitative model for the present-day crustal movement in China. We consider not only the effect of the India-Eurasia collision, the gravitational potential energy difference of the Tibet Plateau, but also the contribution of the shear traction on the bottom of the lithosphere induced by the global mantle convection. The comparison between our results and the velocity field obtained from the GPS observation shows that our model satisfactorily reproduces the general picture of crustal deformation in China. Numerical modeling results reveal that the stress field on the base of the lithosphere induced by the mantle flow is probably a considerable factor that causes the movement and deformation of the lithosphere in continental China with its eflfcet focuing on the Eastern China A numerical research on the small-scale convection with variable viscosity in the upper mantle is introduced in Chapter 6. Based on a two-dimensional model, small-scale convection in the mantle-lithosphere system with variable viscosity is researched by using of finite element method. Variation of viscosity in exponential form with temperature is considered in this paper The results show that if viscosity is strongly temperature-dependent, the upper part of the system does not take a share in the convection and a stagnant lid, which is identified as lithosphere, is formed on the top of system because of low temperature and high viscosity. The calculated surface heat flow, topography and gravity anomaly are associated well with the convection pattern, namely, the regions with high heat flow and uplift correspond to the upwelling flow, and vice versa.In Chapter 7, we give a brief of future research subject: The inversion of lateral density heterogeneity in the mantle by minimizing the viscous dissipation.