Development of a global evapotranspiration algorithm based on MODIS and global meteorology data

The objective of this research is to develop a global remote sensing evapotranspiration (ET) algorithm based on Cleugh et al.'s [Cleugh, H.A., R. Leuning, Q. Mu, S.W. Running (2007) Regional evaporation estimates from flux tower and MODIS satellite data. Remote Sensing of Environment 106, page 285-304- 2007 (doi: 10.1016/j.rse.2006.07.007).] Penman-Monteith based ET (RS-PM). Our algorithm considers both the surface energy partitioning process and environmental constraints on ET. We use ground-based meteorological observations and remote sensing data from the MODerate Resolution Imaging Spectroradiometer (MODIS) to estimate global ET by (1) adding vapor pressure deficit and minimum air temperature constraints on stomatal conductance; (2) using leaf area index as a scalar for estimating canopy conductance; (3) replacing the Normalized Difference Vegetation Index with the Enhanced Vegetation Index thereby also changing the equation for calculation of the vegetation cover fraction (F_C); and (4) adding a calculation of soil evaporation to the previously proposed RS-PM method.We evaluate our algorithm using ET observations at 19 AmeriFlux eddy covariance flux towers. We calculated ET with both our Revised RS-PM algorithm and the RS-PM algorithm using Global Modeling and Assimilation Office (GMAO v. 4.0.0) meteorological data and compared the resulting ET estimates with observations. Results indicate that our Revised RS-PM algorithm substantially reduces the root mean square error (RMSE) of the 8-day latent heat flux (LE) averaged over the 19 towers from 64.6 W/m² (RS-PM algorithm) to 27.3 W/m² (Revised RS-PM) with tower meteorological data, and from 71.9 W/m² to 29.5 W/m² with GMAO meteorological data. The average LE bias of the tower-driven LE estimates to the LE observations changed from 39.9 W/m² to − 5.8 W/m² and from 48.2 W/m² to − 1.3 W/m² driven by GMAO data. The correlation coefficients increased slightly from 0.70 to 0.76 with the use of tower meteorological data. We then apply our Revised RS-PM algorithm to the globe using 0.05° MODIS remote sensing data and reanalysis meteorological data to obtain the annual global ET (MODIS ET) for 2001. As expected, the spatial pattern of the MODIS ET agrees well with that of the MODIS global terrestrial gross and net primary production (MOD17 GPP/NPP), with the highest ET over tropical forests and the lowest ET values in dry areas with short growing seasons. This MODIS ET product provides critical information on the regional and global water cycle and resulting environment changes.     

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.     

A continuous measure of gross primary production for the conterminous United States derived from MODIS and AmeriFlux data

The quantification of carbon fluxes between the terrestrial biosphere and the atmosphere is of scientific importance and also relevant to climate-policy making. Eddy covariance flux towers provide continuous measurements of ecosystem-level exchange of carbon dioxide spanning diurnal, synoptic, seasonal, and interannual time scales. However, these measurements only represent the fluxes at the scale of the tower footprint. Here we used remotely sensed data from the Moderate Resolution Imaging Spectroradiometer (MODIS) to upscale gross primary productivity (GPP) data from eddy covariance flux towers to the continental scale. We first combined GPP and MODIS data for 42 AmeriFlux towers encompassing a wide range of ecosystem and climate types to develop a predictive GPP model using a regression tree approach. The predictive model was trained using observed GPP over the period 2000-2004, and was validated using observed GPP over the period 2005-2006 and leave-one-out cross-validation. Our model predicted GPP fairly well at the site level. We then used the model to estimate GPP for each 1 km × 1 km cell across the U.S. for each 8-day interval over the period from February 2000 to December 2006 using MODIS data. Our GPP estimates provide a spatially and temporally continuous measure of gross primary production for the U.S. that is a highly constrained by eddy covariance flux data. Our study demonstrated that our empirical approach is effective for upscaling eddy flux GPP data to the continental scale and producing continuous GPP estimates across multiple biomes. With these estimates, we then examined the patterns, magnitude, and interannual variability of GPP. We estimated a gross carbon uptake between 6.91 and 7.33 Pg C yr^− 1 for the conterminous U.S. Drought, fires, and hurricanes reduced annual GPP at regional scales and could have a significant impact on the U.S. net ecosystem carbon exchange. The sources of the interannual variability of U.S. GPP were dominated by these extreme climate events and disturbances.     

