doi: /qxxb 气象学报青藏高原犌犔犃犛犛地表反照率产品精度分析陈爱军 1,2 胡慎慎 1,2 卞林根 3 刘玉洁 4 CHEN Aijun 1,2 HUShenshen 1,2 BIANLingen 3 LIU Yujie 4 1. 中国气象局气溶胶与云降水

Size: px

Start display at page:

Download "doi: /qxxb 气象学报青藏高原犌犔犃犛犛地表反照率产品精度分析陈爱军 1,2 胡慎慎 1,2 卞林根 3 刘玉洁 4 CHEN Aijun 1,2 HUShenshen 1,2 BIANLingen 3 LIU Yujie 4 1. 中国气象局气溶胶与云降水"

层苏
5 years ago
Views:

1 doi: /qxxb 气象学报青藏高原犌犔犃犛犛地表反照率产品精度分析陈爱军 1,2 胡慎慎 1,2 卞林根 3 刘玉洁 4 CHEN Aijun 1,2 HUShenshen 1,2 BIANLingen 3 LIU Yujie 4 1. 中国气象局气溶胶与云降水重点开放实验室, 南京, 南京信息工程大学大气物理学院, 南京, 中国气象科学研究院, 北京, 国家卫星气象中心, 北京, 犓犲狔犔犪犫狅狉犪狋狅狉狔犳狅狉犃犲狉狅狊狅犾犆犾狅狌犱犘狉犲犮犻狆犻狋犪狋犻狅狀狅犳犆犺犻狀犪犕犲狋犲狅狉狅犾狅犵犻犮犪犾犃犱犿犻狀犻狊狋狉犪狋犻狅狀, 犖犪狀犼犻狀犵 , 犆犺犻狀犪 2. 犛犮犺狅狅犾狅犳犃狋犿狅狊狆犺犲狉犻犮犘犺狔狊犻犮狊, 犖犪狀犼犻狀犵犝狀犻狏犲狉狊犻狋狔狅犳犐狀犳狅狉犿犪狋犻狅狀犛犮犻犲狀犮犲牔犜犲犮犺狀狅犾狅犵狔, 犖犪狀犼犻狀犵 , 犆犺犻狀犪 3. 犆犺犻狀犲狊犲犃犮犪犱犲犿狔狅犳犕犲狋犲狅狉狅犾狅犵犻犮犪犾犛犮犻犲狀犮犲狊, 犅犲犻犼犻狀犵 , 犆犺犻狀犪 4. 犖犪狋犻狅狀犪犾犛犪狋犲犾犻狋犲犕犲狋犲狅狉狅犾狅犵狔犆犲狀狋犲狉, 犅犲犻犼犻狀犵 , 犆犺犻狀犪收稿, 改回. 陈爱军, 胡慎慎, 卞林根, 刘玉洁青藏高原 GLASS 地表反照率产品精度分析. 气象学报,73(6): 犆犺犲狀犃犻犼狌狀, 犎狌犛犺犲狀狊犺犲狀, 犅犻犪狀犔犻狀犵犲狀, 犔犻狌犢狌犼犻犲犃狀犪狊犲狊犿犲狀狋狅狀狋犺犲犪犮犮狌狉犪犮狔狅犳狋犺犲犌犔犃犛犛犪犾犫犲犱狅狆狉狅犱狌犮狋狊狅狏犲狉狋犺犲犜犻犫犲狋犪狀犘犾犪狋犲犪狌. 犃犮狋犪犕犲狋犲狅狉狅犾狅犵犻犮犪犛犻狀犻犮犪,73(6): 犃犫狊狋狉犪犮狋 TheaccuraciesoftheGLASS(GlobalLAndSurfaceSatelites)albedoproductsin1kmresolutionareassessedusing thegroundmeasurementsofyear2003atthreestationsoverthetibetanplateau.theresultsshowthattheoveralchange trendsoftheglassblack skyandwhite skyalbedosareconsistentwiththeground observedalbedos,whichcouldindicatethe stateofsurfacefairwel;boththesmalfractionofsnowcoverandcloudcoverhavegreatefectontheaccuracyoftheglass albedoproductsandcloudcovermightresultinhighervaluesofglassalbedoproductsthangroundmeasurements;therm SEs(root mean squareerrors)betweentheglassblack skyalbedos,glass white skyalbedosandground measurements areevidentlydecreasedwheneliminatingtheefectofsnowandcloudcover,whichare0.0155and0.019,respectively. 犓犲狔狑狅狉犱狊 GLASS,Albedo,theTibetanPlateau,Accuracyassessment 摘要应用 2003 年青藏高原 3 个站点的地表反照率观测结果, 对比分析了 GLASS(GlobalLAndSurfaceSatelites) 地表反照率 1km 1km 分辨率产品的精度, 结果表明,GLASS 黑空反照率白空反照率与地表反照率地面观测结果的总体变化趋势基本一致, 能够有效地反映实际地表状态的变化 ; 局地积雪和云覆盖对 GLASS 地表反照率产品的精度影响较大, 云覆盖导致 GLASS 地表反照率可能比实际地表反照率高 ; 消除云覆盖和局地积雪的影响后,GLASS 黑空反照率白空反照率与地表反照率地面观测结果的均方根误差显著降低, 分别为和关键词 GLASS, 地表反照率, 青藏高原, 精度分析中图法分类号 P405 1 引言地表反照率定义为地表反射太阳辐射与入射太阳辐射之比, 不仅影响地气系统的辐射能量收支与平衡, 而且对全球气候变化也有重要影响, 因而地表反照率是陆面过程模式及气候模式中的一个重要参数 (Dickinson,1983; 王鸽等,2010) 青藏高原是世界上海拔最高面积最大的高原, 其大范围热动力作资助课题 : 国家自然科学基金 ( ) 江苏省政府留学基金作者简介 : 陈爱军, 主要从事气象卫星资料分析与应用研究 E mail:chenaijun@nuist.edu.cn

