• CN 62-1175/P
  • ISSN 1006-7639
  • 双月刊
  • 中国科技核心期刊
  • 中国学术期刊综合评价数据库统计源期刊
  • 中文科技期刊数据库收录期刊

干旱气象, 2023, 41(4): 531-539 DOI: 10.11755/j.issn.1006-7639(2023)-04-0531

综述

对流起沙的研究进展

梁国豪,1, 毛睿,2,3, SHAO Yaping4, 李晓岚5, 龚道溢6

1.北京师范大学环境演变与自然灾害教育部重点实验室,北京 100875

2.北京师范大学国家安全与应急管理学院,北京 100875

3.中国科学院地球环境研究所黄土与第四纪地质国家重点实验室,陕西 西安 710061

4.德国科隆大学地球物理与气象研究所,科隆 50923

5.中国气象局沈阳大气环境研究所,辽宁 沈阳 110166

6.北京师范大学地表过程与资源生态国家重点实验室,北京 100875

Research progress of convective turbulent dust emission

LIANG Guohao,1, MAO Rui,2,3, SHAO Yaping4, LI Xiaolan5, GONG Daoyi6

1. Key Laboratory of Environmental Change and Natural Disaster of Ministry of Education, Beijing Normal University, Beijing 100875, China

2. School of National Safety And Emergency Management, Beijing Normal University, Beijing 100875, China

3. State Key Laboratory of Loess and Quaternary Geology, Institute of Earth Environment, Chinese Academy of Sciences, Xi’an 710061, China

4. Institute for Geophysics and Meteorology, University of Cologne, Cologne 50923, Germany

5. Shenyang Institute of Atmospheric Environment, China Meteorological Administration, Shenyang 110166, China

6. State Key Laboratory of Earth Surface Processes and Resource Ecology, Beijing Normal University, Beijing 100875, China

通讯作者: 毛睿(1980—),男,甘肃兰州人,教授,主要从事气候变化及其影响研究。E-mail:MR@bnu.edu.cn

责任编辑: 黄小燕;校对:蔡迪花

收稿日期: 2022-10-26   修回日期: 2022-12-30  

基金资助: 环境演变与自然灾害教育部重点实验室项目(2022-KF-04)
黄土与第四纪地质国家重点实验室开放基金项目(SKLLQG-2220)

Received: 2022-10-26   Revised: 2022-12-30  

作者简介 About authors

梁国豪(1996—),男,山东枣庄人,硕士生,主要从事大气沙尘气溶胶研究。E-mail:lgh@mail.bnu.edu.cn

摘要

对流起沙是一种由热力对流湍流直接夹卷沙尘进入大气的起沙机制,由于其频繁发生,长期累积的沙尘排放贡献不可忽视。根据国内外研究现状,总结归纳对流起沙的发生机理和影响因素,介绍目前相关的参数化方案,并从发生条件和起沙通量等方面分析对流起沙与跃移起沙、尘卷风的异同。最后,对对流起沙的野外观测和参数化方案改进等方面提出建议:应加强对流起沙潜在活跃区的观测及比较不同沙源地对流起沙的特征,进一步构建土壤湿度、植被覆盖等影响因素的修正函数以提升模式模拟性能。

关键词: 对流起沙; 跃移起沙; 尘卷风; 研究进展

Abstract

Convective turbulent dust emission (CTDE) is a new dust emission mechanism that thermal convective turbulent directly entrains dust particles into the atmosphere. Due to frequent occurrence, the long-term contribution of CTDE cannot be ignored. In the light of the previous research, the mechanism and influencing factors of CTDE are summarized, the parameterization schemes of CTDE are introduced, and the similarities and differences between CTDE and saltation-bombardment and/or aggregation-disintegration dust emission as well as dust devil are compared in terms of occurrence condition and dust emission flux. Finally this paper provides a reference for studying on CTDE and gives suggestions on field observation and improvement of parameterization for CTDE, including measuring CTDE events in potentially active areas, comparing the characteristics of CTDE between different dust sources, and building further correction functions of soil moisture and vegetation cover to improve the model performance.

Keywords: convective turbulent dust emission; saltation-bombardment and/or aggregation-disintegration dust emission; dust devil; research progress

