青藏高原地区MODIS反照率的精度分析

陈爱军，梁学伟，卞林根，刘玉洁，王飞，朱小祥

(1.南京信息工程大学气象灾害省部共建教育部重点实验室，江苏南京210044;2.南京信息工程大学大气物理学院，江苏南京210044;3.中国气象科学研究院，北京100081;4.国家卫星气象中心，北京100081;5.国家气象中心，北京100081)

青藏高原地区MODIS反照率的精度分析

陈爱军1，2，梁学伟1，2，卞林根3，刘玉洁4，王飞1，2，朱小祥5

(1.南京信息工程大学气象灾害省部共建教育部重点实验室，江苏南京210044;2.南京信息工程大学大气物理学院，江苏南京210044;3.中国气象科学研究院，北京100081;4.国家卫星气象中心，北京100081;5.国家气象中心，北京100081)

应用2002—2004年青藏高原CAMP/Tibet试验期间4个地面站点的反照率观测结果定量分析Terra MODIS 1 km分辨率短波SW波段(0.3～5.0 μm)反照率全反演结果和当量反演结果的精度。对于全反演结果，黑空反照率、白空反照率与地面观测结果的均方根差分别为0.018 7和0.016 8;对于当量反演结果，黑空反照率、白空反照率与地面观测结果的均方根差分别为0.076 6和0.076 1。综合全反演结果和当量反演结果，则黑空反照率、白空反照率与地面观测结果的均方根差分别为0.067 9和0.067 5。当地面观测结果与MODIS反照率当量反演结果均为“无雪”状态时，黑空反照率、白空反照率与地面观测结果的均方根差分别为0.035 2和0.036 4;当地面观测结果为“积雪”状态，MODIS反照率当量反演结果为“无雪”状态时，黑空反照率、白空反照率与地面观测结果的均方根差分别高达0.155 6和0.154 1。

反照率;验证;MODIS;青藏高原

0 引言

地表反照率是地表反射与入射太阳辐射通量密度之比。地表反照率对地表辐射能量收支、地—气相互作用，以及全球气候变化均有重要影响(Dickinson，1983;Henderson-Sellers and Wilson，1983;Sellers et al.，1995)，是气候和陆面过程模式的重要参数之一(Dickinson et al.，1993;Sellers et al.，1996;Bonan，1998;Blackmon et al.，2001;Dai et al.，2003)。气候和陆面过程模式对地表反照率的精度要求很高，要求其绝对精度应达到0.02～0.05(Henderson-Sellers and Wilson，1983;Sellers et al.，1995)。

卫星遥感已成为连续获取大范围乃至全球高时空分辨率和高精度地表反照率的一条重要途径(Pinty et al.，2000a，2000b;Schaaf et al.，2002，2011;Govaerts et al.，2010;Wagner et al.，2010)。早期，卫星遥感反演地表反照率主要采用统计反演方法(Brest and Goward，1987;Wydick et al.，1987;祝昌汉等，1990;Li and Garand，1994;方宗义等，1996)。自Terra MODIS(MODerate resolution Imaging Spectraradiometer)发射以来，美国采用半经验的物理反演方法获得MODIS全球二向反射分布函数BRDF(BidirectionalReflectanceDistribution Function)和反照率产品(Schaaf et al.，2002)。该产品最初仅由Terra MODIS数据生成，反演周期为16 d，空间分辨率为1 km，产品编号为MOD43(Schaaf et al.，2002，2011)。

MOD43采用RossThick-LiSparse-R核驱动的BRDF模型反演BRDF参数并计算黑空反照率BSA(black sky albedo)和白空反照率WSA(white sky albedo)(Schaaf et al.，2002，2011;陈爱军等，2009)。在MOD43反演过程中，如果16 d期间累积的高质量多角度晴空观测总数n≥7，就采用“全反演法(the full inversion)”，由最小二乘法获得拟合多角度观测数据最佳的BRDF参数;如果3≤n＜7，或者全反演时的拟合效果不佳，就采用“当量反演法(the magnitude inversion)”，利用多角度晴空观测数据修正原型(archetype)BRDF参数数据库提供的先验(a priori)信息，得到与实际情况相符的BRDF参数(Schaaf et al.，2002，2011;陈爱军等，2009)。由于观测数据的数量有限，或者质量欠佳，反演过程又依赖于先验信息，当量反演结果质量自然也会受到影响，其精度可能会不及全反演结果。最新的MODIS全球BRDF和反照率产品已联合Terra和Aqua双星MODIS数据生成，反演周期缩短为8 d，空间分辨率也提高到500 m，产品编号为MCD43(Roman et al.，2009;Schaaf et al.，2011)。

