气象学报  2018, Vol. 76 Issue (6): 968-982   PDF    
http://dx.doi.org/10.11676/qxxb2018.051
中国气象学会主办。
0

文章信息

岳治国, 余兴, 刘贵华, 戴进, 朱延年, 徐小红, 惠英, 陈闯. 2018.
YUE Zhiguo, YU Xing, LIU Guihua, DAI Jin, ZHU Yannian, XU Xiaohong, HUI Ying, CHEN Chuang. 2018.
NPP/ⅦRS卫星反演青藏高原夏季对流云微物理特征
NPP/ⅦRS satellite retrieval of summer convective cloud microphysical properties over the Tibetan Plateau
气象学报, 76(6): 968-982.
Acta Meteorologica Sinica, 76(6): 968-982.
http://dx.doi.org/10.11676/qxxb2018.051

文章历史

2018-03-13 收稿
2018-07-22 改回
NPP/ⅦRS卫星反演青藏高原夏季对流云微物理特征
岳治国1,2, 余兴2, 刘贵华2, 戴进2, 朱延年2, 徐小红2, 惠英2, 陈闯2     
1. 陕西省人工影响天气办公室, 西安, 710014;
2. 陕西省气象科学研究所, 西安, 710016
摘要: 青藏高原(下称高原)对东亚大气环流、气候变化及下游灾害性天气形成、发展有重要影响,研究青藏高原云微物理特征有重要意义。但因高原台站稀少,对云微物理研究不充分。NPP(National Polar-orbiting Partnership)卫星ⅦRS(Visible Infrared Imaging Radiometer Suite)传感器包含17个中分辨率通道(750 m)和5个高分辨通道(375 m),具有反演初生小块对流云的优势,能够利用NPP/ⅦRS反演对流云的微物理特征。利用NPP/ⅦRS卫星格点对流云云物理自动反演(Automatic Mapping of Convective Clouds,AMCC)软件对高原地区2013-2017年夏季(6-8月)过境的ⅦRS资料进行了反演,得到了高原对流云的宏、微观物理特征,并计算了这些物理量在0.33°×0.33°格点上的平均值。分析得出如下结论:(1)反演云底温度(Tb)与那曲探空计算抬升凝结温度(TLCL)线性相关,相关系数为0.87,均方根误差为3.0℃。(2)高原对流云宏、微观物理特征为:一是云底冷(Tb为-5℃),云底离地高度为1800-2200 m,云内含水量低;二是云底云凝结核数浓度(NCCN)为200-400个/mg,最大过饱和度(Smax)为0.7%,NCCN少,Smax大,云滴凝结增长速率更快;三是降水启动厚度(D14)小,为1500-2000 m,雅鲁藏布江流域及藏南地区D14约500-1000 m,更加容易形成降水;四是云顶海拔高度为10-13 km,云厚度从南部5000 m逐渐减小到北部2500 m,云厚有限;五是晶化温度高,从中部、南部-30℃到北部-25℃,加之高原Tb < 0℃,使得云内降水粒子以冰相为主。(3)高原对流云的这些微物理特征决定了其降水具有多发、短时、量小、滴大的特点。这些结论进一步深化了对高原夏季对流云的科学认识。
关键词: 青藏高原     ⅦRS     云微物理反演     对流云     云底温度     云凝结核    
NPP/ⅦRS satellite retrieval of summer convective cloud microphysical properties over the Tibetan Plateau
YUE Zhiguo1,2, YU Xing2, LIU Guihua2, DAI Jin2, ZHU Yannian2, XU Xiaohong2, HUI Ying2, CHEN Chuang2     
1. Office of Weather Modification of Shaanxi Province, Xi'an 710014, China;
2. Meteorological Institute of Shaanxi Province, Xi'an, 710016, China
Abstract: The Tibetan Plateau (TP) has important influences on the formation and development of the East Asian atmospheric circulation, climate change and disastrous weathers in its downstream regions. It is very important to understand the microphysical characteristics of clouds at the TP. However, observations and researches of the microphysical characteristics at this area are insufficient because of the sparse observation stations and few cloud physical data. ⅦRS (Visible Infrared Imaging Radiometer Suite) onboard the Suomi NPP (National Polar-orbiting Partnership) satellite, has an improved imaging spectroradiometer with seventeen channels of 750 m moderate resolution and five channels of 375 m image resolution. The high-resolution satellite has an advantage for observing the small or initial convective clouds. Based on the retrieval methodologies of cloud microphysical properties by NPP satellite were proposed, an automated mapping software package named AMCC (Automatic Mapping of Convective Clouds) has been developed to retrieve properties of convective cloud at the scale of satellite swath. The properties of convective clouds are automatically retrieved by AMCC and their values are averaged in 0.33°×0.33° grids via NPP/ⅦRS satellite data over the TP in the summers from 2013 to 2017. The results show that:(1) the temperature of lifting condensation level (TLCL) in Naqu meteorological station and the temperature of cloud base (Tb) from ⅦRS are linearly correlated with the correlation coefficient of 0.87 and the root-mean-square (RMS) error of 3.0℃; (2) convective clouds at the TP have the following microphysical properties. First, the cloud base temperature (Tb) is about -5℃, and the height of the cloud base above the ground (Hb) ranges between 1800-2200 m, and the cloud water content is small. Second, the cloud condensation nuclei concentrations (NCCN) is between 200 and 400 mg-1 with 0.7% of maximum supersaturation (Smax); consequently, the condensation growth rate of water cloud droplet with less NCCN and higher Smax is fast. Third, because the precipitation initiation depth (D14) varies within 1500-2000 m and 500-1000 m at the Yarlung Zangbo River Basin and the Zangnan region, respectively, and the clouds at these areas are more prone to precipitation. Fourth, the height of the cloud top above sea level (Htop) is between 10-13 km, but the cloud depth (Dcld) is rather small, which is about 5000 m at its southern part and gradually reduces to 2500 m at its northern part. Fifth, the glaciation temperature (Tg) ranges -30℃ at its central and southern parts to -25℃ at its northern part, which, combined with the high Tg and the Tb that is less than 0℃, lead to the domination of ice process in the clouds; (3) these microphysical properties of convective cloud at the TP explains why the precipitation is frequent and lasts over a short time with small amount and large rain drops.
Key words: Tibetan Plateau     ⅦRS     Retrieval of cloud microphysical properties     Convective cloud     Cloud base temperature     Cloud condensation nuclei    
1 引言

