集合—变分混合同化中集合更新方法及循环同化预报研究
发布时间:2020-08-04 05:45
【摘要】:集合-变分混合同化方法将集合卡尔曼滤波法和变分法结合,具有计算成本低、便于同化多种观测资料及背景误差协方差具有流依赖性的优点。混合同化中,能描述天气形势演变特征的流依赖的背景误差协方差隐含于集合扰动中,因此集合预报的效果很大程度上影响着集合-变分混合同化的效果。为建立简便可行的、稳定高效的集合-变分循环同化预报系统,本文对集合-变分混合同化中的集合更新方法进行了一系列研究。首先,为讨论集合变换卡尔曼滤波中不同协方差膨胀方案对混合同化及预报效果的影响,开展了中国中东部区域一次连续大范围降水过程的连续10天的循环同化和预报试验。结果表明:几种协方差膨胀方案相对于无协方差膨胀的方案,均有效地提高了混合同化和预报的效果:对于风场,将预报误差协方差投影到集合子空间的方案和采用平均新息协方差信息的方案表现较好;对于温度场、湿度场和降水预报,采用平均新息协方差信息的方案和采用了同化时次前两次集合预报比率的方案最好。其次,针对集合变换卡尔曼滤波在同化过程中易受不同时次观测资料差异的影响而产生较大程度的协方差震荡,从而可能导致系统不稳定的问题,构建了一种基于物理控制变量扰动的混合同化预报系统。该系统通过对分析场的随机扰动并经过模式变换得到扰动量,从而实现集合更新,而集合扰动的误差统计由静态背景误差协方差提供。由于静态背景误差协方差的平衡属性,因此产生的扰动受平衡约束。连续10天的循环同化及预报试验表明,本同化方案效果明显优于三维变分方案,动力场的整体同化和预报效果与集合转换卡尔曼滤波方案基本相当。本方案不受观测资料差异的影响,在没有经任何参数调试情况下,取得了良好同化和预报效果,为混合同化的便捷运行提供了一种稳定可靠的手段。第三,为了在不增加或仅增加很小计算代价和存储成本的前提下,提高目前基于变分同化的快速更新同化系统的性能,本文初步研究了基于时间滞后集合的快速更新混合同化方案。试验结果表明:时间滞后集合充分利用了历史观测,使各集合成员包含了其对应时次的同化及预报信息,与动力模式协调且具有动力增长结构,因而时间滞后集合体现了天气流型即预报误差协方差的流依赖性,对于动力场和热力场的改进作用明显。
【学位授予单位】:南京信息工程大学
【学位级别】:硕士
【学位授予年份】:2015
【分类号】:P456
【图文】:
2.2邋WRFDA同化系统及其混合同化方法逡逑2.2.1邋WRFDA同化系统逡逑WRFDA系统在WRF数值模式中的作用如图2.1,其中蓝色部分为WRFDA同化系统功逡逑能。图2.1中,3表示由WRF或者由WPS/REAL输出结果得到的初猜场(背景场);x?表逡逑示由WPS/REAL得到的侧边界条件;是由WRFDA同化系统同化观测资料后生成的分析逡逑场;V是WRF模式进行预报得到的预报场;扩为经过OBSPROC邋(观测前处理)得到的观逡逑测场。(需要注意的是:PREPBUFR,雷达,辐射和降水资料不需要经过OBSPROC预处理);逡逑7逡逑
合成员的循环同化及预报试验,每12小时同化…次。在前述四个协方差膨胀方案的基础上,逡逑增加了一组无协方差膨胀方案NOIF邋(No邋Inflation)的试验,共设计了五组协方差膨胀方案逡逑对比试验。试验区域如图2.2,水平方向格点数为160x110,垂直方向38层,水平分辨率为逡逑27km,模式层顶气压为50hPa。采用的模式为WRFV3.5.1,模式微物理过程采用WSM6方逡逑案,积云对流参数化采用Kain-Fritsch对流方案,长波辐射方案为RRTM方案,短波辐射方逡逑案为Dudhia方案,行星边界层方案为Mellor-Yamada-Janjic方案,陆面过程米用Noah-MP逡逑方案》侧边界条件由每6小时一次的GFS再分析资料提供。逡逑在试验中,静态背景误差协方差矩阵由NMC方法计算得到(Parrish邋and邋Derber邋1992),逡逑样本为2011年7月1日00时起,连续一个月每隔12小时一次的WRF模拟的12小时预报逡逑场和24小时预报场。Hybrid集合扰动中,初始集合成员是在GFS再分析场的基础上通过逡逑RandomCV方法扰动得到,在之后的同化时次里采用ETKF方法更新,静态背景误差协方差逡逑和集合背景误差协方差的权重各占一半
协方差膨胀系数数值的变化反映了不同协方差膨胀方案对误差协方差估计的趋势。逡逑为直观比较不同ETKF协方差膨胀方法,本章对协方差膨胀系数随各同化时次的变化进行了逡逑分析。从各协方差膨胀系数随各同化时次变化曲线(图2.5)可以看出,WG07方案协方差逡逑膨胀系数的数值最小,TRNK方案协方差膨胀系数要大于WG07方案,WG03和BOWL方逡逑案协方差膨胀系数则明显较大,其屮WG03协方差膨胀系数在前8个同化时次是最大的。逡逑从图中也可以看出,WG07方案和TRNK方案协方差膨胀系数比较稳定,WG07方案数值基逡逑本维持在8上下,TRNK方案协方差膨胀系数经过3次循环后趋于稳定,数值维持在20-30逡逑之间;BOWL方案经过大概5次循环趋于稳定,数值在40-50附近,但是存在波动现象;逡逑WG03协方差膨胀系数经过3次循环后数值在50上下波动。几种方案的协方差膨胀系数随逡逑同化时次变化不同表现说明:(1)由于WG07方案采取了将预报误差投影到集合子空间的逡逑方法以改善采样误差
【学位授予单位】:南京信息工程大学
【学位级别】:硕士
【学位授予年份】:2015
【分类号】:P456
【图文】:
2.2邋WRFDA同化系统及其混合同化方法逡逑2.2.1邋WRFDA同化系统逡逑WRFDA系统在WRF数值模式中的作用如图2.1,其中蓝色部分为WRFDA同化系统功逡逑能。图2.1中,3表示由WRF或者由WPS/REAL输出结果得到的初猜场(背景场);x?表逡逑示由WPS/REAL得到的侧边界条件;是由WRFDA同化系统同化观测资料后生成的分析逡逑场;V是WRF模式进行预报得到的预报场;扩为经过OBSPROC邋(观测前处理)得到的观逡逑测场。(需要注意的是:PREPBUFR,雷达,辐射和降水资料不需要经过OBSPROC预处理);逡逑7逡逑
合成员的循环同化及预报试验,每12小时同化…次。在前述四个协方差膨胀方案的基础上,逡逑增加了一组无协方差膨胀方案NOIF邋(No邋Inflation)的试验,共设计了五组协方差膨胀方案逡逑对比试验。试验区域如图2.2,水平方向格点数为160x110,垂直方向38层,水平分辨率为逡逑27km,模式层顶气压为50hPa。采用的模式为WRFV3.5.1,模式微物理过程采用WSM6方逡逑案,积云对流参数化采用Kain-Fritsch对流方案,长波辐射方案为RRTM方案,短波辐射方逡逑案为Dudhia方案,行星边界层方案为Mellor-Yamada-Janjic方案,陆面过程米用Noah-MP逡逑方案》侧边界条件由每6小时一次的GFS再分析资料提供。逡逑在试验中,静态背景误差协方差矩阵由NMC方法计算得到(Parrish邋and邋Derber邋1992),逡逑样本为2011年7月1日00时起,连续一个月每隔12小时一次的WRF模拟的12小时预报逡逑场和24小时预报场。Hybrid集合扰动中,初始集合成员是在GFS再分析场的基础上通过逡逑RandomCV方法扰动得到,在之后的同化时次里采用ETKF方法更新,静态背景误差协方差逡逑和集合背景误差协方差的权重各占一半
协方差膨胀系数数值的变化反映了不同协方差膨胀方案对误差协方差估计的趋势。逡逑为直观比较不同ETKF协方差膨胀方法,本章对协方差膨胀系数随各同化时次的变化进行了逡逑分析。从各协方差膨胀系数随各同化时次变化曲线(图2.5)可以看出,WG07方案协方差逡逑膨胀系数的数值最小,TRNK方案协方差膨胀系数要大于WG07方案,WG03和BOWL方逡逑案协方差膨胀系数则明显较大,其屮WG03协方差膨胀系数在前8个同化时次是最大的。逡逑从图中也可以看出,WG07方案和TRNK方案协方差膨胀系数比较稳定,WG07方案数值基逡逑本维持在8上下,TRNK方案协方差膨胀系数经过3次循环后趋于稳定,数值维持在20-30逡逑之间;BOWL方案经过大概5次循环趋于稳定,数值在40-50附近,但是存在波动现象;逡逑WG03协方差膨胀系数经过3次循环后数值在50上下波动。几种方案的协方差膨胀系数随逡逑同化时次变化不同表现说明:(1)由于WG07方案采取了将预报误差投影到集合子空间的逡逑方法以改善采样误差
【参考文献】
相关期刊论文 前10条
1 刘梅;俞剑蔚;韩桂荣;张鹏勃;;2011年7月中旬江淮流域持续性强降水的环流特征分析[J];热带气象学报;2014年01期
2 龚建东;;同化技术:数值天气预报突破的关键——以欧洲中期天气预报中心同化技术演进为例[J];气象科技进展;2013年03期
3 熊春晖;张立凤;关吉平;陶恒锐;苏佳佳;;集合—变分数据同化方法的发展与应用[J];地球科学进展;2013年06期
4 徐枝芳;郝民;朱立娟;龚建东;万齐林;陶士伟;薛纪善;陈德辉;;GRAPES_RAFS系统研发[J];气象;2013年04期
5 陈葆德;王晓峰;李泓;张蕾;;快速更新同化预报的关键技术综述[J];气象科技进展;2013年02期
6 沈艳;潘e
本文编号:2780126
本文链接:https://www.wllwen.com/projectlw/qxxlw/2780126.html