Estimating daily gross primary production of maize based only on MODIS WDRVI and shortwave radiation data

Accurate assessment of temporal changes in gross primary production (GPP) is important for carbon budget assessments and evaluating the impact of climate change on crop productivity. The objective of this study was to devise a simple remote sensing-based GPP model to quantify daily GPP of maize. In the model, (1) daily shortwave radiation (SW), derived from the reanalysis data (North American Land Data Assimilation System; NLDAS-2) and (2) smoothed Wide Dynamic Range Vegetation Index (WDRVI) data, derived from Moderate Resolution Imaging Spectroradiometer (MODIS) 250-m observations were used as proxy variables of the incident photosynthetically active radiation (PAR) and the total canopy chlorophyll content, respectively. The model was calibrated and validated by using tower-based CO₂ flux observations over an 8-year period (2001 to 2008) for one rainfed and two irrigated sites planted to maize as part of the Carbon Sequestration Program at the University of Nebraska-Lincoln. The results showed the temporal features of the product SW*WDRVI closely related to the temporal GPP variations in terms of both daily variations and seasonal patterns. The simple GPP model was able to predict the daily GPP values and accumulated GPP values of maize with high accuracy.     

Comparison of multiple models for estimating gross primary production using MODIS and eddy covariance data in Harvard Forest

Gross primary production (GPP) defined as the overall rate of fixation of carbon through the process of vegetation photosynthesis is important for carbon cycle and climate change research. Three models, the Vegetation Photosynthesis Model (VPM), the Temperature and Greenness (TG) model and the Vegetation Index (VI) model have been compared for the estimation of GPP in Harvard Forest from 2003 to 2006 using climate variables acquired by eddy covariance (EC) measurements and Moderate Resolution Imaging Spectroradiometer (MODIS) satellite images. All these models provide more reliable estimates of GPP than that of MODIS GPP product. High Pearsons correlation coefficients r equal to 0.94, 0.92 and 0.90 are observed for the VPM, the TG and the VI model, respectively. Relationships between GPP and land surface temperature (LST, R² = 0.72), and vapor pressure deficit (VPD, R² = 0.45) indicate that climate variables are important for GPP estimation. Due to proper characterization of temperature, water stress and leaf age by three scalars, VPM best follows the seasonal variations of GPP. By incorporation of the MODIS surface reflectance and LST product, the TG model is the most suitable choice for areas without prior knowledge as it is based entirely on remote sensing observations. Results from the VI model demonstrate the possibility of using a single vegetation index for light use efficiency (LUE) estimation in deciduous forest that is of high spatial heterogeneity. The validation and comparison of models will be helpful in development of future GPP models using combinations of climate variables and/or remote sensing observations.     

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.     

Developing a continental-scale measure of gross primary production by combining MODIS and AmeriFlux data through Support Vector Machine approach

