岷江干旱河谷三种蔷薇植物繁殖更新研究

繁殖更新是植物生活史的重要阶段，在退化生态系统中，植物繁殖更新能力往往较差，是植被恢复的限制环节，因而也成为恢复研究重点和核心。本研究选择岷江干旱河谷广泛分布的三种蔷薇：多苞蔷薇（R. multibracteata）、黄蔷薇（R. hugonis）和川滇蔷薇（R. soulieana）为研究对象，通过野外调查，在查明其生长、繁殖更新状况的基础上，采用控制和模拟实验，对种子和幼苗阶段进行了深入研究，综合分析更新潜力，并提出相对应的促进更新和植被恢复措施。主要结论如下： 1）三种蔷薇在岷江干旱河谷广泛分布，生长和繁殖状况良好，结实量大。各生长指标：株高、基径和冠幅，繁殖指标：结实数量、重量和单果重量都具有显著的空间差异性。基径对多苞蔷薇结实量影响最大；而冠幅对黄蔷薇结实量影响最大。海拔和纬度是对蔷薇生长和繁殖影响最大的环境因素，随着海拔和纬度的升高，植株生长更高大，结实量增加；坡度和坡向对其生长和繁殖也有一定影响，随着坡度 和坡向增加，蔷薇生长和结实受到抑制。 2）三种蔷薇在岷江干旱河谷更新现状不佳, 但更新潜力大。活力种子比率低，动物取食以及两年生幼苗的大量死亡是蔷薇更新的主要限制因素。多苞蔷薇和黄蔷薇的结实率低，川滇蔷薇较高。三种蔷薇种子产量大，但种子质量较差，更新具有充足的种源。三种蔷薇都能形成持久种子库，种子库中种子总量大，但有效种子少，黄蔷薇被动物啃食的比例很高，多苞蔷薇和川滇蔷薇也有一部分种子受到动物破坏。三种蔷薇幼苗库组成特征表现为，当年生幼苗所占比例很高，年龄较大幼苗所占比例小。 3）三种蔷薇都具有不同程度休眠，未经处理种子的发芽率极低。黄蔷薇休眠程度最深，为深度生理休眠；多苞蔷薇为中度生理休眠；川滇蔷薇为非深度生理休眠。三种蔷薇种子在形态上发育成熟，种皮具有透水性。蔷薇果果肉和瘦果中含有抑制物质，其浸泡液抑制了油菜种子萌发，果肉抑制作用更强，果肉和瘦果浸泡液的抑制程度分别为：川滇蔷薇>黄蔷薇>多苞蔷薇。切割和硫酸腐蚀提高了川滇蔷薇种子的发芽能力，而对多苞蔷薇和黄蔷薇没有影响。完全去除瘦果果皮和种皮提高了多苞蔷薇种子发芽率，但对黄蔷薇没有影响。赤霉素和烟水对蔷薇种子萌发没有促进作用。三种蔷薇打破休眠所需低温层积时间分别为：黄蔷薇>多苞蔷薇>川滇蔷薇。对于多苞蔷薇和川滇蔷薇，层积前对种子进行硫酸腐蚀或暖温层积能缩短低温层积时间，提高发芽率。对于多苞蔷薇，变温层积中暖温层积和低温层积具有一定的负补性，即延长暖温层积可以缩短种子萌发对低温层积的需要。 4）多苞蔷薇种子形态特征和种子休眠与萌发在不同海拔梯度间存在较大差异。种子采集时间、采集季节和干藏影响多苞蔷薇和川滇蔷薇的种子休眠。多苞蔷薇果实大小、种子大小和千粒重、种皮厚度随海拔升高而增加，而种子饱满率和活力随海拔升高而降低，种子休眠程度也随海拔升高而增加。种皮厚度与种子大小、千粒重成正相关关系，硫酸腐蚀后的种子经过不同时间的低温层积后，种子发芽率与种皮厚度、种子大小、千粒重、海拔成正相关关系。2006 年采集川滇蔷薇和多苞蔷薇种子休眠程度较2005 年低。种子休眠随种子年龄增加而减弱。高温和干旱能减轻多苞蔷薇和川滇蔷薇种子休眠。 5）三种蔷薇的生长和生物量积累在干旱胁迫条件下受到抑制，而生物量分配、叶片形态特征和水分利用特征等都发生了变化。三种蔷薇的根、茎、叶各器官生物量以及总生物量等在干旱胁迫下明显减小，叶片脱落数量增加。在干旱胁迫条件下，较多的生物量分配到地下部分，从而这使R/S 明显增加。比叶面积（SLA）和冠层面积比（LAR）对干旱胁迫的反应不敏感，仅有部分物种在干旱胁迫条件下发生了变化，并且其变化特点在不同年龄幼苗之间有一定差异。干旱胁迫对WUE 的影响在不同物种间存在差异。多苞蔷薇和黄蔷薇的WUE 随着干旱胁迫的增加而增大， 而川滇蔷薇的WUE 则随干旱胁迫增加而减小。在干旱胁迫条件下，多苞蔷薇和黄蔷薇叶片脱落量和生物量减小幅度较川滇蔷薇大，表明其抗旱能力较强。在干旱胁迫条件，三种蔷薇两年生幼苗的生物量减小幅度较当年生幼苗小，表明两年生幼苗的抗旱能力更强。 6）两种植被恢复措施中，幼苗移栽比播种具有更好的植被恢复效果。播种后，蔷薇种子的发芽率较高，但出苗率都很低，即使出苗，幼苗也几乎在一月内全部死亡。 三种微生境条件下（灌木、半灌木和裸地），种子出苗和幼苗成活没有差异。移栽幼苗总体死亡率都比较低，小于20％。特别是两年生幼苗死亡率更低，小于2％。移栽后的幼苗生长状况良好，在整个生长季中，各生长指标不断增加。生境对幼苗的存活率没有显著影响，但对于幼苗的生长和生物量积累有一定影响，裸地更有利于幼苗生长和生物量积累。与当年生幼苗相比，两年生幼苗具有更高的成活率。总之，三种蔷薇在干旱河谷分布广泛、生长繁殖状况良好，结实量大，具有丰富种源，繁殖更新潜力大，但繁殖更新状况不佳；种子散布后动物对种子的取食、种子的深度休眠过程、种子出苗以及当年生幼苗的存活和定居是更新的主要限制环节。水分是影响结实、种子休眠解除和萌发，幼苗存活和定居的最主要的限制因素。在植被恢复中，应在种子成熟季节大量采集种子，在室内打破休眠后进行人工播种，培育两年生幼苗，通过幼苗移栽方式进行植被恢复。川滇蔷薇应栽种在相对湿润的过渡区，而多苞蔷薇和黄蔷薇可以应用于核心区植被恢复。 Regeneration is an important phase in plant life cycle. It has been a key component of ecological restoration in degradation ecosystem in which plants commonly has poor regeneration. In this paper, we investigated the natural growth, propagation and regeneration status of native three rose species, Rosa multibracteata, R. hugonis and R. soulieana, and analysis the limitation in seed germination and seedling establishment stages. Advice on facilitating the use of these plants in restoration based on the results has been proposed. The results were as follows: 1) Three rose plants widely distributed in the dry valley of the Minjiang River, and made a good performance in growth and propagation. There were significant spatial differences in each growth parameter, such as ramet height, basal diameter, crown diameter and propagation parameters including hip number of a clump, hip mass of a clump and a hip mass. Basa diameter was the most important growth parameter influencing fruit number for R. multibracteata and crown diameter was for R. hugoni. Altitude and latitude had the greatest effect on the growth and propagation of rose plants among environmental conditions. Each parameter of growth and propagation increased with the increase of altitude and latitude. In addition, the increase of slope and aspect limited the growth and propagation. 2) Three rose plants had poor natural regeneration, but great regeneration potential. Low seed viability, predation and higher mortality of current year old seedlings were the limitation in regeneration. R. multibracteata and R. hugonis had higher fruiting rates than R. souliean. All three plants produced a great number of seeds, while their viability was poor. Three rose plants had persistent seed banks, with high total seed number but very low viable seed density. Predation was most severe in R. hugonis, and it also existed to some degree in R. multibracteata and R. soulieana. The seedling age-structure was characteristic of current-year seedlings predominating and few older seedlings were observed. 3) Three rose seeds were dormant and untreated seeds germinated with very low germination percentages. The rose seeds had morphological mature embryos, and achenes were permeable. Some inhabit substances existed in hips and achenes for the extracts of hips and achenes inhibited germination of Brassica campestris. The inhibition effect of the extracts of three rose hip and achenes was R. soulieana>R. hugonis>R. multibracteata. Mechanical and H2SO4 scarification increased R. soulieana germination but had no effect on germination of R. hugonis and R. multibracteata seeds. Full removal of pericarp and testa improved the germination of R. multibracteata but did not affect R. hugonis germination. GA3 and smoke water had no positive effect on rose seed germination. The periods of cold stratification required to released seed dormancy was R. hugonis > R. soulieana >R. multibracteata. H2SO4 scarification and warm stratification shortened cold stratification to release dormancy for R. soulieana and R. multibracteata. Warm stratification had complementary effect for cold stratification, i.e. the longer warm stratification seeds received, the shorter cold stratification were required to obtain the same germination percentage. Three rose seeds had different kinds of dormancy; R. hugonis has deep physiological dormancy, R. multibracteata with intermediate physiological dormancy and R. souliean non-deep physiological dormancy. 4）The seeds traits and dormancy of R. multibracteata showed significant difference across altitudes. Year and season of seed collection had significant effect on seed dormancy for both R. souliean and R. multibracteata. Hip size, seed size, seed weight, seed coat thickness and seed dormancy level increased with the increase of the altitude. There were positive relations between seed coat thickness with seed size and seed weight. Germination percentage of seeds treated with H2SO4 scarification following different periods of cold stratification showed positive relation with seed coat thickness, seed size, seed weight and altitude. Seeds of R. souliean and R. multibracteata collected in 2006 had low dormancy level than those collected in 2005. Seed dormancy decreased with increasing seeds age. High temperature and drought were associated with low dormancy level. 5) Seedling growth, the total dry mass and their components of seedlings were reduced, while leaf senescence accelerated under drought stress. More biomass allocation to root system resulted in higher R/S ratio under drought. Water-use efficiency (WUE) of R. multibracteata and R. hugonis increased, while it declined for R. soulieana under drought stress. R. soulieana seedlings had poor drought-resistance capacity it had more senescent leaves, and its reduction of biomass was stronger than two other rose plants under drought. The reduction degree of one year old seedlings under drought stress was slighter than that of current year seedlings. Therefore, one year old seedling was more drought-resistent compared to current year seedlings. 6）Planting seedlings may have better effect in comparison with direct seeding. Most seeds germinated after seeding, but seedling emergence was very low. More than 80 % seedlings from direct seeding died within a months after emergence. Seedling emergence and survival rate did not show difference among microhabitats. Mortality rates of seedlings artificially planted in microhabitats were general lower than 20 %, and the mortality rate of one year old seedlings was lower than 2 %. Each grow parameter including plant height, leaf number and branch number continually increased after planting. Microhabitat type had effect on the growth parameter and biomass production, but it did not influence the seedling survival. Bare land tended to facilitate seedling growth. One year old seedlings had higher survival rate than current year seedlings. In conculsion, the three rose had wide distribution in the dry valley of the Minjiang River. They produced many seeds and had tolerance to drought stress to some degree. But they had poor regeneration in habitats may be caused by predation, seed dormancy，and high mortality in current year seedlings. We recommend that rose plants should be utilized in restoration by planting two-year old seedlings in spring. A large quantity of seeds should be collected artificially in autumn, release seed dormancy in room, and then cultivate two-year old seedlings by seeding in particular container. R. soulieana seedling probably be planted in transition area, and R. multibracteata and R. hugonis can be used in core area of the dry valley of the Minjiang River.