青藏高原(下称高原)是世界上最高的高原, 平均海拔在4000 m以上,对东亚大气环流、气候变化及灾害性天气的形成和发展都有重要影响(Wan, et al, 2017; 夏昕等, 2016; 徐祥德等, 2015)。由于其特殊的热力和动力作用,高原夏季对流和雷暴活动频繁(张鸿发等, 2003)。云的微物理特征(例如云滴大小、云滴浓度、云形成过程等)是人们认识云的重要途径。同时,云微物理特征的描述对天气和气候模式的预测也有重要影响。因此,研究青藏高原地区云微物理特征有着十分重要的意义。

目前,通常对云的观测有遥测和飞机穿云直接探测。地面遥测方式有激光云高仪、毫米波云雷达、激光雷达、红外成像仪等。其中激光云高仪和激光雷达可以准确地获得高分辨的云底离地高度(Hb)和云顶高度,毫米波雷达的波长更接近云粒子尺度,可以探测到云内结构(吴翀等, 2017; 张倩等, 2016)。卫星探测主要有星载云雷达(CloudSat)和激光雷达(Calipso),星载云雷达能测量到卫星轨道下云和气溶胶的垂直分布,但探测范围十分有限(刘建军等, 2017; 汪会等, 2011; 王宝鉴等, 2017; 张晓等, 2015)。对于云凝结核(cloud condensation nuclei, CCN)数浓度(NCCN)、云滴大小、云滴浓度等特征量,通常需要搭载探测仪器的飞机穿云探测(王黎俊等, 2013)。这些机载云物理探测设备和地基设备(常祎等, 2016; 刘黎平等, 2015)价格高、维护和标定过程复杂,在观测台站稀疏的高原地区使用较少。云底高度的观测主要有人工目测法、探空气球法、云高仪观测等。目测法主观性强、观测不连续,数据准确性差;探空气球法对人工依赖较强,费用高;云高仪观测精确,但未普及使用。

气象卫星观测资料具有全球连续覆盖、多时次、获取及时、便利的优点,广泛地应用于各种云微物理特征的反演。最近20年来,卫星的可见光和红外反演技术已经广泛用来反演云的粒子有效半径(re)、液态水路径、光学厚度、云顶温度(Ttop)和相态等参数。

Rosenfeld等(1998)利用NOAA/AVHRR资料反演了对流云温度(T)和云粒子有效半径(re)的分布特征(T-re分析),提出了研究对流云发展和云降水物理过程的方法(戴进等, 2010; 刘贵华等, 2011b)。这一反演方法也应用于MODIS(Rosenfeld, et al, 2011)、FY3/VIRR(刘贵华等, 2011a)等卫星的云微物理特征反演和强对流识别(Rosenfeld, et al, 2008)研究。

2011年10月28日NPP (National Polar-orbiting Partnership)卫星的成功发射,为卫星云微物理反演提供了新机遇。与MODIS(或AVHRR)等传感器相比,NPP卫星VIIRS (Visible Infrared Imaging Radiometer Suite)传感器热红外通道分辨率(375 m)提高了3倍,在精确反演初始及发展阶段的小对流云方面显示出独特的优势(Hillger, et al, 2013; Rosenfeld, et al, 2014b)。VIIRS的5个375 m高分辨率通道中心波长分别为0.64、0.865、1.615、3.745和11.45 μm,缺少独立的12 μm高分辨率通道,无法使用高分辨率分裂窗(即:11 μm与12 μm高分辨率通道的亮温差)计算水汽吸收用于3.745 μm通道辐射亮温反演云粒子有效半径(re)。Rosenfeld等(2014b)发展了针对VIIRS高分辨率资料的云微物理反演方法,使得反演对流云云底温度(Zhu, et al, 2014)和云底上升气流(Zheng, et al, 2015b)成为可能。这是对以前1 km分辨率云微物理反演能力的新突破,也使反演对流云云底NCCN成为可能(Rosenfeld, et al, 2016)。VIIRS资料也应用于气溶胶在对流云中成冰作用的研究(朱延年等, 2015)。

Zhu等(2014)利用VIIRS资料的T-re分析,发展了反演对流云云底温度(Tb)和边界层比湿的新方法。与位于美国南部大平原地面站点2 a云高仪和探空资料比较发现,反演的Tb标准差仅为1.1℃;通过Tb计算的边界层比湿的误差约为10%。以VIIRS资料反演的Tb、云底离地高度(Hb)及地面气温为基础,Zheng等(2015b)发展了一个反演对流云底上升速度(Wb)的方法,地面激光雷达检验反演的Wb标准差为0.41 m/s。Williams等(2002)研究了WbHb成正比的物理意义。Zheng等(2015a)进一步发展了计算Wb的方法,给出了HbWb的经验线性关系。地面和船载设备观测表明,此计算方法普遍适用于边界层对流云,计算的Wb平均绝对误差为21%。基于高分辨率VIIRS资料反演的对流云T-re廓线,Rosenfeld等(2014a)发展了反演对流云底最大过饱和度(Smax)及CCN活化成绝热云滴的方法。通过地基和船载CCN计数器验证,VIIRS卫星反演的NCCN误差为20%—30%(Rosenfeld, et al, 2016)。另外一些外场观测研究也发现,对流云内偏离绝热变化的现象广泛存在(Brenguier, et al, 2011; Min, et al, 2012),这会造成NCCN和绝热云滴数浓度(Nd, a)反演结果的不确定性(Merk, et al, 2016; Miller, et al, 2016)。

基于以上这些VIIRS高分辨率卫星资料在云微物理反演方法上的持续研究,利用自主研发的卫星格点对流云云物理自动反演软件,反演分析2013—2017年6—8月过境高原的VIIRS高分辨率卫星资料,对云物理特征进行格点平均,以获得近5年夏季高原对流云云物理特征的总体状况,部分弥补高原地区缺少云微物理仪器观测的不足,深化对高原区域云物理特征的科学认识。

2 资料和反演方法 2.1 资料

反演使用07时(世界时,下同)前后扫描高原区域的5 min一景的VIIRS资料,由NOAA数据服务网站(CLASS)提供。为了获得高原夏季平均云物理特征,本反演使用2013—2017年6—8月的共计569景VIIRS L1B资料。反演中还用到美国环境预报中心(NCEP)6 h间隔1°×1° FNL(Final Operational Global Analysis)资料的温度、气压和位势高度廓线作为辅助资料。

https://www.bou.class.noaa.gov/saa/products/welcome

2.2 反演方法

