Simplified atmospheric radiative transfer modelling for estimating incident PAR using MODIS atmosphere products

Photosynthetically active radiation (PAR) is an important parameter in the estimation of vegetation growth. A method is presented with which instantaneous PAR can be calculated with high accuracy from Moderate Resolution Imaging Spectroradiometer (MODIS) atmosphere data products. The method, PARcalc, is based on a simplification of the general radiative transfer equations, whereby the atmosphere is modelled in a few layers depending on the opacity of the atmosphere due to clouds and aerosol load. The four major processes of atmospheric attenuation of solar radiation, Rayleigh scattering, absorption by ozone and water vapour, and aerosol scattering, are treated individually, using calibrated data from MODIS imagery. Comparing the calculations to field measurements made in Costa Rica in 2002, absolute errors in the order of 2-3% are obtained, which should be sufficiently accurate for application of the PARcalc method in models of vegetation growth as a source of incident PAR.     

Estimation of absorbed PAR across Scandinavia from satellite measurements: Part I: Incident PAR

The primary productivity of a plant community can be modeled as the product of the amount of photosynthetically active radiation (PAR) absorbed by the canopy and a light use efficiency factor, where the amount of absorbed PAR (APAR) is the product of the fractional absorption and the amount of incident PAR. By implementing a method, PARcalc, using atmospheric data from the Moderate Resolution Imaging Spectroradiometer (MODIS), incident PAR is estimated in this study. In addition, since many PAR datasets are generated by converting shortwave radiation into PAR, the ratio of PAR to shortwave radiation was also investigated. PARcalc models the photosynthetic photon flux density (PPFD) as a product of atmospheric transmittance, the cosine of the Sun zenith angle, and the solar constant. The atmospheric transmittance includes the attenuation of radiation by Rayleigh and aerosol scattering, and absorption by water and ozone. A cloud transmittance factor which is primarily a function of the cloud optical thickness is added in order to cope with cloudy conditions. The model was implemented at two sites in Sweden, Asa and Norunda, where in situ measurements of PPFD were made during the spring and summer of 2004. Modeled time-series were evaluated against the measurements, and daily sums of PPFD were calculated by fitting of a sine function in combination with linear interpolation of the instantaneous estimates from sunrise to sunset. This gave correlation coefficients at Norunda and Asa of 0.80 and 0.77, respectively, when comparing modeled and measured daily insolation. The average relative errors were 24% and 25%. Corresponding figures for five day averages were 0.91 and 0.86; and 9.3% and 11.9%. Instantaneous estimates of PPFD were modeled with correlation coefficients of 0 88-0 93 and average relative errors from 17.0%. These numbers were acquired when using measured values for determining cloudiness; the corresponding figures when the method is fully implemented using satellite data are 0.84 to 0.71 and 24.9%, respectively. The ratio of PAR to shortwave radiation was measured at Norunda 1 Jan to 31 Oct 2004 and was found to vary between 0.27 and 0.48 on a daily basis with an average of 0.43 for the whole period.     

Mapping evapotranspiration using MODIS and MM5 Four-Dimensional Data Assimilation

Evapotranspiration (ET), the sum of evaporation from soil and transpiration from vegetation, is of vital importance in the hydrologic cycle and must be taken into consideration in assessments of the water resources of any region. The MODerate resolution Imaging Spectroradiometer (MODIS) sensor offers a promising opportunity for estimating daily ET with a 1 km spatial resolution, but is hampered by frequent cloud contamination or data gaps from other factors. In this study, 1) a stand-alone ET model was applied and tested during clear or partial cloudy sky conditions using MODIS-based inputs of land surface and atmospheric data and 2) meteorological simulations by using Four-Dimensional Data Assimilation (FDDA) system between MODIS and the 5th Generation Meso-scale Meteorological Model (MM5) was used in cloudy conditions to facilitate continuous daily ET estimates. The MODIS ET algorithm modified from Mu et al. (2007) is based on the Penman-Monteith equation and was applied to predict ET at flux measurement sites. This algorithm considers both the effects of surface energy partitioning processes and environmental constraints on ET. We devised gap-filling approaches for MODIS aerosol and albedo data that were identified as bottlenecks to determine retrieval rates of insolation and ET. MODIS-derived input variables (i.e., meteorological variables and radiation components) for estimating ET showed a good agreement with flux tower observations at each site. The retrieval rate of MODIS ET doubled at four flux measurement sites after gap-filling with negligible compensation was undertaken for accuracy. In spite of the high accuracy of MODIS-derived input variables, MODIS ET showed meaningful errors at the four flux measurement sites. These errors were mainly associated with errors in the estimated canopy conductance. During clear sky conditions, MODIS was used to calculate ET, while the MODIS-MM5 FDDA system provided input variables for the calculation of ET under cloudy sky conditions. The performance of the MODIS-MM5 FDDA system was evaluated by comparing ET based on MODIS, which showed a good agreement with the MODIS ET for various land cover types. Our results indicate that MODIS can be applied to monitor the land surface energy budget and ET with reasonable accuracy and that MODIS-MM5 FDDA has the potential to provide reasonable input data of ET estimation under cloudy conditions.     