Direct Numerical Simulation of a Spatially Evolving Supersonic Turbulent Boundary Layer at Ma=6

Direct numerical simulation is carried out for a spatially evolving supersonic turbulent boundary layer at free-stream Mach number 6. To overcome numerical instability, the seventh-order WENO scheme is used for the convection terms of Navier-Stokes equations, and fine mesh is adopted to minimize numerical dissipation. Compressibilty effects on the near-wall turbulent kinetic energy budget are studied. The cross-stream extended self-similarity and scaling exponents including the near-wall region are studied. In high Mach number flows, the coherence vortex structures are arranged to be smoother and streamwised, and the hair-pin vortices are less likely to occur.

Liquid phase epitaxy growth of AlxGayIn1-x-yPzAs1-zGaAs with direct band gap up to 2.0 eV

III-V pentenary semiconductor AlGaInPAs with a direct band gap of up to 2.0 eV has been grown successfully on GaAs substrates by liquid phase epitaxy;(LPE). With the introduction of the energy bowing parameters of quaternaries, the theoretical calculations agree well with the measured PL peak energy data from pentenary samples.

Direct numerical test of the B-G-K model equation by the DSMC method

In the present paper the rarefied gas how caused by the sudden change of the wall temperature and the Rayleigh problem are simulated by the DSMC method which has been validated by experiments both in global flour field and velocity distribution function level. The comparison of the simulated results with the accurate numerical solutions of the B-G-K model equation shows that near equilibrium the BG-K equation with corrected collision frequency can give accurate result but as farther away from equilibrium the B-G-K equation is not accurate. This is for the first time that the error caused by the B-G-K model equation has been revealed.

