异质性地表反照率遥感产品真实性检验研究现状及挑战

地表反照率直接决定了地表能够吸收到的太阳辐射能量,是研究气候变化、能量平衡的一个关键参数。遥感为大尺度、连续获取地表反照率提供了一种有效的观测手段。但遥感数据本身的精度限制和反演模型的不确定性,使基于卫星数据反演的反照率产品存在误差,而这种误差的存在又会影响产品的进一步应用。正确的认识这种误差有助于提高产品的应用精度,深化其应用的深度和广度。真实性检验就是正确认识卫星反照率产品准确性、稳定性的重要手段,它是卫星产品从生产到应用的桥梁。虽然目前已经开展了大量的真实性检验工作,但即使是针对同一种卫星产品,真实性检验的结果往往并不一致。其根本原因在于验证中所采用的参考值能不能够准确地代表卫星像元尺度的地面真值。地面观测和卫星产品像元之间巨大的尺度差异以及广泛分布的地表异质性,使地面观测并不能直接作为像元尺度真值与卫星产品进行简单的对比。因此真实性检验过程并不是直接的,而是需要经过一系列严格、独立的过程得到像元尺度真值后与产品在一定的空间和时间范围内进行对比。针对目前真实性检验结果准确性和可信度不高等问题,本文尝试从地面实测数据、尺度转换、验证方式、评价方法及验证中存在的问题等几个方面来论述反照率产品的真实性检验现状及挑战。     

Forward a spatio-temporal trend surface for long-term ground-measured albedo upscaling over heterogeneous land surface

Upscaling ground albedo is challenged by the serious discrepancy between the heterogeneity of land surfaces and the small number of ground-based measurements. Conventional ground-based measurements cannot provide sufficient information on the characteristics of surface albedo at the scale of coarse pixels over heterogeneous land surfaces. One method of overcoming this problem is to introduce high-resolution albedo imagery as ancillary information for upscaling. However, due to the low frequency of updating of high-resolution albedo maps, upscaling time series of ground-based albedo measurements is difficult. This paper proposes a method that is based on the idea of conceptual universal scaling methodology for establishing a spatiotemporal trend surface using very few high-resolution images and time series of ground-based measurements for spatial-temporal upscaling of albedo. The construction of the spatiotemporal trend surface incorporates the spatial information provided by auxiliary remote sensing images and the temporal information provided by long time series of ground observations. This approach was illustrated by upscaling ground-based fine-scale albedo measurements to a coarse scale over the core study area in HiWATER. The results indicate that this method can characterize the spatiotemporal variations in surface albedo well, and the overall correlation coefficient was 0.702 during the study period.     

黑河遥感试验中尺度上推研究的进展与前瞻

尺度问题是遥感科学研究的一个关键科学问题,但其理论和方法的发展严重受限于稀缺的多尺度观测数据。黑河生态水文遥感试验(Hi WATER)的核心目标之一是开展多尺度观测以支持尺度转换研究。本文综述了Hi WATER中定点观测的尺度上推研究进展,内容包括:(1)尝试严格定义了空间平均、空间尺度上推、观测足迹、代表性误差、观测真值等概念;(2)介绍了Hi WATER获取的多尺度(单点—像元—区域—流域)生态水文观测数据;(3)发展了基于地统计理论的多尺度采样方法,改进了基于时间稳定性的采样方法;(4)定量评估了辐射、碳通量、土壤水分、地表温度单点观测的代表性误差,实证了异质性地表遥感产品真实性检验的不确定性主要来源于观测的时空代表性;(5)发展了定点观测的尺度上推方法,将克里格方法推广至回归克里格、面到面、不等精度观测等情形,发展了贝叶斯框架下的非线性尺度上推方法,实证了引入遥感观测作为协同信息可显著提高尺度上推的精度。总之,Hi WATER初步形成了从采样设计、多尺度观测、代表性误差的度量、尺度上推新方法到真实性检验的研究框架。     

关键陆表参数遥感产品真实性检验方法研究进展