MODIS反照率产品受到科学界的广泛关注(Zhou et al.，2003;Tian et al.，2004;Myhre et al.，2005;Lawrence and Chase，2007;Morcrette et al.，2008)，尤其是精度。人们采用多种途径对MODIS反照率产品的精度进行分析和评价，例如:将反照率的地面观测结果直接与其比较(Jin et al.，2003b;Roesch，2004;Wang et al.，2004;Stroeve et al.，2005;焦子锑等，2005;Salomon et al.，2006;Samain et al.，2008;Liu et al.，2009;余予等，2010)、将反照率的地面观测结果与高分辨率遥感反演结果融合后再与其比较(Liang et al.，2002;Susaki et al.，2007;Roman et al.，2009)、将反照率机载观测结果与其比较(焦子锑等，2005;Knobelspiesse et al.，2008)，以及应用卫星遥感反演结果与其比较(Jin et al.，2003b)。虽然研究结果表明MODIS反照率能够满足气候和陆面过程模式对反照率的精度要求，但是，上述研究主要针对MODIS反照率全反演结果，很少涉及MODIS反照率当量反演结果，即使涉及，也只是定性分析(Jin et al.，2003a;Stroeve et al.，2005;Liu et al.，2009;Roman et al.，2009)。最近，陈爱军等(2012)应用“对照反演法”对中国地区相同时空条件下MODIS反照率全反演结果与当量反演结果的差异做了定量分析。事实上，全球陆地约一半以上的区域需要通过当量反演获得MODIS反照率(Jin et al.，2003a)。因此，研究MODIS反照率当量反演结果的精度具有十分重要的意义，也有助于正确分析和评价MODIS反照率产品的精度。

青藏高原是世界上平均海拔最高、面积最大、地形最为复杂的高原，素有“地球第三极”之称。青藏高原高大的地形特征及其对太阳辐射的吸收和反射，不仅影响了该地区地—气相互作用过程(尤其是能量与水分循环过程)，也对亚洲季风、东亚大气环流及全球气候变化有着极大的影响(叶笃正和高由禧，1979;章基嘉等，1988)。MODIS反照率无疑为研究青藏高原地区的地表反照率，以及该地区的地—气相互作用过程提供了一条重要途径。然而，由于云和积雪的影响，青藏高原地区反照率当量反演结果的比例甚至超过全反演结果的比例(陈爱军等，2009;Chen et al.，2007a，2007b)。这增加了青藏高原地区MODIS反照率产品精度的不确定性。Wang et al.(2004)用青藏高原改则自动气象站近3 a的观测结果分析了MOD43B3的精度。余予等(2010)利用2007年5月—2008年8月青藏高原纳木错站无雪期晴空条件下的观测结果，分析了由Terra和Aqua双星MODIS数据生成的反照率产品MCD43B3的精度，但上述研究均未涉及当量反演结果。最近，陈爱军等(2012)应用“对照反演法”定量分析了中国地区相同时空条件下MODIS反照率全反演结果和当量反演结果的差异。

本文应用2002—2004年青藏高原地区4个自动气象站的下行和上行太阳总短波辐射观测资料，并经处理后获得反照率地面观测结果，据此分别定量分析Terra MODIS反照率产品的全反演结果、当量反演结果的精度。

1 资料

1.1 地面观测资料

研究所用地面观测资料来自2001—2005年青藏高原CAMP/Tibet国际合作科学观测试验期间D105、Amdo、MS3478和BJ等4个站点的下行和上行太阳总短波辐射观测资料。这4个站点依次由北向南排列(表1)，站点所处位置比较开阔，能够代表藏北高原大范围的地形和下垫面特征，观测资料具有较好的代表性，其中:D105站位于多年冻土地区，地表为淤泥土壤和砂质土壤;Amdo站周围地面覆盖着高约5 cm的高原草甸;MS3478站的地表为含有细石块的砂黏土，夏季不均匀地生长着高约15 cm的草甸;BJ站地面覆盖着高约5 cm、稀疏的高原草甸，土壤由含有碎石的砂壤土构成。另外，CAMP/Tibet试验前所有观测仪器都经过严格的标定和校准(马伟强等，2004;李英和胡泽勇，2006;马耀明等，2006;陈学龙等，2008)，确保观测结果的准确性和可比性。

表1 地面观测站点的基本信息Table 1 The basic information of ground-based stations

