免费文献传递   相关文献

Applicability of daytime downward longwave radiation parameterized models in Changbai Mountains, Northeast China.

大气向下长波辐射参数化模型在长白山地区的适用


利用中国科学院长白山森林生态系统定位站的近地面气象观测数据,分析评价了目前被广泛使用的8个晴天与8个云天大气长波辐射参数化模型的模拟性能.结果表明: 晴天时Satterlund模型最适用,其偏差(BIAS)与均方根误差(RMSE)分别是-23.34和28.55 W·m-2;系数校正后,虽然其参数值变化不大,但其模拟效果有很大提高,BIAS与RMSE分别降低为-6.33和18.08 W·m-2;云天时Jacobs模型最准确,BIAS和RMSE只有0.38和29.29 W·m-2.对模型中大气发射率的敏感性分析表明,大气发射率对水汽压的变化最敏感,对温度的变化不敏感.应用优选模型(晴天和云天)得到的模拟值与观测值的日变化趋势基本一致,但在云量发生突变的节点上模拟效果不太理想.

 

A total of eight clearsky and eight cloudysky parameterized models for estimating daytime downward longwave radiation were evaluated by using the meteorological data measured in the Changbai Mountains region, Northeast China. The results indicated that the Satterlund model performed better in estimating clearsky downward longwave radiation, and the bias (BIAS) and root mean square error (RMSE) were -23.34 and 28.55 W·m-2, respectively. Although the coefficients were not significantly changed, the performance of Satterlund model was significantly improved after the locally calibrated, and the BIAS and RMSE decreased to -6.33 and 18.08 W·m-2, respectively. Jacobs model was found to be best for modeling cloudysky downward longwave radiation and the BIAS and RMSE were 0.38 and 29.29 W·m-2, respectively. Sensitivity analysis showed that the vapor pressure  was the most sensitive variable to the atmospheric emissivity and the temperature was not sensitive to it. The predicted results of the optimal model (clear and cloudysky) were consistent with the measured data, while the performance of these models was
affected by a sudden change of cloudy cover.


全 文 :大气向下长波辐射参数化模型在
长白山地区的适用性∗
李伟斌1,2  吴家兵1  王安志1  关德新1  金昌杰1∗∗
( 1中国科学院沈阳应用生态研究所森林与土壤生态国家重点实验室, 沈阳 110164; 2中国科学院大学, 北京 100049)
摘  要  利用中国科学院长白山森林生态系统定位站的近地面气象观测数据,分析评价了目
前被广泛使用的 8个晴天与 8个云天大气长波辐射参数化模型的模拟性能.结果表明: 晴天
时 Satterlund 模型最适用,其偏差 (BIAS)与均方根误差 (RMSE)分别是 - 23. 34 和 28. 55
W·m-2;系数校正后,虽然其参数值变化不大,但其模拟效果有很大提高,BIAS与 RMSE分别
降低为-6.33和 18.08 W·m-2;云天时 Jacobs 模型最准确,BIAS 和 RMSE 只有 0.38 和 29.29
W·m-2 .对模型中大气发射率的敏感性分析表明,大气发射率对水汽压的变化最敏感,对温度
的变化不敏感.应用优选模型(晴天和云天)得到的模拟值与观测值的日变化趋势基本一致,
但在云量发生突变的节点上模拟效果不太理想.
关键词  长白山; 下行长波辐射; 参数化模型; 辐射模型
文章编号  1001-9332(2015)02-0497-08  中图分类号  P422.3  文献标识码  A
Applicability of daytime downward longwave radiation parameterized models in Changbai
Mountains, Northeast China. LI Wei⁃bin1,2, WU Jia⁃bing1, WANG An⁃zhi1, GUAN De⁃xin1,
JIN Chang⁃jie1 ( 1State Key Laboratory of Forest and Soil Ecology, Institute of Applied Ecology, Chi⁃
nese Academy of Sciences, Shenyang 110164, China; 2University of Chinese Academy of Sciences,
Beijing 100049, China) . ⁃Chin. J. Appl. Ecol., 2015, 26(2): 497-504.
Abstract: A total of eight clear⁃sky and eight cloudy⁃sky parameterized models for estimating day⁃
time downward longwave radiation were evaluated by using the meteorological data measured in the
Changbai Mountains region, Northeast China. The results indicated that the Satterlund model per⁃
formed better in estimating clear⁃sky downward longwave radiation, and the bias (BIAS) and root
mean square error (RMSE) were -23.34 and 28.55 W·m-2, respectively. Although the coeffi⁃
cients were not significantly changed, the performance of Satterlund model was significantly im⁃
proved after the locally calibrated, and the BIAS and RMSE decreased to - 6. 33 and 18. 08
W·m-2, respectively. Jacobs model was found to be best for modeling cloudy⁃sky downward long⁃
wave radiation and the BIAS and RMSE were 0.38 and 29.29 W·m-2, respectively. Sensitivity
analysis showed that the vapor pressure was the most sensitive variable to the atmospheric emissivity
and the temperature was not sensitive to it. The predicted results of the optimal model (clear⁃ and
cloudy⁃sky) were consistent with the measured data, while the performance of these models was
affected by a sudden change of cloudy cover.
Key words: Changbai Mountains; downward longwave radiation; parameterized model; radiation
model.
∗国家自然科学基金项目(31370614,41371064)资助.
∗∗通讯作者. E⁃mail: cjjin@ iae.ac.cn
2014⁃06⁃12收稿,2014⁃11⁃25接受.
    大气下行长波辐射(简称为大气长波辐射)是
