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.

Discussion on the geometric structure of strange attractor

We discuss the transversal heteroclinic cycle formed by hyperbolic periodic pointes of diffeomorphism on the differential manifold. We point out that every possible kind of transversal heteroclinic cycle has the Smalehorse property and the unstable manifolds of hyperbolic periodic points have the closure relation mutually. Therefore the strange attractor may be the closure of unstable manifolds of a countable number of hyperbolic periodic points. The Henon mapping is used as an example to show that the conclusion is reasonable.

The strange attractor of a kind of 2-dimensional map and dynamical properties on it

On the basis of previous works, the strange attractor in real physical systems is discussed. Louwerier attractor is used as an example to illustrate the geometric structure and dynamical properties of strange attractor. Then the strange attractor of a kind of two-dimensional map is analysed. Based on some conditions, it is proved that the closure of the unstable manifolds of hyberbolic fixed point of map is a strange attractor in real physical systems.

The geometric structure and dynamic properties of lauwerier attractor

 Introduction The strange chaotic attractor (ACS) is an important subject in the nonlinear field. On the basis of the theory of transversal heteroclinic cycles, it is suggested that the strange attractor is the closure of the unstable manifolds of countable infinite hyperbolic periodic points. From this point of view some nonlinear phenomena are explained reasonably. 更多还原

Fatigue crack-growth rate in ferrite martensite dual-phase steel

An empirical study is made on the fatigue crack growth rate in ferrite-martensite dual-phase (FMDP) steel. Particular attention is given to the effect of ferrite content in the range of 24.2% to 41.5% where good fatigue resistance was found at 33.8%. Variations in ferrite content did not affect the crack growth rate when plotted against the effective stress intensity factor range  which was assumed to follow a linear relation with the crack tip stress intensity factor range ΔK. A high  corresponds to uniformly distributed small size ferrite and martensite. No other appreciable correlation could be ralated to the microstructure morphology of the FMDP steel. The closure stress intensity factor , however, is affected by the ferrite content with  reaching a maximum value of 0.7. In general, crack growth followed the interphase between the martensite and ferrite.

Dividing the fatigue crack growth process into Stage I and II where the former would be highly sensitive to changes in ΔK and the latter would increase with ΔK depending on the  ratio. The same data when correlated with the strain energy density factor range ΔS showed negligible dependence on mean stress or R ratio for Stage I crack growth. A parameter α involving the ratio of ultimate stress to yield stress, percent reduction of area and R is introduced for Stage II crack growth so that the  data for different R would collapse onto a single curve with a narrow scatter band when plotted against αΔS.

A passive scalar field convected by turbulence

Classical statistical mechanics is applied to the study of a passive scalar field convected by isotropic turbulence. A complete set of independent real parameters and dynamic equations are worked out to describe the dynamic state of the passive scalar field. The corresponding Liouville equation is solved by a perturbation method based upon a Langevin–Fokker–Planck model. The closure problem is treated by a variational approach reported in earlier papers. Two integral equations are obtained for two unknown functions: the scalar variance spectrum F(k) and the effective damping coefficient (k). The appearance of the energy spectrum of the velocity field in the two integral equations represents the coupling of the scalar field with the velocity field. As an application of the theory, the two integral equations are solved to derive the inertial-convective-range spectrum, obtaining F(k)=0.61 −1/3 k−5/3. Here is the dissipation rate of the scalar variance and is the dissipation rate of the energy of the velocity field. This theoretical value of the scalar Kolmogorov constant, 0.61, is in good agreement with experiments.

Nonequilibrium statistical mechanics of one-dimensional turbulence

The method of statistical mechanics is applied to the study of the one-dimensional model of turbulence proposed in an earlier paper. The closure problem is solved by the variational approach which has been developed for the three-dimensional case, yielding two integral equations for two unknown functions. By solving the two integral equations, the Kolmogorov k−5/3 law is derived and the (one-dimensional) Kolmogorov constant Ko is evaluated, obtaining Ko=0.55, which is in good agreement with the result of numerical experiments on one-dimensional turbulence.

Universal equilibrium range of turbulence

An attempt is made to determine the form of F(x), the dimensionless function of universal nature which occurs in the energy spectrum for the universal equilibrium range of fully developed turbulence, by the method of statistical mechanics without introducing any parameter of semiempirical nature. Then, the validity of the variational approach to the closure problem of turbulence theory is tested by applying it to the study of the universal equilbrium range of turbulence.

