基于中国通量网的MODIS短波反照率验证与分析

遥感地表反照率产品的验证与分析是将其应用于环境研究的基础。采用中国通量网的地表实测短波反照率数据对MODIS反照率产品进行对比和分析,针对选取的8个地面站点,提取了2004年的MODIS反照率产品并进行验证。这些站点的植被覆盖情况涵盖了草地、森林和农业用地。结果显示MODIS在多数情况下能提供准确的地表反照率产品。针对各个站点的误差、均方根误差、相关系数分析都显示了这个结果,总体反演误差在0.002左右。较大的误差出现在有冰雪影响的时候,排除受积雪影响的数据,总体均方根误差可达0.028。分析了引起误差的原因并提出了改进意见。     

Validation of the collection 5 MODIS FPAR product in a heterogeneous agricultural landscape in arid Uzbekistan using multitemporal RapidEye imagery

The fraction of photosynthetically active radiation (FPAR) absorbed by a vegetation canopy is an important variable for global vegetation modelling and is operationally available from data of the Terra Moderate Resolution Imaging Spectroradiometer (MODIS) satellite sensor starting from the year 2000. Product validation is ongoing and important for constant product improvement, but few studies have investigated the specific accuracy of MODIS FPAR using in situ measurements and none have focused on agricultural areas. This study therefore presents a validation of the collection 5 MODIS FPAR product in a heterogeneous agricultural landscape in western Uzbekistan. High-resolution FPAR maps were compiled via linear regression between in situ FPAR measurements and the RapidEye normalized difference vegetation index (NDVI) for the 2009 season. The data were aggregated to the MODIS scale for comparison. Data on the percentage cover of agricultural crops per MODIS pixel allowed investigation of the impact of spatial heterogeneity on MODIS FPAR accuracy. Overall, the collection 5 MODIS FPAR overestimated RapidEye FPAR between approximately 6% and 15%. MODIS quality flags, the underlying biome classification and spatial heterogeneity were investigated as potential sources of error. MODIS data quality was very good in all cases. A comparison of the MODIS land-cover product with high-resolution land-use classification revealed a significant misclassification by MODIS. Yet, we found that the overestimation of MODIS FPAR is independent of classification accuracy. The results indicate that the amount of background information, present even in the most homogeneous pixels (～70% crop cover), is most likely the reason for the overestimation. The behaviour of pure pixels could not be investigated due to a lack of appropriate pixels.     

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.     

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

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

Evaluation of the MODIS (MOD10A1) daily snow albedo product over the Greenland ice sheet

This study evaluates the performance of the beta-test MODIS (MOD10A1) daily albedo product using in situ data collected in Greenland during summer 2004. Results indicate the beta-test product tracks the general seasonal variability in albedo but exhibits significant more temporal variability than observed at the stations. This may indicate problems with the cloud detection algorithm, and/or failure of the BRDF model to adequately model the bidirectional reflectance of snow. Comparisons with in situ observations at five automatic weather stations in Greenland indicate an overall RMSE of 0.067 for the Terra instrument and an RMSE of 0.075 on Aqua. The Terra-retrieved-albedo are slightly better correlated with the in situ data than the Aqua retrievals (r = 0.79 versus r = 0.77). Comparisons were also made between the MODIS daily albedo product and the MODIS 16-day albedo product (MOD43B3). Results indicate general correspondence between the two products, with better agreement found using the Terra-retrieved-albedo than the Aqua-retrieved albedo. The reason for the differences in albedo between the Aqua and Terra satellites remains unclear. At the stations examined, both the Terra and Aqua retrievals were made at nearly the same time of the day and therefore the differences in albedo between the satellites cannot be explained by differences in solar illumination. Finally, the albedo derived using MODIS data and the direct estimation algorithm (DEA) was also compared with 2004 Greenland in situ data. Results from this comparison suggest that the DEA performs well as long as the solar zenith angle of the observation is not greater than 70°.     

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.     

Seasonal variability in aerosol optical depth over India: a spatio-temporal analysis using the MODIS aerosol product

