基于对偶数理论的资料同化新方法

针对变分资料同化中目标泛函梯度计算精度不高且复杂等问题, 提出了一种基于对偶数理论的资料同化新方法, 主要优点是: 能避免复杂的伴随模式开发及其逆向积分, 只需在对偶数空间通过正向积分就能同时计算出目标泛函和梯度向量的值. 首先利用对偶数理论把梯度分析过程转换为对偶数空间中目标泛函计算过程, 简单、高效和高精度地获得梯度向量值; 其次结合典型的最优化方法, 给出了非线性物理系统资料同化问题的新求解算法; 最后对Lorenz 63混沌系统、包含开关的不可微物理模型和抛物型偏微分方程分别进行了资料同化数值实验, 结果表明: 新方法能有效和准确地估计出预报模式的初始条件或物理参数值.     

New data assimilation system DNDAS for high-dimensional models

The tangent linear(TL) models and adjoint(AD) models have brought great difficulties for the development of variational data assimilation system. It might be impossible to develop them perfectly without great efforts, either by hand, or by automatic differentiation tools. In order to break these limitations, a new data assimilation system, dual-number data assimilation system(DNDAS), is designed based on the dual-number automatic differentiation principles. We investigate the performance of DNDAS with two different optimization schemes and subsequently give a discussion on whether DNDAS is appropriate for high-dimensional forecast models. The new data assimilation system can avoid the complicated reverse integration of the adjoint model, and it only needs the forward integration in the dual-number space to obtain the cost function and its gradient vector concurrently. To verify the correctness and effectiveness of DNDAS, we implemented DNDAS on a simple ordinary differential model and the Lorenz-63 model with different optimization methods. We then concentrate on the adaptability of DNDAS to the Lorenz-96 model with high-dimensional state variables. The results indicate that whether the system is simple or nonlinear, DNDAS can accurately reconstruct the initial condition for the forecast model and has a strong anti-noise characteristic. Given adequate computing resource, the quasi-Newton optimization method performs better than the conjugate gradient method in DNDAS.     

Impact of GNSS radio occultation bending angle data assimilation in YH4DVAR system

Global Navigation Satellite System （GNSS） radio occultation measurements have been assimilated into the four- dimensional variational data assimilation system （YH4DVAR） using a one-dimensional bending angle operator （GBAO） as a new type of observation. For the sake of verifying the impact of GNSS radio occultation （RO） measurements to the data assimilation system, three experiments have been conducted. The statistical results of the analysis error experiment and forecast skill experiment show that the GNSS RO measurements have an impact on the analysis system. The typhoon forecast experiment shows the impact on the important weather process. They all have a positive impact on the weather forecast. Lastly, we look forward to future work on the observation system simulation experiment （OSSE） to investigate the impact of GNSS RO measurements as a function of observation number, which is an effective method to estimate the saturation of the observation number.     

GPS掩星一维弯曲角算子在四维变分资料同化系统中的实现方法研究

GPS无线电掩星资料的良好特性使其作为十分重要的观测资料在欧洲 及美国数值天气预报中效果明显. 本文设计了一维弯曲角算子在四维变分资料同化系统中的实现方法和质量控制策略. 详细阐述了一维弯曲角正演算子及其切线性/伴随算子在全球四维变分资料同化系统中的实现步骤和过程, 其中折射率计算中采用Rüeger的折射率最佳平均因子的三项表达式, 计算出更加精确的折射率值. 对2012年8月份的GPS无线电掩星资料进行了观测误差统计, 引入了更加符合四维变分资料同化系统实际情况的一维弯曲角观测误差模型和质量控制策略. 同时, 对比了观测误差统计试验中同化GPS无线电掩星资料的模式背景场和未同化GPS无线电掩星资料的背景场的观测误差统计结果, 发现同化GPS无线电掩星资料之后的观测误差统计结果在全球区域正作用十分明显, 尤其是在中低纬度地区, 高度在10 km 以上的区域观测误差减小十分明显. 关键词： GPS无线电掩星 一维弯曲角算子 四维变分同化系统 质量控制