中国东北温带样带区域和景观尺度生态模型的研究

样带是沿全球变化某一驱动因素的主要梯度而设置的由一系列研究站点构成的区域，被认为是研究全球变化与陆地生态系统关系的最有效的途径。而模型研究是全球变化研究中不可或缺的手段。本文即采用模型研究方法研究中国东北温带样带(NECT)区域，试图揭示温带生态系统对于全球变化（尤其是降水）的反应机制。 中国东北温带样带（NECT)位于42°N - 46°N，108°E - 132°E，长约二千多公里，是最早被列入GCTE的四条样带之一，从东到西有明显的湿度梯度，被认为是温带区域研究水分梯度的代表性样带。本文研究主要集中在： 1．NECT中环境数据库的建立,本文采用EIS作为数据管理系统。由于EIS管理空间数据的特点是根据确定的地理坐标来提供空间定位，因而每一环境因子的属性值分布都有确定的地理坐标与其对应，特别适合于样带这种研究区域较大，同时又要求有精确空间定位的区域。NECT环境数据库包括地形、气候、植被、土壤、土地利用、水文、孢粉数据及社会经济等分库、本数据库力图提供各环境因子的各种属性值而代替仅仅提供类型值。 2．NECT中PFTs的划分PFTs的划分被认为是建立DGVM的前提。本文认为PFTs的划分是模型研究中一个尺度上升过程的结果，不同的尺度，不同的研究目标导致不同的PFTs的划分。在NECT区域中，考虑植被对全球变化中降水因子的不同反映机制，采用生活型、高度、耐旱特性、叶子大小、叶子季相、主根深度和木质化程度等指标根据- TWINSPAN和FCLUS进行划分，得到以下9种NECT区域中植被功能类型：常绿针叶树种、落叶针叶树种、落叶阔叶树种、落叶小叶灌木、落叶小叶半灌木、落叶强旱生半灌木、多年生中旱生草本、适应旱生环境的多年生草本和多年生强旱生草本。对NECT区域中PFTs的DCA分析表明降水是控制PFTs在NECT区域中分布的主要环境因子。在代表景观层次的长白山PFTs的划分中，则采用树种有记载的最大寿命、最大胸径、最大树高、各树种生长参数、树种自然分布区内>5℃的有效积温的最小值和最大值、耐阴、耐旱、喜肥特性、树种的扩散更新，就地下种更新和萌条更新能力参数及叶子大小和类型等指标采用上述软件得到的以下PFTs:即不耐荫阔叶树种、耐荫阔叶树种、耐荫针叶树种和不耐荫的阳生针叶树种。 3.NECT中BCM模型的建立和预测 本文认为土壤水是决定SPC系统水分状况的直接指标。而均衡土壤水分剖面代表了土壤水的多年平均状态，因而本文以Watershed模型为基础，模拟NECT区域中任意一点的均衡土壤水分剖面（精度为每经纬网格32×48个点）;然后根据这个均衡水分剖面用计算LAI子模型确定该水分剖面所能支持的LAI;进而根据这个LAI由Biome等模型划分出Biome在NECT中的分布。全球变化的结果将改变区域中任意一点的土壤水分状况，从而影响植被的LAI，进而导致Biome的改变。本模型成功的模拟了LAI和Biome在NECT区域中的分布，利用85-90生长季每月平均的NDVI作相关检验表明除5月份以外，相关系数都>0.7，而5月份也达到0,6457，都达到了极显著的程度。尤为重要的是，模型对于不同植被类型的NDVI与LAI的对应关系有良好的模拟，如针叶林的LAI在相同的NDVI值下明显比阔叶林小，因而模型模拟的LAI在NECT东部针叶林分布区LAI值比针阔混交林明显偏小，而与Spanmen等(1990)所提出的针叶林叶面积指数与NDVI关系非常一致。模型的预测显示：(1) T+20C (PET+15%)，Precipitation+20%，LAI总体上变化不大，且空间变化呈现复杂性，总体上表现出草原植被LAI减少，而森林的LAI增加;Biome层次表现出针阔混交林和矮草原面积扩大，针叶林和森林草原面积减少，其中对于该情形下变化最为明显的是针叶林和森林草原。NECT东部区域发育在沙性土上的植被的LAI明显增加，而科尔沁沙地植被的LAI则维持不变。(2)T+40C (PET+30%)，Precipitation+20%，LAI总体上将减小0.14，但空间分布不均。东部森林区域LA1将维持不变或增加（主要为针叶林），草原植被LAI仍表现出减少趋势;在Biome层次上则表现出草原面积的扩大。对于第一种情形下LAI有增有减的森林革原地区则表现出减小的一致性，总体来说，第二种情形比第一种情形表现出相当的干旱性。从对两种全球变化情形的反应来看，针叶林和森林草原是NECT中对全球变化驱动因子温度和降水的敏感植被类型;丽科尔沁沙地植被表现出相当的稳定性，表明该沙地的敏感性主要是由于人类活动这个因子造成的。 4.NECT中景观层次NPP模型的建立和预测 景观层次之所以成为模型研究中一个独特的层次，是由于地形效应的存在。地形效应对于水、热。营养物质的进行重新分配，从而进一步控制了生态系统的分布。本文选择NECT区域中森林生态系统的代表性分布小流域一二道白河小流域为研究区域。首先，应用Sunlight模型来模拟小流域任意一点所截取的能用于光合作用的太阳辐射能。Sunlight模型充分考虑了由于栅格的坡度、坡向和遮蔽度对可照时间和太阳直射辐射的影响以及坡度和可祝度对太阳散射辐射的影响，并提供了消除大气状况从站点观察数据推测的方法，即太阳直射辐射转换系数Rb和太阳散射辐射转换系数R，结合植被的分布特性，得到IPR在小流域中的分布。结果表明，IPR在小流域中相差不大，与高程呈正相关。进而利用温度修正模型得到温度修正系数，平均为0.446，表明温度对NPP的限制效应比较大;而水分修正系数则通过Topmodel模拟每一栅格的地下水位，由这个地下水位通过前述Waterbalance模型模拟均衡土壤水分剖面，进而求出水分修正系数，平均为0.86，表明该流域水分状况良好，水分状况对NPP的限制性不强。模拟结果显示：海拔1700m的岳桦林平均NPP为7.31t.ha-l年“两实测的NPP值分别为20.19 t.ha-l年 13.45 t.ha-l年‘1和5.15 t.ha-1年”。可见模型模拟的整个NPP的分布趋势与现实NPP分布一致。模型的预测表明：(1)．T+2℃( PET+15%)，P+20%时整个小流域的平均NPP将从14.60 t.ha-1年-1增加到16.65 t.ha_l年，且不同生态系统的增加幅度大致类似。以岳桦林生态系统稍高，增加幅度为15.3%，引起该小流域NPP增加的主要原因是温度状况的改善，而水分状况与变化前一致。(2)T+4℃( PET+30%)，P+20%的变化前景下NPP的增加幅度更大，整个小流域的平均NPP将达到18.63 t.ha-1年，且不同生态系统变化幅度不同．阔叶红松林、暗针叶林和岳桦林生态系统NPP增加幅度将分别达到25.8%，28.3%和32.8%.引起整个小流域NPP增加的主要原因仍旧是温度状况的大幅度提高，平均从变化前的0.477提高到0.610，而水分状况则仍保持相当的稳定性。

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

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.

