Evaluation of ASTER and MODIS land surface temperature and emissivity products using long-term surface longwave radiation observations at SURFRAD sites

Land surface temperature (LST) and emissivity are key parameters in estimating the land surface radiation budget, a major controlling factor of global climate and environmental change. In this study, Terra Advanced Spaceborne Thermal Emission Reflection Radiometer (ASTER) and Aqua MODerate resolution Imaging Spectroradiometer (MODIS) Collection 5 LST and emissivity products are evaluated using long-term ground-based longwave radiation observations collected at six Surface Radiation Budget Network (SURFRAD) sites from 2000 to 2007. LSTs at a spatial resolution of 90 m from 197 ASTER images during 2000-2007 are directly compared to ground observations at the six SURFRAD sites. For nighttime data, ASTER LST has an average bias of 0.1 °C and the average bias is 0.3 °C during daytime. Aqua MODIS LST at 1 km resolution during nighttime retrieved from a split-window algorithm is evaluated from 2002 to 2007. MODIS LST has an average bias of − 0.2 °C. LST heterogeneity (defined as the Standard Deviation, STD, of ASTER LSTs in 1 × 1 km² region, 11 × 11 pixel in total) and instrument calibration error of pyrgeometer are key factors impacting the ASTER and MODIS LST evaluation using ground-based radiation measurements. The heterogeneity of nighttime ASTER LST is 1.2 °C, which accounts for 71% of the STD of the comparison, while the heterogeneity of the daytime LST is 2.4 °C, which accounts for 60% of the STD. Collection 5 broadband emissivity is 0.01 larger than that of MODIS Collection 4 products and ASTER emissivity. It is essential to filter out the abnormal low values of ASTER daily emissivity data in summer time before its application.     

Scaling of land surface temperature using satellite data: A case examination on ASTER and MODIS products over a heterogeneous terrain area

Land surface temperature (LST) is a key parameter in numerous environmental studies. Surface heterogeneity induces uncertainty in pixel-wise LST. Spatial scaling may account for the uncertainty, however, different approaches lead to differences in scaled values. Satellite-retrieved LST may be representative of the pixel-wise LST and useful for scaling analysis, but the limited accuracy of retrieved values adds uncertainty into the scaled values. Based on the Stefan-Boltzmann (S-B) law, this study proposed scaling approaches for LST over flat and relief areas to explore the combined uncertainties in scaling using satellite-retrieved data. To take advantage of simultaneous, multi-resolution observations at coincident nadirs by the Advanced Spaceborne Thermal Emission Reflection Radiometer (ASTER) and the MODerate-resolution Imaging Spectroradiometer (MODIS), LST products from these two sensors were examined for part of the Loess Plateau in China. 90-m ASTER LST data were scaled up to 1 km using the proposed approaches, and variation in the LST was generally reduced after scaling. Amongst the sources of uncertainties, surface heterogeneity (emissivity) and different scaling approaches resulted in very minor differences, with a maximum difference of 0.2 K for the upscaled LST. Terrain features, taken as an areal weighting factor, had negligible effects on the upscaled value. Limited accuracy of the retrieved LST was the major uncertainty. The overall LST increased 0.6 K on average with correction for terrain-induced angular effect and 0.4 K for both angular and adjacency effects over the study area. Accounting for terrain correction in scaling is necessary for rugged areas. With terrain correction, the upscaled ASTER LST achieved an agreement of − 0.1 ± 1.87 K and a root mean square error (RMSE) of 1.87 K overall with the 1-km MODIS LST rectified by Wan et al.'s approach [Wan, Z., Zhang, Y., Zhang Q., Li, Z.-L. (2002b), Validation of the land-surface temperature products retrieved from Terra Moderate Resolution Imaging Spectroradiometer data. Remote Sensing of Environment, 83, 163-180]. Refining the rectification approach resulted in a better agreement of − 0.2 ± 1.57 K and a RMSE of 1.58 K.     

Estimation of subpixel land surface temperature using an endmember index based technique: A case examination on ASTER and MODIS temperature products over a heterogeneous area

