土壤侵蚀危险度评价研究进展

通过对已有土壤侵蚀危险度评价研究的分析,将这些评价方法归纳为土壤抗蚀年限评价法、土壤侵蚀危险度主导因子分级评价法和土壤侵蚀危险度模型评价法3大类,并在对各种评价方法综合分析的基础上,对土壤侵蚀危险度评价的研究方向和发展趋势进行讨论。建议今后加强土壤侵蚀危险度定量评价方法,尤其是模型方法的研究。

侵蚀环境人工灌木林土壤活性有机碳与碳库管理指数演变

采用时空互代法,以黄土丘陵区纸坊沟流域生态恢复过程中不同年限的人工柠条和沙棘林为研究对象,选取坡耕地和天然侧柏林为对照,分析了植被恢复过程中土壤有机碳(TOC)、活性有机碳(LOC)、非活性有机碳(NLOC)及碳库管理指数的演变特征。结果表明,侵蚀环境下的坡耕地由于人为干扰,土壤碳库含量偏低,退耕营造柠条林可以显著增加土壤碳库各组分含量,并随恢复年限呈显著线性关系,25 a时TOC、LOC和NLOC分别较坡耕地增加271%、144%和204%,仅为侧柏林的32%、30%和29%,碳库指数和碳库管理指数较坡耕地明显增加,增幅分别达到144%和108%,仅为侧柏林的28%和43%;不同灌木林对土壤碳库管理的改善作用不同,恢复年限相同的沙棘林土壤碳库组分含量和管理指数明显高于柠条林,坡耕地营造灌木林后土壤经营和管理水平得到了显著改善,土壤系统向着良性方向转变。相关性分析表明有机碳、活性有机碳、非活性有机碳、碳库指数、碳库管理指数与土壤主要肥力因子相关性极其密切,可以作为反映生态恢复过程土壤质量演变的指标。

模拟降雨下初始含水量对砂黄土硝态氮迁移特征的影响

利用室内人工模拟降雨,研究了不同初始含水量砂黄土在降雨条件下入渗-径流、土壤侵蚀,以及NO3--N随径流流失和土壤深层淋溶特征。结果表明,初始含水量对产流时刻的影响在相对含水量为49.4%和76.9%之间存在一个转折点,高初始含水量较低含水量产流提前大约15 min;土壤侵蚀量随着土壤初始含水量的增加而增加,相对含水量为97.1%时,侵蚀泥沙量分别是相对含水量22.9%时的2.8倍,49.4%时的2.3倍,76.9%时的1.5倍。初始含水量高的处理径流初始NO3--N浓度高,随后各处理均衰减很快,10 min左右NO3--N含量趋于雨水本底值;土壤初始含水量越低,NO3--N被淋洗的程度越严重,土壤剖面中NO3--N的浓度峰越深。对于黄土高原坡地砂黄土NO3--N迁移特征来看,按照NO3--N迁移数量,随径流和泥沙流失量比向土壤深层迁移的数量小。说明在降雨条件下,NO3--N主要通过土壤深层淋溶损失,且土壤初始含水量越低其损失越严重。针对黄土高原降水量小、分布集中的特点,采取措施增加入渗,蓄积水分,在一定含水量下施肥,以提高氮肥利用率,降低NO3--N的淋溶。

分布式水文模型的误差分析

分布式水文模型以其具有明确物理意义的参数结构和对空间分异性的全面反映,能够准确详尽地描述和模拟流域内真实的降水径流过程而被广泛需求和关注。在模拟土地利用、土地覆盖、水土流失等各种变化过程的水文响应,面源污染、陆面过程、气候变化影响评价等诸多领域都有广泛的应用。模型的预报精度和误差至关重要,决定了模型的应用和推广。在分析分布式水文模型建立和验证过程的基础上,提出了模型的4类误差来源:被排除在外的因素引起的误差,实测历史记录资料的随机或系统误差,参数误差和模型结构误差,讨论了各类误差的分析与计算方法,为模型的发展和成长提供了依据。

侵蚀环境人工刺槐林土壤水稳性团聚体演变及其养分效应

