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

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

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

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.

小波域地震资料的处理方法（WSSP）研究

The modeling formula based on seismic wavelet can well simulate zero - phase wavelet and hybrid-phase wavelet, and approximate maximal - phase and minimal - phase wavelet in a certain sense. The modeling wavelet can be used as wavelet function after suitable modification item added to meet some conditions. On the basis of the modified Morlet wavelet, the derivative wavelet function has been derived. As a basic wavelet, it can be sued for high resolution frequency - division processing and instantaneous feature extraction, in acoordance with the signal expanding characters in time and scale domains by each wavelet structured. Finally, an application example proves the effectiveness and reasonability of the method. Based on the analysis of SVD (Singular Value Decomposition) filter, by taking wavelet as basic wavelet and combining SVD filter and wavelet transform, a new de - noising method, which is Based on multi - dimension and multi-space de - noising method, is proposed. The implementation of this method is discussed the detail. Theoretical analysis and modeling show that the method has strong capacity of de - noising and keeping attributes of effective wave. It is a good tool for de - noising when the S/N ratio is poor. To give prominence to high frequency information of reflection event of important layer and to take account of other frequency information under processing seismic data, it is difficult for deconvolution filter to realize this goal. A filter from Fourier Transform has some problems for realizing the goal. In this paper, a new method is put forward, that is a method of processing seismic data in frequency division from wavelet transform and reconstruction. In ordinary seismic processing methods for resolution improvement, deconvolution operator has poor part characteristics, thus influencing the operator frequency. In wavelet transform, wavelet function has very good part characteristics. Frequency - division data processing in wavelet transform also brings quite good high resolution data, but it needs more time than deconvolution method does. On the basis of frequency - division processing method in wavelet domain, a new technique is put forward, which involves 1) designing filter operators equivalent to deconvolution operator in time and frequency domains in wavelet transform, 2) obtaining derivative wavelet function that is suitable to high - resolution seismic data processing, and 3) processing high resolution seismic data by deconvolution method in time domain. In the method of producing some instantaneous characteristic signals by using Hilbert transform, Hilbert transform is very sensitive to high - frequency random noise. As a result, even though there exist weak high - frequency noises in seismic signals, the obtained instantaneous characteristics of seismic signals may be still submerged by the noises. One method for having instantaneous characteristics of seismic signals in wavelet domain is put forward, which obtains directly the instantaneous characteristics of seismic signals by taking the characteristics of both the real part (real signals, namely seismic signals) and the imaginary part (the Hilbert transfom of real signals) of wavelet transform. The method has the functions of frequency division and noise removal. What is more, the weak wave whose frequency is lower than that of high - frequency random noise is retained in the obtained instantaneous characteristics of seismic signals, and the weak wave may be seen in instantaneous characteristic sections (such as instantaneous frequency, instantaneous phase and instantaneous amplitude). Impedance inversion is one of tools in the description of oil reservoir. one of methods in impedance inversion is Generalized Linear Inversion. This method has higher precision of inversion. But, this method is sensitive to noise of seismic data, so that error results are got. The description of oil reservoir in researching important geological layer, in order to give prominence to geological characteristics of the important layer, not only high frequency impedance to research thin sand layer, but other frequency impedance are needed. It is difficult for some impedance inversion method to realize the goal. Wavelet transform is very good in denoising and processing in frequency division. Therefore, in the paper, a method of impedance inversion is put forward based on wavelet transform, that is impedance inversion in frequency division from wavelet transform and reconstruction. in this paper, based on wavelet transform, methods of time - frequency analysis is given. Fanally, methods above are in application on real oil field - Sansan oil field.

济阳-昌潍坳陷火成岩油气藏体系形成机理和分布规律研究

