2. 中国气象科学研究院,北京 100081
2. Chinese Academy of Meteorological Sciences, Beijing 100081
在过去几十年中,测量云中结构和水成物的新技术与装备已有很大发展(NRC,2003;雷恒池等,2008;郭学良等,2013)。作为云物理学研究的重要组成部分与重要观测手段,飞机云物理探测特别是机载云物理探测系统应用,因具有穿云直接观测的特点,对于了解大气环境特征、云雾降水形成发展过程、人工影响天气科学指挥作业、天气气候变化与精细化数值预报研究以及飞行安全(如积冰)等方面都起着至关重要的作用,并被国内外广泛地应用于地球物理(大气)科学各类观测研究中以及指导人工影响天气作业实践(范慧君和马培民,1981;NRC,2003;毛节泰和郑国光,2006;IPCC,2007)。
20世纪70年代,国际上就开始广泛使用美国PMS(Particle Measuring Systems)公司的光电探头进行云微物理参数探测(Knollenberg,1981;游来光,1987;魏强等,1997;Lawson and Jensen, 1998;郭学良和郑国光,2010)。这些粒子测量设备作为探测云物理结构最常见仪器设备能够自动完成对粒子的取样、测量和记录,具有取样频率高、取样尺度分辨率高等特点,可对大范围云系进行连续观测。20世纪80年代以后中国气象科学研究院先后主持或参与改装了多架不同型号的云物理观测和人工增雨播云作业飞机。飞机装载各种气象和云物理探测仪器进入云中直接测量,取得了十多个省(区)的云降水微物理结构和大气核的资料(姚展予,2006)。“十五”期间,国家进一步加大对人工影响天气研究试验支持,中国气象科学研究院集合新引进PMS粒子测量系统与自研的空地数传系统、碘化银末端燃烧器,主持改装完成了我国新一代人工增雨综合探测作业飞机夏延-ⅢA设备与技术集成(杨绍忠和陶文有,2002;刘卫国等,2003;王广河和姚展予,2003)。之后国内陆续引进多套粒子测量系统及其他机载云物理探测仪器,先后改装建设了系列化作业飞机和探测飞机,并在我国西北、东北和华北等数个省(市)组织进行了飞机云物理结构探测和人工影响天气作业试验(雷恒池等,2008;郭学良等,2013;王磊等,2014;蔡兆鑫等,2019)。随着人工影响天气业务科学化、规范化、现代化建设加快,人工增雨雪逐步向飞机作业为主的空中平台转变,空地一体化作业对飞机作业平台和配套机载探测能力、播撒能力和指挥调度能力的要求也越来越高。
作为空中催化播撒最主要的平台和云物理探测最直接的手段,高性能的人工影响天气专用飞机及其机载云物理探测能力建设对于促进人工影响天气现代业务发展起着举足轻重的作用。2012年5月国家发展改革委批复立项“新增千亿斤粮食工程:东北区域人工影响天气能力建设项目”(以下简称“东北人影工程”),作为核心建设任务及项目最重要、最关键的技术装备,建设的两架新舟60国家高性能增雨飞机既可满足不同目的作业需要,又能作为空中试验平台开展大气科学等相关试验任务。2014年12月国家发展改革委联合中国气象局共同编制颁布了《全国人工影响天气发展规划(2014—2020年)》,提出大力加强飞机作业力量,加快国家高性能作业飞机建设,集成先进探测设备、催化播撒设备和空地通信设备,建成具有国际先进水平的云物理飞机探测平台。伴随国家和地方人工影响天气能力建设,大量先进机载探测设备投入使用,对深入开展云降水物理科学研究和开拓新的研究领域均发挥出越来越重要作用,并对人工影响天气业务科技水平提高起到了重要推动作用(Guo and Zheng, 2009)。
目前全国人工影响天气飞机平台与机载云物理探测设备建设应用相关技术标准和规范仍存在很多不足。与先进国家同类工作相比,人工影响天气飞机机载云物理探测系统应用中所存在的问题主要体现在进口机载探测设备使用过程中受标定、维护、培训等诸多条件制约,实际业务环境应用中部分探测设备的稳定性、可靠性以及数据质量也逐渐显露出了一些问题。本文针对“东北人影工程”建设的新舟60国家高性能增雨飞机,结合该飞机2018年10月21日在湖北省开展的一次飞行案例,通过飞机平台自身设备与加装的机载任务系统设备对一些关键飞行参数和环境气象要素进行观测对比分析,为后续国家和地方高性能人工影响天气飞机建设机载任务系统集成设计和机载探测数据合理应用提供重要技术支持。
1 飞机观测平台与功能 1.1 机载探测重要作用人工影响天气是建立在云雾降水物理学基础上的一门应用科学技术。监测了解云雾降水形成和发展中各种物理过程及其宏观与微观特征变化,是揭示成云致雨机制、发展人工催化原理与技术方法的基础。作为云降水物理学研究的重要组成部分和空中催化播撒最主要平台,利用人工影响天气专用飞机及其配备的先进机载云物理探测系统进行穿云采样,它是获取云微物理特征和监测云降水演变最直接同时也是最有效的观测手段。在人工影响天气工作中,结合催化剂不同环境条件下的活化特性,充分利用高性能飞机及机载云物理探测设备实时准确地获取飞行位置及其周围环境基本气象要素以及云微物理特征量的时空分布变化,对突破和解决催化条件(表 1)监测识别、精准确定催化要素以及科学检验催化效果等诸项关键技术问题都具有重要支撑作用。
新舟60国家高性能增雨飞机(图 1)是“东北人影工程”的核心建设内容。自2012年8月飞机建设项目正式启动,2015年12月加入业务运行,同时参与系列重大科学试验和气象服务保障任务。该飞机也是我国首次拥有具有适航资格的专业化大中型人工增雨专用飞机。
新舟60增雨飞机机上构建有一套任务系统,将加装的大气探测、催化作业、宏观成像、空地通信等子系统的任务设备高度集成化和自动化(表 2)。在飞机任务系统中,大气探测子系统(亦为粒子测量系统)与宏观成像系统相配合可实时连续测量大气气溶胶、云凝结核、云降水(云滴、雨滴、冰雪晶等)粒子谱分布,云和降水粒子形状或图像,云中液态/固态水含量,大气温压湿和风等环境气象要素监测数据,飞行高度、经纬度、空速、航姿和航向等飞机位置参数,云宏观尺度、云中降水粒子相态、飞机积冰状况和飞机催化作业情况等云宏观影像。另外,大气探测子系统配套的粒子分析与显示系统(PADS),集成了卫星导航定位、云粒子测量系统、气象数据监测系统等多参量数据采集、处理和分析显示一体,用于对大气探测子系统数据的收集、处理、存储,对各类探测仪器工作参数和采样数据的基本显示及数据回放。
在飞机任务系统建设应用中,一方面需要充分考虑各子系统对引用的重要飞行参数与气象要素相同参数的验证(廖捷和任芝花,2011;赵增亮等,2011),并加以优选集成设计,另一方面也需考虑为同一基本参数增加必要备份,避免实际使用中这些基本参数的观测缺失,为探测试验和人工影响天气指挥作业提供基础保障。
2 机载不同观测设备及其性能比较 2.1 飞机导航定位飞机导航定位主要用于确定飞机位置信息,通常包括经度、纬度和海拔高度等重要飞行参数,它是航空管制与飞行安全必不可少的组成部分,同时也是了解云物理特征时空分布和精准识别催化部位并实施播云作业的前提。新舟60增雨飞机目前机上同时存在三套定位源:(1)飞机平台自身全球卫星定位系统(GPS),其天线安装在机头顶部正上方;(2)任务系统大气探测子系统中机载综合气象测量系统(AIMMS-20),其自含GPS全球卫星定位系统;(3)任务系统空地通信子系统加装的中国自主建立的北斗卫星导航系统(BDS),其天线安装在机身前部顶部。此外,飞机飞行高度值还可由飞机平台自身气压高度表来测定。
任务系统加装的加拿大Aventech Research Inc.公司制造的机载综合气象测量系统AIMMS-20(图 2),该系统由舱内、舱外两部分构成。舱内设备有全球定位系统(GPS)模块、惯导测量单元(IMU)、中央处理模块(CPM)。舱外设备包括外挂在左机翼下方的大气数据探头(ADP)和分别安装在左右机翼末端的GPS天线。该设备所含GPS模块具备精密三维定位能力,能够实时获取飞机惯性加速度及角速率,提供飞机位置、速度、姿态(滚转角、俯仰角和偏航角)和精密授时。
新舟60增雨飞机平台自备大气数据系统(ADS)用来测量、计算以及显示机体周围自由气流的大气压力(包括静压和全压)、总温、大气密度、气压高度、升降速度、马赫数、指示空速、真空速、攻角和侧滑角等重要飞行参数及环境气象要素,这些大气数据是飞行控制系统、导航系统、指示系统等不可或缺的信息,ADS的性能和精度也直接关系到飞行安全。在诸多大气数据中,包括全压、静压、总温、攻角和侧滑角等最基本参数直接通过ADS安装在飞机机头外部和机身上的空速管、静压孔、总温传感器和迎角传感器采集获取,其他飞行参数则通过对基本大气数据进一步解算加以获得。
任务系统中的AIMMS-20,除了具备对飞机导航定位功能外,其自带的大气压力、温度、湿度传感器还为登机人员指挥作业提供了环境气象条件实时测量参考。AIMMS-20环境气象要素测量的主要性能参数如表 3所示。
图 3为加装在增雨飞机右机翼下方由美国DMT(Droplet Measurement Technologies Inc.)公司制造的组合式云粒子图像探头(CIP)。该探测设备除了能够实时探测云粒子谱和显示粒子二维成像外,其还自带了对气压、温度、湿度等环境气象要素测量的传感器。通过CIP获取环境气压与温度,可转换得到真空速,一方面为本探头及其他粒子探头提供粒子数浓度计算用,如式(1)所示,另一方面这些要素也可作为AIMMS-20相同环境参数的备份,同时为指挥人员提供催化条件识别的关键参数,用以指导播云作业。
$ N_{\mathrm{T}}=\frac{\sum\limits_{i=1}^{m} n_{i}}{S \cdot T A S \cdot t} $ | (1) |
式中:NT为总的粒子数浓度,单位为个·cm-3;ni为第i通道测得粒子数;m为总测量通道数;S为探头采样面积;TAS为真空速;t为采样时间。
3 飞行参数不同观测对比分析结合2018年10月21日一次飞行案例,以下对新舟60增雨飞机平台GPS、ADS与飞机任务系统BDS、AIMMS-20、CIP等不同观测源获取的飞机位置及真空速等关键飞行参数及环境气象要素进行简单的对比分析。
3.1 飞机定位观测比较图 4为2018年10月21日新舟60增雨飞机(机号B-3726)在湖北省的一次飞行航迹显示。飞机定位信息取自任务系统AIMMS-20。当日12:21:53飞机从襄阳刘集机场起飞,至16:28:47落地,共飞行约4 h。飞机起飞后快速爬升到适宜催化高度执行人工增雨作业,随后下降进入设定目标区,并在该目标区内按300 m一层分层先上后下对云系的上下部分别做探测。探测零度层高度约为3 900 m,飞行中还观测到较强程度的结冰。
飞行高度随时间变化显示,飞机起飞后向上一直爬升,12:39:23达到4 650 m,然后维持在该高度做平直飞行,机上人员开展增雨播撒作业。13:33:54飞机再向上爬升,13:36:30上升到4 939 m并保持高度继续平飞。13:45:24后飞行下降,约6 min后于13:51:20降至3 050 m,然后维持高度水平探测。自14:17:06开始,飞机在设定的目标区内按垂直300 m一层分层盘旋做垂直探测,并向云系上部爬升,最高飞行至4 640 m附近。15:11:00后飞机再次快速下降高度,15:16:54飞行高度降为2 770 m,随后继续每300 m一层分层盘旋向下垂直探测到1 840 m。完成目标区内垂直探测后,飞机于16:19:00返场,并在16:28:47落地。
以任务系统AIMMS-20定位信息为基准,图 5给出了该飞行案例飞机平台GPS与任务系统BDS二者相对AIMMS-20定位经纬度的偏差。虽然上述定位设备机上安装位置有所差异,但总体上看三者经纬度定位较为接近,经纬度定位偏差较小,飞机平台GPS比BDS定位经纬度值更接近AIMMS-20。比较经度定位偏差,飞机平台GPS与AIMMS-20经度最大、最小偏差分别为0.003°、-0.003 1°,平均偏差为0.000 1°,小于10 m;BDS与AIMMS-20经度最大、最小偏差分别为0.003 5°、-0.003 2°,平均偏差也为0.000 1°。纬度定位偏差比较显示,飞机平台GPS和AIMMS-20最大、最小偏差分别为0.001 8°、-0.002 1°,BDS和AIMMS-20最大、最小偏差分别为0.002 5°、-0.002 6°,飞机平台GPS、任务系统BDS与AIMMS-20纬度定位平均偏差均为0°,非常相近。
从飞机自带气压高度表和任务系统BDS分别与AIMMS-20定位的海拔高度偏差(图 6)易看出,飞行中飞机气压高度表观测海拔高度明显低于BDS和AIMMS-20测量值,且随飞行高度升高,气压高度表观测海拔高度偏差也越大,与飞行高度基本呈反相关。这与气压高度表高空气压测量和压高转换计算有关。在2018年10月21日飞行案例中,飞机气压高度表海拔高度值比AIMMS-20探测高度至少低21.4 m,最多时偏低151.2 m(约3%)。
图 6中同时给出了BDS与AIMMS-20探测的海拔高度偏差。相对AIMMS-20,BDS定位高度最大高出43.1 m,最小低55.7 m(约6%)。二者高度偏差总体较小,特别是在飞机飞行高度较高时,二者更为接近。
3.2 真空速观测比较图 7显示了飞机平台ADS与任务系统AIMMS-20、CIP探头对真空速TAS观测比较。TAS是飞机相对于空气的真实速度,其用于粒子测量系统采样体积计算,其测量准确程度直接关系着基于该测量值所有机载探头采样粒子数浓度计算准确度。由图 7可见,飞机平台ADS与任务系统AIMMS-20、CIP探头真空速观测变化趋势相同,飞机平台ADS与AIMMS-20观测非常接近,二者相关系数达到0.99;CIP探测真空速则明显低于ADS和AIMMS-20。此飞行案例中,ADS的TAS比AIMMS-20平均偏快0.8 m·s-1,最高偏快8.0 m·s-1,最低则偏慢10.9 m·s-1,且均发生在起飞或降落阶段。二者差异较大时均发生在AIMMS-20秒记录的瞬时风速出现陡增或陡降的时候,这表明AIMMS-20对环境风速瞬时变化响应要比飞机平台ADS观测更灵敏。同AIMMS-20相比,CIP探测真空速平均偏慢约10 m·s-1,二者相关系数亦为0.99。
另外,任务系统AIMMS-20飞行中真空速偶尔会出现奇异值,如13:31:40、16:18:10机载AIMMS20探测TAS分别为46.1和48.6 m·s-1,而对应时刻前后秒记录的TAS值分别大于100和87 m·s-1,实际飞行中基本不可能发生。因此,出现的明显低于前后秒记录数值应为仪器观测异常需加以剔除。
4 环境气象要素不同观测对比分析 4.1 气压观测比较图 8显示,新舟60增雨飞机任务系统AIMMS-20与CIP探测的环境气压二者变化趋势完全一致。在2018年10月21日飞行中,任务系统AIMMS-20探测环境气压从地面1 015.7 hPa一直变化到飞行高度最高时554.0 hPa。与AIMMS-20比较,CIP探测的环境气压最大偏高6.0 hPa,最小偏低8.6 hPa(不到1%),平均偏高1.9 hPa。二者对环境气压观测总体上非常接近。
图 9显示的是飞机平台ADS与任务系统AIMMS-20、CIP观测的环境气温及相对湿度变化。三者对气温观测变化的趋势相同,且观测值也较接近。飞机平台ADS观测气温总体上偏低于AIMMS-20,该案例中,ADS观测环境气温比AIMMS-20平均偏低1.4℃,最大偏高1.4℃,最小偏低6.0℃,低幅最大时出现在12:15:38,观测显示该时段大气中存在明显逆温,自12:15:33—12:15:38,5 s短暂飞行环境气温从16.5℃上升到了21.5℃,增幅达5℃。任务系统CIP探测的环境气温比AIMMS-20平均偏低0.6℃,最大偏高1.9℃,最小偏低3.7℃。从环境相对湿度变化看,AIMMS-20与CIP湿度变化趋势相近,相关系数为0.9,但二者观测湿度值存在较大差异,CIP环境湿度比AIMMS-20平均偏低8.6%,最大偏高7.2%,最小偏低19.4%。CIP与AIMMS-20探测环境湿度最大值分别为95.7%、100%,最小值分别为57.2%、64.5%。结合云区湿度特征和云降水粒子实际采样分析,飞机入云并长时间探测到云滴粒子,环境空气湿度一般不低于85%,很容易得知AIMMS-20探测环境湿度比CIP更为准确。
结合新舟60国家增雨飞机2018年10月21日飞行案例,通过对该飞机平台GPS、ADS与任务系统加装的BDS、AIMMS-20、CIP等不同观测源获取的关键飞行参数及气象要素观测对比分析,结果表明:
(1) 飞机平台GPS与任务系统BDS、AIMMS-20(GPS模块)三套定位源经纬度定位偏差较小,飞机平台GPS比任务系统BDS定位经纬度更接近AIMMS-20。BDS与AIMMS-20二者探测海拔高度偏差总体较小,飞行高度较高时更接近。飞机平台气压高度表观测海拔高度明显低于任务系统BDS和AIMMS-20测量值,且随飞行高度升高,气压高度表与后两者对海拔高度观测偏差也越大。
(2) 飞机平台ADS与任务系统AIMMS-20、CIP观测的真空速变化趋势相同,飞机平台ADS与AIMMS-20真空速观测值非常接近,ADS观测真空速比AIMMS-20平均偏快0.8 m·s-1,AIMMS-20对环境风速瞬时变化响应要比飞机平台ADS更灵敏;CIP探测真空速明显低于ADS和AIMMS-20,CIP探测真空速比AIMMS-20平均偏慢约10 m·s-1。
(3) 任务系统AIMMS-20和CIP探测环境气压变化趋势完全一致,且环境气压观测值非常接近。飞机平台ADS与任务系统AIMMS-20、CIP观测环境气温变化趋势相同且观测值较接近。ADS观测环境气温比AIMMS-20平均偏低1.4℃,偏低最大时观测显示该时段大气存在明显逆温。CIP探测环境气温比AIMMS-20平均偏低0.6℃。AIMMS-20与CIP湿度变化趋势相近,相关系数为0.9,但二者观测湿度值存在较大差异,CIP环境湿度比AIMMS-20平均偏低8.6%,AIMMS-20探测环境湿度比CIP更为准确。
新舟60增雨飞机不同机载设备观测获取的飞行参数和气象要素存在一定差异。这些差异,一方面由于飞机上不同机载设备安装位置不同所致,另一方面,大气与云结构不均匀性也会导致不同设备观测差异。通过飞机平台与飞机任务系统不同设备观测对比分析,提供的飞行参数和气象要素等基础观测数据,不仅为云物理机载探测数据合理应用提供技术指导,还能够为国家和地方人工影响天气飞机机载设备系统集成设计提供重要技术支持。
蔡兆鑫, 蔡淼, 李培仁, 等, 2019. 大陆性积云不同发展阶段宏观和微观物理特性的飞机观测研究[J]. 大气科学, 43(6): 1191-1203. Cai Z X, Cai M, Li P R, et al, 2019. Aircraft observation research on macro and microphysics characteristics of continental cumulus cloud at different development stages[J]. Chin J Atmos Sci, 43(6): 1191-1203 (in Chinese).
|
范蕙君, 马培民, 1981. 云的卫星探测和飞机探测[J]. 气象, 7(7): 42-44. Fan H J, Ma P M, 1981. Satellite and aircraft detection of clouds[J]. Meteor Mon, 7(7): 42-44 (in Chinese).
|
郭学良, 付丹红, 胡朝霞, 2013. 云降水物理与人工影响天气研究进展(2008—2012年)[J]. 大气科学, 37(2): 351-363. Guo X L, Fu D H, Hu Z X, 2013. Progress in cloud physics, precipitation, and weather modification during 2008-2012[J]. Chin J Atmos Sci, 37(2): 351-363 (in Chinese).
|
郭学良, 郑国光, 2010. 大气物理与人工影响天气[M]. 北京: 气象出版社: 526-527. Guo X L, Zheng G G, 2010. Atmospheric Physics and Weather Modification[M].
Beijing: China Meteorological Press: 526-527 (in Chinese).
|
雷恒池, 洪延超, 赵震, 等, 2008. 近年来云降水物理和人工影响天气研究进展[J]. 大气科学, 32(4): 967-974. Lei H C, Hong Y C, Zhao Z, et al, 2008. Advances in cloud and precipitation physics and weather modification in recent years[J]. Chin J Atmos Sci, 32(4): 967-974 (in Chinese).
|
廖捷, 任芝花, 2011. 我国飞机观测气温和常规高空观测气温的对比分析[J]. 气象, 37(3): 263-269. Liao J, Ren Z H, 2011. Comparative analysis on aircraft and radiosonde temperatures in China[J]. Meteor Mon, 37(3): 263-269 (in Chinese).
|
刘卫国, 苏正军, 王广河, 等, 2003. 新一代机载PMS粒子测量系统及应用[J]. 应用气象学报, 14(S1): 11-18. Liu W G, Su Z J, Wang G H, et al, 2003. Development and application of new-generation airborne particle measuring system[J]. J Appl Meteor Sci, 14(S1): 11-18 (in Chinese).
|
毛节泰, 郑国光, 2006. 对人工影响天气若干问题的探讨[J]. 应用气象学报, 17(5): 643-646. Mao J T, Zheng G G, 2006. Discussions on some weather modification issues[J]. J Appl Meteor Sci, 17(5): 643-646 (in Chinese).
|
王广河, 姚展予, 2003. 人工增雨综合技术研究[J]. 应用气象学报, 14(S1): 1-10. Wang G H, Yao Z Y, 2003. Study of comprehensive techniques for precipitation enhancement[J]. J Appl Meteor Sci, 14(S1): 1-10 (in Chinese).
|
王磊, 李成才, 赵增亮, 等, 2014. 飞机积冰云微物理特征分析及监测技术研究[J]. 气象, 49(2): 196-205. Wang L, Li C C, Zhao Z L, et al, 2014. Microphysical property analysis and detection of air icing clouds[J]. Meteor Mon, 49(2): 196-205 (in Chinese).
|
魏强, 高建春, 钱越英, 1997. 机载粒子测量系统及资料处理[J]. 气象, 23(5): 37-40. Wei Q, Gao J C, Qian Y Y, 1997. Airborne partical measuring system and data processing[J]. Meteor Mon, 23(5): 37-40 (in Chinese).
|
杨绍忠, 陶文有, 2002. 一种新型机载测温仪[J]. 气象, 28(12): 43-46. Yang S Z, Tao W Y, 2002. A new type airborne thermograph[J]. Meteor Mon, 28(12): 43-46 (in Chinese).
|
姚展予, 2006. 中国气象科学研究院人工影响天气研究进展回顾[J]. 应用气象学报, 17(6): 786-795. Yao Z Y, 2006. Review of weather modification research in Chinese Academy of Meteorological Sciences[J]. J Appl Meteor Sci, 17(6): 786-795 (in Chinese).
|
游来光, 1987. 我国云物理和人工影响天气研究工作新进展[J]. 气象科技, (1): 59-64. You L G, 1987. The new progress of cloud physics and weather modification research in China[J]. Meteor Sci Technol, (1): 59-64 (in Chinese).
|
赵增亮, 毛节泰, 王磊, 等, 2011. 一次典型层积云的飞机观测结果及与卫星资料的对比分析[J]. 气象学报, 69(3): 521-527. Zhao Z L, Mao J T, Wang L, et al, 2011. In situ aircraft observations of one typical stratocumulus cloud process compared with the satellite measurements[J]. Acta Meteor Sin, 69(3): 521-527 (in Chinese).
|
Guo X L, Zheng G G, 2009. Advances in weather modification from 1997 to 2007 in China[J]. Adv Atmos Sci, 26(2): 240-252. DOI:10.1007/s00376-009-0240-8
|
IPCC, 2007.Climate Change 2007: The physical science basis[R].Cambridge and New York, NY: Cambridge University Press.
|
Knollenberg R G, 1981.Techniques for probing cloud microstructure[M]//Hobbs P V.Clouds, Their Formation, Optical Properties, and Effects.New York: Academic Press.
|
Lawson R P, Jensen T L, 1998.Improved microphysical observations in mixed phase clouds[C]//Proceedings of the AMS Conference on Cloud Physics.Boston: AMC.
|
NRC (National Research Council), 2003. Critical Issues in Weather Modification Research[J]. Washington D.C.:National Academy Press: 51-53.
|