文章快速检索     高级检索
  高原气象  2016, Vol. 35 Issue (6): 1409-1418  DOI: 10.7522/j.issn.1000-0534.2015.00111
0

引用本文 [复制中英文]

陈爱军, 梁学伟, 卞林根, 等 . 2016. 青藏高原MODIS地表反照率反演结果的空间分布[J]. 高原气象, 35(6): 1409-1418. DOI: 10.7522/j.issn.1000-0534.2015.00111
[复制中文]
CHEN Aijun, LIANG Xuewei, BIAN Lingen, et al . 2016. Spatial Distribution Characteristics of MODIS Land Surface Albedo Inversions over the Qinghai-Xizang Plateau[J]. PLATEAU METEOROLOGY, 35(6): 1409-1418. DOI: 10.7522/j.issn.1000-0534.2015.00111.
[复制英文]

资助项目

国家自然科学基金项目(40875015);安徽省气象局预报员专项项目(kY201614)

作者简介

陈爱军(1972-), 男, 湖北荆门人, 副教授, 主要从事气象卫星资料分析与应用研究.E-mail:chenaijun@nuist.edu.cn

文章历史

收稿日期: 2015-02-27
定稿日期: 2015-12-17
青藏高原MODIS地表反照率反演结果的空间分布
陈爱军1,2, 梁学伟3, 卞林根4, 刘玉洁5     
1. 南京信息工程大学 气象灾害预报预警与评估协同创新中心, 南京 210044;
2. 南京信息工程大学 大气物理学院, 南京 210044;
3. 安徽省气象台, 合肥 230031;
4. 中国气象科学研究院, 北京 100081;
5. 国家卫星气象中心, 北京 100081
摘要: 为了促进MODIS地表反照率产品的应用及其反演算法的改进,应用MODIS地表反照率反演质量数据MCD43B2,并结合数字海拔高程模型DEM(Digital Elevation Model)数据和MODIS地表覆盖类型数据MOD12Q1,统计分析了20032013年青藏高原MODIS地表反照率不同反演结果的空间分布情况,结果表明:(1)MODIS地表反照率全反演结果主要分布在中西部和北部,春季全反演结果的概率超过80%,夏季向西北收缩且概率降低,秋季高概率全反演结果的范围最大、冬季最小。(2)冬季高概率当量反演结果主要分布在青藏高原西北、东南两翼,其他季节主要分布在昌都、那曲、林芝、山南地区和拉萨市。(3)在山南地区中东部、林芝及其与昌都地区交界处构成的“入”字形区域,四季无反演结果的概率都比较高,最高可达100%。(4)全反演结果的概率随海拔的增加呈下降趋势,当量反演结果则相反,无反演结果的概率在各海拔区间都较恒定且不超过10%。(5)在三种主要地表类型中,开放灌木区和裸土稀疏植被区全反演结果的概率高于草地,约为70%;草地当量反演结果的概率最高,约为30%;三者无反演结果的概率大致相当,均不足10%。以上结果表明青藏高原地区MODIS地表反照率反演质量的空间分布具有较为明显的地域特征,并且与海拔和地表覆盖类型存在一定的联系。
关键词: 青藏高原    地表反照率    反演质量    MODIS    空间分布    
1 引言

青藏高原占中国国土面积的1/4,平均海拔超过4500 m,有“世界屋脊”和“第三极”之称,是全球地形最为复杂的高原。青藏高原对大气环流有着重要影响,通过全球动量、 能量和水分循环影响着区域和全球气候变化,对中国及亚洲地区的灾害性天气乃至全球气候变化均有重大影响(叶笃正和高由禧,1979; 章基嘉等,1988)。地表反照率是控制地表辐射能量收支、 影响地气相互作用的一个重要的地表参数(Dickinson et al,1993; Blackmon et al,2001; Dai et al,2003)。获取青藏高原地表辐射和反照率,一直是各种青藏高原大气科学观测试验的一项重要内容(《青藏高原气象科学实验文集》编辑组,1984a1984b1987; 陶诗言等,1999; Zhang et al,2012; 张人禾和徐祥德,2012; 李跃青等,2010,2012; 周长艳等,2012; 李英和胡泽勇,2007; 武荣盛和马耀明,2010; 冯超等,2010)。然而,由于青藏高原独特的地形和气候特点,只有通过卫星遥感才能比较全面、 准确、 及时地了解青藏高原地表反照率的时空分布及其变化情况(方宗义等,1996; 徐兴奎和林朝晖,2002; 陈爱军等,2009)。

