Parametrized BRDF for atmospheric and topographic correction and albedo estimation in Jiangxi rugged terrain,China

A method is presented for bi‐directional reflectance distribution function (BRDF) parametrization for topographic correction and surface reflectance estimation from Landsat Thematic Mapper (TM) over rugged terrain. Following this reflectance, albedo is calculated accurately. BRDF is parametrized using a land‐cover map and Landsat TM to build a BRDF factor to remove the variation of relative solar incident angle and relative sensor viewing angle per pixel. Based on the BRDF factor and radiative transfer model, solar direct radiance correction, sky diffuse radiance and adjacent terrain reflected radiance correction were introduced into the atmospheric‐topographic correction method. Solar direct radiance, sky diffuse radiance and adjacent terrain reflected radiance, as well as atmospheric transmittance and path radiance, are analysed in detail and calculated per pixel using a look‐up table (LUT) with a digital elevation model (DEM). The method is applied to Landsat TM imagery that covers a rugged area in Jiangxi province, China. Results show that atmospheric and topographic correction based on BRDF gives better surface reflectance compared with sole atmospheric correction and two other useful atmospheric‐topographic correction methods. Finally, surface albedo is calculated based on this topography‐corrected reflectance and shows a reasonable accuracy in albedo estimation.     

Variability in surface BRDF at different spatial scales (30 m-500 m) over a mixed agricultural landscape as retrieved from airborne and satellite spectral measurements

Over the past decade, the role of multiangle remote sensing has been central to the development of algorithms for the retrieval of global land surface properties including models of the bidirectional reflectance distribution function (BRDF), albedo, land cover/dynamics, burned area extent, as well as other key surface biophysical quantities impacted by the anisotropic reflectance characteristics of vegetation. In this study, a new retrieval strategy for fine-to-moderate resolution multiangle observations was developed, based on the operational sequence used to retrieve the Moderate Resolution Imaging Spectroradiometer (MODIS) Collection 5 reflectance and BRDF/albedo products. The algorithm makes use of a semiempirical kernel-driven bidirectional reflectance model to provide estimates of intrinsic albedo (i.e., directional-hemispherical reflectance and bihemispherical reflectance), model parameters describing the BRDF, and extensive quality assurance information. The new retrieval strategy was applied to NASA's Cloud Absorption Radiometer (CAR) data acquired during the 2007 Cloud and Land Surface Interaction Campaign (CLASIC) over the well-instrumented Atmospheric Radiation Measurement Program (ARM) Southern Great Plains (SGP) Cloud and Radiation Testbed (CART) site in Oklahoma, USA. For the case analyzed, we obtained ~ 1.6 million individual surface bidirectional reflectance factor (BRF) retrievals, from nadir to 75° off-nadir, and at spatial resolutions ranging from 3 m to 500 m. This unique dataset was used to examine the interaction of the spatial and angular characteristics of a mixed agricultural landscape; and provided the basis for detailed assessments of: (1) the use of a land cover type-specific a priori knowledge in kernel-driven BRDF model inversions; (2) the interaction between surface reflectance anisotropy and instrument spatial resolution; and (3) the uncertainties that arise when sub-pixel differences in the BRDF are aggregated to a moderate resolution satellite pixel. Results offer empirical evidence concerning the influence of scale and spatial heterogeneity in kernel-driven BRDF models; providing potential new insights into the behavior and characteristics of different surface radiative properties related to land/use cover change and vegetation structure.     

Validation of Narrow-band Surface Albedo Retrieved from FY-3C MERSI Satellite Data