2 陈爱军等 : 青藏高原 GLASS 地表反照率产品精度分析 1115 用及地气相互作用过程, 对中国东部和南部以及亚洲地区的灾害性天气, 甚至全球气候变化均有重大影响 ( 叶笃正等,1979) 研究并获得青藏高原地表反照率, 无疑有助于更好地认识和理解青藏高原地气相互作用过程及其热动力作用和气候影响青藏高原地形复杂, 气候条件恶劣, 只有通过卫星遥感反演才能获得时空分布连续的高时空分辨率地表反照率 ( 陈爱军等,2009) 目前, 除美国 Terra 和 Aqua 双星 MODIS (MODerate resolutionima gingspectroradiometer) 全球地表反照率产品 MCD43(Schaaf,etal,2002,2011; 陈爱军等,2009) 欧洲全球地表反照率产品 GlobAlbedo(Lewis,et al,2011) 之外, 北京师范大学全球陆表特征参量产品生成与应用研究重点项目开发的全球变化与地球系统科学研究的参数产品 GLASS(GlobalLAnd SafaceSatelites)( 梁顺林等,2014), 也可以提供青藏高原卫星遥感反演的地表反照率 GLASS 地表反照率产品不仅是中国唯一自主开发的高精度时空连续的全球地表反照率产品, 也是目前国际上时间序列最长的全球地表反照率产品 ( 梁顺林等, 2014), 对于促进青藏高原地气相互作用研究具有重要的应用价值气候和陆面过程模式对地表反照率的精度要求较高, 一般要求其绝对精度在 (Selers, etal,1995) 近年来, 不少研究分析和评价了 GLASS 地表反照率产品的精度, 例如 : 王立钊等 (2014) 利用北美洲和亚洲地区的 5 个全球通量观测网络 (FLUXNET) 站点数据和 LandsatTM 数据验证了 GLASS 反照率产品的精度, 结果表明,GLASS 地表反照率产品的精度较高, 总体误差约为 ; 齐文栋等 (2014) 利用年全球 24 个站点数据和 MODIS 数据对比分析了 GLASS 地表反照率产品的精度, 结果表明 GLASS 地表反照率产品与地面观测数据的均方根误差约为 0.05, 与 MCD43B3 产品的均方根误差约为 0.03; 王健等 (2013) 利用中国和美国地区农田常绿针叶林和混交林等不同地表覆盖类型的 4 个 FLUXNET 站点数据分析了 GLASS 地表反照率产品的精度, 结果表明地表覆盖均一区域的精度较高, 而地表覆盖复杂地区的精度较低然而, 分析和验证青藏高原 GLASS 地表反照率产品精度的研究仍然十分少见青藏高原地形复杂云覆盖度较高, 而且天气变化剧烈, 地表状态变化较为频繁这些因素一方面给青藏高原地表反照率的卫星遥感反演带来挑战, 另一方面也促使用户在分析和应用卫星遥感反演结果时, 必须高度重视其反演质量和精度, 尤其用于气候和陆面过程模式时因此, 将利用 CAMP/Tibet 国际合作科学观测试验期间自动气象站观测的地表辐射数据, 分析和验证 GLASS 地表反照率产品青藏高原地区的精度研究一方面有助于增进人们进一步了解青藏高原 GLASS 地表反照率产品的精度, 另一方面也有助于促进其在青藏高原地气相互作用研究中发挥应有的作用 2 GLASS 地表反照率产品简介目前,GLASS 地表反照率产品的时间范围为年, 其中, 年的 GLASS 地表反照率产品由 AVHRR(Advanced Very High ResolutionRadiometer) 数据生成, 空间分辨率为 , 采用等经纬度投影 ; 年的 GLASS 地表反照率产品由 MODIS 数据生成, 空间分辨率为 1km 1km, 采用正弦投影 (Liu,etal, 2013b) GLASS 地表反照率产品首先采用 AB(Angu larbin) 算法生成中间产品, 然后采用 STF (Statis tics basedtemporalfiltering) 算法合成终级产品 (LiuN F,etal,2011,2013; 刘强等,2012; 梁顺林等,2014) AB 算法以多波段地表反射率与短波波段地表反照率的线性统计关系为基础, 分为两种 : 一种采用大气校正后的多波段地表反射率, 记为 AB1; 另一种采用大气层顶的多波段反射率, 记为 AB2 STF 算法的实质是一种数据融合方法, 对 GLASS 地表反照率中间产品进行时间序列平滑缺失填补和数据融合等处理 (LiuQ,etal,2013; 刘强等,2012; 梁顺林等,2014) 基于 MODIS 数据的 GLASS 地表反照率中间产品有 4 种 : 采用 AB1 算法分别由大气校正后的 Terra 和 Aqua MODIS 数据生成的记为 GLASS02A21 和 GLASS02A22; 采用 AB2 算法分别由 Terra 和 Aqua MODIS 大气顶的反射率数据生成的记为 GLASS02A23 和 GLASS02A24 然后, 采用 STF 算法将这 4 种逐日生成的中间产品每 8d 以 17 d 的时域滤波窗口合成终级产品 GLASS02A06 基于 AVHRR 数据的 GLASS 地表反照率产品由 AVHRR 数据采用 AB1 算法生成, 并采用 STF 算法每 8d 以 33d 的时域滤波窗口合