地面与大气能量交换的关键分量,是影响全球能量
平衡的重要组成部分[1],在天气预报、能量平衡及
全球气候变化等领域的研究中起着非常重要的作
用[2-5] .目前,关于长波辐射的数据主要通过以下 3
种方法获得:直接观测[6]、利用简单的经验模型计
算[2,4,7]以及用 LOWTRAN、MODTRAN、STREAMER
等根据辐射传输方程转换得到的复杂模型计
算[8-10] .由于测量长波辐射的仪器价格较短波辐射
仪器昂贵[11-12],因此目前除了一些国家气象观测站
有零星数据外,再很难找到相关的准确数据.辐射的
应 用 生 态 学 报  2015年 2月  第 26卷  第 2期                                                         
Chinese Journal of Applied Ecology, Feb. 2015, 26(2): 497-504
复杂模型在日常应用中由于其输入数据的短缺而存
在很大限制性[13];相反,简单的经验参数化模型直
接利用大气温度和水汽压数据就能计算出大气长波
辐射值,而大气温度和水汽压数据在各个气象站都
较易观测获得,所以运用简单的参数化模型来量化
研究长波辐射受到很多学者的关注.
国外对长波辐射的量化研究起步较早,已构建
了很多性能较好的模型,这些模型大部分通过建立
长波辐射与环境因子(大气温度、水汽压)间的经验
关系式而得,其仅适用于特定地区(或生态系统)而
不能被广泛应用.为了提高现有模型的适用性,很多
学者对现有模型做了很多改进,如 Prata[14]引入水
汽密度来降低地表温度廓线对模拟结果的影响;Jin
等[15]利用大气探空数据改进了 Brutsaert[16]模型的
主要参数,这些改进都增强了模型的物理基础.国内
的相关研究起步较晚,还没有提出能被广泛利用的
模型.现有的研究大多是对已有模型的改进,如黄妙
芬等[17]利用北京小汤山国家精准农业示范基地的
数据对几个常见的辐射模型进行了改进;余珊珊
等[18]比较分析了一些常用模型在中国黑河流域的
模拟性能,并提出了适用于中国西北干旱半干旱地
区的长波辐射模型.虽然经验模型相对简单实用,但
大多数经验模型是利用某一地区的数据建立,所以
在使用过程中都因研究区域的地形、气候条件差异
而存在局地性.
长白山地区位于我国吉林省东部,具有海拔高、
云量多的独特区域特点,该区具有我国典型的温带
森林生态系统[19-20] .本研究利用中国科学院长白山
森林生态系统定位站的观测数据对最常用的及新近
提出的 8个晴天参数化模型和 8个云天的参数化模
型在该类区域内的适用性进行了分析,得到了最适
合该类区域的晴天与云天的模型以及系数,还分析
了模型对水汽压和大气温度的敏感性及模拟值和观
测值的日变化趋势.
1  研究地区与研究方法
1􀆰 1  研究区概况
试验数据来源于中国科学院长白山森林生态系
统定位站,该站位于中国吉林省安图县二道白河镇
(41°42′45″—42°45′18″ N,127° 33′30″ —128° 16′48″
E),该地海拔 736 m,年均降水量 600 ~ 1000 mm,年
均气温 3.5 ℃,属典型的季风温带大陆山地气候.气
象观测场位于定位站内,周围数十千米内地势平坦.
观测场按国家基本气象站标准设计:南北和东西边
长 35 m×25 m,地面为均质草皮,四周是大范围的天
然阔叶红松林[2] .
1􀆰 2  参数化模型
1􀆰 2􀆰 1晴天参数化模型   有研究表明,地面接收到
的大气长波辐射有近 80%来源于地面以上 500 m内
各层大气向下发射的长波辐射[22],它主要受大气温
度、大气成分(水汽和 CO2 等)和云的影响[23] .在晴
天无云的情况下,长波辐射模型的基本形式为:
R↓lc =εaσTa 4 (1)
式中:R↓lc为晴天大气长波辐射(W·m
-2);εa为大
气发射率;σ 为玻尔兹曼常数(σ = 5.670373 × 10-8
W·m-2·K-4);Ta为近地面大气温度(K).此公式与
斯特藩-玻尔兹曼定律相近,对于温度均匀的黑体
来说,εa为常数,而地球表面的大气在垂直方向上存
在明显的温度变化廓线关系,因此在实际计算中大
气发射率通常是用近地面气温 Ta(K)和水汽压 ea
(hPa)估算而得.
Brunt[24]通过经验数据推出发射率的表达式
为:
εa = 0.65+0.48·ea 0.5 (2)
发射率仅是水汽压的函数. Idso 等[13]认为,发
射率仅为大气温度的函数,并提出计算表达式:
εa = 1-0.261exp[-7.77×10
-4·(273-Ta) 2] (3)
大部分学者认为,发射率是受水汽压和温度综
合影响的,并提出了很多应用广泛的计算表达式.表
1中列出了应用较广泛和新发表的晴天参数化模
型.其中,Prata[14]提出的模型、 Dilley等[25]提出的模
型是基于辐射传输理论推导而来,并经过水汽密度
校正;Brutsaert[16]提出的模型是基于对近似标准大
气 Schwarzschild方程的解析解法,这 3 个模型具有
一定的物理意义,更易被接受应用;Brunt[24]提出的
模型纯属经验模型.Abramowitz 等[26]提出的模型是
利用全球 10个通量观测网点的数据进行验证而得,
可直接应用于云天,由于它只是大气温度与水汽压
的函数,用此模型时就不再需要繁琐的云量测定.
1􀆰 2􀆰 2云天参数化模型  云的存在会增加到达地面
的长波辐射,且这种影响非常复杂,为了得到云天长
波辐射的计算模型,很多研究引入一个参数———云
量.云量一般用目视观测,将其定义为把天空分为 10
份时云所占的比例.由于这种观测存在一定主观性,
因此很难进行各站点之间的比较.本文利用 Craw⁃
ford等[3]提出的估计云量( cf)的方法,其估计公式
为:
cf = 1-Rs / Rso (4)
894 应  用  生  态  学  报                                      26卷
表 1  晴天与云天下行长波辐射参数化模型
Table 1  Downward longwave radiation models for clear⁃sky and cloudy⁃sky
晴天下行长波辐射参数化模型
Clear⁃sky longwave radiation model
文献
Reference
公式
Formula
系数
Coefficient
云天下行长波辐射参数化模型
Cloudy⁃sky longwave radiation model
文献
Reference
公式
Formula
[23] R↓lc =(a1+b1 ea 0.5) σTa 4
a1 = 0.605
b1 = 0.048
[32] R↓l =R↓lc·(1+0.22cf 2.75)
[24]
R↓lc ={(1-a2·exp[ -b2·(273-Ta ) 2 ]}
σTa 4
a2 = 0.261
b2 = 7.77×10-4
[31] R↓l =R↓lc·(1+0.26cf)
[15] R↓lc =a3·(ea / Ta) b3σTa 4
a3 = 1.24
b3 = 0.14
[33] R↓l =R↓lc·(1+0.0496cf 2.45)
[28] R↓lc =a4·[1-exp(-ea
T a / b4)] σTa 4 a4 = 1.08
b4 = 2016
[34] R↓l =R↓lc·(1-cf 4)+0.952·cf 4σTa 4
[13]
R↓lc =(1-[1+ a5·ea / Ta]·exp[ -b5 +c5 a5
·ea / Ta]) 0.5) σTa 4
a5 = 46.5,b5 = 1.20
c5 = 3.0
[3] R↓l =R↓lc·(1-cf)+ cfσTa 4
[25]
R↓lc =a6+b6·(Ta / 273.16) 6 +d6·ea / Ta·
c6 / 25
a6 = 59.38,b6 = 113.7
c6 = 96.96,d6 = 46.5
[27] R↓l =R↓lc·(1+0.242·cf 0.583)
[29] R↓lc =a7·(ea / Ta) b7σTa 4
a7 = 0.576
b7 = 0.202
[27]
R↓l = R↓lc ·(1 - cf 0.671 ) + 0. 990· cf 0.671
σTa 4
[26] R↓lc =a8 ea+b8Ta-c8
a8 = 0.031,b8 = 2.84
c8 = 522.5
[26] R↓lc =a8 ea+b8Ta-c8
R↓lc: 晴天大气长波辐射 Longwave radiation in clear⁃sky (W·m-2); R↓l: 云天大气长波辐射 Longwave radiation in cloudy⁃sky (W·m-2); Ta:
近地面大气温度 Surface air temperature (K); ea: 近地面水汽压 Surface air vapor (hPa); cf: 云量 Fractional cloud cover.
式中:Rs为地表观测到的太阳总辐射(短波);Rso为
下行太阳总辐射的理论值(短波).
应用较广泛的云天参数化模型(表 1)主要有以
下两种形式:
R↓l = R↓lc·(1+mcf n) (5)
R↓l = R↓lc·(1+cf t)+scf tσTa 4 (6)
式中:m、n、s 和 t 是与云量有关的参数;cf为云量.但
也有特例,Abramowitz等[26]提出的模型中就没有考
虑云量,仅仅是水汽压与温度的函数.
1􀆰 3  数据选取及处理
根据数据质量选取定位站 2004年 1月 1日到 9
月 1日 8:00—17:30 的所有气象观测数据.所选数
据包括观测塔 32 m高处测得的太阳总辐射、反射辐
射、大气长波辐射、大气温度以及气压等.其中,太阳
总辐射、反射辐射采用 CM11短波辐射计测量(Kipp
& Zonen , Netherlands) (0.1 W·m-2);大气长波辐
射采用 CNR⁃1辐射表测量(0.1 W·m-2);空气温湿
度测量采用 HMP⁃45C 空气温湿仪 ( Vaisala, Fin⁃
land) ( 0. 001K);大气压采用 CS105 气压表测量
(1 hPa);水汽压是利用大气压以及温度计算得到.
记录时间间隔为 30 min,剔除明显的测量异常值.然
后按照 cf值大小将所有数据分为晴天( cf <0.05)数
据和云天(cf≥0.05)数据[27],最后分别对晴天与云
天情况下的大气长波辐射进行计算分析.
采用软件 Matlab 7.10 进行误差计算、拟合、线
性回归分析;采用 Origin 8.5软件作图.
2  结果与分析
2􀆰 1  晴天参数化模型的适用性
由于经验模型都有局地性,因此在将其应用于
局部地区时需对参数进行校正.本文利用站点的观
测数据,对文中提到的 8 个常用的晴天参数化模型
的计算结果进行比较.由表 2 可以看出,除了 Idso
等[13]提出模型和 Abramowitz 等[26]提出模型的模拟
结果比实测值高外,其他模型都低于实测值,这种差
异主要是由发射率的不同参数化方案造成的;
Brunt[24]提出模型的平均偏差(BIAS)、均方根误差
(RMSE)、平均绝对误差(MAE)、平均相对误差率
(PMRE)值最小,其次是 Satterlund[28]提出模型、Idso
等[13]提出模型,但 Satterlund[28]提出模型模拟值与
实测值的相关性好于其他两个模型(R2 = 0􀆰 897).
Prata[14]提出模型和 Brutsaert[15]提出模型的 BIAS、
RMSE、MAE和 PMRE值最大,但它们的线性相关性
相对较好,其中,Brutsaert[16]提出模型在所有模型中
的线性相关性最好(R2 = 0.905).
    为了进一步比较各个模型的计算结果,对模型
