基于Apache Spark的MODIS海表温度反演方法

为应对海量遥感影像快速计算的需求,通过对影像获取、算法和计算过程优化和改进,提出了一种基于Apache Spark并行计算框架的MODIS海表温度反演方法,实现了海量MODIS遥感影像的海表温度快速反演.应用四轮网络查询请求获取特定的时空范围影像数据,提高影像获取阶段的效率;应用简化算法参数、拟合过程变量改进海表温度劈窗算法,使之适合快速并行计算;应用弹性分布式数据集（RDD）窄依赖关系的优点,避免并行计算中的数据交换延迟.通过单机模式与集群模式对比实验,发现集成了并行计算框架的集群模式影像处理效率约为单机模式的10倍.研究结果表明了融合集群计算技术的海表温度反演过程有效提高了传统单机应用程序的处理效率.     

Identification and analysis of urban surface temperature patterns in Greater Athens, Greece, using MODIS imagery

Thermal infrared images are being acquired by satellites for more than two decades enabling studies of the human-induced Urban Heat Island (UHI) phenomenon. As a result, the requirement of the scientific community for fast and efficient methods for extracting and analyzing the thermal patterns from a vast volume of acquired data has emerged. The present paper proposes an innovative object-based image analysis procedure to extract thermal patterns for the quantitative analysis of satellite-derived Land Surface Temperature (LST) maps. The spatial and thermal attributes associated with these objects are then calculated and used for the analyses of the intensity, the position and the spatial extent of UHIs. A case study was conducted in the Greater Athens Area, Greece. More than 3000 LST images of the area acquired by MODIS sensor over a decade were analyzed. Three daytime hot-spots were identified and studied (Megara, Elefsina-Aspropyrgos and Mesogeia). They were all found to exhibit similar behavior, gradually increasing their maximum temperature during the summer season and reaching their maxima in mid-July. The hot-spots' thermal intensities compared to a suburban area were of 9-10 °C and were found to be highly correlated to their areal extent. During the night-time, Athens center developed a typical UHI spatially coinciding with the dense urban fabric. The nighttime maximum LST peaked (on average) at the end of July, two weeks later than the daytime surface patterns. The mean spatial extent of UHI in Athens was 55.2 km², whilst its mean intensity was 5.6 °C. The proposed automatic extraction process can be customized for other cities and potentially used for comparison of LST patterns and UHI behavior between different cities.     

Satellite-based estimation of surface vapor pressure deficits using MODIS land surface temperature data

Vapor Pressure Deficit (VPD) is a principle mediator of global terrestrial CO₂ uptake and water vapor loss through plant stomata. As such, methods to estimate VPD accurately and efficiently are critical for ecosystem and climate modeling efforts. Based on prior work relating energy partitioning, remotely sensed land surface temperature (LST), and VPD, we developed simple linear models to predict VPD using saturated vapor pressure calculated from MODIS LST at a number of different temporal and spatial resolutions. We developed and assessed the LST-VPD models using three data sets: (1) instantaneous and daytime average ground-based VPD and radiometric temperature from the Soil Moisture Experiments in 2002 (SMEX02); (2) daytime average VPD from AmeriFlux eddy covariance flux tower observations; and (3) estimated daytime average VPD from Global Surface Summary of Day (GSSD) observations. We estimated model parameters for VPD estimation both regionally (MOD11 A2) and globally (MOD11 C2) with RMSE values ranging from .32 to .38 kPa. VPD was overestimated along coastlines and underestimated in arid regions with low vegetation cover. Also, residuals were larger with higher VPDs because of the non-linear function of saturation vapor pressure with LST. Linear relationships were seen at multiple scales and appear useful for estimation purposes within a range of 0 to 2.5 kPa.     

Quality assessment of Landsat surface reflectance products using MODIS data

