Surface pigments,algal biomass profiles,and potential production of the euphotic layer:Relationships reinvestigated in view of remote-sensing applications

Maps of surface chlorophyllous pigment (Chl a + Pheo a) are currently produced from ocean color sensors. Transforming such maps into maps of primary production can be reliably done only by using light-production models in conjuction with additional information about the column-integrated pigment content and its vertical distribution. As a preliminary effort in this direction. $\ticksim 4,000$ vertical profiles pigment (Chl a + Pheo a) determined only in oceanic Case 1 waters have been statistically analyzed. They were scaled according to dimensionless depths (actual depth divided by the depth of the euphotic layer, $Z_e$) and expressed as dimensionless concentrations (actual concentration divided by the mean concentration within the euphotic layer). The depth $Z_e$ generally unknown, was computed with a previously develop bio-optical model. Highly sifnificant relationships were found allowing $\langle C \rangle_tot$, the pigment content of the euphotic layer, to be inferred from the surface concentration, $\bar C_pd$, observed within the layer of one penetration depth. According to their $\bar C_pd$ values (ranging from $0.01 to > 10 mg m^-3$), we categorized the profiles into seven trophic situations and computed a mean vertical profile for each. Between a quasi-uniform profile in eutrophic waters and a profile with a strong deep maximum in oligotrophic waters, the shape evolves rather regularly. The wellmixed cold waters, essentially in the Antarctic zone, have been separately examined. On average, their profiles are featureless, without deep maxima, whatever their trophic state. Averaged values their profiles are featureless, without deep maxima, whatever their trophic state. Averaged values their profiles are featureless, without deep maxima, whatever their trophic state. Averaged values of $ρ$, the ratio of Chl a tp (Chl a + Pheo a), have also been obtained for each trophic category. The energy stored by photosynthesizing algae, once normalized with respect to the integrated chlorophyll biomass $\langle C \rangle _tot $ is proportional to the available photosythetic energy at the surface via a parameter $ψ∗$ which is the cross-section for photosynthesis per unit of areal chlorophyll. By tanking advantage of the relative stability of $ψ∗.$ we can compute primary production from ocean color data acquired from space. For such a computation, inputs are the irradiance field at the ocean surface, the "surface" pigment from which $\langle C \rangle _tot$ can be derived, the mean $ρ value pertinent to the trophic situation as depicted by the $\bar C_pd or $\langle C \rangle _tot$ values, and the cross-section $ψ∗$. Instead of a contant $ψ∗.$ value, the mean profiles can be used; they allow the climatological field of the $ψ∗.$ parameter to be adjusted through the parallel use of a spectral light-production model.

实行三个基本转变促进黄土丘陵区农业和经济持续发展

依据黄土丘陵区几十年治理开发的正反两方面经验 ,认为该地区要真正抓住西部大开发大好机遇 ,迅速改变落后面貌 ,实现农业和经济的快速及可持续发展 ,必须在指导思想和农业生产上实现三个根本转变。即 :变以粮为主的广种薄收为高效的农林果牧综合发展 ;变自给自足的小农经济为以商品生产为主的集约经营 ;变生态经济相悖论为生态经济统一观

Assessing the consequence of land use change on agricultural productivity in China

China's cultivated land has been undergoing dramatic changes along with its rapidly growing economy and population. The impacts of land use transformation on food production at the national scale, however, have been poorly understood due to the lack of detailed spatially explicit agricultural productivity information on cropland change and crop productivity. This study evaluates the effect of the cropland transformation on agricultural productivity by combining the land use data of China for the period of 1990-2000 from TM images and a satellite-based NPP (net primary production) model driven with NOAH/AVHRR data. The cropland area of China has a net increase of 2.79 Mha in the study period, which causes a slightly increased agricultural productivity (6.96 Mt C) at the national level. Although the newly cultivated lands compensated for the loss from urban expansion, but the contribution to production is insignificant because of the low productivity. The decrease in crop production resulting from urban expansion is about twice of that from abandonment of arable lands to forests and grasslands. The productivity of arable lands occupied by urban expansion was 80% higher than that of the newly cultivated lands in the regions with unfavorable natural conditions. Significance of cropland transformation impacts is spatially diverse with the differences in land use change intensity and land productivity across China. The increase in arable land area and yet decline in land quality may reduce the production potential and sustainability of China's agro-ecosystems. (C) 2008 Elsevier B.V. All rights reserved.