整合多尺度地球物理数据建立高分辨率储层模型的方法及应用

Most of the fields in China are in the middle-late development phase or are mature fields. It becomes more and more difficult to develop the remaining oil/gas. Therefore, it is import to enhance oil/gas recovery in order to maintain the production. Fine scale modeling is a key to improve the recovery. Incorporation of geological, seismic and well log data to 3D earth modeling is essential to build such models. In Ken71 field, well log, cross-well seismic and 3D seismic data are available. A key issue is to build 3D earth model with these multi-scales data for oil field development.In this dissertation, studies on sequential Gaussian-Bayesian simulation have been conducted. Its comparison with cokriging and sequential Gaussian simulation has been performed. The realizations generated by sequential Gaussian-Bayesian simulation have higher vertical resolution than those generated by other methods. Less differences between these realization and true case are observed. With field data, it is proved that incorporating well log, cross-well seismic and 3D seismic into 3D fine scale model is reliable. In addition, the advantages of sequential Gaussian-Bayesian simulation and conditions for input data are demonstrated. In Ken71 field, the impedance difference between sandstone and shale is small. It would be difficult to identify sandstone in the reservoir with traditional impedance inversion. After comparisons of different inversion techniques, stochastic hillclimbing inversion was applied. With this method, shale content inversion is performed using 3D seismic data. Then, the inverted results of shale content and well log data are incorporated into 3D models. This demonstrates a procedure to build fine scale models using multi scale seismic data, especially 3D seismic amplitude volume.The models generated through sequential Gaussian-Bayesian simulation have several advantages including: (1) higher vertical resolution compared with 3D inverted acoustic impedance (AI); (2) consistency of lateral variation as 3D inverted AI; (3) more reliability due to integration cross-well seismic data. It is observed that the precision of the model depends on the 3D inversion.