Land surface temperature (LST) is a key parameter in numerous environmental studies. Surface heterogeneity induces uncertainty in estimating subpixel temperature. To take an advantage of simultaneous, multi-resolution observations at coincident nadirs by the Advanced Spaceborne Thermal Emission Reflection Radiometer (ASTER) and the MODerate-resolution Imaging Spectroradiometer (MODIS), LST products from the two sensors were examined for a portion of suburb area in Beijing, China. We selected Soil-Adjusted Vegetation Index (SAVI), Normalized Multi-band Drought Index (NMDI), Normalized Difference Built-up Index (NDBI) and Normalized Difference Water Index (NDWI) as representative remote sensing indices for four land cover types (vegetation, bare soil, impervious and water area), respectively. By using support vector machines, the overall classification accuracy of the four land cover types with inputs of the four remote sensing indices, extracted from ASTER visible near infrared (VNIR) bands and shortwave infrared (SWIR) bands, reached 97.66%, and Kappa coefficient was 0.9632. In order to lower the subpixel temperature estimation error caused by re-sampling of remote sensing data, a disaggregation method for subpixel temperature using the remote sensing endmember index based technique (DisEMI) was established in this study. Firstly, the area ratios and statistical information of endmember remote sensing indices were calculated from ASTER VNIR/SWIR data at 990 m and 90 m resolutions, respectively. Secondly, the relationship between the 990 m resolution MODIS LST and the corresponding input parameters (area ratios and endmember indices at the 990 m resolution) was trained by a genetic algorithm and self-organizing feature map artificial neural network (GA-SOFM-ANN). Finally, the trained models were employed to estimate the 90 m resolution subpixel temperature with inputs of area ratios and endmember indices at the 90 m resolution. ASTER LST product was used for verifying the estimated subpixel temperature, and the verified results indicate that the estimated temperature distribution was basically consistent with that of ASTER LST product. A better agreement was found between temperatures derived by our proposed method (DisEMI) and the ASTER 90 m data (R² = 0.709 and RMSE = 2.702 K).     

ENVISAT/AATSR derived land surface temperature over a heterogeneous region

In this paper a method for evaluating land surface temperature (LST) algorithms over heterogeneous areas is presented. The evaluation was made for a set of 12 algorithms derived by using the split-window (SW) and dual-angle (DA) techniques for estimating sea and land surface temperature (SST and LST) from Advanced Along-Track Scanning Radiometer (AATSR) data. A validation of the proposed algorithms was carried out over a heterogeneous region of Morocco in the framework of the WATERMED (WATer use Efficiency in natural vegetation and agricultural areas by Remote sensing in the MEDiterranean basin) project. AATSR data and in situ measurements over this heterogenous region were compared by implementing a classification based strategy over a higher spatial resolution Landsat image. Three reference classes were considered when performing the classification from the Landsat image. Ground based measurements where then used to assign an effective surface radiometric temperature to each of these three classes. Finally, an averaging procedure based on class proportion was implemented for deriving surface radiometric temperature at the AATSR pixel scale. For this heterogeneous site, the results showed that LST can be obtained with a root mean-square error (RMSE) lower than 1.7 K from the split-window algorithms. Dual-angle algorithms, on the other hand, provided greater RMSE due to the different surfaces observed in the nadir and forward views. The results suggest that to retrieve LST from 1 km pixels over heterogeneous surfaces spatial averaging is required to improve accuracy on temperature estimation.     

Atmospheric correction of optical imagery from MODIS and Reanalysis atmospheric products

