2023年7月17日发(作者：)

气　　　象第４７卷第１期ｏｌ．４７　Ｎｏ．１　Ｖ　　　　　　　　　　　　　　　　　　　年月Ｊａｎｕａｒ２０２１１０２１ｙ　２ＭＥＴＥＯＲＯＬＯＧＩＣＡＬＭＯＮＴＨＬＹ赵萌初，李清泉，沈新勇，等，］气象，：２０２１．１９８４—２０１７年ＩＳＣＣＰＦＨ辐射资料在青藏高原地区的适用性评估［Ｊ．４７（１）１１２３．，ＬｉＱＱ，ＳｈｅｎＸＹ，ｅｔａｌ２０２１．ＡｌｉｃａｂｉｌｉｔｓｓｅｓｓｍｅｎｔｏｆＩＳＣＣＰＦＨｒａｄｉａｔｉｏｎｐｒｏｄｕｃｔｏｖｅｒｔｈｅＴｉｂｅｔａｎＰｌａｔｅａｕｄｕｒＺｈａｏＭＣ，ｐｐｙａ［］，（）：（）Ｊ．ＭｅｔｅｏｒＭｏｎ４７１１１２３ｉｎＣｈｉｎｅｓｅ．ｉｎ９８４－２０１７ｇ１—１９８４２０１７年犐犛犆犆犘犉犎辐射资料在青藏高原地区的适用性评估２１，２１，３４赵萌初１，　李清泉　沈新勇　王庆元１南京信息工程大学气象灾害教育部重点实验室／气候与环境变化国际合作联合实验室／气象灾害预报预警与评估协同创新中心，南京２１００４４中国气象局气候研究开放实验室，北京１０００８１２国家气候中心，珠海），珠海５１９０８２３南方海洋科学与工程广东省实验室（天津３４天津市气象局，０００７４提　要：长序列卫星辐射资料的缺乏一直是制约青藏高原（以下简称高原）辐射长期变化研究的重要原因之一。对国际上最新提供的１以下简称Ｆ长序列卫星辐射资料中的大气顶逸出长波辐射（、到达地面短波辐射９８４—２０１７年ＩＳＣＣＰＦＨ（Ｈ）ＯＬＲ）（、地面向上长波辐射（、到达地面长波辐射（进行分析，评估了ＦＳＷＤ）ＬＷＵ）ＬＷＤ）Ｈ辐射资料在全天气条件下的青藏高原地区的适用性。结果表明：与观测相比，其中ＯＦＨ资料的４种辐射通量气候平均值误差均小于５％，ＬＲ和ＳＷＤ的偏差较小，ＬＷＵ的偏差最大。ＦＨ资料能正确反映高原各辐射通量的冬季增强趋势，ＯＬＲ和ＬＷＤ在各季节的长期变化趋势均与观测一致，在高原地区，ＬＷＵ则呈现虚假的减弱趋势。总体来说，ＦＨ资料的地面短波辐射通量比长波辐射通量适用性好。进一步对长波辐射偏差原因分析显示，气温偏差会增强ＬＷＤ的气候态和长期趋势，而地温偏差对ＬＷＵ的作用与之相反。辐射云和水汽的差异导致最终Ｆ模型、Ｈ资料中的ＬＷＤ气候态和长期趋势较观测略偏弱，ＦＨ资料的计算方案在一定程度上修正了地温偏差造成的ＬＷＵ偏弱。研究结果将为使用ＦＨ辐射资料提供参考依据。关键词：青藏高原，卫星资料，辐射通量，长期变化ＩＳＣＣＰＦＨ，中图分类号：：／犗犐Ｐ４０５　　　　　　文献标志码：Ａ　　　　　　犇１０．７５１９．ｉｓｓｎ．１００００５２６．２０２１．０１．００２ｊＡｌｉｃａｂｉｌｉｔｓｓｅｓｓｍｅｎｔｏｆＩＳＣＣＰＦＨＲａｄｉａｔｉｏｎＰｒｏｄｕｃｔｐｐｙＡｏｖｅｒｔｈｅＴｉｂｅｔａｎＰｌａｔｅａｕＤｕｒｉｎ９８４－２０１７ｇ１１２１２１３４ＺＨＡＯＭｅｎｃｈｕＩＱｉｎｕａｎＨＥＮＸｉｎｏｎｉｎｕａｎ　Ｌ　Ｓｇｇｑｙｇ　ＷＡＮＧＱｇｙ，，，，Ｍ／１ＫｅａｂｏｒａｔｏｒｆＭｅｔｅｏｒｏｌｏｉｃａｌＤｉｓａｓｔｅｒｉｎｉｓｔｒｆＥｄｕｃａｔｉｏｎＪｏｉｎｔＩｎｔｅｒｎａｔｉｏｎａｌＲｅｓｅａｒｃｈＬａｂｏｒａｔｏｒｆｙＬｙｏｇｙｏｙｏ／，ＣｏｌｌａｂｏｒａｔｉｖｅＩｎｎｏｖａｔｉｏｎＣｅｎｔｅｒｏｎＦｏｒｅｃａｓｔａｎｄＥｖａｌｕａｔｉｏｎｏｆＭｅｔｅｏｒｏｌｏｉｃａｌＤｉｓａｓｔｅｒｓＣｌｉｍａｔｅａｎｄＥｎｖｉｒｏｎｍｅｎｔＣｈａｎｅｇｇ，ＮＮａｎｉｎｎｉｖｅｒｓｉｔｆＩｎｆｏｒｍａｔｉｏｎＳｃｉｅｎｃｅａｎｄＴｅｃｈｎｏｌｏａｎｉｎ１００４４ｊｇＵｙｏｇｙｊｇ２，Ｎ，ａｔｉｏｎａｌＣｌｉｍａｔｅＣｅｎｔｒｅＣＭＡ，Ｂｅｉｉｎ０００８１２ＬａｂｏｒａｔｏｒｏｒＣｌｉｍａｔｅＳｔｕｄｉｅｓｊｇ１ｙｆ），３ＳｏｕｔｈｅｒｎＭａｒｉｎｅＳｃｉｅｎｃｅａｎｄＥｎｉｎｅｅｒｉｎｕａｎｄｏｎａｂｏｒａｔｏｒＺｈｕｈａｉＺｈｕｈａｉ５１９０８２ｇｇＧｇｇＬｙ（，Ｔｉａｎｉｎ３０００７４４ＴｉａｎｉｎＭｅｔｅｏｒｏｌｏｉｃａｌＳｅｒｖｉｃｅｊｊｇ：犃犫狊狋狉犪犮狋ＴｈｅｓｈｏｒｔａｅｏｆｌｏｎｔｅｒｍｓａｔｅｌｌｉｔｅｒａｄｉａｔｉｏｎｐｒｏｄｕｃｔｓｈａｓｒｅｓｔｒｉｃｔｅｄｔｈｅｓｔｕｄｆｒａｄｉａｔｉｏｎｓａｔｉｏｇｇｙｏｐｔｅｍｏｒａｌｖａｒｉａｔｉｏｎｏｖｅｒｔｈｅＴｉｂｅｔａｎＰｌａｔｅａｕ（ＴＰ）．Ｔｈｉｓｐａｅｒｐｒｅｓｅｎｔｓｅｖａｌｕａｔｉｏｎｏｆｔｈｅｓｔａｔｅｏｆｔｈｅａｒｔｐｐ）、国家“第二次青藏高原综合科学考察研究”（、国家自然科学基金重大ＸＤＡ２０１００３０４２０１９ＱＺＫＫ０２０８）　中国科学院战略性先导科技专项（项目（）、国家重点研发计划（）和国家自然科学基金项目（）共同资助４１７９０４７１２０１６ＹＦＡ０６０２２００４１９７５０５４、４１９３０９６７２０２０年２月３日收稿；０２０年１０月２０日收修定稿　２赵萌初，主要从事气候变化的研究．：第一作者：Ｅｍａｉｌ７３６８９０２２５＠ｑ．ｃｏｍｑ通讯作者：沈新勇，主要从事中尺度气象学和气候动力学研究．：Ｅｍａｉｌｓｈｅｎｘｕｉｓｔ．ｅｄｕ．ｃｎ＠ｎｙ７卷　　　　　　　　　　　　　　　　　　　　　气　　象　　　　　　　　　　　　　　　　　　　　第４１２ｓａｔｅｌｌｉｔｅｐｒｏｄｕｃｔｓｃａｌｌｅｄＩＳＣＣＰＦＨ（ｈｅｒｅｉｎａｆｔｅｒｔｏｂｅｒｅｆｅｒｅｄａｓＦＨ）ｏｖｅｒｔｈｅＴＰ．ＴｈｅａｌｌｓｋＬＲ，ｙＯ，ｄｏｗｎｗａｒｄｓｈｏｒｔｗａｖｅｒａｄｉａｔｉｏｎ（ＳＷＤ）ｕｗａｒｄｌｏｎｗａｖｅｒａｄｉａｔｉｏｎ（ＬＷＵ）ａｎｄｄｏｗｎｗａｒｄｌｏｎｗａｖｅｒａｄｉｐｇｇａｔｉｏｎ（ＬＷＤ）ｄｕｒｉｎ９８４－２０１７ｗｅｒｅｃｏｍａｒｅｄｗｉｔｈｏｂｓｅｒｖａｔｉｏｎ．Ｔｈｅｒｅｓｕｌｔｓｉｎｄｉｃａｔｅｔｈａｔｔｈｅｄｅｖｉａｔｉｏｎｓｇ１ｐｏｆｒａｄｉａｔｉｏｎｆｌｕｘｃｌｉｍａｔｏｌｏｒｅｗｉｔｈｉｎ５％，ｓｌｉｈｔｅｒｒｏｒｅｘｉｓｔｓｉｎＯＬＲａｎｄＳＷＤｗｈｉｌｅｌａｒｅｂｉａｓｉｓｆｏｕｎｄｇｙａｇｇ，ｉｎＬＷＵ．ＦＨｒａｄｉａｔｉｏｎｆｌｕｘｅｓａｒｅｐｒｏｅｒｌｉｎｃｒｅａｓｉｎｉｎｗｉｎｔｅｒｔｈｅｔｒｅｎｄｓｏｆＯＬＲａｎｄＬＷＤａｒｅｃｏｎｓｉｓｔｅｎｔｐｙｇｗｉｔｈｏｂｓｅｒｖａｔｉｏｎｓｉｎａｌｌｆｏｕｒｓｅａｓｏｎｓｂｕｔＬＷＵｗｅａｋｅｎｓｆａｌｓｅｌ．ＴｈｅＦＨｓｈｏｒｔｗａｖｅｒａｄｉａｔｉｏｎｉｓｇｅｎｅｒａｌｌｙｙｂｅｔｔｅｒｔｈａｎｌｏｎｗａｖｅｏｖｅｒｔｈｅＴＰ．Ｔｈｅｌｏｎｗａｖｅｒａｄｉａｔｉｏｎｅｒｒｏｒｓｃａｎｂｅａｔｔｒｉｂｕｔｅｄｔｏｔｈａｔａｉｒｔｅｍｅｒａｔｕｒｅｇｇｐ，ａｎｄｔｈｅｉｎｆｌｕｅｎｃｅｏｆｓｕｒｆａｃｅｔｅｍｅｒａｔｕｒｅｏｎＬＷＵｉｓｅｒｒｏｒｗｏｕｌｄｏｖｅｒｅｓｔｉｍａｔｅＬＷＤｃｌｉｍａｔｏｌｏｎｄｔｒｅｎｄｐｇｙａ，ｏｏｓｉｔｅ．ＴｈｅｒａｄｉａｔｉｏｎｍｏｄｅｌｃｌｏｕｄａｎｄｍｏｉｓｔｕｒｅｄｉｆｆｅｒｅｎｃｅｓｍａｅｓｕｌｔｉｎａｒｅｃｉａｂｌｅｕｎｄｅｒｅｓｔｉｍａｔｉｏｎｏｆｐｐｙｒｐｐＬＷＤ，ｗｈｉｌｅｔｈｅｃａｌｃｕｌａｔｉｏｎｐｒｏｃｅｓｓｏｆＦＨｐａｒｔｌｏｒｒｅｃｔｓｔｈｅｕｎｄｅｒｅｓｔｉｍａｔｉｏｎｏｆＬＷＵ．ＴｈｉｓｓｔｕｄｏｕｌｄｙｃｙｗｒｏｖｉｄｅｓｏｍｅｒｅｆｅｒｅｎｃｅｓｆｏｒｔｈｅｆｕｔｕｒｅｕｓｅｏｆＦＨｐｒｏｄｕｃｔｓ．ｐ：，，犓犲狅狉犱狊ＴｉｂｅｔａｎＰｌａｔｅａｕ，ｓａｔｅｌｌｉｔｅｐｒｏｄｕｃｔＩＳＣＣＰＦＨ，ｒａｄｉａｔｉｏｎｆｌｕｘｌｏｎｔｅｒｍｖａｒｉａｔｉｏｎｇ狔狑（，；，；；ＲｏｓｓｏｗａｎｄＳｃｈｉｆｆｅｒ１９９１Ｚｈａｎｔａｌ２００４２００６ｇｅ，）。目前这几套主流卫星辐射资料在高Ｃｏｘｅｔａｌ２００６原地区的适用性均弱于平原地区，ＧＥＷＥＸＳＲＢ资料青藏高原平均海拔高度超过４ｋ占我国陆地ｍ，面积的四分之一，是世界上海拔最高的高原，被称为“世界屋脊”“第三极”。高原位于东西风带间，海拔高度变化快，水汽状况受季风影响强，高原动力和热美洲等地气候造成影响（段安力作用将对中国东北、民等，师锐和何光碧，张镇宏等，，２０１８；２０１８；２０１９）辐射加热在高原热力作用中起至关重要的作用（叶，，笃正和高由禧，１９７９；Ｗｕｅｔａｌ１９９７；Ｗａｎｔａｌｇｅ）。对高原辐射的研究将有利于进一步认识高２００８原气候变化。为了定量研究高原辐射，早期的研究主要通过高原试验直接测量地面辐射，但这种方法存在几个缺点，一是站点数量少，难以全面分析高原辐射空间分布，且无法获取大气顶辐射特征；二是持续时间短，不能对辐射量长期变化进行分析（季国良等，马伟强等，）。此外周允华（）、翁笃鸣１９９５；２００４１９８４）和Ｗ（）也利用常规气象要素等（１９８８ａｎｔａｌ２００９ｇｅ和经验公式计算辐射量，进而分析辐射长期变化特征。与前一种方法类似，该方法也存在区域的局限目前高原观测站点大部分位于高原中东部，高原性，西部很难获得长期观测资料。卫星观测的发展，为研究高原辐射变化提供了新研究途径（蒋兴文和李跃清，王艺等，２０１０；２０１６；于涵等，王美蓉等，。国际卫星云项目２０１８；２０１９）（，ＩｎｔｅｒｎａｔｉｏｎａｌＳａｔｅｌｌｉｔｅＣｌｏｕｄＣｌｉｍａｔｏｌｏｒｏｅｃｔｇｙＰｊ自１ＩＳＣＣＰ）９８３年起收集并分析卫星观测辐射数据进而衍生出并推断全球范围云的性质和分布，ＩＳＣＣＰＦ系列、ＧＥＷＥＸＳＲＢ等卫星反演辐射资料在计算过程中对海拔高度影响的忽视，导致到达地面短波辐射相对于观测资料的低估；而ＩＳＣＣＰＦＤ（以下简称Ｆ资料使用的ＴＤ）ＯＶＳ数据在高原地区大气温度的误差，造成ＦＤ资料向下长波辐射的低估（，；，。Ｙａｎｔａｌ２００６ａ２０１０；Ｗａｎｔａｌ２０１２）ｇｅｇｅ云和地球辐射能量系统（２０００年开始，Ｃｌｏｕｄａｎｄ项目ｔｈｅＥａｒｔｈ’ｓＲａｄｉａｎｔＥｎｅｒｓｔｅｍ，ＣＥＲＥＳ）ｇｙＳｙ／／通过搭载在ＴｅｒｒａＡｕａＮＯＡＡ２０卫星上的六种ｑ仪器直接观测大气顶的地球反射短波辐射和射出长波辐射，并通过其他项目（如中分辨率成像光谱仪，，ＭＯＤｅｒａｔｅＲｅｓｏｌｕｔｉｏｎＩｍａｉｎｅｃｔｒｏｒａｄｉｏｍｅｔｅｒＭＯｇｇＳｐ）观测的云量、地面温度以及地表比辐射率等气象ＤＩＳ要素反演地表辐射通量（，Ｗｉｅｌｉｃｋｉｅｔａｌ１９９６；Ｋａｔｏ，）。Ｃｅｔａｌ２０１８ＥＲＥＳ资料相比于前面两套资料在数据观测和处理上有很大进步，但在表征地表净辐射时仍存在较大误差，其误差来源于计算地面辐射和观测站点相比，其到达地时气象数据版本的更替；面短波辐射产品在全球范围总体决定系数为０．６９，但在高原地区仅为０．（，；３５Ｇｕｉｅｔａｌ２０１０ａ２０１０ｂ；Ｊｉａ，）。这说明在利用卫星辐射资料分析高原２０１８ｅｔａｌ地区辐射特征前，需要同时对资料空间分布和时间演变的可靠性进行评估。以下２０１９年末，ＮＡＳＡ完成了对ＩＳＣＣＰＦＨ（资料的延长处理，使其成为目前时间跨度简称ＦＨ）最长的卫星辐射资料，相比于２００９年结束的ＧＥＷＥＸＳＲＢ资料和２０００年开始的ＣＥＲＥＳ资料，３４年的时间长度为研究大气顶及地面辐射特征的长期变化提供了新的方案。ＦＨ资料延长到２０１７年的同引　言１９８４—２０１７年ＩＳＣＣＰＦＨ辐射资料在青藏高原地区的适用性评估　　　　　　　　　第１期　　　　　　　赵萌初等：１３时，采用新的辐射传输模型，改进大气吸收短波辐射的参数化方案，提高了模式垂直分辨率（４３层），改进了与长波辐射关系密切的水汽过程，更新了输入模型的气象要素数据源，并大幅提高输出数据分辨率（Ｓｃｈｍｉｄｔｅｔａｌ，２００６；Ｙｏｕｎｇｅｔａｌ，２０１８）。目前，ＦＨ资料算法和输入数据的改进能否改善其在高原地区的适用性尚未得知。本次研究将１９８４—２０１７年ＦＨ资料与基于观测数据计算的长时间序列辐射量进行对比，从气候态、长期趋势和时间演变三个方面，对ＦＨ资料的长、短波辐射通量在高原地区的适用性进行评估，并进一步探讨误差来源，为未来采用ＦＨ资料开展各种研究和应用提供参考依据。　资料和方法１．１　资　料观测资料采用中国气象局提供的１９８４—２０１７年高原地区常规气象站逐月观测数据，参照Ｄｕａｎｅｔａｌ（２０１１）方法选取７０个数据可靠的站点，站点分布如图１所示。用到的气象要素有气温、地温、水汽压、日照百分率以及采用太阳辐射综合观测模型（Ｈｙｂｉｒｄ）研制的到达地面短波辐射（Ｙａｎｇｅｔａｌ，２００１）。本文采用了１９８４—２０１７年ＦＨ卫星辐射资料的全天气状况下逸出长波辐射（ＯＬＲ）、到达地面的短波辐射（ＳＷＤ）、到达地面的长波辐射（ＬＷＤ）、地面向上长波辐射（ＬＷＵ），以及ＩＳＣＣＰＨ系列气温（ＳＡＴ）、地温（ＳＴ）和地表发射率数据。ＦＨ资料采图１　高原地区７０个气象观测站分布（填色为高原地形高度，白色点为站点分布）Ｆｉｇ．１　Ｄｉｓｔｒｉｂｕｔｉｏｎｏｆ７０ｍｅｔｅｏｒｏｌｏｇｉｃａｌｏｂｓｅｒｖａｔｉｏｎｓｔａｔｉｏｎｓｏｖｅｒｔｈｅＴＰ（ｃｏｌｏｒｅｄ：ｈｅｉｇｈｔｏｆｔｅｒｒａｉｎｅｌｅｖａｔｉｏｎ，ｗｈｉｔｅｐｏｉｎｔ：ｌｏｃａｔｉｏｎｏｆｓｔａｔｉｏｎ）用ＩＳＣＣＰＨ系列１０ｋｍ分辨率辐照度、云和地面气象要素数据（ＦＤ资料分辨率为３０ｋｍ）、ＭＡＣｖ２气溶胶数据集和ｎｎＨＩＲＳ大气温湿数据集（ＦＤ资料分别为ＭＡＣｖ１和ＴＯＶＳ），基于ＮＡＳＡ的ＧＣＭ模式Ｅ版本辐射传输模型计算辐射通量（ＦＤ资料为Ｄ版本），水平分辨率为１°×１°（ＦＤ资料为２．５°×２．５°）。此外，ＮＯＡＡ的ＯＬＲ数据被用来检验ＦＨ资料的ＯＬＲ准确性。１．２　方　法为了检验评估ＦＨ卫星辐射资料，本文首先采用Ｃｒａｗｆｏｒｄ和Ｄｕｃｈｏｎ模型（简称ＣＤ９９模型），根据高原台站观测数据计算ＬＷＤ（ＣｒａｗｆｏｒｄａｎｄＤｕｃｈｏｎ，１９９９）。这个计算方法被广泛使用，且已有研究证明ＣＤ９９模型适用于包括高原在内的高海拔地区（Ｙａｎｇｅｔａｌ，２００６ｂ；Ｌｈｏｍｍｅｅｔａｌ，２００７），计算方法如下：犔犠犇＝εασ犜４ａ（１）εα＝犮犾犳＋（１－犮犾犳）［（犲０．１４２８６１．２４犜（ａ）］２）犮犾犳＝０．７２２３－０．８６３６狀／犖＋０．１４１３（狀／犖）２（３）式中：σ＝５．６７×１０－８Ｗ·ｍ－２·Ｋ－４，为斯蒂芬玻尔兹曼常数，εα为大气发射率，犲为水汽压（单位：ｈＰａ），犜ａ为ＳＡＴ（单位：Ｋ），犮犾犳代表云反射作用，狀为实际日照时数（单位：ｈ），犖为最大可能日照时数（单位：ｈ），狀／犖代表日照百分率。ＬＷＵ包含了地面发射长波辐射以及小部分地面反射的ＬＷＤ，采用以下公式计算：犔犠犝＝εｓσ犜４ｓ＋（１－εｓ）犔犠犇（４）式中：犜ｓ代表ＳＴ（单位：Ｋ），εｓ为地表发射率。已有研究表明，高原中东部εｓ大约为０．９６（Ｇａｏｅｔａｌ，１９９８），因此本文在计算台站观测ＬＷＵ时，εｓ取固定值０．９６。分别采用一元线性回归和皮尔逊相关系数犚计算比较资料线性趋势和时间演变的相关性，平均偏差（ＭＢＥ）和决定系数犚２被用来定量ＦＨ资料的偏差大小。为统一比较卫星和站点观测数据，采用双线性插值将ＦＨ和ＮＯＡＡＯＬＲ资料插值到高原地区的７０个观测站点。２　ＦＨ资料气候态和长期趋势的检验高原气候具有一定的区域差异，该现象在观测站点集中的高原中东部（ＣＥＴＰ）也存在。在研究１７卷　　　　　　　　　　　　　　　　　　　　　气　　象　　　　　　　　　　　　　　　　　　　　第４１４通常将其分为ＣＥＴＰ的辐射及大气热源等特征时，南北两部分进行讨论（谢志昂和段安民，除多２０１７；等，于涵等，）。通过对四种观测辐射通量２０１８；２０１８ＥＯＦ分析也发现，ＥＯＦ第一模态表现为全区一致，第二模态以３图略）。那３°Ｎ为界呈南北反相变化（么，ＦＨ资料在高原南部和北部的适用性是否存在区域差异？为验证本文中以高原全部站点平均的辐我们还以射通量特征进行ＦＨ资料评估的合理性，和北部３３°Ｎ为分界分别计算了高原南部（ＳＴＰ）（站点平均的辐射通量气候态和长期趋势。ＮＴＰ）高原地区多年平均的年平均和季节平均辐射通量以及温度如图２所示，ＦＨ资料中的ＯＬＲ在不同季节均弱于观测，其中夏季偏差较小（约为２Ｗ·ｍ－２），春季偏差最大（约为８Ｗ·ｍ－２）。ＦＨ中的ＳＷＤ在冬季与观测较为接近，偏差在５Ｗ·ｍ－２之内，但夏季ＦＨ和观测之间的误差超过１５Ｗ·ｍ－２。ＳＷＤ虽然和太阳辐射关系密切，但受大气状况影响也十分明观测中ＳＷＤ年循环并非稳定的先增强、后减显，而是在３月和５月大幅增强，弱，７月强度甚至强于并于５月达６月；ＦＨ中的ＳＷＤ从１月开始逐渐增加，图２　１、高原北部（和高原整体（年平均（）和冬季（）、９８４—２０１７年高原南部（ＳＴＰ）ＮＴＰ）ＴＰ）ａｂ春季（）、夏季（）、秋季（）平均的Ｆ）资料辐射通量以及ＳｃｄｅＨ与观测（ＯＢＳＴ和ＳＡＴＦｉ．２　ＡｎｎｕａｌａｎｄｓｅａｓｏｎａｌｍｅａｎｓｏｆｔｈｅＦＨａｎｄｏｂｓｅｒｖｅｄｒａｄｉａｔｉｏｎｆｌｕｘｅｓＳＴａｎｄｇＳＡＴｏｖｅｒｔｈｅＳＴＰ，ＮＴＰａｎｄＴＰｄｕｒｉｎ９８４－２０１７ｇ１（），（）ｗ，（），（），（）ａａｎｎｕａｌｂｉｎｔｅｒｃｓｒｉｎｄｓｕｍｍｅｒｅａｕｔｕｍｎｐｇ１９８４—２０１７年ＩＳＣＣＰＦＨ辐射资料在青藏高原地区的适用性评估　　　　　　　　　第１期　　　　　　　赵萌初等：１５到最大，此后逐渐减弱（图略），这可能与ＦＨ资料计算短波辐射过程中对云量的逐月变化的估计不足有关。ＦＨ中的长波辐射通量及其对应的气温和地温与观测资料有一定差别（图２）：ＬＷＤ存在较大季节差异，夏季辐射强度可达到冬季的１．５倍，ＦＨ资料且误差也随着辐射强度增强低估了高原地区ＬＷＤ，而增大，夏季误差约为１两套资料ＬＷＵ５Ｗ·ｍ－２；的季节变化与ＬＷＤ相似，但ＦＨ资料在冬季存在在夏季存在２６Ｗ·ｍ－２的高估，０Ｗ·ｍ－２的低估。图２中可以发现ＦＨ资料ＬＷＵ偏差与地温偏差十分相似，均为冬季高估而夏季低估，尤其是夏季ＦＨ资料低估了大约１１Ｋ的地温。理论上ＬＷＤ与气温关系密切，但是ＦＨ资料逐月气候态气温均高于观测，ＬＷＤ气候态低于观测。从四个季节看，ＦＨ的Ｏ最大偏差出现ＬＲ在不同季节的偏差都较小，，其他三种辐射通量的夏在春季（－６．６５Ｗ·ｍ－２）季偏差最大。除夏季外其他季节ＦＨ的各辐射通量并以冬、春季误差最ＭＢＥ均不超过１０Ｗ·ｍ－２，小。总体来说，ＦＨ资料的四种辐射通量ＭＢＥ相对虽然Ｆ于其气候态基本上都在５％以内，Ｈ资料最大地温偏差超过１由于计算ＬＷＵ时采用了开１Ｋ，尔文温标，温度的最大ＭＢＥ百分比仍小于４％。因此，ＦＨ资料辐射通量气候态的误差范围较小。Ｈ资料是目前唯一从１９８４年开始并延续到　　Ｆ现在的卫星辐射资料，其重要优势之一就是便于研究高原辐射长期变化趋势，因此我们对ＦＨ资料高原地区长期趋势进行了评估。图３比较了高原地区平均的ＦＨ与观测资料不同季节线性趋势。ＦＨ资料图３　同图２，但为Ｆ）资料辐射通量线性趋势以及ＳＨ和观测（ＯＢＳＴ和ＳＡＴ的线性趋势，Ｆｉ．３　ＳａｍｅａｓＦｉ．２，ｂｕｔｆｏｒｌｉｎｅａｒｔｒｅｎｄｓｏｆｔｈｅａｎｎｕａｌａｎｄｓｅａｓｏｎａｌｍｅａｎｒａｄｉａｔｉｏｎｆｌｕｘｅｓＳＴａｎｄＳＡＴｇｇ７卷　　　　　　　　　　　　　　　　　　　　　气　　象　　　　　　　　　　　　　　　　　　　　第４１６能正确反映冬季高原辐射通量和气温的增强，ＯＬＲ和长波辐射的冬季趋势低于观测约０．６～１．６Ｗ·－１），ｍ－２·（１０ａＳＷＤ和气温较观测偏强。其他季节ＦＨ资料的趋势与观测有较大差别：ＦＨ的ＯＬＲ－１增加趋势仅有冬季显著［）］，１．６Ｗ·ｍ－２·（１０ａ大差异，但观测量的增长明显。ＦＨ资料趋势微弱，观测的春季ＳＷＤ为南部减弱北部增强，而夏、秋季与春季相反，在高原北部增强南部减弱。ＦＨ资料正确反映了秋、冬季高原地区ＳＷＤ长期趋势的空间分布，但未表现出春、夏季高原中部ＳＷＤ的减弱但观测的ＯＬＲ在秋、冬季均显著增加，趋势分别为３．１０和１．５３Ｗ·ｍ－２·（１０ａ）－１；观测的ＳＷＤ在夏季显著减弱，冬季增强，其他季节呈不显著的减弱趋势，而ＦＨ的ＳＷＤ在各个季节明显增强，春季趋势达到５．０１Ｗ·ｍ－２·（１０ａ）－１；基于观测的地面长波辐射在不同季节均显著增强，但ＦＨ的ＬＷＤ增加趋势不显著，ＬＷＵ仅在冬季显著增强，其他季节存在不同程度的减弱趋势。ＦＨ中的ＬＷＵ不变或减弱的趋势显然是存在明显误差的，研究表明在全球加速变暖背景下高原也有显著的增暖趋势（丁一汇和张莉，２００８；段安民等，２０１６），这在图３观测气温和地温的显著增强中也有所体现。ＦＨ资料使用的ｎｎＨＩＲＳ气温与观测表现为一致增暖，因此其对ＬＷＤ长期趋势的表现较好，但ＦＨ的地温在高原地区存在明显的趋势误差，年平均、夏季和秋季线性趋势减弱，且夏季减弱趋势达到－０．４Ｋ·（１０ａ）－１，这可能导致了ＦＨ中的ＬＷＵ趋势计算错误。此外，即使ＦＨ的气温增暖比观测强，ＦＨ中的ＬＷＤ增强依然明显弱于观测。　　由图２和图３可见，尽管ＦＨ资料对高原南、北部辐射通量气候态的高（低）估的量值有一定差异，但ＦＨ与观测资料的强弱关系未发生改变。部分辐射通量的长期趋势在高原南北部有不同，如观测中的秋季ＳＷＤ在ＳＴＰ减弱明显，在ＮＴＰ趋势几乎为０；ＦＨ资料的秋季ＬＷＤ在ＳＴＰ增加，在ＮＴＰ略微减弱。但与气候态类似，ＦＨ资料对同一辐射通量长期趋势的估算在高原南、北部同时偏高或偏低（图３）。上述结果表明ＦＨ资料对高原辐射通量气候态和长期趋势的高（低）估情况在高原南北部一致，但偏差的大小有不同。为进一步探究ＦＨ资料对高原不同区域辐射通量长期趋势的计算偏差，本文还给出了１９８４—２０１７年高原地区观测和ＦＨ资料辐射通量线性趋势的空间分布（图４），填色为ＦＨ资料辐射通量线性趋势，圆点内为观测辐射通量线性趋势。ＯＬＲ的观测量在高原全区呈增长趋势，ＦＨ资料能正确表现出ＯＬＲ在高原东南部增强较北部快的现象，但强度整体弱于观测。春到秋季高原中部ＯＬＲ趋势存在较趋势。ＦＨ资料的ＬＷＤ表现为高原主体增强且四周减弱，这与观测基本吻合，但高原东南角观测ＬＷＤ在各季节增速均超过４Ｗ·ｍ－２·（１０ａ）－１，而ＦＨ资料在这一区域为减弱趋势。与上文分析结果相似，ＦＨ在整个高原地区均存在对ＬＷＤ线性趋势的低估，尤其是观测中夏、秋季高原北部ＬＷＤ增长较快，但ＦＨ资料趋势微弱或有负增长。冬季ＦＨ资料的ＬＷＵ趋势与观测一致，但其他季节仅有高原南北侧３０００ｍ海拔附近有正确的线性趋势。总体来说，ＦＨ资料能正确表现各辐射通量年循环及各季节平均气候态。在１９８４—２０１７年长期趋势方面，ＦＨ中的ＳＷＤ最好，ＯＬＲ和ＬＷＤ次之，ＬＷＵ出现明显误差。ＦＨ资料的地面长波辐射误差的可能原因是其计算时采用的温度数据集存在偏差，本文在第４节对误差的来源和量值做详细分析。３　ＦＨ资料时间演变的检验本节对高原地区去线性趋势后ＦＨ资料时间演变的可靠性进行评估。图５对比了不同季节ＦＨ和观测资料的ＯＬＲ和ＳＷＤ，两套资料ＯＬＲ在各个季节均呈１１分布。春、夏季，ＦＨ的ＳＷＤ在部分站点略大于观测，总体分布较ＯＬＲ略离散，但基本也沿对角线分布。ＦＨ资料年平均ＯＬＲ和观测变化不同步，相关系数仅有０．０２；除年平均外其他季节与观测变化较为一致，相关系数均通过０．０１显著性水平检验（表１）。ＦＨ资料年平均和春季ＳＷＤ与观测相关性不显著，但其他季节相关系数均超过０．５，冬、春季ＦＨ的ＳＷＤ在２０００年附近有明显增强，但观测中没有这一现象（图略）。ＦＨ与观测的ＯＬＲ和ＳＷＤ高原地区相关系数的空间分布见图７，相关性存在明显季节差异，年平均辐射通量的相关性普遍较差，而夏、秋季几乎所有站点相关性能通过０．０１显著性水平检验。ＦＨ与观测ＯＬＲ的相关系数从高原东南部向北递减，说明ＦＨ资料对高原北部ＯＬＲ变化的表现不足；高原中部部分站点的冬、春季ＳＷＤ的相关系数未能通过显著性水平检验，但也均为正相关。１９８４—２０１７年ＩＳＣＣＰＦＨ辐射资料在青藏高原地区的适用性评估　　　　　　　　　第１期　　　　　　　赵萌初等：１７图４　１）、冬季（）、春季（）、夏季（）和秋季（）９８４—２０１７年高原地区年平均（ａｅｉｌｍ～ｐｔ～ｄ～ｈ～～ｑ－１平均观测（圆点）和Ｆ填色）辐射通量线性趋势的空间分布［单位：）］Ｈ资料（Ｗ·Ｍ－２·（１０ａ（由外到内黑色实线分别代表高原１８００和３０００ｍ海拔高度））），），Ｆｉ．４　Ｓａｔｉａｌｄｉｓｔｒｉｂｕｔｉｏｎｏｆｔｈｅｌｉｎｅａｒｔｒｅｎｄｓｏｆａｎｎｕａｌ（ａ－ｄａｎｄｗｉｎｔｅｒ（ｅ－ｈｓｒｉｎｉ－ｌｇｐｐｇ（），）ｍ））ｓｕｍｍｅｒ（ｍ－ｐａｕｔｕｍｎ（ｅａｎｒａｄｉａｔｉｏｎｆｌｕｘｅｓｏｆｏｂｓｅｒｖａｔｉｏｎ（ｄｏｔｓａｎｄＦＨ（ｃｏｌｏｒｅｄ－ｔｑ－１：）］ｏｖｅｒｔｈｅＴＰｄｕｒｉｎ９８４－２０１７［ｕｎｉｔＷ·ｍ－２·（１０ａｇ１（，）Ｂｌａｃｋｓｏｌｉｄｌｉｎｅｓｆｒｏｍｔｈｅｏｕｔｓｉｄｅｔｏｔｈｅｉｎｓｉｄｅｒｅｒｅｓｅｎｔｔｈｅ１８００ａｎｄ３０００ｍｅｌｅｖａｔｉｏｎｓｒｅｓｅｃｔｉｖｅｌｐｐｙ表１　１９８４—２０１７年高原地区平均的去趋势犉犎与观测辐射通量、犛犜和犛犃犜的相关系数，犜犪犫犾犲１　犆狅狉狉犲犾犪狋犻狅狀犮狅犲犳犳犻犮犻犲狀狋狊狅犳犱犲狋狉犲狀犱犲犱狉犪犱犻犪狋犻狅狀犳犾狌狓犲狊犛犜犪狀犱犛犃犜犫犲狋狑犲犲狀犉犎狆狉狅犱狌犮狋犪狀犱狅犫狊犲狉狏犪狋犻狅狀犗犔犚年冬春夏秋０．０２０．４５０．３７０．５７犛犠犇０．１９０．６１犔犠犇００．３５犔犠犝０．５１０．７６０．４９０．５００．４４犛犃犜０．４６０．６００．５２０．５３犛犜０．４４０．６８０．４６０．４７０．４６０．２１０．５７０．６２－０．０５０．３１０．２５０．３２０．２８１和０．０５显著性水平检验。　　　　　　注：和分别代表通过０．：ａ，ｏｔｅｎｄｎｄｉｃａｔｅｔｈａｔｔｈｅｃｏｒｒｅｌａｔｉｏｎｃｏｅｆｆｉｃｉｅｎｔｓｈａｖｅｐａｓｓｅｄｔｈｅｓｉｎｉｆｉｃａｎｃｅｔｅｓｔｓａｔ０．１ａｎｄ０．０５ｌｅｖｅｌｓｒｅｓｅｃｔｉｖｅｌ．　　　　　　ＮｉｇｐｙＤ资料的地面长波辐射在高原地区存在偏差　　Ｆ（，），这种误差在ＦＹａｎｔａｌ２０１０Ｈ资料中似乎未ｇｅ被完全改善。ＦＨ与观测资料长波辐射和温度的散点图（图６）离散更加明显，ＦＨ资料整体弱于观测，７卷　　　　　　　　　　　　　　　　　　　　　气　　象　　　　　　　　　　　　　　　　　　　　第４１８图５　高原地区年平均（，）、冬季（，）、春季（，）、夏季（，）、秋季（，）平均的观测与ＦａｆｂｃｈｄｉｅＨ资料辐射通量时间序列散点图ｇｊ（）（）ａＯＬＲ，ｆＳＷＤ～ｅ～ｊ，），），，），，），，）Ｆｉ．５　Ｓｃａｔｔｅｒｉｎｌｏｔｓｏｆａｎｎｕａｌ（ａｆａｎｄｗｉｎｔｅｒ（ｂｓｒｉｎｃｈｓｕｍｍｅｒ（ｄｉａｕｔｕｍｎ（ｅｇｇｐｐｇ（ｇｊｂｅｔｗｅｅｎｏｂｓｅｒｖｅｄａｎｄＦＨＯＬＲａｎｄＳＷＤ（），（）ａ－ｅＯＬＲｆ－ＷＤｊＳ（）ＡｂｓｃｉｓｓａａｎｄｏｒｄｉｎａｔｅｒｅｒｅｓｅｎｔＦＨｐｒｏｄｕｃｔａｎｄｏｂｓｅｒｖａｔｉｏｎ，ｒｅｓｅｃｔｉｖｅｌｐｐｙ（横坐标代表Ｆ纵坐标代表观测资料）Ｈ资料，图６　同图５，但为平均的观测与ＦＨ资料辐射通量及ＳＴ和ＳＡＴ（），（），（），（）ａＡＴｐｔＳＴ～ｅＬＷＤｆ～～ｊＬＷＵｋ～ｏＳ，（），（）））Ｆｉ．６　Ｓ－ｔａｍｅａｓＦｉ．５ｂｕｔｂｅｔｗｅｅｎｏｂｓｅｒｖｅｄａｎｄＦＨｆｏｒａ－ｅＬＷＤｆ－ｋ－ｏＳＡＴａｎｄ（ＳＴｇｐｇｊＬＷＵ，（１９８４—２０１７年ＩＳＣＣＰＦＨ辐射资料在青藏高原地区的适用性评估　　　　　　　　　第１期　　　　　　　赵萌初等：１９部分变量如年平均和夏季ＬＷＵ以及夏季地温明显偏离了１１线。从ＦＨ与观测资料相关系数的空间分布（图７）也可以看出，与ＯＬＲ相比，ＦＨ中的冬季ＦＨ中的ＬＷＤ误差更加明显。在高原南部，其他季节两套资ＬＷＤ与３～１０个站点无显著相关，料的相关系数超过－０．４；ＦＨ对高原北部ＬＷＤ时春季最差，两套资料相关性间演变的表现冬季最好、还表现为从高原东部向西部减弱的空间分布。上文分析表明ＦＨ中的ＬＷＵ气候态和趋势存在误差，但从相关性来看ＦＨ的ＬＷＵ的时间演变可靠，各相关性最差的季节与观测在高原全区均为正相关，年平均ＬＷＵ也有超过半数站点显著性通过０．０１显著性水平检验。两种资料的高原北部ＬＷＵ相关性略差，春夏季约有１０个站点相关系数不显著。Ｈ和观测资料的去趋势地面长波辐射和温　　从Ｆ度时间序列间相关系数（表１）看出：ＦＨ的ＬＷＤ与观测相关性在年平均和春季约为０，其他季节相关系数最大值也仅有０．时间演变差别较大。Ｆ３５，Ｈ中的ＬＷＤ和气温从１有高—低—高的年９８４年起，这恰好与观测反位相（图略），说明Ｆ代际变化，Ｈ资料ＬＷＤ的年代际变化存在明显偏差。两套资料各季节ＬＷＵ的相关系数均通过０．０１显著性水平检验，时间序列变化基本吻合。将资料间地面长波辐表１），我们发现射相关系数与温度相关系数对比（两套资料的地温显著相关，与地温紧密联系的但是这种温度与辐射间的对应关ＬＷＵ也显著相关。图７　同图４，但为ＦＨ与观测资料辐射通量相关系数（“代表该站点的辐射通量相关系数通过０．×”０１显著性水平检验）Ｆｉ．７　ＳａｍｅａｓＦｉ．４，ｂｕｔｆｏｒｃｏｒｒｅｌａｔｉｏｎｃｏｅｆｆｉｃｉｅｎｔｓｏｆａｎｎｕａｌａｎｄｓｅａｓｏｎａｌｍｅａｎｇｇｒａｄｉａｔｉｏｎｆｌｕｘｅｓｂｅｔｗｅｅｎＦＨａｎｄｏｂｓｅｒｖａｔｉｏｎ（）Ｔｈｅ“×”ｓｍｂｏｌｉｎｄｉｃａｔｅｓｔｈａｔｔｈｅｃｏｒｒｅｌａｔｉｏｎｃｏｅｆｆｉｃｉｅｎｔｏｆｒａｄｉａｔｉｏｎｆｌｕｘｅｓｈａｓｐａｓｓｅｄｔｈｅｓｉｎｉｆｉｃａｎｃｅｔｅｓｔａｔ０．０１ｌｅｖｅｌｙｇ７卷　　　　　　　　　　　　　　　　　　　　　气　　象　　　　　　　　　　　　　　　　　　　　第４２０系在气温和ＬＷＤ之间存在差异，虽然ＦＨ和观测的气温也有较强的相关，但资料间ＬＷＤ相关性却很差。这说明除温度外，可能有其他因素造成了ＬＷＤ和ＬＷＵ时间演变的明显偏差。上文也表明ＦＨ中的地面长波辐射气候平均值和长期趋势也有因此有必要对误差来源进行深入讨论。较大误差，资料的地温和地表发射率数据计算ＬＷＵ。Ｓ１与观测之间的差异代表了ＦＨ温度数据对长波辐射的影响；Ｓ２与观测的差别代表了不同地表发射率对ＬＷＵ的影响；ＦＨ资料和Ｓ１和Ｓ３的差异分别代表不同计算模型和云、水汽差异对ＬＷＤ和ＬＷＵ的影响。首先检验ＦＨ资料温度数据对长波辐射的影响。当使用ＦＨ资料温度数据计算长波辐射时，其他季节较观测ＬＷＤ的夏季气候态与观测一致，偏高６～１０Ｗ·ｍ－２；ＬＷＵ的冬季气候态与观测一致，其他季节偏低明显，春、秋季偏差分别为－２５．６７和－１夏季偏差达到－５８．３５Ｗ·ｍ－２，６．７Ｗ·偏差约为观测气候值的１表２）。图８ｍ－２，４％（ａ中春季是观测的两倍，在夏、秋季Ｓ１的线性趋势在冬、也高于观测，说明ＦＨ资料气温的使用会使ＬＷＤ－１线性趋势偏高１～２Ｗ·ｍ－２·（）。与原Ｆ１０ａＨ４　ＦＨ资料长波辐射误差来源分析长波辐射通量计算的准确性取决于输入数据和辐射模型的质量。在无法直接获取ＦＨ资料长波辐本文采用Ｃ射模型的情况下，Ｄ９９模型（Ｃｒａｗｆｏｒｄ，，将ＦＨ资料的温度和地表发射ａｎｄＤｕｃｈｏｎ１９９９）率数据代入Ｃ并比较计算结果、观Ｄ９９模型中计算，测和原Ｆ深入分析长波辐射Ｈ资料中的辐射通量，误差来源。误差来源可以分为三类：温度、地表发射包括计算模型及大气中云和水汽状况）。率和其他（具体计算方法如下：方案１（是仅使用ＦＳ１）Ｈ资料的温度数据，其他均使用观测数据计算长波辐射；方（）是使用ＦＨ资料的地表发射率数据，其他案２Ｓ２方案３（是使用Ｆ均使用观测数据计算ＬＷＵ；Ｓ３）Ｈ资料相似，Ｓ１计算的ＬＷＵ趋势除冬季外均与观测相反，这说明ＦＨ资料地温趋势偏差是造成ＬＷＵ趋势错误的重要原因。决定系数犚２用来表征不同表３），当温度数据集对长波辐射时间演变的影响（更换ＦＨ温度数据后，Ｓ１与观测ＬＷＤ的犚２各季节表２　１９８４—２０１７年高原地区平均的犉犎与犛１、犛２数据相对于观测的平均偏差犕犅单位：犈（犠·犿－２）犜犪犫犾犲２　犕犅犈狅犳犉犎，犛１，犛２狉犪犱犻犪狋犻狅狀犳犾狌狓狉犲犾犪狋犻狏犲狋狅狅犫狊犲狉狏犪狋犻狅狀狅狏犲狉狋犺犲犜犘：犠·犿－２）狌狀犻狋犱狌狉犻狀９８４－２０１７（犵１犔犠犇犗犅犛年冬春夏秋２６０．３５２１２．９９２５９．７５３２０．６３２５９．８３犔犠犝犛１６．４５７．８７９．５３－０．１９８．９７犉犎－１０．７６－７．５２－６．３０－１８．９５－９．４９犗犅犛３４８．６３３０６．３０３５５．１７４００．２１３４５．９３犉犎－６．７９３．８１－３．２６－２２．４３－７．５１犛１－２３．４６０．４６－２５．６７－５６．７２－１８．３５犛２４．１２３．５１４．３２５．０６４．２５ＯＢＳ代表观测辐射通量的气候平均值。　　　　　　　　　　　　　　注：：ＯｏｔｅＢＳｒｅｒｅｓｅｎｔｓｔｈｅｃｌｉｍａｔｏｌｏｉｃａｌｍｅａｎｏｆｏｂｓｅｒｖａｔｉｏｎ．　　　　　　　　　　　　　　Ｎｐｇ图８　１９８４—２０１７年高原地区年平均和季节平均的观测、ＦＨ、Ｓ１和Ｓ２的辐射通量线性趋势（）（）ａＬＷＤ，ｂＬＷＵＦｉ．８　Ｌｉｎｅａｒｔｒｅｎｄｓｏｆａｎｎｕａｌａｎｄｓｅａｓｏｎａｌｍｅａｎｒａｄｉａｔｉｏｎｆｌｕｘｅｓｏｆｇ，ｏｂｓｅｒｖａｔｉｏｎＦＨ，Ｓ１ａｎｄＳ２ｏｖｅｒｔｈｅＴＰｄｕｒｉｎ９８４－２０１７ｇ１（））ａＬＷＤ，（ｂＬＷＵ１９８４—２０１７年ＩＳＣＣＰＦＨ辐射资料在青藏高原地区的适用性评估　　　　　　　　　第１期　　　　　　　赵萌初等：２１表３　１９８４—２０１７年高原地区平均的犉犎与犛１、犛２数据相对于观测资料的决定系数犚２犜犪犫犾犲３　犇犲狋犲狉犿犻狀犪狋犻狅狀犮狅犲犳犳犻犮犻犲狀狋狊狅犳犉犎，犛１犪狀犱犛２狉犲犾犪狋犻狏犲狋狅狅犫狊犲狉狏犪狋犻狅狀狅狏犲狉狋犺犲犜犘犱狌狉犻狀犵１９８４－２０１７犔犠犇犔犠犝犉犎犛１犉犎犛１犛２年００．０６０．２８０．２３０．９９冬０．１３０．２８０．５８０．５００．９９春００．０６０．２４０．２３０．９９夏０．０８０．２３０．２８０．３１１．００秋０．０６０．２４０．２００．２２１．００均较低，ＬＷＵ的犚２较高，但最高也仅为５０％（冬季），表明温度数据集仍是造成辐射资料时间演变偏差的重要因素。　　对于地表发射率造成的误差，Ｓ２仅替换ＦＨ地表发射率，其他使用观测数据计算ＬＷＵ，与采取固定地表发射率方案计算的观测ＬＷＵ相比，ＦＨ地表辐射率方案气候态略高出３．５～５Ｗ·ｍ－２，偏差百分率低于２％。从图８ｂ中也可以发现地表发射率对趋势的影响很小，可以忽略不计。Ｓ２与观测ＬＷＵ的犚２也接近１００％，说明ＦＨ资料地表辐射率几乎不会影响ＬＷＵ的准确性（表３）。从ＦＨ与Ｓ１和Ｓ３的各季节长波辐射平均偏差发现（表４），ＦＨ资料算法和水汽的误差会导致计算的ＬＷＤ偏低１７Ｗ·ｍ－２左右；ＦＨ地温数据的使用会导致ＬＷＵ偏低２０～５０Ｗ·ｍ－２，但ＦＨ资料计算方案修正了其中６～３０Ｗ·ｍ－２的偏差。ＦＨ资料ＬＷＤ线性趋势小于Ｓ１（图８），表明算法和水汽的偏差会导致计算的ＬＷＤ趋势偏低３～６Ｗ·ｍ－２，但由于卫星资料气温的显著增强，最终ＦＨ资料中的ＬＷＤ增长速度慢于观测。此外，ＦＨ与Ｓ１和Ｓ３的决定系数犚２基本都在０．７以上（表４），说明算法和水汽的偏差未对长波辐射通量的时间演变造成明显影响。表４　１９８４—２０１７年高原地区平均的犉犎资料相对于犛１的犔犠犇、犛３的犔犠犝平均偏差和决定系数犚２犜犪犫犾犲４　犕犅犈犪狀犱犱犲狋犲狉犿犻狀犪狋犻狅狀犮狅犲犳犳犻犮犻犲狀狋狅犳犉犎狆狉狅犱狌犮狋狉犲犾犪狋犻狏犲狋狅犔犠犇狅犳犛１犪狀犱犔犠犝狅犳犛３狅狏犲狉狋犺犲犜犘犱狌狉犻狀犵１９８４－２０１７犔犠犇犔犠犝犕犅犈／（Ｗ·ｍ－２）犚２犕犅犈／（Ｗ·ｍ－２）犚２年－１７．１７０．７１１２．７５０．９３冬－１５．３９０．７８－０．１５０．９５春－１５．８３０．６９１８．３６０．８７夏－１８．６７０．５８２９．９３０．７８秋－１８．４６０．７３６．９２０．９５５　结论与讨论本文从气候态、长期趋势和时间演变三个方面对目前唯一一套３４年的长序列辐射资料ＩＳＣＣＰＦＨ在我国青藏高原地区适用性进行了全面评估，并深入探讨长波辐射通量存在明显偏差的原因，定量分析了不同误差来源对长波辐射的影响程度。主要结论如下：（１）从不同季节气候态来看，ＦＨ资料各辐射通量与观测之间的误差随辐射强度增加而增大，但误差均在５％以内。ＦＨ资料低估了１～８Ｗ·ｍ－２的高原地区ＯＬＲ，对ＳＷＤ的高估在５～２０Ｗ·ｍ－２，地面长波辐射气候态偏差较大，夏季对ＬＷＤ和ＬＷＵ的低估达到３０Ｗ·ｍ－２。ＦＨ资料未能完全反映ＳＷＤ月际变化，对高原地区大气短波辐射过程的表现不足是产生偏差的可能原因。（２）ＦＨ资料各辐射通量线性趋势的可靠性存在明显差异，同一辐射通量趋势的季节差异也较显著。ＦＨ资料能正确反映高原冬季各辐射通量的增强趋势；其ＯＬＲ和ＬＷＤ在各季节线性趋势均和观测一致；春、夏、秋季ＳＷＤ和ＬＷＵ的趋势与观测完全相反；在高原持续增暖的情况下ＬＷＵ的减弱趋势存在明显错误。ＦＨ资料三种长波辐射趋势在高原全区总体偏弱。ＦＨ资料正确反映了各季节ＯＬＲ、ＬＷＤ和秋、冬季ＳＷＤ的线性趋势空间分布，但未表现出春、夏季高原中部ＳＷＤ的减弱趋势；除冬季外其他季节ＦＨ的ＬＷＵ在高原东部减弱，与观测偏差较大。（３）对去线性趋势后ＦＨ资料辐射通量时间演变的分析发现，ＦＨ中的ＬＷＵ与观测在不同季节均有较强相关性，ＯＬＲ和ＳＷＤ也能反映冬、夏季高原地区辐射变化特征。虽ＦＨ的ＬＷＤ与观测有一致气候态和趋势，但其与观测数据的时间演变相关性较差，尤其是在年代际变化方面与观测反位相。ＦＨ资料能基本反映高原全区ＳＷＤ和ＬＷＵ的时间演变，ＦＨ与观测ＯＬＲ在高原主体的相关性强于高原四周，ＬＷＤ的误差在海拔快速增加的高原南坡最为明显。（４）温度偏差不是造成ＦＨ资料地面长波辐射误差的唯一原因。ＦＨ气温数据会造成ＬＷＤ的气候平均值相对于观测偏高６～１０Ｗ·ｍ－２、线性趋势偏高１～２Ｗ·ｍ－２·（１０ａ）－１；算法和水汽的作７卷　　　　　　　　　　　　　　　　　　　　　气　　象　　　　　　　　　　　　　　　　　　　　第４２２用将导致ＬＷＤ相对于观测气候态偏低１５～１９Ｗ·ｍ－２，线性趋势偏低２～４Ｗ·ｍ－２·（１０ａ）－１，此外算法和水汽的误差对ＬＷＤ时间演变的影响大于气温。ＦＨ地表发射率相对于固定地表发射率方案会高估ＬＷＵ３～５Ｗ·ｍ－２，但不会影响ＬＷＵ的长期趋势和时间演变可靠性，ＦＨ地温数据的使用会导致ＬＷＵ偏低２０～５０Ｗ·ｍ－２，但ＦＨ资料计算方案对偏差有一定修正作用。目前对高原辐射长期变化特征的研究仍然欠缺，各辐射资料之间的不一致和长序列资料的匮乏是阻碍研究的重要原因，随着ＩＳＣＣＰＨ系列长序列卫星资料的出现，一定程度上有助于相关研究的推进。本文的分析表明多方面原因造成了ＦＨ资料在高原地区的偏差，因此在使用该资料时选择合适季节或变量会得到更加可靠的结果。除此之外，仍有一些问题尚未明确，例如除了本文评估的四种辐射通量外，大气顶逸出短波辐射和地面反射短波辐射这两个变量很难获得长期观测资料进行评估，这将导致ＦＨ资料大气顶和地面辐射平衡存在不确定性；在明确误差来源后，是否能对资料偏差订正，改善其在高原地区适用性也有待尝试。参考文献除多，洛桑曲珍，林志强，等，２０１８．近３０年青藏高原雪深时空变化特征分析［Ｊ］．气象，４４（２）：２３３２４３．ＣｈｕＤ，ＬｕｏｓａｎｇＱＺ，ＬｉｎＺＱ，ｅｔａｌ，２０１８．ＳｐａｔｉｏｔｅｍｐｏｒａｌｖａｒｉａｔｉｏｎｏｆｓｎｏｗｄｅｐｔｈｏｎＴｉｂｅｔａｎＰｌａｔｅａｕｏｖｅｒｔｈｅｌａｓｔ３０ｙｅａｒｓ［Ｊ］．ＭｅｔｅｏｒＭｏｎ，４４（２）：２３３２４３（ｉｎＣｈｉｎｅｓｅ）．丁一汇，张莉，２００８．青藏高原与中国其他地区气候突变时间的比较［Ｊ］．大气科学，３２（４）：７９４８０５．ＤｉｎｇＹＨ，ＺｈａｎｇＬ，２００８．ＩｎｔｅｒｃｏｍｐａｒｉｓｏｎｏｆｔｈｅｔｉｍｅｆｏｒｃｌｉｍａｔｅａｂｒｕｐｔｃｈａｎｇｅｂｅｔｗｅｅｎｔｈｅＴｉｂｅｔａｎＰｌａｔｅａｕａｎｄｏｔｈｅｒｒｅｇｉｏｎｓｉｎＣｈｉｎａ［Ｊ］．ＣｈｉｎＪＡｔｍｏｓＳｃｉ，３２（４）：７９４８０５（ｉｎＣｈｉｎｅｓｅ）．段安民，肖志祥，王子谦，２０１８．青藏高原冬春积雪和地表热源影响亚洲夏季风的研究进展［Ｊ］．大气科学，４２（４）：７５５７６６．ＤｕａｎＡＭ，ＸｉａｏＺＸ，ＷａｎｇＺＱ，２０１８．ＩｍｐａｃｔｓｏｆｔｈｅＴｉｂｅｔａｎＰｌａｔｅａｕｗｉｎｔｅｒ／ｓｐｒｉｎｇｓｎｏｗｄｅｐｔｈａｎｄｓｕｒｆａｃｅｈｅａｔｓｏｕｒｃｅｏｎＡｓｉａｎｓｕｍｍｅｒｍｏｎｓｏｏｎ：ａｒｅｖｉｅｗ［Ｊ］．ＣｈｉｎＪＡｔｍｏｓＳｃｉ，４２（４）：７５５７６６（ｉｎＣｈｉｎｅｓｅ）．段安民，肖志祥，吴国雄，２０１６．１９７９—２０１４年全球变暖背景下青藏高原气候变化特征［Ｊ］．气候变化研究进展，１２（５）：３７４３８１．ＤｕａｎＡＭ，ＸｉａｏＺＸ，ＷｕＧＸ，２０１６．ＣｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｃｌｉｍａｔｅｃｈａｎｇｅｏｖｅｒｔｈｅＴｉｂｅｔａｎＰｌａｔｅａｕｕｎｄｅｒｔｈｅｇｌｏｂａｌｗａｒｍｉｎｇｄｕｒｉｎｇ１９７９－２０１４［Ｊ］．ＣｌｉｍａｔｅＣｈａｎｇｅＲｅｓ，１２（５）：３７４３８１（ｉｎＣｈｉｎｅｓｅ）．季国良，江灏，吕兰芝，１９９５．青藏高原的长波辐射特征［Ｊ］．高原气象，１４（４）：４５１４５８．ＪｉＧＬ，ＪｉａｎｇＨ，ＬｖＬＺ，１９９５．ＣｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｌｏｎｇｗａｖｅｒａｄｉａｔｉｏｎｏｖｅｒｔｈｅＱｉｎｇｈａｉＸｉｚａｎｇＰｌａｔｅａｕ［Ｊ］．ＰｌａｔｅａｕＭｅｔｅｏｒ，１４（４）：４５１４５８（ｉｎＣｈｉｎｅｓｅ）．蒋兴文，李跃清，２０１０．青藏高原地表辐射的气候特征［Ｊ］．资源科学，３２（１０）：１９３２１９４２．ＪｉａｎｇＸＷ，ＬｉＹＱ，２０１０．ＣｌｉｍａｔｏｌｏｇｉｃａｌｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｓｕｒｆａｃｅｒａｄｉａｔｉｏｎｏｖｅｒｔｈｅＴｉｂｅｔａｎＰｌａｔｅａｕ［Ｊ］．ＲｅｓｏｕｒＳｃｉ，３２（１０）：１９３２１９４２（ｉｎＣｈｉｎｅｓｅ）．马伟强，马耀明，胡泽勇，等，２００４．藏北高原地面辐射收支的初步分析［Ｊ］．高原气象，２３（３）：３４８３５２．ＭａＷＱ，ＭａＹＭ，ＨｕＺＹ，ｅｔａｌ，２００４．ＡｎａｌｙｓｅｓｏｎｓｕｒｆａｃｅｒａｄｉａｔｉｏｎｂｕｄｇｅｔｉｎｎｏｒｔｈｅｒｎＴｉｂｅｔａｎＰｌａｔｅａｕ［Ｊ］．ＰｌａｔｅａｕＭｅｔｅｏｒ，２３（３）：３４８３５２（ｉｎＣｈｉｎｅｓｅ）．师锐，何光碧，２０１８．移出高原后长生命史高原低涡在不同移动路径下的大尺度环流特征及差异［Ｊ］．气象，４４（２）：２１３２２１．ＳｈｉＲ，ＨｅＧＢ，２０１８．Ｃｏｍｐａｒｉｓｏｎａｎａｌｙｓｉｓｏｆｃｉｒｃｕｌａｔｉｏｎｆｅａｔｕｒｅｓｏｆｄｉｆｆｅｒｅｎｔｍｏｖｉｎｇｔｒａｃｋｓｏｆｌｏｎｇｌｉｆｅｃｙｃｌｅｐｌａｔｅａｕｖｏｒｔｅｘｅｓｍｏｖｉｎｇｏｕｔｏｆｔｈｅｐｌａｔｅａｕ［Ｊ］．ＭｅｔｅｏｒＭｏｎ，４４（２）：２１３２２１（ｉｎＣｈｉｎｅｓｅ）．王美蓉，郭栋，钟珊珊，２０１９．多源资料在青藏高原大气热源计算中的适用性分析［Ｊ］．气象，４５（１２）：１７１８１７２６．ＷａｎｇＭＲ，ＧｕｏＤ，ＺｈｏｎｇＳＳ，２０１９．ＣｏｍｐａｒｉｓｉｏｎｏｆｔｈｅｍｕｌｔｉｓｏｕｒｃｅｄａｔａｓｅｔｓｉｎｃａｌｃｕｌａｔｉｏｎｏｆｔｈｅａｔｍｏｓｐｈｅｒｉｃｈｅａｔｓｏｕｒｃｅｓｏｖｅｒｔｈｅＴｉｂｅｔａｎＰｌａｔｅａｕ［Ｊ］．ＭｅｔｅｏｒＭｏｎ，４５（１２）：１７１８１７２６（ｉｎＣｈｉｎｅｓｅ）．王艺，伯癑，王澄海，２０１６．青藏高原中东部云量变化与气温的不对称升高［Ｊ］．高原气象，３５（４）：９０８９１９．ＷａｎｇＹ，ＢｏＹ，ＷａｎｇＣＨ，２０１６．ＲｅｌａｔｉｏｎｓｏｆｃｌｏｕｄａｍｏｕｎｔｔｏａｓｙｍｍｅｔｒｉｃｄｉｕｒｎａｌｔｅｍｐｅｒａｔｕｒｅｃｈａｎｇｅｉｎｃｅｎｔｒａｌａｎｄｅａｓｔｅｒｎＱｉｎｇｈａｉＸｉｚａｎｇＰｌａｔｅａｕ［Ｊ］．ＰｌａｔｅａｕＭｅｔｅｏｒ，３５（４）：９０８９１９（ｉｎＣｈｉｎｅｓｅ）．翁笃鸣，孙治安，史兵，１９８８．中国地表净辐射的气候学研究［Ｊ］．南京气象学院学报，１１（２）：１３２１４３．ＷｅｎｇＤＭ，ＳｕｎＺＡ，ＳｈｉＢ，１９８８．ＡｃｌｉｍａｔｏｌｏｇｉｃａｌｓｔｕｄｙｏｆｔｈｅｓｕｒｆａｃｅｎｅｔｒａｄｉａｔｉｏｎｏｖｅｒＣｈｉｎａ［Ｊ］．ＪＮａｎｊｉｎｇＩｎｓｔＭｅｔｅｏｒ，１１（２）：１３２１４３（ｉｎＣｈｉｎｅｓｅ）．谢志昂，段安民，２０１７．盛夏青藏高原热源与菲律宾海对流活动的联系［Ｊ］．大气科学，４１（４）：８１１８３０．ＸｉｅＺＡ，ＤｕａｎＡＭ，２０１７．ＲｅｌａｔｉｏｎｓｈｉｐｂｅｔｗｅｅｎｔｈｅＴｉｂｅｔａｎＰｌａｔｅａｕｈｅａｔｓｏｕｒｃｅａｎｄｃｏｎｖｅｃｔｉｏｎｏｖｅｒｔｈｅＰｈｉｌｉｐｐｉｎｅＳｅａ［Ｊ］．ＣｈｉｎＪＡｔｍｏｓＳｃｉ，４１（４）：８１１８３０（ｉｎＣｈｉｎｅｓｅ）．叶笃正，高由禧，１９７９．青藏高原气象学［Ｍ］．北京：科学出版社．ＹｅｈＴＣ，ＧａｏＹＸ，１９７９．ＭｅｔｅｏｒｏｌｏｇｙｏｆｔｈｅＱｉｎｇｈａｉＸｉｚａｎｇ（Ｔｉｂｅｔ）Ｐｌａｔｅａｕ［Ｍ］．Ｂｅｉｊｉｎｇ：ＳｃｉｅｎｃｅＰｒｅｓｓ（ｉｎＣｈｉｎｅｓｅ）．于涵，张杰，刘诗梦，２０１８．基于ＣＥＲＥＳ卫星资料的青藏高原有效辐射变化规律［Ｊ］．高原气象，３７（１）：１０６１２２．ＹｕＨ，ＺｈａｎｇＪ，ＬｉｕＳＭ，２０１８．ＴｈｅｖａｒｉａｔｉｏｎｏｆｅｆｆｅｃｔｉｖｅｒａｄｉａｔｉｏｎｉｎＱｉｎｇｈａｉＴｉｂｅｔａｎＰｌａｔｅａｕｂａｓｅｄｏｎｔｈｅＣＥＲＥＳｓａｔｅｌｌｉｔｅｄａｔａ［Ｊ］．ＰｌａｔｅａｕＭｅｔｅｏｒ，３７（１）：１０６１２２（ｉｎＣｈｉｎｅｓｅ）．张镇宏，蔡景就，乔云亭，等，２０１９．青藏高原夏季大气视热源与中国东部降水的关系的年代际变化［Ｊ］．大气科学，４３（５）：９９０１００４．ＺｈａｎｇＺＨ，ＣａｉＪＪ，ＱｉａｏＹＴ，ｅｔａｌ，２０１９．ＩｎｔｅｒｄｅｃａｄａｌｃｈａｎｇｅｉｎｔｈｅｒｅｌａｔｉｏｎｂｅｔｗｅｅｎａｔｍｏｓｐｈｅｒｉｃａｐｐａｒｅｎｔｈｅａｔｓｏｕｒｃｅｓｏｖｅｒＴｉｂｅｔａｎＰｌａｔｅａｕａｎｄｐｒｅｃｉｐｉｔａｔｉｏｎｉｎＥａｓｔｅｒｎＣｈｉｎａｉｎｓｕｍｍｅｒ［Ｊ］．ＣｈｉｎＪＡｔｍｏｓＳｃｉ，４３（５）：９９０１００４（ｉｎＣｈｉｎｅｓｅ）．周允华，１９８４．青藏高原地面长波辐射经验计算方法［Ｊ］．地理学报，５１（２）：１４８１６２．ＺｈｏｕＹＨ，１９８４．Ａｎｅｍｐｉｒｉｃａｌｍｅｔｈｏｄｆｏｒ１９８４—２０１７年ＩＳＣＣＰＦＨ辐射资料在青藏高原地区的适用性评估　　　　　　　　　第１期　　　　　　　赵萌初等：２３ｅｓｔｉｍａｔｉｎｇｓｕｒｆａｃｅｌｏｎｇｗａｖｅｒａｄｉａｔｉｏｎｅｘｃｈａｎｇｅｓｉｎｔｈｅＱｉｎｇｈａｉＸｉｚａｎｇＰｌａｔｅａｕ［Ｊ］．ＡｃｔａＧｅｏｇｒＳｉｎ，５１（２）：１４８１６２（ｉｎＣｈｉｎｅｓｅ）．ＣｏｘＳＪ，ＳｔａｃｋｈｏｕｓｅＪｒＰＷ，ＧｕｐｔａＳＫ，ｅｔａｌ，２００６．ＴｈｅＮＡＳＡ／ＧＥＷＥＸｓｕｒｆａｃｅｒａｄｉａｔｉｏｎｂｕｄｇｅｔｐｒｏｊｅｃｔ：ｏｖｅｒｖｉｅｗａｎｄａｎａｌｙｓｉｓ［Ｃ］／／Ｔｈｅ１２ｔｈＣｏｎｆｅｒｅｎｃｅｏｎＡｔｍｏｓｐｈｅｒｉｃＲａｄｉａｔｉｏｎＭａｄｉｓｏｎ．Ｗｉｓｃｏｎｓｉｎ，ＵＳＡ：ＮＡＳＡ．ＣｒａｗｆｏｒｄＴＭ，ＤｕｃｈｏｎＣＥ，１９９９．Ａｎｉｍｐｒｏｖｅｄｐａｒａｍｅｔｅｒｉｚａｔｉｏｎｆｏｒｅｓｔｉｍａｔｉｎｇｅｆｆｅｃｔｉｖｅａｔｍｏｓｐｈｅｒｉｃｅｍｉｓｓｉｖｉｔｙｆｏｒｕｓｅｉｎｃａｌｃｕｌａｔｉｎｇｄａｙｔｉｍｅｄｏｗｎｗｅｌｌｉｎｇｌｏｎｇｗａｖｅｒａｄｉａｔｉｏｎ［Ｊ］．ＪＡｐｐｌＭｅｔｅｏｒ，３８（４）：４７４４８０．ＤｕａｎＡＭ，ＬｉＦ，ＷａｎｇＭＲ，ｅｔａｌ，２０１１．ＰｅｒｓｉｓｔｅｎｔｗｅａｋｅｎｉｎｇｔｒｅｎｄｉｎｔｈｅｓｐｒｉｎｇｓｅｎｓｉｂｌｅｈｅａｔｓｏｕｒｃｅｏｖｅｒｔｈｅＴｉｂｅｔａｎＰｌａｔｅａｕａｎｄｉｔｓｉｍｐａｃｔｏｎｔｈｅＡｓｉａｎｓｕｍｍｅｒｍｏｎｓｏｏｎ［Ｊ］．ＪＣｌｉｍａｔｅ，２４（２１）：５６７１５６８２．ＧａｏＧ，ＺｈａｎｇＬＹ，ＷｅｎｇＤＭ，１９９８．ＱｉｎｇｈａｉＸｉｚａｎｇｓｕｒｆａｃｅｌｏｎｇｗａｖｅｒａｄｉａｔｉｏｎｐａｒａｍｅｔｅｒｉｚａｔｉｏｎｗｉｔｈｉｔｓｃｌｉｍａｔｉｃｃａｌｃｕｌａｔｉｏｎ［Ｊ］．ＡｃｔａＭｅｔｅｏｒＳｉｎ，１２（２）：１９０１９８．ＧｕｉＳ，ＬｉａｎｇＳＬ，ＬｉＬ，２０１０ａ．Ｅｖａｌｕａｔｉｏｎｏｆｓａｔｅｌｌｉｔｅｅｓｔｉｍａｔｅｄｓｕｒｆａｃｅｌｏｎｇｗａｖｅｒａｄｉａｔｉｏｎｕｓｉｎｇｇｒｏｕｎｄｂａｓｅｄｏｂｓｅｒｖａｔｉｏｎｓ［Ｊ］．ＪＧｅｏｐｈｙｓＲｅｓＡｔｍｏｓ，１１５（Ｄ１８）：Ｄ１８２１４．ＧｕｉＳ，ＬｉａｎｇＳＬ，ＷａｎｇＫＣ，ｅｔａｌ，２０１０ｂ．Ａｓｓｅｓｓｍｅｎｔｏｆｔｈｒｅｅｓａｔｅｌｌｉｔｅｅｓｔｉｍａｔｅｄｌａｎｄｓｕｒｆａｃｅｄｏｗｎｗｅｌｌｉｎｇｓｈｏｒｔｗａｖｅｉｒｒａｄｉａｎｃｅｄａｔａｓｅｔｓ［Ｊ］．ＩＥＥＥＧｅｏｓｃｉＲｅｍｏｔｅＳｅｎｓＬｅｔｔ，７（４）：７７６７８０．ＪｉａＡＬ，ＬｉａｎｇＳＬ，ＪｉａｎｇＢ，ｅｔａｌ，２０１８．Ｃｏｍｐｒｅｈｅｎｓｉｖｅａｓｓｅｓｓｍｅｎｔｏｆｇｌｏｂａｌｓｕｒｆａｃｅｎｅｔｒａｄｉａｔｉｏｎｐｒｏｄｕｃｔｓａｎｄｕｎｃｅｒｔａｉｎｔｙａｎａｌｙｓｉｓ［Ｊ］．ＪＧｅｏｐｈｙｓＲｅｓＡｔｍｏｓ，１２３（４）：１９７０１９８９．ＫａｔｏＳ，ＲｏｓｅＦＧ，ＲｕｔａｎＤＡ，ｅｔａｌ，２０１８．Ｓｕｒｆａｃｅｉｒｒａｄｉａｎｃｅｓｏｆｅｄｉｔｉｏｎ４．０ｃｌｏｕｄｓａｎｄｔｈｅｅａｒｔｈ＇ｓｒａｄｉａｎｔｅｎｅｒｇｙｓｙｓｔｅｍ（ＣＥＲＥＳ）ｅｎｅｒｇｙｂａｌａｎｃｅｄａｎｄｆｉｌｌｅｄ（ＥＢＡＦ）ｄａｔａｐｒｏｄｕｃｔ［Ｊ］．ＪＣｌｉｍａｔｅ，３１（１１）：４５０１４５２７．ＬｈｏｍｍｅＪＰ，ＶａｃｈｅｒＪＪ，ＲｏｃｈｅｔｅａｕＡ，２００７．ＥｓｔｉｍａｔｉｎｇｄｏｗｎｗａｒｄｌｏｎｇｗａｖｅｒａｄｉａｔｉｏｎｏｎｔｈｅＡｎｄｅａｎＡｌｔｉｐｌａｎｏ［Ｊ］．ＡｇｒＦｏｒＭｅｔｅｏｒ，１４５（３／４）：１３９１４８．ＲｏｓｓｏｗＷＢ，ＳｃｈｉｆｆｅｒＲＡ，１９９１．ＩＳＣＣＰｃｌｏｕｄｄａｔａｐｒｏｄｕｃｔｓ［Ｊ］．ＢｕｌｌＡｍｅｒＭｅｔｅｏｒＳｏｃ，７２（１）：２２０．ＳｃｈｍｉｄｔＧＡ，ＲｕｅｄｙＲ，ＨａｎｓｅｎＪＥ，ｅｔａｌ，２００６．ＰｒｅｓｅｎｔｄａｙａｔｍｏｓｐｈｅｒｉｃｓｉｍｕｌａｔｉｏｎｓｕｓｉｎｇＧＩＳＳＭｏｄｅｌＥ：ｃｏｍｐａｒｉｓｏｎｔｏｉｎｓｉｔｕ，ｓａｔｅｌｌｉｔｅ，ａｎｄｒｅａｎａｌｙｓｉｓｄａｔａ［Ｊ］．ＪＣｌｉｍａｔｅ，１９（２）：１５３１９２．ＷａｎｇＢ，ＢａｏＱ，ＨｏｓｋｉｎｓＢ，ｅｔａｌ，２００８．ＴｉｂｅｔａｎＰｌａｔｅａｕｗａｒｍｉｎｇａｎｄｐｒｅｃｉｐｉｔａｔｉｏｎｃｈａｎｇｅｓｉｎＥａｓｔＡｓｉａ［Ｊ］．ＧｅｏｐｈｙｓＲｅｓＬｅｔｔ，３５（１４）：Ｌ１４７０２．ＷａｎｇＭＲ，ＺｈｏｕＳＷ，ＤｕａｎＡＭ，２０１２．ＴｒｅｎｄｉｎｔｈｅａｔｍｏｓｐｈｅｒｉｃｈｅａｔｓｏｕｒｃｅｏｖｅｒｔｈｅｃｅｎｔｒａｌａｎｄｅａｓｔｅｒｎＴｉｂｅｔａｎＰｌａｔｅａｕｄｕｒｉｎｇｒｅｃｅｎｔｄｅｃａｄｅｓ：ｃｏｍｐａｒｉｓｏｎｏｆｏｂｓｅｒｖａｔｉｏｎｓａｎｄｒｅａｎａｌｙｓｉｓｄａｔａ［Ｊ］．ＣｈｉｎＳｃｉＢｕｌｌ，５７（５）：５４８５５７．ＷａｎｇＷＨ，ＬｉａｎｇＳＬ，ＡｕｇｕｓｔｉｎｅＪＡ，２００９．ＥｓｔｉｍａｔｉｎｇｈｉｇｈｓｐａｔｉａｌｒｅｓｏｌｕｔｉｏｎｃｌｅａｒｓｋｙｌａｎｄｓｕｒｆａｃｅｕｐｗｅｌｌｉｎｇｌｏｎｇｗａｖｅｒａｄｉａｔｉｏｎｆｒｏｍＭＯＤＩＳｄａｔａ［Ｊ］．ＩＥＥＥＴｒａｎｓＧｅｏｓｃｉＲｅｍｏｔｅＳｅｎｓ，４７（５）：１５５９１５７０．ＷｉｅｌｉｃｋｉＢＡ，ＢａｒｋｓｔｒｏｍＢＲ，ＨａｒｒｉｓｏｎＥＦ，ｅｔａｌ，１９９６．ＣｌｏｕｄｓａｎｄｔｈｅＥａｒｔｈ’ｓＲａｄｉａｎｔＥｎｅｒｇｙＳｙｓｔｅｍ（ＣＥＲＥＳ）：ａｎｅａｒｔｈｏｂｓｅｒｖｉｎｇｓｙｓｔｅｍｅｘｐｅｒｉｍｅｎｔ［Ｊ］．ＢｕｌｌＡｍＭｅｔｅｏｒｏｌＳｏｃ，７７（５）：８５３８６８．ＷｕＧＸ，ＬｉＷＰ，ＧｕｏＨ，ｅｔａｌ，１９９７．ＳｅｎｓｉｂｌｅｈｅａｔｄｒｉｖｅｎａｉｒｐｕｍｐｏｖｅｒｔｈｅＴｉｂｅｔａｎＰｌａｔｅａｕａｎｄｉｔｓｉｍｐａｃｔｓｏｎｔｈｅＡｓｉａｎｓｕｍｍｅｒｍｏｎｓｏｏｎ［Ｍ］∥ＹｅＤＺ．ＣｏｌｌｅｃｔｉｏｎｓｏｎｔｈｅＭｅｍｏｒｙｏｆＺｈａｏＪｉｕｚｈａｎｇ．Ｂｅｉｊｉｎｇ：ＳｃｉｅｎｃｅＰｒｅｓｓ：１１６１２６．ＹａｎｇＫ，ＨｕａｎｇＧＷ，ＴａｍａｉＮ，２００１．Ａｈｙｂｒｉｄｍｏｄｅｌｆｏｒｅｓｔｉｍａｔｉｎｇｇｌｏｂａｌｓｏｌａｒｒａｄｉａｔｉｏｎ［Ｊ］．ＳｏｌａｒＥｎｅｒｇｙ，７０（１）：１３２２．ＹａｎｇＫ，ＫｏｉｋｅＴ，ＳｔａｃｋｈｏｕｓｅＰ，ｅｔａｌ，２００６ａ．ＡｎａｓｓｅｓｓｍｅｎｔｏｆｓａｔｅｌｌｉｔｅｓｕｒｆａｃｅｒａｄｉａｔｉｏｎｐｒｏｄｕｃｔｓｆｏｒｈｉｇｈｌａｎｄｓｗｉｔｈＴｉｂｅｔｉｎｓｔｒｕｍｅｎｔａｌｄａｔａ［Ｊ］．ＧｅｏｐｈｙｓＲｅｓＬｅｔｔ，３３（２２）：Ｌ２２４０３．ＹａｎｇＫ，ＫｏｉｋｅＴ，ＹｅＢＳ，２００６ｂ．Ｉｍｐｒｏｖｉｎｇｅｓｔｉｍａｔｉｏｎｏｆｈｏｕｒｌｙ，ｄａｉｌｙ，ａｎｄｍｏｎｔｈｌｙｓｏｌａｒｒａｄｉａｔｉｏｎｂｙｉｍｐｏｒｔｉｎｇｇｌｏｂａｌｄａｔａｓｅｔｓ［Ｊ］．ＡｇｒＦｏｒＭｅｔｅｏｒ，１３７（１／２）：４３５５．ＹａｎｇＫ，ＨｅＪ，ＴａｎｇＷＪ，ｅｔａｌ，２０１０．Ｏｎｄｏｗｎｗａｒｄｓｈｏｒｔｗａｖｅａｎｄｌｏｎｇｗａｖｅｒａｄｉａｔｉｏｎｓｏｖｅｒｈｉｇｈａｌｔｉｔｕｄｅｒｅｇｉｏｎｓ：ｏｂｓｅｒｖａｔｉｏｎａｎｄｍｏｄｅｌｉｎｇｉｎｔｈｅＴｉｂｅｔａｎＰｌａｔｅａｕ［Ｊ］．ＡｇｒｉｃＦｏｒＭｅｔｅｏｒ，１５０（１）：３８４６．ＹｏｕｎｇＡＨ，ＫｎａｐｐＫＲ，ＩｎａｍｄａｒＡ，ｅｔａｌ，２０１８．ＴｈｅｉｎｔｅｒｎａｔｉｏｎａｌｓａｔｅｌｌｉｔｅｃｌｏｕｄｃｌｉｍａｔｏｌｏｇｙｐｒｏｊｅｃｔＨｓｅｒｉｅｓｃｌｉｍａｔｅｄａｔａｒｅｃｏｒｄｐｒｏｄｕｃｔ［Ｊ］．ＥａｒｔｈＳｙｓｔｅｍＳｃｉｅｎｃｅＤａｔａ，１０（１）：５８３５９３．ＺｈａｎｇＹＣ，ＲｏｓｓｏｗＷＢ，ＬａｃｉｓＡＡ，ｅｔａｌ，２００４．ＣａｌｃｕｌａｔｉｏｎｏｆｒａｄｉａｔｉｖｅｆｌｕｘｅｓｆｒｏｍｔｈｅｓｕｒｆａｃｅｔｏｔｏｐｏｆａｔｍｏｓｐｈｅｒｅｂａｓｅｄｏｎＩＳＣＣＰａｎｄｏｔｈｅｒｇｌｏｂａｌｄａｔａｓｅｔｓ：ｒｅｆｉｎｅｍｅｎｔｓｏｆｔｈｅｒａｄｉａｔｉｖｅｔｒａｎｓｆｅｒｍｏｄｅｌａｎｄｔｈｅｉｎｐｕｔｄａｔａ［Ｊ］．ＪＧｅｏｐｈｙｓＲｅｓＡｔｍｏｓ，１０９（Ｄ１９）：Ｄ１９１０５．ＺｈａｎｇＹＣ，ＲｏｓｓｏｗＷＢ，ＳｔａｃｋｈｏｕｓｅＪｒＰＷ，２００６．Ｃｏｍｐａｒｉｓｏｎｏｆｄｉｆｆｅｒｅｎｔｇｌｏｂａｌｉｎｆｏｒｍａｔｉｏｎｓｏｕｒｃｅｓｕｓｅｄｉｎｓｕｒｆａｃｅｒａｄｉａｔｉｖｅｆｌｕｘｃａｌｃｕｌａｔｉｏｎ：ｒａｄｉａｔｉｖｅｐｒｏｐｅｒｔｉｅｓｏｆｔｈｅｎｅａｒｓｕｒｆａｃｅａｔｍｏｓｐｈｅｒｅ［Ｊ］．ＪＧｅｏｐｈｙｓＲｅｓＡｔｍｏｓ，１１１（Ｄ１３）：Ｄ１３１０６．