HJ-1A/B星地表反照率生成算法与检验

提出了适合环境与灾害监测预报小卫星-A、B星(简称HJ-1A/B星)CCD相机的大气订正算法,并基于不同地表特性和大气条件下的辐射传输模拟数据,建立HJ-1A/B星的窄波段向宽波段反照率转换的模型.利用多级灰阶靶标实测数据、敦煌检验场实测数据验证了大气订正算法以及转换模型的可靠性和精度,并将HJ-1A/B星影像数据计算的反照率产品与同时相的MODIS反照率产品进行对比分析.结果表明:文章提出的HJ-1A/B星CCD相机大气订正算法可有效校正大气影响;窄波段向宽波段反照率转换模型反演的反照率精度可靠;基于研究成果生成的HJ-1A/B星地表反照率与MODIS反照率产品一致性较好,满足后续遥感数据定量化模型研究的精度需要.     

基于TM影像天山北坡地表反照率反演方法的研究

地表反照率是气候模型和地表能量平衡方程中的重要参数。基于6S模型估算地表反照率,忽略了地形起伏的影响,不同波段组合的地表反照率也有待进一步研究。本文选取干旱区典型流域——天山北坡三工河流域为研究区域,以TM影像为数据源逐步进行地形校正、大气校正等,从而提取窄波段地表反照率。在此基础上,根据亮度、绿度、湿度3个特征变量的物理意义,以各波段能量权重为转换参数对窄波段地表反照率进行组合,实现研究区宽波段地表反照率的反演,得出基于不同波段的物理意义的地表反照率。     

FY-3C中分辨率成像光谱仪数据的窄波段地表反照率验证研究

地表反照率数据对地表能量平衡和全球变化研究具有重要意义。基于2014年FY-3C卫星250 m分辨率的反射率数据和角度数据,选取非洲及北美洲的4个区域作为研究区,采用RossThick-LiSparseR模型作为BRDF(Bidirectional Reflectance Distribution Function)核模型反演了地表窄波段反照率,得到250 m分辨率的4个窄波段黑空、白空反照率。将反演得到的FY-3C地表窄波反照率产品与MODIS反照率产品（MCD43A3）数据进行了交叉验证,结果表明：FY-3C窄波段反照率与对应MODIS窄波段反照率对比的均方根误差在0.01~0.04,平均偏差(MBIAS)为0.09,FY-3C窄波段反照率与对应的MODIS窄波段反照率在可见光波段、近红外波段有较好的一致性。本研究提升了国产风云极轨卫星的应用范围,可为FY-3C地表反照率业务化产品提供算法支撑。     

HJ-1/CCD地表反照率估算及其与NDVI关系分析

参照TM地表反照率反演算法,建立了适用于HJ-1/CCD传感器估算地表反照率的算法。通过应用6S辐射传输模型建立查找表,对覆盖内蒙古自治区锡林浩特市的HJ-1/CCD数据进行大气校正,并根据反照率定义,回归分析得到可见光波段地表反照率;通过与地面实测数据进行对比分析,表明该算法估算得到的地表反照率精度较高,其最大相对误差为14.32%。同时,本文将估算结果与同时期NDVI进行拟合分析,得出地表反照率与NDVI存在较高的负相关关系。     

中国不同土地利用类型分光辐射地表反照率研究

选取2000～2015年遥感反演的地表反照率数据及2000和2015年两期土地利用数据,采用经典统计学方法分析不同土地利用类型分光辐射地表反照率的特征及其年际变化趋势。为了解土地利用类型的反照率特征、认知区域气候或陆面模式中能量模块的相关分光辐射变量的物理过程提供科学依据。结果表明:不同的土地利用类型具有差异明显的地表反照率特征,同种土地利用类型地表反照率的差异甚至超过了不同类别土地利用类型之间的差异,说明了地表反照率巨大的空间异质性;大部分土地利用类型在短波总辐射及分光辐射地表反照率满足:近红外短波可见光,说明短波地表反照率的上限值更大程度上取决于近红外波段的地表反照率;研究时段内,各土地利用类型在3个波段地表反照率分别呈现出了不同的变化趋势,但是大部分土地利用类型分光辐射地表反照率的年际变化速率较小,基本保持稳定。     