In this paper we analyze the differences obtained in the atmospheric correction of optical imagery covering bands located in the Visible and Near Infra-Red (VNIR), Short-Wave Infra-Red (SWIR) and Themal-Infrared (TIR) spectral regions when atmospheric profiles extracted from different sources are used. In particular, three sensors were used, Compact High Resolution Imaging Spectrometer (CHRIS), Advanced Spaceborne Thermal Emission and Reflection radiometer (ASTER) and Landsat5 Thematic Mapper (TM), whereas four atmospheric profiles sources were considered: i) local soundings launched near the sensor overpass time, ii) Moderate Resolution Radiometer (MODIS) atmospheric profiles product (MOD07), iii) Atmospheric Correction Parameter Calculator (ACPC) generated by the National Center for Environmental Prediction (NCEP) and iv) Modified Atmospheric Profiles from Reanalysis Information (MAPRI), which includes data from NCEP and National Center of Atmospheric Research (NCAR) Reanalysis project but interpolated to 34 atmospheric levels and resampled to 0.5° × 0.5°. MODIS aerosol product (MOD04) was also used to extract Aerosol Optical Thickness (AOT) values at 550 nm. Analysis was performed for three test dates (12th July 2003, 18th July 2004 and 13th July 2005) over an agricultural area in Spain. Results showed that air temperature vertical profiles were similar for the four sources, whereas dew point temperature profiles showed significant differences at some particular levels. Atmospheric profiles were used as input to MODTRAN4 radiative transfer code in order to compute atmospheric parameters involved in atmospheric correction, with the aim of retrieving surface reflectances in the case of VNIR and SWIR regions, and Land Surface Temperature (LST) in the case of the TIR region. For the VNIR and SWIR region, significant differences depending on the atmospheric profile used were not found, particularly in the Visible region in which the AOT content is the main parameter involved in the atmospheric correction. In the case of TIR, differences depending on the atmospheric profile used were appreciable, since in this case the main parameter involved in the atmospheric correction is the water vapor content, which depends on the vertical profile. In terms of LST retrieval from ASTER data (2004 test case), all profiles provided satisfactory results compared to the ones obtained when using a local sounding, with errors of 0.3 K for ACPC and MAPRI cases and 0.7 K for MOD07. When retrieving LST from TM data (2005 test case), errors for MOD07 and MAPRI were 0.6 and 0.9 K respectively, whereas ACPC provided an error of 2 K. The results presented in this paper show that the different atmospheric profile sources are useful for accurate atmospheric correction when local soundings are not available. In particular, MOD07 product provides atmospheric information at the highest spatial resolution, 5 km, although its use is limited from 2000 to present, whereas MAPRI provides historical information from 1970 to present, but at lower spatial resolution.     

Disaggregation of MODIS surface temperature over an agricultural area using a time series of Formosat-2 images

The temporal frequency of the thermal data provided by current spaceborne high-resolution imagery systems is inadequate for agricultural applications. As an alternative to the lack of high-resolution observations, kilometric thermal data can be disaggregated using a green (photosynthetically active) vegetation index e.g. NDVI (Normalized Difference Vegetation Index) collected at high resolution. Nevertheless, this approach is only valid in the conditions where vegetation temperature is approximately uniform. To extend the validity domain of the classical approach, a new methodology is developed by representing the temperature difference between photosynthetically and non-photosynthetically active vegetation. In practice, both photosynthetically and non-photosynthetically active vegetation fractions are derived from a time series of Formosat-2 shortwave data, and then included in the disaggregation procedure. The approach is tested over a 16 km by 10 km irrigated cropping area in Mexico during a whole agricultural season. Kilometric MODIS (MODerate resolution Imaging Spectroradiometer) surface temperature is disaggregated at 100 m resolution, and disaggregated temperature is subsequently compared against concurrent ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) data. Statistical results indicate that the new methodology is more robust than the classical one, and is always more accurate when fractional non-photosynthetically active vegetation cover is larger than 0.10. The mean correlation coefficient and slope between disaggregated and ASTER temperature is increased from 0.75 to 0.81 and from 0.60 to 0.77, respectively. The approach is also tested using the MODIS data re-sampled at 2 km resolution. Aggregation reduces errors in MODIS data and consequently increases the disaggregation accuracy.     

Separating surface emissivity and temperature using two-channel spectral indices and emissivity composites and comparison with a vegetation fraction method

