暖温带山地森林生态系统水分问题研究

本文的研究是中国科学院院重大项目“暖温带森林生态系统结构、功能及生产力持续发展”的主要内容之一。作者以详实的第一手资料，从森林小气候及环境特征、森林降水的水文学效应及降水化学、森林的热量平衡及蒸发散、树木个体的水分生理生态学几个方面阐述、分析了暖温带山地森林生态系统主要林分的水分及其相关生态学问题。 在森林小气候及环境特征一章，作者从不同季节的日变化和生长季的月际变化两个视角，以落叶阔叶混交林和油松林为研究对象，考察了林冠上和林下四个不同梯度的风速、气温、湿度、地温的时空动态。 在森林降水的水文学效应和降水化学一章，笔者以1993、1994年试验年度的83次降雨观测资料为基础，分析了暖温带落叶阔叶混交林、辽东栎林、油松林、落叶松林、次生灌丛降水总量与各降水分量的关系，建立了单次降雨与各降雨分量的经验模型，并给出了生长季林冠作用层和林地作用层的水量分配的月际动态。在探讨上述水量关系的同时，作者还分析了前四类林分大气降水及各降水分量中N、K、Ca、S、Mg、P、Al七种元素的浓度及含量变化，就不同树种对上述元素的选择性交换作了探讨，比较了不同林分的降水化学效应差异。 在第四章，作者以落叶阔叶混交林和油松林为研究对象，分析了两类林分在94试验年度生长季辐射平衡、显热通量、潜热通量、蒸发散以及土壤热通量的季节变化和日变化特征。 在树木个体的水分生理生态部分，作者应用压力室一容积技术测定了暖温带落叶阔叶林、油松林和次生灌丛10种主要树种的水分生理指标：日最低水势值、最大膨压时的渗透势、膨压为零时的渗透势、初始质壁分离时渗透水的相对含量、初始质壁分离时的相对含水量、质外体水的相对含量、细胞最大弹性模量，并比较了不同树种间上述指标与抗旱性的关系。此外，作者还应用Li-1600稳态气孔计测定了上述林分中主要树种的日均蒸腾强度的季节动态，并比较了上下两面叶片蒸腾特性的差异。最后，作者采用九种水分生理指标对10种主要树种的抗旱性作了主分量分析，给出了综合性抗旱指标。 在第六章，作者应用热脉冲技术系统地研究了暖温带山地森林主要乔木树种的树干液流的时空变化特征，并应用时序分析方法对上述树种的树液流量变化建立了自回归模型，在此基础上提出了生理惯性指标，给予了生理学解释。

浑善达克退化生态系统恢复研究

　　浑善达克沙地是我国四大沙地之一，地处北方干旱半干旱区，为草原区向荒漠区过渡的地带。长期以来，由于人类不合理的生产活动，加上这里脆弱的生态环境，已引起了严重的土地退化问题：流动沙丘面积由1950s年占沙地总面积的2%增加到1990s的近50%。因此浑善达克沙地成为我国研究土地退化、防治沙尘暴的重点地区，本文从自然科学和人文科学相结合的角度出发，对浑善达克沙地草地退化原因、自然恢复潜力、恢复过程、适宜物种选择，以及社区生存、生产等方面进行了综合研究，得出以下主要结论： 　　 1) 浑善达克沙地土壤种子库中含有大量的种子，在退化草地自然恢复中表现出极大的潜力。这些种子在摆脱人和牲畜干扰的前提下，可以萌发、定居并形成植物群落，使退化草地恢复;当地表达到一定程度的植被覆盖，可有效减少沙尘暴的危害。土壤种子库中的植物种类与地上植被有极显著的相似性(P<0.05)，这是地上植被形成稳定群落的基础。种子库中的物种组成影响植被恢复演替的进程;反之，恢复演替也制约着种子库组成和幼苗建立。 　　2)退化沙地草地围封后，对不同恢复阶段草地的群落学调查表明，该地区自然恢复过程大致分为3个阶段：围封2年的恢复早期，流动沙丘向半固定沙丘转变;围封3-5年的恢复中期，半固定沙丘向固定沙丘转变;围封6年后的恢复后期，为固定沙丘稳定发展阶段。根据生活型及植物种类随恢复演替的变化规律，浑善达克沙地植被演替的总体趋势可归纳为：沙米 (Agriophyllum squarrosum)+ 雾冰藜(Bassia dasyphylla)群落→黄柳(Salix gordejevii)+ 冰草(Agropyron cristatum)群落→褐沙蒿(Artemisia intramongolica) + 冰草群落→沙地榆(Ulmus pumila var. sabulosa)疏林+冰草群落。在围封禁牧下，浑善达克退化沙地草地在较短时间内实现自然恢复，因此制约退化草地恢复演替的关键因素主要是人为和牲畜的干扰，只要排除了这种干扰因素，浑善达克大面积的退化沙地草地完全能够借助自然力实现生态恢复。 　　3) 浑善达克沙地3种生境下84种植物叶片渗透势值和含水量，表现出不同功能型上的差异。总体变化趋势为： 深根系 > 浅根系;灌木 > 乔木 > 草本;分布在湿地和丘间低地的植物叶片渗透势和含水量较高，而生长在沙丘上的植物叶片渗透势较低，需要有发达的根系吸收土壤深层的水。不同植物具有独特的水分利用特性，使它们能共存于同一生态系统中。这些不同植物功能型表现出的植物水分生理生态特性，表明浑善达克天然分布的植物群落发育有完善的利用水资源的能力，能够保证在很大降水波动条件下分布有丰富的植物群落和较高的生物生产力，构成该特殊类型生态系统很强的恢复潜力。另外，浑善达克沙地沙丘的存在是该类生态系统恢复弹性较高的另一重要原因。 　　4)本地种与引入种在生理生态上表现出不同的适应能力。在相似的太阳辐射和叶片温度下，引入种旱柳的叶片水势较高，而净光合速率、气孔导度、水分利用效率则较低。这表明它的光合潜能在改变环境中没有正常发挥。同时，引入种较低的最大光化学效率进一步表明它抵抗环境胁迫的能力较低。当土壤水分可利用程度降低而导致水分竞争时，引入种很可能在竞争中被淘汰。因此，在生态恢复中，应尽量避免引入外来种，大量使用本地种。 　　5) 生态恢复不仅是自然科学问题，更重要的是社会经济问题。为了充分认识当地社区的参与对生态恢复的作用，在实验过程中，调查了当地居民自本实验开展以来的思想观念、经济收入和生产效益等变化。在生态恢复中当地社区的积极参与是保证恢复成功与否的关键因素。从发挥“自然力” 和“以人为本”的指导思想出发，在生态恢复中应注重充分利用自然的力量;在管理方面，要以解决社区居民的生产生活实际需求为目标。只有这样，才能保证生态、社会和经济可持续发展。

