基于非饱和土壤水流模型及地面点观测的土壤湿度数据同化方案39

1、基于非饱和土壤水流模型及地面点观测的土壤湿度数据同化方案 张生雷 谢正辉 中国科学院大气物理研究所2006年8月10日2006年全球变化科学研究生暑期学习交流报告内 容引言陆面数据同化方法非饱和土壤水流问题的预报算子土壤湿度数据同化方案及试验小结与讨论一、引言土壤湿度通过影响陆面和大气之间的水分和能量的平衡，对气候、区域环境变化起着非常重要的作用，受到大气科学、土壤学、农业工程、环境工程和地下水动力学等领域的高度重视；获取土壤湿度的途径： 场地点观测、陆面水文模式模拟和遥感反演来估计。 中国土壤湿度观测的情况以农业气象观测为根底的每旬逢8观测，逢1发报，目前约有600个站；以农气加密观测为根底

2、，逢3观测，逢6发报，约有400个站；2005年新布的100个自动观测站，目前刚安装运行，发报站不到10个，资料尚未经过比照分析。二、陆面数据同化方法直接插入法直接插入法是一种最简单的方法，就是在观测数据点上简单地用观测值替代模型的状态场而不考虑观测误差。该方法的前提假设是认为观测数据绝对准确。 2、Kalman滤波方法考虑如下的随机非线性离散系统： 式中 和 是非线性函数，向量 和 为系统过程噪声和观测噪声序列，假设为零均值白噪声序列。 Kalman滤波的预估校正算法示意图预报局部时间更新1、状态预报 2、误差预报 分析局部观测更新1、计算增益矩阵2、对变量的分析3、对分析变量的误差的分析集

3、合Kalman滤波计算流程图集合Kalman滤波是用一个集合统计的方法来估计Kalman滤波方程组中的分析误差协方差矩阵和背景场误差协方差矩阵；主要思路：先根据背景场和观测值的特征误差分布来对背景场和观测值加以一系列的扰动，然后用这些加上不同扰动的背景场和观测场进行分析，得到一组分析值。然后用这组分析值的差异作为分析误差的统计样本来进行分析误差协方差的估计。对这组分析值作一个短期预报后，也可以得到一组预报值。同样，把这组预报值的差异作为背景误差的统计样本来进行背景误差协方差的估计。三、非饱和土壤水流问题的预报算子连续预报算子 为地表通量， 向下为正，原点取在地表， 为土壤厚度。离散的状态方程土

4、壤参数：四、土壤湿度数据同化方案及试验一、理想同化数值试验方案 数值试验设计为孪生试验，即观测值由模式生成；初估值的误差方差阵取为对角矩阵，即各层间不相关，且对角线上元素给定0.5的标准差；模式误差方差阵取为对角矩阵，且以半小时模拟时间的系统状态的5%作为标准差；观测误差方差阵取为对角矩阵，且以观测值的2%作为标准差 由于是做理想试验,故为了简化和集中于算法研究,这里不使用辐射亮温估计模型,而直接将真实廓线上的值作为相应位置的观测值,从而观测矩阵为 Soil parametersRetrieval conditionSoil typeClay loamDepth100cmSoil moistu

5、re 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的同化(

6、实线)，观测深度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小时引入

7、一次资料 模拟点，观测深度20cm的同化(实线)，观测深度50cm的同化(点线)，观测深度80cm的同化(虚线)，观测点虚线每3小时引入一次资料模拟点，观测深度2、6cm的同化(实线)，观测深度2、80cm的同化(点线)，观测深度50、80cm的同化(虚线)，观测点虚线每3小时引入一次资料 模拟点虚线，观测深度2、6、10cm的同化(实线)，观测深度2、6、10、20cm的同化(点线)，观测深度2、6、10、20、50cm的同化(虚线)，观测点每3小时引入一次资料模拟点，观测深度2cm的同化(实线)，观测深度6cm的同化(点线)，观测深度10cm的同化(虚线)，观测点虚线每24小时引入一次资料

