Bayesian maximum entropy data fusion of field-observed leaf area index (LAI) and Landsat Enhanced Thematic Mapper Plus-derived LAI

Accurate high-resolution leaf area index (LAI) reference maps are necessary for the validation of coarser-resolution satellite-derived LAI products. In this article, we propose an efficient method based on the Bayesian Maximum Entropy (BME) paradigm to combine field observations and Landsat Enhanced Thematic Mapper Plus (ETM+)-derived LAI surfaces in order to produce more accurate LAI reference maps. This method takes into account the uncertainties associated with field observations and with the regression relationship between ETM+-derived LAI and field measurements to perform a non-linear prediction of LAI, the variable of interest. In order to demonstrate the difference by soft data and hard data, we estimate the LAI reference maps by three BME interpolation methods, BME1, BME2, and BME3. BME1 and BME2 perform maximum estimation and mean estimation, respectively, by taking the ETM+-derived LAI as interval soft data and the field LAI measurements as hard data. BME3 is utilized when ETM+-derived LAI surfaces are processed as uniform probability soft data and field measurements are processed as Gaussian probability soft data. Three study sites are selected from the BigFoot project (NASA's Earth Observing System validation programme) (http://www.fsl.orst.edu/larse/bigfoot/index.html). In regard to the mean and standard deviation of LAI surfaces, standard deviation predicted by BME methods has lower values than that derived by ETM+. The mean value of the BME-predicted LAI, which takes into account the uncertainties of field measurements, is lower than that of ETM+-derived LAI at each study site. A comparison with field measurements shows that BME1, BME2, and BME3 have root mean square errors (RMSE) of 0.455, 0.485, and 0.517 and average biases of??0.017,??0.010, and??0.304, respectively. The RMSEs and biases of the predicted LAI surfaces are less when compared to the ETM+-derived LAI, which has the average RMSE and bias of 0.642 and??0.080. When the field measurements are processed as soft data, the predicted LAI by BME3 has more bias than those of the predictions by BME1 and BME2, but has less RMSE than that of the ETM+-derived LAI by 0.125. In summary, BME is capable of incorporating the spatial autocorrelation and the uncertainties in the field LAI measurements into the LAI surface estimation to produce a more accurate LAI surface with less RMSE in validation. The maximum estimation has relatively better accuracy than the mean estimation. The results indicate that the BME is a promising method for fusing point-scale and area-scale data.     

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

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

MODIS spectral signals at a flux tower site: Relationships with high-resolution data, and CO2 flux and light use efficiency measurements

In this study, spectral indices were calculated from single date HyMap (3 m; 126 bands), Hyperion (30 m; 242 bands), ASTER (15/30 m; 9 bands), and a time series of MODIS nadir BRDF-adjusted reflectance (NBAR; 1 km, 7 bands) for a study area surrounding the Tumbarumba flux tower site in eastern Australia. The study involved: a) the calculation of a range of physiologically-based vegetation indices from ASTER, HyMap, Hyperion and MOD43B NBAR imagery over the flux tower site; b) comparison across scales between HyMap, Hyperion and MODIS for the normalized difference water index (NDWI) and the Red-Green ratio; c) analysis of relationships between tower-based flux and light use efficiency (LUE) measurements and seasonal and climatic constraints on growth; and d) examination of relationships between fluxes, LUE and time series of NDVI, NDWI and Red-Green ratio. Strong seasonal patterns of variation were observed in NDWI and Red/Green ratio from MODIS NBAR. Correlations between fine (3 and 30 m) and coarse (1 km) scale indices for a small region around the flux tower site were moderately good for Red/Green ratio, but poor for NDWI. Hymap NDWI values for the understorey canopy were much lower than values for the tree canopy. MODIS NDWI was negatively correlated with CO₂ fluxes during warm and cool seasons. The correlation indicated that surface reflectance, affected by a spectrally bright grassland understorey canopy, was decoupled from growth of trees with access to deep soil moisture. The application of physiologically-based indices at earth observation scale requires careful attention to applicability of band configurations, contribution of vegetation components to reflectance signals, mechanistic relationships between biochemical processes and spectral indexes, and incorporation of ancillary information into any analysis.     