2013年我国启动了“青藏高原地气耦合系统变化及其全球气候效应”重大研究计划(简称青藏高原重大研究计划),该计划的三个核心科学问题均与地气相互作用有关(http://www.nsfc.gov.cn/publish/portal0/zdyjjh/002/info24876.htm),迫切需要青藏高原时空分布连续的高精度、 高时空分辨率地表反照率用于地气相互作用研究。美国EOS(Earth Observing System)/MODIS(MODerate resolution Imaging Spectroradiometer)全球地表反照率产品(编号MCD43),是目前唯一能够准实时提供青藏高原高时空分辨率地表反照率的卫星遥感反演业务化产品(Schaaf et al,20022011),对青藏高原研究具有重要的应用价值。气候和陆面过程模式不仅对地表反照率的精度要求较高,一般绝对精度达到0.02~0.05(Henderson-Sellers and Wilson,1983; Sellers et al,1995),而且要求地表反照率具有良好的时空分辨率和时空分布连续性(http://www.wmo.int/pages/prog/gcos/Publications/gcos-154.pdf)

目前,无论是全球范围(Jin et al,2003a2003b; Stroeve et al,2005; Liu et al,2009; Román et al,200920112013; Wang et al,20122014)还是青藏高原(Wang et al,20042005; 余予等,2010; 陈爱军等,2012a),大量研究都通过多种途径和方式分析和验证了MODIS地表反照率的反演精度。然而,涉及MODIS地表反照率时空分布连续性的研究比较少,尤其是不同反演结果的时空分布情况。不久前,陈爱军等(2016a2016b)利用MODIS地表反照率反演质量数据MCD43B2,研究了20032013年青藏高原MODIS地表反照率反演效率和反演质量情况,本文在此基础上将进一步研究20032013年青藏高原MODIS地表反照率不同反演结果的空间分布情况,正确认识和了解青藏高原MODIS地表反照率不同反演结果的空间分布情况,并为改进青藏高原地表反照率的卫星遥感反演研究提供参考和依据。

2.1 资料选取

所用数据为20032013年MODIS地表反照率反演质量数据MCD43B2(Collection 005)。该数据提供Terra和Aqua双星MODIS数据联合反演地表反照率的反演质量,与MODIS地表反照率数据MCD43B3具有完全对应的时空关系,空间分辨率为1 km,时间间隔为8天(Schaaf et al,2011)。

在MCD43B2中,BRDF_Albedo_Band_Quality的科学数据集(Scientific Data Set,SDS)提供MODIS 1~7波段BRDF参数的反演质量(http://lpdaac.usgs.gov/products/modis_products_table/mc-d43b2)。该科学数据集由500 m分辨率MCD43A2数据对应的SDS转换而来。与1 km分辨率数据对应的4个500 m分辨率像元中,如果至少3个像元是“最佳全反演结果”,质量标记为0; 如果至少3个像元是“全反演结果”,质量标记为1; 如果最多2个像元为“全反演结果”,且最多有1个像元为填充值,质量标记为2; 如果4个像元均为“当量反演结果”,或者最多2个像元为填充值,质量标记为3; 如果至少3个像元为填充值,质量标记为4; 4个像元均为填充值时,质量标记为15。由于质量标记1、 2和4所占比例都较小(陈爱军等,2016a),因此,本文将MCD43B2中质量标记0和1对应的地表反照率,称为全反演结果; 质量标记2和3对应的地表反照率,称为当量反演结果; 质量标记4和填充值15对应的结果称为无反演结果。

MODIS地表反照率数据MCD43采用MODLAND正弦投影(http://modis-land.gsfc.nasa.gov/MODLAND_grid.html),正弦投影编号H24V05、 H25V05、 H25V06、 H26V05和H26V06的5个区域正好覆盖青藏高原。研究所用数据为这5个区域的MCD43B2数据,由美国地质调查局USGS(United States Geological Survey)地球资源观测与科学EROS(Earth Resource Observation and Science)数据中心提供,数据质量完全可以满足科学研究与分析要求。

2.2 方法简介

由于相同反演结果的比例在MODIS 1~7通道的差异较小(陈爱军等,2016a),本文以1~7通道相同反演结果的平均比例统计2003-2013年青藏高原MOIDS地表反照率不同反演结果的平均概率,即:

${{p}_{x}}=\frac{\underset{i=1}{\overset{11}{\mathop{\sum }}}\,{{\left( {{N}_{i}} \right)}_{x}}}{11\times n}\times 100%,$ (1)

其中: x为不同的MODIS地表反照率反演结果,即全反演结果、 当量反演结果或无反演结果;(Ni)x为反演结果x在第i年的总频数; n为统计的反演周期总数; px为在n个反演周期内获得反演结果x的11年平均概率。然后,以不同反演结果的11年平均概率(如无特殊说明,下同)的空间分布,研究青藏高原MODIS地表反照率全反演结果、 当量反演结果和无反演结果的空间分布情况。

3 结果分析 3.1 总体情况

由青藏高原MODIS地表反照率不同反演结果的概率空间分布(图 1)可知,青藏高原上地表反照率全反演结果的概率大多超过60%,一些地区超过80%,例如: 青海海西蒙古族藏族自治州、 海南藏族自治州和西藏阿里地区、 日喀则地区、 那曲地区的绝大部分区域。然而,西藏林芝地区、 山南地区东南部、 那曲地区东南部、 昌都地区西北部,以及青海玉树藏族自治州东部、 果洛藏族自治州的部分区域等少数地区,全反演结果的概率比较低,不足40%,雅鲁藏布江下游至横断山脉的林芝地区甚至不足20%(图 1a)。全反演结果概率较低的区域,当量反演结果或者无反演结果的概率自然较高。例如: 那曲、 昌都、 林芝等三地交界处,当量反演结果的概率超过60%(图 1b); 山南地区东南部、 林芝地区及其与昌都地区的交界处形成的“入”字形区域,无反演结果概率超过40%,最高甚至超过90%(图 1c)。湖泊地区全反演结果的概率较低,当量反演结果的概率较高。例如: 青海境内的青海湖、 西藏境内的色林错湖和纳木错湖,全反演结果的概率只有约30%,当量反演结果的概率约为60%。这主要是因为湖泊上空水汽密度大,容易形成云,卫星过境时无法获得充足的晴空观测数据,无法获得全反演结果。青藏高原当量反演结果和无反演结果概率较高的地区,大多是因为云覆盖而无法获得充足的晴空观测数据。

图 1 青藏高原MODIS地表反照率不同反演结果的概率分布 (a)全反演结果,(b)当量反演结果,(c)无反演结果,字母A-F分别表示青海省海北藏族自治州、 海东地区、 海南藏族自治州、 黄南藏族自治州、 果洛藏族自治州和海西蒙古族藏族自治州(格尔木市代管),网格左上角: 78.32°E,40.00°N, 右上角: 117.47°E,40.00°N; 左下角: 63.86°E,20.00°N,右下角: 95.77°E,20.00°N Figure 1 The distribution of the probability for different kinds of MODIS land surface albedo over the Qinghai-Xizang Plateau. (a)for the full inversions,(b)for the magnitude inversions,(c)for no inversions.Characters A-F in these figures are Haibei Tibetan Autonomous Prefecture,Haidong Prefecture,Hainan Tibetan Autonomous Prefecture,Huangnan Tibetan Autonomous Prefecture,Golog Tibetan Autonomous Prefecture and Haixi Mongolian-Tibetan Autonomous Prefecture(Governed by Golmud City),respectively.The upper left corner of the grid is 78.32°E,40.00°N,the upper right corner is 117.47°E, 40.00°N,the bottom left corner is 63.86°E,20.00°N,the bottom right corner is 95.77°E,20.00°N
3.2 季节特征

本文将第9~19个8天期间(为3月6日至6月1日)划为春季,第20~30个8天期间(为6月2日至8月28日)划为夏季,第31~42个8天期间(为8月29日至12月2日)划为秋季,第43~45个8天期间(为12月326日)和次年第1~8个8天期间(为1月1日至3月5日)划分为冬季,进一步分析青藏高原MODIS地表反照率不同反演结果空间分布的季节特征。

春季除东南部之外,青藏高原其他区域全反演结果的概率超过80%,甚至达到100%。在东南部,包括青海省果洛藏族自治州、 玉树藏族自治州东部以及西藏自治区昌都地区、 林芝地区和山南地区的东部等地区,全反演结果的概率较低,一般不足30%; 林芝地区和山南地区部分区域全反演结果的概率几乎为零,而且这些区域明显呈“入”字形(图 2a)。夏季,西藏全反演结果概率较大的范围整体向西北收缩,日喀则地区,以及那曲和阿里两个地区的东南部,全反演结果的概率降至50%以下,其中部分地区的概率甚至不足10%。春季林芝和山南地区全反演结果概率为零且呈“入”字形的区域,夏季全反演结果的概率仍然为零,且范围呈弥散状扩大(图 2b)。秋季,全反演结果概率较高的范围又向青藏高原东南方向扩展,其范围在四季中最大: 除林芝、 山南地区呈“入”字形区域全反演结果的概率仍然较低之外,其他区域全反演结果的概率基本都超过80%,部分地区甚至接近100%(图 2c)。冬季,全反演结果概率较高的区域,从青藏高原西北、 东南两侧向中部收缩,其范围为四季中最小,只有西藏阿里地区中东部、 那曲地区中东部、 日喀则地区、 拉萨市、 昌都地区东南部、 青海海西蒙古族藏族自治州及海北藏族自治州、 海南藏族自治州的部分区域,全反演结果的概率比较高。林芝和山南地区“入”字形区域全反演结果的概率仍然较低,约为10%(图 2d)。

图 2 青藏高原MODIS地表反照率全反演结果不同季节的概率分布 (a)~(d)分别为春、 夏、 秋、 冬季,其余说明同图 1 Figure 2 The distribution of the probability in different seasons for MODIS land surface albedo full inversions over the Qinghai-Xizang Plateau. (a)~(d)are for spring,summer,autumn and winter,respectively.Others are the same as Fig. 1

对于当量反演结果,春季在青海果洛藏族自治州南部、 玉树藏族自治州南部,以及西藏昌都地区、 那曲地区的东南部、 林芝地区、 拉萨市、 山南地区等地的概率最高超过70%; 其他地区的概率较低,约为10%,甚至为零(图 3a); 夏季,当量反演结果概率较高的区域扩大至日喀则地区,但林芝地区的概率却下降至10%左右,有些区域甚至为零(图 3b); 秋季,当量反演结果的概率基本在10%左右,只有西藏林芝地区,以及那曲、 昌都和青海玉树藏族自治州交界处,当量反演结果的概率略高,约为60%(图 3c); 冬季,除上述区域之外,当量反演结果概率较高的区域还有那曲地区的西北部、 玉树藏族自治州的西北部、 海西蒙古族藏族自治州(格尔木市代管)的西北部、 海西蒙古族藏族自治州的东北部(图 3c)。此外,一年四季湖泊水体区域当量反演结果的概率都比较高。

图 3 青藏高原MODIS地表反照率当量反演结果不同季节的概率分布 (a)~(d)分别为春、 夏、 秋、 冬季,其余说明同图 1 Figure 3 The distribution of the probability in different seasons for MODIS land surface albedo magnitude inversions over the Qinghai-Xizang Plateau. (a)~(d)are for spring,summer,autumn and winter,respectively.Others are the same as Fig. 1

无论春、 夏、 秋、 冬季,青藏高原东南部林芝地区与山南地区东南部的“入”字形区域,无反演结果的概率都较高(图 4)。春季,只有“入”字形及其附近区域无反演结果的概率较高(图 4a)。夏季,不只是“入”字形区域,该区域下方的三角形区域,无反演结果的概率进一步增大,大多超过80%,甚至达到100%(图 4b)。秋季,虽然无反演结果概率较高的区域仍然集中在“入”字形区域,但是范围明显减小(图 4c)。冬季,青藏高原东南部“入”字形区域无反演结果概率较高的区域进一步缩小,但青海海西蒙古族藏族自治州、 西藏阿里地区西北部、 那曲地区西北部,以及青海玉树藏族自治州等地都弥散性地分布着无反演结果概率较高的区域(图 4d)。

图 4 青藏高原MODIS地表反照率无反演结果不同季节的概率分布 (a)~(d)分别为春、 夏、 秋、 冬季,其余说明同图 1 Figure 4 The distribution of the probability in different seasons for MODIS land surface albedo no inversions over the Qinghai-Xizang Plateau. (a)~(d)are for spring,summer,autumn and winter,respectively.Others are the same as Fig. 1

上述分析表明,春季除东南部的部分区域之外的其他区域,都能以较高的概率获得MODIS地表反照率全反演结果; 夏季获得全反演结果概率较高的区域向西北方向收缩; 秋季获得全反演结果概率较高的范围最大;冬季获得全反演结果概率较高的范围最小,从西北、 东南两个方向收缩至青藏高原的中部。无论春夏秋冬,东南部林芝地区与山南地区东南部的“入”字形区域,无反演结果的概率始终都比较高。对于当量反演结果,除冬季分布在青藏高原西北和东南两翼之外,其他季节主要分布在西藏昌都地区、 那曲地区、 林芝地区、 拉萨市和山南地区,其概率一般在50%左右,最高超过70%。

3.3 随海拔的分布特征

应用DEM(Digital Elevation Model)数据分析结果表明,青藏高原海拔介于2.5~6 km的区域约占高原总面积的97%,其中: 4.5~5.0 km海拔的区域所占比例最大,约36%; 其次是5.0~5.5 km海拔的区域,所占比例约为23.5%(表 1)。各海拔区间MODIS地表反照率全反演结果的概率超过50%,其中: 2.5~3.0 km海拔区域的概率最高,接近80%; 5.5~6.0 km海拔区域的概率最低,略高于50%; 4~4.5 km海拔的区域的概率也较低,略高于60%。总体而言,青藏高原MODIS地表反照率全反演结果概率随海拔的增加呈下降趋势。随着海拔的升高,一方面云覆盖的概率加大,获得晴空观测数据的概率减小; 另一方面,短时积雪的概率也会增加。这两方面的因素都会影响MODIS地表反照率的反演质量,无法获得全反演结果。由于全反演结果的概率降低,当量反演结果的概率自然随海拔的升高呈加大的趋势,其中: 5.5~6.0 km海拔区域,当量反演结果的概率最高,约为35%; 4~4.5 km海拔区域的概率也较高,接近30%; 2.5~3.0 km海拔区域的概率最低,约为10%。除海拔介于5.5~6.0 km的区域外,其他海拔区间当量反演结果的概率都低于30%。无反演结果的概率在不同海拔高度区间相对稳定且低于10%,只有海拔介于5.5~6.0 km的区域例外,约为15%(图 5)。

表 1 青藏高原不同海拔区间的面积比例 Table 1 The ratio of the area in different elevation over the Qinghai-Xizang Plateau
图 5 青藏高原不同海拔区间MODIS地表反照率反演结果的概率 Figure 5 The probability of different kinds of MODIS land surface albedo in versions in different ranges of elevation over the Qinghai-Xizang Plateau

总之,青藏高原MODIS地表反照率全反演结果的概率随海拔的升高呈下降趋势,当量反演结果则相反,无反演结果的概率在各海拔区间基本保持稳定。

3.4 随地表类型的分布特征

应用MODIS地表覆盖类型数据MCD12Q1的分析结果表明,青藏高原地表覆盖类型以开放灌木、 草地和裸土稀疏植被为主,其中: 草地占32%,开放灌木占30%,裸土稀疏植被占28%。草地全反演结果的概率较低,略高于60%; 开放灌木区和裸土稀疏植被区全反演结果的概率大致相当,约为70%。青藏高原三种主要的地表覆盖类型中,草地当量反演结果的概率最高,约为30%; 开放灌木区和裸土稀疏植被区当量反演结果的概率比较接近,约为20%; 无反演结果在三种地表类型的概率都不足10%,且相差很小(图 6)。

图 6 青藏高原三种主要地表类型MODIS地表 反照率反演结果的概率 Figure 6 The probability of different kinds of MODIS land surface albedo inversions in three main land surface types over the Qinghai-Xizang Plateau
4 结论

(1) 青藏高原MODIS地表反照率全反演结果主要集中在中西部,包括西藏阿里地区、 日喀则地区、 那曲地区(东南少数区域除外),以及青海除玉树藏族自治州东南部、 果洛藏族自治州之外的区域,平均概率超过60%。青藏高原东南部,包括西藏林芝地区、 山南地区东南部、 那曲地区东南部、 昌都地区西北部,以及青海玉树藏族自治州东部、 果洛藏族自治州的部分区域,全反演结果的概率比较低。全反演结果概率较低的区域,当量反演结果或无反演结果的概率自然较高,尤其西藏山南地区中东部、 林芝地区及其与昌都地区交界处构成的“入”字形区域,无反演结果的概率高达90%以上。

(2) 春季青藏高原西北部MODIS地表反照率全反演结果的概率较高,超过80%,甚至达到100%,这些的区域范围在夏季不仅向西北收缩,而且概率有所降低; 秋季全反演结果概率较高的范围最大、 冬季最小。当量反演结果概率较高的区域在冬季主要分布在青藏高原西北、 东南两翼,其他季节主要分布在西藏昌都地区、 那曲地区、 林芝地区、 拉萨市和山南地区,其概率一般在50%左右,最高超过70%。在四季中,东南部林芝地区与山南地区东南部的“入”字形区域,无反演结果的概率始终都比较高。

(3)青藏高原MODIS地表反照率全反演结果的概率随海拔的升高呈下降趋势,但是都超过50%; 当量反演结果的平均比例随海拔的升高呈上升趋势,最高约为35%; 各海拔区间无反演结果的概率比较稳定,一般不超过10%。

(4) 三种主要地表类型中,开放灌木区和裸土稀疏植被区MODIS地表反照率全反演结果的概率大致相当,约为70%; 草地当量反演结果的概率最高,约为30%; 三者无反演结果的概率大致相当,均不足10%。

总之,青藏高原地区MODIS地表反照率反演质量的空间分布具有明显的地域特征,且与海拔和地表覆盖类型存在一定的联系。大量分析和验证结果表明,MODIS地表反照率全反演结果的精度完全可以满足气候和陆面过程模式的精度要求,当量反演结果则难以满足(Jin et al,2003a2003b; Stroeve et al,2005; Liu et al,2009; Román et al,200920112013; Wang et al,20122014; Wang et al,20042005; 余予等,2010; 陈爱军等,2012a); 云覆盖是导致MODIS地表反照率当量反演结果或者无反演结果的主要原因(陈爱军等,2012b2012c)。因此,对于青藏高原而言,无论是改进MODIS地表反照率反演算法,还是建立卫星遥感反演地表反照率反演新方法,未来的关键在于改进青藏高原中东部地表反照率的反演质量,将这些区域高概率的当量反演结果、 无反演结果转化为全反演结果,尤其东南部林芝地区、 山南地区东南部等“入”字形区域。只有这样,才能获得青藏高原卫星遥感反演的高时空分辨率、 时空分布连续的高精度地表反照率,促进青藏高原重大研究计划和地气相互作用研究。

致谢 感谢江苏省政府留学基金对本文的资助。
参考文献
Blackmon M, Boville B, Bryan F, et al. 2001. The community climate system model[J]. Bull Amer Meteor Soc, 82 (11): 2357–2376. DOI:10.1175/1520-0477(2001)082<2357:TCCSM>2.3.CO;2
Dai Y, Zeng X, Dickinson R E, et al. 2003. The common land model[J]. Bull Amer Meteor Soc, 84 (8): 1013–1023. DOI:10.1175/BAMS-84-8-1013
Dickinson R E, Kennedy P J, Henderson-Sellers A. 1993. Biosphere-Atmosphere Transfer Scheme (BATS) version 1e as coupled to the NCAR community climate model[M]. .
Henderson-Sellers A, Wilson M F. 1983. Surface albedo data for climatic modeling[J]. Rev Geophys, 21 (8): 1743–1778. DOI:10.1029/RG021i008p01743
Jin Y, Schaaf C B, Gao F, et al. 2003a. Consistency of MODIS surface bidirectional reflectance distribution function and albedo retrievals:1[J]. Algorithm performance[J]. J Geophys Res:Atmospheres (19842012), 108 (D5).
Jin Y, 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:Atmospheres (19842012), 108 (D5).
Liu J, Schaaf C, Strahler A, et al. 2009. Validation of Moderate Resolution Imaging Spectroradiometer (MODIS) albedo retrieval algorithm:Dependenceof albedo on solar zenith angle[J]. J Geophys Res:Atmospheres (19842012), 114 (D1).
Román 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 Sensing Environ, 113 (11): 2476–2498. DOI:10.1016/j.rse.2009.07.009
Román M O, Gatebe C K, Schaaf C B, et al. 2011. Variability in surface BRDF at different spatial scales (30~500 m) over a mixed agricultural landscape as retrieved from airborne and satellite spectral measurements[J]. Remote Sensing Environ, 115 (9): 2184–2203. DOI:10.1016/j.rse.2011.04.012
Román M O, Gatebe C K, Shuai Y, et al. 2013. Use of in situ and airborne multi-angle data to assess MODIS-and Landsat-based estimates of directional reflectance and albedo[J]. Geoscience and Remote Sensing, IEEE Transactions on, 51 (3): 1393–1404. DOI:10.1109/TGRS.2013.2243457
Schaaf C B, Gao F, Strahler A H, et al. 2002. First operational BRDF, albedo nadir reflectance products from MODIS[J]. RemoteSensing of Environment, 83 (1): 135–148.
Schaaf C B, Liu J, Gao F, et al. 2011. Aqua and Terra MODIS albedo and reflectance anisotropy products[M]//Land Remote Sensing and Global Environmental Change. Springer New York:549-561. http://cn.bing.com/academic/profile?id=535fe542ff2f016562d14c7688434ddf&encoded=0&v=paper_preview&mkt=zh-cn
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 Sensing Environ, 51 (1): 3–26. DOI:10.1016/0034-4257(94)00061-Q
Stroeve J, Box J E, Gao F, et al. 2005. Accuracy assessment of the MODIS 16-day albedo product for snow:comparisons with Greenland in situ measurements[J]. Remote Sensing Environ, 94 (1): 46–60. DOI:10.1016/j.rse.2004.09.001
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:Atmospheres (19842012), 109 (D5).
Wang K, Wang P, Liu J, et al. 2005. Variation of surface albedo and soil thermal parameters with soil moisture content at a semi-desert site on the western Tibetan Plateau[J]. Bound-Layer Meteor, 116 (1): 117–129. DOI:10.1007/s10546-004-7403-z
Wang Z, Schaaf C B, Chopping M J, et al. 2012. Evaluation of Moderate-resolution Imaging Spectroradiometer (MODIS) snow albedo product (MCD43A) over tundra[J]. Remote Sensing Environ, 117 : 264–280. DOI:10.1016/j.rse.2011.10.002
Wang Z, Schaaf C B, Strahler A H, et al. 2014. Evaluation of MODIS albedo product (MCD43A) over grassland, agriculture and forest surface types during dormant and snow-covered periods[J]. Remote Sensing Environ, 140 (1): 60–77.
Zhang R, Koike T, Xu X, et al. 2012. A ChinaJapan cooperative JICA atmospheric observing network over the Tibetan Plateau (JICA/Tibet Project):An overviews[J]. J Meteor Soc Japan, 90c : 1–16. DOI:10.2151/jmsj.2012-C01
陈爱军, 卞林根, 刘玉洁. 2009. 应用MODIS资料反演青藏高原地区地表反照率[J]. 南京气象学院学报 , 32 (2): 222–229. Chen Aijun, Bian Lingen, Liu Yujjie. 2009. Using MODIS Data to Retrieve Albedo over the Qinghai-Tibet Plateau[J]. J Nanjing Insti Meteor, 32 (2): 222–229.
陈爱军, 梁学伟, 卞林根, 等. 2012a. 青藏高原地区MODIS反照率的精度分析[J]. 大气科学学报 , 35 (6): 664–672. Chen Aijun, Liang Xuewei, Bian Lingen, et al. 2012a. Assessment on the accuracy of MODIS albedos over the Tibetan Plateau[J]. Trans Atmos Sci, 35 (6): 664–672.
陈爱军, 王飞, 卞林根, 等. 2012b. 青藏高原地区MODIS反照率两种反演结果差异的对比分析[J]. 高原气象 , 31 (6): 1479–1487. Chen Aijun, Wang Fei, Bian Lingen, et al. 2012b. Study on difference between two kinds of MODIS albedo over the Qinghai-Xizang Plateau[J]. Plateau Meteor, 31 (6): 1479–1487.
陈爱军, 刘玉洁, 卞林根, 等. 2012c. 中国地区MODIS反照率两种反演结果的比较[J]. 气象学报 , 70 (5): 1119–1127. Chen Aijun, Liu Yujie, Bian Lingen, et al. 2012c. Analysis of the differences the two kinds of MODIS albedos over China[J]. Acta Meteor Sinica, 70 (5): 1119–1127.
陈爱军, 梁学伟, 卞林根, 等. 2016a. 青藏高原MODIS地表反照率反演质量分析[J]. 高原气象 , 35 (2): 277–284. DOI:10.7522/j.issn.1000-0534.2015.00015 Chen Aijun, Liang Xuewei, Bian Lingen, et al. 2016a. Analysis on MODIS albedo retrieval quality over the Qinghai-Xizang Plateau[J]. Plateau Meteor, 35 (2): 277–284. DOI:10.7522/j.issn.1000-0534.2015.00015
陈爱军, 周婵, 卞林根, 等. 2016b. 藏北高原 GlobAlbedo 地表反照率的精度分析[J]. 高原气象 , 35 (4): 887–894. DOI:10.7522/j.issn.1000-0534.2015.00097 Chen Aijun, Zhou Chan, Bian Lingen, et al. 2016b. Assessment on the accuracy of GlobAlbedo over the Northern Qinghai-Xizang Plateau[J]. Plateau Meteor, 35 (4): 887–894. DOI:10.7522/j.issn.1000-0534.2015.00097
方宗义, 刘玉洁, 林曼芸. 1996. 青藏高原地表反照率计算研究[J]. 气象学报 , 54 (5): 580–589. Fang Zongyi, Liu Yujie, Lin Manyun. 1996. The research and calculation of surface albedo over Tibetan Plateau from satellite data[J]. Acta Meteor Sinica, 54 (5): 580–589.
冯超, 古松, 赵亮, 等. 2010. 青藏高原三江源区退化草地生态系统的地表反照率特征[J]. 高原气象 , 29 (1): 70–77. Feng Chao, Gu Song, Zhao Liang, et al. 2010. Albedo characteristics of degraded grassland ecosystem in the source region of three rivers in QinghaiTibetan Plateau[J]. Plateau Meteor, 29 (1): 70–77.
李跃清, 赵兴炳, 邓波. 2010. 2010年夏季西南涡加密观测科学试验[J]. 高原山地气象研究 , 30 (4): 80–84. Li Yueqing, Zhao Xingbing, Deng Bo. 2010. Intensive observation scientific experiment of the southwest vortex in the summer of 2010[J]. Plateau Mount Meteor Res, 30 (4): 80–84.
李跃清, 赵兴炳, 张利红, 等. 2012. 2011年夏季西南涡加密观测科学试验[J]. 高原山地气象研究 , 31 (4): 7–11. Li Yueqing, Zhao Xingbing, Zhang Lihong, et al. 2012. Intensive observation scientific experiment of the southwest vortex in the summer of 2011[J]. Plateau Mount Meteor Res, 31 (4): 7–11.
李英, 胡泽勇. 2007. 藏北高原地表反照率的初步研究[J]. 高原气象 , 25 (6): 1034–1041. Li Ying, Hu Zeyong. 2007. A preliminary study on land-surface albedo in Northern Tibetan Plateau[J]. Plateau Meteor, 25 (6): 1034–1041.
《青藏高原气象科学实验文集》编辑组. 1984a. 青藏高原气象科学实验文集(一)[R]. 北京:科学出版社. Editorial team of Anthologies of Tibetan Plateau meteorological experiment. 1984a. Anthologies of Tibetan Plateau meteorological experiment (1)[R]. Beijing:Science Press.
《青藏高原气象科学实验文集》编辑组. 1984b. 青藏高原气象科学实验文集(二)[R]. 北京:科学出版社. Editorial team of Anthologies of Tibetan Plateau meteorological experiment. 1984b. Anthologies of Tibetan Plateau meteorological experiment (2)[R]. Beijing:Science Press. http://www.oalib.com/references/15910637
《青藏高原气象科学实验文集》编辑组. 1987. 青藏高原气象科学实验文集(三)[R]. 北京:科学出版社. Editorial team of Anthologies of Tibetan Plateau meteorological experiment. 1987.Anthologies of Tibetan Plateau meteorological experiment (3)[R]. Beijing:Science Press.
陶诗言, 陈联寿, 徐祥德, 等. 1999. 第二次青藏高原大气科学试验理论研究进展(一)、(二)、(三)[M]. 北京: 气象出版社. Tao Shiyan, Chen Lianshou, Xu Xiangde, et al. 1999. The theoretical research progress of 2nd atmospheric science experiment on Tibetan Plateau (1-3)[M]. Beijing: Science Press .
武荣盛, 马耀明. 2010. 青藏高原不同地区辐射特征对比分析[J]. 高原气象 , 29 (2): 251–259. Wu Rongsheng, Ma Yaoming. 2010. Comparative analysis on radiation characteristics in different areas over the Tibetan Plateau[J]. Plateau Meteor, 29 (2): 251–259.
徐兴奎, 林朝晖. 2002. 青藏高原地表月平均反照率的遥感反演[J]. 高原气象 , 21 (3): 233–237. Xu Xingkui, Lin Chaohui. 2002. Remote sensing retrieval of surface monthly mean albedo in QinghaiXizang Plateau[J]. Plateau Meteor, 21 (3): 233–237.
叶笃正, 高由禧. 1979. 青藏高原气象学[M]. 北京: 科学出版社, 7-9. Ye Duzheng, Gao Youxi. 1979. Meteorology of the QinghaiXizang Plateau[M]. Beijing: Science Press , 7-9.
余予, 陈洪滨, 夏祥鳌, 等. 2010. 青藏高原纳木错站地表反照率观测与 MODIS 资料的对比分析[J]. 高原气象 , 29 (2): 260–267. Yu Yu, Chen Hongbin, Xia Xiangao, et al. 2010. Comparison of surface albedo measurement with MODIS product at Namco station of Tibetan Plateau[J]. Plateau Meteor, 29 (2): 260–267.
章基嘉, 朱抱真, 朱福康, 等. 1988. 青藏高原气象学进展[M]. 北京: 科学出版社. Zhang Jijia, Zhu Baozhen, Zhu Kangfu, et al. 1988. Advances in Meteorology of the QinghaiXizang Plateau[M]. Beijing: Science Press .
张人禾, 徐祥德. 2012. 青藏高原及东缘新一代大气综合探测系统应用平台——中日合作JICA项目[J]. 中国工程科学 , 14 (9): 102–112. Zhang Renhe, Xu Xiangde. 2012. An applying platform for the new generation of the comprehensive atmospheric observing system over the Tibetan Plateau and its eastern region-a China-Japan Cooperative JICA Project[J]. Eng Sci, 14 (9): 102–112.
周长艳, 张虹娇, 赵兴炳, 等. 2012. 近三十多年青藏高原大气科学试验观测布局综述[J]. 高原山地气象研究 , 32 (1): 81–87. Zhou Changyan, Zhang Hongjiao, Zhao Xingbing, et al. 2012. Overview of the observation network of the main Tibetan Plateau experiments of atmospheric sciences in recent decades[J]. Plateau Mount Meteor Res, 32 (1): 81–87.
Spatial Distribution Characteristics of MODIS Land Surface Albedo Inversions over the Qinghai-Xizang Plateau
CHEN Aijun1,2 , LIANG Xuewei3 , BIAN Lingen4 , LIU Yujie5     
1. Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters, Nanjing University of Information Science & Technology, Nanjing 210044, China;
2. School of Atmospheric Physics, Nanjing University of Information Science & Technology, Nanjing 210044, China;
3. Anhui Meteorological Observatory, Hefei 230031, China;
4. Chinese Academy of Meteorology Science, Beijing 100081, China;
5. National Satellite Meteorology Center, Beijing 100081, China
Abstract: In order to promote the application of MODIS (MODerate Imaging Spectra-radiometer) land surface albedo (LSA) and the improvement of its retrieval algorithm, the spatial distribution of MODIS retrievals over the Qinghai-Xizang Plateau (QXP) for different kind of retrieval quality with the MODIS LSA retrieval quality product MCD43B2 from 2003 to 2013 were statistically analyzed, as well as the digital elevation model (DEM) and the MODIS land surface cover classification product MOD12Q1. The results are as following:(1) The MODIS LSA full inversions (FI) are mainly located on the central and west regions, and the north regions of the QXP. The probability of FI over these regions is above 80% in spring. In summer, not only the area with high probability of FI shrinks to the north-west of the QXP, but also the probability decreases. The area with high probability of FI in autumn is the largest, while the smallest in winter. (2) Main areas with high probability of MODIS LSA magnitude inversions (MI) are distributed over the north-west and the south-east of the QXP in winter, while over Changdu Prefecture, Nagqu Prefecture, Nyingchi Prefecture, Shannan Prefecture and Lhasa City in the other three seasons. (3) The probability of no inversions (NI) for MODIS LSA is the highest over the Chinese character "into"-shaped area, which is formed in the central and east of Shannan Prefecture, Nyingchi Prefecture and its adjacent area of Changdu Prefecture. Moreover, the probability over this area is high over the four seasons and the highest is close to 100%. (4) The general trend for the probability of FI is decreasing as the increase of the elevation, while it is opposite for that of MI. The probability of NI is relatively stable and no more than 10% in each elevation range. (5) Among three main kinds of land cover type over the QXP, the probabilities of FI over the open shrubs and the naked (sparse vegetated) are around 70% and higher than that over the grassland. Accordingly, the probability of MI over the grassland is the highest and 30% or so. The probability of NI for these three kinds of land cover type is about same and no less than 10%. All above reveal that the spatial distribution of MODIS LSA retrieval quality over the QXP indicates obviously regional characteristics and it has a certain relationship with local altitude and surface cover type.
Key Words: The Qinghai-Xizang Plateau    Land surface albedo    Retrieval quality    MODIS    Spatial distribution