8、模拟点，观测深度20cm的同化(实线)，观测深度50cm的同化(点线)，观测深度80cm的同化(虚线)，观测点虚线每24小时引入一次资料二、同化常规土壤湿度场地点观测数据的数值试验方案 图1非饱和土壤水流模型图 、 、 和 分别表示降水、蒸发、地表径流和冠层截流合肥土壤湿度观测站经度：117.2 纬度：31.9时间：1986年1月1日1993年12月31日观测：土壤观测深度1米，共分11层，分别是05，510，1020，2030，3040，4050，5060，6070，7080，8090，90100cm，共296次观测。同化方法：扩展Kalman滤波 土壤湿度同化中使用的土壤参数土壤参数土壤类

9、型 粘壤土饱和土壤湿度0.467饱和导水率 0.000245cm/s指数 b8.02饱和土壤水势 -26.3cm土壤深度、层数3米 31层长沙马坡岭土壤湿度观测站经度： 113.1 纬度： 28.2时间：1985年1月1日1989年12月31日观测：土壤观测深度1米，共分11层，分别是05，510，1020，2030，3040，4050，5060，6070，7080，8090，90100cm，共180次观测；同化方法：集合Kalman滤波同化试验结论基于非饱和土壤水流模型，初步建立了一个利用kalman滤波方法同化近地表土壤湿度观测资料的陆面数据同化方案；同化对模拟结果有一定改进，说明该同化方

10、案具有一定的合理性；近地表观测对同化结果有较大影响；观测层数多少对深层同化结果有一定影响；加大观测频率,可以进一步改善数据同化的效果 五、讨论扩展Kalman滤波同化中模式误差、观测误差问题；集合Kalman滤波同化中样本数及样本的抽取问题；土壤参数及观测数据的准备；陆面数据同化方案中的模式；多种资料的利用。谢 谢 !9H5E2B+x(u$rZoWkThQeMbJ7G4D1z-w*t!qYmVjSgOdLaI6F3B0y)v%s#pXlUiQfNcK8H5E2A+x*u$rZnWkThPeMaJ7G4C1z-w&t!pYmVjRgOdL9I6F3B0y(v%s#oXlUiQfNbK8H5D2A

11、+x*u$qZnWkShPeMaJ7F4C1z)w&t!pYmUjRgOcL9I6E3B+y(v%r#oXlTiQeNbK8G5D2A-x*u$qZnVkShPdMaJ3B0y(v%s#oXlUiQfNbK8H5D2A+x*u$qZnWkShPeMaJ7F4C1z)w&t!pYmUjRgOcL9I6E3B+y(v%r#oXlTiQfNbK8G5D2A-x*u$qZnVkShPdMaJ7F4C0z)w&s!pYmUjRfOcL9H6E3B+y(u%r#oWlTiQeNbJ8G5D1A-x*t$qYnVkSgPdMaI7F4C0z)v&s!pXmUjRfOcK9H6E2B+y(u%rZoWlThQ

12、eNbJ8G4D1A-w*t$qYnVjSgPdLaI7F3C0y)v&s#pXmUiRfNcK9H5E2B+x(u%rZoWkThQeMbJ8G4D1z-w*t!qYnVjSgOdLaI6F3C0y)v%s#pXlUiRfNcK8H5E2A+x(u$rZnWkThPeMbJ7G4C1z-w&t!qYmVjRgOdL9I6F3B0y)v%s#oXlUiQfNcK8H5D2A+x*u$rZnWkShPeMaJ7G4C1z)w&t!pYmVjRgOcL9I6E3B0y(v%r#oXlTiQfNbK8G5D2A-x*u$qZnWkShPdMaJ7F4C1z)w&s!pYmUjRgOcL9H6E3B+

