临界起沙摩阻风速u.docxVIP

临界起沙摩阻风速u 1 阵风锋的临界起沙阈值的确定 风蚀是土壤在风荷载下的空间位移过程，发生在世界干旱和半干旱地区。风蚀产生的气溶胶是目前大气中最主要的气溶胶来源,据估计全球每年粉尘释放量达3 000 Mt。风蚀的产生是2种力共同作用的结果,即:空气动力使得颗粒脱离地表,重力和颗粒间的粘性力阻碍颗粒的运动。前者由代表地面风剪切的摩阻风速u*来表征,后者由风蚀发生所需最小的摩阻风速——临界摩阻风速u*t来表示。u*和u*t之间的平衡受一系列的环境因子控制,如:风速、温度等天气条件,土壤类型,土壤状态,地表微尺度特征和地表大尺度特征等。因而,风蚀表现了较强的空间变化性和时间间歇性。 一个集成的沙尘天气数值预报系统包括气象模块、陆面模块、风蚀模块、输送沉降模块和GIS数据库。其中,风蚀方案是整个集成沙尘天气数值预报系统的核心,对风蚀机理的深刻认识是模拟沙尘天气的基础。风蚀模块利用从气象模块中获取的摩阻风速、陆面模块得到的土壤含水量,以及GIS数据库中获取的地表参数,对风蚀的临界起沙摩阻风速u*t,水平跃移通量Q和沙尘释放率F进行参数化。某一粒径组的水平跃移通量和沙尘释放率均为摩阻风速和临界摩阻风速(或临界风速)的函数。沙尘释放频率由风速超过临界起沙风速的次数决定,而沙尘释放强度则依赖于风速超过临界起沙风速的程度。临界起沙摩阻风速u*t或临界起沙风速ut是风蚀模块中最重要的参数。 早期,人们多利用风洞提供的可控环境研究和确定起沙的临界条件。但是,风洞不能完全准确地反映野外的真实条件。在自然条件下,由于风蚀过程的间歇性和准确获取大气和土壤要素的困难,起沙阈值的确定更为复杂。因此,有必要开展野外实验,采用不同的测量方法和手段,开展不同下垫面的风蚀起沙观测研究,确定临界起沙阈值。野外试验通常给出局地的临界起沙阈值,而沙尘天气数值模拟需要给出模拟区域内临界起沙阈值的时空分布。早期简单的沙尘模式一般采用统一的临界起沙阈值,如Westphal等在研究撒哈拉地区起沙通量时,临界起沙摩阻风速u*t统一采用0.6 m/s。集成的风蚀起沙模式中,通常基于风蚀起沙的物理模型,建立复杂的临界起沙摩阻风速参数化方案,将临界摩阻风速表达为土壤粒径分布、地表粗糙度、土壤湿度等因子的半经验函数。考虑影响临界起沙阈值的因子众多,为减小参数化方案中的不确定性,一些学者基于地面风速和沙尘天气出现频率,利用沙尘暴源区的常规地面气象观测资料,采用统计方法确定不同沙尘源区的临界起沙阈值。 确定临界起沙阈值的不同方案可以互相校验,其结果有利于更加深刻认识起沙机理、校准临界起沙摩阻风速(或临界风速)参数化方案。本文从野外观测、模式参数化方案和统计理论3个方面回顾了近年来临界起沙阈值(临界起沙风速ut或临界起沙摩阻风速u*t)的确定方法,探讨不同方案的优缺点和代表性。针对存在的问题和难点,对今后的研究和发展提出一些建议。 2 风蚀起沙临界速度计算 获取临界起沙风速的野外实验通常需要风速和沙尘两方面的观测。 采用置于地表或距地表几厘米高度处的压电传感器记录沙粒的跃移活动[1,14,15,16,22,23,24],通过分析沙尘颗粒跃移后撞击传感器次数与风速的关系,给出风蚀起沙判别的相对标准。将引起跃移活动的最低风速作为临界起沙风速。这种方法相对简单,但在实际应用中,由于风力和跃移活动的不完全相关,较难给出固定的临界风速代表值。 为了克服上述缺点,Stout等建立了“时间分数相等方法”(time fraction equivalence method)定量计算临界起沙风速。定义沙粒跃移的时间分数γp为给定时间段内,观测点探测到跃移沙粒的时间分数。假设沙粒跃移的时间分数γp与风速超过其临界值的时间分数γu相等,通过迭代方法计算起沙临界风速,具体地: 首先,根据每个采样间隔的沙粒撞击测量p(t)计算沙粒跃移的间歇性因子bp(t): bp(t)=0,若p(t)=0;bp(t)=1,若p(t)0 沙粒跃移的时间分数γp为采样时段内bp(t)观测的平均值: γp=1ΝΝ∑i=1bpi=ˉbp(t)(1)γp=1N∑i=1Nbpi=bp(t)ˉˉˉˉˉˉˉ(1) 同理,先给定临界风速的初始估计值ut,定义风速的间歇性因子bu(t): bu(t)=0,若u(t)ut;bu(t)=1,若u(t)ut 风速超过其给定临界值的时间分数γu为: γu=1ΝΝ∑i=1bui=ˉbu(t)(2)γu=1N∑i=1Nbui=bu(t)ˉˉˉˉˉˉˉ(2) 若γuγp,则增大ut,使得γu减小;若γuγp,则降低ut,使得γu增大。重复迭代,当(γu-γp)满足某一阈值时,得到ut值。这种方法确定的ut和γp受到测量跃移活动仪器灵敏度的影响。Stout分析了平均时间对计算临界风速的影响,指出平均时间接近零时,对应的临