3 1116 犃犮狋犪犕犲狋犲狅狉狅犾狅犵犻犮犪犛犻狀犻犮犪气象学报 2015,73(6) 成终级产品 GLASS02B05 文中所用数据为 GLASS 地表反照率产品 GLASS02A06( 第 1 版 ) 该产品主要提供反照率及数据质量信息, 其中, 反照率包括短波波段 ( μm) 的黑空反照率 BSA(Black SkyAlbedo) 和白空反照率 WSA(White SkyAlbedo); 数据质量信息包括反照率的整体质量有效性不确定性地表覆盖状态 STF 算法时域滤波窗口的时间, 以及输入数据信息 3 地面观测资料及其处理研究所用地面辐射观测数据来自 CAMP/Tibet 国际合作科学观测试验期间 BJ D105 和 NPAM 站自动气象站 2003 年观测的下行和上行太阳总短波辐射资料这 3 个站点所处位置比较开阔, 能够代表藏北高原大范围的地形和典型的下垫面特征 ( 表 1), 其中,BJ 站地面覆盖高约 5cm 稀疏的高原草甸, 土壤由含有碎石的砂壤土构成 ;D105 站位于多年冻土区, 地表为淤泥土壤和砂质土壤 ;NPAM 站地表为含有细石块的砂黏土, 夏季不均匀地生长高约 15cm 的草甸观测试验期间, 所有仪器都经过严格的标定和校准, 确保观测数据的准确性表 1 地面观测站点的基本信息 Table1 Thebasicinformationoftheground basedstations 站名经度 ( E) 纬度 ( N) 海拔高度 (m) 地表覆盖类型 BJ 高原草甸 D 多年冻土区 NPAM 草甸为了获得与研究所用 GLASS 地表反照率产品对应的地面观测结果, 对于这 3 个自动气象站观测的下行和上行太阳短波总辐射资料, 采用如下方法做处理 : 首先, 剔除明显有误的观测数据, 原则是下行太阳短波总辐射不低于上行太阳短波总辐射其次, 计算小时平均地表反照率将 10min 或 30min 采样的下行和上行太阳短波总辐射观测数据按小时平均, 筛选出晴空观测数据, 再计算小时平均的地表反照率采取的晴空标准为下行太阳短波总辐射不低于 600 W/m 2 第 3 步, 计算日平均地表反照率由于 GLASS 地表反照率是当地正午时刻的地表反照率, 为了尽可能地获得与 GLASS 地表反照率对应且具有代表性的日平均地表反照率地面观测结果, 取当地正午时刻前后 1h 的小时平均地表反照率的平均值最后, 计算 8d 平均地表反照率 GLASS 地表反照率终级产品是每 8d 以 17d 的时域滤波窗口合成的平均地表反照率, 例如 : 时间标记为第 9 天的 GLASS 地表反照率终级产品, 实际是第 1 17 天日地表反照率的合成结果 (Liang,etal,2014) 计算 8d 平均地表反照率的地面观测结果时, 亦采用与 GLASS 地表反照率完全相同的时间周期考虑到积雪对地表反照率的影响, 按 17d 中无雪日或积雪日的多数取平均值, 其中积雪日标准为当地正午时刻的地表反照率不小于 0.4 图 1 给出的是 2003 年 BJ 站观测的日平均地表图年 BJ 站反照率地面观测处理结果 Fig.1 GroundmeasurementsofalbedoattheBJstationin2003

