两类滤波方法在估算土壤湿度方面的对比研究-气象学专业毕业论文.docxVIP

摘 摘 要 集合卡曼滤波(EnlF)利用蒙特卡罗方法统计背景场的误差协方差矩阵， 解决了卡曼滤波在数据同化时需要对模型算子和观测算子线性化的缺点，易于使 用，被广泛地应用到陆面数据同化中来估算土壤湿度廓线。 EnKF虽然可以用非线性模式产生集合背景预报场，但是最优分析场是建立 在模式误差间关联为线性且误差为高斯分布的假设上，而实际土壤湿度方程为高 度非线性。同时，土壤湿度又存在固定取值范围；过干和过湿都会造成集合样本 偏斜。为了考查这些不利因素可能的影响，特引入非线性粒子滤波，来定量比较 它们在估算土壤湿度廓线的速度和精度等方面异同。 首先，本论文引入了基于贝叶斯原理计算粒子的后验概率分布的粒子滤波 中采样重要性重采样算法(Sampling importance Resampling，SIR)，它不需模式误 差为正态以及误差为线性关联的假设，并且可以提供误差分布的高阶矩信息；它 在更新背景场样本时每个粒子的权重不同，而集合卡曼滤波粒子的权重都是相同 的。其次，为了评估土壤湿度样本在土壤过干或过湿时出现的偏斜性甚至双模态， 考虑到西北黄土高原陆面过程的重要性，本论文选取黄土土壤的逐渐变干过程束 进行观测模拟同化试验。最后，选用2003年土壤湿度试验(SMEX03)的站点 观测资料进行实际大气观测驱动下同化试验。通过上述实验可以得出： l、估算土壤湿度的EnKF同化方案虽然不满足集合卡曼滤波最优解需要的 条件，但它却能快速地估算出土壤湿度廓线，加之所需粒子样本较少，计算量也 较小：而SIR虽然无此限制，但只有在大样本的条件下才能达到相同的土壤湿度 估算精度，因此无论从速度和计算量来看都比EnKF逊色。 2、两者样本统计分布随同化过程也完全不同，EnKF的样本集合的边缘概 率密度分布始终保持为一个单峰分布，而SIR方法的样本粒子却经历了从一个峰 值到两个再到一个分布的变化过程。 3、EnKF同化方案在改进土壤湿度廓线估计的同时也使计算的地表感热和 潜热通量更接近真值。 基于上述结论，EnKF基本满足陆面数据同化系统对同化算法的要求，可用 于同化表层土壤湿度观测(或微波量温)束估算土壤湿度廓线。 关键词：集合卡曼滤波(EnKF)粒子滤波 关键词：集合卡曼滤波(EnKF)粒子滤波 公用陆面模式(CoLM)土壤湿度 资料同化 Richards方程 AbstractThe Abstract The ensemble Kalman filter(EnKF)calculates background error covariance matrix using Monte-Carlo method and is able to resolve the nonlinearity and discontinuity exist within model operator and observation operator．At the same，it is an easy to use，flexible，and efficient data assimilation algorithm widely used in Soil moisture Assimilation System． Although EnKF can use non-linear models to gain the background forecast field， but it based on the normality approximation of model error and observational error which is the assumption of Gaussian distribution．However，the soil moisture equation is highly nonlinear,and soil moisture is fixed range，The sample be hiighly skewed toward the wet dry ends．In order to examine the possible impact of these negative factors，This paper gives the results from the EnKF compared with those obtained from sequential importance resampling(SIR)particle filter that is of nonlinear filters． First，this thesis introduces the Sampling importance Resampling，It