Comparisons of static, quasi-static and dynamic 3D porous media scale network models for two-phase immiscible flow in porous media

A dynamic 3D pore-scale network model is formulated for investigating the effect of interfacial tension and oil-water viscosity during chemical flooding. The model takes into account both viscous and capillary forces in analyzing the impact of chemical properties on flow behavior or displacement configuration, while the static model with conventional invasion percolation algorithm incorporates the capillary pressure only. From comparisons of simulation results from these models. it indicates that the static pore scale network model can be used successfully when the capillary number is low. With the capillary increases due to the enhancement of water viscosity or decrease of interfacial tension, only the quasi-static and dynamic model can give insight into the displacement mechanisms.

Effects of Subgrid-Scale Modeling on Time Correlations in Large Eddy Simulation

The effects of the unresolved subgrid-scale (SGS) motions on the energy balance of the resolved scales in large eddy simulation (LES) have been investigated actively because modeling the energy transfer between the resolved and unresolved scales is crucial to constructing accurate SGS models. But the subgrid scales not only modify the energy balance, they also contribute to temporal decorrelation of the resolved scales. The importance of this effect in applications including the predictability problem and the evaluation of sound radiation by turbulent flows motivates the present study of the effect of SGS modeling on turbulent time correlations. This paper compares the two-point, two-time Eulerian velocity correlation in isotropic homogeneous turbulence evaluated by direct numerical simulation (DNS) with the correlations evaluated by LES using a standard spectral eddy viscosity. It proves convenient to express the two-point correlations in terms of spatial Fourier decomposition of the velocity field. The LES fields are more coherent than the DNS fields: their time correlations decay more slowly at all resolved scales of motion and both their integral scales and microscales are larger than those of the DNS field. Filtering alone is not responsible for this effect: in the Fourier representation, the time correlations of the filtered DNS field are identical to those of the DNS field itself. The possibility of modeling the decorrelating effects of the unresolved scales of motion by including a random force in the model is briefly discussed. The results could have applications to the problem of computing sound sources in isotropic homogeneous turbulence by LES

Orientation-based Continuum Damage Models for Rocks

A general formulation of the Helmholtz free energy used in thermodynamics of damage process of rocks is derived within a multi-scale framework. Such a physically-based thermodynamic state potential has a hybrid, discrete/continuum, nature in the sense tha

Micro-scale Mechanics of the Surface-Nanocrystalline Al-Alloy Material