研究所用的下行和上行太阳短波总辐射资料来自以上4个站点的自动气象站观测资料。1)D105站:2002—2003年、2004年1月1日—10月14日，时间分辨率均为1 h。2)Amdo站:2002年1月1日—8月21日，时间分辨率为10 min;2004年4月22日—6月10日、8月12日—12月31日，时间分辨率为30 min。3)MS3478站:2002—2003年、2004年1月1日—8月12日，时间分辨率均为1 h。4)BJ站:2002年，时间分辨率为1 h;2003年、2004年6月7日—8月23日，时间分辨率均为10 min。

1.2 MODIS反照率数据

研究采用美国地质调查局USGS(United States Geological Survey)地球资源观测与科学EROS(Earth Resource Observation and Science)数据中心提供的L3级全球1 km正弦投影网格(L3 Global 1 km SIN Grid)Terra MODIS 16 d反照率产品MOD43(V004)。该产品已经过验证，适合科学研究，主要包括:1)MODIS 1—7通道和可见光(0.3～0.7 μm)、近红外(0.7～5.0 μm)、短波(0.3～5.0 μm)3个宽波段的黑空反照率BSA、白空反照率WSA，分别以1 km和0.05°气候网格两种空间分辨率提供。2)MODIS 1—7通道星下点BRDF校正后的反射率(NBAR，Nadir BRDF-Adjusted Reflectances)，空间分辨率为1 km。3)MODIS 1—7波段和可见光、近红外、短波3个宽波段的BRDF参数，空间分辨率为1 km。4)数据质量信息，包括反照率反演质量、地表类型、当地正午太阳天顶角和积雪状态等信息(Schaaf et al.，2002)。

MODIS陆地正弦投影网格(http://landdb1.nascom.nasa.gov/developers/is_tiles/is_grid.html)编号h25v05区正好覆盖CAMP/Tibet试验区，故选用该区域2002—2004年短波SW波段1 km分辨率的黑空反照率BSA和白空反照率WSA进行研究。

2 地面观测数据处理

为了获得与Terra MODIS反照率产品具有相同周期的地面观测结果，需要对自动气象站观测的下行和上行太阳短波总辐射资料进行适当地处理:1)根据下行太阳短波总辐射不低于上行太阳短波总辐射的原则，剔除明显有误的观测数据。2)将10或30 min的采样数据按小时平均，得到小时平均的下行和上行太阳短波总辐射数据。3)以下行太阳短波总辐射不小于600 W/m2为“晴空”标准(马伟强等，2004)，以上行与下行太阳短波总辐射之比计算“晴空”条件的小时平均地表反照率;否则，记为0。4)以非0的小时平均地表反照率计算“晴空”日平均地表反照率。5)以“晴空”日的平均地表反照率计算16 d平均地表反照率。

考虑到积雪对地表反照率的影响较大，结合CAMP/Tibet试验区的地表类型，为了与“积雪”条件下的MODIS反照率反演算法一致(Schaaf et al.，2002，2011;陈爱军等，2009)，计算16 d平均地表反照率时，如果日平均地表反照率低于0.4，则视为“无雪”日反照率，否则视为“积雪”日反照率。如果“积雪”日反照率不少于8 d，就以“积雪”日反照率平均值作为16 d平均地表反照率;反之，以“无雪”日反照率平均值作为16 d平均地表反照率。

图1是根据BJ站2002年第17—32、33—48天的下行和上行太阳短波总辐射资料计算的小时平均、日平均和16 d平均地表反照率。可以看出:第20天的日平均地表反照率小于0.4(约为0.3)，表明当天的反照率观测结果为“无雪”日的反照率;第28天的小时平均反照率虽然接近0.7，但当天的下行太阳短波总辐射均小于600 W/m2，为非“晴空”反照率观测结果，故当天的日平均地表反照率记为0。除第20、28天之外，其余14 d的日平均地表反照率均高于0.4，有的甚至超过0.7，表明均为“积雪”日反照率(图1a)。第17—32天期间，“积雪”日明显多于“无雪”日，应以“积雪”日的反照率平均值作为该16 d的反照率观测结果，约为0.66。从第33天起，BJ站的日平均反照率逐渐减小，第38天时已降至0.4左右，第43—48天期间已稳定维持在0.25左右(图1b)。第33—48天期间，只有第33—37天的日平均地表反照率不低于0.4，“积雪”日为5 d，明显少于“无雪”日，应以11 d“无雪”日的反照率平均值作为该16 d的反照率观测结果，约为0.28。

图1 BJ站2002年第17—32天(a)和第33—48天(b)的反照率地面观测结果Fig.1 Ground-observed albedos at BJ station during the period of(a)days 17—32 and(b)days 33—48 of 2002

3 结果分析

