首先以考虑地形遮蔽的分布式可照时间理论模型结合气象台站观测的日照时间计算了日照百分率,然后对我国第二代静止气象卫星(FY-2G)的总云量遥感影像产品重采样后,根据日照百分率和云量的负相关性,分别建立了基于遥感总云量和观测站云量(总云量、低云量)的日尺度日照百分率估算模型,并以气象站点分布稀疏的新疆区域为例,对两种估算模型使用IDW和Kriging插值方法的模拟效果进行了检验,得出以下结论:1) 实际地形下新疆区域可照时间的空间分布受地形影响较大,全区四季平均可照时间分别为:春季1165 h,夏季1286 h,秋季964 h,冬季823 h。2) 基于遥感云量的条带状重采样方案考虑了日照轨迹和云的区域移动变化,重采样后的云量值与日照百分率的相关性有所提升,相关系数为0.756。3) 建立的单站分季节遥感集成日照百分率模型相关系数冬夏差异明显,夏季最高,春秋次之,且分布形态较为一致,冬季最低,低值主要集中在北疆沿天山一带。4) 从模拟的效果来看,遥感集成日照百分率模型(平均绝对误差为14.8%)要明显优于基于观测站云量模型模拟的结果,由于遥感集成日照百分率充分发挥了卫星在空间上连续观测的优势,空间分布更为连续,通过检验后的该方法可以在站点稀少的西部地区进行业务应用。 The sunshine percentage was calculated using the theoretical model of the distributed exposure time considering the terrain factors and the sunshine hours observed by the meteorological station. After re-sampling the total cloud amount remote sensing image products by the second generation geostationary meteorological satellite (FY-2G) in China, the daily scale sunshine percentage estimation model was established based on remote sensing total cloud amount and observation station cloud amount (total cloud amount and low cloud amount) used the negative correlation between cloud cover and sunshine percentage. The spatial simulations are completed of the two interpolation methods IDW and Kriging in Xinjiang. Finally, simulation results are verified by the test station. The conclusions are as follows: 1) The spatial distribution of the four seasons in Xinjiang has obvious changes. The average seasons in the whole region are 1165 h in spring, 1286 h in summer, 964 h in autumn and 823 h in winter. 2) The resampling scheme based on remote sensing cloud cover takes into account the cloud movement. The correlation was improved after resampling between cloud cover and sunshine percentage, and the correlation coefficient was 0.756. 3) The correlation coefficient of sunshine percentage model based on Remote Sensing has obvious differences in the two seasons in winter and summer; the highest correlation was found in summer; the distribution pat-terns were similar in spring and autumn; the correlation coefficient was the lowest in winter and represented along the Tianshan Mountains in the northern Xinjiang. 4) The sunshine percentage accuracy of the RS integrated model is obviously better than that based on the observation station cloud cover model. The mean absolute error of the sunshine percentage was 14.8%. As that gives full play to the superiority of the continuous observation of satellites in space, the spatial distribution is more continuous. The post-verification method can be used for business applications in the sparsely stationed western region.
陈鹏翔1,2,彭冬梅3,张旭1
1新疆维吾尔自治区气候中心,新疆 乌鲁木齐
2新疆维吾尔自治区政府投资项目评审中心,新疆 乌鲁木齐
3新疆维吾尔自治区兴农网信息中心,新疆 乌鲁木齐
收稿日期:2018年8月23日;录用日期:2018年8月30日;发布日期:2018年9月6日
首先以考虑地形遮蔽的分布式可照时间理论模型结合气象台站观测的日照时间计算了日照百分率,然后对我国第二代静止气象卫星(FY-2G)的总云量遥感影像产品重采样后,根据日照百分率和云量的负相关性,分别建立了基于遥感总云量和观测站云量(总云量、低云量)的日尺度日照百分率估算模型,并以气象站点分布稀疏的新疆区域为例,对两种估算模型使用IDW和Kriging插值方法的模拟效果进行了检验,得出以下结论:1) 实际地形下新疆区域可照时间的空间分布受地形影响较大,全区四季平均可照时间分别为:春季1165 h,夏季1286 h,秋季964 h,冬季823 h。2) 基于遥感云量的条带状重采样方案考虑了日照轨迹和云的区域移动变化,重采样后的云量值与日照百分率的相关性有所提升,相关系数为0.756。3) 建立的单站分季节遥感集成日照百分率模型相关系数冬夏差异明显,夏季最高,春秋次之,且分布形态较为一致,冬季最低,低值主要集中在北疆沿天山一带。4) 从模拟的效果来看,遥感集成日照百分率模型(平均绝对误差为14.8%)要明显优于基于观测站云量模型模拟的结果,由于遥感集成日照百分率充分发挥了卫星在空间上连续观测的优势,空间分布更为连续,通过检验后的该方法可以在站点稀少的西部地区进行业务应用。
关键词 :新疆,遥感,日照时间,重采样
Copyright © 2018 by authors and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
日照百分率是表征气候变化的主要气象要素之一 [
近些年随着空间技术的发展,卫星所提供的资料在空间和时间上弥补了地面观测的局限性,越来越多的学者开始将卫星遥感数据应用于地面太阳辐射有关参量的反演 [
我国也有自主研发的风云卫星云产品,对本国卫星数据质量进行检验和评估,可以为数据进一步改进提供科学依据。风云二号气象卫星(FY-2)是我国自行研制的第一代地球静止轨道气象卫星,与极地轨道气象卫星相辅相成,构成我国气象卫星应用体系。王怀清等 [
新疆地处亚欧大陆腹地,位于我国西北部,面积约占我国陆地面积的六分之一,其四面被高山所环绕,中部有天山横亘,从而形成了“三山夹两盆”的特殊地貌格局。由于新疆远离海洋,使得湿润的海洋气流难以进入,形成了极端干燥的大陆性气候,是世界典型干旱区之一,同时,新疆具有得天独厚的水土光热资源,日照时间长、积温多、昼夜温差大、无霜期长,年太阳能辐射量在全国范围内仅次于西藏,具有丰富的太阳能光热资源 [
本文所用的地面观测资料包括新疆维吾尔自治区境内的102个气象站2015年6月1日-2016年5月31日(以下称为全年)逐日日照时数的观测数据。其中65个气象站具有总云量、低云量的观测数据,因此用于建立日照百分率估算模型,剩余37个气象站作为检验站点。图1出示了新疆区域建模站点和检验站点的间分布。
FY-2G静止卫星主要搭载扫描辐射计和空间环境监测器两个有效载荷,其可见光红外自旋扫描辐射仪(VISSR)具有5个通道,可以获取白天可见光云图、昼夜红外云图和水气分布图 [
基于FY-2G总云量的日照百分率计算主要考虑三个方面:复杂地形下的可照时间的计算、基于FY-2G总云量重采样方案的确定和日照百分率估算模型的建立。首先利用分布式天文可照时数理论模型,结合地形遮蔽的分析计算复杂地形下可照时间,气象站观测的日照时间与对应的地理可照时间的之比值就是该站点的实际日照百分率。然后对FY-2G总云量遥感数据进行重采样分析,确定最优重采样方案,最后利用云量与日照百分率之间良好的负相关关系,使用FY-2G总云量遥感数据建立估算模型得到日照百分率栅格数据。图2给出了计算流程示意图。
日照百分率是台站观测到的日照时间与相应时段可照时间的比值,“可照时间”定义有两种,即天文可照时间和地理可照时间,前者是指不考虑大气影响和地形遮蔽的最大可能日照时间,后者是指考虑地形遮蔽而不考虑大气影响的可能日照时间,某地的天文可照时数可根据有关太阳轨迹方程精确计算得到 [
图1. 新疆气象站点分布图
图2. 基于FY-2G总云量的日照百分率计算流程图
图3. 地理可照时间模型示意图
图4为新疆区域春季(3~5月)、夏季(6~8月)、秋季(9~11月)及冬季(12月、1~2月)可照时间的空间分布。实际地形下新疆各季可照时间的空间分布差异明显,具有明显的季节变化特征。全区四季可照时间平均分别为:春季1165 h,夏季1286 h,秋季964 h,冬季823 h,可照时间以夏季最长,春秋季次之,冬季最短,区域内春秋两季变化差异较夏冬两季平缓,表现出四季分布的不对称性。从夏季和冬季的可照时间空间分布中可以看出有较为显著纬向分布特征,夏季可照时间表现出自北向南的递减,而冬季则表现出相反分布特点。无论哪一个季节,在三大山脉中由地形对可照时间造成的影响都非常明显。
气象站观测到的总云量是指云遮蔽天空视野的成数,其物理意义是“视云量”,也就是云在观测者“视角”中占的比例,与卫星资料中的云量具有不同的物理意义。实际情况下,某一格点周围一定区域范围内的云都有可能对该格点的日照情况产生影响,因此,为了寻求影响日照的最佳云覆盖范围,通常使用重采样分析方法。常用的遥感重采样方法为窗口分析,以对应格点为中心,按东西和南北两个方向上等距的正方形区域为窗口进行采样。Shi等 [
对于重采样后的云量数据,以全年FY-2G总云量与日照百分率相关系数的绝对值为指标来选取最优重采样方案。图6为全年3种网格采样方案下FY-2G总云量与日照百分率的线性相关变化,h代表采样方案中南北方向的网格递增数量,w代表采样方案中东西方向上的网格递增数量,从图中可以看出,随着采样方案中南北向网格的递增整体相关系数呈减小趋势,其中南北向递增1个网格的相关系数整体最大。3种采样方案的相关系数随着东西向网格的递增呈现出先增大后减小的趋势,均在网格递增量为3个时达到顶点,当采样方案中南北向网格递增量取1、东西向网格递增量取3时全年总云量与日照百分率的相关系数最大,为0.756,因此,选取气象站所在位置南北方向各递增1个像元,东西方向各递增3个像元共21个像元为采样区域,其算术平均值作为该点的总云量值。
根据日照百分率与云量之间较好的负相关关系,建立两种日照百分率估算模型,模型I为一元线性
图4. 新疆区域四季可照时间空间分布(a:春季;b:夏季;c:秋季;d:冬季)
图5. 网格点重采样示意图
图6. 不同采样方案的相关系数变化
回归模型,以使用重采样方案后的FY-2G日总云量数据作为日照百分率的影响因子,模型II为二元线性回归模型,以观测站的总云量和低云量同时作为影响日照百分率的因子。
有基于MODIS云量的研究指出 [
按上述确定的最优重采样方案,对全年每日的FY-2G总云量影像数据进行重采样,利用建立的基于遥感云量的日照百分率模型,完成全年每日的日照百分率的估算,同时使用模型II完成基于观测站云量的日照百分率估算,在计算过程中模型系数和常数项的插值方法分别采用IDW和Kriging,从空间分布和模拟误差两个方面来分析不同模型不同插值方案对日照百分率的模拟程度。
图8是使用不同插值方案下的遥感模型与观测站模型日照百分率空间分布对比图,(a)~(d)分别为2015年7月15日,2015年10月15日,2016年1月15日,2016年4月15日遥感集成的日照百分率空间分布,(e)~(h)为对应日期观测站云量模型计算的日照百分率空间分布,1为使用IDW插值方案,2是使用Kriging插值方案,从图中可以看出,不论使用何种插值方案使用遥感集成与观测站云量模拟的日照百分率在空间分布上吻合度都较高,但是遥感集成日照百分率的空间分布更加连续,结果更加客观。
图7. 不同季节模型I的相关系数分布(a:春季,b:夏季,c:秋季,d:冬季)
图8. 遥感集成与观测站插值日照百分率空间分布对比
使用未参与模型运算的37个验证站点全年各月15日的日照百分率观测数据,与两种模型计算的日照百分率数据进行检验。为了直观和全面地评估不同方法对日照百分率的模拟能力,引入泰勒图分析方法,它是由不同模型模拟的结果与实际观测的相关系数及其均方差比值(均方差比值是模拟和观测之间离散程度的比值,1代表与实际观测一致)组成的极坐标图。图9显示了使用不同插值方法的遥感云量集成模型和观测站云量模型模拟的日照百分率与实际观测场的泰勒图分布,圆点代表不同方法模拟的结果,距离原点(黑色实心方块)的越近,表示模拟的效果越好(图中不同方法的表示:ST代表基于观测站云量模型、RS代表遥感云量模型、IDW和Kriging代表不同插值方法)。从图中可以看出,使用观测站云量模型模拟的结果相关系数都在0.6以下,使用遥感云量模型模拟结果相关系数都超过了0.7,其中使用Kriging插值方案的相关系数最高,达到了0.72。从均方差之比(黑色虚线)来看,使用遥感云量集成模型要比观测站云量模型更接近1,这说明使用遥感集成模型模拟的空间变率更加接近实际观测。对于标准化均方根误差来说,使用观测站云量模型模拟的日照百分率与实际观测的均方根误差都在0.8以上,而使用遥感云量模型的均方根误差减少至0.7左右,说明基于遥感集成模型模拟的日照百分率要比基于观测站云量模型模拟的偏差幅度明显降低。对于不同的插值方法,使用IDW和Kriging插值方案的遥感集成日照百分率模型的平均绝对误差为15.1%和14.8%,整体来说小于基于观测云量模拟的日照百分率(平均绝对误差分别为19.09%和19.77%),其中Kriging插值方案的结果要略优于IDW插值方案,最后将本文计算结果
图9. 四种方法模拟的日照百分率与实际观测场的泰勒图分布。(弧度轴表示模拟场和观测场的相关系数,黑色虚线表示均方差之比,圆点到observed的距离(绿色实线)为模拟场与观测场的标准化均方根误差;字母序号代表不同模拟方法,见右上角)
与国内相关研究做了比对,Shi等 [
本文首先计算了考虑地形遮蔽的气象站日照百分率,同时引入FY-2G的总云量遥感影像产品,对遥感总云量数据在一定范围内重采样后,根据云量和日照百分率的负相关性,分别建立了基于FY-2G总云量和观测站总云量、低云量的日尺度日照百分率模型,并使用IDW和Kriging两种插值方法完成新疆区域全年各月15日的日照百分率空间模拟,得出以下结论:
1) 实际地形下新疆四季可照时间的空间分布差异明显,具有明显的季节变化特征。全区四季可照时间平均分别为:春季1165 h,夏季1286 h,秋季964 h,冬季823 h。冬、夏两季的可照时间空间分布有较为显著纬向分布特征。三大山脉中由地形对可照时间造成的影响全年都非常明显。
2) 基于遥感云量的条带状重采样方案考虑了日照轨迹和云的区域移动变化,通过实验得出的最优重采样方案为选取气象站所在位置南北方向各递增1个像元,东西方向各递增3个像元,共21个像元的算术平均值作为该点的总云量,与日照百分率的相关系数为0.756。
3) 本文建立的单站分季节的遥感集成日照百分率模型相关系数在夏季最高,单站相关系数均在0.8以上;春秋次之,北疆相关系数整体要高于南疆,且分布形态较为一致,相关较低的地方主要在南疆中部及阿克苏地区周边;冬季最低,南疆的相关系数普遍高于北疆,相关系数分布最低的地方主要集中在北疆沿天山一带。
4) 使用检验站点对估算结果的检验分析显示,基于遥感集成模型的模拟结果与实际观测值的相关系数、均方差比值以及标准化均方根误差三个指标都要明显优于基于观测站云量的模型模拟的结果,使用遥感集成日照百分率模型的平均绝对误差为14.8%,小于基于观测云量模拟的日照百分率(平均绝对误差分别为19.09%),其中Kriging插值方案的结果要略优于IDW插值方案。同时,使用两种模型两种插值方案模拟的日照百分率空间分布都较为一致,但是,基于遥感集成的日照百分率模型充分发挥了卫星在空间上连续观测的优势,克服了气象站点分布不规律的限制,模拟结果更为客观,该方法可以在站点稀少的西部地区进行业务应用。
国家自然科学基金(41505101)。
陈鹏翔,彭冬梅,张 旭. 基于FY-2G总云量的日照百分率估算及检验 Estimated and Tested of Sunshine Percentage Based on Total Cloud Amount Form FY-2G[J]. 气候变化研究快报, 2018, 07(05): 381-390. https://doi.org/10.12677/CCRL.2018.75041