Mapping incident photosynthetically active radiation from MODIS data over China

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

Statistical evaluation of remotely sensed snow-cover products with constraints from streamflow and SNOTEL measurements

Using streamflow and Snowpack Telemetry (SNOTEL) measurements as constraints, the evaluation of the Moderate Resolution Imaging Spectroradiometer (MODIS) daily and 8-day snow-cover products is carried out using the Upper Rio Grande River Basin as a test site. A time series of the snow areal extent (SAE) of the Upper Rio Grande Basin is retrieved from the MODIS tile h09v05 covering the time period from February 2000 to June 2004 using an automatic Geographic Information System (GIS)-based algorithm developed for this study. Statistical analysis between the streamflow at Otowi (NM) station and the SAE retrieved from the two MODIS snow-cover products shows that there is a statistically significant correlation between the streamflow and SAE for both products. This relationship can be disturbed by heavy rainstorms in the later springtime, especially in May. Correlation analyses show that the MODIS 8-day product has a better correlation (r=−0.404) with streamflow and has less percentage of spurious snowmelt events in wintertime than the MODIS daily product (r=−0.300). Intercomparison of these two products, with the SNOTEL data sets as the ground truth, shows that (1) the MODIS 8-day product has higher classification accuracy for both snow and land; (2) the omission error of misclassifying snow as land is similar for both products, both are low; (3) the MODIS 8-day product has a slightly higher commission error of misclassifying land as snow than the MODIS daily product; and (4) the MODIS daily product has higher omission errors of misclassifying both snow and land as clouds. Clouds are the major cause for reduction of the overall accuracy of the MODIS daily product. Improvement in suppressing clouds in the 8-day product is obvious from this comparison study. The sacrifice is the temporal resolution that is reduced from 1 to 8 days. The significance of the results is that the 8-day product will be more useful in evaluating the streamflow response to the snow-cover extent changes, especially from the long-term point of view considering its lower temporal resolution than the daily product. For clear days, the MODIS daily algorithm works quite well or even better than the MODIS 8-day algorithm.     

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.     

Using variance analysis of multitemporal MODIS images for rice field mapping in Bali Province,Indonesia

Existing methods for rice field classification have some limitations due to the large variety of land covers attributed to rice fields. This study used temporal variance analysis of daily Moderate Resolution Imaging Spectroradiometer (MODIS) satellite images to discriminate rice fields from other land uses. The classification result was then compared with the reference data. Regression analysis showed that regency and district comparisons produced coefficients of determination (R ²) of 0.97490 and 0.92298, whereas the root mean square errors (RMSEs) were 1570.70 and 551.36 ha, respectively. The overall accuracy of the method in this study was 87.91%, with commission and omission errors of 35.45% and 17.68%, respectively. Kappa analysis showed strong agreement between the results of the analysis of the MODIS data using the method developed in this study and the reference data, with a kappa coefficient value of 0.8371. The results of this study indicated that the algorithm for variance analysis of multitemporal MODIS images could potentially be applied for rice field mapping.     

Integrated assessment on multi-temporal and multi-sensor combinations for reducing cloud obscuration of MODIS snow cover products of the Pacific Northwest USA

