1、 Ge -疏松土壤的净扬起量kg/（m2Gt-跃移颗粒被截留量kg/（m2s）。疏松土壤的净扬起量不仅受到地表覆盖的影响,同时也受到地表粗糙度的影响,其表达式如下7: （2） （3） （4）式中,Ce-扬起系数（1/m）;qe-模拟区域最大的扬起土壤量kg/（mCcov-地表残茬覆盖率（%）;Fcov-由团聚体、地表壳和石块覆盖而不会产生扬沙的地表所占百分比;Za-模拟区域空气动力学粗糙度（mm）;u*-模拟摩阻速度（m/s）;u*t-地表临界摩阻速度（m/s）。对于给定的地表状况,只要土壤释放量超过风的输送量,则跃移和蠕移颗粒运动量就会下降。另外,当挟带沙尘的风穿过粗糙地表时,大部分土壤颗粒
2、仍随风运动,但有部分的跃移和蠕移土壤颗粒被地表覆盖物（土块,石块和倒伏残茬等）和直立残茬截留。其截留量计算式为7,8: （5） （6） （7） （8）式中,Ci-直立残茬的截留系数（1/m）;n-单位面积的茎秆数;D-茎秆直径（mm）;Hres-茎秆高度（mm）;Ct-地表覆盖物造成的粗糙地表的截留系数（1/m）;qc-跃移和蠕移的输送能力kg/（m1.2控制部分摩阻速度的大小是决定土壤风蚀发生与否的重要因素。只有当摩阻速度超过临界摩阻速度时,地表才发生土壤风蚀现象。1.2.1摩阻速度摩阻速度主要受模拟区域空气动力学粗糙度的影响。风蚀观测过程中,由于条件有限,不可能在每一模拟区域内安装气象观测
3、站。由于气象观测点与模拟区域之间的环境存在差异,空气动力学粗糙度也就不同,因此在两环境下的摩阻速度也必将存在一定差异。为了便于将气象观测站的风速直接用于模拟区域风蚀研究,对气象观测点的摩阻速度进行了修正。研究表明9,在近地100 m高度内,风速随高度的变化符合对数关系,故气象观测点的摩阻速度可表示为7,10: （9）式中,u-风速观测高度Z处实测风速值（m/s）;Z0-气象观测点的空气动力学粗糙度（mm）。地表无直立残茬时,模拟区域摩阻速度为:u*=u*wZaZ00.067（10） （10）地表有直立残茬时,模拟区域摩阻速度为:（11）1.2.2临界摩阻速度临界摩阻速度的大小体现地表土壤可蚀性
4、程度,其主要受地表粗糙度、石块和地表壳覆盖、残茬覆盖、地表团聚体大小和密度以及地表土壤湿度的影响。在国内已有的风蚀模型中,考虑了地表壳、石块和2 mm的团聚体覆盖对风蚀的影响,没有考虑残茬覆盖,本模型全面考虑了这些因素对土壤风蚀的影响。裸露的地表,临界摩阻速度主要考虑到地表覆盖（壳、石块和2 mm的团聚体）和地表粗糙度的影响,计算式为11,12: （12）（1）覆盖的影响由于残茬覆盖可以进一步减少裸露的地表面积,对其下伏的土壤颗粒起到保护作用,从而增加了地表的抗风蚀能力。临界摩阻速度增加量与地表残茬覆盖率之间的关系为13: （13）式中,u*tR-由残茬覆盖引起的临界摩阻速度变化量（m/s）;
5、Sallaway14通过试验得出裸露地表所占百分比与地表残茬覆盖量之间呈负指数关系,本研究通过田间试验测定,证实了这一结论,并得出适合丰宁地区春麦作物的地表残茬覆盖率Ccov与残茬覆盖量Pres（t/hm2）关系如下: （14）由于河北坝上高原属于半干旱地区,降雨量少,故残茬的腐解过程对地表残茬覆盖量的影响不大,所以本模型中未考虑残茬覆盖的腐解过程,主要考虑耕作对地表残茬覆盖量的影响,借鉴PERFECT中的“残茬耕作子模型”,并根据丰宁坝上的试验测定,标定不同机具作业对残茬覆盖率的影响（如表1所示）,用于定量分析不同的耕作体系对土壤风蚀的影响。（2）土壤含水量的影响土壤湿度直接影响着土壤的抗风
6、蚀能力,土壤在湿润状态下凝聚力强,不易发生风蚀。Saleh和Fryrear15研究表明随着地表含水量的增加,临界摩阻速度增加,其增加量与土壤含水量之间呈线性关系。本研究针对试验区典型砂壤土,通过风洞试验建立了临界摩阻速度增加量与含水量之间的定量关系（如图1所示）: （15）式中,u*tW-由土壤含水量变化引起的临界摩阻速度的变化量（m/s）;H-地表（05 cm）土壤含水量。则经土壤含水量和残茬覆盖修正的地表临界摩阻速度u*t为: （16）表1不同耕作机具对地表残茬覆盖的影响Table 1Effect of tillage machinery on straw mulch机具Machinery
7、铧式犁翻耕Moldboard plough浅松Shallow tillage圆盘耙Disc harrow深松Subsoiling免耕No-tillage冬季休闲Winter fallow覆盖的减少率（%）Coverage decrease100*2630*20*35*为丰宁坝上测试结果;Experiment results in Fengning, Hebei.图1临界摩阻风速与土壤含水量的关系Fig.1Relationship between threshold wind velocity and soil water content1.3影响风蚀的因素1.3.1土壤含水量旱地农田的主要供水
8、项为降水,主要支出项为地表径流和蒸散,连续日土壤水分平衡方程如下所示 （17）式中,swi和swi-1分别为第i日和第i-1的土壤含水量（mm）;P-降水量（mm）;ET-蒸散量（mm）;Q-径流量（mm）;D-渗漏量（mm）。对于土壤风蚀研究来说,土壤-大气界面的土壤湿度是直接影响土壤可蚀性的重要因子。由于降雨大到可以产生径流时,不可能有风蚀出现,所以研究风蚀时可以不考虑地表径流。此外,河北坝上地区年降水量较少,属于半干旱地区,渗漏量约为零,因此降水和土壤蒸散决定着地表的干湿程度。针对保护性耕作的特点,本模型主要考虑地表残茬覆盖对地表蒸发的影响。（1）潜在蒸散量蒸腾蒸发（蒸散）计算一般包括潜
9、在蒸散和实际蒸散,实际蒸散又可以分为地表蒸发和作物蒸腾两部分。本文主要考虑一年一熟地区冬春季节严重的土壤风蚀时期,因此实际蒸散只包含地表蒸发,而且主要决定于大气的物理状况、土壤物理状况和水文特性。本文选择Van Bavel16修正的Penman公式计算潜在蒸散量,主要是因为该方程中不包含经验常数或函数,只要求输入每天的平均气温、最高气温、最低气温、相对湿度、太阳辐射量和风速等气象资料,这些数据可以从气象站获得。而且Jensen通过大量的试验研究证明,该方法最适用于干旱、半干旱地区潜在土壤蒸散量的估计17,表达式为: （18）式中,ETP-潜在蒸散量（mm/d）;L-水汽化潜热（MJ/kg）;-
10、饱和水汽压随温度变化的曲线斜率（kpa/）;es、ea-分别为高度z上的饱和与实际水汽压（kPa）;K-紊流扩散系数kg/（m2kPad）;H-地表辐射能量MJ/（m2d）。（2）土壤蒸发土壤潜在蒸发与潜在蒸散、土壤含水量及地表面的覆盖情况有关。土壤潜在蒸发是在农田潜在蒸散的基础上考虑作物覆盖与残茬覆盖的影响而得。本模型主要估算发生土壤风蚀高峰期冬春季的风蚀土壤流失量,该地区为一年一熟作物种植区,根据田间实际测得在春季地表秸秆覆盖可以减少地面温度12,从而减少了地表蒸发,因此对于保护性耕作,土壤潜在蒸发主要考虑地表残茬覆盖的影响。研究表明,秸秆覆盖可明显地降低土壤水分蒸发,且二者之间呈明显的负
11、相关关系:（19）式中,Pres-地表残茬覆盖量（kg/hm2）。当土壤中水分充足时,地表土壤水分蒸发通常以潜在的蒸发率进行。当土壤中水分不充足时,假设自下层进入最上模拟层的土壤水量超过潜在的土壤蒸发量,则认为土壤实际蒸发等于潜在的土壤蒸发率;如果进入最上模拟层的土壤水量小于潜在的土壤蒸发量,则根据最上一层土壤有效水估计潜在的土壤蒸发率,其表达式为:EA=EPAw0.7EA=EP（Aw/0.7）0Aw 0.7EA= 0Aw 0 （20）其中,EA-实际土壤蒸发（mm/h）;Aw-地表土壤中实际有效水分与最大有效水分之比;w-凋萎含水率（体积%）;f-田间持水率（体积%）;v-地表土壤含水率（体
12、积%）。（3）土壤水分再分配达西定律和质量守恒原理是土壤水分运动所遵循的基本规律,田间土壤水分运动,可近似为一维垂直向的流动,则基于达西定律估算土壤水再分配的流量方程为18-20:Qi= （hi- 1 -hi）KT/D （21）式中,Q-t时间单位面积通过的水流量（m）;h-土壤水头（m）;K-土壤导水率（m/s）;T-时间步长（s）;D-土壤水在邻近土层间流动的距离（m）。根据质量守恒原理,土壤水分平衡方程为:Wi=WIi+Qni （22）Hi=Wi/Ti （23）式中,Wi-第i层土壤湿度（mm）;WIi-第i层初始土壤湿度（mm）;Qni-流入第i层的实际水量（mm）;Hi-第i层土壤含
13、水量（m3/m3）;Ti-第i层层厚（mm）。1.3.2空气动力学粗糙度空气动力学粗糙度主要受垄、地表覆盖（包括土块、石块和倒伏残茬等覆盖）和直立残茬的影响。本试验地是以平作方式种植为主,所以本模型以地表覆盖和直立残茬两影响因素计算空气动力学粗糙度。（1）地表覆盖的影响地表覆盖对空气动力学粗糙度的影响体现在随机粗糙度的大小,随机粗糙度受许多因素的影响作用,比如耕作、降雨量、径流以及土壤质地等。研究表明19,降雨量和耕作是对随机粗糙度影响最大的两因素。本模型也主要考虑这两个因素,由地表覆盖物影响的空气动力学粗糙度为: （24）式中,SZrr-当前的随机粗糙度,是指初始随机粗糙度（SZrro）经累
14、积降雨（Pc）冲刷和覆盖影响的随机粗糙度值,表达式为: （25）A= 91.08 + 765.8Fsil式中,Fsil、Fsand-分别为表土层（05 cm）中粘土和砂土含量。研究表明不同耕作机具对随机粗糙度影响程度不同,表2给出不同耕作机具对随机粗糙度影响程度的经验值4,21,22。表2不同耕作机具对地表粗糙度的影响Table 2Effect of tillage machinery on ground roughness1.00.60.80.2（2）直立残茬的影响直立残茬对空气动力学粗糙度的影响体现在残茬的叶和茎对风的阻碍作用,由于直立残茬的存在将明显提高零风速的高度即增加空气动力学粗糙度
15、。本试验区秋季联合收割机收割小麦一般留茬较低,因此直立残茬对空气动力学粗糙度的影响,主要表现在无叶茎杆对空气动力学粗糙度的影响,其表达式为: （26）其中,A= 28.41 - 3.72 lnD;B=-3.052 + 0.6lnD;C=- 8.33 + 1.54 lnD;n、D、Hres如前所述。由于有直立残茬时,直立残茬对空气动力学粗糙度的影响远大于地面覆盖的影响,所以在本模型中:Za=WZrr无直立残茬时Za=WZv有直立残茬时 （27）2模型初步验证本文利用在河北省丰宁县鱼儿山镇保护性耕作土壤风蚀测试的结果,对所建立的旱地保护性耕作风蚀模型进行了检验。2.1模型运行的准备模型要求输入的数
16、据包括:气象数据、耕作管理措施、土壤参数、模拟区域地理位置和模拟初始土壤含水量。2.1.1气象数据气象数据包括日平均气温、日最高气温、日最低气温、日相对湿度、日太阳辐射、日平均风速、小时平均风速、小时风向和小时平均降雨。上述气象数据均由气象站获得。2.1.2耕作管理措施不同的耕作管理体系是由一套耕作管理措施构成的,例如耕作日期、所用机具及残茬处理等。这些措施影响了地表残茬覆盖量、地表粗糙度、地表土壤团聚体结构及土壤含水量等,因而可根据风蚀模拟结果评价不同耕作体系防止风蚀的效应。这是本模型的主要目的。2.1.3土壤参数土壤参数包括土壤水分常数、土壤机械物理组成、容重和土壤反射率等。（1）土壤水分
17、常数包括凋萎含水量、田间持水量和饱和含水量,田间实测数据如表3所示:表3土壤水分常数（%）Table 3Soil water constant处理方式Treatment凋萎含水率Withered water content饱和含水率Saturated water content田间持水率Field water capacity重量含水率Water content（w/w）体积含水率（v/v）传统耕作Traditional tillage2.353.0131.940.8317.8922.90免耕覆盖No-tillage and mulch2.273.0633.344.9618.0324.34（2
18、）土壤机械物理组成和容重田间实测数据如表4所示:表4土壤机械物理组成和容重Table 4Soil physical characters and bulk density颗粒含量百分数（粒径:mm）Percentage of particles （particle size: mm）容重Bulk density（g/cm3）10.250.250.050.050.010.010.0050.0050.0010.00123.0455.878.353.131.058.561.2321.1251.212.644.212.18.641.44（3）土壤裸地反射率C变化范围为0.10（干燥、有机质含量高、黑色
19、、粘重土壤）到0.30（砂地）,当土壤含水量增加时,C略微减小。所试土壤为细砂壤土,有机质含量低于2%,则C值为0.15。2.1.4模拟区域地理位置及地表状况模拟区域地理位置主要包括模拟区的海拔、经度、纬度及子午线等。试区海拔 1500 m、经度115、纬度41。子午线是基于地球每小时旋转15计算出来的,将地球分成不同时间带,每一时间包含15经度,东边的边界即为标准子午线。对于中心地带从西经90到105的时间带的子午线为90,本初子午线为0根据上述确定方法,本试区子午线为105地表状况主要包括随机粗糙度、地表残茬覆盖量等。目前随机粗糙度多采用计量钉的方法测定。作者基于计量钉的测量原理,采用两相
20、交直线先确定一个水平面,然后测量线上各标记处地表距已知水平面的垂直距离,最后计算出标准差,根据定义其值即为该试区的随机粗糙度。传统耕作和保护性耕作处理试区随机粗糙度实际测量结果分别为14.9 mm和28.5 mm。地表残茬覆盖量测定采用一个1 m1 m的方框,在供测试的地块,按照“S”形,选择10个点,将方框摆放在地表,拣出方框内的全部秸秆（不要土下面的根茬）,将秸秆烘干,称重,10点秸秆称重平均后即为单位面积秸秆覆盖量（t/hm2）。保护性耕作处理试区地表残茬覆盖量实际测量结果为0.61 t/hm2。2.1.5土壤初始含水量以2002年3月23日开始田间风蚀观测试验为起点进行模拟,取当时实测
21、的土壤含水量为土壤水分模拟的初始值。2.2风蚀模型验证根据2002年3月23日至4月30日期间田间风蚀观测及风蚀模拟结果,验证风蚀模型的模拟效果。不同处理土壤流失量逐日预测值与实测值的比较如图2所示,统计关系如图3所示。从图中可以看出保护性耕作与传统耕作的土壤流失量预测值与实测值吻合的比较好,相关系数为0.94,斜率在0.85左右、预测值与实测值比值在1.15左右,均接近于1。3结语针对农田土壤风蚀,建立了保护性耕作土壤风蚀模型,以小时为步长,根据气象数据、地表土壤水分、秸秆残茬覆盖率及地表粗糙度,模拟不同耕作体系下农田土壤风蚀流失量情况,为评价不同保护性耕作措施对土壤风蚀的防治效应提供了依据
22、。国外的农田风蚀预报模型一般考虑秸秆残茬覆盖率及地表粗糙度两个因素,本模型在此基础上增加了可明显增加地表抗蚀性的因素地表土壤水分对风蚀的影响。从模拟效果来看,本文针对保护性耕作技术特点所设计的风蚀模型是成功的。由于只有1 a的试验数据,对模型的验证与评价还只是初步的,尚需要进一步验证。图2土壤流失量预测值与实测值比较Fig.2Contrast between simulated and observed value of soil loss图3土壤流失量模拟值与实测值的关系Fig.3Statistical relationship between simulated and observed
23、value of soill loss参考文献:1L E Wagner. An overview of the wind erosion prediction systemR. Kansas City, Missouri: Westin Crown Center,1996.2J R Williams, C A Jones ,P T Dyke. A modeling approach todeter mining the relationship between erosion and soilproductivity J. Transactions of the ASAE, 1984,27（1
24、）:129-144.3Littleboy M, Silburn D M, Freebairn D M .PERFECT-acom-pute simulation model of productivity erosion runoff functionsto evaluate conservation techniquesR. Brisbane, Australia:Queensland Department of Primary Industries, 1989.4王晓燕.旱地机械化保护性耕作径流与土壤水分平衡模型试验研究D.北京:中国农业大学工程学院, 2000.5L J Hagen,L
25、E Wager ,E L Skidmore. Analytical solutionsand sensitivity analyses for sediment transport in WEPS J.Transactions of the ASAE, 1999,42（6）: 1715-1721.6Stout J. Wind erosion within a simple fieldJ. Trans Am SocAgric Engin, 1990, 33（5）: 1597-1600.7Hagen L J. Crop residue effects on aerodynamic processes andwind erosion J. Theoretical and Applied Climatology,1996, 54: 39-46. 8Hagen L J, D V Armbrust. Aerodynamic roughness andsaltation trapping efficiency of tillage ridges J. Trans AmSoc Agric Engin. 1992, 35（4）: 1179-1184.9翁笃鸣,陈万隆,沈党成,等.小气候和农田小气候M.北京:农业出版社,1981.10朱朝云,丁国栋,杨明远.风沙物理学M.北
