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.     

Validation of MODIS-GPP product at 10 flux sites in northern China

Gross primary production (GPP) is an important variable in studies of the carbon cycle and climate change. The Moderate Resolution Imaging Spectroradiometer (MODIS)-GPP product (MOD17) provides global GPP data for terrestrial ecosystems; however, it is not well validated in China. In this study, an eddy covariance (EC) system observed GPP at 10 sites in northern China and was used to validate MOD17. The results indicated that MOD17 presents a strong bias in the study region due to the meteorological data, MODIS FPAR (fraction of absorbed photosynthetically active radiation) (MOD15), and the model parameters in the MODIS-GPP algorithm, Biome Parameters Look Up Table (BPLUT). Maximum light-use efficiency (?₀) had the strongest impact on the predicted GPP of the MODIS-GPP algorithm. After using the inputs observed in situ and improving parameters in the MODIS-GPP algorithm, the model could explain 85% of the EC-observed GPP of the sites, whereas the MODIS-GPP algorithm without in situ inputs and parameters only explained 26% of EC-observed GPP.     

Scaling Gross Primary Production (GPP) over boreal and deciduous forest landscapes in support of MODIS GPP product validation

The Moderate Resolution Imaging Radiometer (MODIS) is the primary instrument in the NASA Earth Observing System for monitoring the seasonality of global terrestrial vegetation. Estimates of 8-day mean daily gross primary production (GPP) at the 1 km spatial resolution are now operationally produced by the MODIS Land Science Team for the global terrestrial surface using a production efficiency approach. In this study, the 2001 MODIS GPP product was compared with scaled GPP estimates (25 km²) based on ground measurements at two forested sites. The ground-based GPP scaling approach relied on a carbon cycle process model run in a spatially distributed mode. Land cover classification and maximum annual leaf area index, as derived from Landsat ETM+ imagery, were used in model initiation. The model was driven by daily meteorological observations from an eddy covariance flux tower situated at the center of each site. Model simulated GPPs were corroborated with daily GPP estimates from the flux tower. At the hardwood forest site, the MODIS GPP phenology started earlier than was indicated by the scaled GPP, and the summertime GPP from MODIS was generally lower than the scaled GPP values. The fall-off in production at the end of the growing season was similar to the validation data. At the boreal forest site, the GPP phenologies generally agreed because both responded to the strong signal associated with minimum temperature. The midsummer MODIS GPP there was generally higher than the ground-based GPP. The differences between the MODIS GPP products and the ground-based GPPs were driven by differences in the timing of FPAR and the magnitude of light use efficiency as well as by differences in other inputs to the MODIS GPP algorithm—daily incident PAR, minimum temperature, and vapor pressure deficit. Ground-based scaling of GPP has the potential to improve the parameterization of light use efficiency in satellite-based GPP monitoring algorithms.     

Evaluation of Collections 4 and 5 of the MODIS Gross Primary Productivity product and algorithm improvement at a tropical savanna site in northern Australia

In this study, we assessed the accuracy of the MODIS (Moderate Resolution Imaging Spectroradiometer) GPP (gross primary productivity) Collections 4.5, 4.8 and 5 along with Leaf Area Index (LAI), fraction of absorbed Photosynthetically Active Radiation (fPAR), light use efficiency (LUE) and meteorological variables that are used to estimate GPP for a northern Australian savanna site. Results of this study indicated that the MODIS products captured the seasonal variation in GPP, LAI and fPAR well. Using the index of agreement (IOA), it was found that Collections 4.5 and 4.8 (IOA 0.89 respectively) agreed reasonably well with flux tower measurements between 2001 and 2006. It was also found that MODIS Collection 4.5 predicted the dry season GPP well (Relative Predictive Error (RPE) 4.17%, IOA 0.72 and Root Mean Square Error (RMSE) of 1.05 g C m^− 2 day^− 1), whilst Collection 4.8 performed better in capturing wet season dynamics (RPE 1.11%, IOA 0.80 and RMSE of 0.91 g C m^− 2 day^− 1). Although the wet season magnitude of GPP was predicted well by Collection 4.8, an examination of the inputs to the GPP algorithm revealed that MODIS fPAR was too high, but this was compensated by PAR and LUE that was too low. Although LAI and fPAR estimated by Collection 5 were more accurate, GPP for this Collection resulted in a much lower value (RPE 25%) due to errors in other factors. Recalculation of MODIS GPP using site specific input parameters indicated that MODIS fPAR was the main reason for the differences between MODIS and tower derived GPP followed by LUE and meteorological inputs. GPP calculated using all site specific values agreed very well with tower data on an annual basis (IOA 0.94, RPE 6.06% and RMSE 0.83 g C m^− 2 day^− 1) but the early initiation of the growing season calculated by the MODIS algorithm was improved when the vapor pressure deficit (VPD) function was replaced with a soil water deficit function. The results of this study however, reinforce previous findings in water limited regions, like Australia, and incorporation of soil moisture in a LUE model is needed to accurately estimate the productivity.     