MODIS黑空反照率BSA和白空反照率WSA分别代表完全晴空和完全阴天条件下的地表反照率，实际地表反照率要根据天空漫射光比例由两者内插而得(Schaaf et al.，2002，2011;陈爱军等，2009)。由于天空漫射光比例与气溶胶光学厚度有关，获得地面站点的气溶胶光学厚度观测数据又比较困难，故以反照率地面观测结果与站点对应像元的MODIS黑空反照率BSA、白空反照率WSA直接进行比较。如果黑空反照率BSA或者白空反照率WSA中的任意一个与地面观测结果的绝对偏差超过0.05，那么，根据天空漫射光比例计算的实际地表反照率与地面观测结果的绝对偏差也必然超过0.05。

图2为BJ站2003年22个16 d(每年第23个16 d的实际天数不足16 d，不予考虑)的反照率地面观测结果和MODIS黑空反照率BSA、白空反照率WSA。可以看出:1)第1天开始的16 d，反照率地面观测结果约为0.70，明显受到积雪的影响。虽然MODIS反照率为当量反演结果，但黑空反照率BSA、白空反照率WSA均接近0.75，与地面观测结果非常接近。2)第17天开始的16 d，反照率地面观测结果约为0.60，表明仍然受到积雪的影响;但是，MODIS黑空反照率BSA、白空反照率WSA均介于0.20和0.25之间，与地面观测结果的差异十分显著。这可能与观测站点附近的局地积雪有关:第1个16 d期间，站点周边大范围积雪，MODIS能够较好地探测到地面积雪信息，反照率反演结果能够较好地反映地表状态，与地面站点观测结果的差异较小;第2个16 d期间，地面积雪逐渐消融，虽然站点(甚至只是观测仪器)附近存在局部积雪，反照率地面观测结果较大，但MODIS已无法探测到局地积雪信息，反照率反演结果自然偏低，与地面观测结果产生较大的差异。3)在22个16 d中，只有5个16 d的MODIS反照率是全反演结果，即第33、273、289、305、321天开始的16 d;第97天开始的16 d没有MODIS反照率反演结果;通过当量反演获得MODIS反照率的16 d多达16个。这主要与当地常年有云覆盖有关(尤其在雨季)，很难在16 d反演周期内获得不少于7个高质量多角度晴空观测数据并通过全反演成功获得BRDF参数(Chen et al.，2007a，2007b;陈爱军等，2007;陈爱军等，2009)。此外，也与冬春季地面积雪的影响有关。

图3为2003年BJ站反照率观测结果与MODIS黑空反照率BSA、白空反照率WSA的绝对偏差。可以看出:第17天开始的16 d，反照率地面观测结果与MODIS黑空反照率BSA、白空反照率WSA的绝对偏差都非常大，达到0.35以上，其余16 d的绝对偏差都不超过0.1。进一步分析表明，除第33天开始的16 d之外，其余16 d反照率地面观测结果与MODIS黑空反照率BSA、白空反照率WSA的绝对偏差均不超过0.05。

图2 2003年BJ站反照率地面观测结果与Terra MODIS反照率的比较(实心点表示全反演结果;空心点表示当量反演结果)Fig.2 Comparison between ground-observed albedos and Terra MODIS albedos at BJ station in 2003(the solid dots represent the full inversion results，and the hollow dots do the magnitude inversion results)

综合分析图2和图3可以看出:除第33天开始的16 d之外，其余4个通过全反演获得MODIS反照率的16 d中，MODIS黑空反照率BSA、白空反照率WSA与地面观测结果的绝对偏差都小于0.05;16个通过当量反演获得MODIS反照率的16 d中，MODIS黑空反照率BSA、白空反照率WSA与地面观测结果的绝对偏差均小于0.05。也就是说，在获得MODIS反照率反演结果的21个16 d中，有20个与反照率地面观测结果的绝对偏差小于0.05，满足绝对精度不超过0.05的要求，所占比例达到95%以上。总体而言，2003年BJ站MODIS反照率反演结果的精度较高。

2002—2004年，与D105、Amdo、MS3478和BJ站的反照率地面观测结果对应，且站点对应像元的MODIS反照率为全反演结果、当量反演结果的分别有41和138个16 d。结果表明:MODIS黑空反照率BSA全反演结果、白空反照率WSA反照率全反演结果与地面观测结果的均方根差分别为0.018 7和0.016 8，可以满足绝对精度不低于0.02的要求;MODIS黑空反照率BSA当量反演结果、白空反照率WSA当量反演结果与地面观测结果的均方根差则高达0.076 6和0.076 1，无法满足绝对精度不低于0.05的最低要求;如果综合MODIS反照率全反演结果和当量反演结果，则黑空反照率BSA、白空反照率WSA与地面观测结果的均方根差为0.067 9和0.067 5，仍然无法满足绝对精度不低于0.05的最低要求(表2)。总之，MODIS反照率全反演结果的精度较高，完全可以满足气候和陆面过程模式的反照率精度要求;当量反演结果的精度远低于全反演结果的精度，无法满足目前气候和陆面过程模式的反照率精度要求。