采用时空互代法,以典型侵蚀环境纸坊沟流域生态恢复过程中不同年限的人工刺槐林为研究对象,选取坡耕地和天然侧柏林为参照,分析了植被恢复过程中土壤水稳性团聚体变化规律及其与土壤养分状况之间的相互关系。结果表明,侵蚀环境下的坡耕地由于人为耕作干扰,土壤水稳性团聚体含量低下,抗蚀性能较差。营造刺槐林前5 a土壤水稳性团聚体含量较坡耕地显著快速增加,随后增幅变缓,成对数增长,恢复25 a后土壤水稳性团聚体含量已经达到天然侧柏林水平。在植被恢复过程中土壤小粒径的水稳性团聚体逐渐聚集转变形成大粒径的团聚体。相关性分析和回归分析表明,植被恢复过程中>0.25 mm土壤水稳性团聚体与有机碳、全氮、碱解氮、速效钾、容重等相关性达到显著(p<0.05)或极显著水平(p<0.01),与全磷和速效磷相关性较弱。坡耕地退耕营造刺槐林后可以减少人为干扰,增加碳素和氮素供给,提高水稳性团聚体含量,使土壤抗蚀性能提高。

自然侵蚀量和容许土壤流失量与水土流失治理标准

在总结分析国内有关自然侵蚀量、容许土壤流失量研究成果的基础上,探讨自然侵蚀量与容许土壤流失量的关系,讨论制订水土流失治理标准的思路。认为水土流失治理标准的确定有3个参考值:1)标准值,即一定条件下的容许土壤流失量,是水土流失治理至少要达到的目标,且随着对不同土地利用类型的功能需求与可实施的最佳水土保持措施、以及所在水土流失类型区的侵蚀危害与治理约束条件的不同而不同;2)理想值,即正常自然侵蚀状态下的土壤流失量;3)极端值,即土壤流失量为0,不发生水土流失。水土流失治理应是先控制到现时生态环境与社会经济条件下的容许土壤流失量范围内,逐步达到自然正常侵蚀量或制止水土流失的发生;还应引入环境伦理、环境美学及景观设计等理念,最终实现土地的可持续利用、区域生态系统的健康稳定及人与自然的和谐友好发展。

黄土坡面水力学特征参数与土壤侵蚀量间关系研究

通过室内人工降雨模拟试验,研究黄土坡面水流的水力学特征参数与土壤侵蚀量间的关系。结果表明:随着放水流量的增加,出口断面平均流速、雷诺数在逐渐增大,弗劳德数在逐渐减小,坡面流流态也从层流状态逐渐过渡到紊流状态,当放水流量在10.1~19.9 L/min范围内时,坡面水流的流态为层流与紊流之间的混和状态;出口断面的平均流速随着放水流量的增加而增大,二者之间呈良好的幂函数关系;黄土坡面的累计泥沙侵蚀总量随着雷诺数的增加而增大,二者呈幂函数关系;单宽土壤侵蚀产沙量与雷诺数之间拟合关系的好坏受雷诺数取值的影响,当雷诺数的值在临界雷诺数(取500)附近时,单宽土壤侵蚀产沙量与雷诺数之间的负相关性良好。研究结果不仅可以揭示土壤侵蚀水动力过程的内在机制,而且对黄土高原地区生态环境建设具有重要的现实意义。

有限元子域法分析土石混合介质的饱和导水率

在黄土高原水蚀风蚀交错地区,碎石的存在影响着土壤的水力学特性,因此研究土石混合介质的渗流运动对该地区的植被恢复建设具有重要的指导意义。虽然用有限元法分析地下水渗流的工作很多,但对土石混合介质而言,由于有限元法的解题规模迅速增加,其研究甚少。本模型采用子域法矩形单元,假定土石体中石块含量增加不影响土的孔隙率且石块不透水分析了土石混合介质中的渗流过程。结果表明,影响土石混合介质饱和导水率的主要因素是土石介质的含石率;石块的大小基本不影响混合介质饱和导水率;石块分布方式增加水流路径时,会影响混合介质饱和导水率。这一数值分析结论与土石体渗流问题的部分理论结果和试验结果吻合良好。本模型子域法计算效率高,结果较准确,可用于较大规模的土石混合介质渗流分析,还可推广用于土石混合介质的三维分析和多级子域法有限元技术以获得更大规模的解题效果。