Evaluation of satellite based primary production modelling in the semi-arid Sahel

The MODIS (Moderate Resolution Imaging Spectroradiometer) primary productivity products are evaluated against observed Above-ground Net Primary Production (AGNPP) in the semi-arid Senegal 2001. MODIS net primary productivity (NPP) modelling is a light use efficiency (LUE) based approach incorporating constraints on vegetation productivity arising from simulated radiation, water demand and temperature data from NASA's Data Assimilation Office (DAO). Annually integrated MODIS PSN (MOD17A2 net photosynthesis, Collection 4) explains more of the observed biomass variation (r² = 0.77) than MODIS fAPAR (fraction Absorbed Photosynthetically Active Radiation, Collection 4) (r² = 0.72), indicating the effect of including the canopy stress scalar (ε_s) based on DAO data combined with modelled maintenance respiration costs (of leaf and fine roots). Annual MODIS NPP (MOD17A3, Collection 4 (C4) and Collection 4.5 (C4.5)) including growth respiration and live wood maintenance respiration costs and modified DAO input (C4.5) however increases the residual unexplained observed AGNPP variance (C4 NPP; r² = 0.49) (C4.5 NPP; r² = 0.37). The overall quality of the annual NPP MODIS C4 and C4.5 products are moderate for the semi-arid Senegal because of the annual respiration cost modelling and a change in C4.5 biome-specific parameters stored in a Biome Properties Look-Up Table (BPLUT) is the main contributor to the observed discrepancy between C4 and C4.5 NPP. The dynamic range of the values of all MOD17 products was too low when compared to observed AGNPP. An estimate of canopy water stress (SIWSI; Shortwave Infrared Water Stress Index) derived from MODIS channels 2 and 6 and photosynthetically active radiation (PAR) irradiance derived from geostationary METEOSAT data were tested for primary production modelling using a stepwise linear regression analysis. PAR irradiance was combined with MODIS fAPAR into APAR (Absorbed Photosynthetically Active Radiation) explaining 79% of the observed AGNPP variation. Introducing SIWSI significantly increased the explained variance of observed AGNPP (r² = 0.89). MODIS-derived percentage tree cover was tested as a predictor based on the hypothesis that tree cover provides information on differences in respiratory costs between trees and grasses thereby accounting for variations in the LUE conversion efficiency ε. No significant reduction in residual unexplained AGNPP variance was found. Earth observation based derivation of PAR and canopy water stress from SIWSI suggest potential improvements to primary production models in semi-arid biomes that can be implemented in general NPP modelling LUE methodology.     

Estimation of net primary productivity of forests by modified CASA models and remotely sensed data

