西藏林芝地区混合像元MODIS地表温度产品验证

西藏林芝地区地形复杂、土地覆盖类型多样，MODIS地表温度 (land surface temperature，LST) 产品验证面临处理混合像元的难题，为获得与像元尺度 (1 km) 相匹配的地表温度数据，该文提出采用多点同时观测结合面积加权的方法，将该方法应用于验证林芝地区2013年6月10日夜间晴空MODIS/LST产品。结果显示:单点观测对像元的代表性不足，容易低估产品精度 (10个样本均方根误差为2.2 K)，面积加权法可获得综合性更好的地面LST信息，对MODIS/LST产品的精度给出更高的评价 (30个样本均方根误差为1.40 K)。对于地表类型混杂程度高且地势较为平坦的像元，面积加权法的优势更为明显，可将卫星LST产品与地面LST之间的差异由3 K降至1 K以内。

MODIS/LST Product Validation for Mixed Pixels at Linzhi of Tibet

Southeastern part of Tibet is featured with complicated terrain and diverse land cover types. Validation of MODIS/LST product (1-km spatial resolution) in this region is faced with mixed pixel issue. Point-based LST measurements cannot represent the pixel well. To obtain ground LST measurements at pixel scale, traditional method usually depends on the high spatial resolution of thermal images, such as Aster and TM. However, these data are often unavailable due to persistent cloud cover and long repeat cycle. Therefore, a new simple method called area-weighted average (AWA) method is proposed, in which land cover map at high spatial resolution is combined with multi-site field observations to model the hypothetical observations at moderate pixel scale. The assumption of AWA method is that field observations can be shared within the same land cover. The AWA method is applied and analyzed on the case of Linzhi (with an area of 20 km2) which locates in southeastern part of Tibet. First, 5 field stations are set up on 5 typical land covers: Grassland, farmland, floodplain, forest at sunny slope, and forest at shadowy slope. The upward and downward long-wave radiations are measured simultaneously. Then the land cover map at 30 m spatial resolution is derived from TM image using maximum likelihood classification method. For every 1-km MODIS pixel, the fraction of each typical land cover is calculated, and the radiation at MODIS pixel-scale is estimated through area-weighted averaging. The broadband emissivity is calculated using linear combination of narrowband emissivity of MODIS band 31 and 32. Finally, LST at MODIS pixel-scale can be calculated based on Stefan-Boltzmann law. The AWA method is used for validating daily product of MODIS/LST from Terra and Aqua platforms on 10 June 2013 (LST at night is used because it changes slowly both in temporal and spatial domain). Results show that the RMSE of MODIS/LST is below 1.4 K (n=30) when applying the AWA method. If a point-based measurement is used to directly represent a MODIS pixel, the RMSE is more than 2.2 K (n=10), showing a tendency of over-estimation. The error of Aqua LST is slightly greater than that of Terra LST, probably due to a larger sensor view zenith angle during overpass. Superiority of the AWA method is more noticeable for pixels with high land cover heterogeneity and gentle terrain. The difference in LST between satellite and field observations can be decreased from 3 K to 1 K. However, for pixels with homogeneous land covers or with very tough terrains, the advantage of AWA method is limited. To further improve the AWA method, terrain adjustment should be taken into account when extrapolating point-based measurements to the same land cover but from another region, because the slope and aspect will influence the surface energy balance process even when the land cover stays the same. Results also indicate MODIS/LST data at nighttime in Linzhi Area are accurate, which are very meaningful considering the low density of meteorological stations in this area.