PDF (4772KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

梁国豪, 毛睿, SHAO Yaping, 李晓岚, 龚道溢. 对流起沙的研究进展[J]. 干旱气象, 2023, 41(4): 531-539 DOI:10.11755/j.issn.1006-7639(2023)-04-0531

LIANG Guohao, MAO Rui, SHAO Yaping, LI Xiaolan, GONG Daoyi. Research progress of convective turbulent dust emission[J]. Arid Meteorology, 2023, 41(4): 531-539 DOI:10.11755/j.issn.1006-7639(2023)-04-0531

引言

沙尘气溶胶是大气气溶胶的重要组成部分,全球每年约有1 000~3 000 Tg的沙尘气溶胶释放到大气中(Zender et al.,2004),相当于对流层气溶胶总量的50%左右(Cakmur et al.,2006),对区域和全球尺度气候变化有重要影响。沙尘气溶胶通过吸收和散射太阳辐射直接影响地-气系统的辐射平衡(Sokolik and Toon,1996;Heinold et al.,2011),作为云凝结核或冰核影响云微物理过程,进而间接影响降水和云的变化,沉降到冰雪表面,影响冰雪的反照率和融化速度(Bullard et al.,2016)。此外,沙尘气溶胶携带铁、磷等微量营养素沉降到海洋,能大幅度提高浮游植物的初级生产力(Jickells et al.,2005;Mahowald et al.,2005)。因此研究大气沙尘气溶胶的来源和数量对气候变化研究有重要意义。

沙尘排放是大气沙尘气溶胶的来源,沙尘排放机制长期受国内外学者广泛关注。目前已经确认的沙尘排放机制有3种:沙粒轰击、土壤团聚体分裂和空气动力夹带(Shao,2001;Shao,2008)。其中,沙粒轰击和土壤团聚体分裂被认为是最有效和最重要的排放机制,大多数区域和全球沙尘排放模型中主要关注这两种机制(Marticorena and Bergametti,1995;Alfaro and Gomes,2001;Shao,2004;Shao et al.,2010),且二者均伴随跳跃过程,称为跃移起沙(Saltation-bombardment and/or Aggregation-disintegration Dust Emission, SADE)。空气动力夹带的起沙通量通常比跃移起沙小一个量级(Shao et al.,1993;Loosmore and Hunt,2000),在以往的沙尘模拟中该机制的贡献往往被忽略(Laurent et al.,2006;Park et al.,2010),然而空气动力夹带频繁发生,长期累积的排放量十分可观(Koch and Renno,2005;Shao,2008;Klose et al.,2014)。

在晴朗弱风天气条件下,太阳辐射导致的对流湍流对沙尘的直接夹卷是空气动力夹带的主要形式(Klose and Shao,2012;Li et al.,2014),称为对流起沙(Convective Turbulent Dust Emission,CTDE)。2016年内蒙古奈曼旗科尔沁沙地对流起沙通量约占该地区所有沙尘排放的27%(Ju et al.,2018),说明这一机制的长期贡献不可忽视。对流起沙在较长时间尺度上对全球沙尘排放有较大贡献(Chkhetiani et al.,2012),在影响大气辐射收支平衡及海洋初级生产力等方面具有重要意义(Klose,2014)。因此,本文基于近年来对流起沙的相关研究成果,围绕对流起沙机制及参数化方案综述研究进展,从发生条件和起沙通量等方面比较对流起沙和跃移起沙、尘卷风的异同,并对未来对流起沙的观测研究及参数化方案的改进提出展望。

1 对流起沙及其对区域大气环境影响

1.1 对流起沙机理

对流起沙是在晴朗弱风条件下大气对流湍流对沙尘颗粒的直接夹卷。如图1所示,由于平均风速较小,无法启动跃移起沙,而太阳辐射加热地表,浮力对大气湍流的发展起主要作用,形成对流混合层(Hess and Spillane,1990;Balme and Greeley,2006;Ito et al.,2010),旺盛的对流湍流不断将动量输送到地表,产生间歇性的强空气动力抬升,克服土壤粒子的重力和粘性力,直接抬升表面松散细小的沙尘颗粒进入大气(Klose,2014)。

图1

图1   对流起沙示意图(改自Klose and Shao, 2013

Fig.1   Schematic diagram of convective turbulent dust emission(adapted from Klose and Shao, 2013


不同于跃移起沙由强水平风将高空动量输送到地表,对流起沙的动量来源主要依赖湍流的大涡结构,因此湍流活动的强度与对流起沙的强度密切相关,同时由于湍流活动的不确定性,对流起沙的发生也具有很大的随机性。

1.2 对流起沙影响因素

1.2.1 地表要素

土壤颗粒的粒径及分布是影响对流起沙的决定性因素,反映土壤释放沙尘的能力。对流起沙释放的沙尘粒径大致小于20 μm,主要是位于土壤表层细小松散的土壤颗粒。然而,在某些热对流湍流活动很强且平均风较弱情况下却没有观测到对流起沙事件,这是由于地表可以释放的沙尘颗粒数量有限(Ju et al.,2018)。

植被覆盖是影响对流起沙的重要因素,主要体现在对起沙面积和感热通量的影响。植被覆盖的增加一方面会减少潜在起沙区域,同时减少到达地表的太阳辐射,加之植物的蒸腾作用导致潜热通量增加、感热通量减小,不利于对流起沙发生。科尔沁沙地观测站附近的植被覆盖率约0.5,仍然观测到对流起沙事件(Li et al,2014),相比之下跃移起沙在植被覆盖率大于0.2时,沙尘浓度出现的范围明显减少(张钛仁等,2008)。因此植被覆盖对对流起沙的影响需要结合更多观测资料做进一步研究。

土壤湿度增加使得土壤颗粒间粘性力增大,细小沙尘颗粒聚集成大团聚体,从而减少沙尘排放。Klose等(2014)对土壤湿度对对流起沙影响进行了初步校正,但改进效果不理想,原因主要是缺乏观测数据。土壤湿度对对流起沙影响仍需结合更多观测数据做进一步研究,包括不同土壤类型及湿度条件下沙尘数据。

1.2.2 大气湍流特征量

湍流强度受动力湍流和热力湍流共同影响,摩擦速度u*反映动力湍流强度,对流速度尺度w*反映热力湍流强度。热力湍流对对流起沙强度有主要贡献,如2004年春季浑善达克沙地非沙尘天气时起沙通量F为-5.0~5.0 µg·m-2·s-1,白天近地层热力湍流作用强于动力湍流,湍流交换以感热交换为主,正起沙通量可能与对流起沙有关(张宏升等,2007);科尔沁沙地发生对流起沙事件时F为0.0~10.0 µg·m-2·s-1,且伴有强感热通量,Fw*呈正相关(Li et al.,2014)。

感热通量与边界层高度反映了热力湍流的强度,地表感热通量越强,对流活动越旺盛,对流起沙通量越大(Ito et al.,2010;Ju et al.,2018)。同时,随着湍流活动增强,边界层高度增加,伴随着更高的对流起沙通量(Ju et al.,2018)。

1.3 对流起沙通量与大气湍流特征量的关系

利用WRF-Chem模式对2016年科尔沁沙地对流起沙进行数值模拟,起沙方案选择KS12方案(Klose and Shao, 2012),对科尔沁沙尘天气监测与观测实验站(Ju et al.,2018)附近模拟结果进行统计分析,得到对流起沙通量与大气湍流特征量的散点图(图2)。可以看出,对流起沙通量与对流速度尺度呈非线性正相关,即随着对流速度尺度增加起沙通量显著增加[图2(a)],决定系数R2为0.64,通过α=0.05的显著性检验,表明对流起沙主要受热力因素影响,其通量随地表感热通量增加近乎线性增加[图2(c)],决定系数R2为0.43,通过α=0.05的显著性检验,这与Ito等(2010)的研究结果相符,即较强的表面感热通量会导致更强的对流湍流,进而产生更大的起沙通量。随着湍流活动不断增强,边界层高度相应也增高,起沙通量与边界层高度也呈正相关[图2(d)],决定系数R2为0.51,通过α=0.05的显著性检验。对流起沙与摩擦速度之间相关性不明显[图2(b)],决定系数R2为0.17,说明动力因素对对流起沙贡献较小。以上结论与Ju等(2018)对科尔沁沙地的研究结果基本一致。

图2

图2   WRF-Chem模式模拟的2016年科尔沁沙地对流起沙通量与对流速度尺度(a)、摩擦速度(b)、感热通量(c)、边界层高度(d)散点图

Fig.2   Scatter plots between dust emission flux and convective velocity scale (a), friction velocity (b), sensible heat flux (c), planetary boundary layer height (d) simulated by WRF-Chem model at Horqin Sandy Land in 2016


1.4 对流起沙对区域大气沙尘含量影响

图3为2016年6月1—4日WRF-Chem模式模拟的沙尘柱质量浓度空间分布,可以看出科尔沁沙地对流起沙事件的传输过程。模拟结果显示,6月1日[图3(a)]在科尔沁沙地发生了对流起沙,此时沙尘范围较小;2日[图3(b)]沙尘的影响范围逐步扩大;3日[图3(c)]沙尘出现明显的向东传输直至到达日本上空;4日[图3(d)]出现一条明显的沙尘传输带,沙尘从日本本土继续向太平洋方向传输。总体来看,科尔沁地区的对流起沙事件对日本和太平洋地区沙尘具有一定贡献,传输的沙尘柱质量浓度约为1~10 mg·m-2。2016年科尔沁沙地全年对流起沙累积排放沙尘约1.588 Tg,对全球沙尘排放量的贡献约0.05%~0.16%,表明在沙尘数值模拟中应考虑对流起沙的影响。

图3

图3   WRF-Chem模式模拟的2016年6月1日(a)、2日(b)、3日(c)、4日(d)沙尘柱质量浓度空间分布(单位:mg·m-2

(绿色方框区域为科尔沁沙地)

Fig.3   The spatial distribution of mass concentration of dust column simulated by WRF-Chem model on 1 (a), 2 (b), 3 (c), 4 (d) June 2016 (Unit: mg·m-2

(The green box area represents the Horqin Sandy Land)


2 对流起沙与其他起沙机制比较

除对流起沙外还存在其他起沙机制,主要有跃移起沙和尘卷风。跃移起沙主要发生在扬沙、浮尘和沙尘暴等强沙尘事件中,当摩擦速度u*大于临界摩擦速度u*t时,大沙粒或土壤团聚体获得足够的动能后发射并冲击土壤表面,使土壤表面或土壤团聚体裂解产生沙尘颗粒进入大气。尘卷风是一种发生在边界层的垂直涡旋,其形成原因主要是太阳照射造成地表受热不均,在浮力作用下形成热力对流,在一定角动量作用下形成旋转上升的对流涡(段佳鹏等,2013),夹卷沙尘进入大气。

在起沙模型构建方面,Lu和Shao(1999)借鉴韧性金属侵蚀研究成果发展的跃移起沙体积模型已在沙尘模式中广泛应用(Shao,2001;Shao,2004;Shao et al.,2011);赵建华等(20052012)类比于电子跃迁模型,提出了冲击起沙跃迁模型,通过理论推导提出气动起沙和冲击起沙条件。Klose和Shao(2012)基于微物理过程构建了对流起沙模型,目前已应用于塔克拉玛干沙漠和澳大利亚沙漠。Rennó等(1998)基于热力学理论提出尘卷风总起沙量经验公式,并广泛用于尘卷风起沙量的估算;Klose和Shao(2016)认为尘卷风是对流起沙的一部分,可以将对流起沙模型用于尘卷风起沙量的估算。

表1列出对流起沙与其他起沙机制特征比较。可以看出,对流起沙以热力起沙为主(Klose and Shao,2013),发生条件同时满足摩擦速度小于临界摩擦速度(u*<u*t),对流速度尺度w*>0,大气稳定度(w*/u*)一般为3~6(Li et al.,2014)。对流起沙发生频率较高,主要发生在08:00—15:00(北京时,下同),单次对流起沙事件的持续时间一般小于3 h,在夏季最旺盛(Li et al.,2014)。尽管对流起沙事件的水平尺度和垂直尺度尚不清楚,但图3模拟结果表明,对流起沙在一定条件下也可以实现远距离传输。跃移起沙以动力起沙为主(Shao,2008),当摩擦速度大于临界摩擦速度(u*>u*t)时起沙发生,其发生频率较低(李耀辉等,2004),主要发生在春季(元天刚等,2016),影响空间尺度较大(吴焕波等,2018),在我国影响范围西起新疆,东到沿海地区(杨艳等,2012),垂直方向上沙尘可达对流层中上层(姜学恭等,2014),持续时间较长(段伯隆等,2021)。尘卷风与对流起沙相似,以热力起沙为主,发生条件同时满足摩擦速度小于临界摩擦速度(u*<u*t)、对流速度尺度w*>0、大气稳定度(w*/u*)>5(韩超信等,2021),表明尘卷风的发生相比对流起沙需要更强的热对流条件和更光滑的地表条件。尘卷风主要发生在12:00—18:00,且在14:00—15:00最旺盛,单次事件持续时间一般小于30 min,夏季发生较多(栾兆鹏等,2016),一般影响范围较小,垂直高度在几米到几千米不等(Sinclair,1969;Sinclair,1973)。

表1   对流起沙与其他起沙机制的特征比较

Tab.1  Comparison of characteristics between convective turbulent dust emission and other dust emission mechanism

特征对流起沙跃移起沙尘卷风
排放方式热力起沙为主(Klose and Shao,2013动力起沙为主(Shao,2008热力起沙为主(Klose and Shao,2016
发生条件同时满足u*<u*tw*>0、3<(w*/u*)<6
(Li et al.,2014
u*>u*t(Shao,2008同时满足u*<u*tw*>0、(w*/u*)>5
(韩超信等,2021
发生频率高(Li et al.,2014较低(李耀辉等,2004高(Balme and Greeley,2006
主要时段08:00—15:00(Li et al.,2014无明显规律12:00—18:00(栾兆鹏等,2016
单次事件时长大部分小于3 h(Li et al.,2014较长(段伯隆等,2021大部分小于30 min(栾兆鹏等,2016
季节特征夏季>秋季>春季>冬季(Li et al.,2014春季最多(元天刚等,2016夏季较多(栾兆鹏等,2016
水平尺度和
垂直尺度
尚不清楚水平范围较大,近地面到对流层中上层(姜学恭等,2014水平范围较小,几米至几千米不等
(Sinclair,1969;Sinclair,1973
起沙通量0.0~10.0 µg·m-2·s-1
(Li et al.,2014
27.7~206.0 µg·m-2·s-1
(沈志宝等,2003;沈建国等,2008
约0.7×106µg·m-2·s-1
(Han et al.,2016
对总沙尘排放的贡献率27.00%(Ju et al.,201865.17%~76.00%(Ju et al.,2018;刘莹等,201824.00%(刘莹等,2018

新窗口打开| 下载CSV


起沙通量方面,2002年4月敦煌戈壁沙漠的两次弱沙尘暴事件最大起沙通量为27.7 µg·m-2·s-1(沈志宝等,2003);2006年3月26日内蒙古朱日和地区的沙尘暴事件平均起沙通量为60.0 µg·m-2·s-1,最大值为206.0 µg·m-2·s-1(沈建国等,2008);相比之下科尔沁站观测的对流起沙通量约0.0~10.0 µg·m-2·s-1,偶尔达30.0 µg·m-2·s-1(Li et al.,2014),尘卷风的起沙通量为0.7×106 µg·m-2·s-1(Han et al.,2016),远大于对流起沙。可能原因有两点:一是尘卷风特殊的螺旋结构导致大量沙尘颗粒不断聚集和卷吸,沙尘浓度更高;二是计算对流起沙通量考虑的沙尘颗粒粒径小于20 µm,而尘卷风起沙通量的计算并没有粒径限制,大于20 µm的沙尘颗粒也包括在内。

贡献率方面,2006—2015年甘肃省民勤县强沙尘事件对总沙尘排放的贡献率平均为76.00%,其中尘卷风贡献率达24.00%(刘莹等,2018);2016年科尔沁沙地跃移起沙对总沙尘排放的贡献率为65.17%,而对流起沙的贡献率约27.00%(Ju et al.,2018)。综上所述,对流起沙的贡献率与尘卷风相近,而跃移起沙的起沙通量和贡献率都明显高于对流起沙。然而,近几十年我国北方主要沙尘源地强沙尘事件的强度整体呈减弱趋势,且春季减弱趋势最显著(元天刚等,2016;曹晓云等,2021;祖拜旦和依马木,2015;赵霞和张中伟,2011),在强沙尘事件减少的气候背景下对流起沙的长期贡献可能愈加重要。

3 对流起沙参数化方案

3.1 基于观测的经验参数化方案

研究表明,对流起沙与摩擦速度u*的相关性较低,主要受对流速度尺度w*影响(Li et al.,2014),因此Ju等(2018)利用2011—2015年科尔沁沙尘天气监测与观测实验站对流起沙事件的起沙通量F与对流速度尺度w*进行拟合,开发并验证了一种新的对流起沙经验参数化方案(F=2.93w*1.5),该方案对科尔沁地区对流起沙通量估算有很好的适用性。然而,该经验参数化方案对其他沙尘源地的对流起沙模拟可能并不适用,因为对流起沙不仅受对流速度尺度w*影响,土壤粒径、植被覆盖、土壤湿度等地表要素的影响也同样重要,这些影响因素需要在参数化方案中加以考虑。因此,要以对流起沙的物理过程为基础,开发适用性更广的参数化方案。

3.2 基于微物理过程的参数化方案

目前有关对流起沙的参数化方案主要为KS12方案(Klose and Shao,2012),该方案以地表中单个沙尘颗粒在对流湍流中的物理过程为基础,通过分析受力过程构建参数化方案,其特点是使用概率密度函数对空气抬升力和粘性力进行处理,体现了对流起沙过程的随机性(Klose and Shao,2012)。对于给定的空气抬升力Fl和颗粒间粘性力Fi沙尘排放通量F-为:

F-=αN2z-wtm+TFl-Fidz   Fl>Fi+mg,  z>d0              else

式中:z(m)为土壤表面的粘性亚层厚度;αN是与土壤质地有关的经验参数,量纲为m-2wt(m·s-1)为粒子最终速度;m(kg)为粒子质量;T(s)为粒子响应时间;g(m·s-2)重力加速度。对总沙尘排放通量F的计算公式如下:

F=0dmax00FlF-pjFidFipFldFlpAdjδdj

式中:p(F1)和pj (Fi)分别为空气抬升力F1和颗粒间粘性力Fi的概率密度函数;dj(μm)是沙尘粒子j的粒径,dmax为最大粒径,一般取20 μm; pAdj代表粒子的横截面积概率分布函数;δ表示方程的微分。Klose等(2014)基于大涡模拟提出了新的对流速度尺度wm=w*3+19u*31/3,用于F1的计算,p(F1)则采用韦伯分布。Klose和Shao(2012)根据Zimon(1982)的数据对Fi进行计算,p(Fi)近似为对数正态分布。该参数化方案对于F1Fi以及p(F1)pFi的参数化仍需进一步验证和改进,以提高模拟的准确性。

4 当前对流起沙研究重点

4.1 加强野外观测

目前,我国仅在科尔沁沙地开展了对流起沙观测研究,未来需要加强在其他区域(包括塔克拉玛干沙漠周边绿洲区、青藏高原雅江河谷区等)的对流起沙观测工作。塔克拉玛干沙漠周边地区土壤颗粒细腻疏松(Klose and Shao,2012),而青藏高原雅鲁藏布江及其支流河谷分布着大量的流动沙丘(韩永翔等,2004),拥有充足的沙源,同时上述地区中午地表感热通量高,有利于对流起沙发生。2021年7月初在雅江河谷地区多次目击到尘卷风,塔克拉玛干沙漠北缘的肖塘气象站观测显示1992—2011年3—9月尘卷风发生日数占全年总发生日数的90.9%(马明杰等,2019),因此有必要加强相关观测研究。

加强对流起沙野外观测主要包括3个方面:一是对沙尘源区土壤采样,完善不同类型土壤特征信息,包括土壤粒径分布、土壤湿度等;二是对近地面风速、气温等气象要素监测,计算相关大气湍流特征量;三是对沙尘浓度直接观测,进一步分析不同地区对流起沙发生频次、持续时间和空间变化。

4.2 下垫面影响因素校正

以往有关跃移起沙的参数化方案中,临界摩擦速度u*t直接参与沙尘排放通量的计算,首先对理想条件(干燥和裸露土壤)的u*t进行计算,然后考虑表面粗糙度、土壤湿度、含盐量、结皮等下垫面影响对其进行校正(Shao,2008)。由于Klose等(2012)构建的对流起沙方案中并没有使用u*t,需要在计算粘性力时考虑下垫面的影响,因此今后需要更多的实测数据来改进粘性力的修正,并进行敏感性试验。

4.3 完善参数化方案

目前Klose和Shao(2012)的对流起沙参数化方案展现了不错的模拟效果,但仍在某些方面需要进一步完善。首先方案中使用新的对流速度尺度wm=w*3+19u*31/3是大涡模拟结果,引入u*以反映动力湍流的贡献,然而在实际模拟过程中有可能夸大动力湍流的贡献(Ju et al.,2018),尤其在风速较大情况下应对u*加以限制,从而去除过高u*条件下的虚假起沙。另外方案中使用了经验参数αN,其代表土壤起沙能力的大小,直接影响起沙通量的计算,不同地区αN有一定差异,因此基于不同沙源地实测数据将αN构建为土壤类型、土壤湿度等影响要素的函数,将对全球对流起沙模拟提供一定帮助。

5 总结与展望

对流起沙是在白天大气热力对流旺盛且平均风较弱条件下,对流湍流直接将土壤表面松散细小的沙尘颗粒抬升进入大气的过程。相比跃移起沙,对流起沙发生更加频繁,累积排放时间和排放通量更加可观,在较长时间尺度上可能对全球沙尘气溶胶产生重要贡献。

Ju等(2018)通过观测资料拟合得出的经验参数化方案可以很好地模拟当地对流起沙,但对其他区域缺乏适用性。Klose和Shao(2012)以沙尘微物理过程为基础构建的对流起沙参数化方案考虑了空气抬升力和颗粒间粘性力的随机特性,对区域和全球的对流起沙模拟提供重要参考,但仍需进一步完善和发展。

基于当前研究现状,对未来对流起沙研究提出以下建议:(1)在对流起沙活跃地区加强外场观测,收集实时地面沙尘气溶胶观测资料,同时利用卫星观测优势,获取高分辨率的土壤湿度、植被覆盖度、地表粗糙度等资料(朱好和张宏升,2011);(2)研究不同沙源地对流起沙的异同,尤其比较荒漠化草原和沙漠戈壁的对流起沙物理过程,尝试探索其发生机理;(3)利用风洞试验和大涡模拟等试验结果进一步探究不同沙源沙粒粒径对对流起沙的影响,改进对流起沙参数化方案,完善植被覆盖和土壤湿度修正函数,提升模式模拟性能。

参考文献

曹晓云, 肖建设, 乔斌, , 2021.

1961—2019年柴达木盆地沙尘强度时空变化特征

[J]. 干旱气象, 39(1): 46-53.

[本文引用: 1]

段伯隆, 刘新伟, 郭润霞, , 2021.

“3·15”北方强沙尘暴天气成因分析

[J]. 干旱气象, 39(4): 541-553.

[本文引用: 2]

段佳鹏, 韩永翔, 赵天良, , 2013.

尘卷风对沙尘气溶胶的贡献及其与太阳辐射的关系

[J]. 中国环境科学, 33(1): 43-48.

[本文引用: 1]

韩超信, 汤耀国, 韩永翔, , 2021.

中国北方地区尘卷风时空分布的数值模拟

[J]. 干旱区地理, 44(4): 1 003-1 010.

[本文引用: 2]

韩永翔, 奚晓霞, 宋连春, , 2004.

青藏高原沙尘及其可能的气候意义

[J]. 中国沙漠, 24(5): 588-592.

[本文引用: 1]

利用1961-2000年高原91个站的气象实测资料并结合高原沙区分布状况, 分析了近40a青藏高原沙尘暴的时空分布特征。结果表明, 高原具备发生沙尘暴的前提条件, 大片的流动沙丘和大片荒漠化土地为沙尘暴的发生提供了充足的沙源, 年沙尘暴发生频率非常高, 大体有以羌塘高原为中心向东南逐渐减少的趋势。从12月至翌年4月, 沙尘暴发生中心从藏南的雅鲁藏布江上游河谷地区依次逐渐向北扩展到羌塘高原南部、羌塘高原及塔里木盆地南部, 这种季节性摆动与副热带西风急流的位置变化密切相关, 加上高原海拔4000~5000m的高度, 细粒物质被轻松地扬升到西风急流区, 传往遥远的北太平洋地区, 高原成为远程传输最高效的沙尘源地之一。沉降在北太平洋的沙尘, 加强了海洋生物泵的效率, 进而可能对全球气候产生影响。

姜学恭, 陈受钧, 云静波, 2014.

基于CALIPSO资料的沙尘暴过程沙尘垂直结构特征分析

[J]. 气象, 40(3): 269-279.

[本文引用: 2]

李耀辉, 张存杰, 高学杰, 2004.

西北地区大风日数的时空分布特征

[J]. 中国沙漠, 24(6): 715-723.

[本文引用: 2]

利用选取的西北5省(区)以及内蒙古西部(110&#176;E以西)分布均匀的127个气象站1960-2000年41a逐月大风出现日数资料, 分析研究了西北地区大风的空间、时间特征, 并具体分析了大风与沙尘暴的时空关系, 揭示了西北地区大风分布的一些新事实。西北大风天气可划分为较少区(年均大风日数小于10d)、较多区(年均大风日数10~50d)、多发区(年均大风日数50~100d)和频发区(年均大风日数大于100d)。西北地区大部分区域为大风较多区, 占总站数的614%, 大风频发区分布最小; 大风最频繁发生的地方在新疆西北部的阿拉山口, 年平均大风日数超过160d, 平均不到3d就有一次大风天气, 大风日数最少的地方是陕西北部延安, 平均每年发生大风天气的日数不到1d。大风日数空间分布与地形有很大关系, 两山之间的峡谷地带以及高山和青藏高原极易出现大风天气。西北地区多数台站近40a来大风呈减少趋势, 其中新疆西北部、甘肃河西走廊西部和陕西东部等地区减少最为明显, 大风增加的区域主要集中在新疆东北部到青海西部地区, 其代表站年均大风日数从20世纪60年代到80年代中期以后增加了近3倍, 达到190d。总体来讲, 西北地区大风天气最多的季节是春季, 以5月最多, 其次是夏季, 秋、冬季特别是秋季大风最少; 陕西、甘肃中南部夏季大风较多, 青海东南部则夏季最最少, 冬季大风更多一些。进一步分析表明, 西北地区大风频发区与沙尘暴频发区并不完全重合, 例如, 南疆是西北乃至我国沙尘暴最频发区之一, 但是南疆却是西北地区大风的较少区, 年大风日数远少于沙尘暴日数。同样是沙源丰富的沙漠地区, 也都是我国沙尘暴的主要频发区, 但是塔克拉玛干沙漠及其附近地区的沙尘暴日数远多于大风日数, 而巴丹吉林沙漠地区的大风日数却比沙尘暴日数明显偏多。最近40 a西北地区大风与沙尘暴发生次数随时间变化趋势一致, 基本呈线性减少特征, 这说明在下垫面状况不变或变化不大的情况下, 近年来沙尘暴次数减少可能主要是由于大风天气(沙尘暴驱动因子)减少而造成的。大风的时间变化可以决定沙尘暴随时间的变化。

刘莹, 韩永翔, 杨文清, , 2018.

沙尘天气, 尘卷风对沙漠地区起沙量的贡献

[J]. 中国沙漠, 38(6): 1 175-1 179.

[本文引用: 3]

栾兆鹏, 赵天良, 韩永翔, , 2016.

干旱半干旱地区尘卷风研究进展

[J]. 沙漠与绿洲气象, 10(2): 1-8.

[本文引用: 4]

马明杰, 杨兴华, 何清, , 2019.

塔克拉玛干荒漠-绿洲过渡带尘卷风活动特征——以肖塘为例

[J]. 中国沙漠, 39(2): 115-121.

DOI      [本文引用: 1]

利用1992-2011年塔克拉玛干沙漠北缘荒漠-绿洲过渡带肖塘气象站的观测资料,分析了该地区尘卷风的年、月变化规律及其与气象因子的关系。结果表明:(1)1992-2011年尘卷风发生日数总体呈波动递减趋势;尘卷风主要发生在3-9月,占全年总日数的90.9%,其中4-7月占全年总日数的70%左右。(2)尘卷风月发生日数随月平均地表与1.5 m高处温差的增大而线性增加(r=0.875,Pr=0.89,Pr=-0.869,P<0.01)。

沈建国, 孙照渤, 章秋英, , 2008.

干旱草原地区起沙通量的初步研究

[J]. 中国沙漠, 28(6): 1 045-1 049.

[本文引用: 2]

沈志宝, 申彦波, 杜明远, , 2003.

沙尘暴期间戈壁沙地起沙率的观测结果

[J]. 高原气象, 22(6):545-550.

[本文引用: 2]

利用中日合作"风送沙尘的形成、输送机制及其对气候与环境影响(ADEC)的研究"项目敦煌站的沙尘粒子数浓度观测资料,估算了2002年4月敦煌地区两次沙尘天气过程中戈壁沙地的起沙率和起沙量,并分析了地面起沙量(质量通量)的尺度分布。在两次强度较弱、持续时间较短的沙尘天气过程期间,戈壁沙地的平均起沙率为1.58×10<sup>-8</sup>和9.95×10<sup>-9</sup>kg·m<sup>-2</sup>·s<sup>-1</sup>,最大起沙率为2.77×10-8kg·m<sup>-2</sup>·s<sup>-1</sup>。地面沙尘释放量的73%以上为d>5.0μm的大粒子。最后讨论了起沙率与摩擦速度的关系。

吴焕波, 陈强, 谷新波, , 2018.

一次强沙尘天气过程及其对PM10时空分布的影响

[J]. 干旱气象, 36(1): 117-123.

DOI      [本文引用: 1]

2017年5月3&mdash;6日内蒙古出现了一次强沙尘天气,内蒙古各城市AQI值达到500。利用常规气象资料和MICAPS资料,结合HYSPLIT(hybrid single particle lagrangian integrated trajectory)模式分析本次沙尘天气成因。利用PM10小时监测值和AQI资料,从时空两方面评估沙尘天气对途经各城市PM10质量浓度分布的影响程度。结果表明:地面受冷高压和蒙古气旋影响,冷锋快速东移,高空槽和低涡东移发展,使高空强冷空气南下东移,配合有利的层结条件和外来沙尘的长距离输送,导致沙尘天气爆发;沙尘天气对PM10质量浓度的影响表现为前期主要影响西部和东部,中期对全区均有影响,后期主要对东部有影响。另外,中期为污染最严重时段。空间影响可划分为西部、中部、东部和通辽市4类,且从西到东、从北到南有一定的滞后性,滞后时间为1~10 h。

杨艳, 王杰, 田明中, , 2012.

中国沙尘暴分布规律及研究方法分析

[J]. 中国沙漠, 32(2): 465-472.

[本文引用: 1]

在分析中国沙尘暴研究现状的基础上,对其时空分布规律、沙尘源地、移动路径、研究方法和手段进行了系统的总结。结果表明,中国沙尘暴分布与北方的沙漠、沙地、戈壁等荒漠化地区分布一致,存在4大沙尘源区。干枯的湖泊、弃荒土地、沙砾草场可能成为东北地区的主要沙源,应引起注意;沙尘暴移动路径主要有3条,分别归纳了各路径的特征;根据沙尘源区示踪方法的特点和实例,指出了适用条件和应用范围,为今后沙尘暴的研究和防治提供依据和思路。

元天刚, 陈思宇, 康丽泰, , 2016.

1961—2010年中国北方沙尘源区沙尘强度时空分布特征及变化趋势

[J]. 干旱气象, 34(6): 927-935.

DOI      [本文引用: 3]

利用1961&mdash;2010年我国北方沙尘源区134个地面气象站沙尘暴、扬沙和浮尘发生频率逐月资料,结合定义的沙尘指数,系统地分析了中国北方7大主要沙尘源区沙尘强度的时空分布特征及变化趋势。结果表明:西部的塔克拉玛干沙漠沙尘强度最高,多年平均沙尘指数高达235,其次是中部的阿拉善高原和鄂尔多斯高原沙漠群(182),东北部的呼伦贝尔沙地沙尘强度最小,多年平均沙尘指数仅为23。总体来看,近50 a来中国沙尘源区的沙尘强度呈明显减小趋势,沙尘指数在1972、1987和2000年出现突变。其中,中部的阿拉善高原和鄂尔多斯高原沙漠群沙尘强度减小趋势最为显著(-6.3 a-1),其次是西部的塔克拉玛干沙漠(-5.9 a-1)。EOF分析结果表明,中国北方沙尘源区各地的沙尘强度整体变化一致,塔克拉玛干沙漠、阿拉善高原和鄂尔多斯高原沙漠群是沙尘强度变化中心;在此基础上,东西部地区的沙尘强度呈明显的反相变化;此外,中国北方沙尘源区的沙尘强度在1980年代后期发生显著变化,沙尘指数至今处于较低值。

张宏升, 朱好, 彭艳, , 2007.

沙尘天气过程沙地下垫面沙尘通量的获取与分析研究

[J]. 气象学报, 65(5): 744-752.

[本文引用: 1]

张钛仁, 宋振鑫, 王金艳, , 2008.

植被参数变化对沙尘起沙影响机理的数值模拟

[J]. 高原气象, 27(2): 392-400.

[本文引用: 1]

分析了叶面积指数和植被覆盖变化对沙尘暴起沙过程的影响机理。在一个集成的沙尘数值模拟预报系统中, 将模拟区域的叶面积指数分成了0, 0.5, 3.0和卫星反演的实测值4种, 结果表明, 叶面积指数对沙尘的数值模拟有较大的影响, 利用卫星反演的叶面积指数计算沙尘浓度时, 沙尘浓度的变化介于叶面积指数为0~0.5之间。将植被覆盖度数据分成0, 0.2和卫星反演的实测值3个类型。当植被覆盖度都为0时, 沙尘浓度最大; 当植被覆盖度为0.2时, 沙尘浓度出现的范围明显减少。

赵建华, 张强, 隆霄, 2012.

冲击起沙的概念模型

[J]. 中国沙漠, 32(2): 323-330.

[本文引用: 1]

冲击起沙是最重要的起沙方式。鉴于Lu和Shao提出的体积模型的不足,通过分析沙床受到冲击后应力的传递规律,类比于电子的跃迁模型,提出了冲击起沙的跃迁模型来解决冲击起沙问题。该模型的基本思想是,沙粒冲击沙床时,能量会沿冲击点向周围粒子扩散,当某一粒子获得的能量大于跃出所需要的能量时,该粒子即从床面跃出;多余的能量将转化为该粒子的跃出初动能。冲击起沙问题由此变成了一个研究跃出能量的大小和能量传输规律的问题。根据该模型,基于能量传输的吸收系数法和传输方程近似法,给出了两维和三维情形下冲击起沙释放的粒子数、速度分布和冲击坑形状的表达式。两维情形下显示,冲击坑的大小与释放的粒子数和粒径大小正比,而且伴随着冲击坑的形成,可能会发生崩塌现象。吸收系数法显示,两维冲击坑随着吸收系数的增大是先增大,后反而出现了收缩。传输方程近似法显示,两维冲击坑的大小正比于kt。

赵建华, 张强, 袁铁, , 2005.

沙粒启动机制的理论分析

[J]. 中国沙漠, 25(6): 853-862.

[本文引用: 1]

从理论上对沙粒的初始启动与垂直抬升作了理论分析。结果指出,沙粒脱离地面的位置顺序是沙平面上凸拐点处、斜面上外突较大的地方、极大值点、水平面、极小值点以及凹槽的拐点处;提出了沙丘的平化运动、缩移运动与持相运动;气动力起沙的条件是:小粒径沙粒、大风、涡度场(与水平风同符号)、不稳定(层结不稳定与切变不稳定)以及粘性流;冲击起沙的条件是:大粒径沙粒、小风、涡度场(宜与水平风反符号)、稳定大气以及湍流;最后提出了沙尘暴启动的交替循环机制。

赵霞, 张中伟, 2011.

1960—2008年新疆沙雅沙尘天气变化特征及原因分析

[J]. 沙漠与绿洲气象, 5(4): 20-25.

[本文引用: 1]

朱好, 张宏升, 2011.

沙尘释放通量外场观测和参数化研究进展

[J]. 北京大学学报(自然科学版), 47(4): 768-776.

[本文引用: 1]

祖拜旦, 依马木, 2015.

巴楚县1961—2010年沙尘天气的变化趋势分析

[J]. 沙漠与绿洲气象, 9(5): 36-40.

[本文引用: 1]

ALFARO S C, GOMES L, 2001.

Modeling mineral aerosol production by wind erosion: emission intensities and aerosol size distributions in source areas

[J]. Journal of Geophysical Research: Atmospheres, 106(D16): 18 075-18 084.

[本文引用: 1]

BALME M, GREELEY R, 2006.

Dust devils on Earth and Mars

[J]. Reviews of Geophysics, 44(3), RG3003. DOI: 10.1029/2005RG000188.

DOI      [本文引用: 2]

BULLARD J E, BADDOCK M, BRADWELL T, et al, 2016.

High-latitude dust in the Earth system

[J]. Reviews of Geophysics, 54(2): 447-485.

DOI      URL     [本文引用: 1]

CAKMUR R V, MILLER R L, PERLWITZ J, et al, 2006.

Constraining the magnitude of the global dust cycle by minimizing the difference between a model and observations

[J]. Journal of Geophysical Research: Atmospheres, 111(D6), D06207. DOI: 10.1029/2005JD005796.

DOI      [本文引用: 1]

CHKHETIANI O G, GLEDZER E B, ARTAMONOVA M S, et al, 2012.

Dust resuspension under weak wind conditions: Direct observations and model

[J]. Atmospheric Chemistry and Physics, 12(11): 5 147-5 162.

[本文引用: 1]

HAN Y, WANG K, LIU F, et al, 2016.

The contribution of dust devils and dusty plumes to the aerosol budget in western China

[J]. Atmospheric Environment, 126: 21-27.

DOI      URL     [本文引用: 2]

HEINOLD B, TEGEN I, BAUER S, et al, 2011.

Regional modelling of Saharan dust and biomass-burning smoke: Part 2: direct radiative forcing and atmospheric dynamic response

[J]. Tellus B: Chemical and Physical Meteorology, 63(4): 800-813.

DOI      URL     [本文引用: 1]

HESS G D, SPILLANE K T, 1990.

Characteristics of dust devils in Australia

[J]. Journal of Applied Meteorology and Climatology, 29(6): 498-507.

[本文引用: 1]

ITO J, TANAKA R, NIINO H, et al, 2010.

Large eddy simulation of dust devils in a diurnally-evolving convective mixed layer

[J]. Journal of the Meteorological Society of Japan, 88(1): 63-77.

[本文引用: 3]

JICKELLS T D, AN Z S, ANDERSEN K K, et al, 2005.

Global iron connections between desert dust, ocean biogeochemistry, and climate

[J]. Science, 308(5718): 67-71.

DOI      PMID      [本文引用: 1]

The environmental conditions of Earth, including the climate, are determined by physical, chemical, biological, and human interactions that transform and transport materials and energy. This is the "Earth system": a highly complex entity characterized by multiple nonlinear responses and thresholds, with linkages between disparate components. One important part of this system is the iron cycle, in which iron-containing soil dust is transported from land through the atmosphere to the oceans, affecting ocean biogeochemistry and hence having feedback effects on climate and dust production. Here we review the key components of this cycle, identifying critical uncertainties and priorities for future research.

JU T T, LI X L, ZHANG H S, et al, 2018.

Parameterization of dust flux emitted by convective turbulent dust emission (CTDE) over the Horqin Sandy Land area

[J]. Atmospheric Environment, 187: 62-69.

DOI      URL     [本文引用: 12]

KLOSE M, SHAO Y, 2012.

Stochastic parameterization of dust emission and application to convective atmospheric conditions

[J]. Atmospheric Chemistry and Physics, 12(16): 7 309-7 320.

[本文引用: 10]

KLOSE M, SHAO Y, 2013.

Large-eddy simulation of turbulent dust emission

[J]. Aeolian Research, 8: 49-58.

DOI      URL     [本文引用: 4]

KLOSE M, 2014. Convective turbulent dust emission: process, parameterization, and relevance in the Earth system[D]. Cologne: University of Cologne.

[本文引用: 5]

KLOSE M, SHAO Y, LI X, et al, 2014.

Further development of a parameterization for convective turbulent dust emission and evaluation based on field observations

[J]. Journal of Geophysical Research: Atmospheres, 2014, 119(17): 10 441-10 457.

KLOSE M, SHAO Y, 2016.

A numerical study on dust devils with implications to global dust budget estimates

[J]. Aeolian Research, 22: 47-58.

DOI      URL     [本文引用: 2]

KOCH J, RENNO N O, 2005.

The role of convective plumes and vortices on the global aerosol budget

[J]. Geophysical Research Letters, 32(18), L18806. DOI: 10.1029/2005GL023420.

DOI      [本文引用: 1]

LAURENT B, MARTICORENA B, BERGAMETTI G, et al, 2006.

Modeling mineral dust emissions from Chinese and Mongolian deserts

[J]. Global and Planetary Change, 52(1/4): 121-141.

DOI      URL     [本文引用: 1]

LI X L, KLOSE M, SHAO Y, et al, 2014.

Convective turbulent dust emission (CTDE) observed over Horqin Sandy Land area and validation of a CTDE scheme

[J]. Journal of Geophysical Research: Atmospheres, 119(16): 9 980-9 992.

DOI      URL     [本文引用: 13]

LOOSMORE G A, HUNT J R, 2000.

Dust resuspension without saltation

[J]. Journal of Geophysical Research: Atmospheres, 105(D16): 20 663-20 671.

DOI      URL     [本文引用: 1]

LU H, SHAO Y, 1999.

A new model for dust emission by saltation bombardment

[J]. Journal of Geophysical Research: Atmospheres, 104(D14): 16 827-16 842.

[本文引用: 1]

MAHOWALD N M, BAKER A R, BERGAMETTI G, et al, 2005.

Atmospheric global dust cycle and iron inputs to the ocean

[J]. Global Biogeochemical Cycles, 19(4), GB4025. DOI: 10.1029/2004GB002402.

DOI      [本文引用: 1]

MARTICORENA B, BERGAMETTI G, 1995.

Modeling the atmospheric dust cycle. 1. Design of a soil-derived dust emission scheme

[J]. Journal of Geophysical Research: Atmospheres, 100(D8): 16 415-16 430.

[本文引用: 1]

PARK S U, CHOE A, LEE E H, et al, 2010.

The Asian dust aerosol model 2 (ADAM2) with the use of normalized difference vegetation index (NDVI) obtained from the Spot4/vegetation data

[J]. Theoretical and Applied Climatology, 101: 191-208.

DOI      URL     [本文引用: 1]

RENNÓ N O, BURKETT M L, LARKIN M P, 1998.

A simple thermodynamical theory for dust devils

[J]. Journal of the Atmospheric Sciences, 55(21): 3 244-3 252.

DOI      URL     [本文引用: 1]

SHAO Y, RAUPACH M R, FINDLATER P A, 1993.

Effect of saltation bombardment on the entrainment of dust by wind

[J]. Journal of Geophysical Research: Atmospheres, 98(D7): 12 719-12 726.

[本文引用: 1]

SHAO Y, 2001.

A model for mineral dust emission

[J]. Journal of Geophysical Research: Atmospheres, 106(D17): 20 239-20 254.

DOI      URL     [本文引用: 2]

SHAO Y, 2004.

Simplification of a dust emission scheme and comparison with data

[J]. Journal of Geophysical Research: Atmospheres, 109(D10), D10202. DOI:10.1029/2003JD004372.

DOI      URL     [本文引用: 2]

SHAO Y P, 2008. Physics and modelling of wind erosion[M]. Dordrecht: Springer Netherlands.

[本文引用: 6]

SHAO Y, FINK A H, KLOSE M, 2010.

Numerical simulation of a continental-scale Saharan dust event

[J]. Journal of Geophysical Research: Atmospheres, 115(D13), D13205. DOI: 10.1029/2009JD012678.

DOI      URL     [本文引用: 1]

SHAO Y, ISHIZUKA M, MIKAMI M, et al, 2011.

Parameterization of size-resolved dust emission and validation with measurements

[J]. Journal of Geophysical Research: Atmospheres, 116(D8), D08203. DOI: 10.1029/2010JD014527.

DOI      [本文引用: 1]

SINCLAIR P C, 1969.

General characteristics of dust devils

[J]. Journal of Applied Meteorology and Climatology, 8(1): 32-45.

[本文引用: 2]

SINCLAIR P C, 1973.

The lower structure of dust devils

[J]. Journal of the Atmospheric Sciences, 30(8): 1 599-1 619.

DOI      URL     [本文引用: 2]

SOKOLIK I N, TOON O B, 1996.

Direct radiative forcing by anthropogenic airborne mineral aerosols

[J]. Nature, 381: 681-683.

DOI      [本文引用: 1]

ZENDER C S, MILLER R, TEGEN I, 2004.

Quantifying mineral dust mass budgets: terminology, constraints, and current estimates

[J]. Eos, Transactions, American Geophysical Union, 85(48): 509-512.

[本文引用: 1]

ZIMON A D, 1982. Adhesion of dust and powder[M]. New York: Consultants Bureau.

[本文引用: 1]

/