To increase the accuracy of predicting net primary productivity (NPP), in this study, Carnegie–Ames–Stanford Approach (CASA) model was modified by developing new methods to estimate absorbed photosynthetically active radiation or fraction of photosynthetically active radiation (FPAR) and water stress coefficient (WSC). In the modified model, FPAR was derived based on its non-linear relationship with leaf area index. Moreover, WSC was estimated using leaf water potential from soil moisture instead of a traditional evapotranspiration-based method. This study was conducted in Baiyun District area of Guangzhou, China, using Gaofen-1 (GF-1), Landsat 7, and Moderate Resolution Imaging Spectroradiometer (MODIS) satellite images. The predictions from the original and three modified CASA models and MODIS NPP product MOD17A3 were compared with field observations. The results showed that all the CASA-based models led to similar spatial distributions of forest aboveground NPP estimates. Overall, the estimates increased with elevation because the valley bottoms were dominated by developed or urbanized areas whereas the hillslopes and hilltops were largely vegetated. Based on root mean square error (RMSE) and relative RMSE between the observed and predicted values, the CASA model that integrated the modifications of both FPAR and WSC increased the estimation accuracy of NPP by 8.1% over the original one. The increase in accuracy was mainly contributed by the modification of FPAR. This suggested that the modification of FPAR provided greater potential than that of WSC for improving the predictions of CASA model. Compared to the CASA models, MOD17A3 had lower accuracy of aboveground NPP estimates. This study also showed that the fine spatial resolution GF-1 image provided a new source of data used to estimate NPP of forest ecosystems.     

Towards the development of a regional version of MOD17 for the determination of gross and net primary productivity of oil palm trees

Conducting quantitative studies on the carbon balance or productivity of oil palm is important for understanding the role of this ecosystem in global climate change. The MOD17 algorithm is used for processing data from the Moderate Resolution Imaging Spectroradiometer (MODIS) to generate the values of gross primary productivity (GPP) and net primary productivity for input to global carbon cycle modelling. In view of the increasing importance of data on carbon sequestration at regional and national levels, we have studied one important factor affecting the accuracy of the implementation of MOD17 at the sub-global level, namely the database of MODIS land cover (MOD12Q1) used by MOD17. By using a study area of approximately 7 km × 7 km (49 MODIS pixels) in semi-rural Johor in Peninsular Malaysia and using Google Earth 0.75 m resolution images as ground data, we found that the land-cover type for only 16 of these 49 MODIS pixels was correctly identified by MOD12Q1 using its 1 km resolution land-cover database. This leads to errors of 24% to 50% in the maximum light use efficiency, leading to corresponding errors of 24% to 50% in the GPP. We show that by using the Finer Resolution Observation and Monitoring – Global Land Cover (FROM-GLC) land-cover database developed by Gong et al., this particular error can be essentially eliminated, but at the cost of using extra computing resources.     

基于MODIS数据的玉米植被参数估算方法的对比分析

基于实测数据建立了FPAR、LAI的植被指数估算模型（NDVI、RVI、NDWI）,并将其应用于MODIS BRDF数据对德惠地区玉米FPAR、LAI进行估算,然后将MODIS 15A2 FPAR/LAI产品值分别与BRDF估算值、地面实测值进行对比分析。主要得出以下结论：植被指数NDVI、RVI都能较好地用于实测数据和MODIS BRDF数据的FPAR、LAI估算;NDWI虽然在实测数据中估算玉米FPAR、LAI的效果优于NDVI、RVI,但其应用于MODIS BRDF数据估算FPAR、LAI时,效果却较差。BRDF数据估算FPAR与MODIS 15A2 FPAR值的关系因生长时期不同而异,在玉米生长前期,前者高于后者,而生长后期两者却较相近;BRDF估算LAI值一直都高于MODIS 15A2 LAI产品值。生长季前期,MOD15A2 FPAR、LAI值接近实测值,而在后期却高于实测值。通过分析也表明,玉米苗期MODIS 15A2 FPAR数值变化范围较小,产品算法对实际FPAR变化尚不够敏感,这可能是影响MODIS FPAR产品精度的一个原因。     

Using LiDAR and quickbird data to model plant production and quantify uncertainties associated with wetland detection and land cover generalizations