MODIS (Moderate Resolution Imaging Spectroradiometer) snow cover products, of daily, freely available, worldwide spatial extent at medium spatial resolution, have been widely applied in regional snow cover and modeling studies, although high cloud obscuration remains a concern in some applications. In this study, various approaches including daily combination, adjacent temporal deduction, fixed-day combination, flexible multi-day combination, and multi-sensor combination are assessed to remove cloud obscuration while still maintain the temporal and spatial resolutions. The performance of the resultant snow cover maps are quantitatively evaluated against in situ observations at 244 SNOTEL stations over the Pacific Northwest USA during the period of 2006-2008 hydrological years. Results indicate that daily Terra and Aqua MODIS combination and adjacent temporal deduction can reduce cloud obscuration and classification errors although an annual mean of 37% cloud coverage remains. Classification errors in snow-covered months are actually small and tend to underestimate the snow cover. Primary errors of MODIS daily, fixed and flexible multi-day combination products occur during transient months. Flexible multi-day combination is an efficient approach to maintain the balance between temporal resolution and realistic estimation of snow cover extent since it uses two thresholds to control the combination processes. Multi-sensor combinations (daily or multi-day), taking advantage of MODIS high spatial resolution and AMSR-E cloud penetration ability, provide cloud-free products but bring larger image underestimation errors as compared with their MODIS counterparts.     

Toward improved daily snow cover mapping with advanced combination of MODIS and AMSR-E measurements

Taking three snow seasons from November 1 to March 31 of year 2002 to 2005 in northern Xinjiang, China as an example, this study develops a new daily snow cover product (500 m) through combining MODIS daily snow cover data and AMSR-E daily snow water equivalent (SWE) data. By taking advantage of both high spatial resolution of optical data and cloud transparency of passive microwave data, the new daily snow cover product greatly complements the deficiency of MODIS product when cloud cover is present especially for snow cover product on a daily basis and effectively improves daily snow detection accuracy. In our example, the daily snow agreement of the new product with the in situ measurements at 20 stations is 75.4%, which is much higher than the 33.7% of the MODIS daily product in all weather conditions, even a little higher than the 71% of the MODIS 8-day product (cloud cover of ~ 5%). Our results also indicate that i) AMSR-E daily SWE imagery generally agrees with MOD10A1 data in detecting snow cover, with overall agreement of 93.4% and snow agreement of 96.6% in the study area; ii) AMSR-E daily SWE imagery underestimates the snow covered area (SCA) due to its coarse spatial resolution; iii) The new snow cover product can better and effectively capture daily SCA dynamics during the snow seasons, which plays a significant role in reduction, mitigation, and prevention of snow-caused disasters in pastoral areas.     

Analysis of MODIS imagery for detection of clear cuts in the boreal forest in north-west Russia

This study examines the feasibility of using MODIS images (MOD02 products) for the detection and monitoring of forest clear cuts in the boreal forest in north-west Russia. The proposed approach combines three change detection methods, including Change Vector Analysis, Textural Analysis using the coefficient of variation, and Constrained Energy Minimization analysis. For each individual method a series of thresholds was tested in order to obtain an optimal identification of clear cuts. A clear cut detection was only accepted if the change was detected by each individual method. All input parameters needed were derived from a set of reference clear cuts, mapped from 30 m resolution Landsat ETM+ imagery and used also for accuracy assessment. Change assessment was tested with MODIS images of two and of three acquisition dates. Referring to two test sites (Karelia, Komi) the detection omission and commission errors, assessed within a 3 × 3 pixels moving kernel, were at 23% and 8%, and at 21% and 17%, respectively. In terms of detectable clear cut size, a detection accuracy of about 90% can be expected for clear cuts in the size category above 15 ha, which contains the majority of cuts in the region. MODIS therefore provides good capabilities for large scale monitoring of major clear cut activities in the boreal forests of north-western Russia.     

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.     

An application of MODIS data to snow cover monitoring in a pastoral area: A case study in Northern Xinjiang, China

Snow is an important land cover on the earth's surface. It is characterized by its changing nature. Monitoring snow cover extent plays a significant role in dynamic studies and prevention of snow-caused disasters in pastoral areas. Using NASA EOS Terra/MODIS snow cover products and in situ observation data during the four snow seasons from November 1 to March 31 of year 2001 to 2005 in northern Xinjiang area, the accuracy of MODIS snow cover mapping algorithm under varied snow depth and land cover types was analyzed. The overall accuracy of MODIS daily snow cover mapping algorithm in clear sky condition is high at 98.5%; snow agreement reaches 98.2%, and ranges from 77.8% to 100% over the 4-year period for individual sites. Snow depth (SD) is one of the major factors affecting the accuracy of MODIS snow cover maps. MODIS does not identify any snow for SD less than 0.5 cm. The overall accuracy increases with snow depth if SD is equal to or greater than 3 cm, and decreases for SD below 3 cm. Land cover has an important influence in the accuracy of MODIS snow cover maps. The use of MOD10A1 snow cover products is severely affected by cloud cover. The 8-day composite products of MOD10A2 can effectively minimize the effect of cloud cover in most cases. Cloud cover in excess of 10% occurs on 99% of the MOD10A1 products and 14.7% of the MOD10A2 products analyzed during the four snow seasons. User-defined multiple day composite images based on MOD10A1, with flexibilities of selecting composite period, starting and ending date and composite sequence of MOD10A1 products, have an advantage in effectively monitoring snow cover extent for regional snow-caused disasters in pastoral areas.     