遥感产品真实性检验是定量遥感基础研究的重要环节,是客观评价遥感产品精度、稳定性和一致性的重要手段,对于提高遥感定量化水平、推动和扩展遥感产品的深化应用具有重要的意义。本文基于当前陆表参数遥感产品真实性检验研究的相关成果,归纳总结了陆表参数遥感产品真实性检验方法,并就其中的5种真实性检验方法,即基于地面单点测量检验、基于地面多点采样检验、基于高分辨率数据的检验、交叉检验和间接检验,阐述了各方法的特点、局限性、适用条件。文章在最后探讨了陆表关键参数遥感产品真实性检验评价指标和验证中影响精度评价的关键因素。本文对于遥感产品真实性检验工作的开展和真实性检验系统研发具有重要的指导意义。     

MODIS V006和V005全球反照率产品精度对比分析

反照率作为一种非常重要的地表能量平衡、全球变化研究的参数,在众多研究领域中得到了广泛的应用,到目前为止已经有多种全球范围的反照率产品进行业务化生产和发布,针对不同反照率产品质量评价的研究也变得愈加重要。MODLAND团队在MODIS V005反照率产品反演算法的基础上通过改进16天周期内观测数据加权的方法生产出新版本的反照率产品MODIS V006。本文针对MODIS两个版本V005及V006的反照率产品,利用FLUXNET地面站点数据,比较验证两个版本反照率的总体精度以及在不同地表类型条件下的精度差异,同时通过交叉验证的方法分析二者的差异及稳定性。验证结果表明,MODIS V006反照率产品虽然在全反演高质量的数据比例上较V005有所下降,但是在同等条件下V006在提高时间分辨率的同时其精度也有所提高,在不同的地表类型条件下精度也优于V005,且在时间序列分布上具有稳定比例的高质量数据,可以满足大多数应用的精度需求。     

Linear unmixing of MODIS albedo composites to infer subpixel land cover type albedos

We use a linear unmixing approach to test how land use and forestry maps, in combination with the MODIS BRDF/albedo product, can be used to estimate land cover type albedos in boreal regions. Operational land use maps from three test areas in Finland and Canada were used to test the method. The resulting endmember albedo estimates had low standard errors of the mean and were realistic for the main land cover types. The estimated albedos were fairly consistent with albedo measurements conducted with a telescope mast and pure pixel albedos. Problems with the method are the possible errors in the land cover maps, lack of good quality winter MODIS albedo composites and the mismatch between the MODIS pixels and the true observation area. The results emphasize the role of tree species as determinant of forest albedo. Comprehensive spatial and temporal measurements of land cover albedo are usually not possible with in situ mast measurements, and the spatial resolution of MODIS albedo product is often too low to allow direct comparison of pixel albedos and land cover types in areas with heterogeneous vegetation. Hence, and since local forestry maps exist for most temperate and boreal regions, we believe that the proposed method will be useful in estimating average regional land cover type albedos as well as in tracking changes in them.     

Functional relation of land surface albedo with climatological variables: a review on remote sensing techniques and recent research developments

Surface albedo has been documented as one of the Essential Climate Variables (ECV) of the Global Climate Observing System (GCOS) that governs the Earth's Radiation Budget. The availability of surface albedo data is necessary for a comprehensive environmental modelling study. Thus, both temporal and spatial scale issues need to be rectified. This study reports about the availability of surface albedo data through in-situ and remote sensing satellite observations. In this paper, we reviewed the existing models for surface albedo derivation and various initiatives taken by related environmental agencies in order to understand the issues of climate with respect to surface albedo. This investigation evaluated the major activities on albedo-related research specifically for the retrieval methods used to derive the albedo values. Two main existing albedo measurement methods are derived through in-situ measurement and remotely sensed observations. In-situ measurement supported with number of instruments and techniques such aspyrheliometers, pyranometers and Baseline Surface Radiation Network (BSRN) and remotely sensed observations using angularly integrated Bi-directional Reflectance Distribution Function (BRDF) by both geostationary and polar orbit satellites. The investigation results reveals that the temporal and spatial scaling is the major issues when the albedo values are needed for microclimatic study, i.e. high-resolution time-series analyses and at heterogeneity and impervious surface. Thus, an improved technique of albedo retrieval at better spatial and temporal scale is required to fulfil the need for such kind of studies. Amongst many others, there are two downscaling methods that have been identified to be used in resolving the spatial scaling biased issues: Smoothing Filter-based Intensity Modulation (SFIM) and Pixel Block Intensity Modulation (PBIM). The temporal issues can be resolved using the multiple regression techniques of land surface temperature, selected air quality parameters, aerosol and daily skylight.     