在已有卫星云微物理反演研究(Rosenfeld et al, 1998, 2014a, 2014b, 2016; Zheng, et al, 2015a; Zhu et al, 2014, 2015)的基础上,改进云像元识别算法,增加对流云自动识别功能,开发了NPP/VIIRS卫星格点对流云云物理自动反演(Automatic Mapping of Convective Clouds,AMCC)软件。AMCC软件由云像元识别、re反演、Tb反演、对流云识别、T-re分析、NCCN反演和质量控制等模块组成。AMCC反演方法如下:

众所周知,VIIRS接收来自云顶附近云粒子反射的可见光和发射的红外辐射。为了研究云内的信息,用一定空间范围内不同发展高度的对流云簇顶部Tre,近似代替对流云内部不同高度上的Tre,从而得到这个区域对流云内不同高度上的云物理特征,即通过各态历经假定进行时、空转换。这种假设的云物理基础为:在一定空间区域内,大气的温度、湿度、气压、NCCN和上升气流等环境条件相近,湿空气上升冷却饱和凝结形成云,这些云簇的云底高度、云底温度、云底NCCN相近;相似的水汽、NCCN和动力条件下,不同气块按照绝热上升运动在相同高度上形成的云粒子大小相近,而对流云可近似为湿空气块的绝热上升运动,故用对流云簇顶部不同的Tre来近似对流云内不同高度上的Tre。当然,这个空间范围不能太大和太小,范围太大时会包含生成环境不相似的对流云,不满足假设条件;范围太小则每个高度层re像元数太少,也就缺少对这个区域内对流云的代表性(戴进等, 2010; Rosenfeld, et al, 1998)。目前,测试的结果为,使用像元数为96×96(在星下点的分辨率为36 km×36 km)区域内的云像元代表这一区域对流云簇,进行对流云簇的T-re分析,得到这个区域内对流云的平均云物理特征。使用同样的区域和分析方法遍历VIIRS数据,可得到卫星扫描区域的格点化云物理特征。一系列的飞机穿云观测(Andreae, et al, 2004; Prabha, et al, 2011; Rosenfeld et al, 1998, 2000; Braga, et al, 2017)表明,这种假设具有一定的合理性。

由于地面上行辐射会穿过薄云使反演结果失真,所以云掩膜方案的目标是剔除薄云,留下地面上行辐射不能透过的厚云作为研究对象。厚云掩膜方案的具体做法是:利用辐射传输模式,计算了不同地表(陆面或海面)反射率、太阳天顶角(Zsol)、卫星天顶角(Zsat)和太阳与卫星相对方位角(Arel)组合条件下厚云(云的光学厚度为10)的0.64 μm通道最小反射率(ρ0.6_min),建立ρ0.6_min查算表。通过VIIRS每个像元的ZsolZsatArel和下垫面类型,可在ρ0.6_min查算表中得到这个像元的ρ0.6_min。如果像元的0.64 μm通道反射率(ρ0.6h)大于这个像元的ρ0.6_min就被识别为厚云像元。这个云掩膜方案是对以前厚云掩膜方案的重要改进。以前的厚云掩膜方案是简单使用ρ0.6h大于0.4的阈值法,会过度剔除云边界像元,造成反演的云底温度过冷。改进后的方法与太阳、卫星的角度等有关,物理意义更为合理。

这些识别为厚云的像元中会包含有上层为薄云的多层云,这些像元的辐射亮温不能真实反映云顶信息,可通过分裂窗(高分辨11 μm和中分辨12 μm通道的亮温差)的动态阈值剔除这些多层厚云像元(Rosenfeld, et al, 2014b)。根据这些厚云像元3.745 μm高分辨率通道的辐射亮温反演3.745 μm的反射率(ρ3.7h),再利用ρ3.7hZsolZsatArel在已经建立的re查算表(戴进等, 2010; Nakajima, et al, 1990; Rosenfeld et al, 1998, 2014b)中得到re

一般认为,云滴在对流云云底附近核化,随上升气流开始凝结增长和碰并增长。因此,对流云云底像元具有最高的温度和最小的re。发展阶段对流云re随温度递减而增大的特征可用于识别对流云和非对流云(Lensky, et al, 1997)。对识别为对流云的96×96像元区域进行T-re分析,得到对流云T-re廓线和晶化温度(Tg)(戴进等, 2010; Rosenfeld, et al, 1998; 朱延年等, 2015)。选取温度最高云像元的温度作为对流云的Tb (Zhu, et al, 2014),利用FNL温度和位势高度的垂直廓线线性内插得到Hb。当云底re随着高度升高(温度降低)而增大到14 μm(T14)左右时降水过程开始启动(形成雨滴),从TbT14之间的距离称为降水启动厚度(D14),代表了降水发生的难易程度(Freud, et al, 2008; Rosenfeld, et al, 1998; Zhu, et al, 2015)。选取温度最低的像元可得到云顶温度(Ttop)和云顶海拔高度(Htop),Htop减去Hb和地面高度可得到云厚(Dcld)。

选取T-re廓线云底附近re小于20 μm(避免云滴产生明显的碰并过程)的最接近绝热凝结部分反演绝热云滴数浓度Nd, a(Rosenfeld, et al, 2014a)

(1)

式中,LWCa为气块从Tb开始绝热抬升凝结产生的绝热水含量,Mrv, a为绝热云滴质量,ρ为水的密度,rv为云滴平均体积半径。根据观测结果,这里取re =1.08 rv。考虑到实际云底的空气混合,对流云云底NCCN=Nd, a/1.3 (Brenguier, et al, 2011; Freud, et al, 2011)。

云底最大水汽过饱和度(Smax)计算式为

(2)

式中,C是与云底温度和云底气压相关的系数(Pinsky, et al, 2012),Wb为云底上升速度,Wb=0.0009Hb(Wb单位为m/s,Hb单位为m)(Zheng, et al, 2015a)。

AMCC反演对流云特征的格点分辨率约为36 km×36 km (0.33°×0.33°),每个格点的对流云宏微观物理特征量有TbHbWbSmaxNCCND14TtopHtopDcldTg等。卫星资料边缘的Zsat较大,这些区域云特征反演结果的精度较低,尤其在背对太阳一侧且Zsat较大的区域受云阴影的影响较大。为保证反演精度,文中只使用Zsat在-20°—50°(星下点左侧为负)范围的资料进行反演。

3 反演结果 3.1 反演云底温度的验证