The temperature-independent thermal infrared spectral indices (TISI) method is employed for the separation of land surface temperature (LST) and emissivity from surface radiances (atmospherically corrected satellite data). The daytime reflected solar irradiance and the surface emission at ∼3.8 μm have comparable magnitudes. Using surface radiances and a combination of day-night 2-channel TISI ratios, the ∼3.8 μm reflectivity is derived. For implementing the TISI method, coefficients for NOAA 9-16 AVHRR channels are obtained. A numerical analysis with simulated surface radiances shows that for most surface types (showing nearly Lambertian behavior) the achievable accuracy is ∼0.005 for emissivity (AVHRR channel-5) and ∼1.5 K for LST. Data from the European Centre for Medium-Range Weather Forecasts (ECMWF) is used for calculation of atmospheric attenuation. Comparisons are made over a part of central Europe on two different dates (seasons). Clouds pose a major problem to surface observations; hence, monthly emissivity composites are derived. Additionally, using TISI-based monthly composites of emissivities, a normalized difference vegetation index (NDVI)-based method is tuned to the particular study area and the results are intercompared. Once the coefficients are known, the NDVI method is easily implemented but holds well only for vegetated areas. The error of the NDVI-based emissivities (with respect to the TISI results) ranges between −0.038 and 0.032, but for vegetated areas the peak of the error-histogram is at ∼0.002. The algorithm for retrieving emissivity via TISI was validated with synthetic data. Due to the different spatial scales of satellite and surface measurements and the lack of homogeneous areas, which are representative for low-resolution pixels and ground measurements, ground-validation is a daunting task. However, for operational products ground-truth validation is necessary. Therefore, also an approach to identify suitable validation sites for meteorological satellite products in Europe is described.     

Investigating the effects of soil moisture on thermal infrared land surface temperature and emissivity using satellite retrievals and laboratory measurements

This study investigates the effects of soil moisture (SM) on thermal infrared (TIR) land surface emissivity (LSE) using field- and satellite-measurements. Laboratory measurements were used to simulate the effects of rainfall and subsequent surface evaporation on the LSE for two different sand types. The results showed that the LSE returned to the dry equilibrium state within an hour after initial wetting, and during the drying process the SM changes were uncorrelated with changes in LSE. Satellite retrievals of LSE from the Atmospheric Infrared Sounder (AIRS) and Moderate Resolution Imaging Spectroradiometer (MODIS) were examined for an anomalous rainfall event over the Namib Desert in Namibia during April, 2006. The results showed that increases in Advanced Microwave Scanning Radiometer (AMSR-E) derived soil moisture and Tropical Rainfall Measuring Mission (TRMM) rainfall estimates corresponded closely with LSE increases of between 0.08-0.3 at 8.6 µm and up to 0.03 at 11 µm for MODIS v4 and AIRS products. This dependence was lost in the more recent MODIS v5 product which artificially removed the correlation due to a stronger coupling with the split-window algorithm, and is lost in any algorithms that force the LSE to a pre-determined constant as in split-window type algorithms like those planned for use with the NPOESS Visible Infrared Imager Radiometer Suite (VIIRS). Good agreement was found between MODIS land surface temperatures (LSTs) derived from the Temperature Emissivity Separation (TES) and day/night v4 algorithm (MOD11B1 v4), while the split-window dependent products (MOD11B1 v5 and MOD11A1) had cooler mean temperatures on the order of 1-2 K over the Namib Desert for the month of April 2006.     

Ground measurements for the validation of land surface temperatures derived from AATSR and MODIS data

An experimental site was set up in a large, flat and homogeneous area of rice crops for the validation of satellite derived land surface temperature (LST). Experimental campaigns were held in the summers of 2002-2004, when rice crops show full vegetation cover. LSTs were measured radiometrically along transects covering an area of 1 km². A total number of four thermal radiometers were used, which were calibrated and inter-compared through the campaigns. Radiometric temperatures were corrected for emissivity effects using field emissivity and downwelling sky radiance measurements. A database of ground-based LSTs corresponding to morning, cloud-free overpasses of Envisat/Advanced Along-Track Scanning Radiometer (AATSR) and Terra/Moderate Resolution Imaging Spectroradiometer (MODIS) is presented. Ground LSTs ranged from 25 to 32 °C, with uncertainties between ± 0.5 and ± 0.9 °C. The largest part of these uncertainties was due to the spatial variability of surface temperature. The database was used for the validation of LSTs derived from the operational AATSR and MODIS split-window algorithms, which are currently used to generate the LST product in the L2 level data. A quadratic, emissivity dependent split-window equation applicable to both AATSR and MODIS data was checked as well. Although the number of cases analyzed is limited (five concurrences for AATSR and eleven for MODIS), it can be concluded that the split-window algorithms work well, provided that the characteristics of the area are adequately prescribed, either through the classification of the land cover type and the vegetation cover, or with the surface emissivity. In this case, the AATSR LSTs yielded an average error or bias of − 0.9 °C (ground minus algorithm), with a standard deviation of 0.9 °C. The MODIS LST product agreed well with the ground LSTs, with differences comparable or smaller than the uncertainties of the ground measurements for most of the days (bias of + 0.1 °C and standard deviation of 0.6 °C, for cloud-free cases and viewing angles smaller than 60°). The quadratic split-window algorithm resulted in small average errors (+ 0.3 °C for AATSR and 0.0 °C for MODIS), with differences not exceeding ± 1.0 °C for most of the days (standard deviation of 0.9 °C for AATSR and 0.5 °C for MODIS).     