Improvements of the MODIS terrestrial gross and net primary production global data set

MODIS primary production products (MOD17) are the first regular, near-real-time data sets for repeated monitoring of vegetation primary production on vegetated land at 1-km resolution at an 8-day interval. But both the inconsistent spatial resolution between the gridded meteorological data and MODIS pixels, and the cloud-contaminated MODIS FPAR/LAI (MOD15A2) retrievals can introduce considerable errors to Collection4 primary production (denoted as C4 MOD17) results. Here, we aim to rectify these problems through reprocessing key inputs to MODIS primary vegetation productivity algorithm, resulting in improved Collection5 MOD17 (here denoted as C5 MOD17) estimates. This was accomplished by spatial interpolation of the coarse resolution meteorological data input and with temporal filling of cloud-contaminated MOD15A2 data. Furthermore, we modified the Biome Parameter Look-Up Table (BPLUT) based on recent synthesized NPP data and some observed GPP derived from some flux tower measurements to keep up with the improvements in upstream inputs. Because MOD17 is one of the down-stream MODIS land products, the performance of the algorithm can be largely influenced by the uncertainties from upstream inputs, such as land cover, FPAR/LAI, the meteorological data, and algorithm itself. MODIS GPP fits well with GPP derived from 12 flux towers over North America. Globally, the 3-year MOD17 NPP is comparable to the Ecosystem Model-Data Intercomparison (EMDI) NPP data set, and global total MODIS GPP and NPP are inversely related to the observed atmospheric CO₂ growth rates, and MEI index, indicating MOD17 are reliable products. From 2001 to 2003, mean global total GPP and NPP estimated by MODIS are 109.29 Pg C/year and 56.02 Pg C/year, respectively. Based on this research, the improved global MODIS primary production data set is now ready for monitoring ecological conditions, natural resources and environmental changes.     

基于MODIS和AMSR-E遥感数据的土壤水分降尺度研究

微波传感器获得的土壤水分产品空间分辨率一般都很粗,而流域尺度上的研究需要中高分辨率的土壤水分数据。用MODIS逐日地表温度产品MOD11A1和逐日地表反射率产品MOD09GA构建温度-植被指数特征空间,并计算得到TVDI(Temperature Vegetation Dryness Index)指数,它与土壤水分呈负相关关系,能够反映土壤水分的空间分布模式,但并不是真实的土壤水分值。在AMSR-E像元尺度上求得TVDI与土壤水分的负相关系数,进而对VUA AMSR-E土壤水分产品进行降尺度计算得到0.01°分辨率的真实土壤水分值。经NAFE06(The National Airborne Field Experiment 2006)试验地面采样数据验证,降尺度后的土壤水分均方根误差平均值为6.1%。     

Evaluation of MODIS Snow Products in Southwestern Xinjiang Using the Google Earth Engine

Google Earth Engine(GEE) is a cloud\|based geospatial processing platform that can analyze geospatial data to achieve parallel processing of massive remote sensing data on a global scale,providing support for remote sensing big data and large\|area research.MODIS snow cover mapping is a global snow cover product established using MODIS data and has been widely used in regional and global climate and environmental monitoring.In the GEE,millions of remote sensing images are stored,including MODIS daily snow products MOD10A1 V5 data and Landsat data.Taking the three research areas in southwestern Xinjiang as examples,the Landsat stored by the GEE were selected,and the NDSI was used to extract the snow cover as the true value of the land cover to evaluate the MOD10A1 accuracy.The results show that the average overall accuracy of MOD10A1 in the snow cover season in southwestern Xinjiang during the period from 2000 to 2016 is 82%,the average misjudgment rate is 2.9%,and the average missed rate is 58.8%.The overall accuracy of MOD10A1 can reach 98% under the clear sky conditions.The accuracy of MOD10A1 is effected by the terrain conditions and cloud cover in different regions.Therefore,the GEE can quickly and effectively filter high quality cloudless Landsat images,and evaluate the accuracy of the MOD10A1 in the snow area around the global regions,displaying intuitively the misjudgment and missed areas in the form of online maps.Meanwhile,GEE provides the Landsat simple cloud score function to calculate the regional cloud cover,which makes the influence of cloud cover on the MOD10A1 accuracy assessment more regionally representative.     