4 陈爱军等 : 青藏高原 GLASS 地表反照率产品精度分析 1117 反照率 ( 记为 αd) 和 8d 平均地表反照率 ( 记为 ), 其中, 第 1 天开始的 8d 的平均地表反照率由 2002 年第天 2003 年第 1 9 天的日平均地表反照率 αd 计算第 1 25 天的日平均地表反照率均高于 0.4, 为积雪日, 因而第 9 17 天开始的 8d 的平均地表反照率均高于 0.4, 分别为 0.71 和 0.66 第天, 虽然存在无雪日, 但无雪日明显少于积雪日, 故以积雪日日平均地表反照率的平均值 0.59 作为第 25 天开始的 8d 的平均地表反照率相反, 第和 137 天开始的 8d, 虽然也存在积雪日, 但无雪日明显多于积雪日, 应以无雪日日平均地表反照率的平均值作为对应的 8d 平均地表反照率, 因而它们的均小于 0.3 第天开始的 8d, 日平均地表反照率均小于 0.4, 全部为无雪日, 它们的 8d 平均地表反照率 () 与对应的 αd 没有明显的差异, 具有较强的一致性 4 结果分析黑空和白空反照率分别代表太阳辐射完全直射和完全漫射条件下的地表反照率, 实际地表反照率需要由两者根据天空漫射光比例内插获得 (Schaaf, etal,2002,2011; 陈爱军等,2009) 天空漫射光比例与气溶胶光学厚度密切相关, 获得这 3 个地面观测站点的气溶胶光学厚度数据又存在一定的困难, 因此, 将以地面观测站点对应像元的 GLASS 短波黑空和白空反照率直接与地面观测结果进行比较图 2 给出的是 2003 年 BJ 站观测的 8d 平均地表反照率与对应像元的 GLASS 黑空反照率 ( 记为 αbsa) 白空反照率 ( 记为 αwsa) 的比较结果从图中可以看出 : 除第天开始的 4 个 8d 之外,αbsa 和 αwsa 与其余 8d 的平均地表反照率具有较强的一致性第 1 天开始的 8d,αbsa 分别比偏低 0.16 和 0.17; 第 9 天开始的 8d,αbsa αbsa 都比偏低 0.13 不过, 这两个 8d 的 αbsa 和远超过 0.4, 都符合积雪条件, 表明地表确实存在积雪这说明 GLASS 地表反照率能够有效地反映地表积雪状态, 同时说明地面积雪可能导致 GLASS 地表反照率比地面观测结果低第天开始的 8d 平均地表反照率都远超过 0.4, 表明 BJ 站存在积雪然而,GLASS 黑空反照率 (αbsa) 白空反照率 (αwsa) 都小于 0.4, 与的绝对偏差超过 0.3 分析 BJ 站对应像元的 Terra 和 Aqua MODIS 地表反射率, 在以第 17 天为中心的 17d 时域滤波窗口, 大气校正后地表反射率小于 0.4 的天数都不少于 10d; 在第 25 天为中心的 17d 时域滤波窗口,Terra 和 AquaMODIS 地表反射率小于 0.4 的天数均为 13d 也就是说, 第天开始的 8d 的 GLASS 黑空反照率和白空反照率不超过 0.4 完全合理这表明在第天开始的 8d,αbsa 与存在较大偏差的原因, 极有可能是观测站点附近局地积雪的影响 : 虽然局地积雪导致 BJ 站的 8d 平均地表反照率偏高, 但是它很难显著影响站点对应像元的地表辐射, 因而根据卫星观测辐射反演的 GLASS 黑空和白空反照率较低图年 BJ 站 GLASS 地表反照率与地面观测结果的比较 Fig.2 ComparisonbetweentheGLASSalbedosandground observedalbedosatthebjstationin2003