Surface reflectance adjusted for atmospheric effects is a primary input for land cover change detection and for developing many higher level surface geophysical parameters. With the development of automated atmospheric correction algorithms, it is now feasible to produce large quantities of surface reflectance products using Landsat images. Validation of these products requires in situ measurements, which either do not exist or are difficult to obtain for most Landsat images. The surface reflectance products derived using data acquired by the Moderate Resolution Imaging Spectroradiometer (MODIS), however, have been validated more comprehensively. Because the MODIS on the Terra platform and the Landsat 7 are only half an hour apart following the same orbit, and each of the 6 Landsat spectral bands overlaps with a MODIS band, good agreements between MODIS and Landsat surface reflectance values can be considered indicators of the reliability of the Landsat products, while disagreements may suggest potential quality problems that need to be further investigated. Here we develop a system called Landsat-MODIS Consistency Checking System (LMCCS). This system automatically matches Landsat data with MODIS observations acquired on the same date over the same locations and uses them to calculate a set of agreement metrics. To maximize its portability, Java and open-source libraries were used in developing this system, and object-oriented programming (OOP) principles were followed to make it more flexible for future expansion. As a highly automated system designed to run as a stand-alone package or as a component of other Landsat data processing systems, this system can be used to assess the quality of essentially every Landsat surface reflectance image where spatially and temporally matching MODIS data are available. The effectiveness of this system was demonstrated using it to assess preliminary surface reflectance products derived using the Global Land Survey (GLS) Landsat images for the 2000 epoch. As surface reflectance likely will be a standard product for future Landsat missions, the approach developed in this study can be adapted as an operational quality assessment system for those missions.     

MODIS-derived surface temperature of the Great Salt Lake

The surface temperature of Utah's hypersaline Great Salt Lake is examined between 2000 and 2007 using 3345 images from the Moderate Resolution Imaging Spectroradiometer (MODIS) on board the NASA Earth Observing System Terra and Aqua platforms. This study shows the utility of using a multi-year record of the easily accessible and fully processed MODIS thermal imagery to monitor spatial, diurnal, seasonal, and annual variations in the surface water temperature (SWT) of lakes where long-term in situ measurements are rarely available. A cloud-free Terra image is available on average every day during the summer and early fall, every other day during spring and late fall, and every third day during the winter. MODIS-derived lake SWT exhibits a cool bias (~ − 1.5 °C) relative to in situ temperature observations gathered from three buoys and a slowly-moving watercraft.The dominant SWT signal is the annual cycle (with a range of 26 °C and peak temperature in mid-July) while the diurnal range is as large as 4 °C during the spring season. Year-to-year variations in SWT are largest during the fall with over 1 °C anomalously warm (cold) departures from the 8-year monthly medians observed during fall 2001 (2006). The MODIS imagery provides an updated SWT climatology for operational weather forecasting applications (e.g., lake-effect snow storm prediction) as well as for input into operational and research numerical weather prediction models.     

Mapping insect defoliation in Scots pine with MODIS time-series data

Insect damage is a general problem that disturbs the growth of forests, causing economic losses and affecting carbon sequestration. Coarse-resolution data from satellites are potentially useful for national and regional mapping of forest damage, but the accuracy of these methods has not been fully examined. In this study, a method was tested for the mapping of defoliation in Scots pine [Pinus silvestris] forests in southeast Norway caused by the pine sawfly [Neodiprion sertifer], with the use of multi-temporal MODIS 16-day composite vegetation index data and the TIMESAT processing method. The damage mapping method used differences in summer mean values and angles of the seasonal profiles, indicating decreasing foliage density, to identify pixels that represent areas containing forest damage. In addition to 16-day NDVI the Wide Dynamic Range Vegetation Index (WDRVI) was tested. Damage areas were identified by classifying data into pixels representing damaged versus undamaged forest areas using a boolean combination of thresholded parameters. Classification results were evaluated against the change in LAI estimated from airplane LIDAR measurements, as an indicator of defoliation. The damage classifications detected 71% to 82% of the pixels with damage, and had kappa coefficients varying between 0.48 and 0.63, indicating some overestimation. This was due e.g. to failure to include clear-cut areas in the evaluation data. Damage classification with WDRVI only resulted in slight improvement compared to the NDVI. Only weak relationships were found between the LIDAR-estimated defoliation and the change parameters obtained from MODIS. Consequently, mapping of the degree of defoliation from MODIS was abandoned. In conclusion, the damage detection method based on MODIS data was found to be useful for locating insect damage, but not for estimating its intensity. Control of the detected damage areas using high-resolution remote sensing data, aerial survey, or fieldwork is recommended for accurate delineation in operational applications.     

Evaluation of MODIS land surface temperature data to estimate air temperature in different ecosystems over Africa