Tb的反演在VIIRS云微物理反演中起非常关键的作用,Tb的准确与否会直接影响如WbSmaxNCCN等物理量反演的准确性。探空资料计算的抬升凝结高度可以近似看作云底高度。为了验证Tb反演结果的可靠性,利用那曲气象站(31.48°N,92.06°E)2013—2016年夏季06时第三次青藏高原大气科学试验加密探空资料计算抬升凝结温度(TLCL),与时间和位置相对应37个VIIRS反演Tb进行对比,结果见图 1。VIIRS反演Tb与抬升凝结温度(TLCL)较为吻合(相关系数为0.87,均方根误差为3.0℃),说明VIIRS反演高原对流云Tb具有较高的可靠性。

图 1 那曲气象站抬升凝结温度(TLCL)与VIIRS反演云底温度(Tb)的对比 Figure 1 Comparison between the temperature of lifting condensation level (TLCL) in Naqu meteorological station and the temperature of cloud base (Tb) from VIIRS
3.2 反演个例分析

2016年6月19日07时02分的NPP/VIIRS资料(卫星天顶角-20°—50°)的多光谱伪彩色合成(图 2)显示,青藏高原、印度东部和缅甸的大部分区域为对流云团覆盖(Rosenfeld, et al, 1998)。

图 2 2016年6月19日07时02分NPP/VIIRS多光谱RGB伪彩色合成 (红:中分辨率0.6 μm通道可见光反射率,绿:中分辨率2.2 μm通道可见光反射率,蓝:中分辨率11 μm通道辐射亮温,卫星天顶角为-20°-50°,伪彩色的含义参见Rosenfeld等(1998)) Figure 2 VIIRS RGB composite image of 0.6 μm reflectance (red), 2.2 μm reflectance (green) and 11 μm brightness temperature (blue) at medium resolution on 07:02 UTC 19 June 2016. (The satellite zenith angle is between -20° and 50°. The meaning of RGB composite is after Rosenfeld, et al (1998))

VIIRS反演的对流云云底温度(Tb)(图 3a)与FNL地面2 m温度空间线性插值后计算的抬升凝结温度(图 3b)对比可见,反演的海面、陆地和高大地形等不同下垫面类型的Tb稳定、可靠。从图 3a可见,印度、不丹、孟加拉国、缅甸、泰国等地的对流云为暖底云,Tb集中在10—25℃。高原地区为冷底云,青海南部Tb约为-10℃,西藏中部为-5℃左右。拉萨—日喀则一带云底较暖,约为5℃。利用2016年6月19日06时加密探空计算拉萨和那曲气象站的抬升凝结温度分别为5.3和2.1℃,而07时02分VIIRS资料反演这两个站的Tb分别为5.0和2.0℃,两者一致。

图 3 2016年6月19日07时02分VIIRS反演的对流云云底温度(Tb)和FNL地面2-m温度计算的抬升凝结温度(TLCL_FNL)(单位:℃) (a.Tb,b. TLCL_FNL,图中实心彩色小方块的色彩表示TbTLCL_FNL的大小) Figure 3 Distributions of retrieved cloud base temperature (Tb) from VIIRS on 07:02 UTC 19 June 2016 and the temperature of lifting condensation level (TLCL_FNL) calculated by FNL 2-m temperature over the surface (unit:℃) (a. Tb, b. TLCL_FNL, and the color represents the value of the Tb and TLCL_FNL)

从反演的对流云云底NCCN分布(图 4)可见,缅甸中部NCCN最高,部分地区超过了2000个/mg。孟加拉国和印度布拉马普特拉河流域NCCN较高,为1000个/mg左右,不丹和高原偏西部NCCN大部分为500个/mg左右,高原东部NCCN很小(50—200个/mg),和黄河上游河南县、玛曲地区NCCN的观测平均值(8.2—81.3个/cm3)相当(黄庚等,2002),这主要可能与人类活动和污染较少有关。而拉萨附近NCCN为1000个/mg左右。反映出高NCCN与人口、城市和工业密集区一致的特点。受青海北部沙尘区影响,南部出现800—1100个/mg高NCCN区。虽然高原地区没有与此个例对应的NCCN实测值,但反演的NCCN分布基本与人们的认知一致,表明了AMCC具有反演高原对流云微物理特征的能力。

图 4 2016年6月19日07时02分VIIRS反演的云凝结核浓度 (NCCN,单位:个/mg, 其他同图 3) Figure 4 Distribution of cloud condensation nuclei concentration (NCCN, unit:mg-1) retrieved from VIIRS on 07:02 UTC 19 June 2016 (Others are the same as in Fig. 3)
3.3 高原云微物理量统计特征 3.3.1 统计样本数

为了使反演结果具有较好代表性,反演了2013—2017年5个夏季对流云的TbNCCNSmaxWb等物理量,对0.33°×0.33°区域求平均(仅用样本数不少于10的0.33°×0.33°区域),作为格点值。从每个格点上样本数分布(图 5)可见,高原西部格点样本数大部分在20个以上,高原东部格点样本数在80个以上,部分格点超过了100个。这表明夏季高原东部对流活跃,北部和西部对流较弱,与许多观测结果一致(常祎等,2016张晓等, 2015傅云飞等,2016)。

图 5 0.33°×0.33°格点样本数分布 (实心彩色小方块的色彩表示样本个数) Figure 5 Distribution of sample number at 0.33°×0.33°grids (The colors in solid squares represent the sample number)
3.3.2 高原云底特征

图 6为卫星反演云底温度(Tb)和云底离地高度(Hb)5 a平均分布。由图 6a可见,高原Tb约为-5℃,标准差4℃左右。其中柴达木盆地Tb为-2℃,高原边缘约为0℃,标准差为4℃;而印度、尼泊尔和孟加拉等地Tb为20—25℃,标准差2℃,四川盆地西部Tb为18℃,标准差6℃。Hb特征(图 6b)为:高原中部、南部较高,为1800—2200 m;雅鲁藏布大峡谷为1100 m,尼泊尔、不丹一带Hb最低,约为700 m;柴达木盆地Hb最高,为2500 m;四川盆地为1500 m左右。这与徐祥德等(2006)发现高原夏季平均Hb大于1.5 km、高于同期平原和盆地地区Hb几百米的结论一致。张晓等(2015)基于CloudSat卫星资料得到区域(29.5°—35°N,90°—100°E)云底海拔高度在6—7 km;常祎等(2016)使用云高仪观测那曲地区最低层云的Hb平均为2.38 km,与本研究结果大体相当。

