地表土壤风蚀起沙的临界摩擦速度及其影响因素.docxVIP

地表土壤风蚀起沙的临界摩擦速度及其影响因素 1 对地面起沙量的影响 地表土壤的风蚀和沙化作用是一个动态过程，它反映了运动的空气流和土壤颗粒之间的相互作用。研究表明, 这是一个由许多因素控制影响的复杂的物理过程。对于风蚀起沙的诸多影响因子, 国内外已有很多研究。Wodruff等把地表土壤的风蚀起沙量用土壤侵蚀度、 地表粗糙度、 气候因子、 盛行风方向的田地长度 (或宽度) 和植被覆盖等五个因子表示, 建立了半经验性的土壤风蚀方程; Gillette等认为, 地表的沙尘排放量是风速和地表状况 (包括土壤水分和地表植被等) 的函数; Shao等把所有的影响因素概括地分为天气和气候条件 (主要是大风和少量的降水) 、 土壤状态 (包括土壤的矿物成分、 粒子尺度分布特征、 土壤硬度和土壤水分等) 、 地表粗糙元素 (非侵蚀性土壤集合、 植被覆盖度) 以及土壤的利用和管理等四个方面; 贺大良等利用室内风沙环境风洞对影响土壤风蚀起沙的风况 (含沙量和风速大小) 、 土壤表面的覆盖状况 (植被等) 、 地表物质组成和人为因素 (开垦和放牧) 等影响因子进行了初步的模拟试验。已有的许多风蚀起沙模型表明, 地面风蚀起沙量的多少在计算中主要受摩擦速度u(或风速u) 和地表风蚀起沙的临界摩擦速度u*t(或临界风速ut) 所影响, 上述诸多因素正是通过影响u*(或u) 和u*t(或ut) 进而影响地面起沙量的, 其中u*t主要受大风和地表粗糙度等的影响,u*主要受粒子尺度、 土壤水分、 植被覆盖以及土壤硬度等的影响。另外, Shao认为, 地表土壤的粒子尺度分布在风蚀起沙过程中随时间发生变化, 这一点对地面起沙量也有较大影响。 本文根据敦煌地区戈壁和绿洲的观测资料, 对地表土壤风蚀起沙的临界摩擦速度u*t及其在戈壁和绿洲变化特征的异同和风蚀起沙过程中地表土壤的粒子尺度分布p(d) 及其对垂直尘粒通量 (风蚀起沙量) 的影响进行了分析研究。 2 土壤硬度和地表土壤水分含量的测定方法 理论和试验研究表明, 在地表土壤的风蚀起沙过程中, 沙粒的跳跃运动是一个最关键的过程。对跳跃运动起支配作用的是空气动力, 而空气动力则与u2*成正比, 因此, 地表土壤的风蚀起沙与否必将对应于一个临界摩擦速度u*t, 当u*u*t时, 空气动力大于沙粒的重力和内部粘性力的合力, 跳跃运动不能发生, 地面会因此而起沙; 当u*u*t时, 空气动力小于微粒的重力和内部粘性力的合力, 跳跃运动不能发生, 地面也不会因此而起沙。 图1是根据Gillette等的起沙模型得到的垂直尘粒通量F随u*t的变化, 其中摩擦速度取定值1.0 m·s-1。由图1可见, 对一定的u*,F随u*t的增大而线性减小,u*t越小, 地面起沙量越大,u*t越大, 地面起沙量越小。 临界摩擦速度表征的是地表对风蚀起沙的阻碍能力, 主要受粒子尺度d、 土壤水分含量w、 植被覆盖度λ以及土壤硬度等的影响。其中粒子尺度d对粒子所受的重力、 空气动力和内部粘性力均有影响; 土壤水分含量w取决于降水、 蒸发以及土壤的持水性等; 植被覆盖度λ一方面影响地表粗糙度进而影响u*, 另一方面可使植被根部附近的土壤粒子结合的更紧, 不易因风蚀而脱离地表; 土壤硬度在一般情况下随时间变化非常缓慢, 但当土壤被耕作时, 其硬度将发生很大变化, 对此, 目前还没有一个较好的确定方法。Shao等在考虑d,w和λ等影响因素的基础上, 提出了一个计算u*t的方案, 即: u*t=Η(w)R(λ)√AΝ(σpgd+ερd),(1)u*t=H(w)R(λ)AN(σpgd+ερd)?????????????√,(1) 其中AN和ε是经验常数, 分别近似取0.0123和3×10-4kg·s-2;σP是土壤微粒密度 (2650 kg·m-3) 和空气密度 (1.23 kg·m-3) 的比值;g是重力加速度, 近似取9.81 m·s-2;ρ为空气密度。 H(w) 表征地表土壤水分对风蚀起沙的阻碍作用, 其表达式根据Fecan等的研究得到: Η(w)={1√1+a(w-w′)bw≤w′ww′,(2) 其中a,b为经验常数,w′是土壤水分是否具有明显阻碍作用的临界值,a,b和w′的大小均决定于土壤类型。 R(λ) 表征植被覆盖对风蚀起沙的阻碍作用, 其表达式根据Raupach的研究得到, 即: R(λ)={1√(1-mσλ)(1+mβλ)λ=0λ0,(3) 其中σ是植被的根部面积与叶面面积之比, 经验值取1.45;β是单个植被元素的拖曳系数与没有植被的地表的拖曳系数之比, 经验值取202;m是一个1的常数, 经验值取0.16;λ是植被的切面积指数, 决定于植被覆盖分数f, 由以下经验公式确定: λ={0-0.35ln(1-f)f=1f1.(4) 图2是依