基于主动微波观测资料土壤含水量廓线估计数据同化.pptxVIP

内 容引言陆面数据同化方法非饱和土壤水流问题的预报算子土壤湿度数据同化方案及试验小结与讨论一、引言土壤湿度通过影响陆面和大气之间的水分和能量的平衡，对气候、区域环境变化起着非常重要的作用，受到大气科学、土壤学、农业工程、环境工程和地下水动力学等领域的高度重视；获取土壤湿度的途径： 场地点观测、陆面水文模式模拟和遥感反演来估计。 中国土壤湿度观测的情况以农业气象观测为基础的每旬逢8观测，逢1发报，目前约有600个站；以农气加密观测为基础，逢3观测，逢6发报，约有400个站；2005年新布的100个自动观测站，目前刚安装运行，发报站不到10个，资料尚未经过对比分析。二、陆面数据同化方法直接插入法直接插入法是一种最简单的方法，就是在观测数据点上简单地用观测值替代模型的状态场而不考虑观测误差。该方法的前提假设是认为观测数据绝对准确。 2、Kalman滤波方法考虑如下的随机非线性离散系统： 式中 和 是非线性函数，向量 和 为系统过程噪声和观测噪声序列，假设为零均值白噪声序列。 Kalman滤波的预估—校正算法示意图预报部分（时间更新）1、状态预报 2、误差预报 分析部分（观测更新）1、计算增益矩阵2、对变量的分析3、对分析变量的误差的分析集合Kalman滤波计算流程图集合Kalman滤波是用一个集合统计的方法来估计Kalman滤波方程组中的分析误差协方差矩阵和背景场误差协方差矩阵；主要思路：先根据背景场和观测值的特征误差分布来对背景场和观测值加以一系列的扰动，然后用这些加上不同扰动的背景场和观测场进行分析，得到一组分析值。然后用这组分析值的差异作为分析误差的统计样本来进行分析误差协方差的估计。对这组分析值作一个短期预报后，也可以得到一组预报值。同样，把这组预报值的差异作为背景误差的统计样本来进行背景误差协方差的估计。三、非饱和土壤水流问题的预报算子连续预报算子 为地表通量， 向下为正，原点取在地表， 为土壤厚度。离散的状态方程土壤参数：四、土壤湿度数据同化方案及试验（一）、理想同化数值试验方案 数值试验设计为孪生试验，即观测值由模式生成；初估值的误差方差阵取为对角矩阵，即各层间不相关，且对角线上元素给定0.5的标准差；模式误差方差阵取为对角矩阵，且以半小时模拟时间的系统状态的5%作为标准差；观测误差方差阵取为对角矩阵，且以观测值的2%作为标准差． 由于是做理想试验,故为了简化和集中于算法研究,这里不使用辐射亮温估计模型,而直接将”真实”廓线上的值作为相应位置的观测值,从而观测矩阵为 Soil parametersRetrieval conditionSoil typeClay loamDepth100cmSoil moisture at saturation0.476Number of nodes50Hydraulic conductivity at saturation0.00025cm/sBoundary condition0.5cm/dayExponent b8.52Initial condition0.40Matric potential at saturation-63.0cmBad initial guess0.35扩展Kalman滤波同化试验模拟（点），观测深度2cm的同化(实线)，观测深度6cm的同化(点线)，观测深度10cm的同化(虚线)，观测（点虚线）．每3小时引入一次资料模拟（点），观测深度20cm的同化(实线)，观测深度50cm的同化(点线)，观测深度80cm的同化(虚线)，观测（点虚线）．每3小时引入一次资料模拟（点），观测深度2、6cm的同化(实线)，观测深度2、80cm的同化(点线)，观测深度50、80cm的同化(虚线)，观测（点虚线）．每3小时引入一次资料观测深度2、6、10cm的同化(实线)，观测深度2、6、10、20cm的同化(点线)，观测深度2、6、10、20、50cm的同化(虚线)，观测（点虚线）．每3小时引入一次资料每24小时引入一次资料每24小时引入一次资料直接插入法同化试验模拟（点），观测深度2cm的同化(实线)，观测深度6cm的同化(点线)，观测深度10cm的同化(虚线)，观测（点虚线）．每3小时引入一次资料 模拟（点），观测深度20cm的同化(实线)，观测深度50cm的同化(点线)，观测深度80cm的同化(虚线)，观测（点虚线）．每3小时引入一次资料模拟（点），观测深度2、6cm的同化(实线)，观测深度2、80cm的同化(点线)，观测深度50、80cm的同化(虚线)，观测（点虚线）．每3小时引入一次资料 模拟（点虚线），观测深度2、6、10cm的同化(实线)，观测深度2、6、10、20cm的同化(点线)，观测深度2、6、10、20、50cm的同化(虚线)，观测（点）