景观格局对土壤侵蚀的影响研究-以岷江上游的小流域为例

岷江上游地区生态环境脆弱，加之不合理的土地利用，土壤侵蚀问题日趋严重。准确定量土壤侵蚀、产沙模数的空间分布，从而分析土地利用及其空间配置对于土壤侵蚀过程是否合理？对减少长江上游的输沙量具有重要意义，是一个严峻的研究课题。本文以GIS为平台，建立泥沙输移分布模型，模拟流域的侵蚀模数、泥沙输移比和产沙模数的空间分布；结合模型模拟结果，分析不同土地利用/覆被方式的侵蚀和产沙特征，并以景观生态学的“格局-过程”理论为基础，通过景观指数——景观空间负荷对比指数，探讨土地利用/覆被的空间格局对土壤侵蚀过程的影响。泥沙输移分布模型模拟的产沙量与实测值比较，结果十分理想。模拟结果表明：两个流域大部分区域属于微度到中度侵蚀，平均侵蚀模数略低于长江上游的平均值；两个流域的平均泥沙输移比均远小于长江上游的平均值；河流两侧1 km范围内是流域主要的产沙区域；不同土地利用/覆被方式的侵蚀模数不同——裸岩>居民点>草地 >农田>灌木林>林地；侵蚀广泛分布在草地、林地和灌木林地上，占流域侵蚀总量的80%以上；灌木林和林地是流域主要的产沙源，占总产沙量的60%以上。景观空间负荷对比指数的计算结果表明：两个流域土地利用/覆被随坡度的空间配置有利于水土保持，但随相对距离、相对高度和运移距离的空间配置不利于水土保持；两个流域“源”大“汇”小的土地利用/覆被组成方式有利于水土保持，且镇江关流域的组成方式更优；景观空间负荷对比指数能够比较全面客观地反映景观格局对对土壤侵蚀过程的影响：两个流域的景观格局均较为合理；且镇江关流域的景观格局优于黑水流域。

森林流域单元水文模型的构建与应用

目前全球缺水、水污染、洪涝灾害以及水土流失仍然非常严重，尤其在我国北方地区。流域水文模型可用来进行不同需水管理的情景分析，为解决我国水问题提供科学依据。分布式水文模型是流域水文模型的发展方向，具有显著特点：1）应用前景广泛，不仅可以模拟流域水文过程，还可以协助模拟泥沙或污染物的运移过程，为水利工程设计、水土保持、环境保护等领域提供技术支持；2）能够预测流域土地利用或气候变化下的流域水文响应过程变化，为管理部门提供决策支持；3）模型所需要的参数全部具有物理意义，可通过实际测量确定，适合模拟实测系列较短或是无观测流域的水文过程；4）对于目前国际水文界的前沿问题—水文尺度转换提供了一种有效的解决途径。 然而分布式水文模型还不完善，如1)真实性问题。对一些水文过程和边界条件还不确定。2)尺度转换问题。目前很少考虑尺度对参数有效性的影响。3)检验问题。还无法判断对有些难以测量的水文状态变量的模拟正确与否。4)计算时间和数据存储的问题。有些分布式水文模型虽然具有很强的水文物理基础和完善的模型结构，但是计算时间过长和(或)数据存储过大，难以应用。上述问题的核心就是对分布式水文模型的核心—单元水文模型的研究不够，需要为进一步完善单元水文模型进行研究。 本文采用饱和入渗理论、Saint-Venant方程、Richards方程、Penman-Monteith方程等等构建了以有限差分法求解的适用于森林流域的单元水文模型，并通过实验室模拟试验和坡地径流场资料进行了验证，主要结论为： 通过不同坡度和不同雨强下的室内坡面产汇流实验模拟，表明：该模型模拟的坡面流和壤中流过程与实测过程基本一致，峰现时间、径流历时、峰值流量、出流总量模拟值与实测值的相对误差均较小，基本小于10%。模型的模拟精度较高，实用性较强，为深入研究壤中流机制和改进流域降雨-径流模型提供了理论依据。 通过坡地径流观测场实测资料的验证，表明：该模型模拟的坡面流过程精度较高，累计流量的精度更高于小时过程的精度，离差系数、效率系数、确定系数均较理想，具有应用价值，有助于改善分布式水文模型在森林流域的模拟效果。