Improving continuity of MODIS terrestrial photosynthesis products using an interpolation scheme for cloudy pixels

The Moderate Imaging Spectroradiometer (MODIS) sensors onboard the NASA Terra and Aqua satellites provide the means for frequent measurement and monitoring of the status and seasonal variability in global vegetation phenology and productivity. However, while MODIS reflectance data are often interrupted by clouds, terrestrial processes like photosynthesis are continuous, so MODIS photosynthesis data must be able to cope with cloudy pixels. We developed cloud‐correction algorithms to improve retrievals of the MODIS photosynthesis product (PSNnet) corresponding to clear sky conditions by proposing four alternative cloud‐correction algorithms, which have different levels of complexity and correct errors associated with cloudy‐pixel surface reflectance. The cloud‐correction algorithms were applied at four weather stations, two fluxtower sites and the Pacific Northwest (PNW) region of the USA to test a range of cloud climatologies. Application of the cloud‐correction algorithms increased the magnitude of both daily and annual MODIS PSNnet results. Our results indicate that the proposed cloud correction methods improve the current MODIS PSNnet product considerably at both site and regional scales and weekly to annual time steps for areas subjected to frequent cloud cover. The corrections can be applied as a post‐processing interpolation of PSNnet, and do not require reprocessing of the MOD17A2 algorithm.     

TM与MODIS植被供水指数反演及对比分析

为了探讨多源卫星遥感数据针对同一地区同一时刻旱情所得植被供水指数的差异,选择我国华北地区旱灾发生频率较高、影响较广的河北地区为研究区,针对近年来干旱监测应用较为广泛的Landsat TM/ETM+和EOS MODIS数据,分别进行植被供水指数的提取,并进行两者之间的对比分析,得出以下结论:① TM VSWI (Vegetation Supply Water Index) 与MODIS VSWI之间数值上存在一定的差别,变化范围在-0.51~0.20之间。其中负值主要集中在城镇、裸地及水体地表;② 在植被覆盖区,TM的平均VSWI大于MODIS的,但两者差别不大;③ 在各种植被覆盖度条件下,TM VSWI与MODIS VSWI差值的最小值、最大值和均值均表现出随着植被覆盖度的增加,其值逐渐增大的特点。该结论可以为两种遥感数据源干旱监测的差异分析及综合应用提供重要的参考依据。     

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

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

MODIS time-series imagery for forest disturbance detection and quantification of patch size effects

The Moderate Resolution Imaging Spectroradiometer (MODIS) 250 m single day surface reflectance (MOD09GQK) and 16-day composite gridded vegetation index data (MOD13Q1) were used to detect forest harvest disturbance between 2000 and 2004 in northern Maine. A MODIS multi-date Normalized Difference Vegetation Index (NDVI) forest change detection map was developed from each MODIS data set. A Landsat TM/ETM+ change detection map was developed as a reference to assess the effect of disturbed forest patch size on classification accuracy (agreement) and disturbed area estimates of MODIS. The MODIS single day and 16-day composite data showed no significant difference in overall classification accuracies. However, the 16-day NDVI change detection map had marginally higher overall classification accuracy (at 85%), but had significantly lower detection accuracy related to disturbed patch size than the single day NDVI change detection map. The 16-day composite NDVI data achieved 69% detection accuracy and the single day NDVI achieved 76% when the disturbed patch size was greater than 20 ha. The detection accuracy increased to approximately 90% for both data sets when the patch size exceeded 50 ha. The R² (range 0.6 to 0.9) and slope (range 0.5 to 0.9) of regression lines between Landsat and MODIS data (based on forest disturbance percent of township) increased with the mean disturbed patch size of each township. The 95% confidence intervals of forest disturbance percent estimate for each township were narrow with less than 1% of each township at the mean MODIS forest disturbance level.     