Tuning the Energy Level of Organic Sensitizers for High-Performance Dye-Sensitized Solar Cells

Cost-effective organic sensitizers will play a pivotal role in the future large-scale production and application of dye-sensitized solar cells. Here we report two new organic D-pi-A dyes featuring electron-rich 3,4-ethylenedioxythiophene- and 2,2'-bis(3,4-ethylenedioxythiophene)-conjugated linkers, showing a remarkable red-shifting of photocurrent action spectra compared with their thiophene and bithiophene counterparts. On the basis of the 3-f{5'-[N,N-bis(9,9-dimethylfluorene-2-yl)phenyl]-2,2'-bis(3,4-ethylenedioxythiophene)-5-yl}2-cyanoacrylic acid dye, we have set a new efficiency record of 7.6% for solvent-free dye-sensitized solar cells based on metal-free organic sensitizers. Importantly, the cell exhibits an excellent stability, keeping over 92% of its initial efficiency after 1000 h accelerated tests under full sunlight soaking at 60 degrees C. This achievement will considerably encourage further design and exploration of metal-free organic dyes for higher performance dye-sensitized solar cells.

Observation and theoretical analysis for the sand-waves migration in the North Gulf of South China Sea

According to the environmental characteristic of the north gulf of South China Sea, a quasi-3D mechanics model has been built for simulating the small scale sand-waves migration in the seas of southwest of Hainan Island. Based on the submarine micro-geomorphic data induced by multi-beam system and hydrographic survey record, the migrations of the sand-waves in the study area are predicted. The results show that calculation is consistent with the observation data in the groove of sand ridge, but not well in the crest of sand ridge. It is indicated that the mechanics model should be used to predict the migration of the small scale sand-waves which are dominated by bed load in the seas. This paper is very meaningful to project the route of submarine pipeline.

Robust Guaranteed Cost Control for Yaw Control of Helicopter

针对具有时变不确定性且不确定性界为椭球的线性系统提出了一种新的具有自适应机制的鲁棒保性能控制器设计方法。首先,引入一个具有可由自适应律在线调整的可调参数的目标模型,通过该参数来保证由目标模型与被控模型所获得的误差系统渐近稳定。结合保证目标模型稳定性的设计,最终形成保证闭环系统稳定且控制器增益仿射依赖于可调参数的鲁棒保性能跟踪控制器。应用于安装在试验平台上的小型直升机航向控制中,仿真试验表明了该方法的有效性。

模型直升机航向自适应保性能控制

提出了一种新的模型直升机航向控制算法。针对具有模型不确定性的直升机航向线性模型,提出了一种具有自适应机制的最优保性能控制器。该控制策略通过引入自适应机制降低固定增益控制器所固有的保守性,并且控制器的反馈增益应用线性矩阵不等式(LMIs)方法解得。理论分析和数字仿真表明所设计的控制器具有良好的鲁棒稳定性能。

柴达木盆地地震勘探关键技术研究

The rugged surface topography determined the seismic data acquisition construction conditions and the seismic wave explosive and receiver quality in Qaidam Basin. This dissertation systematically researched the seismic acquisition, imaging process and the attribute analysis techniques of complicated oil and gas reservoir. The main research achievements and cognitions are as follows: 1. Through the stimulation effects research and analysis from the aspect of lithologic water-containing differences, it’s specific that stable hydrous sand layer can effectively enhance the stimulation effects combined with the corresponding field tests. The seismic data S/N ratio has been improved due to the combination explosive stimulation. Through the fold number and maximum offset analyses of target horizon, the complicated geometry has been optimized and the S/N ratio of seismic data has been improved, which made an important basis for improvement of 3D seismic data. 2. It has been proved that the first arrival refraction static correction method under the model constraint of fine surface survey is suitable to the Qaidam Basin of western areas by the real seismic data processing. Although the refraction horizon of near surface has some changes in a certain extent, it’s steady basically. The refraction horizon can be continuously traced in sections, so it’s qualified for the refraction static correction method on the whole. 3. The research is based on the curved-ray pre-stack time migration techniques of rough topography, and improved the imaging precision of complex areas. This techniques adopted the constant and variable velocity scanning mode and enhanced the velocity analysis precision. The 3D pre-stack time migration techniques reasonably solved the imaging and velocity multiple solutions problems of steep-dip faults and the intersections of horizontal layers. What’s more, fine velocity analysis and mute are very important to enhance the imaging precision of the seismic data in complicated Wunan areas. 4. The 3D seismic data edge-preserving processing methods have been realized due to the image process techniques. Because this method uses the large range filter, it can attenuate the noise maximally. The faults, break points, lithologic pinchout points and lithologic body of small scale such as river will not be influenced by blur because of the edge-preserving characterization of the method which is really an effective assistant technique of low S/N ratio seismic data attribute analysis. 5. The use of spectral decomposition technique can effectively identify the reservoirs. The special geology body which will not be identified (or without obvious characters) in the seismic profile may be found through the details changes of different frequencies in the amplitude profiles.