土壤干湿交替对玉米生长速度及其耗水量的影响

在盆栽条件下研究了干湿交替对玉米生长速率、叶片水势、渗透势、气孔导度、相对生长速率和耗水量的影响。结果表明玉米在 3～ 7叶期经历土壤水分缓慢亏缺 ,再进行复水的干湿交替后玉米叶片渗透调节能力明显增加 ,叶片生长表现出补偿效应 ,每次干湿交替后生长速率迅速下降的叶水势趋于下降 ,气孔导度对土壤水分变化非常敏感 ,并在干旱—复水过程中具有后效作用 ,蒸腾耗水量随干—湿交替而具有下降趋势 ,初步证明可在节水灌溉条件下人为控制不同生育时期的供水时间形成干湿交替 ,促进渗透调节能力增强和补偿生长来实现作物高产、高效、优质的目的

Electro-Osmotic Flow and Mixing in Heterogeneous Microchannels

Analytical and numerical studies of secondary electro-osmotic flow EOF and its mixing in microchannels with heterogeneous zeta potentials are carried out in the present work. The secondary EOFs are analyzed by solving the Stokes equation with heterogeneous slip velocity boundary conditions. The analytical results obtained are compared with the direct numerical simulation of the Navier-Stokes equations. The secondary EOFs could transport scalar in larger areas and increase the scalar gradients, which significantly improve the mixing rate of scalars. It is shown that the heterogeneous zeta potentials could generate complex flow patterns and be used to enhance scalar mixing.

Comparative studies of AOT and NaDEHP by fluorescence technique and z-potential (ζ) measurements

The aggregation behaviors of two surfactants with the same hydrophobic tail, sodium bis(2-ethylhexyl)sulfosuccinate (AOT) and sodium bis(2-ethylhexyl)phosphate (NaDEHP), have been investigated by the fluorescence technique and z-potential (ζ) measurements. Five fine peaks of the pyrene molecule fluorescence spectroscopy appear in the surfactant solution, and the micropolarity at which pyrene locates is monitored from the intensity ratio of the first (I1) and the third peak (I3). A wide peak around 475 nm, the emission spectra of the excimer of pyrene molecules, is observed in the NaDEHP solution, while this is not found for the AOT system. The value of I1/I3 decreases in a more limited concentration range for the AOT system than for NaDEHP, indicating that small aggregates can be more easily formed by NaDEHP molecules. The z-potential results for the aggregates formed by the two surfactants show that the interaction between AOT and PVP is stronger than that between NaDEHP and PVP.

A Lennard-Jones Embedded-Atom Potential and Its Application to the Study of Melting

The nearest-neighbour Lennard-Jones potential from the embedded-atom method is extended to a form that includes more than nearest neighbours. The model has been applied to study melting with molecular dynamics. The calculated melting point, fractional volume change on melting, heat of fusion and linear coefficients of thermal expansion are in good agreement with experimental data. We have found that the second and third neighbours influence the melting point distinctly.

Dynamical symmetry of screened Coulomb potential and isotropic harmonic oscillator

