植被-侵蚀动力学模型参数确定及在黄土高原应用.docVIP

植被-侵蚀动力学模型参数的确定及在黄土高原的应用 王费新，王兆印1，2 * 1，清华大学水利系，北京 100084；2，国际泥沙研究培训中心，北京 100044 植被侵蚀动力学是研究流域植被与侵蚀在人类活动影响下演变规律的一门新的边缘学科。由植被-侵蚀动力学模型可绘制植被-侵蚀状态图，以预测停止人为干扰后植被及侵蚀变化趋势，评价水土保持措施成效，提出小流域优化治理方略。植被-侵蚀动力学模型4个参数a、c、b、f是流域重要特征参数，其值的确定是进行植被-侵蚀动力学模拟及绘制植被-侵蚀状态图的前提和基础。它们依赖于流域气候、土壤和地形地貌。本文以黄土高原典型小流域为例，采用试算法确定了它们的取值，通过相关分析发现参数a与P、T，参数c、b、f与D50、S之间存在明显的相关关系，并得到了相关关系式。将其用于吕二沟小流域植被-侵蚀动力学模拟，结果与实际情况符合较好。绘制了各参数的取值分布图。最后以燕儿沟小流域为例，将以上相关关系式用于植被-侵蚀状态图的绘制，及小流域自然状况、植被-侵蚀状态及水土保持效果的评价。 植被侵蚀动力学；黄土高原；模型参数；小流域；相关性S157 文献标识码：A 文章编号：1672-2175（2006）06-1366-06王兆印等[1]于2003年创建了一个新的边缘科学——植被侵蚀动力学，用于研究流域植被与侵蚀在人类活动影响下演变规律，建立了植被侵蚀动力学模型，并成功应用于云南小江及其子流域植被发育及土壤侵蚀演变过程的模拟。此后，王兆印等[2-5]又将其推广到北京西山地区、甘肃定西安家沟及山西离石王家沟等小流域，模拟结果与小流域实测资料符合情况较好；并在此基础上开发了植被-侵蚀状态图这一评价工具，用于探讨植被和侵蚀在没有或停止人类干扰作用后的演变趋势，该工具能较客观地评价各典型小流域植被侵蚀状况和水土保持成效；此外还初步探讨了植被侵蚀动力学参数与流域自然特征的关系。 植被-侵蚀动力学模型中的4个参数是植被-侵蚀动力学理论的重要组成部分。由这4个参数可以模拟和预测小流域在一定水土保持措施下植被发育和土壤侵蚀变化过程；可以做出植被-侵蚀状态图，用于分析小流域自然特性及植被-侵蚀现状，评价各种水土保持措施效果，提出优化治理方略，指导小流域展开综合治理工作[4]。这4个参数取值依赖于流域气候、土壤和地形地貌，是地区的重要特征值，与自然应力和人类活动无关[4]。由于问题的复杂性，目前它们主要采用试算法得到，限制了植被-侵蚀动力学在资料相对缺乏的小流域和地区的应用。本文拟以黄土高原为例，选取典型小流域，通过试算法确定这4个参数的取值，并通过相关性分析得到它们与流域气候、地形地貌及土壤特性之间的定量化关系式，完善植被-侵蚀动力学理论，并推广植被-侵蚀动力学模型和植被-侵蚀状态图在黄土高原及其它水土流失地区的应用。 植被侵蚀动力学模型由耦合的植被-侵蚀动力学方程组组成[2]，其形式如下： 其中V为植被覆盖度，E为侵蚀率，为发生在时刻的致死生态应力，代表人类活动对植被覆盖度的影响，为人类活动对侵蚀率的影响，a、c、b、f即为植被侵蚀动力学模型的4个特性参数，其中，a代表植被涵养水分及繁殖作用下植被自我增长作用，c代表侵蚀对植被的破坏作用，b代表土壤侵蚀自身的促进作用，f代表植被对侵蚀的控制作用。 利用植被-侵蚀动力学方程组可以做出植被-侵蚀状态图[4]。令方程组(1)中各生态应力项为零，得到齐次微分方程组： 其中：，。分别令方程组(2)中及，得到植被覆盖度变化率和侵蚀模数变化率为零的2条直线： 在，的坐标平面上绘制这2条直线及植被-侵蚀状态点，即可得到植被-侵蚀状态图。这2条直线将V-E平面划分成3个区，分别为A、B、C区或者A、D、C区。不同区植被和侵蚀模数变化趋势各不相同：A区，，在没有人类干扰情况下，处于该区的小流域会自动向着侵蚀模数增大、植被覆盖度降低、生态环境恶化的方向发展，为不良状态区；C区与A区相反，，，处于该区的小流域自动向着侵蚀模数降低、植被完善、生态环境改善的方向发展，为良性状态区；B区，，植被和侵蚀模数都增长，为不稳定的过渡区，经过一定时期可能会进入A区或C区；D区，，植被和侵蚀模数都减小，亦为不稳定的过渡区 图1 王家沟植被-侵蚀状态图(●-治理前；▲-治理后). Fig. 1 Vegetation-erosion chart for the Wangjiagou Watershed (●- before reclamation project; ▲- after reclamation project). 图2黄土高原中径变化及典型小流域所在县市分布图(改自文献[7])Fig. 2 Medial diameters variation and the cities d