Initiation and development of water film by seepage

When water seeps upwards through a saturated soil layer, the soil layer may become instability and water films occur and develop. Water film serves as a natural sliding surface because of its very small friction. Accordingly, debris flow may happen. To investigate this phenomenon, a pseudo-three-phase media is presented first. Then discontinuity method is used to analyze the expansion velocity of water film. Finally, perturbation method is used to analyze the case that a water flow is forced to seep upwards through the soil layer while the movement of the skeleton may be neglected relative to that of water. The theoretical evolutions of pore pressure gradient, effective stress, water velocity, the porosity and the eroded fine grains are obtained. It can be seen clearly that with the erosion and re-deposited of fine grains, permeability at some positions in the soil layer becomes smaller and smaller and, the pore pressure gradient becomes bigger and bigger, while the effective stress becomes smaller and smaller. When the effective stress equals zero, e.f. liquefaction, the water film occurs. It is shown also that once a water film occurs, it will be expanded in a speed of (U)(t)/(1 - E >).

A study on water film in saturated sand

Water film can serve as a sliding surface and cause landslides on gentle slopes. The development of "water film" in saturated sand is analyzed numerically and theoretically based on a quasi-three-phase model. It is shown that stable water films initiate and grow if the choking state (where the fluid velocity decreases to near zero) remains steady in a liquefied sand column. Discontinuity can occur in pore water velocity, grain velocity and pore pressure after the initiation of a water film. However, the discontinuity and water film can disappear once the choking state is changed. The key to the formation of water film is the choking in the sand column caused by eroded fine grains.

Plasma and electrode emissions from a 1 kW hydrogen-nitrogen arcjet thruster

Arc root behavior affects the energy transfer and nozzle erosion in an arcjet thruster. To investigate the development of arc root attachment in 1 kW class N2 and H2-N2 arcjet thrusters from the time of ignition to the stably working condition, a kinetic series of end-on view images of the nozzle obtained by a high-speed video camera was analyzed. The addition of hydrogen leads to higher arc voltage levels and the determining factor for the mode of arc root attachment was found to be the nozzle temperature. At lower nozzle temperatures, constricted type attachment with unstable motions of the arc root was observed, while a fully diffused and stable arc root was observed at elevated nozzle temperatures.

Sand storms: CFD analysis of Reynolds stress and collision stress of particles near sand bed

Sand storm is a serious environmental threat to humans. Sand particles are transported by saltation and suspension, causing soil erosion in one place and deposition in another. In order to prevent and predict sand storms, the causes and the manners of particle motions must be studied in detail. In this paper a standard k-epsilon model is used for the gas phase simulation and the discrete element method (DEM) is used to predict the movements of particles using an in-house procedure. The data are summarized in an Eulerian-Eulerian regime after simulation to get the statistical particle Reynolds stress and particle collision stress. The results show that for the current case the Reynolds stress and the air shear stress predominate in the region 20-250 mm above the initial sand bed surface. However, in the region below 3 mm, the collision stress must be taken into account in predicting particle movement. (C) 2010 Chinese Society of Particuology and Institute of Process Engineering, Chinese Academy of Sciences. Published by Elsevier B.V. All rights reserved.

中国卧龙自然保护区不同海拔川滇高山栎(Quercus aquifolioides)群体的遗传变异