化学反应流中分子扩散的二点封闭MCA模型

在化学反应流的概率密度函数(PDF)方法中，对流项和化学反应项都是封闭的，但分子扩散项必须模拟．现有的分子扩散模型都是唯象的，需要引入外加参数，并难以通过一些基本物理过程的检验．本文发展了随机映射逼近(mapping closure approximation，MCA)方法，解析地从控制方程导出一个封闭的分子扩散模型．该方法考虑两点联合概率密度函数方程，引入空间特征尺度，因此解决了以往映射封闭方法中分子扩散速率无法确定的问题．数值模拟表明该方法能用于预测标量扩散的速度，以及概率密度函数和条件平均扩散等统计量．

Subgrid scale fluid velocity timescales seen by inertial particles in large-eddy simulation of particle-laden turbulence

Large-eddy simulation (LES) has emerged as a promising tool for simulating turbulent flows in general and, in recent years,has also been applied to the particle-laden turbulence with some success (Kassinos et al., 2007). The motion of inertial particles is much more complicated than fluid elements, and therefore, LES of turbulent flow laden with inertial particles encounters new challenges. In the conventional LES, only large-scale eddies are explicitly resolved and the effects of unresolved, small or subgrid scale (SGS) eddies on the large-scale eddies are modeled. The SGS turbulent flow field is not available. The effects of SGS turbulent velocity field on particle motion have been studied by Wang and Squires (1996), Armenio et al. (1999), Yamamoto et al. (2001), Shotorban and Mashayek (2006a,b), Fede and Simonin (2006), Berrouk et al. (2007), Bini and Jones (2008), and Pozorski and Apte (2009), amongst others. One contemporary method to include the effects of SGS eddies on inertial particle motions is to introduce a stochastic differential equation (SDE), that is, a Langevin stochastic equation to model the SGS fluid velocity seen by inertial particles (Fede et al., 2006; Shotorban and Mashayek, 2006a; Shotorban and Mashayek, 2006b; Berrouk et al., 2007; Bini and Jones, 2008; Pozorski and Apte, 2009).However, the accuracy of such a Langevin equation model depends primarily on the prescription of the SGS fluid velocity autocorrelation time seen by an inertial particle or the inertial particle–SGS eddy interaction timescale (denoted by $\delt T_{Lp}$ and a second model constant in the diffusion term which controls the intensity of the random force received by an inertial particle (denoted by C_0, see Eq. (7)). From the theoretical point of view, dTLp differs significantly from the Lagrangian fluid velocity correlation time (Reeks, 1977; Wang and Stock, 1993), and this carries the essential nonlinearity in the statistical modeling of particle motion. dTLp and C0 may depend on the filter width and particle Stokes number even for a given turbulent flow. In previous studies, dTLp is modeled either by the fluid SGS Lagrangian timescale (Fede et al., 2006; Shotorban and Mashayek, 2006b; Pozorski and Apte, 2009; Bini and Jones, 2008) or by a simple extension of the timescale obtained from the full flow field (Berrouk et al., 2007). In this work, we shall study the subtle and on-monotonic dependence of $\delt T_{Lp}$ on the filter width and particle Stokes number using a flow field obtained from Direct Numerical Simulation (DNS). We then propose an empirical closure model for $\delta T_{Lp}$. Finally, the model is validated against LES of particle-laden turbulence in predicting single-particle statistics such as particle kinetic energy. As a first step, we consider the particle motion under the one-way coupling assumption in isotropic turbulent flow and neglect the gravitational settling effect. The one-way coupling assumption is only valid for low particle mass loading.

Budgets and dynamics of nitrogen and phosphorus in a shallow, hypereutrophic lake in China

Budgets and dynamics of nitrogen and phosphorus in Lake Donghu were investigated from Oct. 1997 to Sept. 1999. The water residence time was estimated to be 89 days in 1997-1998 and 124 days in 1998-1999. The total external loadings were 53 g N m(-2) yr(-1) and 3.2 g P m(-2) yr(-1) in 1997-1998, and 42 g N m(-2) yr(-1) and 3.1 g P m(-2) yr(-1) in 1998-1999. On average, about 80% of nitrogen and phosphorus input was from sewage outlets, while the rest was from land runoff and precipitation. Ammonium ion was the most abundant form of inorganic nitrogen in the sewage. The nutrient output was mainly through water outflow and fish catch. The percentages of nutrients in fish were estimated to be 7.8%-11.2% for nitrogen and 47.6%- 49.6% for phosphorus. Lake Donghu has a very high nutrient retention (63% for nitrogen and 79% for phosphorus) mainly due to its closure and long water residence time. Sedimentation is an important nutrient retention mechanism in this lake. Using mass balance method, we estimated that denitrification of Lake Donghu involves about 50% of the retained nitrogen. Lake Donghu is rich in inorganic nitrogen and phosphorus and showed great seasonal variation.

Domains in ferroelectric nanodots

Almost free-standing single crystal mesoscale and nanoscale dots of ferroelectric BaTiO(3) have been made by direct focused ion beam patterning of bulk single crystal material. The domain structures which appear in these single crystal dots, after cooling through the Curie temperature, were observed to form into quadrants, with each quadrant consisting of fine 90 degrees stripe domains. The reason that these rather complex domain configurations form is uncertain, but we consider and discuss three possibilities for their genesis: first, that the quadrant features initially form to facilitate field-closure, but then develop 90 degrees shape compensating stripe domains in order to accommodate disclination stresses; second, that they are the result of the impingement of domain packets which nucleate at the sidewalls of the dots forming "Forsbergh" patterns (essentially the result of phase transition kinetics); and third, that 90 degrees domains form to conserve the shape of the nanodot as it is cooled through the Curie temperature but arrange into quadrant packets in order to minimize the energy associated with uncompensated surface charges (thus representing an equilibrium state). While the third model is the preferred one, we note that the second and third models are not mutually exclusive.

黄土丘陵区封禁对侵蚀土壤微生物生物量的影响

采用时空互代法,以典型侵蚀环境纸坊沟流域不同封禁年限的狼牙刺群落和杂灌群落为研究对象,选取放牧地和天然次生林为参照,分析了生态恢复过程中土壤微生物生物量、呼吸强度、代谢商及理化性质的演变特征。结果表明,封禁后土壤理化性质明显改善;微生物生物量随封禁年限的延长变化显著,阳坡随封禁年限增加土壤微生物生物量逐渐增加,25 a后微生物生物量碳、氮、磷较封禁前分别增加252%、161%和174%,但显著低于天然侧柏林,仅为其39.0%、41.8%和53.7%;阴坡封禁前10年微生物生物量迅速增加,随后增加幅度减缓,呈波动式缓慢上升趋势;封禁25 a后微生物生物量碳、氮、磷分别增加108%、93%和102%,但明显低于天然杂灌丛群落和辽东栎林,仅为辽东栎林的54.4%、49.1%和40.1%。土壤呼吸强度在封禁5 a后增大明显,且随着年限增加逐渐上升,阳坡25 a时达到最大值,而阴坡15 a时达到最大值,随后开始有所下降,25 a后降至最低点,但仍显著高于放牧地,相同封禁年限的土壤呼吸强度阴坡明显高于阳坡。qCO2随着封禁进程逐渐降低,25 a后达到最低值。相关性分析显示微生物生物量碳、氮、磷、呼吸强度、qCO2与土壤养...

川西亚高山针叶林主要针叶树种年轮生态学研究

米亚罗地区是四川西部较为典型的亚高山针叶林区域之一。为建立该地区主要针叶树种岷江冷杉、云杉、紫果云杉和红杉的年轮宽度年表资料，了解不同海拔高度岷江冷杉原始林和不同恢复过程的人工针叶林及次生混交林树木径向生长规律，结合样地调查，用生长锥钻取了树木芯样做年轮生态学分析。芯样经过标准化程序固定和打磨抛光后，用WinDENDRO图像分析系统测量年轮宽度序列，用COFECHA程序交叉定年和控制测量数据质量，用ARSTAN程序建立了4个主要针叶树种的地区年表和不同海拔高度岷江冷杉林及人工针叶林和次生混交林针叶树的样地年表。 4个主要针叶树种年轮宽度年表的平均敏感度低于0.2，而其晚材宽度年表都具相对较高的平均敏感度。早材宽度与年轮总宽度标准化年表间的相关系数均在0.9以上；晚材宽度与年轮总宽度标准化年表间的相关系数则种间差异较大，红杉的最高，岷江冷杉的最低。岷江冷杉晚材宽度与年轮总宽度的相关性从1970年以后明显下降，而其他种的相关系数则随时间变化较小。树种之间标准化年表显著正相关，而云杉与紫果云杉和红杉与岷江冷杉之间相关系数明显较高。年表序列的第1主分量表达了4个树种树木共同径向生长变化格局；第2至第4主分量分别表达了云杉属和冷杉属、常绿针叶树种和落叶针叶树种以及云杉和紫果云杉树木径向生长变化差异。 不同海拔高度的8个岷江冷杉样地年轮宽度年表序列敏感度大体上随海拔高度升高而降低。各样地早材宽度与年轮总宽度年表之间的相关系数均在0.9以上，且随海拔高度变化不大；晚材宽度与年轮总宽度之间的相关系数随海拔高度的变化较大，并有随海拔升高而降低的趋势。样地年表序列之间相关系数差异很大，高海拔样地年表间多为显著正相关；低海拔样地年表间的相关系数变化不一；高海拔和低海拔样地年表之间相关性较差，且多不显著。样地年表的第1主分量能解释年表序列总方差的37.5%，反映了不同海拔高度岷江冷杉林木共同的径向生长变化格局；第2和第3主分量分别解释总方差的24.5%和18.2%，表现出明显的高海拔和低海拔样地树木间不同的径向生长变化，除一些样地例外，它们一般与低海拔样地年表有正相关，与高海拔样地年表有负相关。在那些另外的样地，海拔以外的其他因素可能也影响了树木径向生长变化。不同海拔高度样地林木的生长抑制和生长释放频率在不同时期表现出较大的差异，表明了不同的干扰历史和林木补充时间。 人工针叶林和次生混交林各样地林木早材宽度与其年轮总宽度年表之间相关系数均高达0.9以上；晚材宽度与年轮总宽度年表之间也都显著正相关，但人工针叶林样地的明显较高。样地年表序列之间的相关关系表现为，林分起源和经营管理相似的样地年表之间的相关系数明显较高，如人工针叶林与人工针叶林尽管树种不同，但样地年表之间显著正相关，而与次生混交林样地年表间关系不显著；反之亦然。综合比较各项生长参数及不同时期的树木径向生长速率，人工针叶林树木的胸径增长至少在40年以内是优于次生混交林的同种（或不同种）针叶树的。不同样地林木生长释放和生长抑制及人工针叶林树木胸高断面积增长分析表明，除严重的人为干扰外，林分郁闭后林木密度过大是造成高频率生长抑制的主要原因，在林分发育的适当时期通过抚育间伐等措施调控林分密度，是保证林木胸高断面积在一定时期内保持较高的连年增长的关键。日本落叶松作为引进的树种，在海拔3100 m左右种植表现良好，近30年来各项生长指标均高于林龄相近的云杉人工林，因此，适当用其作为川西亚高山针叶林采伐迹地快速恢复是合理的。 Miyaluo area is one of the typical regions covered by subalpine coniferous forests in western Sichuan province of southwestern China. To develop the regional tree-ring width chronology series for the dominant conifers such as Abies faxoniana, Picea asperata, P. purpurea and Larix potaninii, and to understand the radial growth patterns of conifers in Abies faxoniana natural forest stands at different altitudes, and in coniferous plantations and natural regenerated mixed stands in their different restoring processes as well, increment cores were sampled in the field together with conventional plots investigations for dendroecological analyses. After the increment cores being prepared according to standard procedures, the ring widths (total-ring and intra-ring widths) were measured with a WinDENDRO image-analysis system, and the measured tree-ring sequences were crossdated and quality-controlled with the software COFECHA. Using the software ARSTAN, we developed tree-ring width based chronology series of the four dominant conifers, eight site-specific Abies faxoniana chronologies, and seven site-specific chronologies of conifers in coniferous plantations and natural regenerated mixed stands. Mean sensitivities for total ring width chronologies of the four sampled dominant conifers were all below 0.2, while those for the latewood width chronologies of the same species were relatively much higher. Correlation coefficients between standard earlywood and total ring width chronologies of the four conifers were all above 0.9, but those between standard latewood and total ring width chronologies exhibited differences among species, with the coefficient of Larix potaninii the highest and that of Abies faxoniana the lowest. Correlation coefficients between latewood and total ring width of A. faxoniana obviously decreased from 1920-1970 for successive 50-year segments with 10-years lag analyses, though the same for the other three species changed unnoticeably with time. Tree-ring standard chronologies among species showed significant positive correlations, with the correlation coefficients between chronologies of Picea asperata and P. purpurea, and of Larix potaninii and Abies faxoniana relatively much higher. The first principal component of tree-ring chronologies represented the common radial growth patterns of the four conifers in Miyaluo area. The second, third and fourth PCs expressed the differences in radial growth responses for the genus Picea and Abies, for the evergreen and deciduous confers, and for the two species of the genus Picea, respectively. In general, mean sensitivities of the eight Abies faxoniana site-specific tree-ring width chronologies decreased with increasing altitude. The correlation coefficients between earlywood and total ring width chronologies for all sites reached 0.9, which did not change much with altitude; but those between latewood and total ring width chronologies diversified, with a decreasing tendency from lower altitudinal sites to higher altitudinal sites. Correlation coefficients among site chronologies varied considerably, with significant positive correlations among higher site chronologies, mixed results among lower site chronologies, and poor and insignificant correlations between chronologies of higher site and lower site. The first PC, which represents 37.5% of the total variance, reflected a common radial growth response at sites of different altitudes, and it showed a tendency of explaining more variance with increasing altitude. The second and the third PCs contributed to 24.5% and 18.2% of the total variance, respectively, exhibiting distinctive differences in radial growth responses between low- and high-altitudinal sites. With some exceptions, the radial growth represented by the second and third PCs had a positive correlation with that at the low-altitudinal sites and a negative correlation with that at the high-altitudinal sites. For those exceptional sites, factors other than altitude might also play a role in affecting tree-ring growth variations. Trees in stands of different altitudes showed great differences in frequencies of growth suppressions and releases through times, suggesting different disturbance histories and periods when trees recruiting to the canopy. Correlation coefficients between earlywood and total ring width chronologies for all sites of coniferous plantations and natural regenerated mixed stands were also above 0.9; and the same between latewood width and total ring width chronologies all positively correlated, too, with the coefficients of the coniferous plantations obviously much higher. Correlations among site chronologies showed that the coefficients among sites with similar stand origin and management regimes were much higher than those among sites with different stand origin and management regimes. For example, significant positive correlations were found for chronologies among different coniferous plantations, irrespective of species differences; while insignificant correlations between chronologies of the same conifer from a coniferous plantation and a natural regenerated mixed stand, and vise versa. Integrative comparisons of different tree growth parameters and radial growth rates at different stages indicated that the diameter at breast height (DBH) increments for trees in coniferous plantations were faster than those for trees of the same (or different) species in the natural regenerated mixed stands, at least within their first 40 years of stand development. Analyses of growth releases and suppressions, and basal area increments of trees in different stands demonstrated that over-dense individuals after canopy closure was the main factor resulting in high frequencies of radial growth suppressions, with some exceptions of severe man-made disturbances. Therefore, to ensure a continuous basal area current annual increment in certain periods, tree density controlling through thinning in due time during the stand development process are necessary. It should be mentioned that, as an introduced conifer to Miyaluo area, Larix kaempferi grew quite well at altitude of ca. 3100 m after planting in 1970s. In their near 30 years of stand development, Larix kaempferi trees exhibited faster growth in various parameters than Picea asperata trees of the similar stand age did. Thus we think it reasonable to use Larix kaempferi as a fast restoring species at appropriate sites of cutting blanks of subalpine coniferous forests in western Sichuan.

遮荫与外源脱落酸喷施对粗枝云杉(Picea asperata Mast.)不同种群抗旱性的影响

本研究通过粗枝云杉不同种群进行的温室半控制试验，采用植物生态学、生理学和生物化学的研究方法，系统地研究了粗枝云杉不同种群抗旱性的生长、形态、生理和生化机理，并结合有关研究进行综合分析，得出主要研究结论如下： 1.粗枝云杉对干旱胁迫的综合反应 粗枝云杉在干旱胁迫下的适应机制为：（1）相对生长速率及植株结构的调整：干旱胁迫下虽然植株相对生长速率显著降低，且有相对较多的生物量向根部分配，但并未发现细根/总根比增加。（2）粗枝云杉对干旱胁迫的光合作用表现为：干旱胁迫显著地降低了控制的理想条件下的气体交换，但干旱胁迫对PSII最大光化学效率（Fv/Fm）没有影响，表明干旱并未影响到光合机构。（3）干旱还影响了很多生理生化过程，包括渗透调解物质（游离脯氨酸）、膜脂过氧化产物、脱落酸（ABA）含量的增加，以及保护酶活性的升高。这些结果证明植物遭受干旱胁迫后发生了一系列的形态、生理和生化响应，这些变化能提高干旱时期植物的存活和生长能力。 2.粗枝云杉不同种群对干旱胁迫反应的种群差异 粗枝云杉三个种群－干旱种群（四川丹巴和甘肃迭部）和湿润种群（四川黑水）对干旱适应不同，这种不同应归因于它们采用的用水策略不同：在水分良好和干旱胁迫条件下，受试种群在相对生长速率和水分利用效率（WUE）方面都表现出显著的种群间差异。与湿润种群相比，干旱种群在两种水分条件下有更高的WUE。粗枝云杉不同种群的碳同位素组分（δ13C）只在干旱胁迫下有显著差异，并且这种差异在水分良好时比干旱胁迫条件下小，说明生理响应和干旱适应性之间的关系受植物内部抗旱机制和外部环境条件（如水分可利用性）或两者互作效应的影响。这些结果说明干旱种群和湿润种群所采用的用水策略不同。干旱种群有更强的抗旱能力，采用的是节水型的用水策略，而湿润种群抗旱能力较弱，采用的是耗水型的用水策略。 3. 遮荫对粗枝云杉不同种群抗旱性影响 干旱胁迫显著降低了全光条件下叶相对含水量（RWC）、相对生长速率、气体交换参数、PSII的有效量子产量（Y），提高了非光化学猝灭效率（qN）、水分利用效率、脯氨酸（PRO）积累、脱落酸（ABA）含量及保护酶活性。然而这种变化在遮荫条件下不明显。我们得出结论适度遮荫降低了干旱对植物的胁迫作用。另一方面，在干旱条件下，与湿润种群相比，干旱种群抗旱性更强，表现在干旱种群净光合速率与单位重量上叶氮含量（Nmass）降低较少。另外，干旱种群表现出更为敏感的气孔导度，更高的热耗散能力（qN）能力、用水效率、ABA积累、保护酶活性，以及更低的总用水量、相对生长速率。这一结果表明这两种群采用不同的生理策略对干旱和遮荫做出反应。许多生长和生理反应差异与这两个种群原产地气候条件相适应。 4. 外源脱落酸(ABA)喷施对粗枝云杉不同种群抗旱性影响 外源ABA喷施在干旱和水分良好条件下均不同程度地提高了根/茎比，表明根和茎对ABA敏感程度不同。实验结果还表明，外源ABA喷施对这两个种群在干旱胁迫期间影响不同。干旱胁迫期间，伴随着ABA喷施，湿润种群净光合速率(A)显著降低，而干旱种群净光合速率变化不明显。另一方面，外源ABA喷施显著提高了干旱条件下干旱种群的单位叶面积重（LMA）、根/茎比、细根/总根（Ft）比、水分利用效率（WUE）、ABA含量, 以及保护酶活性。然而，外源ABA喷施对湿润种群的上述测定指标没有显著影响。这一结果表明干旱种群对外源ABA喷施更为敏感, 反应在更大的气孔导度降低，更高的生物量可塑性，及更高的水分利用效率、ABA含量和保护酶活性。综上所述，我们得出结论，粗枝云杉对外源ABA敏感性因种群的不同而不同。该研究结果可为两个明显不同种群在适应分化方面提供强有力的证据。 Arid or semi-arid land covers more than half of China's land territory. In arid systems, severe shortages of soil water often coincide with periods of high temperatures and high solar radiation, producing multiple stresses on plant performance. Protection from high radiation loads in shaded microenvironments during drought may compensate for a loss of productivity due to reduced irradiance when water is available. Additionally, ABA, a well-known stress-inducible plant hormone, has long been studied as a potential mediator for induction of drought tolerance in plants. Picea asperata Mast., which is one of the most important tree species used for the production of pulp wood and timber, is a prime reforestation species in western China. In this experiment, different population of P. asperata were used as experiment material to study the adaptability to drought stress and population differences in adaptabiliy, and the effects of shade and exogenous abscisic acid (ABA) application on the drought tolerance. Our results cold provide a strong theoretical evidence and scientific direction for the afforestation, and rehabilitation of ecosystem in the arid and semi-arid area, and provide a strong evidence for adaptive differentiation of different populations, and so may be used as criteria for species selection and tree improvement. The results are as follows: 1. A large set of parallel response to drought stress Drought stress caused pronounced inhibition of the growth and increased relatively dry matter allocation into the root; drought stress also caused pronounced inhibition of photosynthesis, while drought showed no effects on the maximal quantum yield of PSII photochemistry (Fv/Fm) in dark-adapted leaves, indicating that drought had no effects on the primary photochemistry of PSII. However, in light-adapted leaves, drought reduced the quantum yield of PSII electron transport (Y) and increased the non-photochemical quenching (qN). Drought also affected many physiological and biochemical processes, including increases in superoxide dismutase (SOD), ascorbate peroxidase (APX) activities, malondialdehyde and ABA content. These results demonstrate that there are a large set of parallel changes in the morphological, physiological and biochemical responses when plants are exposed to drought stress; these changes may enhance the capability of plants to survive and grow during drought periods. 2. Difference in adaptation to drought stress between contrasting populations of Picea asperata There were significant population differences in growth, dry matter allocation and water use efficiency. Compared with the wet climate population (Heishui), the dry climate population (Dan ba and Jiebu) showed higher LMA, fine root/total root ratio and water use efficiency under drought-stressed treatments. The results suggested that there were different water-use strategies between the dry population and the wet population. The dry climate population with higher drought tolerance may employ a conservative water-use strategy, whereas the wet climate population with lower drought tolerance may employ a prodigal water-use strategy. These variations in drought responses may be used as criteria for species selection and tree improvement. 3. The effects of shade on the drought tolerance For both populations tested, drought resulted in lower needle relative water content (RWC), relative growth rate (RGR), gas exchange parameters and effective PSII quantum yield (Y), and higher non-photochemical quenching (qN), water use efficiency (WUE), proline (PRO) and abscisic acid (ABA) accumulation, superoxide dismutase (SOD), ascorbate peroxidase (APX) activities as well as malondialdehyde (MDA) levels and electrolyte leakage in sun plants, whereas these changes were not significant in shade plants. Our study results implied that shade, applied together with drought, ameliorated the detrimental effects of drought. On the other hand, compared with the wet climate population, the dry climate population was more tolerant to drought in the sun treatment, as indicated by less decreases in A and mass-based leaf nitrogen content (Nmass), more responsive stomata, greater capacity for non-radiative dissipation of excitation energy as heat (analysed by qN), and higher WUE，higher level of antioxidant enzyme activities，higher ABA accumulation as well as lower MDA content and electrolyte leakage. Many of the differences in growth and physiological responses reported here are consistent with the climatic differences between the locations of the populations of P. asperata. 4. The effects of exogenous abscisic acid (ABA) application on the drought tolerance For both populations tested, exogenous ABA application increased root/shoot ratio (Rs) under well-watered and drought-stressed conditions, indicating that there was differential sensitivity to ABA in the roots and shoots. However, it appeared that ABA application affected the two P. asperata populations very differently during drought. CO2 assimilation rate (A) was significantly decreased in the wet climate population, but only to a minor extent in the dry climate population following ABA application during soil drying. On the other hand, ABA application significantly decreased stomatal conductance (gs), transpiration rate (E) and malondialdehyde (MDA) content, and significantly increased leaf mass per area (LMA), Rs, fine root/total root ratio (Ft), water use efficiency (WUE), ABA contents, superoxide dismutase (SOD), ascorbate peroxidase (APX) and catalase (CAT) activities under drought condition in the dry climate population, whereas ABA application did not significantly affect these parameters in the wet population plants. The results clearly demonstrated that the dry climate population was more responsive to ABA application than the wet climate population, as indicated by the strong stomata closure and by greater plasticity of LMA and biomass allocation, as well as by higher WUE, ABA content and anti-oxidative capacity to defense against oxidative stress, possibly predominantly by APX. We concluded that sensitivity to exogenous ABA application is population dependent in P. asperata. Our results provide strong evidence for adaptive differentiation between populations of P. asperata.