Assessment of multi-resolution and multi-sensor data for urban surface temperature retrieval

Data from three thermal sensors with different spatial resolution were assessed for urban surface temperature retrieval over the Yokohama City, Japan. The sensors are Thermal Airborne Broadband Imager (TABI), Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) and MODerate resolution Imaging Spectroradiometer (MODIS). Two algorithms were developed for land surface temperature (LST) retrieval from TABI image and ASTER thermal infrared (TIR) channels 13 and 14. In addition, ASTER LST and MODIS LST products were also collected. All the LST images were assessed by analyzing the relationship between LST and normalized difference vegetation index (NDVI) and by spatial distributions of LST profiles, derived from typical transects over the LST images. In this study, a strong negative relationship between LST and NDVI has been demonstrated although the degree of correlation between NDVI and LST varies slightly among the different LST images. Cross-validation among the LST images retrieved from the three thermal sensors of different spatial resolutions indicates that the LST images retrieved from the 2 channel ASTER data and a single band TABI thermal image using our developed algorithms are reliable. The LST images retrieved from the three sensors should have different potential to urban environmental studies. The MODIS thermal sensor can be used for the synoptic overview of an urban area and for studying urban thermal environment. The ASTER, with its TIR subsystem of 90-m resolution, allows for a more accurate determination of thermal patterns and properties of urban land use/land cover types, and hence, a more accurate determination of the LST. In consideration of the high heterogeneity of urban environment, the TABI thermal image, with a high spatial resolution of 2 m, can be used for rendering and assessing complex urban thermal patterns and detailed distribution of LST at the individual house level more accurately.     

Influence of lake morphology and clarity on water surface temperature as measured by EOS ASTER

Relationships between lake morphometric parameters and nighttime lake surface temperatures were investigated in North American temperate lakes using the ASTER kinetic temperature (AST08) product. Nighttime ASTER kinetic temperature measurements were found to be a good analogue for nighttime surface temperatures. Linear regression between ASTER and buoy-measured temperatures in a test lake were better during the evening (R² = 0.98) than the day (R² = 0.90), presumably due to the greater influence of radiation and latent heat fluxes during daylight hours. Nighttime lake surface temperatures measured in three ASTER scenes were significantly correlated to logarithm of lake area, maximum lake depth, Secchi depth (a measure of lake clarity) and lake order (a measure of lake connection with surface drainage), during October and November. Nighttime lake surface temperatures were significantly correlated only with lake area in July. We hypothesize that morphology was more strongly related to surface temperature in the fall months due to lake turnover during that season. This study suggests that satellite derived thermal data may be useful for calculation of lake heat budgets and evaporation rates, provided surface temperatures are measured in well-mixed lakes.     

Characterising groundwater use by vegetation using a surface energy balance model and satellite observations of land surface temperature

This study presents a novel ‘model-data’ approach to detect groundwater-dependent vegetation (GDV), through differences in modelled and observed land surface temperatures (LST) in space and time. Vegetation groundwater use is inferred where modelled LST exceeds observed LST by more than a threshold determined from consideration of systematic and random errors in model and observations. Modelled LST was derived from a surface energy balance model and LST observations were obtained from Terra-MODIS thermal imagery. The model-data approach, applied in the Condamine River Catchment, Queensland, Australia, identified GDV coincident to existing mapping. GDV were found to use groundwater up to 48% of the time and for as many as 56 consecutive days. Under driest of conditions, groundwater was estimated to contribute up to 0.2 mm h⁻¹ to total ET for GDV. The ability to both detect the location and water-use dynamics of GDV is a significant advancement on previous remote-sensing GDV methods.     