5 1118 犃犮狋犪犕犲狋犲狅狉狅犾狅犵犻犮犪犛犻狀犻犮犪气象学报 2015,73(6) 图 3 给出的是 2003 年 D105 站观测的 8d 平均地表反照率 () 与站点对应像元 GLASS 地表反照率的对比结果从图中可以看出 :GLASS 黑空反照率 (αbsa) 白空反照率 (αwsa) 与地面观测的 8d 平均地表反照率 () 都在 0.4 以下, 且整体变化趋势基本一致, 年平均偏差均为 0.025, 只是冬季 ( 即第天第天 ) 存在一定的偏差,αbsa 与的最大偏差分别为和图年 D105 站 GLASS 地表反照率与地面观测结果的比较 Fig.3 ComparisonbetweentheGLASSalbedosandground observedalbedosatd105stationin2003 图 4 给出的是 2003 年 NPAM 站观测的地表反照率与站点对应像元 GLASS 地表反照率的对比第 1 9 天开始的 8d,NPAM 站对应像元的 GLASS 黑空反照率 (αbsa) 白空反照率 (αwsa) 和地面观测结果 () 存在显著差异 ; 第天开始的 8d, αbsa 和也存在一定的差异, 最大偏差分别为和 0.084; 其余 8d 的 αbsa 和的一致性都非常好分析以 NPAM 站对应像元及其相邻像元大气校正后的地表反射率, 在第 1 天开始的 8d 期间,GLASS 地表反照率终级产品时域滤波窗口的 Terra 和 Aqua MODIS 数据有云覆盖的天数均为 11d 类似地, 在第 9 天开始的 8d,GLASS 地表反照率终级产品时域滤波窗口的 Terra MODIS 数据有云覆盖的天数为 9d,Aqua MODIS 数据有云覆盖的天数为 10 d 因此, 第 1 9 天开始的 8 d GLASS 黑空反照率 (αbsa) 白空反照率 (αwsa) 都超过图年 NPAM 站 GLASS 地表反照率与地面观测结果的比较 Fig.4 ComparisonbetweentheGLASSalbedosandground observedalbedosatthenpamstationin2003

