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.     

基于ART模型的MODIS积雪反照率反演研究

积雪反照率是研究局地或全球的能量收支平衡和气候变化中的重要参数,遥感反演为积雪反照率的获取提供了便利的手段。积雪反照率大小主要取决于积雪的自身物理属性(雪粒径、形状和污染物等因子)以及天气状况,遥感反演反照率大多基于双向反射模型(BRDF),积雪BRDF模型常使用积雪辐射传输模型获得。采用考虑了雪粒径、粒子形状以及污染物影响的渐进辐射传输理论(ART)模型,建立了MODIS积雪反照率反演算法,得到了MODIS 8d合成积雪反照率产品。将此算法应用于具有均一积雪地表的格陵兰岛地区,并使用GC-Net实测数据进行了验证,反演的总均方根误差(RMSE)为0.018,相关系数(r)为0.83,结果表明考虑了积雪特性的ART模型能够较好地反演积雪反照率,而且反演需要的参数较少。     

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.     

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.     

An assessment of the performance of two snow kernels in characterizing snow scattering properties

The kernel-driven RossThick-LiSparseReciprocal (RTLSR) bidirectional reflectance distribution function (BRDF) model has been widely used in the quantitative remote sensing community. However, the performance of this model is challenged when modelling the optical scattering properties of pure snow surfaces. Recently, two snow kernels have been developed to improve the snow anisotropic reflectance in the kernel-driven RTLSR model framework. However, the performances of these two snow kernels must be assessed to identify their potential applications. Therefore, we assess the performances of these two kernels using various BRDF data sources. Our findings demonstrate their differences in several aspects. (1) These two kernels differ in characterizing the variability in BRDF shape as a function of the solar zenith angle (SZA). As the SZA increases, the shape of snow kernel derived by the asymptotic radiative transfer (ART) model (hereinafter named the ART method) changes from a dome shape to bowl shape, which agrees well with the simulation data of the bicontinuous photon tracking (bic-PT) model. The shape of the snow kernel proposed by Qu et al. (hereinafter named the Qu method) based on the Rahman-Pinty-Verstraete (RPV) model maintains a bowl shape for all SZAs. These differences in the kernel performances affect their abilities to fit snow BRDF data with different SZAs. (2) The corrected RTLSR models, with their respective snow kernels, are generally able to model the forward-scattering properties of snow surfaces compared with the original RTLSR model. However, the ART method performs better in capturing the BRDF variations in snow surfaces than the Qu method. This assessment provides an improved understanding of the performance of these two snow kernels and, thus, suggests further applications for the ART snow kernel in the kernel-driven BRDF model framework in the near future.     

Influence of land surface effects on MODIS aerosol retrieval using the BAER method over Korea

As satellite receiving signals are affected by complex radiative transfer processes in the atmosphere and on land surfaces, aerosol retrieval over land from space requires the ability to determine surface reflectance from the remote measurements. To use the Bremen Aerosol Retrieval (BAER) method for aerosol optical thickness (AOT) retrieval over land at a spatial scale of 1×1 km² from Moderate Resolution Imaging Spectroradiometer (MODIS) data, a linear mixing model with a vegetation index was used to calculate surface reflectances. As the vegetation index is affected by the aerosol present in the atmosphere, an empirical linear relationship between short wavelength infrared (SWIR) channel reflectance and visible reflectance was estimated to calculate a modified aerosol free vegetation index (AFRI) value. Based on a modified AFRI obtained from MODIS SWIR channel reflectance, an improved linear mixing model was applied for aerosol retrieval. A comparison of results between calculated and apparent surface reflectance was satisfactory, with a linear fit slope above 0.94, correlation coefficients above 0.84, and standard deviation below 0.008 for the study area. These results can therefore be used for improved aerosol retrieval over land by the BAER method with MODIS Level 1 data.     

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.     

Algorithm for retrieval of aerosol optical properties over the ocean from the Geostationary Ocean Color Imager

