AbstractIn different studies, land surface temperature (LST) observations have been assimilated into the variational data assimilation (VDA) approaches to estimate turbulent heat fluxes. The VDA methods yield accurate turbulent heat fluxes, but they need an adjoint model, which is difficult to derive and code. They also cannot directly calculate the uncertainty of their estimates. To overcome the abovementioned drawbacks, this study assimilates LST data from Geostationary Operational Environmental Satellite into the ensemble Kalman smoother (EnKS) data assimilation system to estimate turbulent heat fluxes. EnKS does not need to derive the adjoint term and directly generates statistical information on the accuracy of its predictions. It uses the heat diffusion equation to simulate LST. EnKS with the state augmentation approach finds the optimal values for the unknown parameters (i.e., evaporative fraction and neutral bulk heat transfer coefficient, CHN) by minimizing the misfit between LST observations from Geostationary Operational Environmental Satellite and LST estimations from the heat diffusion equation. The augmented EnKS scheme is tested over six Ameriflux sites with a wide range of hydrological and vegetative conditions. The results show that EnKS can predict not only the model parameters and turbulent heat fluxes but also their uncertainties over a variety of land surface conditions. Compared to the variational method, EnKS yields suboptimal turbulent heat fluxes. However, suboptimality of EnKS is small, and its results are comparable to those of the VDA method. Overall, EnKS is a feasible and reliable method for estimation of turbulent heat fluxes.