分析反照率地面观测结果和MODIS反照率当量反演结果的积雪状态可知，两者均为“积雪”状态的只有BJ站2003年第1天开始的16 d和D105站2004年第17天开始的16 d。此外，地面观测结果

图3 2003年BJ站反照率地面观测结果与Terra MODIS黑空反照率、白空反照率的绝对偏差Fig.3 The absolute biases between ground-observed albedos and Terra MODIS black-sky albedos(white-sky albedos)at BJ station in 2003

表2 反照率地面观测结果与MODIS反照率的均方根差Table 2 Root-mean-square errors between ground-observed albedos and MODIS albedos

为“积雪”状态、MODIS反照率当量反演结果为“无雪”状态的现象比较常见，没有出现地面观测结果为“无雪”状态、MODIS反照率当量反演结果为“积雪”状态的现象。如果按照获得反照率地面观测结果的积雪状态，将其分为“无雪”组和“积雪”组，“无雪”组地面观测结果与对应的黑空反照率BSA、白空反照率WSA的均方根差为0.035 2和0.036 4，而“积雪”组的均方根差则为0.155 6和0.154 1(表3)。如果综合MODIS反照率全反演结果和“无雪”组的当量反演结果，则地面观测结果与黑空反照率BSA、白空反照率WSA的均方根差为0.032 0和0.032 8。以上分析结果，一方面说明地面观测站点附近的局地积雪对反照率地面观测结果有影响，不利于正确分析和评价MODIS反照率当量反演结果的精度，另一方面也说明在大范围长期无雪区，MODIS反照率当量反演结果精度较高，可以满足地表反照率的绝对精度要求。当然，如何应用地面观测数据正确地分析和评价MODIS反照率当量反演结果的精度，仍然有待深入研究。

表3 不同积雪状态的反照率地面观测结果与MODIS反照率当量反演结果的均方根差Table 3 Root-mean-square errors between ground-observed albedos of different snow conditions and MODIS albedos derived with the magnitude inversion

4 结论

地表反照率是影响地表辐射能量收支的一个重要参数，而且气候和陆面过程模式对其精度要求很高。本文应用青藏高原CAMP/Tibet试验期间D105、Amdo、MS3478和BJ等4个站点的自动气象站观测的下行和上行太阳总短波辐射资料，经过处理后得到16 d平均的反照率地面观测结果，与美国Terra MODIS 1 km分辨率反照率产品MOD43提供的短波SW波段黑空反照率BSA、白空反照率WSA直接进行比较，定量分析了青藏高原地区MODIS反照率全反演结果和当量反演结果的精度，结果表明:

1)MODIS反照率全反演结果能够满足绝对精度不低于0.02的要求，黑空反照率BSA、白空反照率WSA与地面观测结果的均方根差分别为0.018 7和0.016 8;MODIS反照率当量反演结果的精度明显低于全反演结果，无法满足绝对精度不低于0.05的最低要求，黑空反照率BSA、白空反照率WSA与地面观测结果的均方根差分别为0.076 6和0.076 1;综合MODIS反照率全反演结果和当量反演结果，则黑空反照率BSA、白空反照率WSA与地面观测结果的均方根差分别为0.067 9和0.067 5，无法满足绝对精度不低于0.05的要求。

2)地面观测站点附近的局地积雪会影响反照率地面观测结果，不利于正确分析和评价MODIS反照率当量反演结果的精度。如果地面观测结果为“无雪”状态，它们与MODIS黑空反照率BSA当量反演结果、白空反照率WSA当量反演结果的均方根差分别为0.035 2和0.036 4;如果地面观测结果为“积雪”状态，它们与MODIS黑空反照率BSA当量反演结果、白空反照率WSA当量反演结果的均方根差分别为0.155 6和0.154 1。

