第47卷第2期气象Vol.47No.2 2021年2月METEOROLOGICAL MONTHLY February2021
孟琦,2021.考虑气压的晴空大气下行长波辐射参数化方案及其适用性分析[J].气象，47(2):183191.Meng Q,2021.Parame-terizations of downward longwave radiation in clear sky considering air pressure and their applicability analysis[J].Meteor Mon,47(2)：183-191(in Chine).
考虑气压的晴空大气下行长波辐射参数化方案
及其适用性分析*
孟琦12
1中国科学院青藏高原研究所，北京100085
2中国科学院大学，北京100049
提要：基于Prata晴空大气下行长波辐射参数化方案，针对其在高原地区及可降水量较小地区理论精度较差的缺点，考虑不同高度的大气柱发射的长波辐射量不同，提出了三种考虑气压的大气下行长波辐射参数化方案，通过全球ERA-5再分析数据进行最小二乘拟合确定经验常数并在全球不同区域对其适用性进行了分析。模拟结果表明，考虑气压的参数化方案有效地改善了可降水量减小时Prata方案大气发射率收
敛过快的不足，三种新方案在青藏高原安多地区和南美洲圣路易斯地区的下行长波辐射模拟结果相对于观测值和ERA-5再分析数据的平均偏差和均方根误差均更小，模拟精度相对Prata方案均有提高，适用性更强。
关键词：晴空下行长波辐射，参数化方案，气压，适用性分析
中图分类号：P413文献标志码：A DOI：10.7519/j.issn.1000-0526.2021.02.005
Parameterizations of Downward Longwave Radiation in Clear Sky Considering Air Pressure and Their Applicability Analysis
MENG Qi1，
1Insttute of Tibetan Plateau Rearch,Chine Academy of Sciences,Beijing100085
2Universty of Chine Academy of Sciences,Beijing100049
Abstract:Bad on Prata parameterization of clear sky downward longwave radiation and in view of its poor theoretical accuracy in plateau areas and areas with small precipitable water vapor,three new parameterizations of clear sky downward longwave radiation are propod in this paper,in whic
h air pressure is considered becau longwave radiation from different length of atmosphere column is different.Empirical constants are obtained by least square fitting of global ERA-5reanalysis data and their applicability is analyzed in different regions of the world.The simulation results show that parameterizations with consideration of airpressurecane f ectivelyimprovetheconvergenceofatmosphericemissivityofPratascheme whenthe precipitable water vapor can be reduced.The mean bias errors and the root mean square errors of the simulation values of the downward longwave radiation of the new schemes in Amdo Area and San Luis from the obrvation values are obviously smaller than that of Prata?s parameterization scheme.Compared with Prata?s scheme,the simulation accuracies of the new schemes are improved and the applicabilities are better.徐向前元帅
Key words:clear sky downward longwave radiation,parameterization,air pressure,applicability analysis
国家重点基础研发计划(2015CB953701)资助
2020年2月11日收稿；2020年8月19日收修定稿
作者：孟琦，主要从事定量遥感及水文气象方向的研究.E-mail:mengqi@itpcas.ac
184气象第47卷
引言
大气发射的下行长波辐射是全球辐射能量收支平衡中一个非常重要的成分，也是地表辐射模型中辐射预算的关键组成部分，准确计算大气下行长波辐射有利于了解影响地表能量平衡的因素(Zhu et al,2017)。同时，大气下行长波辐射研究在预报夜间霜冻、雾、温度变化、云量及改善区域天气预报、确定辐射冷却率、衡量气候变化和计算全球增暖等方面也发挥着重要作用(彭丽春等，2015)。
目前获得晴空下行长波辐射的方法主要有三种：地面气象站辐射计直接观测、根据复杂的大气辐射传输模型精确计算及通过经验公式计算获得。辐射计直接观测能方便准确地得到下行长波辐射值，但地面观测站点分布不均，因此大部分地区不具备直接观测的条件，目前只有部分国家气象观测站能获取相关数据(李伟斌等，2015)。在已知整层大气的垂直分布特性的条件下利用探空数据，通过复杂的大气辐射传输模型，可以获得下行长波辐射的理论值，然而由于探空数据不易获取，使用上述辐射传输模型计算时，常受到输入数据缺乏的限制(周允华，1984)。通过经验公式进行计算具有模型简单、观测数据(近地面气温、气压、水汽压等)易获得的优势，因此具有很强的实用性，在许多领域得到了广泛的应用(余珊珊等，2011；Alados et al,2012)。例如Zhu et al(2017)在青藏高原地区对Brutsaert (1975)、Garratt(1992)、Prata(1996)和Swinbank (1963)等13个参数化方案进行了局地校正，获得了
适用于该地区的参数化方案。近年来，下行长波辐射参数化研究取得了丰硕成果，例如Prata(1996)提出的参数化方案有效地改善了极端气候条件下的模拟精度，物理意义更加严谨，且适用性更强。Yu et al(2013)利用卫星热红外数据提出了具有高分辨率优势的参数化方案，后来又基于MODIS数据对该方案进行了改进，将模拟精度和适用性进一步提高(Yu et al,2019)。Zhou et al(2019)利用中分辨率成像光谱仪数据及地面实测资料，基于多元自适应回归方法进行了晴空大气下行长波辐射计算，并在全球范围7个站点进行了验证。随着遥感科学的发展，气象观测的时空分辨率不断提高，下行长波辐射参数化方案模拟精度和适用性也得到了改善，但目前基于大气参数的参数化方案通常不具备广泛的适用性，尤其是在高海拔或干旱地区，通常存在较大的偏差，且误差会随着海拔的升高而增大，制约了参数化方案的适用性(Gubler et al,2012)。狗子图片
我国幅员辽阔，地形与气候条件较为复杂，利用上述参数化方案对晴空下行长波辐射进行理论计算的精确度难以满足要求。以青藏高原为例，由于其地处高寒地区，地面有效辐射和辐射平衡具有大致相同的量值，因此，下行长波辐射的精确估计对于高原地面辐射平衡和加热场的研究具有重要意义(周允华,1984；黄妙芬等，2005)。其次，入射长波辐射是青藏高原地区冰川融化的能量来源，利用准确的下行长波辐射值，有助于获得准确的模拟蒸散和冰川排放的结果，这对于冰川灾害预防、绿洲农业用水甚至西藏地区社会经济的整体发展都是必不可少的(Zhuetal，2017；Ohmura，2001；Sedlarand Hock，2009)。此外，研究高寒地区晴空下行长波辐射对光能的利用具有指导性意义，可为新的清洁
能源的开发提供积极广阔的前景。同时青藏高原对我国的气候变化有着重要的影响，精确地定量研究该地区的晴空下行长波辐射，对研究我国天气气候变化也有着重要意义。然而晴空下行长波辐射数据的缺乏，是能量平衡、气候模拟等研究常常面临的困难(余珊珊等，2011；闵敏和吴晓，2019)。
针对目前晴空大气下行长波辐射参数化方案在青藏高原及其他干旱地区适用性较差的问题，本研究基于Prata的参数化方案，考虑到入射到地面的长波辐射来自整层大气，对于不同海拔的地区来说，大气柱的高度不同，下行长波辐射值也不同，由此引进气压因素，根据ERA-5的再分析数据，通过多元函数拟合，获得了新形式的晴空大气发射率参数化方案，进而根据站点资料在全球不同海拔地区进行了模拟，与站点资料的下行长波辐射的观测值进行了对比，分析了新参数化方案在青藏高原及其他地区的适用性。
1物理背景
入射到地面的大气下行长波辐射是其上整层大气发射的，与大气层温度和湿度的垂直分布、厚度密切相关，理论上可表示为垂直分层大气的积分形式。Steffan-Boltzman公式是近似计算长波辐射的常用公式，但该公式只适用于黑体，对于非黑体的大气来说，则需要考虑大气发射率。根据该公式，在晴空条
第2期孟琦：考虑气压的晴空大气下行长波辐射参数化方案及其适用性分析185
件下，大气下行长波辐射计算公式可近似表示为：
LW d=e犜(1)式中LW a为晴空大气下行长波辐射，单位为W•m-2；a为晴空大气发射率犜通常取近地面大气温度代替整层大气的发射等效温度，单位为K；。为玻尔兹曼常数，其值为5.67X10-8W•m^2•。
由于晴空大气发射率无法通过直接测量获得，Brunt(1932)基于热传导与辐射传递之间的相似性，提出了晴空大气发射率与近地面水汽压相关的经验公式(Crawford and Duchon,1999):
e=犪+犫槡(2)式中：e为近地面水汽压，单位为hPa；a、为经验常数，值约在0.34-0.66，值约在0.033-0.127。由于经验常数因地区而异，给应用带来了不便。Swinbank(1963)指出下行长波辐射与水汽压之间并无关系，而仅仅与温度的平方有关，且Brunt (1932)经验公式有效的原因是水汽压与温度的正相关关系，进而提出了仅与近地面气温相关的经验公式(盛裴轩等，2013)：
E a=aT a(3)式中：常数a=5.31X10-14/ff；T a为近地面气温，单位为K
由于上式的常数不随地点而变，观测到的和估计的晴空大气下行长波辐射之间的相关系数为0.99，因而被广泛采用，但在高海拔地区，因水汽吸收的压力效应，需做高度校正(盛裴轩等，2013)。Brutsaert(1975)基于Schwarzschild方程和标准大气假设，提出了物理意义更严格的参数化方案：
从式中可以看出当水汽压为0时，该参数化方案得到的晴空大气发射率为0，但实际上，即使水汽含量
为0，具有一定温度的大气同样能发射长波辐射，因此晴空大气发射率不应为0，故此结果并不能完全适用于计算晴空大气发射率。Prata(1996)指出，基于对单个气体吸收带的窄带模型和宽带模型的了解，可以用修正的指数带模型近似表示整个长波频谱，同时应考虑对水汽路径不均匀性加以修正，据此将辐射传输过程进行了高度简化近似，考虑了温度、湿度的垂直影响，引入了大气可降水量，提出了以下形式的参数化方案：
良=1-(1+w)e Wa1+a2ror](5)式中：犪1，2和犿是经验常数，犪1和犪2通过最小二乘拟合的方法获得，分别取1.2和3.0，犿取1/2。狑是大气可降水量，单位为g-cm-2，可以根据以下方式进行计算：
狑=犮犲(6)
T a
吴世民
翡翠吊兰式中：犮是一个常数，标准大气压假设下其取值为46.5。该参数化方案具有以下良好特性：当大气可降水量趋近于0，即大气接近为干空气时，只要具有一定温度，大气依然发出下行长波辐射，此时晴空大气发射率是一个常数，且该参数化方案估算的晴空大气发射率不会超过黑体的发射率。需要说明，随着海拔高度的升高，此参数化方案结果与观测值的均方根误差会增大，例如在海拔为3580m时，均方根误差达到长波辐射值的15%(Prata,1996)。鉴于此，本研究认为大气柱的厚度会对大气发射率产生影响，而长波辐射参数化方案在高寒地区的应用有望结合气压加以改进。
2引入气压的参数化方案
由于Prata参数化方案具有干空气条件下，下行长波辐射依然存在和晴空大气发射率不超过黑体发射率的优点，本研究在Prata参数化方案的基础上，引入了气压因素，提出了以下三种形式的晴空大气发射率的参数化计算方案：
S1-(1+狑)(1+”犪皿(7)
1-(1+狑))1+了)—〔狑犫+/+犳狑z鬥犮(8)
e〜1—(9)式中：Z=(狆/狆0)狆为当地大气压，狆为标准大气压，狑是大气可降水量，犪、b、c、d、f、g均为经验常数。这里分别称式(7)〜式(9)为方案1、方案2和方案3，上述三种方案中当气压趋近于0时，大气可降水量狑趋近于0，则此时晴空大气发射率也趋近于0，且引入气压项后，晴空大气发射率依然不会超过黑体发射率。三种方案均保留了Prata方案在极限条件下的优良特性，但表达形式不同，大气发射率模拟结果也不相同。方案1参照了Prata进行水汽改正时的方法，在指数部分和系数部分同时添加气压项，方案2和方案3在方案1的基础上考虑了可降水量与气压的交互影响，因此理论上极限情况下表现更好。以上三种参数化方案中的经验常数通过使用2018年每月1日ERA-5全球再分析数据进行最小二乘拟合的方法获得。再分析数据包含全球范
186
气 象
第47卷
围的近地面气温、露点温度、气压、大气可降水量以 及晴空大气下行长波辐射值，时间分辨率为1 h，空
间分辨率为0.25。。为了充分考虑不同气压条件下 长波辐射的差异，本研究对拟合数据进行了筛选，将
每小时全球各点的气压值均从小到大排序并形成十 个区间，在每个气压区间中随机挑选100个气压值, 对应100个地面点，拟合样本容量为288000个。作
为示例，图1给出了经筛选后2018年1月1日00 时全球范围内用于拟合的数据对应的地面点分布,
红色的点表示此时经过筛选的站点。
从图1中可以看出，青藏高原地区等高海拔地
区数据覆盖程度较高，因此晴空数据筛选结果能够 更好地保证高海拔地区拟合结果的准确性，这也正
是对再分析数据进行筛选的目的。在此基础上经过
最小二乘拟合处理后，得到三种考虑气压影响的晴
空下行长波辐射参数化方案的经验系数见表1。
schemes considering air pressure
表1考虑气压的参数化方案经验系数拟合结果
Table 1 Empirical coefficients of parameterized
方案a
犫C
犱
寄的组词是什么/
犵
方案1
3．214
2. 6010 548/
//
方案21．8350 5581 712
0 706— 0.851
0. 381
方案3
0. 179
0 934
0 5980 613
0567 0 130
为了研究上述三种参数化方案对晴空大气发射 率的模拟效果，本文将对应时间ERA-5再分析数据
中获得的所有站点的晴空大气发射率作为参照（即 实际值），分别求得三种参数化方案模拟结果的平均
偏差、均方根误差及与参照大气发射率的相关系数
并做散点密度图，同时将Prata 模型的模拟结果作 为对照，比较不同的参数化方案的模拟精度。
图2a 〜2d 分别给出了本研究提出的三种新参 数化方案及Prata 参数化方案在全球按海拔均匀分
布的多个随机站点（包含平原地区和高原地区）的晴
空大气发射率散点密度，模拟时间为2018年每月1 日，时间分辨率为1 h，样本容量为288000。图中横 轴为根据ERA-5再分析数据中晴空下行长波辐射 得到的晴空大气发射率，根据式（1）有：
LW d 丁 4O l a
(10)
纵轴为考虑气压的参数化方案与Prata 参数化方 案，分别通过ERA-5气压与可降水量数据计算得到 的晴空大气发射率理论值，颜色表示散点密度，颜色 越亮则该区域散点密度越大。从图中可以看出，考
虑气压影响的三种参数化方案得到的晴空大气发射
率理论值与实际值基本上表现出斜率为1的线性关 系，而对于Prata 参数化方案而言，由于其在大气可 降水量不断减小时大气发射率逐渐收敛于0. 665左
右，因此对于大气发射率小于该极限值的情况无法
进行模拟，这也解释了该参数化方案应用于大气可 降水量较小的地区时而造成精度较低的问题。三种
考虑气压的参数化方案与Prata 参数化方案晴空大 气发射率模拟结果相较于ERA-5数据的平均偏差
（MBE ）、均方根误差（RMSE ）及相关系数（R ）计算
结果见表2。
表2新参数化方案与Prata 参数化方案晴空大气
发射率相对于ERA-5数据的MBE 、RMSE 和R
Table 2 Mean bias error , root mean square error and correlation coefficient of new parameterizations
andPrataparameterizationcorrespondingto
氧平衡ERA-5 reanalysis data
统计量方案1
方案2
方案3
Prata 参数化方案
MBE
0 0017. 73X10 —5
4 36X10—5
0 068RMSE 0 0420 0380 0380 088R
0 900
0 920
0 919
0 862
从表2中可以看到，相对于ERA-5再分析数据，
方案2和方案3的MBE 较小，其中方案3的MBE 仅
为4. 36 X 10-5 , Prata 参数化方案的MBE 较大，为
0.068；对于RMSE 和R 而言，方案1、方案2及方案3
亦均显著优于Prata 参数化方案，其中方案2和方案3 精度更高，相关性更好。因此，本研究所提出的三种
参数化方案相比Prata 的参数化方案能更精确地对 晴空大气发射率进行模拟， 尤其是在高海拔大气可 降水量较小的区域，理论上能够弥补Prata 参数化 方案过早收敛的缺陷。
3适用性分析
上一节中晴空大气发射率参数化方案的拟合数
据是经过筛选的全球数据，因此理论上三种新的参 数化方案可对全球范围具有不同气压和可降水量条 件地区的晴空下行长波辐射进行模拟。 为了进一步
分析三种新的参数化方案在青藏高原及其他地区的
适用性，本研究对青藏高原及其他地区的晴空下行长 波辐射进行了理论计算，以位于青藏高原的安多县
（海拔5 200 m ）和位于阿根廷圣路易斯（海拔500 m ） 的两个站点为例，本文给出了利用站点数据对新方
案与Prata 方案的模拟效果进行比较的结果。
第2期孟琦：考虑气压的晴空大气下行长波辐射参数化方案及其适用性分析187
对安多地区晴空下行长波辐射的模拟使用了该
地区1998年5月11日至9月16日的站点数据，包
括气压、近地面气温、相对湿度等观测值。本研究通 过计算云量的方法获取晴空时间(Zhu et al,2017):
式中：狀表示云量，S m 表示实测太阳下行短波辐射 值,S ^r 为晴空时太阳下行短波辐射值，二者单位均
为W ・m -2，后者可以通过理论模型计算获得。由 于夜间太阳下行短波辐射减弱为0 W ・m -2，因此
上述方法无法获取夜晚晴空时间，本研究仅取太阳
图1 2018年1月1日00时拟合数据对应
纲手的耐力测试的地面点分布
Fig. 1 Ground points corresponding to the data
ud for fitting at 00 : 00 BT 1 January 2018
0.4 0.5 0.6 0.7 0.8 0.9 1.0
(C )
(a)
0.9
0.8
0.7
输卸炭轻蕭庇鼬
44
咪®
如菸
K
卑
1.00.6根据下行长波辐射计算得到的晴空大气发射率0.6
根据下行长波辐射计算得到的晴空大气发射率
0.9
(d)
图2方案1(a)，方案2(b)，方案3(c)和Prata 方案(d)晴空大气发射率模拟结果散点密度
Fig. 2 Simulation results of clear sky atmospheric emissivity of Scheme 1 (a ),
Scheme 2 (b ) , Scheme 3 (c ) and Prata's Scheme (d )
天顶角小于70°且狀值小于0.05时作为晴空。用于 计算的大气可降水量无法通过观测获得，目前常用 的推算大气可降水量的方法包括探空资料计算、地
面气象资料推算、地基GPS 探测资料反演等。由于 目前探空资料和GPS 数据较为缺乏，本研究中采用
站点资料，根据经验公式计算青藏高原地区大气可
降水量(Yang et al,2010)：
狑=0. 004 93 犚ge (26・23—犜6)
(12)
式中:RH a 为站点的相对湿度；狑为大气可降水量,五灵丸