地震岩性与物性联合反演技术研究和应用

This paper studies how to more effectively invert seismic data and predict reservoir under complicated sedimentary environment, complex rock physical relationships and fewer drills in offshore areas of China. Based on rock physical and seismic amplitude-preserving process, and according to depositional system and laws of hydrocarbon reservoir, in the light of feature of seismic inversion methods present applied, series methods were studied. A joint inversion technology for complex geological condition had been presented, at the same time the process and method system for reservoir prediction had been established. This method consists four key parts. 1)We presented the new conception called generalized wave impedance, established corresponding inversion process, and provided technical means for joint inversion lithology and petrophysical on complex geological condition. 2)At the aspect of high-resolution nonlinear seismic wave impedance joint inversion, this method used a multistage nonlinear seismic convolution model rather than conventional primary structure Robinson seismic convolution model, and used Caianiello neural network implement inversion. Based on the definition of multistage positive and negative wavelet, it adopted both deterministic and statistical physical mechanism, direct inversion and indirect inversion. It integrated geological knowledge, rock physical theory, well data, and seismic data, and improved the resolution and anti-noise ability of wave impedence inversion. 3)At the aspect of high-resolution nonlinear reservoir physical property joint inversion, this method used nonlinear rock physical model which introduced convolution model into the relationship between wave impedance and porosity/clay. Through multistage decomposition, it handles separately the large- and small-scale components of the impedance-porosity/clay relationships to achieve more accurate rock physical relationships. By means of bidirectional edge detection with wavelets, it uses the Caianiello neural network to finish statistical inversion with combined applications of model-based and deconvolution-based methods. The resulted joint inversion scheme can integrate seismic data, well data, rock physical theory, and geological knowledge for estimation of high-resolution petrophysical parameters. 4)At the aspect of risk assessment of lateral reservoir prediction, this method integrated the seismic lithology identification, petrophysical prediction, multi-scale decomposition of petrophysical parameters, P- and H-spectra, and the match relationship of data got from seismics, well logging and geology. It could describe the complexity of medium preferably. Through applications of the joint inversion of seismic data for lithologic and petrophysical parameters in several selected target areas, the resulted high-resolution lithologic and petrophysical sections(impedance, porosity, clay) show that the joint inversion can significantly improve the spatial description of reservoirs in data sets involving complex deposits. It proved the validity and practicality of this method adequately.

金川镍矿岩体移动与开采稳定性研究