图 6 0.33°×0.33°格点平均云底温度(Tb)(单位:℃)和云底离地高度(Hb)(单位:m) (a.Tb,b.Hb,图中实心彩色小方块的色彩表示平均TbHb值) Figure 6 Distributions of average cloud base temperature (Tb, unit:℃) and height of cloud base above the ground (Hb, unit:m) in 0.33°×0.33° grids (a. Tb, b. Hb, and the colors represent the value of average Tb and Hb)

由此可见,夏季高原Tb约为-5℃,海拨高度约为6000 m(气压约为550 hPa),云底水汽混合比约为4.6 g/kg;而内陆四川盆地为18℃,海拨高度约2000 m(气压约为850 hPa),水汽混合比约为15 g/kg。即高原对流云水汽含量仅为平原地区的1/3,且云的离地高度较高,云的主体均在0℃层高度以上发展,云物理过程以冷云过程为主。相对于暖底对流云,冷底对流云所凝结的水量较少,会造成云中的水成物少,释放潜热少,对流不稳定能量相对较小,雨量、雨强都不大,与傅云飞等(2016)高原深厚弱对流降水平均雨强约1.6 mm/h结论一致。

3.3.3 高原云凝结核浓度(NCCN)和最大过饱和度(Smax)特征

由于总体受人类活动影响较小,高原云底平均NCCN为200—400个/mg,标准差约200个/mg(图 7),Smax在0.7%左右(图 8)。拉萨周边的NCCN为500个/mg左右,标准差约260个/mg,Smax在0.6%左右。雅鲁藏布江流域及藏南地区为NCCN低值区,为100—200个/mg,标准差100—200个/mg,Smax为0.8%左右。青海柴达木盆地NCCN较高,为500—800个/mg,标准差200—800个/mg,可能与区域扬沙或沙尘输送有关。沿着印度—尼泊尔—高原从南向北,随着地势从低到高,NCCN从1000个/mg左右逐渐减小到200个/mg左右。

图 7 平均云凝结核浓度(NCCN)分布 (单位:个/mg,其他同图 6) Figure 7 Distribution of average cloud condensation nuclei concentration (NCCN, unit:mg-1; Others are the same as in Fig. 6)
图 8 平均云底最大过饱和度(Smax)分布 (单位:%),其他同图 Figure 8 Distribution of average cloud base maximum supersaturation (Smax, unit:%; Others are the same as in Fig. 6)

高原云底过饱和度明显高于周边平原和盆地(图 8)。高原中西部Smax为0.8%,北部和东部Smax略低,在0.6%左右;四川盆地西部Smax约为0.3%。暖云云滴凝结增长方程为,其中r为云滴半径,t为时间,G1在给定环境中为常数,S为过饱和度(约翰·M.华莱士等, 2008)。由此可知,r2S成正比,与中国平原地区相比,高原地区大的S使云滴具有更快凝结增长速率,云滴更易通过凝结增长变大,进而更容易开始碰并增长、形成降水。

3.3.4 高原云底上升速度(Wb)

高原对流云云底平均上升速度(Wb)≥1.5 m/s(图 9),雅鲁藏布江流域及藏南地区Wb比高原其他地区略小,约1 m/s。青海柴达木盆地Wb较高,为2—2.3 m/s。高原南侧印度一带的Wb明显低于高原,约为1 m/s。高原云底上升速度较周围大的原因可能是与云形成初期凝华释放的凝华潜热较大有关,因为冷云底再加上较低的环境温度,在云形成初期,尽管凝结量小,但凝华量较大(王宏等, 2002)。徐祥德等(2006)利用声雷达观测发现高原热泡上升速度异常,垂直速度达1 m/s。傅云飞等(2016)发现高原暖顶云降水回波中顶高超过6 km(环境温度低于0℃)的比例高,猜测高原浅对流云内上升气流强,将环境0℃以下的水相粒子和暖空气携带至此高度以上,在低于0℃环境中形成一暖块。图 9Wb分布部分验证了这个猜测。

图 9 平均云底上升速度 (Wb, 单位:m/s;其他同图 6) Figure 9 Distribution of average cloud base updraft (Wb, unit:m/s; others are the same as in Fig. 6)
3.3.5 降水启动厚度(D14)

从对流云平均降水启动厚度(D14)的分布(图 10)可见,高原上D14为1500—2000 m,与平原上4000—5000 m (Zhu, et al, 2015)相比,更容易形成降水。雅鲁藏布江大峡谷及藏南地区D14值明显低于高原其他地区,为500—1000 m,与洁净海洋上对流云相当,表明该区域极易形成降水。青海柴达木盆地D14高,为1500—2000 m。印度及云南北部、四川盆地D14明显高于高原,大部分都超过了2000 m,与NCCN高值区对应(图 7),反映出空气污染对云降水的影响,与D14高值区相对应对流云更难形成降水。

图 10 平均降水启动厚度 (D14,单位:m;其他同图 6) Figure 10 Distribution of average precipitation initialization depth (D14, unit:m; others are the same as in Fig. 6)
3.3.6 高原云顶海拔高度(Htop)和云厚(Dcld)

图 11为高原对流云平均云顶海拔高度(Htop)的分布。高原上Htop为10—13 km;雅鲁藏布江流域及藏南地区Htop值略低于高原其他地区,约10 km;柴达木盆地Htop约为8 km;高原南侧印度一带Htop为6—8 km。那曲5 a平均Htop为10.58 km,C波段连续波雷达和Ka波段毫米波云雷达观测那曲地区的2014年7—8月平均Htop分别为11.41和11.62 km,FY-2E卫星TBB数据计算的Htop只有9.18 km(常祎等, 2016),同时段VIIRS观测那曲26个个例的平均Htop为11.8 km。可见,VIIRS反演的Htop与雷达观测非常接近。

图 11 平均云顶海拔高度 (Htop,单位:km;其他同图 6) Figure 11 Distribution of average cloud top altitude (Htop, unit:km; others are the same as in Fig. 6)

对流云的厚度(Dcld)从高原南部的5000 m左右逐渐减小到北部的2500 m左右(图 12)。表明高原南部对流云发展旺盛,北部对流云发展较弱、较浅薄。有限云厚加之云内含水量低,使得高原上对流云发生强降水和灾害性天气过程概率大为降低,因降水启动厚度较薄,容易形成降水,因此,高原上降水频繁、雨量雨强小,且多以短时阵性降水为主。常祎等(2016)研究指出高原主要以短时阵性降水为主,持续时间小于1 h,远小于同时期中国中东部平原地区。这与高原上云底、云厚、降水启动厚度特征有关。