Estimation of high-spatial resolution clear-sky longwave downward and net radiation over land surfaces from MODIS data

Surface downwelling longwave radiation (LWDN) and surface net longwave radiation (LWNT) are two components in the surface radiation budget. In this study, we developed new linear and nonlinear models using a hybrid method to derive instantaneous clear-sky LWDN over land from the Moderate Resolution Imaging Spectroradiometer (MODIS) TOA radiance at 1 km spatial resolution. The hybrid method is based on extensive radiative transfer simulation (physical) and statistical analysis (statistical). Linear and nonlinear models were derived at 5 sensor view zenith angles (0°, 15°, 30°, 45°, and 60°) to estimated LWDN using channels 27-29 and 31-34. Separate models were developed for daytime and nighttime observations. Surface pressure effect was considered by incorporating elevation in the models. The linear LWDN models account for more than 92% of variations of the simulated data sets, with standard errors less than 16.27 W/m² for all sensor view zenith angles. The nonlinear LWDN models explain more than 93% of variations, with standard errors less than 15.20 W/m². The linear and nonlinear LWDN models were applied to both Terra and Aqua TOA radiance and validated using ground data from six SURFRAD sites. The nonlinear models outperform the linear models at five sites. The averaged root mean squared errors (RMSE) of the nonlinear models are 17.60 W/m² (Terra) and 16.17 W/m² (Aqua), with averaged RMSE ~ 2.5 W/m² smaller than that of the linear models. LWNT was estimated using the nonlinear LWDN models and the artificial neural network (ANN) model method that predicts surface upwelling longwave radiation. LWNT was also validated using the same six SURFRAD sites. The averaged RMSEs are 17.72 (Terra) and 16.88 (Aqua) W/m²; the averaged biases are − 2.08 (Terra) and 1.99 (Aqua) W/m². The LWNT RMSEs are less than 20 W/m² for both Terra and Aqua observations at all sites.     

Development and validation of a technique for merging satellite derived aerosol optical depth from SeaWiFS and MODIS

Merging time series of satellite derived aerosol products from independent missions can support aerosol science by combining in a consistent way temporally overlapping data sets and by increasing data coverage. A merging technique applied to satellite aerosol optical depth τ_a is presented and tested with SeaWiFS and MODIS-Aqua data. The technique relies on least squares fitting of the available τ_a spectra onto a linear or second-order polynomial relation between log-transformed τ_a and wavelengths. First, the sensor specific products are compared with field observations collected by a sun-photometer installed on the Acqua Alta Oceanographic Tower in the northern Adriatic Sea. Mean absolute percentage differences are approximately 21% at 412 and 443 nm, and increase with wavelength, with large overestimates in the red and near-infrared bands. The mean absolute differences are typically 0.04. When inter-compared, the 2 satellite products agree well, with mean absolute percentage differences lower than 20% at all wavelengths and little bias. The results of the comparison of the merger outputs with the field data are well in line with the validation results of the sensor specific products, and are comparable for the various merging procedures. The benefits of merging in terms of data coverage are briefly illustrated for the Mediterranean basin.     

Assessing the coupling between surface albedo derived from MODIS and the fraction of diffuse skylight over spatially-characterized landscapes

In this effort, the MODerate Resolution Imaging Spectroradiometer (MODIS) (Collection V005) Bidirectional Reflectance Distribution Function (BRDF)/Albedo algorithm is used to retrieve instantaneous surface albedo at a point in time and under specific atmospheric conditions. These retrievals are then used to study the role that the fraction of diffuse skylight plays under realistic scenarios of anisotropic diffuse illumination and multiple scattering between the surface and atmosphere. Simulations of the sky radiance using the MODTRAN®5.1 radiative transfer model were performed under different aerosol optical properties, illumination conditions, and surface characteristics to describe these effects on surface albedo retrievals from MODIS. This technique was examined using a validation scheme over four measurement sites with varied aerosol levels and landscapes, ranging from croplands to tundra ecosystems, and over extended time periods. Furthermore, a series of geostatistical analyses were performed to examine the types of spatial patterns observed at each measurement site. In particular, Enhanced Thematic Mapper Plus (ETM+) retrievals of surface albedo were acquired to analyze the change in variogram model parameters as a function of increased window-size. Results were then used to assess the degree to which a given point measurement is able to capture the intrinsic variability at the scale of MODIS observations. Assessments of MODIS instantaneous albedos that account for anisotropic multiple scattering, over snow-free and snow-covered lands and at all diurnal solar zenith angles, show a slight improvement over the albedo formulations that treat the downwelling diffuse radiation as isotropic. Comparisons with field measurements show biases improving by 0.004-0.013 absolute units (root-mean-squared error) or 0.1%-2.0% relative error.     