Jiyang & Changwei depressions are two neighboring depressions in Bahai Bay Basin, the famous oil rich basin in East China. The exploration activities in the past 40 years has proved that, within the basins, there exists not only plentiful sandstone hydrocarbon reservoirs (conventional), but also abundant special reservoirs as igneous rock, mudstone and conglomerate ones which have been knowing as the unconventional in the past, and with the prospecting activity is getting more and more detailed, the unconventional reservoirs are also getting more and more important for further resources, among which, the igneous lithological reservoir be of significance as a new research and exploration area. The purpose of this paper is, with the historical researches and data as base, the System Theory, Practice Theory and Modern Comprehensive Petroleum Geology Theory as guide, the theoretical and practice break through as the goal, and the existing problems in the past as the break through direction, to explore and establish a valid reservoir formation and distribution models for igneous strata in the profile of the eastern faulted basins. After investigating the distribution of the igneous rocks and review the history of the igneous rocks reservoirs in basins, the author focused on the following issues and correspondingly the following progresses have been made: 1.Come to a new basin evolution and structure model named "Combined-Basin-bodies Model" for Jiyang even Eastern faulted basins based on the study on the origin and evolution of Jiyang & Changwei basins, depending on this model, every faulted basin in the Bo-hai Bay Basin is consisted of three Basin-Bodies including the Lower (Mesozoic), Middle (Early Tertiary) and the Upper (Late Tertiary) Bodies, each evolved in different geo-stress setting and with different basin trend, shape and igneous-sedimentary buildings system, and from this one to next one, the basin experienced a kind of process named "shape changing" and "Style changing". 2. Supposed a serious of new realizations as follows (1) There were "multi-level magma sources" including Upper mantel and the Lower, Middle and even the Upper Shell magma Chambers in the historical Magma Processes in the basins; (2) There were "multi-magma accessing or pass" from the first level (Mantel faults) to the second, third and fourth levels (that is the different levels of fault in the basin sediment strata) worked in the geo-historical and magma processes; (3) Three tectonic magma cycles and more periods have been recognized those are matched with the "Basin -body-Model" and (4)The geo-historical magma processes were non-homogeneous in time and space scale and so the magma rocks distributed in "zones" or "belts". 3. The study of magma process's effect on basin petroleum conditions have been made and the following new conclusions were reached: (1) the eruptive rocks were tend to be matched with the "caped source rock", and the magma process were favorable to the maturing of the source rocks. (2) The magma process were fruitful to the accumulation of the non-hydrocarbon reservoirs however a over magma process may damage the grade of resource rock; (3) Eruptive activity provided a fruitful environment for the formation of such new reservoir rocks as "co-eruptive turbidity sandstones" and "thermal water carbonate rocks" and the intrusive process can lead to the origin of "metamorphism rock reservoir"; (4) even if the intrusive process may cause the cap rock broken, the late Tertiary intrusive rocks may indeed provide the lateral seal and act as the cap rock locally even regionally. All above progresses are valuable for reconstructing the magma-sedimentary process history and enriching the theory system of modem petroleum geology. 4. A systematic classification system has been provided and the dominating factors for the origin and distribution of igneous rock reservoirs have been worked out based on the systematic case studies, which are as follows: (1) The classification is given based on multi-factors as the origin type, litho-phase, type of reservoir pore, reservoir ability etc., (2) Each type of reservoir was characterized in a detailed way; (3) There are 7 factors dominated the intrusive reservoir's characteristics including depth of intrusion, litho-facies of surrounding rocks, thickness of intrusive rock, intrusive facies, frequency and size of the working faults, shape and tectonic deformation of rock, erosion strength of the rock and the time of the intrusion ect., in the contrast, 4 factors are for eruptive rocks as volcanic facies, frequency and size of the working faults, strength of erosion and the thermal water processing. 5. Several new concept including "reservoir litho-facies", "composite-volcanic facies" and "reservoir system" ect. Were suggested, based on which the following models were established: (1) A seven reservoir belts model for a intrusive unit profile and further more, (2) a three layers cubic model consisted of three layer as "metamorphic roe layer", "marginal layer" and "the core"; (3) A five zones vertical reservoir sequence model consisted of five litho-facies named A, B, C, D and E for a original lava unit and furthermore three models respectively for a erosion, subsidence and faulted lava unit; (4) A composite volcanic face model for a lava cone or a composite cone that is consisted of three facies as "crater and nearby face", "middle slope" and "far slope", among which, the middle slope face is the most potential reservoir area and producible for oil & gas. 6. The concept of "igneous reservoir" was redefined as the igneous, and then a new concept of "igneous reservoir system" was supposed which means the reservoir system consisted of igneous and associated non-igneous reservoirs, with non-hydrocarbon reservoir included. 7. The origin and distribution of igneous reservoir system were probed and generalized for the exploration applications, and origin models of the main reservoir sub-systems have been established including those of igneous, related non-igneous and non-hydrocarbon. For intrusive rocks, two reservoir formation models have been suggested, one is called "Original or Primary Model", and the another one is "Secondary Model"; Similarly, the eruptive rock reservoirs were divided in three types including "Highly Produced", "Moderately Produced" and "Lowly Produced" and accordingly their formation models were given off; the related non-igneous reservoir system was considered combination of eight reservoirs, among which some ones like the Above Anticline Trap are highly produced; Also, the non-hydrocarbon. Trap system including five kinds of traps was discussed. 8. The concept models for four reservoir systems were suggested, which include the intrusive system consisted of 7 kinds of traps, the land eruptive system with 6 traps, the under water eruptive system including 6 kinds of traps and the non-hydrocarbon system combined by 5 kinds of traps. In this part, the techniques for exploration of igneous reservoir system were also generalized and probed, and based on which and the geological progresses of this paper, the potential resources and distributions of every reservoir system was evaluated and about 186 millions of reserves and eight most potential non-hydrocarbon areas were predicted and outlined. The author believe that the igneous reservoir system is a very important exploration area and its study is only in its early stage, the framework of this paper should be filled with more detailed studies, and only along way, the exploration of igneous reservoir system can go into it's really effective stage.