Direct numerical simulation of transition and turbulence in compressible mixing layer

The three-dimensional compressible Navier-Stokes equations are approximated by a fifth order upwind compact and a sixth order symmetrical compact difference relations combined with three-stage Ronge-Kutta method. The computed results are presented for convective Mach number Mc = 0.8 and Re = 200 with initial data which have equal and opposite oblique waves. From the computed results we can see the variation of coherent structures with time integration and full process of instability, formation of Lambda-vortices, double horseshoe vortices and mushroom structures. The large structures break into small and smaller vortex structures. Finally, the movement of small structure becomes dominant, and flow field turns into turbulence. It is noted that production of small vortex structures is combined with turning of symmetrical structures to unsymmetrical ones. It is shown in the present computation that the flow field turns into turbulence directly from initial instability and there is not vortex pairing in process of transition. It means that for large convective Mach number the transition mechanism for compressible mixing layer differs from that in incompressible mixing layer.

Direct Numerical Simulation of Passive Scalar in Decaying Compressible Turbulence

The passive scalars in the decaying compressible turbulence with the initial Reynolds number (defined by Taylor scale and RMS velocity) Re=72, the initial turbulent Mach numbers (defined by RMS velocity and mean sound speed) Mt=0.2-0.9, and the Schmidt numbers of passive scalar Sc=2-10 are numerically simulated by using a 7th order upwind difference scheme and 8th order group velocity control scheme. The computed results are validated with different numerical methods and different mesh sizes. The Batchelor scaling with k(-1) range is found in scalar spectra. The passive scalar spectra decay faster with the increasing turbulent Mach number. The extended self-similarity (ESS) is found in the passive scalar of compressible turbulence.