It is shown that for the screened Coulomb potential and isotropic harmonic oscillator, there exists an infinite number of closed orbits for suitable angular momentum values. At the aphelion (perihelion) points of classical orbits, an extended Runge-Lenz vector for the screened Coulomb potential and an extended quadrupole tensor for the screened isotropic harmonic oscillator are still conserved. For the screened two-dimensional (2D) Coulomb potential and isotropic harmonic oscillator, the dynamical symmetries SO3 and SU(2) are still preserved at the aphelion (perihelion) points of classical orbits, respectively. For the screened 3D Coulomb potential, the dynamical symmetry SO4 is also preserved at the aphelion (perihelion) points of classical orbits. But for the screened 3D isotropic harmonic oscillator, the dynamical symmetry SU(2) is only preserved at the aphelion (perihelion) points of classical orbits in the eigencoordinate system. For the screened Coulomb potential and isotropic harmonic oscillator, only the energy (but not angular momentum) raising and lowering operators can be constructed from a factorization of the radial Schrodinger equation.

A new N-body potential and its application

Based on the embedded atom method (EAM) proposed by Daw and Baskes and Johnson's model, this paper constructs a new N-body potential for bcc crystal Mo. The procedure of constructing the new N-body potential can be applied to other metals. The dislocation emission from a crack tip has been simulated successfully using molecular dynamics method, the result is in good agreement with the elastic solution.

The positive floating potential of high-temperature particles in plasma

It has been predicted that the floating potential of particles in plasma may become positive when the particle surface temperature is high enough, but, to our knowledge, no positive floating potential has been obtained yet. In the present paper the floating potential theory of high-temperature particles in plasma is developed to cover the positive potential range for the first time, and a general approximate analytical formula for the positive floating potential with a thin plasma sheath and subsonic plasma flow is derived from the new model recently proposed by the authors. The results show that when the floating potential is positive, the net flux of charge incident on the particle approaches a constant similar to the 'electron saturation' phenomena in the case of the electric probes.

A new model for the floating potential of fine particles in plasma

In the plasma processing of ultrafine particles of material, the heat transfer and force are considerably affected by particle charging. In this communication a new model, including thermal electron emission and incorporating the effect of electric field near the particle surface, is developed for metallic spherical particles under the condition of a thin plasma sheath. Based on this model, the particle floating potential, and thus the heat transfer and force, can be detemined more accurately and more realistically than previously.

Potential energy model for turbulent entrainment across density interface in a linearly stratified fluid

A potential energy model is developed for turbulent entrainment in the absence of mean shear in a linearly stratified fluid. The relation between the entrainment distance D and the time t and the relation between dimensionless entrainment rate E and the local Richardson number are obtained. An experiment is made for examination. The experimental results are in good agreement with the model, in which the dimensionless entrainment distance D is given by DBAR = A(i)(SBAR)-1/4(fBAR)1/2(tBAR)1/8, where A(i) is the proportional coefficient, S is the dimensionless stroke, fBAR is the dimensionless frequency of the grid oscillation, tBAR the dimensionless time.

Viscosity-induced mode splitting and potential energy criterion for mode stability

This paper points out that viscosity can induce mode splitting in a uniform infinite cylinder of an incompressible fluid with self-gravitation, and that the potential energy criterion cannot be appropriate to all normal modes obtained, i.e., there will be stable modes with negative potential energy (<0). Therefore the condition >0 is not necessary, although sufficient, for the stability of a mode in an incompressible static fluid or magnetohydrodynamics (MHD) system, which is a correction of both Hare's [Philos. Mag. 8, 1305 (1959)] and Chandrasekhar's [Hydrodynamic and Hydromagnetic Stability (Oxford U.P., Oxford, 1961), p. 604] stability criterion for a mode. These results can also be extended to compressible systems with a polytropic exponent.

High-Strength Composite Fibers: Realizing True Potential of Carbon Nanotubes in Polymer Matrix through Continuous Reticulate Architecture and Molecular Level Couplings

Carbon nanotubes have unprecedented mechanical properties as defect-free nanoscale building blocks, but their potential has not been fully realized in composite materials due to weakness at the interfaces. Here we demonstrate that through load-transfer-favored three-dimensional architecture and molecular level couplings with polymer chains, true potential of CNTs can be realized in composites as Initially envisioned. Composite fibers with reticulate nanotube architectures show order of magnitude improvement in strength compared to randomly dispersed short CNT reinforced composites reported before. The molecular level couplings between nanotubes and polymer chains results in drastic differences in the properties of thermoset and thermoplastic composite fibers, which indicate that conventional macroscopic composite theory falls to explain the overall hybrid behavior at nanoscale.

Splitting the fast and slow motions in molecular dynamics simulations based on the change of cold potential well bottom

Abstract. The atomic motion is coupled by the fast and slow components due to the high frequency vibration of atoms and the low frequency deformation of atomic lattice, respectively. A two-step approximate method was presented to determine the atomic slow motion. The first step is based on the change of the location of the cold potential well bottom and the second step is based on the average of the appropriate slow velocities of the surrounding atoms. The simple tensions of one-dimensional atoms and two-dimensional atoms were performed with the full molecular dynamics simulations. The conjugate gradient method was employed to determine the corresponding location of cold potential well bottom. Results show that our two-step approximate method is appropriate to determine the atomic slow motion under the low strain rate loading. This splitting method may be helpful to develop more efficient molecular modeling methods and simulations pertinent to realistic loading conditions of materials.