Remote sensing is a potentially powerful technology with which to extrapolate eddy covariance-based gross primary production (GPP) to continental scales. In support of this concept, we used meteorological and flux data from the AmeriFlux network and Support Vector Machine (SVM), an inductive machine learning technique, to develop and apply a predictive GPP model for the conterminous U.S. In the following four-step process, we first trained the SVM to predict flux-based GPP from 33 AmeriFlux sites between 2000 and 2003 using three remotely-sensed variables (land surface temperature, enhanced vegetation index (EVI), and land cover) and one ground-measured variable (incident shortwave radiation). Second, we evaluated model performance by predicting GPP for 24 available AmeriFlux sites in 2004. In this independent evaluation, the SVM predicted GPP with a root mean squared error (RMSE) of 1.87 gC/m²/day and an R² of 0.71. Based on annual total GPP at 15 AmeriFlux sites for which the number of 8-day averages in 2004 was no less than 67% (30 out of a possible 45), annual SVM GPP prediction error was 32.1% for non-forest ecosystems and 22.2% for forest ecosystems, while the standard Moderate Resolution Imaging Spectroradiometer GPP product (MOD17) had an error of 50.3% for non-forest ecosystems and 21.5% for forest ecosystems, suggesting that the regionally tuned SVM performed better than the standard global MOD17 GPP for non-forest ecosystems but had similar performance for forest ecosystems. The most important explanatory factor for GPP prediction was EVI, removal of which increased GPP RMSE by 0.85 gC/m²/day in a cross-validation experiment. Third, using the SVM driven by remote sensing data including incident shortwave radiation, we predicted 2004 conterminous U.S. GPP and found that results were consistent with expected spatial and temporal patterns. Finally, as an illustration of SVM GPP for ecological applications, we estimated maximum light use efficiency (e_max), one of the most important factors for standard light use efficiency models, for the conterminous U.S. by integrating the 2004 SVM GPP with the MOD17 GPP algorithm. We found that e_max varied from ∼ 0.86 gC/MJ in grasslands to ∼ 1.56 gC/MJ in deciduous forests, while MOD17 e_max was 0.68 gC/MJ for grasslands and 1.16 gC/MJ for deciduous forests, suggesting that refinements of MOD17 e_max may be beneficial.     

Remote estimation of gross primary production in maize and support for a new paradigm based on total crop chlorophyll content

The accurate quantification of gross primary production (GPP) in crops is important for regional and global studies of carbon budgets. Because of the observed close relationship between GPP and total canopy chlorophyll content in crops, vegetation indices related to chlorophyll can be used as a proxy of GPP. In this study, we justified the approach, tested the performance of several widely used chlorophyll-related vegetation indices in estimating total chlorophyll content and GPP in maize based on spectral data collected at a close range, 6 meters above the top of the canopy, over a period of eight years (2001 to 2008). The results show that GPP can be accurately estimated with chlorophyll-related indices that use near infra-red and either green or the red edge range of the spectrum. These indices provide the best approximation of the widely variable GPP in maize under both irrigated and rainfed conditions.     

Modeling gross primary production of alpine ecosystems in the Tibetan Plateau using MODIS images and climate data

The eddy covariance technique provides measurements of net ecosystem exchange (NEE) of CO₂ between the atmosphere and terrestrial ecosystems, which is widely used to estimate ecosystem respiration and gross primary production (GPP) at a number of CO₂ eddy flux tower sites. In this paper, canopy-level maximum light use efficiency, a key parameter in the satellite-based Vegetation Photosynthesis Model (VPM), was estimated by using the observed CO₂ flux data and photosynthetically active radiation (PAR) data from eddy flux tower sites in an alpine swamp ecosystem, an alpine shrub ecosystem and an alpine meadow ecosystem in Qinghai-Tibetan Plateau, China. The VPM model uses two improved vegetation indices (Enhanced Vegetation Index (EVI), Land Surface Water Index (LSWI)) derived from the Moderate Resolution Imaging Spectral radiometer (MODIS) data and climate data at the flux tower sites, and estimated the seasonal dynamics of GPP of the three alpine grassland ecosystems in Qinghai-Tibetan Plateau. The seasonal dynamics of GPP predicted by the VPM model agreed well with estimated GPP from eddy flux towers. These results demonstrated the potential of the satellite-driven VPM model for scaling-up GPP of alpine grassland ecosystems, a key component for the study of the carbon cycle at regional and global scales.     

Evaluation and comparison of gross primary production estimates for the Northern Great Plains grasslands