Regional scale evapotranspiration estimation using satellite derived albedo and surface temperature

In the studies reteted to surface energy balance, satellite data provides important inputs for estimating regional surface albedo and evapotranspiration. The paper describes the use of satellite data in determining the surface emissivity over heterogeneous a’reas by taking Normalized Difference Vegetation Index (NDVI) as modulating parameter at pixel resolution. The estimated emissivity values have been used to find the surface temperature at the pixel scale. Landsat-TM-visible, NIR, TIR bands data and some ground meteorological data have been used in an energy balance model for estimating surface albedo and evapotranspiration. The ET values derived from the model are in good agreement with the values obtained with. ‘CENTURY MODEL’ and ground observations over the area, suggesting the possible use of this approach fot regional scale studies on evapotranspiration.     

遥感估算地表蒸散发真实性检验研究进展

地表蒸散发是连接土壤—植被—大气连续体的纽带,结合遥感技术估算地表蒸散发已成为获取区域乃至全球尺度时空连续地表蒸散发量的有效手段。由于遥感估算地表蒸散发容易受到地表空间异质性和近地层气象条件复杂性的影响,在模型机理与变量参数化方案、输入数据和时间尺度扩展等方面存在不确定性,影响了其准确度的提高和应用范围的拓展,因此需要开展真实性检验。本文综述了当前遥感估算地表蒸散发(包括植被蒸腾和土壤蒸发)真实性检验研究的相关成果,重点归纳并总结了应用于遥感估算地表蒸散发真实性检验的直接检验法和间接检法的主要原理、适用性和优缺点,在此基础上阐述了当前遥感估算地表蒸散发真实性检验研究所面临的挑战。分析表明：由于地表空间异质性的普遍存在,遥感估算地表蒸散发真实性检验研究在理论和方法方面还受到诸多挑战,今后应打破地表蒸散发遥感产品真实性检验局限在均匀地表的传统思路,发展非均匀地表遥感估算地表蒸散发真实性检验的理论框架,包括地表水热状况空间异质性的度量、非均匀地表验证场的优化布设、非均匀下垫面地表蒸散发的多尺度观测试验、卫星像元/区域尺度地表蒸散发相对真值的获取、验证过程中的不确定性分析以及遥感估算地表蒸散发的实证研究等,并构建一个多源、多尺度、多方法、多层次的真实性检验技术流程,以期把遥感估算地表蒸散发真实性检验作为突破口,提升相应遥感产品的应用水平,推动定量遥感科学的发展。     

融合Landsat ETM+和MODIS数据估算高时空分辨率地表短波反照率

提出一种通过融合高空间低时间分辨率、低空间高时间分辨率地表短波反照率,来估算高时空分辨率地表短波反照率的方法。首先,利用Landsat ETM+数据,通过窄波段到宽波段的转换得到一景或多景空间分辨率较高的ETM+蓝天空短波反照率;然后,在MODIS短波反照率产品基础上,以天空光比例因子为权重,得到空间分辨率较低的MODIS蓝天空短波反照率;最后,利用STARFM(Spatial and Temporal Adaptive Reflectance Fusion Model)模型融合ETM+短波反照率的空间变化信息和MODIS短波反照率的时间变化信息,得到高时空分辨率的地表短波反照率。针对STARFM模型在异质性区域估算精度降低的问题,通过以MODIS反照率影像各像元的端元(各地类)反照率取代MODIS像元反照率来提取时空变化等信息参与STARFM模型的融合过程,达到提高异质性区域估算精度的目的。结果显示,直接利用STARFM模型估算得到的高空间分辨率地表短波反照率处在合理的精度范围内(RMSE0.02),用改进后的STARFM模型估算得到的异质性区域短波反照率和真实ETM+短波反照率间的相关系数增大。     

Monitoring land surface albedo and vegetation dynamics using high spatial and temporal resolution synthetic time series from Landsat and the MODIS BRDF/NBAR/albedo product