Regionally adaptable dNBR-based algorithm for burned area mapping from MODIS data

Recent advances in instrument design have led to considerable improvements in wildfire mapping at regional and global scales. Global and regional active fire and burned area products are currently available from various satellite sensors. While only global products can provide consistent assessments of fire activity at the global, hemispherical or continental scales, the efficiency of their performance differs in various ecosystems. The available regional products are hard-coded to the specifics of a given ecosystem (e.g. boreal forest) and their mapping accuracy drops dramatically outside the intended area. We present a regionally adaptable semi-automated approach to mapping burned area using Moderate Resolution Imaging Spectroradiometer (MODIS) data. This is a flexible remote sensing/GIS-based algorithm which allows for easy modification of algorithm parameterization to adapt it to the regional specifics of fire occurrence in the biome or region of interest. The algorithm is based on Normalized Burned Ratio differencing (dNBR) and therefore retains the variability of spectral response of the area affected by fire and has the potential to be used beyond binary burned/unburned mapping for the first-order characterization of fire impacts from remotely sensed data. The algorithm inputs the MODIS Surface Reflectance 8-Day Composite product (MOD09A1) and the MODIS Active Fire product (MOD14) and outputs yearly maps of burned area with dNBR values and beginning and ending dates of mapping as the attributive information. Comparison of this product with high resolution burn scar information from Landsat ETM+ imagery and fire perimeter data shows high levels of accuracy in reporting burned area across different ecosystems. We evaluated algorithm performance in boreal forests of Central Siberia, Mediterranean-type ecosystems of California, and sagebrush steppe of the Great Basin region of the US. In each ecosystem the MODIS burned area estimates were within 15% of the estimates produced by the high resolution base with the R² between 0.87 and 0.99. In addition, the spatial accuracy of large burn scars in the boreal forests of Central Siberia was also high with Kappa values ranging between 0.76 and 0.79.     

On the upstream inputs into the MODIS primary productivity products using biometric data from oil palm plantations

This study evaluated the influence of upstream inputs into the Moderate Resolution Imaging Spectroradiometer (MODIS) primary productivity products, termed the MOD17, at tropical oil palm plantations (Elaeis guineensis Jacq.). Evaluation of MOD17 using oil palm plantations as test sites is ideal because the plantations are cultivated on large areas which are comparable with the size of MODIS pixels. It is difficult to find test sites covered by other single species in a whole pixel. The upstream inputs studied included (1) MODIS land cover, (2) the National Centers for Environmental Prediction–Department of Energy (NCEP-DOE) Reanalysis 2 meteorological data set, (3) MODIS leaf area index/fraction of photosynthetically active radiation (LAI/fPAR), and (4) MODIS maximum light-use efficiency (maximum LUE). Oil palm biometric and local meteorological data were utilized as ground data. Furthermore, scaling up oil palm LAI and fPAR from plot scale to regional scale (Peninsular Malaysia) was done empirically by correlating oil palm LAI derived from the hemispherical photography technique with radiance information from the Disaster Monitoring Constellation 2 satellite (UK-DMC 2). The upscaled LAI/fPAR developed in this study was used to evaluate the MODIS LAI/fPAR. The results showed that the MODIS land-cover product has an overall accuracy of 78.8% when compared to the Peninsular Malaysia land-use map produced by the Department of Agriculture, Malaysia. Regarding the NCEP-DOE Reanalysis 2 data set, vapour pressure deficit (VPD) and photosynthetically active radiation (PAR) contain large uncertainties in our study area. However, MODIS LAI and fPAR were correlated relatively well with the upscaled LAI (R² = 0.50) and the upscaled fPAR (R² = 0.60), respectively. The constant values of maximum LUE for croplands and evergreen broadleaf forest ecosystems are lower than the maximum LUE of oil palm. The relative predictive error assessment showed that the MOD17 net primary productivity (NPP) overestimated oil palm NPP derived from biometric methods by 142–204%. We replaced the upstream inputs of MOD17 by the local inputs for estimating oil palm GPP and NPP in Peninsular Malaysia. This was done by (1) assigning maximum LUE for oil palm plantations as a constant at 1.68 g C m^?2 day^?1, (2) utilizing meteorological data from local meteorological stations, and (3) using the upscaled fPAR of oil palm plantations. The amount of oil palm GPP and NPP for Peninsular Malaysia in 2010 were estimated to be ~0.09 Pg C year^?1 (or equivalent to ~0.33 Pg CO₂ year^?1) and ~0.03 Pg C year^?1 (~0.11 Pg CO₂ year^?1), respectively, indicating that oil palm plantations in Peninsular Malaysia can play an important role in global carbon sequestration. In the future there is likely to be a demand for MODIS GPP and NPP products that are more accurate than those currently generated by MOD17. We recommend future developments of the MOD17 processing system to allow improvements in the upstream input parameters, in the manner described in this article, both for global processing and for the production of more accurate values for GPP and NPP at regional and local scales.     

