_Article_  
**Comparison of Machine Learning Regression Algorithms for Cotton Leaf Area Index Retrieval Using Sentinel-2 Spectral Bands**  
**[PERSON] 1,2**, **[PERSON] 1,*, [PERSON] 1,2, [PERSON] 1,3 and [PERSON] 1,2**  
1 Key Laboratory of Digital Earth Sciences, Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, Beijing 100101, China; [EMAIL_ADDRESS] (H.M.); [EMAIL_ADDRESS] (F.J.); [EMAIL_ADDRESS] (Q.Z.); [EMAIL_ADDRESS] (H.F)  
2 University of Chinese Academy of Sciences, Beijing 100049, China  
3 College of Geomatics and Geoinformation, Guilin University of Technology, Guilin 541004, China  
* Correspondence: [EMAIL_ADDRESS]; Tel.: +86-010-6486-9473  
**Abstract:** Leaf area index (LAI) is a crucial crop biophysical parameter that has been widely used in a variety of fields. Five state-of-the-art machine learning regression algorithms (MLRAs), namely, artificial neural network (ANN), support vector regression (SVR), Gaussian process regression (GPR), random forest (RF) and gradient boosting regression tree (GBRT), have been used in the retrieval of cotton LAI with Sentinel-2 spectral bands. The performances of the five machine learning models are compared for better applications of MLRAs in remote sensing, since challenging problems remain in the selection of MLRAs for crop LAI retrieval, as well as the decision as to the optimal number for the training sample size and spectral bands to different MLRAs. A comprehensive evaluation was employed with respect to model accuracy, computational efficiency, sensitivity to training sample size and sensitivity to spectral bands. We conducted the comparison of five MLRAs in an agricultural area of Northwest China over three cotton seasons with the corresponding field campaigns for modeling and validation. Results show that the GBRT model outperforms the other models with respect to model accuracy in average (\(\overline{R^{2}}\) = 0.854, \(\overline{RMSE}\) = 0.674 and \(\overline{MAE}\) = 0.456). SVR achieves the best performance in computational efficiency, which means it is fast to train, and to validate that it has great potentials to deliver near-real-time operational products for crop management. As for sensitivity to training sample size, GBRT behaves as the most robust model, and provides the best model accuracy on the average among the variations of training sample size, compared with other models (\(\overline{R^{2}}\) = 0.884, \(\overline{RMSE}\) = 0.615 and \(\overline{MAE}\) = 0.452). Spectral bands sensitivity analysis with dCor (distance correlation), combined with the backward elimination approach, indicates that SVR, GPR and RF provide relatively robust performance to the spectral bands, while ANN outperforms the other models in terms of model accuracy on the average among the reduction of spectral bands (\(\overline{R^{2}}\) = 0.881, \(\overline{RMSE}\) = 0.625 and \(\overline{MAE}\) = 0.480). A comprehensive evaluation indicates that GBRT is an appealing alternative for cotton LAI retrieval, except for its computational efficiency. Despite the different performance of the ML models, all models exhibited considerable potential for cotton LAI retrieval, which could offer accurate crop parameters information timely and accurately for crop fields management and agricultural production decisions.**  
**Keywords:** leaf area index (LAI); machine learning; Sentinel-2; sensitivity analysis; training sample size; spectral bands