“…Variational data assimilation (VDA) methods estimate the key unknowns of the surface energy balance (SEB) equations (i.e., neutral bulk heat transfer coefficient, C HN , and evaporative fraction, EF) by assimilating LST observations into the heat diffusion or force-restore equations (Abdolghafoorian et al, 2017;Bateni, Entekhabi, & Jeng, 2013;Bateni, Entekhabi, & Castelli, 2013;Bateni et al, 2014;Bateni & Entekhabi, 2012a;Bateni & Liang, 2012;Crow & Kustas, 2005;He et al, 2018;Qin et al, 2007;Sini et al, 2008;Xu et al, 2014Xu et al, , 2016Xu, Bateni, et al, 2015;Xu, He, et al, 2019). The unknown parameters of the VDA approaches (i.e., C HN and EF) are obtained by minimizing the misfit between the LST observations and estimations from the force-restore or heat diffusion equation.…”