Two spatially-explicit estimates of gross primary production (GPP) are available for the Northern Great Plains. An empirical piecewise regression (PWR) GPP model was developed from flux tower measurements to map carbon flux across the region. The Moderate Resolution Imaging Spectrometer (MODIS) GPP model is a process-based model that uses flux tower data to calibrate its parameters. Verification and comparison of the regional PWR GPP and the global MODIS GPP are important for the modeling of grassland carbon flux. This study compared GPP estimates from PWR and MODIS models with five towers in the grasslands. Among them, PWR GPP and MODIS GPP showed a good agreement with tower-based GPP at three towers. The global MODIS GPP, however, did not agree well with tower-based GPP at two other towers, probably because of the insensitivity of MODIS model to regional ecosystem and climate change and extreme soil moisture conditions. Cross-validation indicated that the PWR model is relatively robust for predicting regional grassland GPP. However, the PWR model should include a wide variety of flux tower data as the training data sets to obtain more accurate results.In addition, GPP maps based on the PWR and MODIS models were compared for the entire region. In the northwest and south, PWR GPP was much higher than MODIS GPP. These areas were characterized by the higher water holding capacity with a lower proportion of C₄ grasses in the northwest and a higher proportion of C₄ grasses in the south. In the central and southeastern regions, PWR GPP was much lower than MODIS GPP under complicated conditions with generally mixed C₃/C₄ grasses. The analysis indicated that the global MODIS GPP model has some limitations on detecting moisture stress, which may have been caused by the facts that C₃ and C₄ grasses are not distinguished, water stress is driven by vapor pressure deficit (VPD) from coarse meteorological data, and MODIS land cover data are unable to differentiate the sub-pixel cropland components.     

Scaling and assessment of GPP from MODIS using a combination of airborne lidar and eddy covariance measurements over jack pine forests

Understanding the influence of within-pixel land cover heterogeneity is essential for the extrapolation of measured and modeled CO₂ fluxes from the canopy to regional scales using remote sensing. Airborne light detection and ranging (lidar) was used to estimate spatial and temporal variations of gross primary production (GPP) across a jack pine chronosequence of four sites in Saskatchewan, Canada for comparison with the Moderate Resolution Imaging Spectroradiometer (MODIS) GPP product. This study utilizes high resolution canopy structural information obtained from airborne lidar to bridge gaps in spatial representation between plot, eddy covariance (EC), and MODIS estimates of vegetation GPP. First we investigate linkages between canopy structure obtained from measurements and light response curves at a jack pine chronosequence during the growing season of 2004. Second, we use the measured canopy height and foliage cover inputs to create a structure-based GPP model (GPP_Landsberg) which was tested in 2005. The GPP model is then run using lidar data (GPP_Lidar) and compared with eight-day cumulative MODIS GPP (GPP_MODIS) and EC observations (GPP_EC). Finally, we apply the lidar GPP model at spatial resolutions of 1 m to 1000 m to examine the influence of within-pixel heterogeneity and scaling (or pixel aggregation) on GPP_Lidar. When compared over eight-day cumulative periods throughout the 2005 growing season, the standard deviation of differences between GPP_lidar and GPP_MODIS were less than differences between either of them and GPP_EC at all sites. As might be expected, the differences between pixel aggregated GPP estimates are most pronounced at sites with the highest levels of spatial canopy heterogeneity. The results of this study demonstrate one method for using lidar to scale between eddy covariance flux towers and coarse resolution remote sensing pixels using a structure-based Landsberg light curve model.     

Exploring the potential of MODIS EVI for modeling gross primary production across African ecosystems