The estimation of near surface air temperature (Ta) is useful for a wide range of applications such as agriculture, climate related diseases and climate change studies. Air temperature is commonly obtained from synoptic measurements in weather stations. In Africa, the spatial distribution of weather stations is often limited and the dissemination of temperature data is variable, therefore limiting their use for real-time applications. Compensation for this paucity of information may be obtained by using satellite-based methods. However, the derivation of near surface air temperature (Ta), from the land surface temperature (Ts) derived from satellite is far from straight forward. Some studies have tried to derive maximum Ta from satellites through regression analysis but the accuracy obtained is quite variable according to the study. The main objective of this study was to explore the possibility of retrieving high-resolution Ta data from the Moderate Resolution Imaging Spectroradiometer (MODIS) Ts products over different ecosystems in Africa. First, comparisons between night MODIS Ts data with minimum Ta showed that MODIS nighttime products provide a good estimation of minimum Ta over different ecosystems (with (ΔTs − Ta) centered at 0 °C, a mean absolute error (MAE) = 1.73 °C and a standard deviation = 2.4 °C). Secondly, comparisons between day MODIS Ts data with maximum Ta showed that (ΔTs − Ta) strongly varies according to the seasonality, the ecosystems, the solar radiation, and cloud-cover. Two factors proposed in the literature to retrieve maximum Ta from Ts, i.e. the Normalized Difference Vegetation Index (NDVI) and the Solar Zenith Angle (SZA), were analyzed. No strong relationship between (ΔTs − Ta) and (i) NDVI and (ii) SZA was observed, therefore requiring further research on robust methods to retrieve maximum Ta.     

Mapping incident photosynthetically active radiation from MODIS data over China

Photosynthetically active radiation (PAR) is a key input parameter for almost all terrestrial ecosystem models, but the spatial resolution of current PAR products is too coarse to satisfy regional application requirements. In this paper, we present an operational system for PAR retrieval from MODIS data that is based on an idea proposed by [Liang, S., Zheng, T., Liu, R., Fang, H., Tsay, S. -C., & Running, S. (2006). Estimation of incident photosynthetically active radiation from Moderate Resolution Imaging Spectrometer data. Journal of Geophysical Research, 111, D15208. doi:10.1029/2005JD006730]. However, the operational system for PAR retrieval described here contains several improvements. The algorithm utilizes MODIS 1B data combining MODIS land surface products and BRDF model parameters products to directly estimate diffuse PAR, direct PAR and total PAR. Times-series data interpolation removes the noise and cloud contamination of land surface reflectance. PAR is retrieved by searching look-up tables calculated using a radiative transfer model. The system can automatically process MODIS 1B data to generate instantaneous and daily PAR. The instantaneous PAR products are compared with observational data from seven ChinaFLUX stations, and daily total PAR estimates are compared with those estimates of global radiation from 98 meteorological stations over China. The results indicate that this approach can produce reasonable PAR estimates, although this method overestimates PAR for low values of PAR.     

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.     

Heat waves measured with MODIS land surface temperature data predict changes in avian community structure

Heat waves are expected to become more frequent and severe as climate changes, with unknown consequences for biodiversity. We sought to identify ecologically-relevant broad-scale indicators of heat waves based on MODIS land surface temperature (LST) and interpolated air temperature data and assess their associations with avian community structure. Specifically, we asked which data source, time periods, and heat wave indices best predicted changes in avian abundance and species richness. Using mixed effects models, we analyzed associations between these indices and data from the North American Breeding Bird Survey in the central United States between 2000 and 2007 in four ecoregions and five migratory and nesting species groups. We then quantified avian responses to scenarios of severe, but commonly-occurring early, late, and summer-long heat waves. Indices based on MODIS LST data, rather than interpolated air temperatures, were more predictive of avian community structure. Avian communities were more related to 8-day LST exceedances (positive anomalies only); and were generally more sensitive to summer-long heat waves. Across the region, abundance, and to a lesser extent, species richness, declined following heat waves. Among the ecoregions, relationships were most consistently negative in the southern and montane ecoregions, but were positive in a more humid northern ecoregion. Among migratory groups, permanent resident species were the most sensitive, declining in abundance following a summer-long heat wave by 19% and 13% in the montane and southern ecoregions, respectively. Ground-nesting species, which declined in the south by 12% following a late summer heat wave, were more sensitive than avifauna overall. These results demonstrate the value of MODIS LST data for measuring ecologically-relevant heat waves across large regions. Ecologically, these findings highlight the importance of extreme events for avian biodiversity and the considerable variation in response to environmental change associated with different functional groups and geographic regions. The magnitude of the relationships between avian abundance and heat waves reported here raises concerns about the impacts of more frequent and severe heat waves in a warming climate.     

Determination of sea surface temperature at large observation angles using an angular and emissivity-dependent split-window equation