Aerosols are one of the key components of climate systems. They absorb and scatter both solar and terrestrial radiation and produce strong surface as well as atmospheric radiative forcing effects. Aerosol climatology includes the measurement of light extinction by aerosol scattering and absorption, by procedures such as aerosol optical depth (AOD), angstrom exponent (α), single-scattering albedo (ω), and size distribution. This article analyses the dynamics of seasonal AOD over the Indian subcontinent from 2001 to 2009 using the MODIS level 2 data set. The analysis carried out for winter, pre-monsoon, monsoon, and post-monsoon seasons is based on 8 days’ composite AOD data for selected months representative of each season. The spatial variability of AOD has been shown to be 0.47 μm and 0.66 μm for fine- and coarse-mode aerosols, respectively, which illustrates the principle of relative difference. The dynamics of seasonally averaged AOD over the period under study represent an increasing tendency from 0.20 to 0.37 at 0.47 μm and from 0.16 to 0.26 at 0.66 μm during winter (2003–2009), whereas AOD in the pre-monsoon season ranged from 0.24 to 0.16 at 0.47 μm and from 0.24 to 0.16 at 0.66 μm (2005–2009). The monsoon season yielded an AOD of less than 0.15 throughout the study period, and the post-monsoon season recorded an increasing tendency from 0.18 to 0.29 at 0.47 μm and from 0.16 to 0.19 at 0.66 μm (2005–2009), reflecting a similar trend to that of the winter AOD curve. The spatial distribution of AOD shows that the northern part of India – especially the Indo-Gangetic plain – remains most affected by high AOD throughout the year. Such high AOD can be attributed to increasing anthropogenic emission due to an ever-increasing population, and urban, industrial, and other economic activities causing high concentrations of fine-mode organic and inorganic aerosol particles, along with coarse soil and mineral dust over the Indo-Gangetic plain.     

An assessment of the MODIS collection 5 leaf area index product for a region of mixed coniferous forest

Canopy leaf area index (LAI), defined as the single-sided leaf area per unit ground area, is a quantitative measure of canopy foliar area. LAI is a controlling biophysical property of vegetation function, and quantifying LAI is thus vital for understanding energy, carbon and water fluxes between the land surface and the atmosphere. LAI is routinely available from Earth Observation (EO) instruments such as MODIS. However EO-derived estimates of LAI require validation before they are utilised by the ecosystem modelling community. Previous validation work on the MODIS collection 4 (c4) product suggested considerable error especially in forested biomes, and as a result significant modification of the MODIS LAI algorithm has been made for the most recent collection 5 (c5). As a result of these changes the current MODIS LAI product has not been widely validated. We present a validation of the MODIS c5 LAI product over a 121 km² area of mixed coniferous forest in Oregon, USA, based on detailed ground measurements which we have upscaled using high resolution EO data. Our analysis suggests that c5 shows a much more realistic temporal LAI dynamic over c4 values for the site we examined. We find improved spatial consistency between the MODIS c5 LAI product and upscaled in situ measurements. However results also suggest that the c5 LAI product underestimates the upper range of upscaled in situ LAI measurements.     

MODIS、MISR与POLDER 3种全球地表反照率卫星反演产品的比较与分析

对MODIS、MISR和POLDER 3种由多角度卫星观测反演得到的全球地表反照率数据(无冰雪覆盖区域)短波波段(SW,0.3~5 μm)与可见光波段(VIS,0.3~0.7 μm)的黑空地表反照率(DHR)进行了相互比较。3种产品6年平均的全球均值存在显著差异,其值从大到小依次为POLDER\,MISR和MODIS。3种产品的纬向平均在35°N以北区域表现出较大的差异。3种产品彼此之间相关性比较高,其中MODIS与MISR产品的相关性最强,MISR与POLDER产品的相关性最低,短波波段的相关系数(r) 分别为0.939与0.911。3种产品在可见光波段的相关性大于短波波段。在不同地表类型上,3种产品表现出了大致相似的差异,表明其对地表类型并不敏感。对气溶胶的分析表明:MODIS与MISR的550 nm气溶胶光学厚度(AOD)较为相似,其差异不足以解释DHR的差别;但是POLDER的865 nm AOD明显小于MISR,因此可以认为是由于POLDER的AOD估算偏低,导致了POLDER的DHR值大于MODIS与MISR。     

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).     

藏北那曲地区地面与MODIS反演的地表反照率对比分析