总之，MODIS反照率全反演结果的精度较高，能够满足气候和陆面过程模式对地表反照率绝对精度的要求。虽然MODIS反照率当量反演结果的精度低于全反演结果，某些特殊条件下(如长期大范围无雪)的反照率当量反演结果的精度仍然较高。由于局地积雪对地表反照率的影响，提高卫星遥感积雪的空间分辨率，有助于提高积雪判识精度，并提高MODIS反照率反演结果的精度。最新的由Terra和Aqua双星MODIS数据生成的全球反照率产品MCD43，不仅空间分辨率提高至500 m，反演周期也缩短为8 d(Roman et al.，2009;Schaaf et al.，2011)，相信有利于减少积雪对反照率反演结果的影响，进一步提高MODIS反照率的精度。

致谢:本文所用CAMP/Tibet试验地面辐射观测数据由中国科学院青藏高原研究所马耀明研究员提供，深表感谢!

陈爱军，卞林根，刘玉洁，等.2007.MODIS反演多云地区地表反照率的一种新方法及其反演试验［J］.大气科学，31(5):855-862.

陈爱军，卞林根，刘玉洁，等.2009.应用MODIS数据反演青藏高原地区地表反照率［J］.南京气象学院学报，32(2):222-229.

陈爱军，刘玉洁，卞林根，等.2012.中国地区MODIS反照率两种反演结果的比较［J］.气象学报，70(5):1119-1127.

陈学龙，马耀明，李茂善，等.2008.藏北地区近地层大气和土壤特征量分析［J］.高原气象，27(5):941-948.

方宗义，刘玉洁，林曼芸.1996.青藏高原地表反照率计算研究［J］.气象学报，54(5):580-589.

焦子锑，王锦地，谢里欧，等.2005.地面和机载多角度观测数据的反照率反演及对MODIS反照率产品的初步验证［J］.遥感学报，9(1):64-72.

李英，胡泽勇.2006.藏北高原地表反照率的初步研究［J］.高原气象，25(6):1034-1041.

马伟强，马耀明，王介民，等.2004.藏北高原地面辐射收支的初步分析［J］.高原气象，23(3):348-352.

马耀明，姚檀栋，王介民.2006.青藏高原能量和水循环试验研究——GAME/Tibet与CAMP/Tibet研究进展［J］.高原气象，25(2):344-351.

叶笃正，高由禧.1979.青藏高原气象学［M］.北京:科学出版社:7-9.

余予，陈洪滨，夏祥鳌，等.2010.青藏高原纳木错站地表反照率观测与MODIS资料的对比分析［J］.高原气象，29(2):260-267.

章基嘉，朱抱真，朱福康，等.1988.青藏高原气象学进展［M］.北京:科学出版社:14-89.

祝昌汉，朱福康，刘玉洁.1990.青藏高原晴空行星反照率与地表反照率的研究［J］.科学通报，35(20):1563-1565.

Blackmon M，Boville B，Bryan F，et al.2001.The Community Climate System Model［J］.Bull Amer Meteor Soc，82(11):2357-2376.

Bonan G B.1998.The land surface climatology of the NCAR land sur-face model(LSM 1.0)coupled to the NCAR Community Climate Model(CCM3)［J］.J Climate，11:1307-1326.

Brest C L，Goward S N.1987.Deriving surface albedo measurements from narrow band satellite data［J］.Int J Remote Sens，8(3):351-367.

Chen Aijun，Bian Lingen，Liu Yujie，et al.2007a.Retrieving summer albedo over the Tibetan plateau with MODIS data［C］//Proceedings of SPIE，MIPPR 2007:Remote Sensing and GIS Data Processing and Applications;and Innovative Multispectral Technology and Applications.Wuhan，China.

Chen Aijun，Bian Lingen，Liu Yujie.2007b.Deriving albedo over cloudy areas with composite inversion［C］//Proceedings of SPIE，Geoinformatics 2007，Geospatial Information Science.Nanjing，China.

Dai Y，Zeng X，Dickinson R E，et al.2003.The common land model［J］.Bull Amer Meteor Soc，84(8):1013-1023.

Dickinson R E.1983.Land surface processes and climate:Surface albedos and energy balance［J］.Adv Geophys，25:305-353.

Dickinson R E，Henderson-Sellers A，Kennedy P J.1993.Biosphere-Atmosphere Transfer Scheme(BATS)Version 1e as coupled to the NCAR community model［R］//NCAR Tech Note NCAR/TN-387+STR，National Center for Atmosphere Research，Boulder，Colorado.

Govaerts Y M，Wagner S，Lattanzio A，et al.2010.Joint retrieval of surface reflectance and aerosol optical depth from MSG/SEVIRI observations with an optimal estimation approach:1.Theory［J］.J Geophys Res，115，D02203.doi:10.1029/2009JD011779.

Henderson-Sellers A，Wilson M F.1983.Surface albedo data for climatic modeling［J］.Rev Geophys，21:1743-1778.