海拔梯度造成的环境异质性，如崎岖的地形、复杂的植被结构以及花期延迟等可能会极大地影响到物种的形态和遗传变异格局。理解物种形态和遗传变异的海拔格局对于物种多样性的管理和保护是非常重要的。尽管植物群体遗传学是一个飞速发展的研究领域，然而与海拔相关的形态变异、遗传变异及群体间遗传差异的研究却很少。到目前为止，还不清楚遗传变异与海拔之间是否必然的相关性。 川滇高山栎是一种重要的生态和经济型树种，广泛分布于中国西南的四川、西藏、贵州和云南省的高海拔地区，在保持水土、调节气候方面起着十分重要的作用。尽管主要受阳光限制而仅分布于阳坡，但其海拔梯度范围较大，表明川滇高山栎对不同的环境具有很强的适应性。本文通过叶型及生理响应、微卫星分子标记和扩增性片段长度多态性方法，试图探索川滇高山栎叶沿海拔梯度的形态和生理响应及其沿海拔梯度的遗传变异格局，为川滇高山栎的保护和利用提供进一步的遗传学理论依据和技术指导。 对叶形、含氮量及碳同位素的试验结果表明，平均比叶面积、气孔密度、气孔长度和气孔指数等气孔参数随海拔的升高呈非线性变化。在海拔大于2800 m时，川滇高山栎的比叶面积、气孔长度和气孔指数都随海拔升高而降低，但是在海拔小于2800 m时，这些指标都随海拔的升高而增大。相对而言，单位叶面积的含氮量和碳同位素则表现出相反的变化模式。另外，比叶面积是决定碳同位素沿海拔梯度变化的最重要参数。本研究结果表明，海拔2800 m附近是川滇高山栎生长和发育的最适地带，在这里生长的植物叶片厚度更薄、气孔更大、叶碳同位素值更小。 利用六对微卫星引物对五个不同海拔川滇高山栎群体遗传多样性进行研究，结果表明，群体内表现出较高的遗传多样性，平均每位点等位基因数11.33个，平均期望杂合度达0.820。群体间差异较小，分化仅为6.6％。聚类分析也并没有显示出明显的海拔格局。然而低频率等位基因却与海拔呈显著性正相关(R2=0.97, P < 0.01)，表明在高海拔处，川滇高山栎以更多的稀有基因来适应恶劣的环境条件。本试验结果表明由海拔梯度形成的选择性压力对川滇高山栎群体的遗传变异影响并不明显。 为了进一步探讨川滇高山栎群体遗传变异与海拔之间的相互关系，我们还对其进行了扩增性片段长度多态性分析。结果表明：(1)随海拔的升高(从群体WL2到群体WL5)，群体内遗传变异降低，而群体间遗传差异增加；(2)低海拔群体WL1表现出最低的遗传变异性(HE = 0.181)，同时与其余四个群体间呈现出最大的遗传差异性(平均FST = 0.0596)；(3)在除去低海拔群体WL1后，Mantel检测表明群体间遗传距离与海拔距离之间表现出正相关性。另外，研究结果还表明，遗传变异受生境条件(过度的湿热环境)及人为干扰(火烧、砍伐和放牧)的影响，这一点至少在低海拔群体WL1上发生了作用。 通过叶形态、生理及DNA分子水平的研究，结果表明叶形态特征和碳同位素与海拔紧密相关，与海拔之间呈非线性变化，海拔2,800 m附近是川滇高山栎生长和发育的最适地带。海拔梯度在一定程度上会影响到川滇高山栎群体的遗传变异结构，但在这样一个狭窄的地理分布区域里，这种影响并不足以导致群体间较大的遗传分化。同时生境条件及人为干扰也是影响遗传变异的限制性因子，不容忽视。 Altitudinal gradients impose heterogeneous environmental conditions, such as rugged topography, a complex pattern of vegetation and flowering delay, and they likely furthermore markedly affect the morphological and genetic variation pattern of a species. Understanding altitudinal pattern of morphological and genetic variation at a species is important for the management and conservation of species diversity. Although plant population genetics is a fast growing field of research, there are only few recent investigations, which analyzed the genetic differentiation and changes of intra-population variation along altitudinal gradients. At present, it is still unclear whether there are some common patterns of morphological and genetic variation with altitude. Quercus aquifolioides Rehder & E.H. Wilson, which is an important ecological and economical endemic woody plant species, is widely distributed in the Yunnan and Sichuan provinces, Southwest China. Its large range of habitat across different altitudes implies strong adaptation to different environments, although it is mainly restricted to sunny, south facing slopes. It plays a very important role in preventing soil erosion, soil water loss and regulating climate, as well as in retaining ecological stability. In this paper, we tried to understand the altitudinal pattern of morphological and genetic variation along altitudinal gradients through the experiments of leaf morphological and physiological responses, microsatellite analysis and AFLP markers. In leaf morphological and physiological responses experiment, we measured leaf morphology, nitrogen content and carbon isotope composition (as an indicator of water use efficiency) of Q. aquifolioides along an altitudinal gradient. We found that these leaf morphological and physiological responses to altitudinal gradients were non-linear with increasing altitude. Specific leaf area, stomatal length and index increased with increasing altitude below 2,800 m, but decreased with increasing altitude above 2,800 m. In contrast, leaf nitrogen content per unit area and carbon isotope composition showed opposite change patterns. Specific leaf area seemed to be the most important parameter that determined the carbon isotope composition along the altitudinal gradient. Our results suggest that near 2,800 m in altitude could be the optimum zone for growth and development of Q. aquifolioides, and highlight the importance of the influence of altitude in research on plant physiological ecology. Genetic variation and differentiation were investigated among five natural populations of Q. aquifolioides occurring along an altitudinal gradient that varied from 2,000 to 3,600 m above sea level in the Wolong Natural Reserve of China, by analyzing variation at six microsatellite loci. The results showed that the populations were characterized by relatively high intra-population variation with the average number of alleles equaling 11.33 per locus and the average expected heterozygosity (HE) being 0.779. The amount of genetic variation varied only little among populations, which suggests that the influence of altitude factors on microsatellite variation is limited. However, there is a significantly positive correlation between altitude and the number of low-frequency alleles (R2=0.97, P < 0.01), which indicates that Q. aquifolioides from high altitudes has more unique variation, possibly enabling adaptation to severe conditions. F statistics showed the presence of a slight deficiency of heterozygosity (FIS=0.136) and a low level of differentiation among populations (FST=0.066). The result of the cluster analysis demonstrates that the grouping of populations does not correspond to the altitude of the populations. Based on the available data, it is likely that the selective forces related to altitude are not strong enough to significantly differentiate the populations of Q. aquifolioides in terms of microsatellite variation. To further elucidate genetic variation pattern of Q. aquifolioides populations under sub-alpine environments, genetic variation and differentiation were investigated along altitudinal gradients using AFLP markers. The altitudinal populations with an average altitude interval of 400 m, i.e. WL1, WL2, WL3, WL4 and WL5, correspond to the altitudes 2,000, 2,400, 2,800, 3,200 and 3,600 m, respectively. Our results were as follows: (i) decreasing genetic variation (ranging from 0.253 to 0.210) and increasing genetic differentiation with altitude were obtained from the WL2 to the WL5 population; (ii) the WL1 population showed the lowest genetic variation (HE = 0.181) and the highest genetic differentiation (average FST = 0.0596) with the other four populations; (iii) the positive correlation was obtained using Mantel tests between genetic and altitude distances except for the WL1 population. Our results suggest that altitudinal gradients may have influenced the genetic variation pattern of Q. aquifolioides populations to some extent. In addition, habitat environments (unfavorable wet and hot conditions) and human disturbances (burning, grazing and felling) were possible influencing factors, especially to the low-altitude WL1 population. The present study shows that there were close correlations between morphological features and carbon isotope composition in our data. This indicates that a coordinated plant response modified these parameters simultaneously across different altitudes. Around 2,800 m altitude there seems to be an optimum zone for growth and development of Q. aquifolioides, as indicated by thinner leaves, larger stomata and more negative d13C values. All available evidence indicates altitudinal gradients may have influenced the genetic variation pattern of Q. aquifolioides to some extent. Decreasing genetic variation and increasing genetic differentiation with altitude was obtained except for the WL1 population. And the environment of habitats and human disturbances were also contributing factors, which impact genetic variation pattern, especially to the low-altitude WL1 population.