Based on the microscopic observations and measurements, the mechanical behavior of the surface-nanocrystallized Al-alloy material at microscale is investigated experimentally and theoretically. In the experimental research, the compressive stress-strain curves and the hardness depth curves are measured. In the theoretical simulation, based on the material microstructure characteristics and the experimental features of the compression and indentation, the microstructure cell models are developed and the strain gradient plasticity theory is adopted. The material compressive stress-strain curves and the hardness depth curves-are predicted and simulated. Through comparison of the experimental results with the simulation results, the material and model parameters are determined.

Toughness Of Ni/Al2O3 Interfaces As Dependent On Micron-Scale Plasticity And Atomistic-Scale Separation

Ceramic/metal interfaces were studied that fail by atomistic separation accompanied by plastic dissipation in the metal. The macroscopic toughness of the specific Ni alloy/Al2O3 interface considered is typically on the order of ten times the atomistic work of separation in mode I and even higher if combinations of mode I and mode II act on the interface. Inputs to the computational model of interface toughness are: (i) strain gradient plasticity applied to the Ni alloy with a length parameter determined by an indentation test, and (ii) a potential characterizing mixed mode separation of the interface fit to atomistic results. The roles of the several length parameters in the strain gradient plasticity are determined for indentation and crack growth. One of the parameters is shown to be of dominant importance, thus establishing that indentation can be used to measure the relevant length parameter. Recent results for separation of Ni/Al2O3 interfaces computed by atomistic methods are reviewed, including a set of results computed for mixed mode separation. An approximate potential fit to these results is characterized by the work of separation, the peak separation stress for normal separation and the traction-displacement relation in pure shearing of the interface. With these inputs, the model for steady-state crack growth is used to compute the toughness of the interface under mode I and under the full range of mode mix. The effect of interface strength and the work of separation on macroscopic toughness is computed. Fundamental implications for plasticity-enhanced toughness emerge.

Peeling Experiments Of Ductile Thin Films Along Ceramic Substrates - Critical Assessment Of Analytical Models

Two types of peeling experiments are performed in the present research. One is for the Al film/Al2O3 substrate system with an adhesive layer between the film and the substrate. The other one is for the Cu film/Al2O3 substrate system without adhesive layer between the film and the substrate, and the Cu films are electroplated onto the Al2O3 substrates. For the case with adhesive layer, two kinds of adhesives are selected, which are all the mixtures of epoxy and polyimide with mass ratios 1:1.5 and 1:1, respectively. The relationships between energy release rate, the film thickness and the adhesive layer thickness are measured during the steady-state peeling process. The effects of the adhesive layer on the energy release rate are analyzed. Using the experimental results, several analytical criteria for the steady-state peeling based on the bending model and on the two-dimensional finite element analysis model are critically assessed. Through assessment of analytical models, we find that the cohesive zone criterion based on the beam bend model is suitable for a weak interface strength case and it describes a macroscale fracture process zone case, while the two-dimensional finite element model is effective to both the strong interface and weak interface, and it describes a small-scale fracture process zone case. (C) 2007 Elsevier Ltd. All rights reserved.

Effects Of Subgrid-Scale Modeling On Lagrangian Statistics In Large-Eddy Simulation

The application of large-eddy simulation (LES) to turbulent transport processes requires accurate prediction of the Lagrangian statistics of flow fields. However, in most existing SGS models, no explicit consideration is given to Lagrangian statistics. In this paper, we focus on the effects of SGS modeling on Lagrangian statistics in LES ranging from statistics determining single-particle dispersion to those of pair dispersion and multiparticle dispersion. Lagrangian statistics in homogeneous isotropic turbulence are extracted from direct numerical simulation (DNS) and the LES with a spectral eddy-viscosity model. For the case of longtime single-particle dispersion, it is shown that, compared to DNS, LES overpredicts the time scale of the Lagrangian velocity correlation but underpredicts the Lagrangian velocity fluctuation. These two effects tend to cancel one another leading to an accurate prediction of the longtime turbulent dispersion coefficient. Unlike the single-particle dispersion, LES tends to underestimate significantly the rate of relative dispersion of particle pairs and multiple-particles, when initial separation distances are less than the minimum resolved scale due to the lack of subgrid fluctuations. The overprediction of LES on the time scale of the Lagrangian velocity correlation is further confirmed by a theoretical analysis using a turbulence closure theory.