An aerosol retrieval algorithm for the first Geostationary Ocean Color Imager (GOCI) to be launched in March 2010 onboard the Communication, Ocean, and Meteorological Satellite (COMS) is presented. The algorithm retrieves aerosol optical depth (AOD), fine-mode fraction (FMF), and aerosol type in 500 m × 500 m resolution. All the products are retrieved over clear water which is defined by surface reflectance ratio between 640 nm and 860 nm (SRR) less or equal to 2.5, while only AOD is retrieved over turbid water (SRR > 2.5) due to high surface reflectance. To develop optimized algorithm for the target area of GOCI, optical properties of aerosol are analyzed from extensive observation of AERONET sunphotometers to generate lookup table. Surface reflectance of turbid water is determined from 30-day composite of Rayleigh- and gas corrected reflectance. By applying the present algorithm to MODIS top-of-the atmosphere reflectance, three different aerosol cases dominated by anthropogenic aerosol contains black carbon (BC), dust, and non-absorbing aerosol are analyzed to test the algorithm. The algorithm retrieves AOD, and size information together with aerosol type which are consistent with results inferred by RGB image in a qualitative way. The comparison of the retrieved AOD with those of MODIS collection 5 and AERONET sunphotometer observations shows reliable results. Especially, the application of turbid water algorithm significantly increases the accuracy in retrieving AOD at Anmyon station. The sensitivity study between MODIS and GOCI instruments in terms of relative sensitivity and scattering angle shows promising applicability of the present algorithm to future GOCI measurements.     

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.     

Retrieval of aerosol optical thickness from HJ-1 CCD data based on MODIS-derived surface reflectance

The successfully launched Huanjing-1 (HJ-1) satellite by China in 2008 provides a new source of data for monitoring the environment. In this article, we develop a new algorithm for retrieving the aerosol optical thickness (AOT) using HJ-1 charge-coupled device (CCD) data with the assistance of the Moderate Resolution Imaging Spectroradiometer (MODIS) surface reflectance data and the bidirectional reflectance distribution function (BRDF) data products. This algorithm is then used to retrieve AOT in a delta region of the Yangtze River. The retrieved results are assessed for their accuracy by comparison with ground-measured data using sun photometers. Comparison of such derived AOT with in situ AOT measured using sun photometers indicates a root mean squared error (RMSE) of 0.123, and their regression relation has a correlation coefficient of 0.896 that is statistically significant at the 0.01 level. Such a relatively high level of retrieval accuracy suggests that HJ-1 CCD data can be used competently and effectively to retrieve AOT with the assistance of MODIS products that are used to construct the surface reflectance model. This study successfully demonstrates the feasibility of synergistically retrieving AOT from data acquired by different sensors. The lower dependence on data from a sole source means that the retrieval is less restrictive by data availability.     

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.     

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.     

A new climatology of aerosols in the middle and upper stratosphere by alternative analysis of GOMOS observations during 2002–2006

Recent works have shown how the study of stratospheric background aerosol (i.e. in periods uninfluenced by major volcanic eruption) seems more complex as it is now performed by more accurate means. We propose a re-analysis of Global Ozone Monitoring by Occultation of Stars (GOMOS) level 1b data for the period August 2002–July 2006, using the LPC2E processor algorithm, which was developed for the retrieval of aerosol extinction in the middle and upper stratosphere. The main differences with regard to the ‘official’ algorithm are the correction of chromatic scintillation, the spectral domain (which has been restricted to the 400–700 nm region), the use of a Differential Optical Absorption Spectroscopy (DOAS) method for species retrieval, and the use of a fourth-order polynomial to reproduce the wavelength dependence of extinction. Since GOMOS observations are performed using stars of different magnitude and colour, discrepancy in signal-to-noise ratio between several profiles exists, and a data selection concerning standard deviation of aerosol extinction and other parameters becomes necessary. In the middle stratosphere, aerosol extinction profiles obtained with the LPC2E processor seem to be in better agreement with the SAGE III observations and sparse balloon-borne measurements than the ‘official products’. We present global coverage of the 500 nm extinction values from around 15–60 km, and the wavelength dependence in the 400–675 nm spectral range which gives information about the nature of the particles. The well-known layer of liquid aerosols can be observed in the lower stratosphere, where the value of extinction is greater for blue than for red wavelengths, as is typical for small droplets. In the middle stratosphere, relatively high extinction values are found, probably due to the presence of solid particles above 30 km at all latitudes. The presence of soot and interplanetary material in the middle atmosphere is discussed, as well as seasonal patterns common to the several years of analysis, such as the stratospheric cleansing of aerosols above 30 km during polar winters.     