13、y(v%r#oWlTiQeNbK8G5D1A-x*t$qZnVkSgPdMaI7F4C0z)v&s!pXmUjRfOcL9H6E2B+y(u%r#oWlThQeNbJ8G5D1A-w*t$qYnVkSgPdLaI7F3C0z)v&s#pXmUiRfOcK9H5E2B+x(u%rZoWkThQeMbJ8G4D1z-w*t!qYnVjSgPdLaI6F3C0y)v&s#pXlUiRfNcK9H5E2A+x(u$rZoWkThPeMbJ7G4D1z-w&t!qYmVjSgOdL9I6F3B0y)v%s#oXlUiQfNcK8H5E2A+x*u$rZnWkThPeMaJ7G4C1z-w&t!pYmVj

14、RgOdL9I6E3B0y(v%s#oXlTiQfNbK8H5D2A-x*u$qZnWkShPdMaJ7F4C1z)w&s!pYmUjRgOcL9I6E3B+y(v%r#oXlTiQeNbK8G5D2A-x*t$qZnVkShPdMaI7F4C0z)w&s!pXmUjRfOcL9H6E2B+y(u%r#oWlThQeNbJ8G5D1A-x*t$qYnVkSgPdMaI7F3C0z)v&s!pXmUiRfOcK9H6E2B+x(u%rZoWlThQeMbJ8G4D1A-w*t!qYnVjSgPdLaI6F3C0%r#oWlTiQeNbJ8G5D1A-x*t$qYnVkSgPdMaI7F3C0z)

15、v&s!pXmUiRfOcK9H6E2B+x(u%rZoWlThQeMbJ8G4D1A-w*t!qYnVjSgPdLaI7F3C0y)v&s#pXmUiRfNcK9H5E2B+x(u$rZoWkThQeMbJ7G4D1z-w*t!qYmVjSgOdLaI6F3B0y)v%s#pXlUiQfNcK8H5E2A+x(u$rZnWkThPeMbJ7G4C1z-w&t!qYmVjRgOdL9I6F3B0y(v%s#oXlUiQfNbK8H5D2A+x*u$qZnWkShPeMaJ7F4C1z)w&t!pYmUjRgOcL9I6E3B0y(v%r#oXlTiQfNbK8G5D2A-x*u$qZnVkSh

16、PdMaJ7F4C0z)w&s!pYmUjRfOcL9H6E3B+y(u%r#oWlTiQeNbJ8G5D1A-x*t$qZnVkSgPdMaI7F4C0z)v&s!pXmUjRfOcK9H6E2B+y(u%rZoWlThQeNbJ8G4D1A-w*t$qYnVjSgPdLaI7F3C0y)v&s#pXmUiRfOcK9H5E2B+x(u%rZoWkThQeMbJ8G4D1z-w*t!qYnVjSgOdLaI6F3C0y)v%s#pXlUiRfNcK8H5E2A+x(u$rZnWkThPeMbJ3C0z)v&s#pXmUiRfOcK9H5E2B+x(u%rZoWkThQeMbJ8G4D1z-w

17、*t!qYnVjSgOdLaI6F3C0y)v%s#pXlUiRfNcK8H5E2A+x(u$rZoWkThPeMbJ7G4D1z-w&t!qYmVjSgOdL9I6F3B0y)v%s#oXlUiQfNcK8H5D2A+x*u$rZnWkShPeMaJ7G4C1z)w&t!pYmVjRgOcL9I6E3B0y(v%s#oXlTiQfNbK8H5D2A-x*u$qZnWkShPdMaJ7F4C1z)w&s!pYmUjRgOcL9H6E3B+y(v%r#oWlTiQeNbK8G5D1A-x*t$qZnVkShPdMaI7F4C0z)w&s!pXmUjRfOcL9H6E2B+y(u%r#oWlThQ

18、eNbJ8G5D1A-w*t$qYnVkSgPdLaI7F3C0z)v&s#pXmUiRfOcK9H6E2B+x(u%rZoWlThQeMbJ8G4D1A-w*t!qYnVjSgPdLaI6F3C0y)v&WlThQeNbJ8G5D1A-w*t$qYnVkSgPdLaI7F3C0z)v&s!pXmUiRfOcK9H6E2B+x(u%rZoWlThQeMbJ8G4D1A-w*t!qYnVjSgPdLaI6F3C0y)v&s#pXlUiRfNcK9H5E2A+x(u$rZoWkThQeMbJ7G4D1z-w*t!qYmVjSgOdLaI6F3B0y)v%s#pXlUiQfNcK8H5E2A+x*u

