This study evaluated the seasonal productivity of a prairie grassland (Mattheis Ranch, in Alberta, Canada) using a combination of remote sensing, eddy covariance, and field sampling collected in 2012-2013. A primary objective was to evaluate different ways of parameterizing the light-use efficiency (LUE) model for assessing net ecosystem fluxes at two sites with contrasting productivity. Three variations on the NDVI (Normalized Difference Vegetation Index), differing by formula and footprint, were derived: (1) a narrow-band NDVI (NDVI 680,800 , derived from mobile field spectrometer readings); (2) a broad-band proxy NDVI (derived from an automated optical phenology station consisting of broad-band radiometers); and (3) a satellite NDVI (derived from MODIS AQUA and TERRA sensors). Harvested biomass, net CO 2 flux, and NDVI values were compared to provide a basis for assessing seasonal ecosystem productivity and gap filling of tower flux data. All three NDVIs provided good estimates of dry green biomass and were able to clearly show seasonal changes in vegetation growth and senescence, confirming their utility as metrics of productivity. When relating fluxes and optical measurements, temporal aggregation periods were considered to determine the impact of aggregation on model accuracy. NDVI values from the different methods were also calibrated against f APAR green (the fraction of photosynthetically active radiation absorbed by green vegetation) values to parameterize the APAR green (absorbed PAR) term of the LUE (light use efficiency) model for comparison with measured fluxes. While efficiency was assumed to be constant in the model, this analysis revealed hysteresis in the seasonal relationships between fluxes and optical measurements, suggesting a slight change in efficiency between the first and second half of the growing season. Consequently, the best results were obtained by splitting the data into two stages, a greening phase and a senescence phase, and applying separate fits to these two periods. By incorporating the dynamic irradiance regime, the model based on APAR green rather than NDVI best captured the high variability of the fluxes and provided a more realistic depiction of missing fluxes. The strong correlations between these optical measurements and independently measured fluxes demonstrate the utility of integrating optical with flux measurements for gap filling, and provide a foundation for using remote sensing to extrapolate from the flux tower to larger regions (upscaling) for regional analysis of net carbon uptake by grassland ecosystems.