This paper proposes an angular and emissivity-dependent split-window equation that permits the determination of the sea surface temperature (SST) to a reasonable level of accuracy for any observation angle, including large viewing angles at the image edges of satellite sensors with wide swaths. This is the case of the MODIS radiometer both on EOS Terra/Aqua platforms, with observation angles of up to 65° at the surface, for which the split-window equation has been developed in this study. The algorithm takes into account the angular dependence of both the atmospheric correction (due to the increase of the atmospheric optical path with angle) and the emissivity correction (since sea surface emissivity (SSE) decreases with observation angle). Angular-dependent coefficients have been estimated for the atmospheric terms, and also an explicit dependence on the SSE has been included in the algorithm, as this parameter has values different to a blackbody surface for off-nadir angles, the SSEs also being dependent on surface wind speed. The proposed algorithm requires as input data at-sensor brightness temperatures for the split-window bands (31 and 32 of MODIS), the observation angle at each pixel, an estimate of the water vapor content (which is provided by the MODIS MOD07/MYD07 products) and accurate SSE values for both channels. The preliminary results show a good agreement between SSTs estimated by the proposed equation for off-nadir viewings of MODIS-Terra images and in situ SST measurements, with a root-mean square error (RMSE) of about ± 0.3 K, for which the MODIS SST product gives an RMSE larger than ± 0.7 K.     

Mapping PAR using MODIS atmosphere products

Instantaneous PAR (Photosynthetically Active Radiation), computed from atmospheric parameters from individual images from the MODIS sensors aboard the Terra and Aqua satellite platforms, is combined to derive daily integrated PAR and mapped to a local coordinate system. Compared to field observations, the daily integrated PAR values were shown to have average errors in the order of 5-8%, with individual estimation errors as high as 21%, but monthly averages showed much better correspondence with observations yielding averaged absolute errors of around 5%. The error appears to be mostly related to uncertainties in the MODIS aerosol retrieval accuracy. This accuracy and the medium spatial resolution of the PAR map compare very favourably to other sources of PAR data and make this a useful product in the improved assessment of vegetation dynamics.     

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.     

Validating MODIS land surface temperature products using long-term nighttime ground measurements

The Moderate Resolution Imaging Spectroradiometer (MODIS), onboard the NASA Terra and Aqua Earth Observing System satellites, provides multiple land surface temperature (LST) products on a daily basis. However, these products have not been adequately validated. This paper presents preliminary results of validating two MODIS Terra daily LST products, MOD11_L2 (version 4) and MOD07_L2 (version 4), using the FLUXNET and Carbon Europe Integrated Project (CarboEurope-IP) long-term ground measurements over eight vegetated sites. Since ground-measured LSTs were only available over one fixed point in each validation site, the study was carefully designed to mitigate the scale mismatch issue by using nighttime ground measurements concurrent to more than 1800 MODIS Terra overpasses.The preliminary results show that MOD11_L2 LSTs have smaller absolute biases and root mean squared errors (RMSE) than those of MOD07_L2 LSTs in most cases. The match of MOD11_L2 LSTs with ground measurements in the Brookings, Audubon, Canaan Valley, and Black Hills sites is good, yielding absolute biases less than 0.8 °C and RMSEs less than 1.7 °C. In the Fort Peck, Hainich, Tharandt, and Bondville sites, MOD11_L2 LSTs were underestimated by 2-3 °C. Biases in MOD11_L2 LSTs correlate to those in MOD07_L2 LSTs. Since the MOD07_L2 LST product is one of the input parameters to the MOD11_L2 LST algorithm, biases in MOD11_L2 LSTs may be influenced by biases in MOD07_L2 LSTs. The errors in both products depend weakly on sensor view zenith angle but are independent of surface air temperature, humidity, wind speed, and soil moisture.     

Spectral shape-based temporal compositing algorithms for MODIS surface reflectance data

Spectral similarity metrics have previously been used to select representative spectra from a class for use in spectral mixture modeling. Since the tasks of spectral selection for spectral mixture modeling and spectral selection for temporal compositing are similar, these metrics may have utility for temporal compositing. This paper explores the use of two spectral similarity metrics, endmember average root mean square error (EAR) and minimum average spectral angle (MASA), for constructing temporal composites. EAR and MASA compositing algorithms were compared against four previously used algorithms, including maximum NDVI, minimum view zenith angle, minimum blue, and median red. A total of 10 different algorithms were used to create 16-day composites of Moderate Resolution Imaging Spectroradiometer (MODIS) surface reflectance data over a 6-year period. Algorithm performance was assessed based on short-term temporal variability in spectral reflectance and in a selection of indices, both within a southwestern California study area and within five land-cover class subsets. EAR compositing produced the lowest variability for 4 out of 7 MODIS bands, as measured by the root mean square of time series residuals. MASA or EAR compositing produced the lowest root mean square residual values for all of the tested indices. To assess how compositing algorithms might affect remote sensing correlations with biophysical variables, correlations between indices calculated from different composites and live fuel moisture were compared. Correlations between indices and live fuel moisture were higher for shape-based composites compared with the standard composites.     