Seasonal vegetation phenology can significantly alter surface albedo which in turn affects the global energy balance and the albedo warming/cooling feedbacks that impact climate change. To monitor and quantify the surface dynamics of heterogeneous landscapes, high temporal and spatial resolution synthetic time series of albedo and the enhanced vegetation index (EVI) were generated from the 500 m Moderate Resolution Imaging Spectroradiometer (MODIS) operational Collection V006 daily BRDF/NBAR/albedo products and 30 m Landsat 5 albedo and near-nadir reflectance data through the use of the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM). The traditional Landsat Albedo (Shuai et al., 2011) makes use of the MODIS BRDF/Albedo products (MCD43) by assigning appropriate BRDFs from coincident MODIS products to each Landsat image to generate a 30 m Landsat albedo product for that acquisition date. The available cloud free Landsat 5 albedos (due to clouds, generated every 16 days at best) were used in conjunction with the daily MODIS albedos to determine the appropriate 30 m albedos for the intervening daily time steps in this study. These enhanced daily 30 m spatial resolution synthetic time series were then used to track albedo and vegetation phenology dynamics over three Ameriflux tower sites (Harvard Forest in 2007, Santa Rita in 2011 and Walker Branch in 2005). These Ameriflux sites were chosen as they are all quite nearby new towers coming on line for the National Ecological Observatory Network (NEON), and thus represent locations which will be served by spatially paired albedo measures in the near future. The availability of data from the NEON towers will greatly expand the sources of tower albedometer data available for evaluation of satellite products. At these three Ameriflux tower sites the synthetic time series of broadband shortwave albedos were evaluated using the tower albedo measurements with a Root Mean Square Error (RMSE) less than 0.013 and a bias within the range of ±0.006. These synthetic time series provide much greater spatial detail than the 500 m gridded MODIS data, especially over more heterogeneous surfaces, which improves the efforts to characterize and monitor the spatial variation across species and communities. The mean of the difference between maximum and minimum synthetic time series of albedo within the MODIS pixels over a subset of satellite data of Harvard Forest (16 km by 14 km) was as high as 0.2 during the snow-covered period and reduced to around 0.1 during the snow-free period. Similarly, we have used STARFM to also couple MODIS Nadir BRDF Adjusted Reflectances (NBAR) values with Landsat 5 reflectances to generate daily synthetic times series of NBAR and thus Enhanced Vegetation Index (NBAR-EVI) at a 30 m resolution. While normally STARFM is used with directional reflectances, the use of the view angle corrected daily MODIS NBAR values will provide more consistent time series. These synthetic times series of EVI are shown to capture seasonal vegetation dynamics with finer spatial and temporal details, especially over heterogeneous land surfaces.     

面向对象的多尺度加权联合稀疏表示的高空间分辨率遥感影像分类

针对高空间分辨率遥感影像中的地物具有多尺度特性,以及各个尺度的对象特征对地物分类精度的影响具有较强的尺度效性,并结合面向对象影像分析方法和多尺度联合稀疏表示方法在高空间分辨率遥感影像分类中的各自优点,提出了一种面向对象的多尺度加权稀疏表示的高空间分辨率遥感影像分类算法。首先,采用多尺度分割算法获得多尺度分割结果并提取对象的多尺度特征;然后,根据影像对象的多尺度分割质量测度计算各尺度的对象权重,构建面向对象的多尺度加权联合稀疏表示模型;最后,采用2个国产GF-2高空间分辨率遥感数据集和1个高光谱-高空间分辨率航空遥感数据集（WashingtonD.C.数据）验证该算法的有效性。试验结果表明,与SVM、像素级稀疏表示、单尺度和多尺度对象级稀疏表示和深度学习等算法相比较,本文算法获得了较高的OA和Kappa分类精度,提高了各个尺度地物的分类精度,有效抑止了地物分类结果中的椒盐噪声现象,同时保持大尺度地物的区域性和小尺度地物的细节信息。     

A new calibration of the effective scattering albedo and soil roughness parameters in the SMOS SM retrieval algorithm