Surface albedo is one of the driving factors in surface radiant energy balance and surface-atmosphere interaction.It is widely used in surface energy balance, medium and long-term weather forecasting and global change research.This study aims to validate the surface albedo retrieved from FY-3C MERSI. This paper selected four regions in Africa and North America as study areas to validate the retrieved albedo from the reflectance data and angle data of FY-3C MERSI at 250 m resolution in 2014. The semi-empirical kernel-driven BRDF(bidirectional reflectance distribution function) model RossThick-LiSparseR and least squares fitting method were used to calculate the parameter of BRDF. Then four narrow-band black-sky albedos and four narrow-band white-sky albedos can be obtained by angle integration. Finally, the cross-validation of FY-3C surface narrow-band albedo products with MODIS albedo products (MCD43A3) was carried out. The results show that theRoot Mean Square Error(RMSE) between the FY-3C narrow-band albedo and the corresponding MODIS narrow-band albedo is in the range of 0.01 ~ 0.04, and the Mean Bias (MBIAS) is 0.09. FY-3C narrow-band albedo has good consistency with the corresponding MODIS narrow-band albedo in the visible and near-infrared bands. So, the methodologyof using the BRDF model to invert the surface albedo of FY-3C medium resolution imaging spectrometer data is feasible and reliable. The further improvement of the inversion accuracy of FY3C-MERSI surface albedo also depends on the improvement of basic data processing quality, including image geometric correction, calibration, and strict data quality control.     

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地表反照率业务化产品提供算法支撑。     

Regularized kernel-based BRDF model inversion method for ill-posed land surface parameter retrieval

In this paper, we consider the direct solution of the kernel-based bidirectional reflectance distribution function (BRDF) models for the retrieval of land surface albedos. This is an ill-posed problem due to nonuniqueness of the solution and the instability induced by error/noise and small singular values of the linearized system or the linear BRDF model. A robust inversion algorithm is critical for the BRDF/albedo retrieval from the limited number of satellite observations. We propose a promising algorithm for resolving this kind of ill-posed problem encountered in BRDF model inversion using remote sensing data.New techniques for robust estimation of BRDF model parameters are needed to cope with the scarcity of the number of observations. We are reminded by Cornelius Lanczos' dictum: “Lack of information cannot be remedied by mathematical trickery.” Thus identifying a priori information or appropriate constraints, and the embedding of the information or constraints into the regularization algorithm, are pivotal elements of a retrieval algorithm. We develop a regularization method, which is called the numerically truncated singular value decomposition (NTSVD). The method is based on the spectrum of the linear driven kernel, and the a priori information/constraint is based on the minimization of the l₂ norm of the parameters vector. The regularization algorithm is tested using field data as well as satellite data. Numerical experiments with a subset of measurements for each site demonstrate the robustness of the algorithm.     

Assessment of biases in MODIS surface reflectance due to Lambertian approximation

Using MODIS data and the AERONET-based Surface Reflectance Validation Network (ASRVN), this work studies errors of MODIS atmospheric correction caused by the Lambertian approximation. On one hand, this approximation greatly simplifies the radiative transfer model, reduces the size of the look-up tables, and makes operational algorithm faster. On the other hand, uncompensated atmospheric scattering caused by Lambertian model systematically biases the results. For example, for a typical bowl-shaped bidirectional reflectance distribution function (BRDF), the derived reflectance is underestimated at high solar or view zenith angles, where BRDF is high, and is overestimated at low zenith angles where BRDF is low. The magnitude of biases grows with the amount of scattering in the atmosphere, i.e., at shorter wavelengths and at higher aerosol concentration. The slope of regression of Lambertian surface reflectance vs. ASRVN bidirectional reflectance factor (BRF) is about 0.85 in the red and 0.6 in the green bands. This error propagates into the MODIS BRDF/albedo algorithm, slightly reducing the magnitude of overall reflectance and anisotropy of BRDF. This results in a small negative bias of spectral surface albedo. An assessment for the GSFC (Greenbelt, USA) validation site shows the albedo reduction by 0.004 in the near infrared, 0.005 in the red, and 0.008 in the green MODIS bands.     

Land surface albedo retrieval via kernel-based BRDF modeling: I. Statistical inversion method and model comparison