波谱响应函数和波宽对地表温度反演的影响

分析了波谱响应函数和波宽对红外数据反演地表温度的影响。结果表明,波谱响应函数对温度反演的影响与传感器波段设置相关,误差随波宽变大而增加。Landsat、CBERS和环境卫星等宽波段红外数据的通道辐亮度和温度的关系,用有效波长比中心波长代替积分效果更好,但都存在一定误差,针对这些传感器数据的单通道算法须考虑波谱响应的差别。窄波段的MDIS地表温度反演,有效波长比中心波长代替积分效果更好,能满足精度要求。Terra和Aqua的MODIS波谱响应函数差别很小,可相互结合提高地表温度的反演精度。     

基于MODIS数据的黑河流域土壤热惯量反演研究

热惯量法在土壤水分反演中有着广泛的应用。以MODIS数据为基础,选用真实热惯量模型,反演得到了黑河流域的土壤热惯量,为进一步研究流域土壤水分提供可靠的方法和数据。利用地面实测数据对模型参数及反演结果进行了验证,并分析了地表昼夜最大温差、地表反照率及土壤热惯量的季节性变化规律,同时对比了真实热惯量模型与表观热惯量模型反演结果与土壤水分的相关性。结果表明：地表温差、地表反照率及土壤热惯量都具有明显的季节性变化特征;真实热惯量模型相对于表观热惯量模型更有利于土壤水分的反演,且具有广泛的适用性。     

基于宽波段和窄波段植被指数的草地LAI反演对比研究

叶面积指数是一个重要的植被生理生态参数,为探讨不同植被指数反演叶面积指数的可行性,基于同空间分辨率不同光谱分辨率的HJ\|1B CCD1和Hyperion遥感影像数据,以内蒙古自治区赤峰市克斯克腾旗贡格尔草原为研究对象,选取几种常见宽波段植被指数和高光谱窄波段植被指数并结合4种常用回归模型,比较分析了不同植被指数反演叶面积指数的精度。结果表明:对于全部植被指数而言,PVI、MSAVI等综合考虑了土壤、环境等因素的植被指数较传统植被指数NDVI、RVI反演草地LAI精度更高。通过对比发现,在反演草地LAI方面,窄波段植被指数比宽波段植被指数表现出明显的优势。其中,窄波段垂直植被指数PVI验证模型的确定性系数R²为0.65,均方根误差RMSE为0.15,说明实测LAI和模拟LAI值之间具有较好的变化一致性。最后基于Hyperion影像和窄波段垂直植被指数PVI的估算模型生成研究区叶面积指数空间分布图。     

关于地表反照率遥感反演的几个问题

分析了地表反照率对陆面辐射能收支以及区域和全球气候的影响,强调了地表反照率是遥感反演陆面参数时的第一重要参数,地表反照率或多波段遥感中不同谱段的地表反射率的准确反演常常是准确估算其它陆面参数如植被和土地利用/土地覆盖等状况的先决条件。在对当前关于反照率的概念及容易混淆的术语进行阐述和说明的基础上,简述了遥感反演地表反照率的步骤和主要难点的解决方法,进而对常用陆面过程模式计算地表反照率的过程作了分析,并将其结果与MODIS有关产品进行了比较,强调了遥感与陆面过程模式和气候模式的结合。     

一个从MODIS数据同时反演地表温度和发射率的神经网络算法

MODIS的三个热红外波段29、31、32建立了三个辐射传输方程,这三个方程包含了5个未知数(大气平均作用温度、地表温度和三个波段的发射率)。用JPL提供的大约160种地物的波谱数据对MODIS三个波段(29/31/32)发射率之间的关系和用MODTRAN4对大气透过率和大气水汽含量之间关系进行模拟分析。分析结果表明地球物理参数之间存在着大量的潜在信息。由于潜在的信息难以严格地用数学表达式来描述,因此神经网络是非常适合被用来解这种病态反演问题。利用辐射传输模型(RM)和神经网络(NN)反演分析表明神经网络能够被用来精确地同时从MODIS数据中反演地表温度和发射率。地表温度的平均反演误差在0.4°C以下;波段29/31/32发射率平均反演误差都在0.008以下。     