One of the most frequently applied methods for integrating controls on primary production through satellite data is the light use efficiency (LUE) approach, which links vegetation gross or net primary productivity (GPP or NPP) to remotely sensed estimates of absorbed photosynthetically active radiation (APAR). Eddy covariance towers provide continuous measurements of carbon flux, presenting an opportunity for evaluation of satellite estimates of GPP. Here we investigate relationships between eddy covariance estimated GPP, environmental variables derived from flux towers, Moderate Resolution Imaging Spectroradiometer (MODIS) enhanced vegetation index (EVI) and GPP across African savanna ecosystems. MODIS GPP was found to underestimate GPP at the majority of sites, particularly at sites in the Sahel. EVI was found to correlate well with estimated GPP on a site-by-site basis. Combining EVI with tower-measured PAR and evaporative fraction (EF, a measure of water sufficiency) improved the direct relationship between GPP and EVI at the majority of the sites. The slope of this relationship was strongly related to site peak leaf area index (LAI). These results are promising for the extension of GPP through the use of remote sensing data to a regional or even continental scale.     

Comparison of MODIS gross primary production estimates for forests across the U.S.A. with those generated by a simple process model, 3-PGS

Since 2000, NASA's Moderate Resolution Imaging Spectro-radiometer (MODIS) has provided 1 × 1 km estimates of 8-day gross primary production (GPP). The MODIS algorithm computes GPP as a simple function of absorbed photosynthetically active radiation and a regionally assigned light-use conversion efficiency (LUE) that is reduced if temperature or atmospheric vapor pressure deficits are suboptimal. We compared MODIS-derived GPP estimates for forested areas across the United States of America (U.S.A.) with those generated by the 3-PGS (Physiological Principles Predicting Growth using Satellite data) model, the latter of which considers spatial variation in available soil water storage capacity (ASWC) and nitrogen content. We expected seasonal and annual MODIS GPP values to be in close agreement with those derived from the 3-PGS model in regions with adequate precipitation, soil water storage, and moderately fertile soils. 3-PGS was initially run with STATSGO-derived soils information provided by the Oak Ridge National Laboratory. The analysis was expanded to include sensitivity analyses with ASWC set at 50, 100, 300, and 400 mm to identify areas within nine major ecoregions where drought might prove to be a major limitation on GPP. The majority of forests across the U.S.A. were relatively insensitive to large variations in ASW storage. In areas where ASWC was assumed < 200 mm and average annual rainfall was < 100 mm yr^− 1, GPP was predicted to be reduced by > 60%. There was generally good agreement (within 20%) between MODIS and 3-PGS estimates of forest GPP across the U.S.A. GPP predicted by the MODIS model was generally higher in ecoregions with substantial drought and with relatively low soil fertility. The latter, which influences LUE, was more than twice as important as soil drought.     

Comparison of MODIS, eddy covariance determined and physiologically modelled gross primary production (GPP) in a Douglas-fir forest stand

Quantification of the magnitude of net terrestrial carbon (C) uptake, and how it varies inter-annually, is an important question with future potential sequestration influenced by both increased atmospheric CO₂ and changing climate. However the assessment of differences in measured and modeled C accumulation is a challenging task due to the significant fine scale variation occurring in terrestrial productivity due to soil, climate and vegetation characteristics as well as difficulties in measuring carbon accumulation over large spatial areas. The Moderate Resolution Imaging Spectroradiometer (MODIS) offers a means of monitoring gross primary production (GPP), both spatially and temporally, routinely from space. However it is critical to compare and contrast the temporal dynamics of the C and water fluxes with those measured from ground-based networks, or estimated using physiological models. In this paper, using a number of approaches, our objective is to determine if any systematic biases exists in either the MODIS, or the modeled estimates of fluxes, relative to the measurements made over an evergreen, needleleaf temperate rainforest on Vancouver Island, Canada. Results indicate that 8-day GPP as predicted with a simple physiological model (3PGS), forced using local meteorology and canopy characteristics, matched measured fluxes very well (r² = 0.86, p < 0.001) with no significant difference between eddy covariance (EC) and modeled GPP (p < 0.001). In addition, modeled water supply closely matched measured relative available soil water content at the site. Using canopy characteristics from the MODIS fraction of photosynthetically active radiation (fPAR) algorithm, slightly reduced the correspondence of the predictions due to a large number of unsuccessful retrievals (83%) due to sun angle, snow and cloud. Predictions of GPP based on the MODIS GPP algorithm, forced using local meteorology and canopy characteristics, were also highly correlated with EC measurements (r² = 0.89, p < 0.001) however these estimates were biased under predicting GPP. Estimates of GPP based on the most recent MODIS reprocessing (collection 4.5) remained highly correlated (r² = 0.88, p < 0.001) yet were also the most biased with the estimates being 30% less than the EC-measured GPP. Most of the variance in GPP at the site was explained by the absorbed photosynthetically active radiation. We also compared the nighttime respiration as measured over 2 years at the site with the minimum 8-day MODIS land surface temperature and found a significant relationship (r² = 0.57), similar to other studies.     