The two major issues in mining industry are work safety and protection of ground environment when carrying on underground mining activities. Cut-and-fill mining method is increasingly applied in China owing to its advantages of controlling ground pressure and protecting the ground environment effectively. However, some cut-and-fill mines such as Jinchuan nickel mine which has big ore body, broken rock mass and high geostress have unique characteristics on the law of ground pressure and rock mass movement that distinguish from other mining methods. There are still many problems unknown and it is necessary for the further analysis. In this dissertation, vast field survey, geology trenching and relative data analysis are carried out. The distribution of ground fissures and the correlation of the fissures with the location of underground ore body is presented. Using of monitoring data by three-dimension fissure meter and GPS in Jinchuan Deposit Ⅱ, the rule of the surface deformation and the reason of ground fissures generation are analyzed. It is shown that the stress redistribution in surrounding rocks resulting from mining, the existence of the void space underground and the influence of on-going mining activities are three main reasons for the occurrence of ground fissures. Based on actual section planes of No.1 ore body, a large-scale 3D model is established. By this model, the complete process of excavation and filling is simulated and the law of rock mass movement and stability caused by Cut-and-fill Mining is studied. According to simulation results, it is concluded that the deformation of ground surface is still going on developing; the region of subsidence on the ground surface is similar with a circle; the area on the hanging wall side is larger than one on the lower wall side; the contour plots show the centre of subsidence lay on the hanging wall side and the position is near the ore body boundary of 1150m and 1250m where ore body is the thickest. Along strike-line of Jinchuan Deposit Ⅱ, the deformation at the middle of filling body is larger than that in the two sides. Because of the irregular ore body, stress concentrates at the boundary of ore body. With the process of excavation and filling, the high stress release and the stress focus disappear on the hanging wall side. The cut-and-fill mechanism is studied based on monitoring data and numerical simulation. The functions of filling body are discussed. In this dissertation, it is concluded that the stress of filling body is just 2MPa, but the stress of surrounding rock mass is 20MPa. We study the surface movement influenced by the elastic modulus of backfill. The minimal value of the elastic modulus of backfill which can guarantee the safety production of cut-and-fill mine is obtained. Finally, based on the real survey results of the horizontal ore layer and numerical simulation, it is indicated that the horizontal ore layer has destroyed. Key words: cut-and-filling mining, 3D numerical simulation, field monitoring, rock mass movement, cut-and-filling mechanism, the elastic modulus of backfill, the horizontal ore layer

油气勘探中的弱磁异常提取与区域场分析

Petroleum and natural gas is an important strategic resources. The short of the reserves will block the development of economy and threaten the safety of nation, along with the main oil fields of our country coming to the height of power and splendor of the exploitation and exploration. Therefore, it makes a great sense to inaugurate new explorative field and increase the reserves of petroleum and natural gas. Magnetic exploration is a main method of geophysics exploration. the developing observation apparatus and the perfect processing method provide wide space for magnetic exploration in these years. The method of magnetic bright spot is an application of magnetic exploration. The vertical migration of the hydrocarbon changes physical and chemical environment above the hydrocarbon reservoir, the new environment make tervalent iron translate into bivalent iron, that produce small scale magnetic anomaly, that is magnetic bright spot. The method of magnetic bright spot explores oil and gas field by the relation between the hydrocarbon and magnetic anomaly. This paper systemically research to pick-up and identify magnetic bright spot combining an oil field item, then point out advantaged area. In order to test the result, the author use the seismic information to superpose the magnetic bright spot, that prove the magnetic bright spot is reliable. then, the author complete a software to pick and identify the magnetic bright spot. The magnetic basement is very important to research forming and evolvement of the basin, especially, it is a crucial parameter of exploring residual basin in the research on pre-Cenozoic residual. This paper put forward a new method to inverse the interface of the magnetic layer on the basis of previous work, that is the method of separation of magnetic field step by step. The theory of this method is to translate the result of magnetic layer fluctuation to the result of magnetization density change, and the magnetic layer is flat, the paper choose thickness of magnetic layer as unit thickness, and define magnetic layer as a unit-thickness layer in order to convenient calculation, at the same time, define the variational magnetization density as equivalent magnetic density. Then we translate the relation between magnetic field and layer fluctuation to the relation between magnetic field and equivalent magnetic density, then, we can obtain the layer fluctuation through calculating equivalent magnetic density. Contrast to conventional parker method, model experimentation and example checkout prove this method is effective. The merit of this method is to avoid flat result in a strongly fluctuant area because of using a uniform average depth, the result of this method is closer to the fact, and this method is to inverse equivalent magnetic density, then translate equivalent magnetic density to layer fluctuation, this lays a foundation to inverse variational magnetic density in the landscape orientation and portrait.

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

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

复杂介质边界元法数值模拟研究