Development of an approach for generation of temporally complete daily nadir MODIS reflectance time series

Consistent, spatially and temporally complete reflectance time series are required for reliable terrestrial monitoring. The Moderate Resolution Imaging Spectroradiometer (MODIS), like other polar-orbiting wide field of view satellite sensors, can provide global observations on a nearly daily basis, but the sparseness of valid observations due to cloud, residual atmospheric effects, and sensor anomalies, may result in gaps in the derived reflectance time series. This paper presents an approach for the generation of temporally complete daily MODIS 500 m nadir view BRDF-adjusted reflectance (NBAR) time series. The research is illustrated and assessed quantitatively using two years of cloud and snow screened, daily MODIS Terra and Aqua reflectance data at four sites in Africa, and demonstrated for phenology monitoring using NBAR derived NDVI time series. The components of the approach include: 1) an outlier detection algorithm to remove residual anomalous daily observations undetected in the upstream processing, 2) the dynamic generation of NBAR time series on a daily basis when seven or more observations are available for a day under consideration over a 16-day period, and 3) the means to gap fill the NBAR time series where less than 7 observations are available. The MODIS Ross-Thick/Li-Sparse-Reciprocal BRDF model is used with a rolling approach whereby a 16-day BRDF inversion window is moved on a daily overlapping basis to provide more reliable outlier detection and daily NBAR. NBAR gap filling in periods of missing observations is investigated using static land cover specific archetype BRDF parameters and using BRDF parameters defined adaptively from the temporally closest 16-day periods with 7 or more observations. Scaling factor estimators using ordinary least squares (OLS) and median-based robust least squares regression are investigated, and the robust method is demonstrated to provide on average temporally more coherent gap filled NBAR values. For regions with persistent clouds, the utility of the adaptive NBAR gap filling method is demonstrated to be severely limited due to the decreased likelihood that the surface BRDF at each gap can be described reliably. The reliability of the NBAR gap filling methodology is evaluated statistically using a cross-validation approach. For the small number of study site considered, the adaptive method is shown to provide more accurate results than the archetype method when there are more than an average of ~ 4-5 observations per 16-day window, or when a gap day is on average less than about 30 days from a 16-day period with 7 or more observations. The resulting gap free daily NBAR time series and derived daily NBAR NDVI generated by the approach is shown to capture phenological variations in a coherent temporally consistent manner, suggesting that it is a fruitful avenue for future research and validation.     

Validation of MODIS derived aerosol optical depth over the Yangtze River Delta in China

MODIS derived aerosol optical depths (AODs) at 550 nm are compared with sunphotometer CE318 measurements at 7 sites located at Yangtze River Delta (YRD) in China from July to October, 2007. The evaluation result indicates that MODIS AODs (Collection 5, C005) are in good agreement with those from CE318 in dense vegetation regions, but show more differences in those regions with complex underlying surface (such as at lake water and urban surface sites). Reasons for these differences are discussed after removing cases with significant errors caused by validation scheme. The final validation result shows that MODIS AODs are in good agreement with CE318 with a correlation coefficient of 0.85 and RMS of 0.15. 90% of MODIS cases fall in the range of Δτ = ± 0.05 ± 0.20τ, indicating MODIS aerosol retrieval algorithm, aerosol models and surface reflectance estimate are generally suitably reasonable for aerosol retrieval in YRD. However, MODIS AODs show a systemic errors with fitted line of y = 0.75x + 0.13, indicating underestimation of AOD when aerosol loadings are high. Aerosol models and surface reflectance estimations are dominant sources of MODIS aerosol retrieval errors.     

Aerosol mapping over land with imaging spectroscopy using spectral autocorrelation

