层质量新通量垂直分布的数值模拟跟沙粒跃动特征参数的.pdfVIP

  风蚀沙流层质量通量垂直分布的数值模拟 及沙粒跃动特征参数的估算1   2 李万元 ，吕世华，沈志宝 （中科院寒区旱区环境与工程研究所西部气候环境与灾害重点实验室，甘肃 兰州， 730000） 摘 要 风蚀地表上跃动沙流层的垂直质量通量结构与沙丘形成直接相关，也与沙流危害程度密不可分。Shao （2005 ）建立的沙流层相似模型能够很好地模拟出裸露、干燥、疏松、单一沙地地表上沙流层的质量通 量垂直分布廓线，Liu et al. （2006 ）也针对该类地表作过一系列风洞试验，试验结果为相似模型的检验和 修正提供了实践依据。本文将基于试验数据和模拟结果的对比分析，结合模型参数的物理意义和前人的 相关定量表征方法，对沙粒跃离地表的平均初速度和地表粗糙度分别进行修正，从而改进模型对沙流层 质量通量垂直分布廓线的模拟能力；然后利用改进模型对地表沙流通量、沙流层的特征高度和平均高度 等重要参量进行模拟分析，以揭示其对摩擦速度、沙粒粒径和地表粗糙度的响应方式及其机制；最后依 据修正后的沙粒跃离地表的平均初速度合理估算出空气在垂直方向上对跃升沙粒的平均作用力，并详细 剖析摩擦速度、沙粒粒径和地表粗糙度对初速度和作用力的影响方式，这些工作对研究沙粒跃动的特征 轨迹、合理估算沙粒撞击地表的平均速度以及空气对沙粒跃动过程的平均作用皆具理论和实践意义。 关键词：跃动沙流层；质量通量垂直结构；相似模型；沙流特征和平均高度；沙粒跃离地表的初速度； 空气对跃升沙粒的垂直作用力 1 引言 风蚀地表附近形成的“沙尘云团”中，绝大部分组分粒子（约占 75 ％以上）是跃动着 的沙粒[1-3] 。跃动沙粒所致的危害并非仅仅通过总沙流通量的作用来实现[1,4-11] ，同时也取 决于沙流中质量通量的垂直分布状况。其实，后者正是无数单个沙粒跃动轨迹的总体表现， 越来越多的科学家已将之视为研究风吹沙流的关键环节[3,12-15] 。风蚀地表上沙流层的垂直                                                                  1 本文的英文版已发表于 Boundary-Layer Meteorology (2008) 127: 313-332 。 2 E-mail: ywl@lzb.ac.cn.   648   结构也与沙丘的形成直接相关。在野外观测或风洞试验中，总沙流通量一般不能直接获得， 只能从不同高度上测得的沙粒质量和风速廓线出发，先给出表征质量通量垂直分布廓线的 经验性函数，然后将该函数对高度积分，积分结果被视为总沙流通量[2,9,16] ；可见，在这一 实践意义上，沙流质量通量的垂直分布也理应受到足够的重视。 多年来，大量的风洞试验和野外观测使人们对风蚀沙流质量通量的垂直分布状况已有 了很深的认识。实践证明，风蚀过程中约 90%的粒子是集中于近地面 30-87 厘米的范围内 进行输送的，60-80%集中于 5-10 厘米的范围内，平均传输高度约为 63 厘米[2-3,14] 。科学界 已达成广泛共识的是，在地表的跃动沙流中，特别当风速处于中等或较大水平时，沙粒质 量通量随高度增加而呈指数递减，且递减率随风速和沙粒的增大而减小，沙流的平均高度 却随风速和沙粒的增大而增加[2-3,9,13,16,18] 。也有一些科学家发现，沙流通量或沙粒浓度会随 高度的增加以幂函数形式递减[2,14,18-19] 。Chen et al. （1996）却发现塔克拉玛干沙漠地表上 的沙流通量垂直分布廓线既不能用指数函数也不能用幂函数来描述[8,15,18] 。Dong et al.[3,14] 在研究风蚀地表纵向（顺风向）尺度对沙流通量廓线的影响时发现，在沙粒组成的地表上， 沙流通量先是随高度增加而增加，然