GISECOlecture5ext01.ppt

地理信息系统在资源环境研究中的应用 宋现锋 中国科学院研究生院资源与环境学院 Phone: 010-8825-6462 Email: xfsong@ 4 USLE参数的其他标定方法 1. R的简单计算方法 E的数值计算需要详细的降雨过程，即分析自记雨量纸。问题所在？ 由于计算E的困难，R的简易计算方法寻求使用常规降雨资料得到R参数，省去动能E的计算。 针对不同时间分辨率的降雨数据，建立区域性的经验方程。 A) 次降雨R值的计算 国内外研究表明，雨量(P)与动能(E)之间存在很好的线性关系。建立EI30与PI30的相关关系，通过PI30求得EI30。 R = 1.70 (P I30 / 100) – 0.136 (I30 10 mm/h) 相关系数r=0.991; 剩余标准差s=0.2; 相对误差8.0%; n=168 R = 2.35(P I30 / 100) –0.523 (I30 = 10mm/h) 相关系数r=0.989; 剩余标准差s=4.8; 相对误差13.7%; n=766 R为次降雨侵蚀力(100 m ton cm / hm2 h)，P为次降雨量(mm)，I30为最大30分钟降雨强度(mm/h) B) 年降雨R的计算 国际上，有两种途径：一是采用次降雨R值简易计算的基本形式，通过年雨量(Py)或汛期雨量(Pf)与年最大30分钟雨强I30的乘积形式，估算年R值；二是采用月雨量与年雨量模比系数的方法估算R值。 在北方干旱、半干旱地区，雨强因子对土壤流失量的影响比雨量因子大得多，年雨量/月雨量都未明显包含强度概念，这与EI30结构的物理意义不相接近。 R = 1.67 (P10 I60 / 100)0.953 P10为一年内次降雨量大于10mm的总和(mm) R = 0.272(P I60 / 100)1.205, 缺乏P10资料时。 R为年侵蚀力(100 m ton cm / hm2 h a) P为年降雨总量(mm) I60为年最大60分钟降雨量(mm) 使用北京等12站在1956-1983年系列降雨资料，进行相关统计分析得到上述公式。平均误差为15%左右。 C) 多年降雨R的计算 类似年R的估算方法，应用北京等25个站多年平均资料，选择四个特征值建立多年平均降雨侵蚀力的估算公式： R = 0.032 P100.616 I601.394 R = 0.009 P0.564 I601.155 I14400.560 缺乏P10资料时 R为年侵蚀力(100 m ton cm / hm2 h a) P10为次降雨大于10mm的年平均降雨量(mm) I60为年最大60分钟降雨量(mm) I1440为年最大1440分钟降雨量(mm) D) 西非洲年降雨R的计算 R = A x P R为年侵蚀力(173.6 m ton cm / hm2 h a) ??? P为年平均降雨总量(mm) (10年) A – 0.5 in most cases ?0.05 0.6 near the sea (40km) 0.2~0.3 in tropical mountain areas 0.1 in Mediterranean mountain areas E) 中国R值的时空分布 表格参见 景可p75-83 2 我国的K值获取： 一是应用诺谟图。土壤质地分类系统的转换: 苏制--美制(美制计算基础小于2.0mm, 苏制计算基础小于1.0mm)。具体方法：图解法 – 先在半对数纸上画出苏制的土壤颗粒大小分配曲线，然后查出美制某一粒径的百分含量。 二是实测K值。张宪奎在东北黑土地区的实测发现，诺谟图求出的K值较实测值普遍小30%。 思考题： 苏制与美制土壤质地的转换方法？ 提示： 累计曲线，滤波 中国土壤可侵蚀性分布 3 作物管理因子 C 1) 作物管理因子 C影响因素 植被的抗侵蚀作用需要考虑垂直分层的作用，如：冠层、枯落层、根系等。 南方 vs 北方 冠层截留降水作用 不同森林类型的林冠截留率 枯枝落叶层的截留降水作用 根系的截留降水作用 林下灌木草本层的截流率 (六盘山) RUSLE/MUSLE 作物管理因子的C因子计算方法？ 1 采用了植被生长模型模拟不同阶段植被的生物量，冠层高度、盖度，植被残茬等。 2 遥感估算， NDVI 混合象元分解（植被、裸地、岩石、水面） 5 USLE/MUSLE量纲与数据流程 5.1 量纲与单位转换常数 公制与英美制度量衡 长度