图 12 平均云厚 (Dcld,单位:m;其他同图 6) Figure 12 Average depth between cloud base and cloud top (Dcld, unit:m; others are the same as in Fig. 6)
3.3.7 晶化温度(Tg)和晶化厚度(DTg)

晶化温度(Tg)即云粒子几乎全部冰晶化的温度(Rosenfeld et al, 1998, 2011; 朱延年等, 2015)。晶化厚度(DTg)指从云底到Tg之间的厚度。高原中部、南部Tg低于-30℃,北部及青海境内Tg为-25℃左右(图 13),这可能是受北部沙尘气溶胶成冰作用影响。DTgDcld的分布基本一致(图 1412),高原南部DTg在4000 m左右,北部在2000 m左右,高原外围DTg大于6000 m。可见高原对流云比平原地区更容易成冰,加之高原Tb约为-5℃,使得云内降水粒子以冰相为主。当降水发生时,造成降水粒子掉出云底时仍以固态为主,及地时才完全融化,因而来不及破碎,呈现出滴大数浓度小的特点。Chen等(2013)观测高原夏季雨滴谱的谱宽大于同纬度平原地区,常祎等(2016)观测到易出现较大的雨滴和冰相粒子,用VIIRS反演微物理特征可以解释物理成因。

图 13 平均晶化温度 (Tg,单位:℃;其他同图 6) Figure 13 Average glaciation temperature (Tg, unit:℃; others are the same as in Fig. 6)
图 14 平均晶化厚度(DTg,单位:m)分布 (其他同图 6) Figure 14 Average depth between cloud base and glaciation (DTg, unit:m; others are the same as in Fig. 6)
4 结论与讨论

利用基于VIIRS高分辨率资料的格点对流云云物理自动反演软件(AMCC),反演高原5个夏季对流云微物理特征,得出如下结论:

(1) 通过VIIRS反演云底温度与那曲站探空计算抬升凝结温度对比分析,两者线性相关,相关系数为0.87,均方根误差为3.0℃。表明了AMCC反演Tb的可用性。

(2) 高原对流云云宏微观物理特征为,一是云底冷,Tb为-5℃,Hb为1800—2200 m,云内含水量小;二是云底NCCN为200—400个/mg,Smax为0.7%,NCCN少、Smax大,云滴凝结增长速率更快;三是D14小,为1500—2000 m,雅鲁藏布江流域及藏南地区D14约500—1000 m,更加容易形成降水;四是Htop为10—13 km、Dcld有限,从南部5000 m逐渐减小到北部2500 m;五是Tg高,从中部、南部-30℃到北部-25℃,加之高原Tb<0℃,使得云内降水粒子以冰相为主。

(3) 高原对流云微物理特征决定了其降水具有多发、短时、量小、滴大的特点。另外,有限云厚和云内含水量低,使得发生强降水和灾害性天气概率小。

利用NPP/VIIRS反演高原地区对流云多个云物理量,据此得到云物理特征,既能补充高原观测不足,又能发挥卫星多通道资源优势,是值得继续深入研究的工作。同时需要指出的是,文中AMCC反演还存在着不足,一是卫星反演云底高度没有用云高仪校验,而是用探空资料计算的抬升凝结高度检验。按照气块法理论,仅当对流云是从地面经充分混合上升凝结而成时,抬升凝结高度才能表征云底高度,否则代表性就较差。这一点从VIIRS反演的Tb与那曲站的TLCL对比结果可以看出。二是有些物理量还没有经过地基和飞机探测校验,如:云底NCCN、降水启动厚度、晶化温度等。三是自动化反演中的质量控制也存在不完善之处。这些都需要在今后工作中,加强外场观测,用更多中外的地基和机载观测资料校验反演结果,据此不断发展反演算法,逐步提高反演准确性,才能更好发挥卫星探测优势。