The land surface albedo is a key parameter influencing the climate near the ground. Therefore, it must be determined with sufficient accuracy. In this paper, a statistical inversion method is presented in support of the application of kernel-based Bi-directional Reflectance Distribution Function (BRDF) models for the calculation of the surface albedo. The method provides the best linear unbiased estimations (BLUE) of the BRDF model coefficients for an arbitrary number of available angular measurements. When the number of measurements exceeds the number of the estimated coefficients, the QR decomposition method is proposed to improve the ill-conditional features of the inversion matrix. In other cases, the singular value decomposition (SVD) method is suggested. The proposed inversion method is innovative in that it provides confidence intervals for each of the BRDF model coefficients with a prescribed significance expressed by a probability level. Five candidate kernel-driven BRDF models were used in the present simulation study: Li-Sparse, Roujean, Li-Sparse-Wanner, Li-Dense and Walthall. A ground-based reflectance measurement data set including 11 surface types forms the background for the inversion experiments. The results show a strong dependence on the solar zenith angle (SZA) and on the land cover type (LCT) for all candidate models. Owing to this, none model could be recommend in a general manner. The Li-Sparse and the Li-Sparse-Wanner models performed the best for the grass and wheat LCT, while the Roujean model appeared as a favorite for the pine and deciduous forests. The implementation of the confidence interval technique shows that the BRDF model coefficients can be retrieved with an uncertainty of 20-30%, and somewhat greater in the case of forest. The measured angular reflectance curves lie, as a rule, within the uncertainty bands related to the 5% significance level (95% probability). The corresponding albedo estimates can be characterized by an absolute uncertainty of 1-2% in the visible band and 5-10% in the near infrared band, or by 10-30% in relative terms. The reflectance measurements at low SZA values are preferable for BRDF model inversion for the grassland and crop, while medium range of SZA seems to provide more information on forest features. For the majority of LCT, the results of BRDF model inversion seem to be less reliable when considering multi-angular measurements for various SZA than for a single SZA.     

Remote sensing of woody shrub cover in desert grasslands using MISR with a geometric-optical canopy reflectance model

A new method is described for the retrieval of fractional cover of large woody plants (shrubs) at the landscape scale using moderate resolution multi-angle remote sensing data from the Multiangle Imaging SpectroRadiometer (MISR) and a hybrid geometric-optical (GO) canopy reflectance model. Remote sensing from space is the only feasible method for regularly mapping woody shrub cover over large areas, an important application because extensive woody shrub encroachment into former grasslands has been seen in arid and semi-arid grasslands around the world during the last 150 years. The major difficulty in applying GO models in desert grasslands is the spatially dynamic nature of the combined soil and understory background reflectance: the background is important and cannot be modeled as either a Lambertian scatterer or by using a fixed bidirectional reflectance distribution function (BRDF). Candidate predictors of the background BRDF at the Sun-target-MISR angular sampling configurations included the volume scattering kernel weight from a Li-Ross BRDF model; diffuse brightness (ρ0) from the Modified Rahman-Pinty-Verstraete (MRPV) BRDF model; other Li-Ross kernel weights (isotropic, geometric); and MISR near-nadir bidirectional reflectance factors (BRFs) in the blue, green, and near infra-red bands. The best method was multiple regression on the weights of a kernel-driven model and MISR nadir camera blue, green, and near infra-red bidirectional reflectance factors. The results of forward modeling BRFs for a 5.25 km² area in the USDA, ARS Jornada Experimental Range using the Simple Geometric Model (SGM) with this background showed good agreement with the MISR data in both shape and magnitude, with only minor spatial discrepancies. The simulations were shown to be accurate in terms of both absolute value and reflectance anisotropy over all 9 MISR views and for a wide range of canopy configurations (r² = 0.78, RMSE = 0.013, N = 3969). Inversion of the SGM allowed estimation of fractional shrub cover with a root mean square error (RMSE) of 0.03 but a relatively weak correlation (r² = 0.19) with the reference data (shrub cover estimated from high resolution IKONOS panchromatic imagery). The map of retrieved fractional shrub cover was an approximate spatial match to the reference map. Deviations reflect the first-order approximation of the understory BRDF in the MISR viewing plane; errors in the shrub statistics; and the 12 month lag between the two data sets.     

The importance of instrument pointing accuracy for surface bidirectional reflectance distribution function mapping