Spatiotemporal data from satellite remote sensing and surface meteorology networks have made it possible to continuously monitor global plant production, and to identify global trends associated with land cover/use and climate change. Gross primary production (GPP) and net primary production (NPP) are routinely derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) onboard satellites Terra and Aqua, and estimates generally agree with independent measurements at validation sites across the globe. However, the accuracy of GPP and NPP estimates in some regions may be limited by the quality of model input variables and heterogeneity at fine spatial scales. We developed new methods for deriving model inputs (i.e., land cover, leaf area, and photosynthetically active radiation absorbed by plant canopies) from airborne laser altimetry (LiDAR) and Quickbird multispectral data at resolutions ranging from about 30 m to 1 km. In addition, LiDAR-derived biomass was used as a means for computing carbon-use efficiency. Spatial variables were used with temporal data from ground-based monitoring stations to compute a six-year GPP and NPP time series for a 3600 ha study site in the Great Lakes region of North America. Model results compared favorably with independent observations from a 400 m flux tower and a process-based ecosystem model (BIOME-BGC), but only after removing vapor pressure deficit as a constraint on photosynthesis from the MODIS global algorithm. Fine-resolution inputs captured more of the spatial variability, but estimates were similar to coarse-resolution data when integrated across the entire landscape. Failure to account for wetlands had little impact on landscape-scale estimates, because vegetation structure, composition, and conversion efficiencies were similar to upland plant communities. Plant productivity estimates were noticeably improved using LiDAR-derived variables, while uncertainties associated with land cover generalizations and wetlands in this largely forested landscape were considered less important.     

Reprocessing the MODIS Leaf Area Index products for land surface and climate modelling

Land surface and climate modelling requires continuous and consistent Leaf Area Index (LAI). High spatiotemporal resolution and long-time record data are more in demand nowadays and will continue to be in the future. MODIS LAI products meet these requirements to some degree. However, due to the presence of cloud and seasonal snow cover, the instrument problems and the uncertainties of retrieval algorithm, the current MODIS LAI products are spatially and temporally discontinuous and inconsistent, which limits their application in land surface and climate modelling. To improve the MODIS LAI products on a global scale, we considered the characteristics of the MODIS LAI data and made the best use of quality control (QC) information, and developed an integrated two-step method to derive the improved MODIS LAI products effectively and efficiently on a global scale. First, we used the modified temporal spatial filter (mTSF) method taking advantage of background values and QC information at each pixel to do a simple data assimilation for relatively low quality data. Then we applied the post processing-TIMESAT (A software package to analyze time-series of satellite sensor data) Savitzky-Golay (SG) filter to get the final result. We implemented the method to 10 years of the MODIS Collection 5 LAI data. In comparison with the LAI reference maps and the MODIS LAI data, our results showed that the improved MODIS LAI data are closer to the LAI reference maps in magnitude and also more continuous and consistent in both time-series and spatial domains. In addition, simple statistics were used to evaluate the differences between the MODIS LAI and the improved MODIS LAI.     

Evaluating evapotranspiration and water-use efficiency of terrestrial ecosystems in the conterminous United States using MODIS and AmeriFlux data

In this study, we used the remotely-sensed data from the Moderate Resolution Imaging Spectrometer (MODIS), meteorological and eddy flux data and an artificial neural networks (ANNs) technique to develop a daily evapotranspiration (ET) product for the period of 2004-2005 for the conterminous U.S. We then estimated and analyzed the regional water-use efficiency (WUE) based on the developed ET and MODIS gross primary production (GPP) for the region. We first trained the ANNs to predict evapotranspiration fraction (EF) based on the data at 28 AmeriFlux sites between 2003 and 2005. Five remotely-sensed variables including land surface temperature (LST), normalized difference vegetation index (NDVI), normalized difference water index (NDWI), leaf area index (LAI) and photosynthetically active radiation (PAR) and ground-measured air temperature and wind velocity were used. The daily ET was calculated by multiplying net radiation flux derived from remote sensing products with EF. We then evaluated the model performance by comparing modeled ET with the data at 24 AmeriFlux sites in 2006. We found that the ANNs predicted daily ET well (R² = 0.52-0.86). The ANNs were applied to predict the spatial and temporal distributions of daily ET for the conterminous U.S. in 2004 and 2005. The ecosystem WUE for the conterminous U.S. from 2004 to 2005 was calculated using MODIS GPP products (MOD17) and the estimated ET. We found that all ecosystems' WUE-drought relationships showed a two-stage pattern. Specifically, WUE increased when the intensity of drought was moderate; WUE tended to decrease under severe drought. These findings are consistent with the observations that WUE does not monotonously increase in response to water stress. Our study suggests a new water-use efficiency mechanism should be considered in ecosystem modeling. In addition, this study provides a high spatial and temporal resolution ET dataset, an important product for climate change and hydrological cycling studies for the MODIS era.     

