天气学分析诊断 2.2.2 平滑滤波与尺度分离 Barnes滤波.pptVIP

南京信息工程大学大气科学系 王咏青 7.具体计算步骤： （1）输入测站（或格点）经纬度、要素值等初始资料； （2）选择两组合适的滤波常数C1、G1和C2、G2； （3）分别取两组参数用公式（3）、（5）或（9）或计算所有点的低通滤波值； （4）用公式（13）求出每个点的带通滤波值； （5）输出结果。 南京信息工程大学大气科学系 王咏青 参考文献 Maddox R. A. An objective technique for separating macroscale and mesoscale feature in meteorological data. Mon. Wea. Rev. , 1980:1108-1121 陈忠明.常规资料获取中尺度信息的一种简便方法. 气象 陈忠明.气象场中尺度带通滤波研究. 气象学报.1992 段旭等.一次滇西南秋季暴雨的中尺度分析与诊断.高原气象.2003 南京信息工程大学大气科学系 王咏青 南京信息工程大学大气科学系 王咏青 Barnes客观分析及带通滤波法分离气象场中的中尺度信息 大尺度系统可以提供中尺度系统形成的条件和环境场，而中尺度系统才是直接造成暴雨等灾害的天气系统，通过中小尺度系统内的对流活动，暴雨过程又可对大尺度背景场起反馈作用，改变原有的风场、湿度场和大气层结等。 1. 引言 南京信息工程大学大气科学系 王咏青 为了研究不同尺度系统之间的相互作用，可以对空间场进行尺度分离。 尺度分离的方法有多种，例如：谐波分析，即对半球或某一纬度带地区，把一种气象场分解成许多谐波，然后讨论各谐波的贡献和相互之间的作用。适用于半球或整个纬圈。 南京信息工程大学大气科学系 王咏青 对于区域的尺度分离，Barnes S. L. 于1964年根据高斯(Gaussian)权重函数发展了一个中尺度滤波法，后又于 1973、1975、1980、1996年发表论文对此方法进行修正提高。 人们称其为Barnes 中尺度空间滤波法，目前，在MM5、ARPS、RAMS等世界著名的中尺度模式中均使用这种方法。 南京信息工程大学大气科学系 王咏青 2. 原理 设P等压面上某一气象要素的观测值为 ，k=1,2,3,……，k为测站序号， 分别为测站的坐标。根据Barnes滤波原理,由观测值 确定的低通滤波初值场为 (1) 南京信息工程大学大气科学系 王咏青 其中，Wk为权重函数 （2） 连续函数 C 为滤波常数 为测站 到 的距离 M为参加点 处滤波的资料样本数 南京信息工程大学大气科学系 王咏青 南京信息工程大学大气科学系 王咏青 如果 (x,y)点的气象观测值为F(x,y)=Asinkx（假设沿y方向是均匀分布的），其中A是振幅，k是波数，k=2π/λ，λ为波长。经过一系列数学变换，可以得到： =exp(-4π2c/ λ2)Asinkx =R0F(x,y) 南京信息工程大学大气科学系 王咏青 R0为响应函数： 其中：k 为波数； 为波长； C 为滤波常数 （3） 南京信息工程大学大气科学系 王咏青 响应函数R0表示，经过加权平均分析后所得到的值与观测值之间的响应程度，或者滤波后的波振幅与原来波振幅之比。 如果R0 =1，表示此波在滤波过程中完全保留下来，不受削弱；如果R0 =0，则完全滤去。一般R0 在0-1之间。 对于不同的c值，R0 和λ之间的关系见下图： 南京信息工程大学大气科学系 王咏青 其关系示意图如下： 南京信息工程大学大气科学系 王咏青 （1）对波长很长的波，响应函数趋于1，即： （2）C 增大时可以相当好地滤去极小尺度的波，即某些噪音。 从图中可见，响应函数有如下的性质： 南京信息工程大学大气科学系 王咏青 对获取的初值场 做进一步订正： （5） 其中： （6） G 为另一个滤波常数，取值通常为0.2～0.4 3. 改进方案之一 修正后的权重函数 南京信息工程大学大气科学系 王咏青 测站处要素值与初值的差值可直接得到： （4） 南京信息工程大学大气科学系 王