An investigation into the magnitude of error introduced into bidirectional reflectance distribution functions (BRDFs) due to unqualified attitude perturbations and uncertainties in the pointing accuracy of the sensor is described. For a given perturbation in the sensing geometry the magnitude of the error in the measured BRDF depends on the shape of the BRDF, i.e. on the Sun-target-sensor geometry and on the surface cover type. Using typical BRDF data sets measured in the First International Satellite Land Surface Climatology Project (ISLSCP) Field Experiment (FIFE) it is shown that an error of 1° in the sensing geometry could result in approximately 1-5 per cent error in the BRDF. It is suggested that these errors could be reduced by sensing multi-angle data from stabilized platforms and/or by incorporating facilities for recording the sensor viewing geometry.     

Models for surface reflection of radiance and polarized radiance: Comparison with airborne multi-angle photopolarimetric measurements and implications for modeling top-of-atmosphere measurements

In this paper, we investigate the surface-atmosphere radiative interaction in application to the problem of aerosol satellite remote sensing over land. First, we test different models of the Bidirectional Reflectance and Polarization Distribution Function (BRDF and BPDF) for bare soil and vegetation surfaces using multi-angle, multi-spectral photopolarimetric airborne measurements of the Research Scanning Polarimeter (RSP). Then, we investigate the performance of different models of BRDF and BPDF for modeling top-of-atmosphere measurements. We have found that different BRDF models can describe the RSP measurements equally well. However, for soil surfaces, the different BRDF models show a different dependence on illumination geometry (solar zenith and azimuth angles), as well as a different dependence on viewing angle outside the range of RSP measurements. This implies that different models describe the surface-atmosphere interaction differently, leading for soil surfaces to differences in the top-of-atmosphere reflectance up to 4-5%, whereas at surface level the models agree within 2% for RSP illumination and measurement geometry. For vegetation, the different BRDF models show more similar dependence on illumination geometry, meaning that, in general, the differences in top-of-atmosphere reflectances are smaller than the differences in surface total reflectances. For the BPDF, we compare the empirical model of Nadal and Breon (1999) and the model developed by Maignan et al. (2009) with a newly developed model. The latter model compares better with RSP measurements. It was shown that, though all models have essentially different angular profiles at different illumination and viewing geometries, the difference of the top-of-atmosphere degree of linear polarization is less or is of the same order as the degree of linear polarization difference at the surface level taken at RSP illumination and measurement geometry. For the considered models, it can be up to 0.015 but is mostly below 0.005.     

Forward and inverse modelling of canopy directional reflectance using a neural network

This article explores the use of artificial neural networks for both forward and inverse canopy modelling. The forward neural modelling paradigm involved training a network for predicting the bidirectional reflectance distribution function (BRDF) of a canopy given the density of the trees, their height, crown shape, viewing, and illumination geometry. The neural network model was able to predict the BRDF of unseen canopy sites with 90% accuracy. Analysis of the signal captured by the model indicates that the canopy structural parameters, and illumination and viewing geometry, are essential for predicting the BRDF of vegetated surfaces. The inverse neural network model involved learning the underlying relationship between canopy structural parameters and their corresponding bidirectional reflectance. The inversion results show that the R2 between the network predicted canopy parameters and the actual canopy parameters was 0.85 for density and 0.75 for both the crown shape and the height parameters. The results of both forward and inverse modelling suggest that neural networks can model accurately the BRDF of vegetated canopies.     

Multi-temporal MODIS-Landsat data fusion for relative radiometric normalization, gap filling, and prediction of Landsat data

A semi-physical fusion approach that uses the MODIS BRDF/Albedo land surface characterization product and Landsat ETM+ data to predict ETM+ reflectance on the same, an antecedent, or subsequent date is presented. The method may be used for ETM+ cloud/cloud shadow and SLC-off gap filling and for relative radiometric normalization. It is demonstrated over three study sites, one in Africa and two in the U.S. (Oregon and Idaho) that were selected to encompass a range of land cover land use types and temporal variations in solar illumination, land cover, land use, and phenology. Specifically, the 30 m ETM+ spectral reflectance is predicted for a desired date as the product of observed ETM+ reflectance and the ratio of the 500 m surface reflectance modeled using the MODIS BRDF spectral model parameters and the sun-sensor geometry on the predicted and observed Landsat dates. The difference between the predicted and observed ETM+ reflectance (prediction residual) is compared with the difference between the ETM+ reflectance observed on the two dates (temporal residual) and with respect to the MODIS BRDF model parameter quality. For all three scenes, and all but the shortest wavelength band, the mean prediction residual is smaller than the mean temporal residual, by up to a factor of three. The accuracy is typically higher at ETM+ pixel locations where the MODIS BRDF model parameters are derived using the best quality inversions. The method is most accurate for the ETM+ near-infrared (NIR) band; mean NIR prediction residuals are 9%, 12% and 14% of the mean NIR scene reflectance of the African, Oregon and Idaho sites respectively. The developed fusion approach may be applied to any high spatial resolution satellite data, does not require any tuning parameters and so may be automated, is applied on a per-pixel basis and is unaffected by the presence of missing or contaminated neighboring Landsat pixels, accommodates for temporal variations due to surface changes (e.g., phenological, land cover/land use variations) observable at the 500 m MODIS BRDF/Albedo product resolution, and allows for future improvements through BRDF model refinement and error assessment.     