An analysis of subgrid-resolved scale interactions with use of results from direct numerical simulations

Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics have been investigated. It is observed that subgrid-subgrid interaction dominates the turbulent dynamics when cut-off wave number locates in the energy-containing range while resolved-subgrid interaction dominates if it is in the dissipation range; By decomposing the subgrid energy transfer and nonlinear interaction into 'forward' and 'backward' groups according to the sign of triadic interaction, we find that individually each group has very large contribution, but the net of them is much smaller, implying that tremendous cancellation happens between these two groups.

Direct Numerical Simulation of Supersonic Turbulent Boundary Layer Flow

Direct numerical simulations of a spatially evolving supersonic flat-plate turbulent boundary layer flow with free Mach number M = 2.25 and Reynolds number Re = 365000/in are performed. The transition process from laminar to turbulent flow is obtained by solving the three-dimensional compressible Navier-Stokes, equations, using high-order accurate difference schemes. The obtained statistical results agree well with the experimental and theoretical data. From the numerical results it can be seen that the transition process under the considered conditions is the process which skips the Tolimien-Schlichting instability and the second instability through the instability of high gradient shear layer and becomes of laminar flow breakdown. This means that the transition process is a bypass-type transition process. The spanwise asymmetry of the disturbance locally upstream imposed is important to induce the bypass-type transition. Furthermore, with increasing the time disturbance frequency the transition will delay. When the time disturbance frequency is large enough, the transition will disappear.

可压剪切层线性稳定性特征直接数值模拟=Direct Simulation of the Linear Instabilities in Compressible Shear Layers

基于伪随机数生成技术促生白噪声扰动,以高精度迎风/对称紧致混合差分算法求解二维/三维非定常可压Navier-Stokes方程,揭示了可压自由剪切层初始剪切过程中扰动的线性演化特征,以及该过程对扰动波数和方向的内在选择性.验证了所用算法的有效性,表明线性理论同数值模拟相结合是可压剪切层研究的合理途径之一.