Effect of surface properties on the narrow to broadband spectral relationship in clear sky satellite observations

Several computational experiments have been conducted to estimate the difference between clear sky spectral narrowband (0.5–0.7 μm) and broadband (0.3–2.5 μm) planetary albedo for three cases of wavelength-independent surface albedo and four cases of surface wavelength-dependent (snow, dry sand, meadow, water) albedo. The spectral interval of (0.5–0.7 μm) was selected to approximate the bulk of the VISSR visible channel on the COES satellites and Channel 1 of the AVHRR on the NOAA operational satellites. Different atmospheric conditions and solar zenith angles have been simulated. It was demonstrated that the relationship between the spectral narrowband and broadband planetary albedo depends primarily on the assumptions made about the magnitude and wavelength dependence of the surface albedo and less on the atmospheric conditions. Future attempts to parametrize the conversion from narrowband to broadband spectral observations should account for the surface type.     

Narrowband-to-broadband albedo conversion for glacier ice and snow: equations based on modeling and ranges of validity of the equations

In this paper, we propose equations for narrowband-to-broadband (NTB) albedo conversion for glacier ice and snow for four types of satellite sensors: thematic mapper (TM), advanced very high resolution radiometer (AVHRR), moderate resolution imaging spectroradiometer (MODIS), and multi-angle imaging spectroradiometer (MISR). We do this on the basis of spectral albedos and incident spectral irradiances generated with radiative-transfer models of the (sub-)surface (a two-stream model) and the atmosphere, respectively. First, we establish equations for reference values of atmospheric components and the surface elevation. These equations describe measurements with root-mean-square differences of ∼0.016. We then show that the “reference equations” also perform well when total ozone and aerosol optical depth are changed with respect to the reference. The negative effect of humidity and elevation variations on the performance of the equations can be eliminated by adding a correction term. We argue that narrowband albedos are much less sensitive to variations in the incident spectral irradiance than broadband albedos. Hence, our conclusions about the effects of variations in atmospheric composition and elevation are also valid for equations for NTB conversion proposed in other papers.     

Narrowband to broadband conversions of land surface albedo: II. Validation

In the first paper of this series, we developed narrowband to broadband albedo conversion formulae for a series of sensors. These formulae were determined based on extensive radiative transfer simulations under different surface and atmospheric conditions. However, it is important to validate the simulation results using independent measurement data. In this paper, the validation results for three broadband albedos (total-shortwave, -visible and -near-IR albedos) using ground measurement of several cover types on five different days at Beltsville, MD are presented. Results show that the conversion formulae in the previous paper are very accurate and the average residual standard errors of the resulting broadband albedos for most sensors are around 0.02, which meets the required accuracy for land surface modeling.     

Validation of MISR land surface broadband albedo

Land surface broadband albedo is a critical variable for many scientific applications. Due to the scarcity of spectral albedo measurements of the Earth's surface environments, it is useful to construct broadband albedo from spectral albedo data obtained by multi‐angle satellite observations. The Multi‐angle Imaging SpectroRadiometer (MISR) onboard NASA's Earth Observing System (EOS) Terra satellite provides land surface albedo products from multi‐angular observations; however, the products have not been comprehensively validated. We convert MISR spectral albedos to total shortwave albedos and validate them using ground measurements at different validation sites. For most surface types, a published narrowband to broadband conversion formula was used, but a new conversion formula for snow and ice covered sites is developed in this study where the spectral range of the instrument is different. Several comparisons are made: (1) between MISR directional‐hemispherical reflectance (DHR) or albedo and MODIS (Moderate Resolution Imaging Spectroradiometer) DHR; and (2) between MISR spectral DHR and bi‐hemispherical reflectance (BHR). The results show that: (1) both the value and the temporal trends of the MISR shortwave albedo and the ground measured shortwave albedo are in good agreement, with the exception of the snow and ice sites; (2) the MISR DHR conforms well to MODIS DHR; and (3) the values of MISR DHR and BHR are nearly identical.     