Comparison of atmospheric aerosol climatologies over southwestern Spain derived from AERONET and MODIS

The southwestern area of Spain, by its geographical and climatological conditions, is a key location for the characterization of atmospheric aerosol properties. The present study is aimed at evaluating the reliability of satellite-based aerosol climatologies, as inferred from level 2 standard aerosol products such as the Terra-MODIS (Moderate Resolution Imaging Radiometer) MOD04 aerosol product, with an application over this region during the period 2000-2008.This evaluation is carried out by means of comparison with ground-based data from the AERONET station of El Arenosillo (Spain, 37.1N, 6.7W), which has been providing continuous data since 2000. The focus of this paper is the climatology of two aerosol optical parameters: the aerosol optical depth (AOD) and the Ångström exponent.AERONET ground-based measurements give an annual mean value of 0.16 ± 0.12 and a median of 0.12 for the AOD, and a mean value of 1.20 ± 0.47 for the Ångström exponent. The seasonal pattern is characterized by two maxima, the most important maximum occurs in summer months, and the other one in late-winter/early-spring. Lowest values appear in fall and winter, however, a local minimum is observed in July which is only detected with the long-term data series.The mean climatological AOD based on AERONET exhibits complex seasonal patterns (i.e. with multiple local extrema), which are not always captured by MODIS-based climatology. MODIS only reproduces low values of the AOD in winter and high values in summer, as well as the local minimum of July which is sharper when using over-land retrievals. The time series of the AOD retrieved from MODIS both over land and ocean are in relatively good agreement with the ground-based measurements, with a monthly overestimation of about 30% on average, and higher differences in spring. Seasonal patterns from MODIS are better reproduced over land than over ocean. The agreement between daily AERONET and MODIS, as assessed by linear regression, gives correlation coefficients above 80% and an intercept bias below 0.03.     

In-situ measured spectral directional emissivity of snow and ice in the 8–14 μm atmospheric window

The directional emissivity of snow and ice surfaces in the 8–14 μm thermal infrared (TIR) atmospheric window was determined from spectral radiances obtained by field measurements using a portable Fourier transform infrared spectrometer in conjunction with snow pit work. The dependence of the directional emissivity on the surface snow type (grain size and shape) was examined. We obtained emissivity spectra for five different surface types, i.e., fine dendrite snow, medium granular snow, coarse grain snow, welded sun crust snow, and smooth bare ice. The derived emissivities show a distinct spectral contrast at wavelengths λ = 10.5–12.5 μm which is enhanced with increasing the snow grain size. For example, emissivities at both 10.5 μm and 12.5 μm for the nadir angle were 0.997 and 0.984 for the fine dendrite snow, 0.996 and 0.974 for the medium granular snow, 0.995 and 0.971 for the coarse grain snow, 0.992 and 0.968 for the sun crust, and 0.993 and 0.949 for the bare ice, respectively. In addition, the spectral contrast exhibits a strong angular dependence, particularly for the coarser snow and bare ice, e.g., the emissivity at λ = 12.5 μm for the off-nadir angle of 75° reaches down to 0.927, 0.896, and 0.709 for the coarse grain snow, sun crust, and bare ice cases, respectively. The angular dependent emissivity spectra of the bare ice were quite consistent with the spectra predicted by the Fresnel reflectance theory. The observed results firmly demonstrate that the directional emissivity of snow in the TIR can vary depending upon the surface snow type. The high variability of the spectral emissivity of snow also suggests the possibility to discriminate between snow and ice types from space using the brightness temperature difference in the atmospheric window.     

A new model of gross primary productivity for North American ecosystems based solely on the enhanced vegetation index and land surface temperature from MODIS

