基于雷达反射率反演降水率的一维%2b三维变分同化研究及试验.pdfVIP

基于雷达反射率反演降水率的一维+三维变分同化研究及试验 1 2 2 1 张诚忠 薛纪善 张林 万齐林 1、 中国气象局广州热带海洋气象研究所 2 、中国气象科学研究院 1、引言 根据雷达的探测特性，多普勒雷达在中尺度观测上具有其独特的优势：时间与空间分辨 率高，三维的云与降水、风场等信息，因此许多国内外研究者通过同化雷达资料改善短时临 近降水预报，取得良好效果。但也存在一些问题，如同化雷达资料后模式预报对降水的落区、 量级的改进主要体现在前 3 小时，对随后的降水预报改善不大。为了解决这问题，国内外的 学者对同化天气雷达资料做了大量研究，提出了一些方案，使同化后的温度湿度场保持协调， 因而模式对这些信息的记忆更长久，改善短时降水预报。 天气雷达观测可以提供有关云与降水以及风场的空间分布信息，同时还可以估计凝结加热 的空间分布。但多数情况下，这些信息并不是模式的直接基本预报量，在同化前常需要利用统 计或物理模型做一定的反演。它们被同化进模式的方式，有 3 种不同的方式：第一种方式先将 雷达资料反演为模式能接受的参数，再使模式接受这写参数。例如通过雷达反射率估计出降水 信息并进一步推算降水引起的凝结加热，将降水潜热信息采用牛顿松弛逼近（nugding）强迫模 式接受，即所谓的 latent Heat Nudging，（Jones and Macpherson, 1997）。这类方案根据模 式的湿物理过程，即通过假设模式预报降水与观测量相当情况下，将过程产生的潜热量倾向增 加到热力倾向预报方程中，以订正修正模式的温度倾向预报。英国气象局（UKMO）开展基于雷 达的云分析、与降水分析，将这类方案应用于业务，并取得一定的效果。在国内，刘红亚等（2009） 通过 Z - qr 关系和定常暖云方案,反演了雨水、云水 、水汽和垂直速度，再利用牛顿松弛方 法将云雨信息引入数值模式，同化后降水预报有一定的改进。第二种方式是雷达反射率的直接 同化。该类方案将反射率作为观测量直接进入应用于三维（3Dvar）或四维（4Dvar）变分同化。 在MM5 和WRF中的雷达同化模块（xiao 2006.）中，控制变量为空气中的总含水量（含汽态与 水凝物）qt,观测量为雷达反射率。通过暖云过程导出雨水、云水和水汽。模式的雨水可模拟出 模式的反射率，再与观测的反射率比较。通过极小化计算（目标函数及其梯度的计算）得到模 式初始时刻的云雨信息。第三种方式为通过一维+三维或四维等（1Dvar + 3Dvar（4Dvar））两 个步骤完成雷达反射率的同化。这是由反演到直接同化的过渡方案，即采用简单的一维变分方 案将雷达观测资料反演为模式的预报变量，再用三维或四维变分同化方案将反演出来的值同化 进模式。该方案目前已在多个数值预报中心实现业务运行。例如，法国气象中心的 AROME系统 通过Bayesian方法从反射率廓线中估计出大气的一维比湿廓线，进而再经三维变分方法同化湿 度廓线。再如欧洲中期预报中心 （ECMWF） 数值预报中心以雷达（加上其他资料）导出的地面 瞬时降水率，通过1Dvar得到大气可降水量，通过4DVar同化到模式中。 尽管上述的雷达反射率同化方案不尽相同，但这些方案的基本科学问题仍归结为如何建立 正确反映雷达观测与分析变量联系的观测算子。通常大气中雷达回波往往伴随着云、对流和降 水，因此长期以来许多研究者在观测算子的建立中除了考虑动力过程之外，也越来越关注湿物 理过程的处理及其对同化分析结果的影响作用（Fillion et al 1997;Bauer et al 2006;Lopez.et al 2006）。Fillion 通过分别利用 Kuo-Anthes 和 Arakawa 对流方案为观测算子对降水率进行一 维变分同化，结果认为有助于提高湿度场的分析质量和资料稀缺地区的预报准确率。近年来， 复杂的观测算子包括线性的湿物理过程如大尺度凝结方案、对流参数化方案的观测算子已应 用到降水率、卫星资料的同化中，分析结果表明对分析场和预报场有着显著的改善 （Bauer et al 2006）。 目前国内关于将湿物理过程引入观测算子的同化研究已逐步开展。数值预报模式的湿物理 过程一般包含网格尺度的显式云及降水方案与对流参数化方案，但至今国内的研究均只考虑其 中的一项，限制了对各类天气形势的可用性。本研究拟通过以