On the Earth's surface energy exchange determination from ERS satellite ATSR data: Part 2. Short-wave radiation

This is the second in a series of papers which discusses determination of the Earth's surface energy exchange from ERS satellite Along-Track Scanning Radiometer (ATSR) data. The paper concentrates on short-wave radiation on sea and land surfaces. In this paper, three methods were used to determine solar irradiance by using ERS ATSR-2 data. We referred to them as 'D scheme', 'T scheme' and 'O scheme'. Intercomparisons of the three schemes were carried out. The schemes were applied to the land and sea areas. The visible and near-infrared reflectances were derived from ERS-2 ATSR-2 spectral bands by using the atmospheric radiative transfer model developed by Xue and Cracknell. The narrowband reflectances are combined into a measure of surface albedo by use of a weighted averaged scheme. The schemes were applied to the land and sea areas in UK and deforestation area in Brazil. The D scheme can give solar spectral irradiance but can not give broadband solar irradiance because of the wavelength limit of sensor visible bands. The T scheme and O scheme can give good broadband solar irradiance but can not give solar spectral irradiance. The O scheme was developed by Oberhuber, which was used to create climatological datasets for GCMs. The O scheme also includes the effects of humidity and surface temperature. The O scheme is better used for daily or monthly averaged solar radiation. The other two schemes can also be developed to determine the hourly or daily solar irradiance. The results show that it is now possible to derive longterm surface solar irradiance from ATSR-2 data which can be useful in climate and hydrological studies. However, our current analysis is restricted to a small range of conditions and needs to be extended to a larger dataset.     

Narrowband to broadband conversion of Landsat TM glacier albedos

In this paper we present an empirical relationship between the broadband glacier albedo (alpha) and the narrowband glacier albedos in Landsat TM bands 2 and 4 (alpha2 and alpha4, respectively). The relationship was established on the basis of multiple linear regression analysis of 112 ground-based simultaneous measurements of alpha, alpha2 and alpha4 made at 32 sites on the tongue of the Morteratschgletscher, Switzerland. The measurements were carried out over a representative set of glacier surface types ranging from completely debris-covered glacier ice (alpha=0.08)to dry snow (alpha=0.86). The regression model explains more than 99% of the variance of the broadband albedo and the root-mean-square value of the residuals is only 0.009. The relationship enables users of Landsat TM data to make an accurate estimate of the broadband albedo on the basis of narrowband albedos without having to classify the glacier surface.     

Albedo of vegetated surfaces: its variability with differing irradiances

The albedo of four vegetated surfaces was investigated to derive its variability with differing distributions of the irradiance. The results are based on measured values of the spectral biconical reflectance factor, which are combined with calculated spectral irradiances for low and high atmospheric turbidity. The solar zenith angle is varied from 0° to 80°. The derived spectral albedos are then integrated with respect to wavelength in order to achieve the albedo. It is found that the variability of the albedo with respect to the atmospheric turbidity is less than 0.01 in nearly all cases. The variability of the albedo with respect to the solar elevation angle, however, is larger than 0.02 in many cases. For solar elevation angles from 20° to 60°, the variability of the albedo of the four surfaces can be represented by a mean curve which fits the individual variabilities with an accuracy of 0.015.     

Arctic tundra albedo and its estimation from spectral hemispheric reflectance