Validation of collection of 6 MODIS/Terra and MODIS/Aqua gross primary production in an alpine meadow of the Northern Tibetan Plateau

Moderate Resolution Imaging Spectroradiometer (MODIS) continuously monitors gross primary production (GPP), which is an extremely important component of carbon cycling, at the global scale. Uncertainties about MODIS GPP limit our ability to accurately quantify GPP at the regional scales. The Collection 6 MODIS/Terra and MODIS/Aqua GPP products (i.e. MOD17A2H and MYD17A2H) were compared with the estimated GPP (GPP_EC) by eddy covariance measurements in an alpine meadow in the Northern Tibetan Plateau during three consecutive growing seasons of 2005–2007. The Collection 6 MODIS/Terra and MODIS/Aqua fractional photosynthetically active radiation (FPAR) products (i.e. MOD15A2H and MYD15A2H) were also validated. The MOD17A2H and MYD17A2H products tended to overestimate GPP_EC by 2.17% and 7.35% in 2005–2007, respectively, although these differences were not significant. The MOD15A2H and MYD15A2H products also tended to overestimate ground-based FPAR (FPAR_G) by 20.31% and 24.73% in 2005–2007, respectively. The overestimation of FPAR resulted in about 17.51–23.97% overestimation of GPP_EC. The default maximum light-use efficiency (ε_max) of 0.86 g C MJ^?1 only underestimated the ground-based ε_max (0.88 g C MJ^?1) by 2.27%, which in turn resulted in about 2.13–2.72% underestimation of GPP_EC. The meteorology data errors only caused about 0.48–1.06% underestimation of GPP_EC. Therefore, although MODIS Collection 6 GPP had a very high accuracy, the input parameters had relative greater errors in the alpine meadow of the Northern Tibetan Plateau. The differences between MODIS GPP and GPP_EC mainly resulted from FPAR, followed by ε_max and meteorological data.     

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.     

Validation of the MODIS snow product and cloud mask using student and NWS cooperative station observations in the Lower Great Lakes Region

NASA's Moderate Resolution Imaging Spectroradiometer (MODIS) snow product (MOD10) creates automated daily, 8-day composite and monthly regional and global snow cover maps. In this study, the MOD10 daily swath imagery (MOD10_L2) and the MODIS cloud mask (MOD35) were validated in the Lower Great Lakes Region, specifically the area to the east of Lake Michigan. Validation of the MOD10_L2 snow product, MOD35 cloud mask and the MOD10_L2 Liberal Cloud Mask was performed using field observations from K-12 student GLOBE (Global Learning and Observations to Benefit the Environment) and SATELLITES (Students And Teachers Evaluating Local Landscapes to Interpret The Earth from Space) programs. Student data consisted of field observations of snow depth, snow water equivalency, cloud type, and total cloud cover. In addition, observations from the National Weather Service (NWS) Cooperative Observing Stations were used. Student observations were taken during field campaigns in the winter of 2001-2002, a winter with very little snow in the Great Lakes region, and the winters of 2000-2001 and 2002-2003, which had significant snow cover. Validation of the MOD10_L2 version 4 snow product with student observations produced an accuracy of 92% while comparison with the NWS stations produced an accuracy of 86%. The higher NWS error appears to come from forested areas. Twenty-five and fifty percent of the errors observed by the students and NWS stations, respectively, occurred when there was only a trace of snow. In addition, 82% of the MODIS cloud masked pixels were identified as either overcast or broken by the student observers while 74% of the pixels the MODIS cloud mask identified as cloudless were identified as clear, isolated or scattered cloud cover by the student observers. The experimental Liberal Cloud Mask eliminated some common errors associated with the MOD35 cloud mask, however, it was found to omit significant cloud cover.     

