陆面过程是地球表层系统中的关键过程，准确地估计陆表状态对提高地球系统模式性能具有重要的作用。陆面数据同化系统通过有机融合多源对地观测数据和陆面模式的优势，提高了陆表状态变量的可观测性与可预测性。由于陆面模式中参数存在较大不确定性，且同化系统中缺乏有效的模型误差和观测误差的估计方法，这些问题阻碍了同化系统的进一步发展和推广。因此，亟需在数据同化系统中实现高效且协调地同时估计多类参数和误差。

针对上述问题，清华大学地学系卢麾教授课题组此前发展了基于集合卡尔曼滤波（EnKF）的双循环同化算法，解决了误差参数、模型参数、观测算子参数同时估算的难题。但是，双循环算法原型参数估计效率较低，限制了它在大尺度同化研究中的应用潜力。为解决这一问题，研究组进一步基于自动微分的扩展卡尔曼滤波（EKF-AD），通过降低模型集合在参数估计中的成本，实现了数据同化系统多参数的高效估计。为验证该算法，研究设计了两组实验，实验1单点同化实验（DA_EKF和DA_EnKF），实验2青藏高原区域的同化实验（DA_TP）和模拟实验（OL_TP）。

图1 基于EKF和EnKF的双循环LDAS参数优化时间对比。

图2. Shiquanhe（a）、Naqu（b）、Maqu（c）和Pali（d）区域土壤水分模拟值、估计值、和观测值的对比散点图。

图3. DA_TP（a）和OL_TP（b）生产的蒸散发与PML_V2蒸散发产品的对比散点图。

研究结果表明，基于EKF-AD的双循环同化算法，参数估计时间缩短近60倍（图1）。在青藏高原同化实验中，土壤水分的估计误差明显小于默认参数模拟实验（图2）；此外，通过验证发现同化土壤水分对蒸散发的模拟也有一定的积极作用（图3）。

图4 沙土含量((a1)和(b1))、黏土含量((a2)和(b2))以及孔隙度((a3)和(b3))默认值(左图)与估计值(右图)的空间分布图。

图5 在（a）Naqu和（b）Maqu验证网区域上的SMAP、SMAP_corr与站点观测的土壤水分时间序列对比图和散点图。

双循环算法可同时估计模型参数、观测算子参数和误差参数。通过对比分析，发现相较于陆面模式默认的土壤参数，同化系统估计的土壤参数更为合理（图4）；此外，观测算子参数的估计值可以直接用于校正SMAP土壤水分产品的系统性偏差（图5）。

研究基于EKF-AD的双循环算法在青藏高原同化实验中展现了良好的同化性能，为水文、生态、陆面数据同化提供了新的思路。相关研究以“Quick estimation of parameters for the land surface data assimilation system and its influence based on the extended Kalman filter and automatic differentiation”（基于扩展卡尔曼滤波和自动微分技术对陆面数据同化系统参数的快速估计及其影响）为题发表于Science China Earth Sciences。此外，课题组还利用双循环数据同化系统实现了通过同化卫星地表温度观测改进土壤水分模拟。研究成果以“Improving Surface Soil Moisture Estimation Through Assimilating Satellite Land Surface Temperature With a Linear SM-LST Relationship”为题发表于IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing。

清华大学地学系博士后田佳鑫为论文第一作者，卢麾教授为论文通讯作者。合作作者来自于清华大学地学系、中科院地理科学与资源研究所、西南大学。研究得到了第二次青藏高原综合科学考察研究（2019QZKK0206）、国家重点研发计划项目（2022YFC3002901）和国家自然科学基金项目（42271491）支持。

论文信息：

Tian, J., Lu, H., Yang, K. et al. Quick estimation of parameters for the land surface data assimilation system and its influence based on the extended Kalman filter and automatic differentiation. Science China Earth Sciences. 66, 2546–2562 (2023). https://doi.org/10.1007/s11430-022-1180-8

Tian, J., Lu, H., Yang, K. et al. Improving Surface Soil Moisture Estimation Through Assimilating Satellite Land Surface Temperature with a Linear SM-LST Relationship. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing. 16，7777-7790 (2023). http://doi.org/10.1109/JSTARS.2023.3305888

供稿：田佳鑫

编辑：王佳音

审核：卢麾 张强