Comparison and validation of MODIS and VEGETATION global LAI products over four BigFoot sites in North America

A new set of recently developed leaf area index (LAI) algorithms has been employed for producing a global LAI dataset at 1 km resolution and in time-steps of 10 days, using data from the Satellite pour l'observation de la terre (SPOT) VEGETATION (VGT) sensor. In this paper, this new LAI product is compared with the global MODIS Collection 4 LAI product over four validation sites in North America. The accuracy of both LAI products is assessed against seven high resolution ETM+ LAI maps derived from field measurements in 2000, 2001, and 2003. Both products were closely matched outside growing season. The MODIS product tended to be more variable than the VGT product during the summer period when the LAI was maximum. VGT and ETM+ LAI maps agreed well at three out of the four sites. The median relative absolute error of the VGT LAI product varied from 24% to 75% at 1 km scale and it ranged from 34% to 88% for the MODIS LAI product. The importance of correcting field measurements for the clumping effect is illustrated at the deciduous broadleaf forest site (HARV). Inclusion of the sub-pixel land cover information improved the quality of LAI estimates for the prairie grassland KONZ site. Further improvement of the global VGT LAI product is suggested by production and inclusion of pixel-specific global foliage clumping index and forest background reflectance maps that would serve as an input into the VGT LAI algorithms.     

Evaluation of MODIS gross primary productivity and land cover products for the humid tropics using oil palm trees in Peninsular Malaysia and Google Earth imagery

Conducting quantitative studies on the carbon balance or productivity of oil palm is important in understanding the role of this ecosystem in global climate change. In this study, we evaluated the accuracy of MODIS (Moderate Resolution Imaging Spectroradiometer) annual gross primary productivity (GPP) (the product termed MOD-17) and its upstream products, especially the MODIS land cover product (the product termed MOD-12). We used high-resolution Google Earth images to classify the land cover classes and their percentage cover within each 1 km spatial resolution MODIS pixel. We used field-based annual GPP for 2006 to estimate GPP for each pixel based on percentage cover. Both land cover and GPP were then compared to MODIS land cover and GPP products. The results show that for pure pixels that are 100% covered by mature oil palm trees, the RMSE (root mean square error) between MODIS and field-based annual GPP is 18%, and that this is increased to 27% for pixels containing mostly oil palm. Overall, for an area of about 42 km² the RMSE is 26%. We conclude that land cover classification (at 1 km resolution) is one of the main factors for the discrepancy between MODIS and field-based GPP. We also conclude that the accuracy of the MODIS GPP product could be improved significantly by using higher-resolution land cover maps.     

Comparisons of land cover and LAI estimates derived from ETM+ and MODIS for four sites in North America: a quality assessment of 2000/2001 provisional MODIS products