Comparison of land surface emissivity and radiometric temperature derived from MODIS and ASTER sensors

This study compared surface emissivity and radiometric temperature retrievals derived from data collected with the MODerate resolution Imaging Spectroradiometer (MODIS) and Advanced Spaceborne Thermal Emission Reflection Radiometer (ASTER) sensors, onboard the NASA's Earth Observation System (EOS)-TERRA satellite. Two study sites were selected: a semi-arid area located in northern Chihuahuan desert, USA, and a Savannah landscape located in central Africa. Atmospheric corrections were performed using the MODTRAN 4 atmospheric radiative transfer code along with atmospheric profiles generated by the National Center for Environmental Predictions (NCEP). Atmospheric radiative properties were derived from MODTRAN 4 calculations according to the sensor swaths, which yielded different strategies from one sensor to the other. The MODIS estimates were then computed using a designed Temperature-Independent Spectral Indices of Emissivity (TISIE) method. The ASTER estimates were derived using the Temperature Emissivity Separation (TES) algorithm. The MODIS and ASTER radiometric temperature retrievals were in good agreement when the atmospheric corrections were similar, with differences lower than 0.9 K. The emissivity estimates were compared for MODIS/ASTER matching bands at 8.5 and 11 μm. It was shown that the retrievals agreed well, with RMSD ranging from 0.005 to 0.015, and biases ranging from −0.01 to 0.005. At 8.5 μm, the ranges of emissivities from both sensors were very similar. At 11 μm, however, the ranges of MODIS values were broader than those of the ASTER estimates. The larger MODIS values were ascribed to the gray body problem of the TES algorithm, whereas the lower MODIS values were not consistent with field references. Finally, we assessed the combined effects of spatial variability and sensor resolution. It was shown that for the study areas we considered, these effects were not critical.     

Mapping plant functional types from MODIS data using multisource evidential reasoning

Reliable information about the geographic distribution and abundance of major plant functional types (PFTs) around the world is increasingly needed for global change research. Using remote sensing techniques to map PFTs is a relatively recent field of research. This paper presents a method to map PFTs from the Moderate Resolution Imaging Spectroradiometer (MODIS) data using a multisource evidential reasoning (ER) algorithm. The method first utilizes a suite of improved and standard MODIS products to generate evidence measures for each PFT class. The multiple lines of evidence computed from input data are then combined using Dempster's Rule of combination. Finally, a decision rule based on maximum support is used to make classification decisions. The proposed method was tested over the states of Illinois, Indiana, Iowa, and North Dakota, USA where crops dominate. The Cropland Data Layer (CDL) data provided by the United States Department of Agriculture were employed to validate our new PFT maps and the current MODIS PFT product. Our preliminary results suggest that multisource data fusion is a promising approach to improve the mapping of PFTs. For several major PFT classes such as crop, trees, and grass and shrub, the PFT maps generated with the ER method provide greater spatial details compared to the MODIS PFT. The overall accuracies increased for all the four states, with the biggest improvement occurring in Iowa from 51% (MODIS) to 64% (ER). The overall kappa statistic also increased for all the four states, with the biggest improvement occurring in Iowa from 0.03 (MODIS) to 0.38 (ER). The paper concludes with a discussion of several methodological issues pertaining to the further improvement of the ER approach.     

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.     

Percent surface water estimation from MODIS BRDF 16-day image composites

Four models for deriving percent surface water estimates were developed for use with MODIS16-day Bidirectional Reflectance Distribution Function (BRDF) corrected composite images. The models allow intra-annual surface water estimates to be produced with 1 km spatial resolution and an 8-day temporal resolution when applied to image composites from sensors on both the Aqua and Terra platforms. The surface water models are conceptually simple, relying on widely used indices (NDVI, NDWI, and tasseled cap), but computationally intensive. The models differ in the time and effort required to produce or acquire the inputs necessary for model training. The models were applied and tested in Yukon Flats National Wildlife Refuge, an area with varied surface water types including ponds, fens, and the Yukon River and its tributaries. Resulting accuracies peaked with an R² of approximately 0.625, and model accuracies were higher for pixels with higher percentages of water.     

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

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