Tea plant (Camellia sinensis) is a major global crop consumed as a drink after water. Quantifying carbon flux, specifically the net ecosystem exchange (NEE), in tea plantations is essential for determining carbon sequestration and ecosystem carbon balance. The Eddy covariance (EC) system is widely used for continuous monitoring of carbon flux but high costs associated with installation and maintenance limit its widespread adoption. In addition, EC flux data is often discarded due to malfunction of instruments caused by adverse weather conditions. Therefore, additional approaches for estimating NEE are necessary to overcome these challenges and ensure accurate NEE measurement. For this purpose, three standalone tree-based machine learning (ML) models were used for NEE estimation using EC flux data collected from tea ecosystem located in subtropical region (Danyang county of Zhenjiang city) of China. To address the accuracy limitations inherent in standalone ML models, the ensemble mechanism based on voting regressor method was proposed. In addition, k-fold cross-validation based on early stopping process was also used to enhance the performance of standalone ML models. Based on visual plots (scatter diagram, heatMap, Taylor diagram) and performing indices (root-mean-square error (RMSE), determination coefficient (R2), mean absolute error (MAE), mean absolute percentage error (MAPE), Nash–Sutcliffe efficiency (NSE), correlation coefficient (r), Kling Gupta Efficiency (KGE) and index of agreement (d)), the findings indicated that non-linear ensemble-generalized regression neural network (NLE-GRNN) significantly improved standalone ML model's results. In current study, the highest NSE, r and d in case of standalone ML model (DT) achieved 0.49, 0.73 and 0.75 respectively while our proposed NLE-GRNN model improved 48 % in NSE value (NSE = 0.97), 25 % in r value (r = 0.98) and 24 % in d value (d = 0.99). Likewise, NLE-GRNN significantly reduce errors (MAE, MAPE and RMSE) and provides NEE estimate closet to the observed value. The impact of climatic variables on NEE using shapley additive explanations (SHAP) analysis revealed that Rg (solar radiation) and Tair (air temperature) were the prime factors controlling NEE variation in the tea ecosystem. Considering the high accuracy and stability of the studied ML models, it is recommended to apply developed ensemble ML model (NLE-GRNN) for significant improvement of NEE estimate in the tea biomes or other ecosystems.