The MODIS land science team produces a number of standard products, including land cover and leaf area index (LAI). Critical to the success of MODIS and other sensor products is an independent evaluation of product quality. In that context, we describe a study using field data and Landsat ETM+ to map land cover and LAI at four 49-km² sites in North America containing agricultural cropland (AGRO), prairie grassland (KONZ), boreal needleleaf forest, and temperate mixed forest. The purpose was to: (1) develop accurate maps of land cover, based on the MODIS IGBP (International Geosphere-Biosphere Programme) land cover classification scheme; (2) derive continuous surfaces of LAI that capture the mean and variability of the LAI field measurements; and (3) conduct initial MODIS validation exercises to assess the quality of early (i.e., provisional) MODIS products. ETM+ land cover maps varied in overall accuracy from 81% to 95%. The boreal forest was the most spatially complex, had the greatest number of classes, and the lowest accuracy. The intensive agricultural cropland had the simplest spatial structure, the least number of classes, and the highest overall accuracy. At each site, mapped LAI patterns generally followed patterns of land cover across the site. Predicted versus observed LAI indicated a high degree of correspondence between field-based measures and ETM+ predictions of LAI. Direct comparisons of ETM+ land cover maps with Collection 3 MODIS cover maps revealed several important distinctions and similarities. One obvious difference was associated with image/map resolution. ETM+ captured much of the spatial complexity of land cover at the sites. In contrast, the relatively coarse resolution of MODIS did not allow for that level of spatial detail. Over the extent of all sites, the greatest difference was an overprediction by MODIS of evergreen needleleaf forest cover at the boreal forest site, which consisted largely of open shrubland, woody savanna, and savanna. At the agricultural, temperate mixed forest, and prairie grassland sites, ETM+ and MODIS cover estimates were similar. Collection 3 MODIS-based LAI estimates were considerably higher (up to 4 m² m⁻²) than those based on ETM+ LAI at each site. There are numerous probable reasons for this, the most important being the algorithms' sensitivity to MODIS reflectance calibration, its use of a prelaunch AVHRR-based land cover map, and its apparent reliance on mainly red and near-IR reflectance. Samples of Collection 4 LAI products were examined and found to consist of significantly improved LAI predictions for KONZ, and to some extent for AGRO, but not for the other two sites. In this study, we demonstrate that MODIS reflectance data are highly correlated with LAI across three study sites, with relationships increasing in strength from 500 to 1000 m spatial resolution, when shortwave-infrared bands are included.     

Evaluation of MODIS-based estimates of water-use efficiency in Amazonia

This study aimed to assess the spatial-temporal patterns of water-use efficiency (WUE) obtained through MODIS gross primary productivity (GPP) and evapotranspiration (ET) products (MOD17 for GPP and MOD16 for ET) in the Upper Tapajos and Curua-Una River basins, located in the oriental flank of the Amazon region, and to validate the results with flux tower measurements within the Large-Scale Biosphere-Atmosphere Experiment in Amazonia (LBA) project. The spatial variation of WUE was primarily related to the larger presence of forested areas in the Upper Tapajos River basin (western part) compared with the Curua-Una River basin (eastern part), which is situated within the so-called arc of deforestation. Temporally, WUE showed a pronounced seasonal pattern, varying with the dry and wet seasons in the region. A decrease of ~3% in WUE was observed during the dry season, which was related to the low water availability and increased vapour pressure deficit during the dry period, which induces stomatal closure, leading to a decline in the photosynthetic rate. Comparison between the WUE estimates obtained by MODIS data and observations from the LBA towers showed an average error of 17%, varying between ~12% and ~28% for the different sites. MODIS WUE depends on the accuracy of both GPP and ET estimation. In this sense, we highlight that improvements in both MODIS GPP and ET products are necessary and should focus on reducing the uncertainties related to the biophysical vegetation parameters and meteorological data that serve as input information in the algorithms.     

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

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

Estimation of gross primary production in Moso bamboo forest based on light-use efficiency derived from MODIS reflectance data