Sensitivity of vegetation phenology detection to the temporal resolution of satellite data

Vegetation phenology derived from satellite data has increasingly received attention for applications in environmental monitoring and modelling. The accuracy of phenological estimates, however, is unknown at the regional and global level because field validation data are insufficient. To assess the accuracy of satellite‐derived phenology, this study investigates the sensitivity of phenology detection to both the temporal resolution of sampling and the number of consecutive missing values (usually representing cloud cover) in the time series of satellite data. To do this, time series of daily vegetation index data for various ecosystems are modelled and simulated using data from Moderate‐Resolution Imaging Spectroradiometer (MODIS) data. The annual temporal data are then fitted using piecewise logistic functions, which are employed to calculate curvature change rate for detecting phenological transition dates. The results show that vegetation phenology can be estimated with a high precision from time series with temporal resolutions of 6–16 days even if daily data contains some uncertainties. If the temporal resolution is no coarser than 16 days for time series sampled using an average composite, the absolute errors are less than 3 days. On the other hand, the phase shift of temporal sampling is shown to have limited impacts on phenology detection. However, the accuracy of phenology detection may be reduced greatly if missing values in the time series of 16‐day MODIS data occur around the onsets of phenological transition dates. Even so, the probability that the absolute error in phenological estimates is greater than 5 days is less than 4% when only one period is missing in the time series of 16‐day data during vegetation growing seasons; this probability increases to 20% if there are two consecutive missing values.     

Sensitivity of near infrared total water vapour estimate to calibration errors

Analysis of satellite data to estimate the precipitable water (also called the columnar water vapour) amount often leads to systematic errors in deduced precipitable water (PW). The causes of systematic errors are likely to be instrumental calibration errors rather than variability of atmospheric parameters. We use the MODTRAN 4.0 radiative transfer code to model effects of various calibration errors on the Multi-spectral Thermal Imager (MTI) daytime total water vapour estimate. From the considered sources of calibration errors (spectral band centre error, spectral bandwidth error and radiometric calibration error) the radiometric calibration error has the largest influence on the accuracy of total water vapour estimate. When the radiometric calibration error between 1% and 5% is combined with the estimated spectral band centre error of 1 nm and the bandwidth error of 0.5 nm, the total systematic error of the columnar water vapour estimate is expected to be between 8% and 26%. The accuracy of the retrieved PW using the MTI imagery over the NASA Stennis site and Oklahoma DOE (Department of Energy) ARM (Atmospheric Radiation Measurement program) site is about 17%, well within the estimated range due to calibration errors. A similarity between the MTI and the MODIS (Moderate Resolution Imaging Spectral-Radiometer) bands used for water vapour estimate suggests that a similar error analysis may be valid for the MODIS sensor. However, the narrow band instruments (with bandwidth around 10 nm) are much more sensitive to the band centre calibration error.     

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.     

Characterization of seasonal variation of forest canopy in a temperate deciduous broadleaf forest, using daily MODIS data

In this paper, we present an improved procedure for collecting no or little atmosphere- and snow-contaminated observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor. The resultant time series of daily MODIS data of a temperate deciduous broadleaf forest (the Bartlett Experimental Forest) in 2004 show strong seasonal dynamics of surface reflectance of green, near infrared and shortwave infrared bands, and clearly delineate leaf phenology and length of plant growing season. We also estimate the fractions of photosynthetically active radiation (PAR) absorbed by vegetation canopy (FAPAR_canopy), leaf (FAPAR_leaf), and chlorophyll (FAPAR_chl), respectively, using a coupled leaf-canopy radiative transfer model (PROSAIL-2) and daily MODIS data. The Markov Chain Monte Carlo (MCMC) method (the Metropolis algorithm) is used for model inversion, which provides probability distributions of the retrieved variables. A two-step procedure is used to estimate the fractions of absorbed PAR: (1) to retrieve biophysical and biochemical variables from MODIS images using the PROSAIL-2 model; and (2) to calculate the fractions with the estimated model variables from the first step. Inversion and forward simulations of the PROSAIL-2 model are carried out for the temperate deciduous broadleaf forest during day of year (DOY) 184 to 201 in 2005. The reproduced reflectance values from the PROSAIL-2 model agree well with the observed MODIS reflectance for the five spectral bands (green, red, NIR₁, NIR₂, and SWIR₁). The estimated leaf area index, leaf dry matter, leaf chlorophyll content and FAPAR_canopy values are close to field measurements at the site. The results also showed significant differences between FAPAR_canopy and FAPAR_chl at the site. Our results show that MODIS imagery provides important information on biophysical and biochemical variables at both leaf and canopy levels.     

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.     

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.     