We present the results of a field experiment in which the nearly complete bidirectional reflectance distribution function of Alaskan arctic tundra sites early in the growing season is measured by the PARABOLA instrument. The spectral hemispheric reflectances were computed by angular integration of these measurements for three wavebands: red (650-670nm), near-infrared (810-840nm) and shortwave infrared (1620-1690 nm). Total albedo was then estimated by weighting the spectral hemispheric reflectances by the fraction of total solar irradiance in three broadband spectral regions (300-700, 700-1300 and 1300-4000nm) and representing each spectral region by the narrowband PARABOLA measurements. These calculations resulted in albedo estimates with a mean relative error of 15.7 per cent as compared to pyranometer measured albedo. Since vegetation reflectance varies significantly over each of the three broadband regions, additional reflectance weighting factors were computed from a combination of high spectral resolution canopy reflectance data and corresponding computed spectral solar irradiance. This additional reflectance weighting resulted in a reduction in the mean relative error to 7.5 per cent relative to pyranometer measured albedo. It is noted that the three spectral bands of the PARABOLA instrument data reported here are similar to those of the spectral wavebands planned for future Advanced Very High Resolution Radiometer (AVHRR) sensors on National Oceanic and Atmospheric Administration (NOAA) satellites. Therefore the results and techniques presented here may be useful for future global albedo estimation utilizing AVHRR sensors. The analysis presented here may also be applied to albedo estimation from satellite sensors with higher spectral resolution and more complete spectral coverage, such as the future orbiting MODIS sensor, in which the errors of spectral reflectance weighting will be reduced considerably due to a more complete sampling of the reflected spectrum.     

An O2-A Band Atmospheric Correction Algorithm for Tower-based Platform based on Look-up Table

Tower-based spectral observation is an important connecting bridge between flux sites and satellite remote sensing data，and the effect of atmospheric absorption and scattering between horizontal surface and tower-based platform on the atmospheric absorption band such as O2-A is difficult to ignore.Firstly，the influence of atmospheric radiation transfer on the up-welling radiance and down-welling irradiance of the tower-based platform is analyzed，and the atmospheric correction method of based on upward and downward transmittance is established，that is，the influence of the upwelling radiance and down-welling irradiance is corrected by the direct transmittance and the total transmittance.Secondly，using the simulation data of MODTRAN model，the influence of AOD550 and radiative transfer path length on atmospheric transmittance is quantitatively analyzed，and the LUT of AOD550 is established based on the ratio of down-welling irradiance of near-infrared and red bands and solar zenith angle，as well as the upward and downward atmospheric transmittance LUT based on the AOD550 and the radiative transfer path length.Finally，using the canopy spectral data of different growth stages observed by the tower-based platform，the difference of the apparent reflectance between the inside and outside of the O₂-A band absorption line before and after atmospheric correction was analyzed.The results show that the atmospheric correction method based on LUT of AOD550 and radiative transfer path length proposed in this paper can better correct the influence of upwelling radiance and down-welling on the O₂-A absorption band of the tower-based platform，and provides important method support for applications such as SIF observation on the tower platform.     

An algorithm for the retrieval of 30-m snow-free albedo from Landsat surface reflectance and MODIS BRDF

We present a new methodology to generate 30-m resolution land surface albedo using Landsat surface reflectance and anisotropy information from concurrent MODIS 500-m observations. Albedo information at fine spatial resolution is particularly useful for quantifying climate impacts associated with land use change and ecosystem disturbance. The derived white-sky and black-sky spectral albedos may be used to estimate actual spectral albedos by taking into account the proportion of direct and diffuse solar radiation arriving at the ground. A further spectral-to-broadband conversion based on extensive radiative transfer simulations is applied to produce the broadband albedos at visible, near infrared, and shortwave regimes. The accuracy of this approach has been evaluated using 270 Landsat scenes covering six field stations supported by the SURFace RADiation Budget Network (SURFRAD) and Atmospheric Radiation Measurement Southern Great Plains (ARM/SGP) network. Comparison with field measurements shows that Landsat 30-m snow-free shortwave albedos from all seasons generally achieve an absolute accuracy of ±0.02-0.05 for these validation sites during available clear days in 2003-2005, with a root mean square error less than 0.03 and a bias less than 0.02. This level of accuracy has been regarded as sufficient for driving global and regional climate models. The Landsat-based retrievals have also been compared to the operational 16-day MODIS albedo produced every 8-days from MODIS on Terra and Aqua (MCD43A). The Landsat albedo provides more detailed landscape texture, and achieves better agreement (correlation and dynamic range) with in-situ data at the validation stations, particularly when the stations include a heterogeneous mix of surface covers.