Assessing the contribution of Moso bamboo (Phyllostachys pubescens) forest to forest ecosystem carbon storage requires accurate estimation of gross primary production (GPP). Based on measurements of light-use efficiency (LUE), defined as the ratio of measured GPP to photosynthetically active radiation (PAR), from the eddy covariance flux tower, the linear regression model and partial least squares regression model were used for estimation of LUE using the Moderate-Resolution Imaging Spectroradiometer (MODIS) reflectance data. GPP estimates were then calculated by the product of LUE estimates and PAR (named the LUE-PAR model), which was compared with GPP from the GPP algorithm designed for the MODIS sensor aboard the Aqua and Terra platforms (MOD17A2 model) and the EC-LUE model. The results revealed the PLS model performed better than the linear regression model in LUE estimation but had lager uncertainties in high and low LUE values. GPP estimates driven by a MODIS-based radiation product with high spatial resolution was more accurate than those driven by Modern-Era Retrospective Analysis for Research and Applications (MERRA) radiation product from the NASA’s Global Modelling and Assimilation Office data set. The LUE-PAR model had the highest accuracy than the other two LUE models. The GPP values derived from the EC-LUE model driven by photosynthetically active radiation (PAR) from MERRA and maximum LUE from the EC data were overestimated due to the overestimation in MERRA radiation product. The GPP values derived from the MOD17A2 model driven by PAR from the MERRA and maximum LUE from the biome properties look-up table were underestimated due to underestimation in the maximum LUE of Moso bamboo forest. This study implied that the LUE-PAR model driven by LUE estimates from the PLS model and PAR from MERRA is a superior approach in improving GPP simulations, and PAR products with high spatial resolution and accurate species-specific maximum LUE are necessary for the LUE models in estimating GPP at regional scale.     

Land cover mapping in support of LAI and FPAR retrievals from EOS-MODIS and MISR: Classification methods and sensitivities to errors

Land cover maps are used widely to parameterize the biophysical properties of plant canopies in models that describe terrestrial biogeochemical processes. In this paper, we describe the use of supervised classification algorithms to generate land cover maps that characterize the vegetation types required for Leaf Area Index (LAI) and Fraction of Photosynthetically Active Radiation (FPAR) retrievals from MODIS and MISR. As part of this analysis, we examine the sensitivity of remote sensing-based retrievals of LAI and FPAR to land cover information used to parameterize vegetation canopy radiative transfer models. Specifically, a decision tree classification algorithm is used to generate a land cover map of North America from Advanced Very High Resolution Radiometer (AVHRR) data with 1 km spatial resolution using a six-biome classification scheme. To do this, a time series of normalized difference vegetation index data from the AVHRR is used in association with extensive site-based training data compiled using Landsat Thematic Mapper (TM) and ancillary map sources. Accuracy assessment of the map produced via decision tree classification yields a cross-validated map accuracy of 73%. Results comparing LAI and FPAR retrievals using maps from different sources show that disagreement in land cover labels generally do not translate into strong disagreement in LAI and FPAR maps. Further, the main source of disagreement in LAI and FPAR maps can be attributed to specific biome classes that are characterized by a continuum of fractional cover and canopy structure.     

Downscaling real-time vegetation dynamics by fusing multi-temporal MODIS and Landsat NDVI in topographically complex terrain

Canopy phenology is an important factor driving seasonal patterns of water and carbon exchange between land surface and atmosphere. Recent developments of real-time global satellite products (e.g., MODIS) provide the potential to assimilate dynamic canopy measurements with spatially distributed process-based ecohydrological models. However, global satellite products usually are provided with relatively coarse spatial resolutions, averaging out important spatial heterogeneity of both terrain and vegetation. Therefore, bias can result from lumped representation of ecological and hydrological processes especially in topographically complex terrain. Successful downscaling of canopy phenology to high spatial resolution would be indispensable for catchment-scale distributed ecohydrological modeling, aiming at understanding complex patterns of water, carbon and nutrient cycling in mountainous watersheds. Two downscaling approaches are developed in this study to overcome this issue by fusing multi-temporal MODIS and Landsat TM data in conjunction with topographic information to estimate high spatio-temporal resolution biophysical parameters over complex terrain. MODIS FPAR (fraction of absorbed photosynthetically active radiation) is used to provide medium spatial resolution phenology, while the variability of vegetation within a MODIS pixel is characterized by Landsat NDVI. The algorithms depend on the scale-invariant linear relationship between FPAR and NDVI, which is verified in this study. Downscaled vegetation dynamics are successfully validated both temporally and spatially with ground-based continuous FPAR and leaf area index measurements. Topographic correction during the downscaling process has a limited effect on downscaled FPAR products except for the period around the winter solstice in the study area.