Jin Yufang，Schaaf C B，Woodcock C E，et al.2003a.Consistency of MODIS surface bidirectional reflectance distribution function and albedo retrievals:1.Algorithm performance［J］.J Geophys Res，108(D5)，4158.doi:10.1029/2002JD002803.

Jin Yufang，Schaaf C B，Woodcock C E，et al.2003b.Consistency of MODIS surface bidirectional reflectance distribution function and albedo retrievals:2.Validation［J］.J Geophys Res，108(D5)，4159.doi:10.1029/2002JD002804.

Knobelspiesse K D，Cairns B，Roman M，et al.2008.Surface BRDF estimation from an aircraft compared to MODIS and ground estimates at the Southern Great Plains site［J］.J Geophys Res，113，D20105.doi:10.1029/2008JD010062.

Lawrence P J，Chase T N.2007.Representing a new MODIS consistent land surface in the Community Land Model(CLM3.0)［J］.J Geophys Res，112，G01023.doi:10.1029/2006JG000168.

Li Zhanqing，Garand L.1994.Estimation of surface albedo from space:A parameterization for global application［J］.J Geophys Res，99(D4):8335-8350.

Liang S，Fang H，Chen J，et al.2002.Validating MODIS land surface reflectance and albedo products:Methods and preliminary results［J］.Remote Sens Environ，83(1):149-162.doi:10.1016/S0034-4257(02)00092-5.

Liu J，Schaaf C，Strahler A，et al.2009.Validation of Moderate Resolution Imaging Spectroradiometer(MODIS)albedo retrieval algorithms:Dependence of albedo on solar zenith angle［J］.J Geophys Res，114，D01106.doi:10.1029/2008JD009969.

Morcrette J-L，Barker H W，Cole J N S，et al.2008.Impact of a New Radiation Package，McRad，in the ECMWF integrated forecasting system［J］.MonWeaRev，136:4773-4798.doi:10.1175/2008MWR2363.1.

Myhre G，Kvalevag M M，Schaaf C B.2005.Radiative forcing due to anthropogenic vegetation change based on MODIS surface albedo dataset［J］.GeophysResLett，32，L21410.doi:10.1029/2005GL024004.

Pinty B，Roveda F，Verstraete M M，et al.2000a.Surface albedo retrieval from Meteosat，1，Theory［J］.J Geophys Res，105(D14):18099-18112.

Pinty B，Roveda F，Verstraete M M，et al.2000b.Surface albedo retrieval from Meteosat，2，Applications［J］.J Geophys Res，105(D14):18113-18134.

Roesch A，Schaaf C B，Gao F.2004.Use of Moderate-Resolution Imaging Spectroradiometer bidirectional reflectance distribution function products to enhance simulated surface albedo［J］.J Geophys Res，109，D12105.doi:10.1029/2004JD004552.

Roman M O，Schaaf C B，Woodcock C E，et al.2009.The MODIS(Collection V005)BRDF/albedo product:Assessment of spatial representativeness over forested landscapes［J］.Remote Sens Environ，113:2476-2498.

Salomon J G，Schaaf C B，Strahler A H，et al.2006.Validation of the MODIS bidirectional reflectance distribution function and albedo retrievals using combined observations from the Aqua and Terra platforms［J］.IEEE Trans Geosci Remote Sens，44:1555-1565.doi:10.1109/TGRS.2006.871564.

Samain O，Kergoat L，Hiernaux P，et al.2008.Analysis of the in-situ and MODIS albedo variability at multiple time scales in the Sahel［J］.J Geophys Res，113，D14119.doi:10.1029/2007JD009174.

Schaaf，C B，Gao F，Strahler A H，et al.2002.First operational BRDF，albedo and nadir reflectance products from MODIS［J］.Remote Sens Environ，83(1):135-148.

Schaaf C B，Liu Jicheng，Gao F，et al.2011.MODIS albedo and reflectance anisotropy products from Aqua and Terra［C］//Ramachandran B，Justice C，Abrams M.Land remote sensing and global environmental change:NASA’s earth observing system and the science of ASTER and MODIS，remote sensing and digital image processing series，11，Springer-Cerlag:873.

Sellers P J，Meeson B W，Hall F G，et al.1995.Remote sensing of the land surface for studies of global change:Models—algorithms—experiments［J］.Remote Sens Environ，51(1):3-26.

Sellers P J，Randall D A，Collatz G J，et al.1996.A revised land surface parameterization(SiB2)for atmospheric GCMs，part I，Model formulation［J］.J Climate，9(4):676-705.