19、$rZnWkThPeMaJ7G4C1z-w&t!pYmVjRgOdL9I6E3B0y(v%s#oXlUiQfNbK8H5D2A+x*u$qZnWkShPeMaJ7F4C1z)w&t!pYmUjRgOcL9I6E3B+y(v%r#oXlTiQeNbK8G5D2A-x*t$qZnVkShPdMaJ7F4C0z)w&s!pYmUjRfOcL9H6E3B+y(u%r#oWhPeMaJ7F4C1z)w&t!pYmUjRgOcL9I6E3B+y(v%r#oXlTiQeNbK8G5D2A-x*u$qZnVkShPdMaJ7F4C0z)w&s!pYmUjRfOcL9H6E3B+y(u%r#oWlTiQeNbJ

20、8G5D1A-x*t$qYnVkSgPdMaI7F3C0z)v&s!pXmUjRfOcK9H6E2B+y(u%rZoWlThQeNbJ8G4D1A-w*t$qYnVjSgPdLaI7F3C0y)v&s#pXmUiRfNcK9H5E2B+x(u$rZoWkThQeMbJ8G4D1z-w*t!qYnVjSgOdLaI6F3C0y)v%s#pXlUiRfNcK8H5E2A+x(u$rZnWkThPeMbJ7G4C1z-w&t!qYmVjRgOdL9I6F3B0y(v%s#oXlUiQfNcK8H5D2A+x*u$rZnWkShPeMaJ7G4C1z)w&t!pYmVjRgOcL9I6E3B0y(v%

21、r#oXlTiQfNbK8G5D2A-x*u$qZnVkShPdMaJ7F4C)v%s#oXlUiQfNcK8H5D2A+x*u$rZnWkShPeMaJ7G4C1z)w&t!pYmVjRgOcL9I6E3B0y(v%r#oXlTiQfNbK8G5D2A-x*u$qZnWkShPdMaJ7F4C1z)w&s!pYmUjRgOcL9H6E3B+y(v%r#oWlTiQeNbK8G5D1A-x*t$qZnVkSgPdMaI7F4C0z)v&s!pXmUjRfOcL9H6E2B+y(u%r#oWlThQeNbJ8G5D1A-w*t$qYnVkSgPdLaI7F3C0z)v&s#pXmUiRfOcK9

22、H5E2B+x(u%rZoWkThQeMbJ8G4D1A-w*t!qYnVjSgPdLaI6F3C0y)v&s#pXlUiRfNcK9H5E2A+x(u$rZoWkThPeMbJ7G4D1z-w&t!qYmVjSgOdL9I6F3B0y)v%s#oXlUiQfNcK8H5E2A+x*u$rZnWkThPeMaJ7G4C1z-w&t!pYmVjRgOdL9I6E3B0y(v%s#oXlTiQfNG4D1z-w&t!qYmVjSgOdLaI6F3B0y)v%s#pXlUiQfNcK8H5E2A+x*u$rZnWkThPeMaJ7G4C1z-w&t!pYmVjRgOdL9I6E3B0y(v%s#oX

23、lTiQfNbK8H5D2A-x*u$qZnWkShPeMaJ7F4C1z)w&t!pYmUjRgOcL9I6E3B+y(v%r#oXlTiQeNbK8G5D2A-x*t$qZnVkShPdMaI7F4C0z)w&s!pXmUjRfOcL9H6E3B+y(u%r#oWlTiQeNbJ8G5D1A-x*t$qYnVkScL9I6E3B+y(v%r#oXlTiQeNbK8G5D2A-x*t$qZnVkShPdMaI7F4C0z)w&s!pYmUjRfOcL9H6E3B+y(u%r#oWlTiQeNbJ8G5D1A-x*t$qYnVkSgPdMaI7F3C0z)v&s!pXmUiRfOcK9H6E2B+x(u%rZoWlThQeNbJ8G4D1A-w*t$qYnVjSgPdLaI7F3C0y)v&s#pXmUiRfNcK9H5E2B+x(u$rZoWkThQeMbJ7G4D1z-w*t!qYmVjSgOdLaI6F3C0y)v%s#pXlUiRfNcK8H5E2A+x(u$rZnWkThPeMbJ7G4C1z-w&t!qYmVjRgOdL9