Arc Root Motions In An Argon-Hydrogen Direct-Current Plasma Torch At Reduced Pressure

Arc root motions in generating dc argon-hydrogen plasma at reduced pressure are optically observed using a high-speed video camera. The time resolved angular position of the arc root attachment point is measured and analysed. The arc root movement is characterized as a chaotic and jumping motion along the circular direction on the anode surface.

Direct Numerical Simulation Of Hypersonic Boundary-Layer Transition Over A Blunt Cone

Direct numerical simulation of transition How over a blunt cone with a freestream Mach number of 6, Reynolds number of 10,000 based on the nose radius, and a 1-deg angle of attack is performed by using a seventh-order weighted essentially nonoscillatory scheme for the convection terms of the Navier-Stokes equations, together with an eighth-order central finite difference scheme for the viscous terms. The wall blow-and-suction perturbations, including random perturbation and multifrequency perturbation, are used to trigger the transition. The maximum amplitude of the wall-normal velocity disturbance is set to 1% of the freestream velocity. The obtained transition locations on the cone surface agree well with each other far both cases. Transition onset is located at about 500 times the nose radius in the leeward section and 750 times the nose radius in the windward section. The frequency spectrum of velocity and pressure fluctuations at different streamwise locations are analyzed and compared with the linear stability theory. The second-mode disturbance wave is deemed to be the dominating disturbance because the growth rate of the second mode is much higher than the first mode. The reason why transition in the leeward section occurs earlier than that in the windward section is analyzed. It is not because of higher local growth rate of disturbance waves in the leeward section, but because the growth start location of the dominating second-mode wave in the leeward section is much earlier than that in the windward section.

Direct Numerical Simulation Of Three-Dimensional Richtmyer-Meshkov Instability

Direct numerical simulation (DNS) is used to study flow characteristics after interaction of a planar shock with a spherical media interface in each side of which the density is different. This interfacial instability is known as the Richtmyer-Meshkov (R-M) instability. The compressible Navier-Stoke equations are discretized with group velocity control (GVC) modified fourth order accurate compact difference scheme. Three-dimensional numerical simulations are performed for R-M instability installed passing a shock through a spherical interface. Based on numerical results the characteristics of 3D R-M instability are analysed. The evaluation for distortion of the interface, the deformation of the incident shock wave and effects of refraction, reflection and diffraction are presented. The effects of the interfacial instability on produced vorticity and mixing is discussed.

Direct dynamical test for deterministic chaos and optimal embedding of a chaotic time-series

We propose here a local exponential divergence plot which is capable of providing an alternative means of characterizing a complex time series. The suggested plot defines a time-dependent exponent and a ''plus'' exponent. Based on their changes with the embedding dimension and delay time, a criterion for estimating simultaneously the minimal acceptable embedding dimension, the proper delay time, and the largest Lyapunov exponent has been obtained. When redefining the time-dependent exponent LAMBDA(k) curves on a series of shells, we have found that whether a linear envelope to the LAMBDA(k) curves exists can serve as a direct dynamical method of distinguishing chaos from noise.

Direct dynamical test for deterministic chaos

We present a direct and dynamical method to distinguish low-dimensional deterministic chaos from noise. We define a series of time-dependent curves which are closely related to the largest Lyapunov exponent. For a chaotic time series, there exists an envelope to the time-dependent curves, while for a white noise or a noise with the same power spectrum as that of a chaotic time series, the envelope cannot be defined. When a noise is added to a chaotic time series, the envelope is eventually destroyed with the increasing of the amplitude of the noise.

Upwind compact method and direct numerical-simulation of driven flow in a square cavity

A high-order accurate finite-difference scheme, the upwind compact method, is proposed. The 2-D unsteady incompressible Navier-Stokes equations are solved in primitive variables. The nonlinear convection terms in the governing equations are approximated by using upwind biased compact difference, and other spatial derivative terms are discretized by using the fourth-order compact difference. The upwind compact method is used to solve the driven flow in a square cavity. Solutions are obtained for Reynolds numbers as high as 10000. When Re less than or equal to 5000, the results agree well with those in literature. When Re = 7500 and Re = 10000, there is no convergence to a steady laminar solution, and the flow becomes unsteady and periodic.