This study focuses on the calibration of the effective vegetation scattering albedo (ω) and surface soil roughness parameters (H_R, and N_Rp, p = H,V) in the Soil Moisture (SM) retrieval from L-band passive microwave observations using the L-band Microwave Emission of the Biosphere (L-MEB) model. In the current Soil Moisture and Ocean Salinity (SMOS) Level 2 (L2), v620, and Level 3 (L3), v300, SM retrieval algorithms, low vegetated areas are parameterized by ω = 0 and H_R = 0.1, whereas values of ω = 0.06 − 0.08 and H_R = 0.3 are used for forests. Several parameterizations of the vegetation and soil roughness parameters (ω, H_R and N_Rp, p = H,V) were tested in this study, treating SMOS SM retrievals as homogeneous over each pixel instead of retrieving SM over a representative fraction of the pixel, as implemented in the operational SMOS L2 and L3 algorithms. Globally-constant values of ω = 0.10, H_R = 0.4 and N_Rp = −1 (p = H,V) were found to yield SM retrievals that compared best with in situ SM data measured at many sites worldwide from the International Soil Moisture Network (ISMN). The calibration was repeated for collections of in situ sites classified in different land cover categories based on the International Geosphere-Biosphere Programme (IGBP) scheme. Depending on the IGBP land cover class, values of ω and H_R varied, respectively, in the range 0.08–0.12 and 0.1–0.5. A validation exercise based on in situ measurements confirmed that using either a global or an IGBP-based calibration, there was an improvement in the accuracy of the SM retrievals compared to the SMOS L3 SM product considering all statistical metrics (R = 0.61, bias = −0.019 m³ m⁻³, ubRMSE = 0.062 m³ m⁻³ for the IGBP-based calibration; against R = 0.54, bias = −0.034 m³ m⁻³ and ubRMSE = 0.070 m³ m⁻³ for the SMOS L3 SM product). This result is a key step in the calibration of the roughness and vegetation parameters in the operational SMOS retrieval algorithm. The approach presented here is the core of a new forthcoming SMOS optimized SM product.     

Impact of the spatial resolution on the energy balance components on an open-canopy olive orchard

The recent technical improvements in the sensors used to acquire images from land surfaces has made possible to assess the performance of the energy balance models using unprecedented spatial resolutions. Thus, the objective of this work is to evaluate the response of the different energy balance components obtained from METRIC model as a function of the input pixel size. Very high spatial resolution airborne images (≈50 cm) on three dates over olive orchards were used to aggregate different spatial resolutions, ranging from 5 m to 1 km. This study represents the first time that METRIC model has been run with such high spatial resolution imagery in heterogeneous agricultural systems, evaluating the effects caused by its aggregation into coarser pixel sizes. Net radiation and soil heat flux showed a near insensitive behavior to spatial resolution changes, reflecting that the emissivity and albedo respond linearly to pixel aggregation. However, greater discrepancies were obtained for sensible (up to 17%) and latent (up to 23%) heat fluxes at spatial resolutions coarser than 30 × 30 m due to the aggregation of non-linear components, and to the inclusion of non-agricultural areas in such aggregation. Results obtained confirm the good performance of METRIC model when used with high spatial resolution imagery, whereas they warn of some major errors in crop evapotranspiration estimation when medium or large scales are used.     

The influence of wind speed on infrared temperature in impervious surface areas based on in situ measurement data

Wind perturbations can cause a relatively rapid decay in infrared temperature, thus resulting in abnormal spatial patterns of infrared temperature in urban areas and the subsequent reduction in the reliability of infrared temperature measurements. To increase the reliability of such measurements, the effects of wind speed must be evaluated and removed. However, studies on the quantitative estimation of wind speed effects on infrared temperature are limited. In this study, in situ infrared temperature measurements and synchronous meteorological data were used to evaluate the influence of wind speed on in situ infrared temperature measurements of impervious surfaces. Five different impervious surfaces were selected in this study. The technical schemes are proposed for quantitative estimation of wind speed effects: (1) the residual-based method from the diurnal temperature cycle model was proposed to estimate the infrared temperature decay (ITD) due to wind fluctuations; (2) quantile regression method was introduced to define the relationship between wind speed fluctuations and the ITD; and (3) An improved probabilistic prediction interval as well as a ratio method were developed to estimate the magnitude and duration of the ITD. The results indicated that relative extreme wind speed (EWS) was significantly correlated with the range of ITD over 5-min intervals; the hourly decay rate and impact duration of ITD varied with changes in relative wind speed and impervious surface type; and the impact duration of ITD increased with an increase in the relative EWS and lasted more than 1.3 h for the studied impervious surfaces. The above findlings provide us a guidance for in situ measurement of infrared temperature and could be utilized for correcting thermal infrared images.     

