用变分法修正黑河试验数据的初步研究＊.PDF

第 一 ， 卷 撇 期 大 气 科 学 功 n 年 3 月 ChnieseJournalofAtmosphereSciences VMoal,r2.6N20o0,22 用变分法修正黑河试验数据的初步研究* 马 刚 王 强 邱崇践 P+ A (兰州大学大气科学系. “州 刀的00 摘 要 为了对黑河试验观Nil资料进行订正.利用一个一维边界层垂直扩散模式和以拟牛 顿极小化算法兰基础的最优化方法建立了 个变分同化系统，其中、目标函数相刘十控制变 量的梯度用扰4.法求取 并利用模式资料和黑河试验若T时刻的观测资料迸行了 些数值试 验_试验证明，将初始场变量与模式误差系数同时作为控制变量加以调整 能够有效调整初 始场观测资料 并且 ‘用变分同化方法能够较线性插仇 (时间 ‘空间)方法更有效 合理的 修正黑河试验观侧资料场的种种缺陷 关健词:变分同化;边界层 ‘黑河试验，扰动法 1 引言 黑河试验从}990年6月至1992年10月 对欧亚大陆腹地典型的干早气候地区黑 河流域的沙漠、戈壁和绿洲不同下垫面的辐射过程，边界层过程、绿洲生物过程，沙尘 辐射特征和干HT地区降水特征进行了观测.获得1大量宝贵的资料 但由于各种原因. 试验所获得的观测资料在时间 卜有不少间断，同时，一些资料也存在较大误差 从而影 响了对这些资料的进 一步应用 在边界层内 各种要素场的时空变化都很剧烈，插ff分 析和平滑的方法会产生较大误差 另一方面，边界层大气也遵循一定的物理规律，即给 定初始场后 积分适当的大气模式就可以模拟出它们的演变过程 这就是说，有可能利 用数值模式提供的信息来补充和修正观测资料的缺陷。这就是当前数值天气预报中流行 的变分四维同化朗基本思想 它既可以充分利用已有的物理规律 (数f直模式) 又能合 理的使用时间演变的历史资料所包含的信息 此时目标函数被定义为模式输出与观测的 距离，因此变分四维同化过程也可以看作是将模式输出与实际观fl相拟合 显然 这种 方法也可以用 1-.解决黑河试验 (HEIFE)资料的订正和插补工作 本文将利用一个一 维边界层数值模1 用变分同化方法对黑河观测资料进行订正 Z 理论概述 2.1 变分四维同化方法 变分四维同化方法的基本原理是通过修正模式输人量 (控制变量〕，使在同化时段 2000-01-14收到.2(g)]-05-t4收到修改稿 国家自然科学纂金资助项 日4955249 大 气 科 学 26卷 内的模式输吕量 (或是它们的某些函数)与相应观测之间距离最小 因此，我们将日标 函数写成最价二乘形式，将其定义为模式预报场与相应观测场之间的距离 于是，变分 同化方法便构成了一个约束条件 (数值模式)下的泛函极小化问题、通过计算出目标函 数相对于控制变量的梯度，然后运用某种极小化算法 (如共扼梯度法、拟牛顿方法〕， 使其达到极小点 一个经空间离散化后的大气边界层数值模式方程可以写为 卜面的常微分方程组 dX F(X,p (1) dt 其中.X是,,;维向量，表示模式输出，包含所有模式变量在所有格点上的值;p是模式 误差参数 同时，定义表征模式输出与相应观测之间距离的目标函数 J=y{C[X(t)]一。(t,.)obiWI(，){Cx[(t)]一，(，)6s}, (2) 这里，。是同化时段内有观测场的时间次数，川t,).b、是观测场、并满足以t,,)= C[X(t,)J,C表示模式输出场X(t,)到观测场v(t)的映射，这是因为观测场并不一定都 是模式变量 W,(t)是一个NxN维的权重系数矩阵 表明目标函数中