6 陈爱军等 : 青藏高原 GLASS 地表反照率产品精度分析 , 显然受到了云覆盖的影响至于第天开始的 8d, 可能与 NPAM 站夏季非均匀生长的草甸有关, 导致地表反照率地面观测结果偏低, 而 GLASS 地表反照率偏高只有获得更多的观测数据, 才能进一步确定第 1 9 天和第天开始的 8dNPAM 站对应像元的 GLASS 地表反照率与地面观测结果存在差异的原因分析 BJ D105 和 NPAM 三站对应像元的 GLASS 黑空反照率 (αbsa) 白空反照率 (αwsa) 与地面观测结果的均方根误差 (σbsa) 和 (σwsa), 结果表明 ( 表 2), 除 D105 站对应像元的黑空反照率 (αbsa) 白空反照率 (αwsa) 与地面观测结果 8d 平均地表反照率 () 的均方根误差都小于 0.02 之外,BJ 和 NPAM 站则都超过 0.07 综合这 3 个站对应像元的 GLASS 地表反照率, 分析黑空反照率 (αbsa) 白空反照率 (αwsa) 与地面观测结果 () 的均方根误差, 则前者为 , 后者为表 2 GLASS 地表反照率与地面观测结果的均方根误差 Table2 Theroot mean squareerrorsbetweentheglass albedosandground observedalbedos 站点 BJ D105 NPAM 全部 σbsa σwsa 考虑到 GLASS 地表反照率可能受到云覆盖的影响而偏高, 地面观测结果可能因为局地积雪覆盖而偏高, 如果只考虑 GLASS 地表反照率地表反照率地面观测结果中低于 0.4 的周期, 除 D105 站没有任何变化外,BJ 站对应像元的 GLASS 黑空反照率 (αbsa) 白空反照率 (αwsa) 与地面观测结果 8d 平均地表反照率 () 的均方根误差显著减小, 分别为和 ;NPAM 站也显著减小,αbsa 与的均方根误差分别为和同样, 综合 3 个站数据的分析结果显示, 此时 αbsa 与的均方根误差分别为和 ( 表 3) 表 3 GLASS 地表反照率和地表反照率地面观测结果均小于 0.4 时的均方根误差 Table3 Theroot mean squareerorsbetweentheglass albedosandground observedalbedosinthecasethatboth theresultsoftheglassalbedosandground observed albedosaresmalerthan0.4 BJ D105 NPAM 全部 σbsa σwsa 小结地表反照率对地表辐射能量收支有着重要影响, 是气候和陆面过程模式的重要参数 GLASS 地表反照率产品是目前中国唯一自主开发的高精度时空连续一致的全球地表反照率产品, 对于研究青藏高原地表辐射能量收支和地气相互作用具有重要的应用价值应用青藏高原 BJ D105 和 NPAM 三站 2003 年自动气象站观测的下行和上行太阳总短波辐射资料, 经过处理后得到与 GLASS 地表反照率产品对应的地表反照率地面观测结果, 并与站点对应像元的 GLASS 黑空反照率白空反照率进行比较, 分析了青藏高原地区 GLASS 地表反照率产品的精度, 结果表明 : (1)GLASS 黑空反照率白空反照率与地表反照率地面观测结果的总体变化趋势基本一致, 能够有效地反映实际地表状态的变化情况由于云覆盖的影响,GLASS 地表反照率可能高估实际的地表反照率同时, 积雪可能导致 GLASS 地表反照率低于地面观测结果 (2)GLASS 黑空反照率白空反照率在 D105 站与地表反照率地面观测结果的一致性较好, 均方根误差不超过 0.02 BJ 站和 NPAM 站因为云覆盖或者局地积雪的影响,GLASS 黑空反照率白空反照率与地面观测结果的均方根误差均超过 0.05 如果消除这两方面因素的影响, 则显著降低,BJ 站对应像元的 GLASS 黑空反照率白空反照率与地面观测结果的均方根误差分别为和 ; NPAM 站的均方根误差分别为和 ; 综合 3 站数据的分析结果, 均方根误差分别为和 (3) 由于局地积雪和云覆盖对 GLASS 地表反照率产品的精度影响较大如何正确地分析和评价 GLASS 地表反照率产品在局地积雪和云覆盖状态下的精度, 还需要进一步分析和研究致谢 : 感谢中国科学院青藏高原研究所马耀明研究员提供 CAMP/Tibet 试验地面辐射观测数据参考文献陈爱军, 卞林根, 刘玉洁等应用 MODIS 数据反演青藏高原地区地表反照率. 南京气象学院学报,32(2): Chen AJ,BianLG,LiuYJ,etal.2009.UsingMODISdatatore trievealbedooverthe Qinghai TibetPlateau.J NanjingInst