A major impetus to study the rough surface and complex structure in near surface model is because accuracy of seismic observation and geophysical prospecting can be improved. Wave theory study about fluid-satuated porous media has important significance for some scientific problems, such as explore underground resources, study of earth's internal structure, and structure response of multi-phase porous soil under dynamic and seismic effect. Seismic wave numerical modeling is one of the effective methods which understand seismic propagation rules in complex media. As a numerical simulation method, boundary element methods had been widely used in seismic wave field study. This paper mainly studies randomly rough surface scattering which used some approximation solutions based on boundary element method. In addition, I developed a boundary element solution for fluid saturated porous media. In this paper, we used boundary element methods which based on integral expression of wave equation to study the free rough surface scattering effects of Kirchhoff approximation method, Perturbation approximation method, Rytov approximation method and Born series approximation method. Gaussian spectrum model of randomly rough surfaces was chosen as the benchmark model. The approximation methods result were compared with exact results which obtained by boundary element methods, we study that the above approximation methods were applicable how rough surfaces and it is founded that this depends on and ( here is the wavenumber of the incident field, is the RMS height and is the surface correlation length ). In general, Kirchhoff approximation which ignores multiple scatterings between any two surface points has been considered valid for the large-scale roughness components. Perturbation theory based on Taylor series expansion is valid for the small-scale roughness components, as and are .Tests with the Gaussian topographies show that the Rytov approximation methods improves the Kirchhoff approximation in both amplitude and phase but at the cost of an extra treatment of transformation for the wave fields. The realistic methods for the multiscale surfaces come with the Born series approximation and the second-order Born series approximation might be sufficient to guarantee the accuracy of randomly rough surfaces. It could be an appropriate choice that a complex rough surface can be divided into large-, medium-, and small-scale roughness components with their scattering features be studied by the Kirchhoff or Rytov phase approximations, the Born series approximation, and the perturbation theory, respectively. For this purpose, it is important to select appropriate parameters that separate these different scale roughness components to guarantee the divided surfaces satisfy the physical assumptions of the used approximations, respectively. In addition, in this paper, the boundary element methods are used for solving the porous elastic wave propagation and carry out the numerical simulation. Based on the fluid-saturated porous model, this paper analyses and presents the dynamic equation of elastic wave propagation and boundary integral equation formulation of fluid saturated porous media in frequency domain. The fundamental solutions of the elastic wave equations are obtained according to the similarity between thermoelasticity and poroelasticity. At last, the numerical simulation of the elastic wave propagation in the two-phase isotropic media is carried out by using the boundary element method. The results show that a slow quasi P-wave can be seen in both solid and fluid wave-field synthetic seismograms. The boundary element method is effective and feasible.

特大地震激发地球自由振荡的数值模拟研究