的模拟值与观测值进行回归分析.由图 1 可以看出,
Satterlund[28]提出模型的各模拟点相对集中,回归线
接近直线,即模拟值与观测值的一致性较好,其
RMSE、MAE都小于 30 W·m-2,PMRE为 10.1%,说
明模拟效果最好. Idso 等[13]提出模型、Kruk 等[29]提
出模型和 Brunt[24]提出模型的模拟结果偏离观测值
较大;Prata[14]提出模型、Dilley等[25]提出模型的模
9942期                        李伟斌等: 大气向下长波辐射参数化模型在长白山地区的适用性           
表 2  晴天参数化模型的统计结果
Table 2  Summary statistics of clear⁃sky parameterized models
模型
Model
BIAS
(W·m-2)
RMSE
(W·m-2)
MAE
(W·m-2)
PMRE
(%)
线性回归分析 Linear regression analysis
a b R2
文献
Reference
Brunt -18.84 28.17 20.35 8.0 0.69 52.18 0.867 [23]
Idso, et al. 24.99 32.00 28.78 14.0 0.81 68.72 0.849 [24]
Brutsaert -69.40 71.35 69.40 30.9 0.81 -25.71 0.905 [15]
Satterlund -23.34 28.55 23.34 10.1 0.87 5.54 0.897 [28]
Prata -72.58 76.83 72.58 30.9 0.55 29.71 0.878 [13]
Dilley, et al. -45.92 52.91 45.92 18.6 0.54 59.44 0.848 [25]
Kruk, et al. -48.60 51.96 48.60 22.7 1.08 -67.20 0.904 [29]
Abramowitz, et al. 31.66 40.43 35.68 18.1 0.58 127.10 0.835 [26]
BIAS: 偏差 Bias; RMSE: 均方根误差 Root mean square error; MAE: 平均绝对误差 Mean average error; PMRE: 平均相对误差 Root mean square
error; a: 斜率 Slope; b: 截距 Intercept. 下同 The same below.
图 1  晴天各模型模拟下行长波辐射值与观测值的比较
Fig.1  Comparisons between the measured and modeled downward longwave radiation (DLR) in clear⁃sky conditions.
拟值与观测值的一致性最差.
2􀆰 2  晴天参数化模型的系数校正及敏感性分析
为了将模拟误差降到最低,采用整个观测数据
对模拟效果相对较好的 Satterlund[28] 模型、 Kruk
等[29]模型和 Brunt[24]模型进行系数拟合校正后,模
型的模拟性能都有很大提高,模拟值与实测值间的
一致性也有很大改善(表 3).其中,Satterlund[28]提
出模型的拟合效果最好,Brunt[24]模型和 Kruk 等[29]
模型次之,经过系数校正,这 3个模型的 RMSE值都
<20 W·m-2;BIAS 的绝对值均<10 W·m-2,R2均
>0􀆰 9.Kruk等[29]模型的系数变化较敏感,在校正前
后变化较大;Satterlund[28]模型和 Brunt[24]模型中系
数 a变化很小,系数 b变化相对较大.
通过对拟合后的Satterlund[28]模型进行敏感性
表 3  晴天参数化模型系数校正后的模拟结果
Table 3  Summary of locally calibrated coefficients for the better clear⁃sky parameterized models
模型
Model
校正后系数
Calibrated coefficient
BIAS
(W·m-2)
RMSE
(W·m-2)
MAE
(W·m-2)
PMRE
(%)
线性回归分析
Linear regression analysis
a b R2
文献
Reference
Brunt a= 0.599,b= 0.155 0.90 17.13 10.32 4.7 0.91 21.34 0.890 [23]
Satterlund a= 1.188,b= 1.734 -6.33 18.08 9.68 4.2 0.99 -5.72 0.900 [28]
Kruk, et al. a= 1.268,b= 0.095 -8.91 18.89 10.42 4.6 0.97 -1.63 0.900 [29]
Rlc:晴天大气长波辐射 Longwave radiation in clear⁃sky (W·m-2); Rl:云天大气长波辐射 Longwave radiation in cloudy⁃sky (W·m-2); Ta:近地
面大气温度 Surface air temperature (K); ea: 近地面水汽压 Surface air vapor (hPa); cf: 云量 Fractional cloud cover.
005 应  用  生  态  学  报                                      26卷
图 2  晴天大气发射率随水汽压(ea)和大气温度(T)的变化
Fig.2  Variations of atmospheric emissivity in clear⁃sky changed
with water vapor pressure and air temperature.
分析,可以理解输入不同变量后晴天长波辐射变化
程度的大小.研究期间,该地大气温度与水汽压值分
别为 247.58~301.84 K和 0.04~2.01 hPa;平均温度
和平均水汽压值分别为 279.97 K 和 0.59 hPa,标准
差分别为 12.28 K和 0.51 hPa.
在整个观测期间水汽压(或大气温度)的变化
范围内,使大气温度(或水汽压)在其标准差范围内
变化,得到晴天大气发射率随水汽压和大气温度的
变化规律.由图 2可见,大气发射率随着温度的增加
而降低,降低幅度较小;随着水汽压的增大而增加,
增加趋势呈抛物线形.说明大气发射率对水汽压变
化的响应比对大气温度变化的响应更敏感,且当水
汽压越小时越敏感,主要原因在于空气中水汽的吸
收与发射作用,这与复杂的空气成分中水汽对长波
辐射影响最大[30]这一结果相吻合.
2􀆰 3  云天参数化模型的适用性
对文中 8个云天参数化模型进行计算比较.其
中,模型中的晴天大气长波辐射值由系数校正过的
Satterlund[28]模型计算而得.从表 4 可以看出,除了
没有考虑云量的 Abramowitz等[26]模型之外,其他模
型的性能都较好,说明 cf在云天长波辐射模拟中起
着非常重要的作用.模拟效果最好的是 Jacobs[31]
模型,其 BIAS 只有 0. 38 W·m-2,RMSE 为 29􀆰 29
W·m-2,在研究期内其模拟值的线性回归线与直线
y= x几乎重合(图 3),即模拟值与实测值具有非常
好的一致性;其次是 Duarte 等[27]模型,其 BIAS、
RMSE都稍大,分别为 5.52 和 29.88 W·m-2,但其
模拟值与实测值的一致性也较好;其他模型的模拟
结果距观测值偏离均较大.
2􀆰 4  大气长波辐射的平均日变化
选取未用于模型校正的 2004 年 8 月观测数据
并求其日平均,然后将模拟效果最好的晴天 (Satter⁃
lun[28])与云天(Jacobs[31] )大气长波辐射模型相结
合来模拟该月的日平均长波辐射,并与观测值进行
比较.由图 4可以看出,平均长波辐射虽然没有明显
的日变化特征,但也可以看到大气长波辐射日变化
的基本趋势,即在早晨时为最低值,随着大气温度的
升高迅速增大,于午间达到最大,随后开始下降,这
主要是由于长白山地区云量多这一特点所致.总体
来看,模型模拟值与观测值的变化趋势基本一致,但
由于该区域云量变化较快,所以在云量发生突变的
节点上模拟效果不太理想,这将是今后研究该类区
域大气长波辐射需要关注的问题.
表 4  云天参数化模型的统计结果
Table 4  Summary statistics of cloudy⁃sky parameterized models
模型
Model
BIAS
(W·m-2)
RMSE
(W·m-2)
MAE
(W·m-2)
PMRE
(%)
线性回归分析
Linear regression analysis
a b R2
文献
Reference
Maykut, et al. -17.90 34.06 27.35 10.0 0.90 9.58 0.826 [32]
Jacobs, et al. 0.38 29.29 22.50 8.7 0.97 7.67 0.838 [31]
Sugita, et al. -30.92 42.84 33.67 11.7 0.81 21.28 0.812 [33]
Konzelmann, et al. -14.48 36.74 28.77 10.8 0.85 28.71 0.764 [34]
Crawford, et al. 22.35 40.99 29.17 11.9 0.87 58.65 0.762 [3]
Duarte, et al. 5.52 29.88 23.99 9.2 0.97 13.23 0.837 [27]
Duarte, et al. 30.59 45.25 34.67 14.2 0.84 75.59 0.768 [27]
Abramowitz, et al. -8.12 45.20 38.29 14.5 0.42 155.12 0.672 [26]
1052期                        李伟斌等: 大气向下长波辐射参数化模型在长白山地区的适用性           
图 3  云天各模型模拟下行长波辐射值与观测值的比较
Fig.3  Comparisons between the measured and modeled downward longwave radiation (DRL) in cloudy⁃sky conditions.
图 4  2004年 8月日平均大气长波辐射
Fig. 4   Mean downward longwave radiation ( DLR) during
August 2004.
Ⅰ: 观测值 Measured value; Ⅱ: 模拟值 Modeled value.
3  结    论
本文利用中国科学院长白山森林生态系统定位
站的实地观测数据对现有被广泛应用的 8个晴天大
气长波辐射参数化模型和 8个云天大气长波辐射参
数化模型进行比较与适用性分析,得到了适用于长
白山地区的长波辐射模型与系数,并对其进行敏感
性分析以及模拟值与观测值的日变化趋势分析.
通过对晴天的 8 个参数化模型对比分析发现,
在系数校正前,Brunt[24]模型的 BIAS、RMSE、MAE
和 PMRE 值最小,其次是 Satterlund[28] 模型、 Idso
等[13]模型,但 Satterlund[28]模型的线性相关性较其
他两个模型好(R2 = 0.897)、且其模拟值与实测值的
一致性更好,所以综合评价 Satterlund[28]模型具有
较好的适宜性.对模型的系数校正以后,其模拟效果
有了进一步提高,其中,Satterlund[28]模型的拟合效
果最好,Brunt[24]模型和 Kruk 等[29]模型次之,它们
的 RMSE值均<20 W·m-2;BIAS 的绝对值也均<10
W·m-2,R2也均>0.9.Kruk 等[29]模型的系数变化较
敏感,校正前后变化较大; Satterlund[28] 模型和
Brunt[24]模型中系数 a 变化很小、系数 b 变化相对
较大.
对系数校正后的 Satterlund[28]模型进行敏感性
分析发现,大气发射率对温度的敏感性较小;对水汽
压的敏感性较大,且在水汽压越小时越敏感.
在有云情况下,Jacobs[31]模型和 Duarte 等[27]
模型的模拟效果非常好,可直接在长白山地区应
用而不需系数校正,其 BIAS 分别为 0. 38 和 5. 52
W·m-2;RMSE分别为 29.29和 29.88 W·m-2;研究
期间,模拟值与实测值的一致性也非常好,R2分别为
0􀆰 838与 0.837.其他几个模型的模拟效果相差不大,
而没有考虑云量变化的 Abramowitz等[26]模型的性能
最差,其模拟值与实测值间的差异较明显,说明云量
在云天大气长波辐射模拟中起着非常重要的作用.
大气长波辐射虽然没有明显的日变化特征,但
模拟值与观测值的变化趋势基本一致;只是由于该
区域云量较多且变化较快,所以在云量发生变化的
节点上模拟效果不太理想,这将是今后研究亟待解
决的一个问题.
基于本文的分析,Satterlund[28]模型在晴天具有
较好的模拟效果,经系数校正后可用于相似于长白
205 应  用  生  态  学  报                                      26卷
山地区的区域.Jacobs[31]模型和 Duarte 等[27]模型在
云天有很好的适宜性,可直接用于该类区域.
参考文献
[1]  He Q (何   清), Miao Q⁃L (缪启龙), Li S (李
帅), et al. Characteristic of long wave radiation over the
Taklimakan Desert Hinterland. Plateau Meteorology (高
原气象), 2009, 28(3): 642-646 (in Chinese)
[2]   Sridhar V, Elliott RL. On the development of a simple
downwelling longwave radiation scheme. Agricultural and
Forest Meteorology, 2002, 112: 237-243
[3]  Crawford TM, Duchon CE. An improved parameteriza⁃
tion for estimating effective atmospheric emissivity for
use in calculating daytime downwelling longwave radia⁃
tion. Journal of Applied Meteorology, 1999, 38: 474 -
480
[4]   Choi M, Jacobs JM, Kustas WP. Assessment of clear
and cloudy sky parameterizations for daily downwelling
longwave radiation over different land surfaces in Flori⁃
da, USA. Geophysical Research Letters, 2008, 35:
L20402
[5]  Li H⁃D (李辉东), Guan D⁃X (关德新), Jin C⁃J (金
昌杰), et al. Measurement and estimation methods and
research progress of snow evaporation in forests. Chinese
Journal of Applied Ecology (应用生态学报), 2011, 24
(12): 3603-3609 (in Chinese)
[6]  Wang KC, Liang SL. Global atmospheric downward
longwave radiation over land surface under all⁃sky condi⁃
tions from 1973 to 2008. Journal of Geophysical Re⁃
search: Atmospheres, 2009, 114: D19101
[7]  Kjaersgaard JH, Plauborg FL, Hansen S. Comparison of
models for calculating daytime long⁃wave irradiance
using long term data set. Agricultural and Forest Meteo⁃
rology, 2007, 143: 49-63
[8]  Kneizys FX, Shettle EP, Abreu LW, et al. User’ s
Guide to LOWTRAN7, AFGL⁃TR⁃88⁃0177. Environ⁃
mental Research Papers No. 1010, Air Force Geopysics
Laboratory, Hanscom AFB, MF, 1988
[9]  Snell HE, Anderson GP, Wang J, et al. Validation of
FASE ( FASCODE for the environment) and MODT⁃
RAN3: Updates and comparisons with clear⁃sky
measurements. Proceedings of SPIE⁃The International
Society for, Optical Engineering, Paris, France, 1995:
194-204
[10]  Key JR, Schweiger AJ. Tools for atmospheric radiative
transfer: STREAMER and FLUXNET. Computers & Geo⁃
sciences, 1998, 24: 443-451
[11]  Francis JA. A method to derive downwelling longwave
fluxes at the Arctic surface from TIROS operational verti⁃
cal sounder data. Journal of Geophysical Research:
Atmospheres, 1997, 102: 1795-1806
[12]  Wu R⁃S (武荣盛), Ma Y⁃M (马耀明). Comparative
analyses on radiation characteristics in different areas
over the Tibetan Plateau. Plateau Meteorology (高原气
象), 2010, 29(2): 251-259 (in Chinese)
[13]  Idso SB, Jackson RD. Thermal radiation from the atmos⁃
phere. Journal of Geophysical Research, 1969, 74:
5397-5403
[14]  Prata A. A new long⁃wave formula for estimating down⁃
ward clear⁃sky radiation at the surface. Quarterly Journal
of the Royal Meteorological Society, 1996, 122: 1127-
1151
[15]  Jin X, Barber D, Papakyriakou T. A new clear⁃sky
downward longwave radiative flux parameterization for
Arctic areas based on rawinsonde data. Journal of Geo⁃
physical Research: Atmospheres, 2006, 111: 1-7
[16]  Brutsaert W. On a derivable formula for long⁃wave radia⁃
tion from clear skies. Water Resources Research, 1975,
11: 742-744
[17]  Huang M⁃F (黄妙芬), Xing X⁃F (邢旭峰), Zhu Q⁃J
(朱启疆), et al. Improvement of models for the estima⁃
tion of downward atmospheric long⁃wave radiation for
retrieving net radiation fluxes based on remote sensing.
Geographical Research (地理研究), 2005, 24 ( 5):
757-765 (in Chinese)
[18]  Yu S⁃S (余珊珊), Xin X⁃Z (辛晓洲), Liu Q⁃H (柳
钦火). Comparison of atmospheric downward longwave
radiation parameterizations. Advances in Earth Science
(地球科学进展), 2011, 26(7): 751-762 ( in Chi⁃
nese)
[19]  Jin C⁃J (金昌杰), Guan D⁃X (关德新), Zhu T⁃Y
(朱廷耀). Spectral characteristics of solar radiation in
broadleaved Korean pine forest in Changbai Mountain.
Chinese Journal of Applied Ecology (应用生态学报),
2000, 11(1): 19-21 (in Chinese)
[20]  Zhang X⁃J (张新建), Yuan F⁃H (袁凤辉), Chen N⁃
N (陈妮娜), et al. Energy balance and evapotranspira⁃
tion in broad⁃leaved Korean pine forest in Changbai
Mountains. Chinese Journal of Applied Ecology (应用生
态学报), 2011, 22(2): 607-613 (in Chinese)
[21]  Zhang M (张   弥), Guan D⁃X (关德新), Han S⁃J
(韩士杰), et al. Climatic dynamics of broadleaved Ko⁃
rean pine forest in Changbai Mountain during the last 22
years. Chinese Journal of Ecology (生态学杂志),
2005, 24(9): 1007-1012 (in Chinese)
[22]  Schmetz J. Towards a surface radiation climatology: Re⁃
trieval of downward irradiances from satellites. Atmos⁃
pheric Research, 1989, 23: 287-321
[23]  Ji C⁃R (吉春容), Zou C (邹  陈), Fan Z⁃A (范子
昂), et al. Radiation characteristics of Picea schrenki⁃
ana forest area in middle Tianshan Mountains of China.
Journal of Arid Land Resources and Environment (干旱
区资源与环境), 2011, 25(12): 82-88 (in Chinese)
[24]  Brunt D. Notes on radiation in the atmosphere. Quarterly
Journal of the Royal Meteorological Society, 1932, 58:
389-420
[25]  Dilley A, O’ Brien D. Estimating downward clear sky
long⁃wave irradiance at the surface from screen tempera⁃
ture and precipitable water. Quarterly Journal of the
Royal Meteorological Society, 1998, 124: 1391-1401
[26]  Abramowitz G, Pouyanné L, Ajami H. On the informa⁃
tion content of surface meteorology for downward atmos⁃
pheric long⁃wave radiation synthesis. Geophysical Re⁃
3052期                        李伟斌等: 大气向下长波辐射参数化模型在长白山地区的适用性           
search Letters, 2012, 39: L04808
[27]  Duarte HF, Dias NL, Maggiotto SR. Assessing daytime
downward longwave radiation estimates for clear and
cloudy skies in Southern Brazil. Agricultural and Forest
Meteorology, 2006, 139: 171-181
[28]  Satterlund DR. An improved equation for estimating
long⁃wave radiation from the atmosphere. Water Re⁃
sources Research, 1979, 15: 1649-1650
[29]  Kruk NS, Vendrame ÍF, Da Rocha HR, et al. Down⁃
ward longwave radiation estimates for clear and all⁃sky
conditions in the Sertãozinho region of Säo Paulo, Brazil.
Theoretical and Applied Climatology, 2010, 99: 115 -
123
[30]   Viúdez⁃Mora A, Calb􀆩 J, Gonz􀅡lez J, et al. Modeling
atmospheric longwave radiation at the surface under
cloudless skies. Journal of Geophysical Research: Atmos⁃
pheres, 2009, 114: 1-12
[31]  Jacobs JD. Radiation climate of Broughton Island. Ener⁃
gy budget studies in relation to fast⁃ice breakup proces⁃
ses in Davis Strait Occas. Paper No. 26, Institute of
Arctic and Alpine Research, University of Colorado,
Boulder, United States, 1978: 105-120
[32]  Maykut G, Church P. Radiation climate of Barrow,
Alaska, 1962-1966. Journal of Applied Meteorology,
1973, 12: 620-628
[33]  Sugita M, Brutsaert W. Cloud effect in the estimation of
instantaneous downward longwave radiation. Water Re⁃
sources Research, 1993, 29: 599-605
[34]  Konzelmann T, Van de Wal RSW, Greuell W, et al.
Parameterization of global and longwave incoming radia⁃
tion for the Greenland ice sheet. Global and Planetary
Change, 1994, 9: 143-164
作者简介  李伟斌,男,1990 年生,硕士研究生.主要从事森
林生态系统水文气象功能与全球变化研究. E⁃mail: lwbin12
@ 126.com
责任编辑  杨  弘
405 应  用  生  态  学  报                                      26卷