内蒙古草地植被群落生物量的梯度特征研究

Relationship between biology and environment is always the theme of ecology. Transect is becoming one of the important methods in studies on relationship between global change and terrestrial ecosystems, especially for analysis of its driving factors. Inner Mongolia Grassland is the most important in China Grassland Transect brought forward by Yu GR. In this study, changes in grassland community biomass along gradients of weather conditions in Inner Mongolia was researched by the method of transect. Methods of regression about biomass were also compared. The transect was set from Eerguna county to Alashan county (38° 07' 35" ~50° 12' 20" N, 101° 55' 25" -120° 20' 46" E) in Inner Mongolia, China. The sample sites were mainly chosen along the gradient of grassland type, meadow steppe-* typical steppe-*desert steppe-*steppification desert-^desert. The study was carried out when grassland community biomass got the peak in August or September, 2003 and 2004. And data of 49 sample sites was gotten, which included biomass, mean annual temperature, annual precipitation, accumulated temperature above zero, annual hours of sunshine and other statistical and descriptive data. The aboveground biomass was harvested, and the belowground biomass was obtained by coring (30 cm deep). Then all the biomass samples were dried within (80 + 5) °C in oven and weighted. The conclusion is as follows: 1) From the northeast to the southwest in Inner Mongolia, along the gradient of grassland type, meadow steppe-*typical steppe-*desert steppe-*steppification desert-* desert, the cover degree of vegetation community reduces. 2) By unitary regression analysis, biomass is negatively correlated with mean annual temperature, s^CTC accumulated temperature, ^10°C accumulated temperature and annual hours of sunshine, among which mean annual temperature is crucial, and positively with mean annual precipitation and mean annual relative humidity, and the correlation coefficient between biomass and mean annual relative humidity is higher. Altitude doesn't act on it evidently. Result of multiple regression analysis indicates that as the primary restrictive factor, precipitation affects biomass through complicated way on large scale, and its impaction is certainly important. Along the gradient of grassland type, total biomass reduces. The proportion of aboveground biomass to total biomass reduces and then increases after desert steppe. The trend of below ground biomass's proportion to total biomass is adverse to that of aboveground biomass. 3) Precipitation is not always the only driving factor along the transect for below-/aboveground biomass ratio of different vegetation type composed by different species, and distribution of temperature and precipitation is more important, which is much different among climatic regions, so that the trend of below-/aboveground biomass ratio along the grassland transect may change much through the circumscription of semiarid region and arid region. 4) Among reproductive allocation of aboveground biomass, only the proportion of stem in total biomass notably correlates to the given parameters. Stem/leaf biomass ratio decreases when longitude and latitude increase, caloric variables decrease, and variables about water increase from desert to meadow steppe. The change trends are good modeled by logarithm or binomial equations. 5) 0'-10 cm belowground biomass highly correlates to environmental parameters, whose proportion to total biomass changes most distinctly and increases along the gradient from the west to the east. The deeper belowground biomass responses to the environmental change on the adverse trend but not so sensitively as the surface layer. Because the change value of 0~10 cm belowground biomass is always more than that of below 10 cm along the gradient, the deference between them is balanced by aboveground biomass's change by the resource allocation equilibrium hypothesis.

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

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.