基于森林模型参数先验知识估算高分辨率叶面积指数

目前,估算高分辨率叶面积指数LAI（Leaf Area Index）的常用方法是采用大量地面测量数据和遥感数据建立统计模型,再用统计模型估算LAI。然而,与农田地面测量实验相比,森林地面测量实验获取的观测数据更加有限,这使得基于统计模型的森林高分辨率LAI的估算精度低,难以满足应用需求。为此,本文提出一种基于森林模型参数先验知识、使用森林研究区少量的LAI地面测量数据和归一化植被指数NDVI数据估算森林高分辨率LAI的方法。首先,获取全球20个森林实验区的LAI地面测量数据和NDVI数据,建立LAI-NDVI统计模型并提取森林模型参数的先验知识。然后,以一个新的森林站点Concepción作为研究区,将该研究区的数据分为建模数据和验证数据两个部分。使用研究区有限的建模数据对森林模型参数先验知识进行本地化校正得到优化模型,优化模型用于估算森林高分辨率LAI,使用验证数据评价LAI的估算精度。同时,选取了Camerons站点、Gnangara站点、Hirsikangas站点评价本文方法的LAI估算精度。使用地面测量LAI验证基于森林模型参数先验知识估算高分辨率LAI的结果精度,经验证4个森林站点的均方根误差分别为0.6680,0.4449,0.2863,0.5755。研究结果表明：在仅有少量观测数据时,采用本方法能有效地提高森林高分辨率LAI的估算精度。因此,本方法可为森林高分辨率LAI的遥感估算提供参考。     

利用中国通量网和MODIS数据评估MISR地表反照率产品

地表反照率是反映地表能量平衡的重要参数。本文通过中国陆地生态系统通量观测研究网络的实测反照率和MODIS的地表反照率产品对MISR的短波反照率数据进行验证和分析:提取了中国通量网中的8个站点的数据和对应的MODIS、MISR的反照率产品用于验证。验证的结果显示,在多数站点,MISR短波反照率能与地面数据相吻合,大部分的反演误差都集中在0.04以内;MISR与MODIS短波反照率的吻合度更高,总体的误差为0.018,均方根误差在0.04左右。总的来说,MISR地表反照率产品具有较高的反演质量。     

CBERS-02B多光谱数据在城市不透水面 估算中的可用性研究

以厦门岛为研究区,以CBERS-02B的CCD影像为数据源,采用基于可变端元的线性光谱混合模型估算了城市不 透水面组分含量,并探讨了该方法的实现过程与优势。通过端元评估确定了研究区的４个典型端元,即高反射不透水 面、低反射不透水面、高反射土壤和植被。在此基础上,以高、低反射不透水面端元的组分含量对城市不透水面含量 进行估算。精度评价结果显示：基于可变端元的方法要优于一般带全约束法;而在混合像元分解过程中加入全色波段 （band5）有助于提高模型估算精度,使得在像元尺度的精度与采用Landsat的已有报道相近,而在土地利用单元尺度实 现了对城市不透水面的无偏估计。研究实例也表明,尽管目前CBERS-02B数据在辐射定标和地理定位等方面还有待改 进,通过采用适当的处理过程和技术手段,依然能利用该数据对城市不透水面进行有效估算。     

利用AMSR-E数据反演华北平原冬小麦单散射反照率

首先,基于冬小麦不同生育期的地面实测参数,构建了组成冬小麦冠层的、包括不同尺寸和含水量的介电散射体模拟数据库,并在此基础上建立冬小麦单散射反照率和光学厚度分别在C(6.925 GHz)和X(10.65 GHz)波段之间的依赖关系。然后,根据一阶参数化模型推导得到的微波植被指数MVIs(Microwave Vegetation Indices)的物理表达式,结合AMSR-E被动微波亮温数据,反演了华北平原地区冬小麦不同生育期的单散射反照率。与MODIS日归一化差异植被指数NDVI的对比结果显示:冬小麦单散射反照率与NDVI随时间的变化趋势大致相同,但在冬小麦的抽穗期到乳熟期,NDVI呈现饱和趋势,而单散射反照率对小麦的生长变化仍旧比较敏感,在指示冬小麦生长方面具有一定优势。