Many current models of ecosystem carbon exchange based on remote sensing, such as the MODIS product termed MOD17, still require considerable input from ground based meteorological measurements and look up tables based on vegetation type. Since these data are often not available at the same spatial scale as the remote sensing imagery, they can introduce substantial errors into the carbon exchange estimates. Here we present further development of a gross primary production (GPP) model based entirely on remote sensing data. In contrast to an earlier model based only on the enhanced vegetation index (EVI), this model, termed the Temperature and Greenness (TG) model, also includes the land surface temperature (LST) product from MODIS. In addition to its obvious relationship to vegetation temperature, LST was correlated with vapor pressure deficit and photosynthetically active radiation. Combination of EVI and LST in the model substantially improved the correlation between predicted and measured GPP at 11 eddy correlation flux towers in a wide range of vegetation types across North America. In many cases, the TG model provided substantially better predictions of GPP than did the MODIS GPP product. However, both models resulted in poor predictions for sparse shrub habitats where solar angle effects on remote sensing indices were large. Although it may be possible to improve the MODIS GPP product through improved parameterization, our results suggest that simpler models based entirely on remote sensing can provide equally good predictions of GPP.     

A simple retrieval method for land surface temperature and fraction of water surface determination from satellite microwave brightness temperatures in sub-arctic areas

A strong linear relationship is found between Special Sensor Microwave/Imager (SSM/I) microwave (19 and 37 GHz) surface emissivities at horizontal and vertical polarizations over snow- and ice-free land surfaces. This allows retrieving the land surface emissivity and temperature from satellite microwave brightness temperatures after atmospheric corrections. Over the Canadian sub-arctic continental area, we show that the main factor modifying the emissivity is the fraction of water surface (FWS) within a pixel. Accordingly, a map of the fraction of water surface across the Canadian landmass is derived, given a correspondence within 6% as compared to the 1 km² Canadian National Topographic Database of water-covered areas. The microwave-derived surface temperatures are compared to synchronous in situ air and ground surface temperatures and also with independent satellite IR measurements over areas without snow or ice. Root mean square differences range between 2° and 3.5°, with mean bias error of the order of 1-3°. Better results are always obtained with the 37 GHz channel rather than with the 19 GHz channel. Over dense vegetation, the microwave-derived surface temperature is closer to the air temperature (at surface level) than to the ground temperature. The proposed simple retrieval algorithm, not sensitive to cloud cover, appears very useful for monitoring summer interannual or seasonal trends of the fraction of surface water, as well as the daily land surface temperature variation, which are very important parameters in environmental change analysis.     

Multitemporal-patch ensemble inversion of coupled surface-atmosphere radiative transfer models for land surface characterization

In the last few years, encouraging results using radiative transfer model inversion techniques were obtained for land biophysical variables retrieval. However, the inversion of radiative transfer models is a severely ill-posed problem that may lead to significant uncertainties in the biophysical variables estimates. Improvement of performances of the inversion process requires more information to be exploited including better radiative transfer models, exploitation of proper prior information on the distribution of canopy and atmosphere variables, knowledge of uncertainties in satellite measurements, as well as possible spatial and temporal constraints. In this study we focus on the use of coupled atmosphere-surface radiative transfer models (SMAC + SAIL + PROSPECT) to estimate some key biophysical variables from top of atmosphere canopy reflectance data. The inversion is achieved over an ensemble of pixels belonging to a spatial window where aerosol properties are supposed to be constant, and over a temporal window of few days where vegetation state is assumed not to vary. The ensemble inversion scheme accounting for the spatial and temporal constraints is described. Top of atmosphere reflectance observations are simulated for 13 bands within the visible and near infrared domains. The coupled model is inverted with a variational method implementation dedicated to solve very large inverse problems. It is based on the use of the adjoint model and a Quasi-Newton optimisation technique with BFGS update.The multitemporal-patch inversion approach exploiting the spatial and temporal constraints is compared to the classical instantaneous-local inversion applied on single pixel and date. The ‘ensemble’ approach shows significant performance improvements when retrieving aerosol optical thickness τ₅₅₀ and some canopy characteristics (LAI, LAI × C_ab and ALA). Conclusions are drawn on the interest of such approaches, and perspectives are given, with due attention to their applicability within operational algorithms.