散射辐射的准确估算对于评价其对陆地生态系统碳交换的影响具有重要意义.基于我国中亚热带江西千烟洲气象观测场2012年3月1日—2013年2月28日散射辐射实际观测数据对目前常用的5个散射辐射分解模型(Reindl-1、Reindl-2、Reindl-3、Boland、BRL)的模拟结果进行验证.结果表明: 在30 min尺度上,虽然5个模型在总体上都可以较好地模拟该地区的散射辐射,但模型模拟效果随着晴空指数(kt)的升高而显著降低.特别是当kt>0.75时,各模型已无法模拟该地区散射辐射.从散射辐射季节变化的模拟来看,5个模型能够很好地模拟大多数月份的散射辐射.5个模型年尺度散射辐射模拟值与观测值的相对偏差最高为7.1%(BRL),最低为0.04%(Reindl-1),平均为3.6%.在全年辐射最强、温度最高和降水偏少的夏季,5个模型的模拟值均出现了过高估计.以7月为例,散射辐射被高估14.5%~28.2%,平均高估21.2%.这可能与高kt条件下散射辐射的估算方法有关,这种不确定性需要在模型应用中做进一步深入评价.根据验证结果并考虑模拟精度和输入变量的易获取性,5个模型的模拟效果依次为BRL>Reindl-3>Reindl-2>Reindl-1>Boland.
Accurate estimation of diffuse radiation is of great significance for evaluating its effect on terrestrial ecosystem carbon exchange. Based on the observed diffuse radiation data in the meteorological observation field in midsubtropical Qianyanzhou, Jiangxi, China from March 1, 2012 to February 28, 2013, the simulated results of five widely used diffuse radiation decomposition models (Reindl-1, Reindl-2, Reindl-3, Boland, BRL) were validated. The results indicated that, on the 30 min scale, all of the five models could well simulate the diffuse radiation of this area overall. But the effect of models decreased significantly with the rising of clearness index (kt). Especially when kt>0.75, each model was unable to simulate diffuse radiation in the region. Regarding the simulation of seasonal change of diffuse radiation, the five models could simulate diffuse radiation well in most months. Relative deviation between simulated and observed values of yearly diffuse radiation of five models had a maximum of 7.1% (BRL), a minimum of 0.04% (Reindl-1), and an average of 3.6%. The simulated values of the five models appeared to be overestimated in the summer when radiation was strongest, temperature was highest, and precipitation was relatively low. For example, in July, the diffuse radiation was overestimated by 14.5%-28.2%, 21.2% on average. This was primarily due to the method of estimating diffuse radiation under the condition of high kt. The uncertainty requires further evaluation in the model application. Considering the results of validation, simulation precision and the accessibility of input variables, the order of the simulation performance of five models was BRL>Reindl-3>Reindl-2>Reindl-1>Boland.
全 文 :江西千烟洲地区散射辐射的模拟与验证∗
韩佳音1,2 李胜功1 张雷明1∗∗ 温学发1 李庆康1 王辉民1
( 1中国科学院地理科学与资源研究所生态系统网络观测与模拟重点实验室, 北京 100101; 2中国科学院大学, 北京 100049)
摘 要 散射辐射的准确估算对于评价其对陆地生态系统碳交换的影响具有重要意义.基于
我国中亚热带江西千烟洲气象观测场 2012 年 3 月 1 日—2013 年 2 月 28 日散射辐射实际观
测数据对目前常用的 5个散射辐射分解模型(Reindl⁃1、Reindl⁃2、Reindl⁃3、Boland、BRL)的模
拟结果进行验证.结果表明: 在 30 min尺度上,虽然 5个模型在总体上都可以较好地模拟该地
区的散射辐射,但模型模拟效果随着晴空指数(kt)的升高而显著降低.特别是当 kt >0.75 时,
各模型已无法模拟该地区散射辐射.从散射辐射季节变化的模拟来看,5 个模型能够很好地模
拟大多数月份的散射辐射.5个模型年尺度散射辐射模拟值与观测值的相对偏差最高为 7.1%
(BRL),最低为 0.04%(Reindl⁃1),平均为 3.6%.在全年辐射最强、温度最高和降水偏少的夏
季,5个模型的模拟值均出现了过高估计.以 7 月为例,散射辐射被高估 14.5% ~ 28.2%,平均
高估 21.2%.这可能与高 kt条件下散射辐射的估算方法有关,这种不确定性需要在模型应用中
做进一步深入评价.根据验证结果并考虑模拟精度和输入变量的易获取性,5 个模型的模拟效
果依次为 BRL>Reindl⁃3>Reindl⁃2>Reindl⁃1>Boland.
关键词 总辐射; 散射辐射; 晴空指数; 散射分数; 中亚热带; 季节性干旱
文章编号 1001-9332(2015)10-2991-09 中图分类号 P414.5, P422.3 文献标识码 A
Simulation and validation of diffuse radiation in Qianyanzhou area, Jiangxi, China. HAN Jia⁃
yin1,2, LI Sheng⁃gong1, ZHANG Lei⁃ming1, WEN Xue⁃fa1, LI Qing⁃kang1, WANG Hui⁃min1 (1Key
Laboratory of Ecosystem Network Observation and Modeling, Institute of Geographic Sciences and Natu⁃
ral Resources Research, Chinese Academy of Sciences, Beijing 100101, China; 2University of Chinese
Academy of Sciences, Beijing 100049, China) . ⁃Chin. J. Appl. Ecol., 2015, 26(10): 2991-2999.
Abstract: Accurate estimation of diffuse radiation is of great significance for evaluating its effect on
terrestrial ecosystem carbon exchange. Based on the observed diffuse radiation data in the meteoro⁃
logical observation field in mid⁃subtropical Qianyanzhou, Jiangxi, China from March 1, 2012 to
February 28, 2013, the simulated results of five widely used diffuse radiation decomposition models
(Reindl⁃1, Reindl⁃2, Reindl⁃3, Boland, BRL) were validated. The results indicated that, on the
30 min scale, all of the five models could well simulate the diffuse radiation of this area overall. But
the effect of models decreased significantly with the rising of clearness index (kt). Especially when
kt>0.75, each model was unable to simulate diffuse radiation in the region. Regarding the simula⁃
tion of seasonal change of diffuse radiation, the five models could simulate diffuse radiation well in
most months. Relative deviation between simulated and observed values of yearly diffuse radiation of
five models had a maximum of 7.1% (BRL), a minimum of 0.04% (Reindl⁃1), and an average of
3.6%. The simulated values of the five models appeared to be overestimated in the summer when ra⁃
diation was strongest, temperature was highest, and precipitation was relatively low. For example,
in July, the diffuse radiation was overestimated by 14.5%-28.2%, 21.2% on average. This was pri⁃
marily due to the method of estimating diffuse radiation under the condition of high kt . The uncer⁃
tainty requires further evaluation in the model application. Considering the results of validation,
simulation precision and the accessibility of input variables, the order of the simulation performance
of five models was BRL>Reindl⁃3>Reindl⁃2>Reindl⁃1>Boland.
Key words: global radiation; diffuse radiation; clearness index; diffuse fraction; mid⁃subtropical
area; seasonal drought.
∗国家自然科学基金项目(31170422)和国家自然科学基金委员会重大项目(31290221)资助.
∗∗通讯作者. E⁃mail: zhanglm@ igsnrr.ac.cn
2015⁃01⁃14收稿,2015⁃06⁃11接受.
应 用 生 态 学 报 2015年 10月 第 26卷 第 10期
Chinese Journal of Applied Ecology, Oct. 2015, 26(10): 2991-2999
绿色植物进行光合作用离不开太阳辐射,太阳
辐射强度和组成的动态变化对其有重要影响[1-2] .太
阳辐射经过地球大气时,其中一部分会被吸收、反射
或散射,到达地面的太阳辐射可以分为直接辐射和
散射辐射[3] .近年来散射辐射在陆地生态系统生产
力中的作用逐步引起人们的关注,一些研究发现,增
加散射辐射对增加陆地碳汇有积极影响[4-6] .但也有
研究表明,生态系统在散射辐射比例提高的情况下
CO2吸收量减少[7] .增加散射辐射对陆地生态系统净
碳吸收的影响目前还不明确,因此阐明散射辐射对
陆地生态系统碳汇的作用在未来的研究中有重要
意义.
总辐射作为一个重要的环境要素在全球范围内
开展了长期、连续的观测.相对而言,散射辐射的观
测站点少且观测时间不长[8-9] .为了准确评价散射辐
射在不同时间尺度上的变化特征及其对陆地生态系
统碳水循环过程的影响,需要长期并且连续的散射
辐射数据.目前,相关研究基本上采用模型模拟的方
法对散射辐射进行估算.常用的散射辐射模型大体
可以按模型结构分为参数模型和分解模型两大
类[10],前者一般需要输入大气状况的详细信息,如
云的类型、数量、分布、日照百分率、大气浑浊度、降
水量等;后者通常只需要输入总辐射、太阳高度角等
较易获得的观测数据或理论计算数据.分解模型相
比参数模型输入的变量少,所以被较广泛地应用并
且其估算结果较为可靠.例如,在模型应用方面的文
献中,Gu等[11]基于分解模型模拟了加拿大地区的
散射辐射,发现散射辐射的增加可以显著提高北方
森林生态系统的 CO2吸收能力;Ren 等[12]利用分解
模型分析了中国地区 1981—2010 年散射辐射的时
空变异性;Ridley 等[13]验证结果表明,分解模型中
的 Boland模型[14-15]在欧洲、非洲、澳洲和亚洲都表
现出来良好的性能,并将其发展为一个全球性的散
射辐射计算模型.
江西千烟洲地区位于我国的中亚热带,这里的
人工针叶林是我国南方人工林的代表性区域,同时
也是我国重要的陆地碳汇分布区[16] .何学兆等[17]基
于 MAESTRA模型分析了光合有效辐射总量及其散
射辐射比例变化对该地区森林总初级生产力
(GPP)的影响.李登秋等[18]基于 BEPS 模型分析了
太阳辐射变化对该地区 GPP 的影响.但这些研究均
是利用散射辐射模型模拟得到的散射辐射变化特
征,缺少实际观测数据的验证.因此,在该地区开展
散射辐射与其他要素的平行观测不仅有利于准确刻
画散射辐射的动态变化特征,也将为进一步研究散
射辐射对森林生态系统碳水交换过程提供数据支持
和科学依据.
本研究利用千烟洲地区散射辐射实际观测数
据,对目前常用的散射辐射模型进行验证,并通过模
型之间的对比分析评价其在该地区的适用情况,筛
选确定适合于该地区的散射辐射模型,可为准确评
价该地区长期散射辐射变化对中亚热带人工针叶林
生态系统的碳水循环过程提供重要的科学依据.
1 研究地区与研究方法
1 1 研究地点概况
本研究在中国科学院千烟洲红壤丘陵综合开发
试验站 (以下简称千烟洲站) 进行. 千烟洲站
(26°44′48″ N,115°04′13″ E)位于江西省吉安市泰
和县灌溪镇境内.站内海拔多在 100 m 左右,相对高
度 20~50 m,属于典型的红壤丘陵区.研究区属典型
亚热带季风气候,年平均气温 17.8 ℃, 0 ℃以上活
动积温 6523 ℃,10 ℃以上活动积温 6015 ℃;多年
平均降水量 1360 mm;年日照时数 1785 h,年日照百
分率 44%,太阳年总辐射量 4349 MJ·m-2;无霜期
290 d.夏季易受到副热带高压的控制而出现高温少
雨的季节性干旱.
主要土壤类型有红壤、水稻土、潮土、草甸土等,
成土母质多为红色砂岩、砂砾岩或泥岩,以及河流冲
积物.植被属中亚热带常绿阔叶林带,但原生植被已
破坏殆尽,现主要为人工林或草、灌丛次生植被,主
要树种有马尾松(Pinus massoniana)、湿地松(Pinus
elliottii)、 杉 木 ( Cunninghamia lanceolata )、 木 荷
(Schima crenata)、柑橘(Citrus spp.)等[19] .
1 2 观测仪器
本研究的散射辐射观测在千烟洲站的综合气象
观测场内进行.观测区域内场地开阔,下垫面平坦、
均质,周围没有树木或建筑物遮挡,利用 SPN1
(Delta⁃T Devices Ltd, 英国)和 RSR3(华益瑞, 中
国)开展散射辐射的长期、连续观测.观测仪器安装
在距离地面 1.5 m高度.测定期间,人工检查和清洁
传感器表面以减少测定误差.
SPN1仪器的测定原理是在测定探头上安装了
多孔的半球形遮光罩,开孔和非开孔的表面积相同,
每个开孔下均有辐射传感器.利用散射辐射无方向
性传播的特点,观测时测定到的最大和最小辐射即
分别为总辐射和散射辐射,可同时测定这两个要
素[20] .
2992 应 用 生 态 学 报 26卷
RSR3由辐射测定探头 LI200X(Campbell Scien⁃
tific, Inc., 美国)和自动旋转遮光臂两部分组成.在
测定时,遮光臂每 10 s 旋转一次,每个测定周期内
选取最大值和最小值分别作为总辐射和散射辐射.
相对上述两种仪器而言, CMP11 ( Kipp &
Zonen, 荷兰)以其精度高和稳定性好而得到广泛应
用[21] .之前的研究结果表明,RSR3 测定的散射辐射
与 CMP11 和 SPN1 的测定结果之间具有很好的一
致性,但其观测的总辐射存在系统偏差[22] .图 1 比
较了同一观测期间 SPN1、RSR3与 CMP11平行观测
的散射辐射,从中可以看出,SPN1 测定的散射辐射
相对 CMP11测定的结果有较好的一致性,符合测定
标准可以使用.因此,在后续的分析和计算中均采用
SPN1的测定结果.
1 3 数据选取
本研究选取了 2012 年 3 月 1 日—2013 年 2 月
28日 SPN1 的总辐射和散射辐射观测数据,其中
2012年 11月 12日 8:25—8:28 的 4 条分钟数据缺
失,2012年 11 月 26 日 14:24—15:16 的 53 条分钟
数据缺失,对以上缺失数据运用平均变化法(滑动
窗口为 15 min)进行插补,可得到完整的一年数据.
为处理方便,本研究取 30 min 的平均值以便于与涡
度相关通量塔上同一时间的温度、相对湿度和降雨
量 30 min尺度数据对应.
1 4 质量控制
对观测数据进行了质量控制[23] .因为太阳高度
角在日出和日落时较小,此时仪器容易出现余弦误
差而使测量值偏大.通常认为剔除太阳高度角在
6° ~10°以下的值都为可接受范围[24],所以在后续处
图 1 SPN1、RSR3与 CMP11测定的散射辐射结果对比
Fig. 1 Comparison of diffuse radiation measured by SPN1,
RSR3 and CMP11.
虚线为 1 ∶ 1线,实线为拟合的趋势线 The dash line was 1:1 line, and
the solid line was fitting trend line.
理时将太阳高度角小于 7°时的测量值全部剔除.
降水也是影响散射辐射仪器准确测量的一个因
素,根据降雨量数据将降水前 1 h和降水后 2 h数据
剔除.
总辐射、散射辐射、地外大气上界辐射分别记作
IG、ID、IE .理论上 IG不能超过 IE,ID不能超过 IG .在阴
天和晴朗条件下,根据 Reindl等[25]提出的两个约束
条件将以下数据剔除: IG / IE < 0. 2 且 ID / IG < 0. 9;
IG / IE>0.6且 ID / IG>0.8
对所有 17520 条数据进行质量控制后得到
5859条有效数据.
1 5 模型介绍
在分解模型中有两个较为重要的物理量:晴空
指数(kt)和散射分数(kd) [26] .kt = IG / IE,其范围是
0~1.kd = ID / IG,其范围是 0~ 1.前者表示天空的晴朗
和通透程度,后者反映散射辐射在总辐射中的变化,
二者之间存在一定的关联关系.如何准确构建基于
kt估算 kd的函数关系,进而通过总辐射计算散射辐
射是散射辐射模型模拟精确的关键.
1 5 1 Reindl模型 Reindl 模型因形式简单而运用
较为广泛.Reindl 等[25]利用北美和欧洲的 5 个站点
观测得到的 22000条小时尺度数据,从 28 个潜在预
测变量中总结出 4个对散射辐射有显著影响的预测
变量,具体表现形式有 3 种,分别记作 Reindl⁃1、Re⁃
indl⁃2和 Reindl⁃3:
Reindl⁃1:
kd = 1.02-0.248kt (kt≤0.3)
kd = 1.45-1.67kt (0.3<kt<0.78)
kd = 0.147 (kt≥0.78)
Reindl⁃2:
kd = 1.02-0.254kt+0.0123sinβ (kt≤0.3)
kd = 1.4-1.749kt+0.177sinβ (0.3<kt<0.78)
kd = 0.486kt-0.182sinβ (kt≥0.78)
Reindl⁃3:
kd = 1 - 0. 232kt + 0. 0239sinβ - 6. 82 × 10
-4 Ta +
0 0195Φ (kt≤0.3)
kd = 1. 329 - 1. 761kt + 0. 267sinβ- 3. 57 × 10
-3 Ta +
0 106Φ (0.3<kt<0.78)
kd = 0.426kt -0.256sinβ+3.49×10
-3Ta +0.0734Φ
(kt≥0.78)
式中:β 为太阳高度角(°);Ta为环境温度(℃);Φ
为相对湿度(%).
1 5 2 Boland 模型 Boland 等[14]利用 Australia 的
Victoria站点的数据构建了一个简单的指数模型,其
399210期 韩佳音等: 江西千烟洲地区散射辐射的模拟与验证
具体形式如下:
kd =
1
1+exp[7.997(kt-0.586]
1 5 3 BRL 模型 Ridley 等[13]在 Boland 模型的基
础上又引入了其他的输入变量,具体形式如下:
kd = 1 / [ 1 + exp ( - 5. 38 + 6. 63kt + 0. 006AST -
0 007β+1.75K t+1.31ψ)]
ψ=
(kt+1+kt-1) / 2
kt+1
kt-1
ì
î
í
ïï
ïï
日出<t<日落
t=日出
t=日落
式中:AST为当地太阳时,其计算方法见文献[27];
K t为日尺度晴空指数,即一日内地表水平面接收到
总辐射之和与一日内地外大气上界水平面接收到太
阳辐射之和的比值.
1 6 模型评价指标
决定系数( coefficient of determination, R2)、平
均偏差(mean bias error, MBE)和均方根误差( root
mean square error, RMSE)计算公式如下:
R2 =
[∑(Ycal - Ycm)(Yobs - Yom)] 2
∑(Ycal - Ycm) 2∑(Yobs - Yom) 2
MBE =
∑(Yobs - Ycal)
n
RMSE =
∑(Yobs - Ycal) 2
n
式中:Ycal为模拟值;Ycm为模拟值的平均值;Yobs为实
测值;Yom为实测值的平均值;n为样本数量.
2 结果与分析
2 1 不同散射辐射模型的验证
分解模型的核心是构建和刻画散射分数( kd)
与晴空指数( kt)的函数关系.基于实际观测的 30
min散射辐射( ID)和总辐射( IG)计算得到 kd,同时
利用观测的总辐射( IG)和估算的地外大气上界辐射
( IE)计算得到 kt .总体而言,5 个模型在不同程度上
均能较好地模拟 kd和 kt的关系,但不同模型之间存
在明显差异.对只考虑 kt的单因子模型而言,Reindl⁃
1和 Boland 模型的模拟结果只是一条折线或曲线,
显然对 kd和 kt关系的变异解释能力不足,从而在一
定程度上限制了这类模型准确模拟散射辐射的能
力.相对于单因子模型,Reindl⁃2模型引入了太阳高
度角因子,Reindl⁃3 模型在 Reindl⁃2 模型的基础上
增加了温度和相对湿度的影响. BRL 模型则在 Bo⁃
land模型的基础上考虑了当地太阳时和日尺度晴空
指数因子.从这些多因子模型的模拟效果来看,对 kd
和 kt关系的模拟能力得到明显提高(图 2).
依据kt区分不同的天空晴朗情况,进一步分析
图 2 不同散射辐射模型的验证
Fig.2 Validation of different diffuse radiation models.
空心圆和实心圆分别表示观测和模拟的散射辐射 Solid and empty
circles denoted observed and modeled diffuse radiation, respectively.
4992 应 用 生 态 学 报 26卷
了不同模型的模拟效果.采用 Yousif 等[28]的方法进
行划分,0 ≤ kt<0.2为阴天;0.2≤kt<0.6为多云;0.6
≤kt<0.75为晴朗;0.75≤kt≤1 为非常晴朗.虽然在
总体上,不同模型的模拟结果可以反映出散射辐射
的实际变化,但是在不同的晴空条件下,模型模拟的
效果存在显著差异(表 1).不同模型均在 kt较低(阴
天到多云)时的模拟效果最好.例如,当 kt <0.2 时,
不同模型均可以解释 99.9%以上的散射辐射变化.
但是随着 kt的进一步增大,模型的模拟能力显著下
降.当 kt≥0.75 时,多个模型已经无法准确反映散射
辐射的实际变化.在不同模型之间,特别是单因子模
型,由于只考虑了 kt的影响,在 kt较大时,其模拟效
果显著低于多因子模型.相对而言,BRL 模型在高 kt
条件下(kt≥0.75)的模拟效果在一定程度上相对
较好.
2 2 散射辐射的季节变化及其模拟
由图 3 可以看出,千烟洲地区总辐射的波动很
大,并且受到了降水的明显影响,7—8 月的总辐射
最高,在此期间的日尺度晴空指数(K t)通常为全年
表 1 不同散射辐射模型性能评价指标
Table 1 Performance evaluation index of different diffuse
radiation models
模型
Model
kt n R2 平均偏差
MBE
(W·m-2)
均方根
误差
RMSE
(W·m-2)
Reindl⁃1 [0, 0.2) 1726 0.99 -0.48 0.03
[0.2, 0.6) 2683 0.84 -4.83 1.52
[0.6, 0.75) 1401 0.46 -0.43 1.50
[0.75, 1] 49 0.00 47.28 11.61
[0, 1] 5859 0.84 -1.30 0.64
Reindl⁃2 [0, 0.2) 1726 0.99 -0.89 0.04
[0.2, 0.6) 2683 0.84 -21.34 1.51
[0.6, 0.75) 1401 0.46 -21.66 1.58
[0.75, 1] 49 0.08 3.77 9.18
[0, 1] 5859 0.84 -12.90 0.65
Reindl⁃3 [0, 0.2) 1726 0.99 -0.64 0.03
[0.2, 0.6) 2683 0.87 -8.00 1.34
[0.6, 0.75) 1401 0.54 -0.46 1.36
[0.75, 1] 49 0.16 32.33 9.65
[0, 1] 5859 0.87 -2.41 0.57
Boland [0, 0.2) 1726 0.99 0.57 0.04
[0.2, 0.6) 2683 0.83 -17.79 1.55
[0.6, 0.75) 1401 0.46 -6.24 1.50
[0.75, 1] 49 0.00 31.41 10.78
[0, 1] 5859 0.83 -6.65 0.64
BRL [0, 0.2) 1726 0.99 -0.24 0.03
[0.2, 0.6) 2683 0.86 -31.76 1.50
[0.6, 0.75) 1401 0.59 -16.90 1.36
[0.75, 1] 49 0.25 14.22 8.24
[0, 1] 5859 0.87 -14.29 0.60
图 3 千烟洲地区气候要素季节变化
Fig.3 Seasonal variations of climate factors ( global radiation,
temperature, precipitation) and daily clearness index (Kt ) in
Qianyanzhou area (2012⁃03⁃01—2013⁃02⁃28).
Ⅰ: 日尺度晴空指数 Kt; Ⅱ: 总辐射 Global radiation; Ⅲ: 温度 Tem⁃
peraturp; Ⅳ: 降水 Precipitation. 虚线为年分界线 The dash line was
the boundary between the two years.
最高.气温呈单峰趋势,7月达到一年中的最高值.全
年降水为 1337.6 mm,与常年降水量相当,但季节分
配不均,特别是在总辐射最强、气温最高的 7 月,降
水明显偏少于其他月份.
由图 4 可以看出,研究区散射辐射的季节变化
总体上为单峰变化特征.但是在夏季,特别是 7 月散
射辐射明显降低,这与总辐射的季节变化趋势并不
完全一致.其中实测值选用了 SPN1 和 RSR3 的测定
结果,其测定结果一致,散射辐射都在 7 月出现降
低,这说明模型的过高估计不是因为观测仪器的问
题,可能是因为天气较为晴朗,空气中的云或水汽对
总辐射影响较小导致散射辐射减少.月尺度散射辐
射的模拟值和实测值的对比表明,5 个模型在总体
上都可以反映出大多数时段散射辐射的季节变化,
但所有模型模拟的散射辐射均表现为单峰变化,在
7 月达到最高,这导致模型模拟值在夏季显著高于
实测值.
进一步的统计表明,在年总量上,散射辐射的实
际测量值为 2258.18 MJ·m-2,该地区年尺度散射辐
射占总辐射的比例为 54.6%,说明散射辐射的量值
在全年占有很大比重.模型 Reindl⁃1、Reindl⁃2、Re⁃
indl⁃3、Boland 和 BRL 的估算结果分别为 2267. 9、
2400.4、2287.2 、2325.5 和 2418.8 MJ·m-2,模拟值
与实测值的相对偏差分别为0.04%、6.3%、1 3%、
599210期 韩佳音等: 江西千烟洲地区散射辐射的模拟与验证
图 4 各模型在月尺度上模拟值与实测值的比较
Fig.4 Comparison between modeled and observed monthly dif⁃
fuse radiation of different models.
图 5 各模型在 2012 年 7 月日尺度上模拟值与实测值的
比较
Fig.5 Comparison between modeled and observed daily diffuse
radiation of different models in July of 2012.
3 0%、7.1%,平均为 3.6%.但是在偏差明显的 7 月,
5个模型模拟的平均值与测量值的偏差达到 51 5
MJ·m-2,其中相对偏差最高为 28. 2%,最低为
14 5%,平均为 21.2%.
针对散射辐射模拟值在夏季出现偏高估算的现
象,选取偏高最为明显的 7 月进一步分析.图 5 表
明,由于同期降水较少,7月的总辐射大部分处于较
高水平,仅在 7月中旬和下旬出现了两次短暂的降
低.与全年日均 K t值 0.34 相比,7 月的 K t明显偏高,
其平均值为 0.52,高于全年平均 52.9%.从模拟值和
实测值的对比来看,偏差较大的时段大多出现在 K t
较高的情况下.选取 2012 年 7 月 2、11、15 日分别作
为晴天、多云、阴天的代表性天气.在晴朗条件下,5
个模型都显著高估散射辐射,而随着 kt的降低,多云
和阴天条件下的模拟值与实测值逐渐趋于一致(图
6).由于该地区 7月的晴天较多(K t较高),在这种条
件下各模型不容易准确模拟出散射辐射的实际变
化,从而导致该月份的模拟值显著高于实测值.
图 6 各模型在不同天气下模拟值与实测值的比较
Fig.6 Comparison between modeled and observed diffuse radia⁃
tion of different models under different weather conditions in July
of 2012.
3 讨 论
3 1 不同散射辐射模型的评价
从 30 min尺度数据的模拟结果来看,不同模型
的模拟效果一致.此外,不同模型在季节变化上的模
拟也是一致的.总体来说,5 个模型都可以模拟出该
地区的散射辐射变化.但是不同模型之间的模拟效
果存在一定差异.首先,从模型结构上来说,只考虑
kt的单因子模型比多因子模型表现较差:单因子模
型模拟的 kd与 kt的关系只能为一条曲线或折线,而
多因子模型模拟的结果为在一定区域内有规则聚集
的点.其次,在只考虑 kt的单因子模型中,分段模型
(Reindl⁃1)略优于整体模型(Boland);在多因子模
型中,整体模型(BRL)优于分段模型(Reindl⁃2 和
Reindl⁃3).虽然多因子模型的模拟效果总体上优于
单因子模型,但在实际应用中需要更多的输入变量.
在模型计算过程中,以往的研究只保留了太阳
高度角 7°以上的模拟结果,针对太阳高度角较低情
6992 应 用 生 态 学 报 26卷
况下的模型模拟能力,本研究基于 Reindl⁃2 模型对
此进行了探讨.分析表明,太阳高度角小于 7°的观测
记录共 802 条,散射辐射值总量为 19.70 MJ·m-2,
仅占全年散射辐射量(2282.91 MJ·m-2)的 0.9%.
Reindl⁃2的模拟值与观测值依然表现出显著的相关
性,R2、MBE、RMSE 分别为 0. 8289、4 97 W·m-2、
0 81 W·m-2,与 7°以上的研究结果基本相同.由于
本研究站点位于亚热带区域,较低太阳高度角出现
的频率相对高纬度地区偏少,因此对于散射辐射模
型在高纬度地区,特别是太阳高度角处于较低情况
下的模拟效果,还需要结合实际的观测数据进行评
价和分析.
考虑 R2、MBE 和 RMSE 等评价指标,5 个模型
的模拟效果排序为:BRL>Reindl⁃3>Reindl⁃2>Reindl⁃
1>Boland.在千烟洲地区 BRL的模拟效果最好,其次
为 Reindl⁃3,若是在缺少环境温度和相对湿度观测
的情况下,则 Reindl⁃2 也能给出较好的模拟结果.
Reindl⁃1和 Boland可用于快速、粗略估计散射辐射
变化趋势的情形,在精度要求较高的情况下不推荐
使用.
3 2 模拟的不确定性
从表 1 可以看出,模拟的最大不确定性来自于
高 kt情况,这一不确定性在季节尺度上主要表现为
该地区夏季散射辐射的高估.本研究选取的 5 个模
型中,在高 kt条件下 kd均表现为 kt的正相关函数,而
kt又与总辐射正相关.因此,在晴天条件较多并且辐
射强烈的夏季,模拟的散射辐射与总辐射表现出了
一致的变化(图 4).但实测数据显示,SPN1 和 RSR3
两台仪器测定的 7 月散射辐射均出现了明显降低,
这可能与随着晴空条件增多和降水偏少天气出现,
大气中云量或颗粒物的含量减少有关,从而导致散
射辐射的降低.因此散射辐射出现了与总辐射不同
的变化趋势.
当 K t<0.55时,模拟值与实测值的差值基本在 0
左右波动,说明此时模拟值与实测值较为吻合;当
K t>0.55后,模拟值与实测值的差值快速升高,说明
此时容易出现过高估计.K t是由日尺度总辐射计算
出的指标,与其有正相关关系,所以在日尺度总辐射
达到 15 MJ·m-2后,容易出现模型的过高估计.当温
度大约高于 30 ℃后,容易出现模型的过高估计.当相
对湿度小于 0.7 时,较易出现过高估计,而随着相对
湿度的增加,模拟值与实测值逐渐趋于吻合(图 7).
图 7 日尺度散射辐射模拟值和实测值差值与气候要素的关系
Fig.7 Relationship between difference of simulated values and measured values of daily diffuse radiation and climate factors.
利用聚合平均的方法,将 365个日尺度散射辐射模拟值与实测值的差值分别以各气候要素值从小到大排序,平均分为 20组,每组由 18个数据
取平均组成,余下的 5个数据取平均作为第 21组 Using the average of the aggregation method, sort 365 differences between simulated values and ob⁃
served values of daily diffuse radiation from the smallest. It was divided into 20 groups. Each group contains 18 data on average. Took the average of the
rest 5 data as 21th group.
799210期 韩佳音等: 江西千烟洲地区散射辐射的模拟与验证
千烟洲地区地处亚热带,多年平均降水量达到
1360 mm以上,其降水充沛.但是夏季容易出现干旱
少雨的情况[16],导致月尺度上 7 月的总辐射量达到
最高但降水稀少,从而使该月的 K t也较高.例如所选
用数据在 7 月的日尺度总辐射平均值达到 20 9
MJ·m-2,7 月总降水量为 76.6 mm,仅占全年降水
量的 5.7%;并且此时通常平均温度(30.0 ℃)高、相
对湿度(0.66)较低.这些条件均较易导致模型高估
的情况.这种差异在模型应用中需要注意,特别是在
夏季辐射量较高和天空云量较少的条件下需要做进
一步深入评价,有必要基于观测数据对模型进行优
化.但目前夏季高 kt条件下的观测数据还十分有限,
客观上需要进一步积累观测数据以满足模型优化和
验证的需要.
参考文献
[1] Şen Z. Atmospheric environment and renewable ener⁃
gy / / Şen Z, ed. Solar Energy Fundamentals and Model⁃
ing Techniques: Atmosphere, Environment, Climate
Change and Renewable Energy. Berlin: Springer, 2008:
21-45
[2] Szeicz G. Solar radiation for plant growth. Journal of Ap⁃
plied Ecology, 1974, 11: 617-636
[3] Harrouni S. Fractal classification of typical meteorologi⁃
cal days from global solar irradiance: Application to five
sites of different climates / / Badescu V, ed. Modeling
Solar Radiation at the Earth’ s Surface. Berlin: Sprin⁃
ger, 2008: 29-54
[4] Gu L, Baldocchi DD, Wofsy SC, et al. Response of a
deciduous forest to the Mount Pinatubo eruption: En⁃
hanced photosynthesis. Science, 2003, 299: 2035-2038
[5] Alton PB, North PR, Los SO. The impact of diffuse
sunlight on canopy light⁃use efficiency, gross photosyn⁃
thetic product and net ecosystem exchange in three forest
biomes. Global Change Biology, 2007, 13: 776-787
[6] Mercado LM, Bellouin N, Sitch S, et al. Impact of
changes in diffuse radiation on the global land carbon
sink. Nature, 2009, 458: 1014-1017
[7] Alton PB. Reduced carbon sequestration in terrestrial
ecosystems under overcast skies compared to clear skies.
Agricultural and Forest Meteorology, 2008, 148: 1641-
1653
[8] Muneer T, Younes S, Munawwar S. Discourses on solar
radiation modeling. Renewable and Sustainable Energy
Reviews, 2007, 11: 551-602
[9] Chen R, Ersi K, Yang J, et al. Validation of five global
radiation models with measured daily data in China.
Energy Conversion and Management, 2004, 45: 1759-
1769
[10] Ahmad MJ, Tiwari GN. Solar radiation models: A re⁃
view. International Journal of Energy Research, 2011,
35: 271-290
[11] Gu L, Fuentes JD, Shugart HH, et al. Responses of net
ecosystem exchanges of carbon dioxide to changes in
cloudiness: Results from two North American deciduous
forests. Journal of Geophysical Research: Atmospheres,
1999, 104: 31421-31434
[12] Ren XL, He HL, Zhang L, et al. Spatiotemporal varia⁃
bility analysis of diffuse radiation in China during 1981-
2010. Annales Geophysicae, 2013, 31: 277-289
[13] Ridley B, Boland J, Lauret P. Modelling of diffuse solar
fraction with multiple predictors. Renewable Energy,
2010, 35: 478-483
[14] Boland J, Scott L, Luther M. Modelling the diffuse frac⁃
tion of global solar radiation on a horizontal surface. En⁃
vironmetrics, 2001, 12: 103-116
[15] Boland J, Ridley B, Brown B. Models of diffuse solar
radiation. Renewable Energy, 2008, 33: 575-584
[16] Sun X, Wen X, Yu G, et al. Seasonal drought effects
on carbon sequestration of a mid⁃subtropical planted
forest of southeastern China. Science in China Series D:
Earth Sciences, 2006, 49: 110-118
[17] He X⁃Z (何学兆), Zhou T (周 涛), Jia G⁃S (贾根
锁), et al. Modeled effects of changes in the amount
and diffuse fraction of PAR on forest GPP. Journal of
Natural Resources (自然资源学报), 2011, 26(4):
619-634 (in Chinese)
[18] Li D⁃Q (李登秋), Zhou Y⁃L (周艳莲), Ju W⁃M (居
为民), et al. Modelling the effects of changes in solar
radiation on gross primary production in subtropical ever⁃
green needle⁃leaf plantations. Chinese Journal of Plant
Ecology (植物生态学报), 2014, 38(3): 219- 230
(in Chinese)
[19] The South Mountain Area Comprehensive Scientific Ex⁃
pedition, Chinese Academy of Sciences (中国科学院南
方山区综合科学考察队). Red Soil Hilly Development
and Governance: Comprehensive Experimental Study
Development Governance in Qianyanzhou. Beijing: Sci⁃
ence Press, 1989 (in Chinese)
[20] Delta⁃T Devices Ltd. SPN1 Technical Fact Sheet [EB /
OL]. (2014⁃01⁃27) [2014⁃03⁃25]. http: / / www.delta⁃
t.co.uk / product⁃support⁃material.asp
[21] Zhu Z⁃L (朱治林), Sun X⁃M (孙晓敏), Yu G⁃R (于
贵瑞), et al. Radiometers performance attenuation and
data correction in long⁃term observation of total radiation
and photosynthetically active radiation in typical forest
ecosystems in China. Chinese Journal of Applied Ecology
(应用生态学报), 2011, 22(11): 2954 - 2962 ( in
Chinese)
[22] Han J⁃Y (韩佳音), Li S⁃G (李胜功), Zhang L⁃M
(张雷明), et al. Comparative observation of diffuse ra⁃
diation in Qianyanzhou during the spring of 2012. Chi⁃
nese Journal of Applied Ecology (应用生态学报),
8992 应 用 生 态 学 报 26卷
2015, 26(3): 697-703 (in Chinese)
[23] Foken T, Wichura B. Tools for quality assessment of
surface⁃based flux measurements. Agricultural and Forest
Meteorology, 1996, 78: 83-105
[24] Younes S, Claywell R, Muneer T. Quality control of so⁃
lar radiation data: Present status and proposed new ap⁃
proaches. Energy, 2005, 30: 1533-1549
[25] Reindl DT, Beckman WA, Duffie JA. Diffuse fraction
correlations. Solar Energy, 1990, 45: 1-7
[26] Tovar⁃Pescador J. Modelling the statistical properties of
solar radiation and proposal of a technique based on bolt⁃
zmann statistics / / Badescu V, ed. Modeling Solar Radi⁃
ation at the Earth’ s Surface. Berlin: Springer, 2008:
55-91
[27] Şen Z. Solar radiation deterministic models / / Şen Z, ed.
Solar Energy Fundamentals and Modeling Techniques:
Atmosphere, Environment, Climate Change and Renew⁃
able Energy. Berlin: Sringer, 2008: 72-74
[28] Yousif C, Quecedo GO, Santos JB. Comparison of solar
radiation in Marsaxlokk, Malta and Valladolid, Spain.
Renewable Energy, 2013, 49: 203-206
作者简介 韩佳音,男,1989年生,博士研究生. 主要从事生
态系统生态学研究. E⁃mail: hanjy.12s@ igsnrr.ac.cn
责任编辑 杨 弘
封 面 说 明
图片由上海海洋大学与黑龙江水产研究所联合培养研究生王臣 2015 年 6 月 7 日拍摄于松花
江流域汤原县竹帘镇江畔.汤原县位于黑龙江省东北部,三江平原西部边缘,北依小兴安岭千里群
山,南邻松花江、汤旺河两大水系.汤原县属中温带大陆季风气候,四季分明,天气多变,夏季暖热多
雨,冬季寒冷干燥,雨热同季,冬夏气温变幅很大.全年无霜期 136.6 d,年有效积温平均值为 2559.5
℃,年平均降水量 536.8 mm,有三分之二集中于夏秋两季.汤旺河与松花江在此交汇,鱼类资源极其
丰富,经实际调查,采集到鱼类分别隶属于 8目 12科 46种,该地也曾是大麻哈鱼的重要产卵地,近
年来逐渐遭到破坏.为恢复大麻哈鱼资源,黑龙江水产研究所在此建立大麻哈鱼鱼类保护和繁育工
作站,2013—2015年连续 3年在此开展大麻哈鱼增殖放流活动,对大麻哈鱼鱼类资源的恢复具有
重大意义.
999210期 韩佳音等: 江西千烟洲地区散射辐射的模拟与验证