7 1120 犃犮狋犪犕犲狋犲狅狉狅犾狅犵犻犮犪犛犻狀犻犮犪气象学报 2015,73(6) Meteor,32(2): (inchinese) 梁顺林, 张晓通, 肖志强等全球陆表特征参量 (GLASS) 产品 : 算法验证与分析. 北京 : 高等教育出版社,33 72.LiangS L,ZhangX T,XiaoZ Q,etal.2014.GlobalLandSurface Satelite(GLASS)Products:ValidationandAnalysis.Beijing: HigherEducationPress,33 72 (inchinese) 刘强, 瞿瑛, 王立钊等.2012.GLASS 陆表反照率产品使用手册. 北京 : 北京师范大学.LiuQ,QuY,WangLZ,etal GLASS GlobalLandSurfaceBroadbandAlbedoProduct:User's Guide.Beijing:BeijingNormalUniversity(inChinese) 齐文栋, 刘强, 洪友堂种反演算法的地表反照率遥感产品对比分析. 遥感学报,18(3): QiW D,LiuQ,Hong Y T.2014.Comparisonanalysisbasedondiferentinverseal gorithmsofsurfacealbedoproducts.jremotesens,18(3): (inchinese) 王鸽, 韩琳地表反照率研究进展. 高原山地气象研究,30 (2):79 83.WangG,HanL.2010.Progressintheresearchof surfacealbedo.plateau Mountain MeteorRes,30(2):79 83 (inchinese) 王健, 孙林使用 FLUXNET 数据验证 GLASS 地表反照率. 山东建筑大学学报,28(1):50 54.WangJ,Sun L GLASSsurfacealbedovalidationbyusingFLUXNET data.j ShandongJianzhuUniv,28(1):50 54 (inchinese) 王立钊, 郑学昌, 孙林等利用 LandsatTM 数据和地面观测数据验证 GLASS 反照率产品. 遥感学报,18(3): WangL Z,Zheng X C,Sun L,etal.2014.Validationof GLASSalbedoproductthrough LandsatTM dataandground measurements.jremotesens,18(3): (inchinese) 叶笃正, 高由禧青藏高原气象学. 北京 : 科学出版社, YeD Z,GaoY X.1979.TibetanPlateau Meteorology.Bei jing:sciencepress, (inchinese) DickinsonRE.1983.Landsurfaceprocessesandclimate surfaceal bedosandenergybalance.advgeophys,25: LewisP,BrockmannC,DanneO,etal.2011.GlobAlbedo:Algo rithmtheoreticalbasisdocument,323pp.availablefrom htp: 1.pdf LiangSL,ZhangX T,XiaoZQ,etal.2014.Shortwavealbedo GlobalLand Surface Satelite (GLASS)Products.Springer BriefsinEarthSciences.Berlin:SpringerInternationalPublish ing,33 72,doi: / _3 LiuNF,LiuQ,WangLZ,etal.2011.Atemporalfilteringalgo rithmtoreconstructdailyalbedoseriesbasedonglassalbedo product Proceedingsof2011IEEE InternationalGeoscience andremotesensingsymposium.vancouver,bc:ieee, LiuNF,LiuQ,WangLZ,etal.2013.Astatistics basedtemporal filteralgorithmto mapspatiotemporalycontinuousshortwave albedofrom MODISdata.HydrolEarthSystSci,17(6): LiuQ,WangLZ,QuY,etal.2013.Preliminaryevaluationofthe long term GLASSalbedoproduct.IntJDigitEarth,6(1):69 95 SchaafCB,GaoF,StrahlerA H,etal.2002.Firstoperational BRDF,albedoandnadirreflectanceproductsfrom MODIS.Re motesensenviron,83(1 2): SchaafCB,LiuJC,GaoF,etal.2011.MODISalbedoandreflec tanceanisotropyproductsfrom AquaandTerra Ramachand ranb,justicec O,Abrams M J.LandRemoteSensingand GlobalEnvironmentalChange:NASA'sEarth ObservingSys tem andthescienceofasterand MODIS,RemoteSensing and DigitalImage Processing SeriesVolume,11.New York: Springer Verlag, SelersPJ,MeesonB W,HalFG,etal.1995.Remotesensingof theland surfaceforstudies ofglobalchange: Models algo rithms experiments.remotesensenviron,51(1):3 26