Regional mapping of gross light-use efficiency using MODIS spectral indices

Direct estimation of photosynthetic light-use efficiency (LUE) from space would be of significant benefit to LUE-based models which use inputs from remote sensing to estimate terrestrial productivity. The Photochemical Reflectance Index (PRI) has shown promise in tracking LUE at the leaf- to small canopy levels, but its use at regional to global scales still remains a challenge. In this study, we used different formulations of PRI calculated from the MODIS ocean band centered at 531 nm and a set of alternative reference bands at 488, 551, and 678 nm to explore the relationship between PRI and LUE where LUE was measured at eight eddy covariance flux towers located in the boreal forest of Saskatchewan, Canada. The magnitude and variability of LUE was significantly lower at the times when useful MODIS ocean band images were available (i.e. around midday under clear-sky conditions) relative to the rest of the growing season. PRI₆₇₈ (reference band at 678 nm) showed the strongest relationship (r² = 0.70) with LUE_90a (i.e. 90-minute mean LUE calculated using Absorbed Photosynthetically Active Radiation, APAR), but only when all sites were combined. Overall, the relationships between the MODIS PRIs and LUE_90a were always stronger for observations closer to the backscatter direction and there were no significant differences in the strength of the correlations whether LUE was calculated based on incident PAR or on APAR. Predictions of ecosystem photosynthesis at the time of the MODIS overpasses were significantly improved by multiplying either PAR or APAR by MODIS PRI (r² improved from 0.09 to 0.44 and 0.54 depending on the PRI formulation).We used our PRI-LUE model to create a regional LUE_90a map for the three cover types covering 47,500 km² around the flux sites. The MODIS PRI-derived LUE_90a map appeared to capture more realistic spatial heterogeneity of LUE across the landscape compared to a daily LUE map derived using the look-up table in the MODIS GPP (MOD17) algorithm. While our LUE map is only a snapshot of minimum regional LUE_90a values, with appropriate gap-filling methods it could be used to improve regional-scale monitoring of GPP. Moreover, the strong relationship between midday and daily LUE on clear days (r² = 0.93) indicates that instantaneous MODIS observations of LUE_90a could be used to estimate daily LUE. Finally, pixel shadow fraction from the 5-Scale geometric-optical model was closely related to both MODIS PRI and tower-derived LUE suggesting that differences in stand leaf area and in diffuse illumination among flux sites play a role in the relationship we observed between LUE and MODIS PRI.     

Field work and statistical analyses for enhanced interpretation of satellite fire data

Because their broad spatial and temporal coverage, satellites provide the main source of fire data for Amazonia. A key to the application of these tools for environmental studies is the appropriate interpretation of the data they provide. To enhance the interpretation of satellite fire data for this region, we collected ground-based data on fires in 2001 and 2002 using a simple and passive method, and statistically related these data to corresponding estimates from AVHRR and MODIS fire products using error matrices. Multiple methods of analyses from simple to complex produced qualitatively similar results. Total accuracies for both fire products were very high (> 99%) and dominated by accurate (> 99%) non-fire detection. Kappa statistics and fire-detection accuracies were substantially lower, with omission errors higher than commission errors. Results calculated using several different sets of spatial-matching parameters of analysis showed that Kappa was 1-10.6% for AVHRR, and 0-1.4% for MODIS. User's accuracy for fires was 0-40% for AVHRR and 3-100% for MODIS. Producer's accuracy for fires was 0-8% for AVHRR and 0-1% for MODIS. Statistical evaluations of potential explanatory factors showed that fire size and sampling time were dominant factors for low accuracies. Results from this study indicate that current satellite fire products are providing a limited sample of the fire activity in the region, and that ground-based analyses can substantially contribute to the interpretation of these products.