Large earthquakes, such as the Chile earthquake in 1960 and the Sumatra-Andaman earthquake on Dec 26, 2004 in Indonesia, have generated the Earth’s free oscillations. The eigenfrequencies of the Earth’s free oscillations are closely related to the Earth’s internal structures. The conventional methods, which mainly focus on calculating the eigenfrequecies by analytical ways, and the analysis on observations can not easily study the whole processes from earthquake occurrence to the Earth’s free oscillation inspired. Therefore, we try to use numerical method incorporated with large-scale parallel computing to study on the Earth’s free oscillations excited by giant earthquakes. We first give a review of researches and developments of the Earth’s free oscillation, and basical theories under spherical coordinate system. We then give a review of the numerical simulation of seismic wave propagation and basical theories of spectral element method to simulate global seismic wave propagation. As a first step to study the Earth’s free oscillations, we use a finite element method to simulate the propagation of elastic waves and the generation of oscillations of the chime bell of Marquis Yi of Zeng, by striking different parts of the bell, which possesses the oval crosssection. The bronze chime bells of Marquis Yi of Zeng are precious cultural relics of China. The bells have a two-tone acoustic characteristic, i.e., striking different parts of the bell generates different tones. By analysis of the vibration in the bell and the spectrum analysis, we further help the understanding of the mechanism of two-tone acoustic characteristics of the chime bell of Marquis Yi of Zeng. The preliminary calculations have clearly shown that two different modes of oscillation can be generated by striking different parts of the bell, and indicate that finite element numerical simulation of the processes of wave propagation and two-tone generation of the chime bell of Marquis Yi of Zeng is feasible. These analyses provide a new quantitative and visual way to explain the mystery of the two-tone acoustic characteristics. The method suggested by this study can be applied to simulate free oscillations excited by great earthquakes with complex Earth structure. Taking into account of such large-scale structure of the Earth, small-scale low-precision numerical simulation can not simply meet the requirement. The increasing capacity in high-performance parallel computing and progress on fully numerical solutions for seismic wave fields in realistic three-dimensional spherical models, Spectral element method and high-performance parallel computing were incorporated to simulate the seismic wave propagation processes in the Earth’s interior, without the effects of the Earth’s gravitational potential. The numerical simulation shows that, the results of the toroidal modes of our calculation agree well with the theoretical values, although the accuracy of our results is much limited, the calculated peaks are little distorted due to three-dimensional effects. There exist much great differences between our calculated values of spheroidal modes and theoretical values, because we don’t consider the effect the Earth’ gravitation in numerical model, which leads our values are smaller than the theoretical values. When , is much smaller, the effect of the Earth’s gravitation make the periods of spheroidal modes become shorter. However, we now can not consider effects of the Earth’s gravitational potential into the numerical model to simulate the spheroidal oscillations, but those results still demonstrate that, the numerical simulation of the Earth’s free oscillation is very feasible. We make the numerical simulation on processes of the Earth’s free oscillations under spherically symmetric Earth model using different special source mechanisms. The results quantitatively show that Earth’s free oscillations excited by different earthquakes are different, and oscillations at different locations are different for free oscillation excited by the same earthquake. We also explore how the Earth’s medium attenuation will take effects on the Earth’s free oscillations, and take comparisons with the observations. The medium attenuation can make influences on the Earth’s free oscillations, though the effects on lower-frequency fundamental oscillations are weak. At last, taking 2008 Wenchuan earthquake for example, we employ spectral element method incorporated with large-scale parallel computing technology to investigate the characteristics of seismic wave propagation excited by Wenchuan earthquake. We calculate synthetic seismograms with one-point source model and three-point source model respectively. Full 3-D visualization of the numerical results displays the profile of the seismic wave propagation with respect to time. The three-point source, which was proposed by the latest investigations through field observation and reverse estimation, can better demonstrate the spatial and temporal characteristics of the source rupture processes than one-point source. Primary results show that those synthetic signals calculated from three-point source agree well with the observations. This can further reveal that the source rupturing process of Wenchuan earthquake is a multi-rupture process, which is composed by at least three or more stages of rupture processes. In conclusion, the numerical simulation can not only solve some problems concluding the Earth’s ellipticity and anisotropy, which can be easily solved by conventional methods, but also finally solve the problems concluding topography model and lateral heterogeneity. We will try to find a way to fully implement self-gravitation in spectral element method in future, and do our best to continue researching the Earth’s free oscillations using the numerical simulations to see how the Earth’ lateral heterogeneous will affect the Earth’s free oscillations. These will make it possible to bring modal spectral data increasingly to bear on furthering our understanding of the Earth’s three-dimensional structure.

南海北部大陆边缘盆地深部地热特征用地球动力学演化

This report is a conclusion of the major research outcome during my post-doctoral residence of research and work. Its content covers the researches of the deep thermal characteristic and dynamics evolution beneath the northern margin basin of South China Sea. In this report, the each other action and effect between lithosphere ad mantle convection were regarded by the combine of deep and shallow study, subdivision from whole to part, and pay equal attention to determine the nature and fixed quantity. The investigative method we used in this report is geothermal and gravity methods. By the help of geological model and geophysics modeling, we calculated lithosphere thermal structure, rheology structure and mantle convection. Firstly, the report introduces concisely the purpose and the previous achievement to this research. Then, it analyzed the characteristic of heat flow on South China Sea. The structure of deep temperature and thermal has been calculated in some models of heat generation and conduction. The rock rheology structure also was computed by the relationship between temperature and viscosity. All these calculations were finished under the guidelines of combine with geology and geophysics. Meanwhile, the fields both deep mantle convection and small scale upper mantle convection are computed. Beside, the density and temperature disorder resulted by mantle convection were also computed with the convection field. After these, the report bring the contribution of local field of mantle convection, thermal construct and effective viscosity beneath the northern margin basin of South China Sea. And, base on the tectonic background and evolution feature, this report discussion the evolution mechanism of south China Sea and its northern margin basin. The end of this report, the main conclusion of this research was summarized and brings out.