A new global tree-cover percentage map using MODIS data

Global tree-cover percentage is an important parameter used to understand the global environment. However, the available global percentage tree-cover products are few. Producing a new global-scale data set facilitates comparison analysis among maps. Our study was undertaken to map tree-cover percentage on a global scale with better accuracy than previous studies. In this study, we estimated the tree-cover percentage on a global scale at a pixel size of 500 m using a modified supervised regression tree algorithm from the Moderate Resolution Imaging Spectroradiometer (MODIS) data of 2008. Training data were derived from high-resolution images displayed in Google Earth and created by linear mixture simulation. The estimation model was modified to fit the reference data, which were randomly collected from the global area. The estimation result was validated at 1106 random sample points. The root mean square error (RMSE) between estimated and observed tree cover was 13.8%. The produced map was also compared with existing global data sets. The RMSE of our result was better than that of two existing global percentage tree-cover data sets (about 1.8% and 9.2% lower than Vegetation Continuous Fields MOD44B data set and Global Map – Percent Tree Cover data set, respectively). In our result and all other data sets, the accuracy was lower for forests than for other areas. When the produced map was compared with the Vegetation Continuous Fields MOD44B data set, the low agreement between them was concentrated in closed forests of Colombia, India, and Indonesia, and closed to open forests of Central Africa.     

Assessment of MODIS LAI retrievals over soybean crop in Southern Brazil

In the present paper we have looked into the excessive occurrence of 255 standard fill value retrievals in Collection 4 MODIS LAI product over soybean areas from crop year 2001/2002 to 2004/2005, in Southern Brazil. The 255 standard fill value indicates that no leaf area index (LAI) retrieval was possible for the considered pixel. Time series of eight‐day composite LAI images (MOD15A2) and 16‐day composite NDVI images (MOD13Q1) were both compared with a soybean reference map derived from multitemporal Landsat images. The Land Cover Type 3 product (MOD12Q1) was also analysed to verify if the occurrence of those retrievals was related to misclassification of the broadleaf crops biome. Results indicated that the 255 standard fill value retrievals in Collection 4 LAI product were mainly related to soybean areas during peak growing season and occurred in every crop year we have studied. Eventual misclassification in the biome map was not the cause of those retrievals in the Collection 4 MODIS LAI product.     

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.     

Multi-index-based soil moisture estimation using MODIS images

The suitability of using Moderate Resolution Imaging Spectroradiometer (MODIS) images for surface soil moisture estimation to investigate the importance of soil moisture in different applications, such as agriculture, hydrology, meteorology and natural disaster management, is evaluated in this study. Soil moisture field measurements and MODIS images of relevant dates have been acquired. Normalized Difference Vegetation Index (NDVI), Enhanced Vegetation Index (EVI) and Normalized Difference Water Index (NDWI) are calculated from MODIS images. In addition, MODIS Land Surface Temperature (LST) data (MOD11A1) are used in this analysis. Four different soil moisture estimation models, which are based on NDVI–LST, EVI–LST, NDVI–LST–NDWI and EVI–LST–NDWI, are developed and their accuracies are assessed. Statistical analysis shows that replacing EVI with NDVI in the model that is based on LST and NDVI increases the accuracy of soil moisture estimation. Accuracy evaluation of soil moisture estimation using check points shows that the model based on LST, EVI and NDWI values gives a higher accuracy than that based on LST and EVI values. It is concluded that the model based on the three indices is a suitable model to estimate soil moisture through MODIS imagery.     

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.