殐殐 27 1 Vol.27 No ImagingScienceandPhotochemistry Jan.,2009 檭檭檭檭檭檭檭檭檭檭檭檭檭檭殐殐 *O! "# #B) 1,CDE 1, 2, 1 (1., ;2., ) $ %:3 J = 2J4

殐殐 27 1 Vol.27 No ImagingScienceandPhotochemistry Jan.,2009 檭檭檭檭檭檭檭檭檭檭檭檭檭檭殐殐 *O! # #B) 1,CDE 1, 2, 1 (1., ;2., ) $ %:3 J = 2J4 殐殐 27 1 Vol.27 No.1 2009 1 ImagingScienceandPhotochemistry Jan.,2009 檭檭檭檭檭檭檭檭檭檭檭檭檭檭殐殐 *O! "# #B) 1,CDE 1, 2, 1 (1., 200092;2., 200092) $ %:3 J = 2J4 ',V T, _` * C?. + SU L 3 H4 E '.< 4 893?, V 9 P>?.

doi: /qxxb 气象学报 青藏高原犌犔犃犛犛地表反照率产品精度分析 陈爱军 1,2 胡慎慎 1,2 卞林根 3 刘玉洁 4 CHEN Aijun 1,2 HUShenshen 1,2 BIANLingen 3 LIU Yujie 4 1. 中国气象局气溶胶与云降水

doi: /qxxb 气象学报青藏高原犌犔犃犛犛地表反照率产品精度分析陈爱军 1,2 胡慎慎 1,2 卞林根 3 刘玉洁 4 CHEN Aijun 1,2 HUShenshen 1,2 BIANLingen 3 LIU Yujie 4 1. 中国气象局气溶胶与云降水