The objective of this paper was to evaluate the accuracy of two advanced blending algorithms, Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) and Enhanced Spatial and Temporal Adaptive Reflectance Fusion Model (ESTARFM) to downscale Moderate Resolution Imaging Spectroradiometer (MODIS) indices to the spatial resolution of Landsat. We tested two approaches: (i) “Index-then-Blend” (IB); and (ii) “Blend-then-Index” (BI) when simulating nine indices, which are widely used for vegetation studies, environmental moisture assessment and standing water identification. Landsat-like indices, generated using both IB and BI, were simulated on 45 dates in total from three sites. The outputs were then compared with indices calculated from observed Landsat data and pixel-to-pixel accuracy of each simulation was assessed by calculating the: (i) bias; (ii) R²; and (iii) Root Mean Square Deviation (RMSD). The IB approach produced higher accuracies than the BI approach for both blending algorithms for all nine indices at all three sites. We also found that the relative performance of the STARFM and ESTARFM algorithms depended on the spatial and temporal variances of the Landsat-MODIS input indices. Our study suggests that the IB approach should be implemented for blending of environmental indices, as it was: (i) less computationally expensive due to blending single indices rather than multiple bands; (ii) more accurate due to less error propagation; and (iii) less sensitive to the choice of algorithm.