利用藏北高原那曲地区反照率地面观测资料分析了其日内、月均和季节变化特点,在此基础上,与同期的MODIS/Terra反演结果进行了对比分析。结果表明:藏北那曲地区晴天地表反照率日内变化明显,主要表现为早晚高,变幅大,中午低,变幅小的"U"形变化特点。早晨太阳高度角低,反照率高,日内随着太阳高度角的增加反照率逐渐降低,下午14:00~15:00反照率达到日内最小值,之后随着太阳高度角的降低,反照率上升明显。夏季日内最高反照率出现在早晨8:00,其他季节则基本上在下午18:00。太阳高度角同样是影响地表反照率季节性变化的主要原因,两者呈极显著的反相关关系,相关系数为-0.91。冬季平均地表反照率最高,为0.28,其次是春秋两季,均为0.23,夏季最低,为0.19。MODIS/Terra在上午12:00左右过境时反演的地表反照率与地面观测值之间存在很好的一致性,两者平均值都为0.22,相对误差9.60%,绝对误差和均方根误差(RMSE)均为0.02,而在13:00左右过境时,卫星反演值较观测值存在系统性偏小特点,平均偏小14.29%,相对误差为16.45%,绝对误差0.04,均方根误差0.05。此外,MODIS反演的地表反照率比地面观测值日际波动大。若不考虑积雪,冬夏两季地表反照率的空间差异小,而春秋两季空间差异较大,这主要与研究区地面植被类型及其季节性空间分布特征有关。     

Evaluation of MODIS water vapour products over China using radiosonde data

Radiosonde data collected from 83 stations in China from January to December 2012 were used to evaluate Moderate Resolution Imaging Spectroradiometer (MODIS) near-infrared (NIR) and thermal infrared (IR) total precipitable water vapour (PWV) products. The results indicate that MODIS NIR PWV products shows better agreement with radiosonde data than with IR PWV products, with the correlation coefficients up to 0.95. The root mean square errors (RMSEs) of NIR PWV range from 2 to 8 mm with different stations, which shows significant regional differences over China. The mean RMSE is about 5.03 mm (~35%) with a positive deviation of 2.56 mm (~18%), indicating the occurrence of a slight overestimation. Moreover, MODIS IR PWV during night-time has a better agreement with radiosonde PWV than that during daytime. The mean RMSE of IR PWV during daytime was ~6.02 mm (~42%), with a positive deviation of 1.54 mm (~11%). The mean RMSE of IR PWV during night-time was ~5.81 mm (~40%), with a negative deviation of approximately ?0.04 mm (~0.25%). Both the NIR and IR PWV products during daytime tend to be higher than radiosonde PWV.     

A simple correction method for the MODIS surface reflectance product over typical inland waters in China

The Moderate Resolution Imaging Spectroradiometer (MODIS) has the advantage of providing continuous, global, near-daily spatial measurements, and has greatly aided in understanding physical, optical, and biological processes in the global ocean biosphere. However, little research has been implemented for the remote-sensing monitoring of global inland waters. One important factor is that there is no operational atmospheric correction method designed for global inland waters. The MODIS surface reflectance product (MOD09) provides surface reflectance data for land at the global scale, but it does not offer accurate atmospheric correction over inland waters because of the constraints of its primary correction algorithm. The purpose of this article is to provide a simple and operational correction method for the MOD09 product to retrieve the water-leaving reflectance for large inland waters larger than 25 km². The correction method is based on an analysis of additive noises in MOD09 data over inland waters and on the adoption of two assumptions. Field-measured data collected in three typical inland waters in China were used to assess the performance of the correction method to ensure its applicability for waters in different conditions. The results show acceptable agreement with field data over the three inland waterbodies, with a mean relative error of 17.1% in visible bands. Our study demonstrates that the MOD09 correction method is moderately accurate when compared with the optimal method for specific waterbodies, but it has the potential for use in operational data-processing systems to derive water-leaving reflectance data from MOD09 data over inland waters in a variety of conditions and large regions.     

利用MODIS卫星资料分析北京地区地表反照率时空分布及变化特征

利用2000～ 2005 年MOD IS 地表双向反照率(MOD43B3) , 植被指数(MOD13A 2) 和地表覆盖类型(MOD12Q 1) 资料, 分析了北京及周边地区地表反照率的时空分布, 计算了2000～ 2005 年平均地表反照率的时间变化。结果表明, 北京城区年平均地表反照率为0. 12, 山区森林(0. 11) 明显小于平原地区(0. 15) , 而永定河流域反照率较大(0. 18) , 这主要是因为永定河流域植被覆盖度较低(植被指数低) , 山区地表反照率季节性变化不明显。2000～ 2005 年6 年的统计回归显示北京平原地区反照率呈略有下降。     