Stroeve J，Box J E，Gao Feng，et al.2005.Accuracy assessment of the MODIS 16-day albedo product for snow:Comparisons with Greenland in situ measurements［J］.Remote Sens Environ，94(1):46-60.

Susaki J，Yasuoka Y，Kajiwara K，et al.2007.Validation of MODIS albedo products of paddy fields in Japan［J］.IEEE Trans Geosci Remote Sens，45(1):206-217.doi:10.1109/TGRS.2006.882266.

Tian Y，Dickinson R E，Zhou L，et al.2004.Land boundary conditions from MODIS data and consequences for the albedo of a climate model［J］.GeophysResLett，31，L05504.doi:10.1029/2003GL019104.

Wagner S C，Govaerts Y M，Lattanzio A.2010.Joint retrieval of surface reflectance and aerosol optical depth from MSG/SEVIRI observations with an optimal estimation approach:2.Implementation and evaluation［J］.JGeophysRes，115，D02204.doi:10.1029/2009JD011780.

Wang K，Liu J，Zhou X，et al.2004.Validation of the MODIS global land surface albedo product using ground measurements in a semidesert region on the Tibetan Plateau［J］.J Geophys Res，109，D05107.doi:10.1029/2003JD004229.

Wydick J E，Davis P A，Gruber A.1987.Estimation of broadband planetary albedo from operational narrowband satellite measurements［R］//NOAA Technical Report NESDIS 27.US Department of Commerce，National Oceanic and Atmospheric Administration and National Environmental Satellite，Data＆Information Administration.

Zhou L，Dickinson R E，Tian Y，et al.2003.Comparison of seasonal and spatial variations of albedos from Moderate-Resolution Imaging Spectroradiometer(MODIS)and Common Land Model［J］.J Geophys Res，108(D15)，4488.doi:10.1029/2002JD003326.

Assessment on the accuracy of MODIS albedos over the Tibetan Plateau

CHEN Ai-jun1，2，LIANG Xue-wei1，2，BIAN Lin-gen3，LIU Yu-jie4，WANG Fei1，2，ZHU Xiao-xiang5

(1.Key Laboratory of Meteorological Disaster of Ministry of Education，NUIST，Nanjing 210044，China;2.School of Atmospheric Physics，NUIST，Nanjing 210044，China;3.Chinese Academy of Meteorological Sciences，Beijing 100081，China;4.National Satellite Meteorological Center，Beijing 100081，China;5.National Meteorological Center，Beijing 100081，China)

Using ground-observed albedos from four stations over the Tibetan Plateau from 2002 to 2004 in the CAMP/Tibet experiments，the accuracies of Terra MODIS(MODerate-resolution Imaging Spectroradiometer)shortwave bands(0.3—5.0 μm)albedos in 1 km resolution derived with the full inversion and the magnitude inversion are individually quantitatively assessed.When MODIS albedos are derived with the full inversion，the root-mean-square errors(RMSEs)between the ground-observed albedos and the black-sky albedos(BSA)，and the white-sky albedos(WSA)are 0.018 7 and 0.016 8，respectively;when MODIS albedos are derived with the magnitude inversion，they are 0.076 6 and 0.076 1，respectively.And，when MODIS albedos derived with the full inversion and the magnitude inversion are all included，the RMSEs are 0.067 9 and 0.067 5，respectively.When the ground-observed albedos and MODIS albedos derived with the magnitude inversion are both from the no-snow condition，RMSEs between the ground-observed albedos and BSA，and WSA derived with the magnitude inversion are 0.035 2 and 0.036 4，respectively.However，when the ground-observed albedos are fromthe snowy condition，while MODIS albedos derived with the magnitude inversion are from the no-snow condition，RMSEs between the ground-observed albedos and BSA，and WSA derived with the magnitude inversion are 0.155 6 and 0.154 1.

albedo;validation;MODIS;Tibetan Plateau

P407

A

1674-7097(2012)06-0664-09

2012-05-06;改回日期:2012-09-23

国家自然科学基金资助项目(40875015);江苏高校优势学科建设工程资助项目(PAPD)

陈爱军(1972—)，男，湖北荆门人，博士，副教授，研究方向为气象卫星资料分析与应用，chenaijun@nuist.edu.cn.

陈爱军，梁学伟，卞林根，等.2012.青藏高原地区MODIS反照率的精度分析［J］.大气科学学报，35(6):664-672.

Chen Ai-jun，Liang Xue-wei，Bian Lingen，et al.2012.Assessment on the accuracy of MODIS albedos over the Tibetan Plateau［J］.Trans Atmos Sci，35(6):664-672.(in Chinese)

(责任编辑:倪东鸿)