致谢: 本研究使用了NOAA CLASS的VIIRS L1B数据、NCEP FNL数据和中国气象局高原试验加密探空资料,在此一并表示感谢。
参考文献
常祎, 郭学良. 2016. 青藏高原那曲地区夏季对流云结构及雨滴谱分布日变化特征. 科学通报, 61(15): 1706–1720. Chang Y, Guo X L. 2016. Characteristics of convective cloud and precipitation during summer time at Naqu over Tibetan Plateau. Chinese Sci Bull, 61(15): 1706–1720. (in Chinese)
戴进, 余兴, 刘贵华, 等. 2010. 一次暴雨过程中云微物理特征的卫星反演分析. 气象学报, 68(3): 387–397. Dai J, Yu X, Liu G H, et al. 2010. Analyses of satellite retrieval microphysical properties of a rainstorm in the northern part of Shaanxi. Acta Meteor Sinica, 68(3): 387–397. (in Chinese)
傅云飞, 潘晓, 刘国胜, 等. 2016. 基于云亮温和降水回波顶高度分类的夏季青藏高原降水研究. 大气科学, 40(1): 102–120. Fu Y F, Pan X, Liu G S, et al. 2016. Characteristics of precipitation based on cloud brightness temperatures and storm tops in summer Tibetan Plateau. Chinese J Atmos Sci, 40(1): 102–120. (in Chinese)
黄庚, 李淑日, 德力格尔, 等. 2002. 黄河上游云凝结核观测研究. 气象, 28(10): 45–49. Huang G, Li S R, Deligeri, et al. 2002. Observational analysis of content of condensation nuclei in the atmosphere in the upper reaches of Huanghe River. Meteor Mon, 28(10): 45–49. (in Chinese)
刘贵华, 余兴, 师春香, 等. 2011a. FY-3A/VIRR反演云微物理特征及与TERRA/MODIS反演结果的比较. 高原气象, 30(2): 461–470. Liu G H, Yu X, Shi C X, et al. 2011a. Comparisons of microphysical property of cloud retrieved from FY-3A/VIRR and TERRA/MODIS. Plateau Meteor, 30(2): 461–470. (in Chinese)
刘贵华, 余兴, 戴进, 等. 2011b. 地形云人工增雨条件卫星探测反演个例分析. 气象学报, 69(2): 363–369. Liu G H, Yu X, Dai J, et al. 2011b. A case study of the conditions for topographic cloud seeding based on the retrieval of satellite measurements. Acta Meteor Sinica, 69(2): 363–369. (in Chinese)
刘建军, 陈葆德. 2017. 基于CloudSat卫星资料的青藏高原云系发生频率及其结构. 高原气象, 36(3): 632–642. Liu J J, Chen B D. 2017. Cloud occurrence frequency and structure over the Qinghai-Tibetan Plateau from CloudSat observation. Plateau Meteor, 36(3): 632–642. (in Chinese)
刘黎平, 郑佳锋, 阮征, 等. 2015. 2014年青藏高原云和降水多种雷达综合观测试验及云特征初步分析结果. 气象学报, 73(4): 635–647. Liu L P, Zheng J F, Ruan Z, et al. 2015. The preliminary analyses of the cloud properties over the Tibetan Plateau from the field experiments in clouds precipitation with the vavious radars. Acta Meteor Sinica, 73(4): 635–647. (in Chinese)
王宝鉴, 黄玉霞, 魏栋, 等. 2017. TRMM卫星对青藏高原东坡一次大暴雨强降水结构的研究. 气象学报, 75(6): 966–980. Wang B J, Huang Y X, Wei D, et al. 2017. Structure analysis of heavy precipitation over the eastern slope of the Tibetan Plateau based on TRMM data. Acta Meteor Sinica, 75(6): 966–980. (in Chinese)
王宏, 雷恒池, 德力格尔, 等. 2002. 黄河上游地区强对流云特征的模拟分析. 气候与环境研究, 7(4): 397–408. Wang H, Lei H C, Deligeri, et al. 2002. A numerical simulation of characteristics of convective cloud at the upper reaches of the Yellow River. Climatic Environ Res, 7(4): 397–408. (in Chinese)
汪会, 罗亚丽, 张人禾. 2011. 用CloudSat/CALIPSO资料分析亚洲季风区和青藏高原地区云的季节变化特征. 大气科学, 35(6): 1117–1131. Wang H, Luo Y L, Zhang R H. 2011. Analyzing seasonal variation of clouds over the Asian monsoon regions and the Tibetan Plateau region using CloudSat/CALIPSO data. Chinese J Atmos Sci, 35(6): 1117–1131. DOI:10.3878/j.issn.1006-9895.2011.06.11 (in Chinese)
王黎俊, 银燕, 姚展予, 等. 2013. 三江源地区秋季一次层积云飞机人工增雨催化试验的微物理响应. 气象学报, 71(5): 925–939. Wang L J, Yin Y, Yao Z Y, et al. 2013. Microphysical responses as seen in a stratocumulus aircraft seeding experiment in autumn over the Sanjiangyuan national nature reserve. Acta Meteor Sinica, 71(5): 925–939. (in Chinese)
吴翀, 刘黎平, 翟晓春. 2017. Ka波段固态发射机体制云雷达和激光云高仪探测青藏高原夏季云底能力和效果对比分析. 大气科学, 41(4): 659–672. Wu C, Liu L P, Zhai X C. 2017. The comparison of cloud base observations with Ka-band solid-state transmitter-based millimeter wave cloud radar and ceilometer in summer over Tibetan Plateau. Chinese J Atmos Sci, 41(4): 659–672. (in Chinese)
夏昕, 任荣彩, 吴国雄, 等. 2016. 青藏高原周边对流层顶的时空分布、热力成因及动力效应分析. 气象学报, 74(4): 525–541. Xia X, Ren R C, Wu G X, et al. 2016. An analysis on the spatio-temporal variations and dynamic effects of the tropopause and the related stratosphere-troposphere coupling surrounding the Tibetan Plateau area. Acta Meteor Sinica, 74(4): 525–541. (in Chinese)
徐祥德, 陈联寿. 2006. 青藏高原大气科学试验研究进展. 应用气象学报, 17(6): 756–772. Xu X D, Chen L S. 2006. Advances of the study on Tibetan Plateau experiment of atmospheric sciences. J Appl Meteor Sci, 17(6): 756–772. DOI:10.3969/j.issn.1001-7313.2006.06.013 (in Chinese)
徐祥德, 赵天良, 施晓晖, 等. 2015. 青藏高原热力强迫对中国东部降水和水汽输送的调制作用. 气象学报, 73(1): 20–35. Xu X D, Zhao T L, Shi X H, et al. 2015. A study of the role of the Tibetan Plateau's thermal forcing in modulating rainband and moisture transport in eastern China. Acta Meteor Sinica, 73(1): 20–35. (in Chinese)
约翰·M.华莱士, 彼得·V.霍布斯. 2008.大气科学.何金海, 王振会, 银燕, 译. 2版.北京: 科学出版社, 233 pp. Wallace J M, Hobbs P V. 2008. Atmospheric Science. He J H, Wang Z H, Yin Y, trans. 2nd ed. Beijing: Science Press, 233pp (in Chinese)
张鸿发, 郭三刚, 张义军, 等. 2003. 青藏高原强对流雷暴云分布特征. 高原气象, 22(6): 558–564. Zhang H F, Guo S G, Zhang Y J, et al. 2003. Distribution characteristic of severe convective thunderstorm cloud over Qinghai-Xizang Plateau. Plateau Meteor, 22(6): 558–564. DOI:10.3321/j.issn:1000-0534.2003.06.005 (in Chinese)
张倩, 宋小全, 刘金涛, 等. 2016. 激光云高仪探测青藏高原夏季云底高度的研究. 光电子·激光, 27(4): 406–412. Zhang Q, Song X Q, Liu J T, et al. 2016. Observation of cloud base height with ceilometer in Tibetan Plateau during summer. J Optoelectro·Laser, 27(4): 406–412. (in Chinese)
张晓, 段克勤, 石培宏. 2015. 基于CloudSat卫星资料分析青藏高原东部夏季云的垂直结构. 大气科学, 39(6): 1073–1080. Zhang X, Duan K Q, Shi P H. 2015. Cloud vertical profiles from CloudSat data over the eastern Tibetan Plateau during summer. Chinese J Atmos Sci, 39(6): 1073–1080. (in Chinese)
朱延年, 余兴, 徐小红, 等. 2015. 利用卫星分析对流云成冰能力与繁生机制及气溶胶影响. 高原气象, 34(6): 1758–1764. Zhu Y N, Yu X, Xu X H, et al. 2015. Glaciation and ice multiplication of convective clouds and their dependence on aerosol investigated by satellites. Plateau Meteor, 34(6): 1758–1764. (in Chinese)
Andreae M O, Rosenfeld D, Artaxo P, et al. 2004. Smoking rain clouds over the Amazon. Science, 303(5662): 1337–1342. DOI:10.1126/science.1092779
Braga R C, Rosenfeld D, Weigel R, et al. 2017. Further evidence for CCN aerosol concentrations determining the height of warm rain and ice initiation in convective clouds over the Amazon basin. Atmos Chem Phys, 17(23): 14433–14456. DOI:10.5194/acp-17-14433-2017
Brenguier J L, Burnet F, Geoffroy O. 2011. Cloud optical thickness and liquid water path:Does the k coefficient vary with droplet concentration?. Atmos Chem Phys, 11(18): 9771–9786. DOI:10.5194/acp-11-9771-2011
Chen B J, Yang J, Pu J P. 2013. Statistical characteristics of raindrop size distribution in the Meiyu season observed in eastern China. J Meteor Soc Japan, 91(2): 215–227. DOI:10.2151/jmsj.2013-208
Freud E, Rosenfeld D, Kulkarni J R. 2011. Resolving both entrainment-mixing and number of activated CCN in deep convective clouds. Atmos Chem Phys, 11(24): 12887–12900. DOI:10.5194/acp-11-12887-2011
Freud E, Ström J, Rosenfeld D, et al. 2008. Anthropogenic aerosol effects on convective cloud microphysical properties in southern Sweden. Tellus B, 60(2): 286–297. DOI:10.1111/j.1600-0889.2007.00337.x
Hillger D, Kopp T, Lee T, et al. 2013. First-light imagery from suomi NPP ⅦRS. Bull Amer Meteor Soc, 94(7): 1019–1029. DOI:10.1175/BAMS-D-12-00097.1
Lensky I M, Rosenfeld D. 1997. Estimation of precipitation area and rain intensity based on the microphysical properties retrieved from NOAA AVHRR data. J Appl Meteor, 36(3): 234–242. DOI:10.1175/1520-0450(1997)036<0234:EOPAAR>2.0.CO;2
Merk D, Deneke H, Pospichal B, et al. 2016. Investigation of the adiabatic assumption for estimating cloud micro-and macrophysical properties from satellite and ground observations. Atmos Chem Phys, 16(2): 933–952. DOI:10.5194/acp-16-933-2016
Miller D J, Zhang Z B, Ackerman A S, et al. 2016. The impact of cloud vertical profile on liquid water path retrieval based on the bispectral method:A theoretical study based on large-eddy simulations of shallow marine boundary layer clouds. J Geophys Res, 121(8): 4122–4141.
Min Q, Joseph E, Lin Y, et al. 2012. Comparison of MODIS cloud microphysical properties with in-situ measurements over the Southeast Pacific. Atmos Chem Phys, 12(23): 11261–11273. DOI:10.5194/acp-12-11261-2012
Nakajima T, King M D. 1990. Determination of the optical thickness and effective particle radius of clouds from reflected solar radiation measurements.Part Ⅰ:Theory. J Atmos Sci, 47(15): 1878–1893. DOI:10.1175/1520-0469(1990)047<1878:DOTOTA>2.0.CO;2
Pinsky M, Khain A, Mazin I, et al. 2012. Analytical estimation of droplet concentration at cloud base. J Geophys Res, 117(D18): D18211.
Prabha T V, Khain A, Maheshkumar R S, et al. 2011. Microphysics of premonsoon and monsoon clouds as seen from in situ measurements during the cloud aerosol interaction and precipitation enhancement experiment (CAIPEEX). J Atmos Sci, 68(9): 1882–1901. DOI:10.1175/2011JAS3707.1
Rosenfeld D, Lensky I M. 1998. Satellite-based insights into precipitation formation processes in continental and maritime convective clouds. Bull Amer Meteor Soc, 79(11): 2457–2476. DOI:10.1175/1520-0477(1998)079<2457:SBIIPF>2.0.CO;2
Rosenfeld D, Woodley W L. 2000. Deep convective clouds with sustained supercooled liquid water down to -37.5℃. Nature, 405(6785): 440–442. DOI:10.1038/35013030
Rosenfeld D, Woodley W L, Lerner A, et al. 2008. Satellite detection of severe convective storms by their retrieved vertical profiles of cloud particle effective radius and thermodynamic phase. J Geophys Res, 113(D4): D04208. DOI:10.1029/2007JD008600
Rosenfeld D, Yu X, Liu G H, et al. 2011. Glaciation temperatures of convective clouds ingesting desert dust, air pollution and smoke from forest fires. Geophys Res Lett, 38(21): L21804.
Rosenfeld D, Fischman B, Zheng Y T, et al. 2014a. Combined satellite and radar retrievals of drop concentration and CCN at convective cloud base. Geophys Res Lett, 41(9): 3259–3265. DOI:10.1002/2014GL059453
Rosenfeld D, Liu G, Yu X, et al. 2014b. High-resolution (375 m) cloud microstructure as seen from the NPP/ⅦRS satellite imager. Atmos Chem Phys, 14(5): 2479–2496. DOI:10.5194/acp-14-2479-2014
Rosenfeld D, Zheng Y T, Hashimshoni E, et al. 2016. Satellite retrieval of cloud condensation nuclei concentrations by using clouds as CCN chambers. Proc Natl Acad Sci USA, 113(21): 5828–5834. DOI:10.1073/pnas.1514044113
Wan B C, Gao Z Q, Chen F, et al. 2017. Impact of Tibetan Plateau surface heating on persistent extreme precipitation events in southeastern China. Mon Wea Rev, 145(9): 3485–3505. DOI:10.1175/MWR-D-17-0061.1
Williams E, Stanfill S. 2002. The physical origin of the land-ocean contrast in lightning activity. C R Phys, 3(10): 1277–1292. DOI:10.1016/S1631-0705(02)01407-X
Zheng Y T, Rosenfeld D. 2015a. Linear relation between convective cloud base height and updrafts and application to satellite retrievals. Geophys Res Lett, 42(15): 6485–6491. DOI:10.1002/2015GL064809
Zheng Y T, Rosenfeld D, Li Z Q. 2015b. Satellite inference of thermals and cloud-base updraft speeds based on retrieved surface and cloud-base temperatures. J Atmos Sci, 72(6): 2411–2428. DOI:10.1175/JAS-D-14-0283.1
Zhu Y N, Rosenfeld D, Yu X, et al. 2014. Satellite retrieval of convective cloud base temperature based on the NPP/ⅦRS Imager. Geophys Res Lett, 41(4): 1308–1313. DOI:10.1002/2013GL058970
Zhu Y N, Rosenfeld D, Yu X, et al. 2015. Separating aerosol microphysical effects and satellite measurement artifacts of the relationships between warm rain onset height and aerosol optical depth. J Geophys Res, 120(15): 7726–7736.