Using NASA's World Wind virtual globe for interactive internet visualization of the global MODIS burned area product

Three‐dimensional virtual globes are radically changing the way geographic information is perceived by the public. This article describes how NASA World Wind, an open source virtual globe, is currently being used for visualization of the MODIS burned area product. The procedures adopted for converting the product into a format compatible with World Wind, as well as the spatial generalization of these data at different scales, are described. Directions to instructions on how to obtain the MODIS burned area product visualization imagery and use it in World Wind are included. This article highlights the potential benefits of integrating the visualization capability of virtual globes into the next generation of remotely sensed product internet analysis and distribution systems.     

Development and evaluation of a cloud-gap-filled MODIS daily snow-cover product

The utility of the Moderate Resolution Imaging Spectroradiometer (MODIS) snow-cover products is limited by cloud cover which causes gaps in the daily snow-cover map products. We describe a cloud-gap-filled (CGF) daily snow-cover map using a simple algorithm to track cloud persistence, to account for the uncertainty created by the age of the snow observation. Developed from the 0.05° resolution climate-modeling grid daily snow-cover product, MOD10C1, each grid cell of the CGF map provides a cloud-persistence count (CPC) that tells whether the current or a prior day was used to make the snow decision. Percentage of grid cells “observable” is shown to increase dramatically when prior days are considered. The effectiveness of the CGF product is evaluated by conducting a suite of data assimilation experiments using the community Noah land surface model in the NASA Land Information System (LIS) framework. The Noah model forecasts of snow conditions, such as snow-water equivalent (SWE), are updated based on the observations of snow cover which are obtained either from the MOD10C1 standard product or the new CGF product. The assimilation integrations using the CGF maps provide domain-averaged bias improvement of ~11%, whereas such improvement using the standard MOD10C1 maps is ~3%. These improvements suggest that the Noah model underestimates SWE and snow depth fields, and that the assimilation integrations contribute to correcting this systematic error. We conclude that the gap-filling strategy is an effective approach for increasing cloud-free observations of snow cover.     

Prototyping of Land-SAF leaf area index algorithm with VEGETATION and MODIS data over Europe

The Satellite Application Facility on Land Surface Analysis (Land-SAF) aims to provide land surface variables for the meteorological and environmental science communities from EUMETSAT satellites. This study assesses the performance of a simplified (i.e. random distribution of vegetation is assumed) version of the Land-SAF algorithm for the estimation of Leaf Area Index (LAI) when prototyped with VEGETATION (processed in CYCLOPES program) and MODIS reflectances. The prototype estimates of LAI are evaluated both by comparison with validated CYCLOPES and MODIS LAI products derived from the same sensors and directly through comparison with ground-based estimates. Emphasis is given on evaluating the impact of the algorithm and input data on LAI retrieval discrepancies. Analysis is achieved over Europe for the 2000-2003 period. The results demonstrate the capacity of the Land-SAF algorithm to retrieve consistent LAI estimates from multiple optical sensors even when their reflectances present systematic differences. High spatial and temporal consistencies between Land-SAF prototype estimates and existing LAI products are found. The differences between Land-SAF and CYCLOPES LAI are lower than their uncertainties (RMSE (relative RMSE) within 0.4 (30%)). Land-SAF prototype estimates and MODIS LAI show larger discrepancies mainly due to differences in the vegetation structure representation and algorithm assumptions (RMSE ranging from 0.2 (30%) up to 0.8 (40%)). Land-SAF prototype provides higher LAI values than MODIS for herbaceous canopies (i.e. shrubs, grasses and crops) and lower values for woody biomes (i.e. savannas and forests). Direct validation indicates that LAI estimates from prototyping of the Land-SAF algorithm with CYCLOPES and MODIS reflectances achieve similar performances (differences with ground measurements are lower than 0.5 LAI units in 60% and 50% of the cases, respectively) as CYCLOPES and MODIS LAI products. Results from this prototyping exercise appear useful for improved retrieval of LAI and constitute a step forward for refinement, validation and consolidation of the Land-SAF algorithm.     

Validation of the MODIS reflectance product under UK conditions