Comparison of three models for predicting gross primary production across and within forested ecoregions in the contiguous United States

Gross primary production (GPP), the photosynthetic uptake of carbon, is an important variable in the global carbon cycle. Although continuous measurements of GPP are being collected from a network of micro-meteorological towers, each site represents a small area with records available for only a limited period. As a result, GPP is commonly modeled over forested landscapes as a function of climatic and soil variables, often supplemented with satellite-derived estimates of the vegetation's light-absorbing properties. Since the late 1990s, a number of models have been developed to provide seasonal and annual estimates of GPP across much of the Earth. Each model, however, contains different underlying assumptions and requires different amounts of data. As a result, predictions vary, sometimes significantly. In this paper we compare modeled estimates of GPP for forested areas across the U.S.A. derived from: NASA's MODIS Product (MOD17); the C-Fix model using SPOT-VGT satellite-derived vegetation data; and the Physiological Principles Predicting Growth from Satellites (3-PGS) model, a process-based model that requires information on both climate and soil properties. The models predicted average ecoregion values of forest GPP between 9.8 and 14.1 MgC ha^− 1 y^− 1 across the United States. 3-PGS predicted the lowest values while the C-Fix model, which included a CO₂ fertilization factor, produced the highest estimates. In the western part of the country, estimates of GPP within a given ecoregion varied by as much as 50%, whereas in the northeast, where topography and climate are less extreme, variation in GPP was less than 10%. Within ecoregions, 3PGS predicted the most variation, reflecting its sensitivity to variation in soil properties. We conclude that where model predictions disagree, an opportunity is presented to evaluate underlying assumptions through sensitivity analyses, additional data collection and where more detailed study is warranted.     

Predicting riparian evapotranspiration from MODIS vegetation indices and meteorological data

A vegetation index (VI) model for predicting evapotranspiration (ET) from data from the Moderate Resolution Imaging Spectrometer (MODIS) on the EOS-1 Terra satellite and ground meteorological data was developed for riparian vegetation along the Middle Rio Grande River in New Mexico. Ground ET measurements obtained from eddy covariance towers at four riparian sites were correlated with MODIS VIs, MODIS land surface temperatures (LSTs), and ground micrometeorological data over four years. Sites included two saltcedar (Tamarix ramosissima) and two Rio Grande cottonwood (Populus deltoides ssp. Wislizennii) dominated stands. The Enhanced Vegetation Index (EVI) was more closely correlated (r=0.76) with ET than the Normalized Difference Vegetation Index (NDVI; r=0.68) for ET data combined over sites and species. Air temperature (T_a) measured over the canopy from towers was the meteorological variable that was most closely correlated with ET (r=0.82). MODIS LST data at 1- and 5-km resolutions were too coarse to accurately measure the radiant surface temperature within the narrow riparian corridor; hence, energy balance methods for estimating ET using MODIS LSTs were not successful. On the other hand, a multivariate regression equation for predicting ET from EVI and T_a had an r²=0.82 across sites, species, and years. The equation was similar to VI-ET models developed for crop species. The finding that ET predictions did not require species-specific equations is significant, inasmuch as these are mixed vegetation zones that cannot be easily mapped at the species level.     