基于数据机理的植被叶面积指数遥感反演研究

定量获取地表植被高精度时序及空间覆盖的叶面积指数（Leaf Area Index, LAI）是生态监测及农业生产应用的重要研究内容。通过使用Moderate Resolution Imaging Spectroradiometer（MODIS）植被冠层多角度观测MOD09GA数据及叶面积指数MOD15A2数据,发展了一种参数化的叶面积指数遥感反演方法并完成了必要的检验分析。研究使用基于辐射传输理论的RossThick LiSparse Reciprocal（RTLSR）核驱动模型及Scattering by Arbitrarily Inclined Leaves with Hotspot（SAILH）模型进行植被冠层辐射特征的提取,使用Anisotropic Index （ANIX）异质性指数作为指示植被冠层二向反射分布Bidirectional Reflectance Distribution Function（BRDF）的辅助特征信息,发展了基于数据机理（Data-Based Mechanistic, DBM）的植被叶面积指数建模和估算方法。通过必要的林地、农作物、草地植被实验区反演及数值分析可得知：①时间序列多角度遥感观测数据结合数据机理的叶面积指数估算方法,可实现模型参数的时序动态更新,改进叶面积指数估算结果的时序完整性及精度。②异质性指数可以用做指示植被冠层二向反射分布特征信息,可降低因观测数据几何条件差异所导致的反演结果不确定情况,同时能够补充植被时序生长过程表现的植被结构变化等动态特征。经研究实践,可将算法应用于时空尺度的叶面积指数估算,并能够为生态、农业应用提供植被的高精度遥感监测指标。     

SIERRA: A new approach to atmospheric and topographic corrections for hyperspectral imagery

In mountainous areas, slope and altitude variations modulate the airborne sensed hyperspectral radiance image. A new algorithm, SIERRA, has been developed for atmospheric, relief and BRDF corrections in order to extract the surface reflectance in the form of bi-hemispherical albedo that does not depend on solar incidence and observation angles. The forward modeling efforts focus on the estimation of diffuse irradiance and upwelling diffuse radiance, and on the formulation of BRDF effects. The inversion scheme consists of four steps, that go deeper and deeper into the phenomena's complexity.To validate the model, reflectance images are assessed from radiance images simulated with different radiative transfer codes or forward models: MODTRAN4 in the case of homogeneous and flat ground, AMARTIS and SIERRA forward models for heterogeneous and mountainous cases. The surface reflectance is retrieved with a 5% relative error under standard acquisition conditions.SIERRA is applied to HyMap data acquired over the hilly landscape near Calanas, Spain. The hypercube reflectances are compared with those obtained using ATCOR4 and COCHISE. The benefit of relief correction is clearly demonstrated.     

Variations in reflectance with seasonality and viewing geometry: implications for semi-arid vegetation mapping with MISR data

