Long-term global land surface fractional vegetation cover (FVC) products are essential for various applications. Currently, several global FVC products have been generated from medium spatial resolution remote sensing data. However, validation results indicate that there are inconsistencies and spatial and temporal discontinuities in the current FVC products. Therefore, the Global LAnd Surface Satellite (GLASS) FVC product algorithm using general regression neural networks (GRNNs), which achieves an FVC estimation accuracy comparable to that of the GEOV1 FVC product with much improved spatial and temporal continuities, was developed. However, the computational efficiency of the GRNNs method is low and unsatisfactory for generating the long-term GLASS FVC product. Therefore, the objective of this study was to discover an alternative algorithm for generating the GLASS FVC product that has both an accuracy comparable to that of the GRNNs method and adequate computational efficiency. Four commonly used machine learning methods, back-propagation neural networks (BPNNs), GRNNs, support vector regression (SVR), and multivariate adaptive regression splines (MARS), were evaluated. After comparing its performance of training accuracy and computational efficiency with the other three methods, the MARS model was preliminarily selected as the most suitable algorithm for generating the GLASS FVC product. Direct validation results indicated that the performance of the MARS model (R 2 = 0.836, RMSE = 0.1488) was comparable to that of the GRNNs method (R 2 = 0.8353, RMSE = 0.1495), and the global land surface FVC generated from the MARS model had good spatial and temporal consistency with that generated from the GRNNs method. Furthermore, the computational efficiency of MARS was much higher than that of the GRNNs method. Therefore, the MARS model is a suitable algorithm for generating the GLASS FVC product from Moderate Resolution Imaging Spectroradiometer (MODIS) data.