A new method for aerosol retrieval over land is proposed that makes explicit use of the contiguous, high-resolution spectral coverage of imaging spectrometers. The method is labelled Aerosol Retrieval by Interrelated Abundances (ARIA) and is based on unmixing of the short-wave infrared sensor signal by region-specific endmembers, assuming low aerosol radiative influence in this spectral region. Derived endmember abundances are transferred to the visible part of the spectrum in order to approximate surface reflectance where aerosol influence is generally strongest. Spectral autocorrelation of surface spectra is a precondition for ARIA and demonstrated using a reference spectrum database. The re-mixed surface reflectance is used as input quantity for the inversion of aerosol optical depth τ_a at 0.55 µm wavelength on a pixel basis. Except for the choice of endmembers and the atmospheric vertical profile, no a priori assumptions on the image scene are required. The potential of the presented method for aerosol retrieval is demonstrated for an Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) scene, collected in California in 2000. Comparisons with existing aerosol retrieval methods showed encouraging results in terms of achieved spatial smoothness and degree of uncertainty of aerosol optical depth across the scene.     

Enhanced algorithm for estimating snow depth from geostationary satellites

Observations in the visible and infrared spectral bands from the Imager instrument onboard Geostationary Operational Environmental Satellite (GOES) have been used to derive snow depth. The technique makes use of correlation between depth of the snow pack and satellite-derived subpixel fractional snow cover. Previous efforts to infer snow depth from satellite data with this technique were focused on grasslands and croplands, where the snow depth/snow fraction relationship is most pronounced. In this paper we improve the retrieval algorithm to extend snow depth estimates to forested areas. The enhanced algorithm accounts for the tree cover fraction and for the type of forest, deciduous or coniferous.The developed technique was used to derive maps of snow depth over mid-latitude areas of North America during winter seasons of 2003-2004 and 2004-2005. Satellite-based snow depth maps were produced daily at 4 km spatial resolution. To validate the retrievals we compared them with surface observations of snow depth and with the snow depth analysis prepared at the NOAA National Operational Hydrological Remote Sensing Center (NOHRSC). The estimated retrieval error was about 30% for snow depths below 30 cm and increased to 50% for snow depths ranging from 30 to 50 cm. Snow depth retrievals were limited to scenes with less than 80% deciduous forest cover fraction and less than 50% needle leaf forest cover.     

Aerosol optical depth spatio-temporal characterization over the Canadian BOREAS domain

A complete set of Advanced Very High Resolution Radiometer (AVHRR) data (75 images) is used to retrieve aerosol optical depth (AOD) over dense vegetation and over lake water in the visible AVHRR channel. The present approach for remote sensing of aerosols from the National Oceanic and Atmospheric Administration (NOAA)-11 AVHRR sensor is based on the detection of atmospherically dominated signals over dark surface covers such as dense dark vegetation (DDV). Such targets were identified using the reflective portion of the middle-wave AVHRR channel 3 signal. When a fixed DDV surface reflectance is subtracted from the observed reflectance after correction for all other atmospheric effects, the remaining part, which is due to aerosols, is inverted to derive aerosol optical thickness using a look-up table (LUT) similar to that used in water surface inversion. The algorithm was applied to the daily afternoon NOAA-11 AVHRR (1?km×1?km) data acquired from the end of May to mid-August 1994 over the Canadian 1000?km×1000?km Boreal Ecosystem Atmosphere Study (BOREAS) domain. A validation analysis using five ground-based Sun photometers within the studied area shows the good performance of the retrieval algorithm. The approach allows detailed analysis of the AOD spatio-temporal behaviour at the regional scale useful for climate and transport model validation.     

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.     

Robust estimation of BRDF model parameters

The effect of the Bidirectional Reflectance Distribution Function (BRDF) is one of the most important factors in correcting the reflectance obtained from remotely sensed data. Estimation of BRDF model parameters can be deteriorated by various factors; contamination of the observations by undetected subresolution clouds or snow patches, inconsistent atmospheric correction in multiangular time series due to uncertainties in the atmospheric parameters, slight variations of the surface condition during a period of observation, for example due to soil moisture changes, diurnal effects on vegetation structure, and geolocation errors [Lucht and Roujean, 2000]. In the present paper, parameter estimation robustness is examined using Bidirectional Reflectance Factor (BRF) data measured for paddy fields in Japan. We compare both the M-estimator and the least median of squares (LMedS) methods for robust parameter estimation to the ordinary least squares method (LSM). In experiments, simulated data that were produced by adding noises to the data measured on the ground surface were used. Experimental results demonstrate that if a robust estimation is sought, the LMedS method can be adopted for the robust estimation of a BRDF model parameter.