The purpose of this article is to understand the effect of multi-temporal multi-angle data on vegetation community type mapping in desert regions. Based on data from the multi-angle imaging spectroradiometer (MISR), a set of 46 multi-temporal classification experiments were carried out in the Jornada Experimental Range in New Mexico, USA. Besides multi-angle observations, bidirectional reflectance distribution function (BRDF) model parameters were also used as input data for the classifications. The experiments used two widely accepted BRDF models, the Rahman–Pinty–Verstraete (RPV) model and the Ross-thin Li-sparse reciprocal (RTnLS) model. The experiments show that multi-temporal multi-angle classifications can yield a more accurate mapping than multi-temporal nadir classifications, and multi-temporal BRDF model parameters combined with a single nadir image can provide an accuracy roughly the same as all multi-temporal multi-angle observations for the vegetation mapping. These findings opened not only a path of reducing data dimensionality for multi-temporal multi-angle classifications, but also a way of merging products of both MISR and moderate resolution imaging spectroradiometer (MODIS) to improve semi-arid vegetation mapping.     

Remote sensing of aerosols optical thickness over various sites using SeaWiFS or VEGETATION and ground measurements

A key problem in aerosol retrieval is to distinguish between surface and atmospheric contributions to the variability in the satellite signal. A major contribution in the surface-related variability is caused by the non-Lambertian nature of the Earth surface reflectance and the fact that the illumination/observation geometry varies considerably between successive observations of the same area (with a polar orbiting sensor). In principle, if the surface boundary condition can be specified with sufficient accuracy by means of a bidirectional reflectance distribution function (BRDF), the two contributions can be unfolded and aerosol information retrieved. This approach has been tested using combined datasets made of satellite measured “top of atmosphere” (TOA) radiance and corresponding ground estimation of the aerosol optical thickness. Studying a time series of data, taking into account geometrical conditions and assuming the ground BRDF to be constant during the time period, a coupled surface/atmosphere model was used to investigate the retrieval of aerosol optical thickness (AOT) over several sites. By fitting a subset of satellite observations associated with ground photometer data, a best fit of BRDF model parameters could be determined. This surface characterization is then used to reduce the model unknowns to AOT only and thereby to permit its retrieval from the satellite data alone, by means of a simple inversion process. The study was conducted on three European AERONET sites and using satellite data from both the VEGETATION and Sea viewing Wide Field of view (SeaWiFS) sensors. In all cases, the AOT retrieved from satellite was in good agreement with the measurements.     

Analysis of broadband surface BRDFs derived from TOA SW CERES measurements for surfaces classified by the IGBP land cover

Most studies on the reflectance properties of the Earth's surface are addressed estimating the bidirectional reflectance distribution function (BRDF) of high spatial resolution and high spectral resolution satellite measurements. This article assesses the development of broadband (BB) BRDFs from radiances corresponding to large footprints classified according to the International Geosphere-Biosphere Programme (IGBP) land-cover classification. Top-of-atmosphere (TOA) shortwave (SW) CERES (Clouds and the Earth's Radiant Energy System) measurements are employed to invert the bidirectional reflectance factor (BRF) Rahman–Pinty–Verstraete (RPV) model for regions identified with the same IGBP type. The inversion of this non-linear parametric model is optimized to improve the computation efficiency and merged into a radiative transfer model to correct the surface radiances for the atmospheric effect. Analysis of the nature of the reflectance field simulated for several regions selected for every IGBP type determines whether the creation of general BRF models for surfaces defined by the same IGBP land cover is feasible. According to the results gathered in this study, the BB BRDFs for regions classified by the IGBP classification show values for the coefficient of variation (CV) between 3.5% and 44.1%. Consequently, the high differences achieved in the reflectance fields discourage the creation of BRDFs based on the IGBP land types.     

Efficient hyper-parameter determination for regularised linear BRDF parameter retrieval

ABSTRACT

Linear kernel driven models of the surface Bidirectional Reflectance Distribution Function (BRDF) are valuable tools for exploiting Earth observation data acquired at different sun–sensor geometries. Here we present a method that efficiently determines linear BRDF model weights using Tikhonov smoothing where the smoothing parameter λ is determined via a Generalized Singular Value Decomposition with the root mean square error prescribed depending on the MODIS band. We applied this method to twenty-six different deciduous broadleaf sites across an entire year using the MODIS Terra and Aqua reflectance data products. Kernel weights and white sky albedo derived from this GSVD method were generally consistent with those provided by the MCD43 data products. The GSVD derived results had less sample variability compared to the MCD43 data products, attributable to the assumed smoothness between kernel weights in the Tikhonov smoothing method. The GSVD technique consistently outperforms MCD43 in the reconstruction of observed MODIS reflectance data, of which retrievals from this method will do a better job of estimating albedo and normalizing data to specified geometries.     