A second-order dynamic subgrid-scale stress model

A second-order dynamic model based on the general relation between the subgrid-scale stress and the velocity gradient tensors was proposed. A priori test of the second-order model was made using moderate resolution direct numerical simulation date at high Reynolds number ( Taylor microscale Reynolds number R-lambda = 102 similar to 216) for homogeneous, isotropic forced flow, decaying flow, and homogeneous rotating flow. Numerical testing shows that the second-order dynamic model significantly improves the correlation coefficient when compared to the first-order dynamic models.

Statistical analysis of synthetic earthquake catalogs generated by models with various levels of fault zone disorder

The stress release model, a stochastic version of the elastic rebound theory, is applied to the large events from four synthetic earthquake catalogs generated by models with various levels of disorder in distribution of fault zone strength (Ben-Zion, 1996) They include models with uniform properties (U), a Parkfield-type asperity (A), fractal brittle properties (F), and multi-size-scale heterogeneities (M). The results show that the degree of regularity or predictability in the assumed fault properties, based on both the Akaike information criterion and simulations, follows the order U, F, A, and M, which is in good agreement with that obtained by pattern recognition techniques applied to the full set of synthetic data. Data simulated from the best fitting stress release models reproduce, both visually and in distributional terms, the main features of the original catalogs. The differences in character and the quality of prediction between the four cases are shown to be dependent on two main aspects: the parameter controlling the sensitivity to departures from the mean stress level and the frequency-magnitude distribution, which differs substantially between the four cases. In particular, it is shown that the predictability of the data is strongly affected by the form of frequency-magnitude distribution, being greatly reduced if a pure Gutenburg-Richter form is assumed to hold out to high magnitudes.

Subgrid-scale contributions to Lagrangian time correlations in isotropic turbulence

The application of large-eddy simulation (LES) to particle-laden turbulence raises such a fundamental question as whether the LES with a subgrid scale (SGS) model can correctly predict Lagrangian time correlations (LTCs). Most of the currently existing SGS models are constructed based on the energy budget equations. Therefore, they are able to correctly predict energy spectra, but they may not ensure the correct prediction on the LTCs. Previous researches investigated the effect of the SGS modeling on the Eulerian time correlations. This paper is devoted to study the LTCs in LES. A direct numerical simulation (DNS) and the LES with a spectral eddy viscosity model are performed for isotropic turbulence and the LTCs are calculated using the passive vector method. Both a priori and a posteriori tests are carried out. It is observed that the subgrid-scale contributions to the LTCs cannot be simply ignored and the LES overpredicts the LTCs than the DNS. It is concluded from the straining hypothesis that an accurate prediction of enstrophy spectra is most critical to the prediction of the LTCs.

Modelling the long-term sustainability of indigenous hunting in Manu National Park, Peru: landscape-scale management implications for Amazonia

P> Widespread hunting throughout Amazonia threatens the persistence of large primates and other vertebrates. Most studies have used models of limited validity and restricted spatial and temporal scales to assess the sustainability. We use human-demographi

PHYSICAL MODEL FOR SMALL SCALE SANDWAVES MIGRATION IN NORTH GULF OF SOUTH CHINA SEA

The aim of this paper is to investigate the mechanism of small scale sand-wave migration. 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 wave migration. The calculation results are shown to be consistent with the observed data in the trough of sand ridge. Considering the effect of environmental actions and sand wave features, we develop an effective formula to predict sand-wave migration. It is indicated that the physical models should be used to predict the migration of the small scale sand-wave, which is rarely dominated by wave activity.