Predicting gross primary production from the enhanced vegetation index and photosynthetically active radiation: Evaluation and calibration

The approach of using primarily satellite observations to estimate ecosystem gross primary production (GPP) without resorting to interpolation of many surface observations has recently shown promising results. Previous work has shown that the remote sensing based greenness and radiation (GR) model can give accurate GPP estimates in crops. However, the feasibility of its application and the model calibration to other ecosystems remain unknown. With the enhanced vegetation index (EVI) derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) images and the surface based estimates of photosynthetically active radiation (PAR), we provide an analysis of the GR model for estimating monthly GPP using flux measurements at fifteen sites, representing a wide range of ecosystems with various canopy structures and climate characteristics. Results demonstrate that the GR model can provide better estimates of GPP than that of the temperature and greenness (TG) model for the overall data classified as non-forest (NF), deciduous forest (DF) and evergreen forest (EF) sites. Calibration of the GR model is also conducted and has shown reasonable results for all sites with a root mean square error of 47.18 g C/m²/month. Different coefficients acquired for the three plant functional types indicate that there are shifts of importance among various factors that determine the monthly vegetation GPP. The analysis firstly shows the potential use of the GR model in estimating GPP across biomes while it also points to the needs of further considerations in future operational applications.     

Use of remotely sensed and ancillary data for estimating forest gross primary productivity in Italy

The current paper describes the development and testing of a procedure which can use widely available remotely sensed and ancillary data to assess large-scale patterns of forest productivity in Italy. To reach this objective a straightforward model (C-Fix) was applied which is based on the relationship between photosynthetically active radiation absorbed by plant canopies and relevant gross primary productivity (GPP). The original C-Fix methodology was improved by using more abundant ancillary information and more efficient techniques for NDVI data processing. In particular, two extraction methods were applied to NDVI data, derived from two sensors (NOAA-AVHRR and SPOT-VGT) to feed C-Fix. The accuracy of the model outputs was assessed through comparison with annual and monthly values of forest GPP derived from eight eddy covariance flux towers. The results obtained indicated the superiority of SPOT-VGT over NOAA-AVHRR data and a higher efficiency of the more advanced NDVI extraction method. Globally, the procedure was proved to be of easy and objective implementation and allowed the evaluation of mean productivity levels of existing forests on the national scale.     

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.     

Deriving maximal light use efficiency from coordinated flux measurements and satellite data for regional gross primary production modeling

Remote sensing models based on light use efficiency (LUE) provide promising tools for monitoring spatial and temporal variation of gross primary production (GPP) at regional scale. In most of current LUE-based models, maximal LUE (ε_max) heavily relies on land cover types and is considered as a constant, rather than a variable for a certain vegetation type or even entire eco-region. However, species composition and plant functional types are often highly heterogeneous in a given land cover class; therefore, spatial heterogeneity of ε_max must be fully considered in GPP modeling, so that a single cover type does not equate to a single ε_max value. A spatial dataset of ε_max accurately represents the spatial heterogeneity of maximal light use would be of significant beneficial to regional GPP models. Here, we developed a spatial dataset of ε_max by integrating eddy covariance flux measurements from 14 field sites in a network of coordinated observation across northern China and satellite derived indices such as enhanced vegetation index (EVI) and visible albedo to simulate regional distribution of GPP. This dynamic modeling method recognizes the spatial heterogeneity of ε_max and reduces the uncertainties in mixed pixels. Further, we simulated GPP with the spatial dataset of ε_max generated above. Both ε_max and growing season GPP show complex patterns over northern China that reflect influences of humidity, green vegetation fractions, and land use intensity. “Green spots” such as oasis meadow and alpine forests in dryland and “brown spots” such as build-up and heavily degraded vegetation in the east are clearly captured by the simulation. The correlation between simulated GPP and EC measured GPP indicate that the simulated GPP from this new approach is well matched with flux-measured GPP. Those results have demonstrated the importance of considering ε_max as both a spatially and temporally variable values in GPP modeling.