Aerosol optical depth retrieval over snow using AATSR data

Aerosol observations over the Arctic are important because of the effects of aerosols on Arctic climate, such as their direct and indirect effects on the Earth's radiation balance and on snow albedo. Although information on aerosol properties is available from ground-based measurements, passive remote sensing using satellite measurements would offer the advantage of large spatial coverage with good temporal resolution, even though, due to light limitations, this is only available during the Arctic summer. However, aerosol optical depth (AOD) retrieval over the Arctic region is a great challenge due to the high reflectance of snow and ice and due to the high solar zenith angle. In this article, we describe a retrieval algorithm using Advanced Along-Track Scanning Radiometer (AATSR) data, a radiometer flying on the European Space Agency (ESA) Environmental Satellite (ENVISAT), which offers two views (near nadir and at 55° forward) at seven wavelengths in the visible thermal-infrared (VIS-TIR). The main idea of the Dual-View Multi-Spectral (DVMS) approach is to use the dual view to separate contributions to reflectance measured at the top of the atmosphere (TOA) due to atmospheric aerosol and the underlying surface. The algorithm uses an analytical snow bidirectional reflectance distribution function (BRDF) model for the estimation of the ratio of snow reflectances in the nadir and forward views, as well as an estimate of the atmospheric contribution to TOA reflectance obtained using the dark pixel method over the adjacent ocean surface, assuming that this value applies over nearby land surfaces in the absence of significant sources across the coastline. An iteration involving all four AATSR wavebands in the visible near-infrared (VIS-NIR) is used to retrieve the relevant information. The method is illustrated for AATSR overpasses over Greenland with clear sky in April 2009. Comparison of the retrieved AOD with AErosol Robotic Network (AERONET) data shows a correlation coefficient of 0.75. The AODs retrieved from AATSR using the DVMS approach and those obtained from AERONET data show similar temporal trends, but the AERONET results are more variable and the highest AOD values are mostly missed by the DVMS approach. Limitations of the DVMS method are discussed. The pure-snow BRDF model needs further correction in order to obtain a better estimation for mixtures of snow and ice.     

Large area mapping of southwestern forest crown cover, canopy height, and biomass using the NASA Multiangle Imaging Spectro-Radiometer

A rapid canopy reflectance model inversion experiment was performed using multi-angle reflectance data from the NASA Multi-angle Imaging Spectro-Radiometer (MISR) on the Earth Observing System Terra satellite, with the goal of obtaining measures of forest fractional crown cover, mean canopy height, and aboveground woody biomass for large parts of south-eastern Arizona and southern New Mexico (> 200,000 km²). MISR red band bidirectional reflectance estimates in nine views mapped to a 250 m grid were used to adjust the Simple Geometric-optical Model (SGM). The soil-understory background signal was partly decoupled a priori by developing regression relationships with the nadir camera blue, green, and near-infrared reflectance data and the isotropic, geometric, and volume scattering kernel weights of the LiSparse–RossThin kernel-driven bidirectional reflectance distribution function (BRDF) model adjusted against MISR red band data. The SGM's mean crown radius and crown shape parameters were adjusted using the Praxis optimization algorithm, allowing retrieval of fractional crown cover and mean canopy height, and estimation of aboveground woody biomass by linear rescaling of the dot product of cover and height. Retrieved distributions of crown cover, mean canopy height, and aboveground woody biomass for forested areas showed good matches with maps from the United States Department of Agriculture (USDA) Forest Service, with R² values of 0.78, 0.69, and 0.81, and absolute mean errors of 0.10, 2.2 m, and 4.5 tons acre^- 1 (10.1 Mg ha^- 1), respectively, after filtering for high root mean square error (RMSE) on model fitting, the effects of topographic shading, and the removal of a small number of outliers. This is the first use of data from the MISR instrument to produce maps of crown cover, canopy height, and woody biomass over a large area by seeking to exploit the structural effects of canopies reflected in the observed anisotropy patterns in these explicitly multiangle data.