Surface reflectance obtained from remote-sensing data is the main input to almost all remote-sensing applications. The availability and special characteristics of Moderate Resolution Imaging Spectroradiometer (MODIS) products have led to their use worldwide. Validation of the MODIS reflectance product is then crucial to provid information on uncertainty in the reflectance data, and in other MODIS products and in the applied surface–atmosphere models. Compact Airborne Spectrographic Imager (CASI) and Système Pour l'Observation de la Terre (SPOT) data, collected during the Network for Calibration and Validation in Earth Observation (NCAVEO) 2006 Field Campaign, were applied to validate daily MODIS reflectance data over a site in the southern UK. The difference in the view geometry of at-nadir CASI and SPOT data and off-nadir MODIS data was dealt with using a semi-empirical bidirectional reflectance distribution function (BRDF) model. The validation results showed that for our particular study site, the absolute errors in the MODIS reflectance product were too large to allow the albedo data to be used directly in climate models. The errors were mainly related to the uncertainties in the MODIS atmospheric variables, the BRDF model, and undetected clouds and cloud shadows. More generally, the study highlights the extreme difficulty of achieving pixel-level validation of coarse spatial resolution satellite sensor data in an environment in which the atmosphere is constantly changing, and in which the landscape is characterized by high space–time heterogeneity.     

Assessment of the MODIS LAI product for Australian ecosystems

The leaf area index (LAI) product from the Moderate Resolution Imaging Spectroradiometer (MODIS) is important for monitoring and modelling global change and terrestrial dynamics at many scales. The algorithm relies on spectral reflectances and a six biome land cover classification. Evaluation of the specific behaviour and performance of the product for regions of the globe such as Australia are needed to assist with product refinement and validation. We made an assessment of Collection 4 of the MODIS LAI product using four approaches: (a) assessment against a continental scale Structural Classification of Australian Vegetation (SCAV); (b) assessment against a continental scale land use classification (LUC); (c) assessment against historical field-based measurement of LAI collected prior to the Terra Mission; and (d) direct comparison of MODIS LAI with coincident field measurements of LAI, mostly from hemispherical photography. The MODIS LAI product produced a wide variety of geographically and structurally specific temporal response profiles between different classes and even for sub-groups within classes of the SCAV. Historical and concurrent field measurements indicated that MODIS LAI was giving reasonable estimates for LAI for most cover types and land use types, but that major overestimation of LAI occurs in some eastern Australian open forests and woodlands. The six biome structural land cover classification showed some significant deviations in class allocation compared to the SCAV particularly where grasslands are allocated to shrubland, savanna woodlands are allocated to shrubland, savanna and broadleaf forest, and open forests are allocated to savanna and broadleaf forest. The land cover and LAI products could benefit from some additional examination of Australian data addressing the structural representation of Eucalypt canopies in the “space of canopy realisation” for savanna and broadleaf forest classes.     

Evaluation of the MODIS LAI product using independent lidar-derived LAI: A case study in mixed conifer forest

This study presents an alternative assessment of the MODIS LAI product for a 58,000 ha evergreen needleleaf forest located in the western Rocky Mountain range in northern Idaho by using lidar data to model (R² = 0.86, RMSE = 0.76) and map LAI at higher resolution across a large number of MODIS pixels in their entirety. Moderate resolution (30 m) lidar-based LAI estimates were aggregated to the resolution of the 1-km MODIS LAI product and compared to temporally-coincident MODIS retrievals. Differences in the MODIS and lidar-derived values of LAI were grouped and analyzed by several different factors, including MODIS retrieval algorithm, sun/sensor geometry, and sub-pixel heterogeneity in both vegetation and terrain characteristics. Of particular interest is the disparity in the results when MODIS LAI was analyzed according to algorithm retrieval class. We observed relatively good agreement between lidar-derived and MODIS LAI values for pixels retrieved with the main RT algorithm without saturation for LAI LAI ≤ 4. Moreover, for the entire range of LAI values, considerable overestimation of LAI (relative to lidar-derived LAI) occurred when either the main RT with saturation or back-up algorithm retrievals were used to populate the composite product regardless of sub-pixel vegetation structural complexity or sun/sensor geometry. These results are significant because algorithm retrievals based on the main radiative transfer algorithm with or without saturation are characterized as suitable for validation and subsequent ecosystem modeling, yet the magnitude of difference appears to be specific to retrieval quality class and vegetation structural characteristics.