利用观测资料结合遥感数据,基于半阶导数法估算青藏高原多年冻土区唐古拉和西大滩两个观测点2010—2012年的地表土壤热通量,并与“观测值”(利用10 cm土壤热通量和5 cm土壤温度观测数据计算得到的地表土壤热通量)及已有的3种遥感估算方法(GLEAM、PM-RS和组合法)进行对比分析。结果表明:1)半阶导数法估算的地表土壤热通量(G0)与观测G0具有较好的一致性,估算精度高;2)与已有的GLEAM、PM-RS和组合法相比,半阶导数法表现更好,估算的两个站点G0的相关系数(R)分别为0.85和0.81,偏差(Bias)分别为-0.70和0.08 W/m2,均方根误差(RMSE)分别为5.78和8.79 W/m2,平均绝对误差(MAE)分别为4.22和6.81 W/m2。相较于其他方法,半阶导数法的优势在于物理机制明确,参数设置简洁,仅需单层地表温度值和土壤热属性数据作为模型输入。