An enhanced ensemble Kalman filter scheme incorporating model error in sequential coupling between flow and geomechanics

In this work, we construct a new methodology for enhancing the predictive accuracy of sequential methods for coupling flow and geomechanics while preserving low computational cost. The new computational approach is developed within the framework of the fixed-stress split algorithm procedure in conjunction with data assimilation based on the ensemble Kalman filter (EnKF). In this context, we identify the high-fidelity model with the two-way formulation where additional source term appears in the flow equation containing the time derivative of total mean stress. The iterative scheme is then interlaced with data assimilation steps, which also incorporate the modeling error inherent to the EnKF framework. Such a procedure gives rise to an “enhanced one-way formulation,” exhibiting substantial improvement in accuracy compared with the classical one-way method. The governing equations are discretized by mixed finite elements, and numerical simulation of a 2D slab problem between injection and production wells illustrate the tremendous achievement of the method proposed herein.     

基于雷达资料的上海地区暴雨面雨量计算及应用

根据上海地区小时雷达资料的定量降水估测(quantitative precipitation estimation,QPE)网格产品和经过质量控制的自动站雨量资料进行暴雨面雨量计算方法研究,计算方法为:将暴雨区内自动站周围9个网格点的QPE平均值作为自动站的雷达估测雨量,计算各测站雨量与雷达估测雨量的差值,用克里金插值方法将差值场插值到雷达降水估算产品相同的网格点上,再将网格点上残差插值数据加上网格点上雷达估测雨量,得到各网格点用自动站订正后的雷达估测小时雨量资料。文中24个典型暴雨各测站订正后QPE与实测雨量的平均绝对误差比订正前减小了27%;两个典型暴雨主要降水阶段经国家气象站资料订正后QPE与实测雨量的平均误差比订正前减小了33%～39%,暴雨过程总雨量经国家气象站资料订正后减小了34%～59%。根据以上方法得到的网格点雨量计算和绘制上海地区2007—2015年24个典型暴雨降水区域的小时面雨量值和图、过程面雨量值和图、行政区和水利片面雨量值和图。最后基于.NET Framework 4.0基础架构软件开发平台,使用Microsoft Visual Studio 2012和Microsoft Visual Studio软件开发工具制作"基于雷达资料暴雨面雨量自动化计算查询系统",供业务科研人员实时查询和计算。     

上海地区春季最高气温预报失败案例分析

本文分析了2015年3月17—18日上海地区连续两天发生最高气温预报失误的天气背景,并使用当日实况观测和业务预报使用的数值模式资料,剖析预报失败的原因,分析表明:对天空状况的误判是导致17日预报失败的主要原因,且东南风预报偏强更进一步增大了预报误差;冷空气影响时间的判断失误是导致18日预报失败的主要原因。从模式预报的实时检验、预报的跳跃性和不确定性角度分析了预报中存在的问题:预报员应重视本地和上游实况,从传统对单一确定性模式预报的依赖向多模式多起报时次及能提供概率预报和不确定性信息的业务集合预报的分析思路转型。此外,还需加强对集合预报的系统性检验、评估及数值预报释用产品的开发,增加包含不确定性信息的公众天气预报发布形式。     

GRAPES区域集合预报模式的初值扰动增长特征

基于GRAPES-REPS（Global and Regional Assimilation and Prediction Enhanced System-Regional Ensemble Prediction System）区域集合预报模式和集合变换卡尔曼滤波（Ensemble Transform Kalman Filter,ETKF）初值扰动方法,对2015年6月1~15日10 km与15 km水平分辨率分别进行集合预报试验,通过分析ETKF初值扰动分量、初值扰动方差准确率、动能谱、扰动能量演变、日变化及集合离散度、均方根误差等特征,揭示GRAPES-REPS区域集合预报ETKF初值扰动结构及增长特征。结果表明:（1）ETKF初值扰动方案产生的扰动能够保持所有正交、不相关方向的误差方差,且ETKF分量α参数值及放大因子具有较好的稳定性。（2）ETKF初值扰动方法生成的扰动场以大尺度扰动为主,扰动结构及能量具有随流型依赖特征,低层以内能扰动为主,高层以动能扰动为主,且集合扰动可以有效捕捉预报误差的结构。（3）GRAPES区域集合预报初值扰动总能量和集合离散度随预报时效的延长均呈发展趋势,但离散度增长率小于均方根误差增长率,即集合预报总体存在集合离散度不足的问题。（4）水平分辨率提高可以增加中高层大尺度扰动波谱能量,明显改进等压面及近地面风场及温度场的集合预报效果。值得指出的是,GRAPES-REPS区域集合预报低层内能扰动能量存在明显的日变化特征,特别是青藏高原地区更加显著,需要进一步研究青藏高原初值扰动结构的合理性。     

GRAPES区域集合预报尺度混合初始扰动构造的新方案

集合预报初始扰动能否准确反映预报误差的结构特征是决定区域集合预报质量的关键因素之一。本文针对GRAPES区域数值预报模式,发展设计了一种基于资料同化思想的混合尺度初始扰动构造新方案。该方案以全球大尺度信息为背景场,区域模式预报作为观测资料,借助GRAPES三维变分同化系统,将高质量的全球大尺度信息与区域模式预报中质量较高的中小尺度信息有效融合,构造混合尺度区域集合预报初始扰动,并通过个例试验和批量试验,比较分析了新方案和原区域集合预报的性能。试验结果表明,基于资料同化构造的初始扰动能够有效融合全球大尺度信息和中小尺度天气系统的信息,其降水概率预报更具参考价值。总体上看,区域集合预报混合初始扰动新方案能够较好地改进区域集合预报质量,尤其是对高度场和温度场效果更为显著,但对风场的集合预报性能影响略小。     

Impact of Perturbation Schemes on the Ensemble Prediction in a Coupled Lorenz Model

Based on a simple coupled Lorenz model, we investigate how to assess a suitable initial perturbation scheme for ensemble forecasting in a multiscale system involving slow dynamics and fast dynamics. Four initial perturbation approaches are used in the ensemble forecasting experiments: the random perturbation(RP), the bred vector(BV), the ensemble transform Kalman filter(ETKF), and the nonlinear local Lyapunov vector(NLLV) methods. Results show that,regardless of the method used, the ensemble ave...     

两种空气质量数值模式的应用评估与集合改进研究

利用2018—2019年国控站观测,评估CAMx和CMAQ模式对广东珠海主要污染物时空分布与演变特征的预报能力,并引入多元线性回归和随机森林方法对预报结果进行集成,探究不同集合方法的改进能力。结果表明：CMAQ在各污染物浓度季节-日变化方面明显优于CAMx,但两者存在明显系统偏差,并对多数污染物(除O₃之外)的昼夜和空间变化的模拟能力仍存在明显缺陷。例如,CMAQ合理地还原了CO、PM_2.5、PM₁₀、SO₂、O₃和NO₂的季节变化,相关系数介于0.72～0.84,但NMB分别达到-0.58、-0.18、-0.30、1.52,-0.16和-0.20,RMSE分别达到0.40 mg/m³、6.86、16.02、10.71、25.05和10.21μg/m³。同时,基于不同污染物构建的两种集合方法均有效移除了系统偏差,加强了CMAQ的模拟优势,并且随机森林方法明显优于多元线性回归,但两者均对模式缺陷无明显改进。进一步分...     

Impact of Coastal Radar Observability on the Forecast of the Track and Rainfall of Typhoon Morakot(2009)Using WRF-based Ensemble Kalman Filter Data Assimilation

This study explored the impact of coastal radar observability on the forecast of the track and rainfall of Typhoon Morakot(2009)using a WRF-based ensemble Kalman filter(EnKF)data assimilation(DA)system.The results showed that the performance of radar EnKF DA was quite sensitive to the number of radars being assimilated and the DA timing relative to the landfall of the tropical cyclone(TC).It was found that assimilating radial velocity(Vr)data from all the four operational radars during the 6 h immediately before TC landfall was quite important for the track and rainfall forecasts after the TC made landfall.The TC track forecast error could be decreased by about 43% and the 24-h rainfall forecast skill could be almost tripled.Assimilating Vr data from a single radar outperformed the experiment without DA, though with less improvement compared to the multiple-radar DA experiment.Different forecast performances were obtained by assimilating different radars, which was closely related to the first-time wind analysis increment, the location of moisture transport, the quasi-stationary rainband, and the local convergence line.However, only assimilating Vr data when the TC was farther away from making landfall might worsen TC track and rainfall forecasts.Besides, this work also demonstrated that Vr data from multiple radars, instead of a single radar, should be used for verification to obtain a more reliable assessment of the EnKF performance.     

Calibrating a Land Parcel Cellular Automaton (LP-CA) for urban growth simulation based on ensemble learning

The reliability of raster cellular automaton (CA) models for fine-scale land change simulations has been increasingly questioned, because regular pixels/grids cannot precisely represent irregular geographical entities and their interactions. Vector CA models can address these deficiencies due to the ability of the vector data structure to represent realistic urban entities. This study presents a new land parcel cellular automaton (LP-CA) model for simulating urban land changes. The innovation of this model is the use of ensemble learning method for automatic calibration. The proposed model is applied in Shenzhen, China. The experimental results indicate that bagging-Naïve Bayes yields the highest calibration accuracy among a set of selected classifiers. The assessment of neighborhood sensitivity suggests that the LP-CA model achieves the highest simulation accuracy with neighbor radius r = 2. The calibrated LP-CA is used to project future urban land use changes in Shenzhen, and the results are found to be consistent with those specified in the official city plan.     

集合方法在月动力预报信息提取中的应用

本工作将集合方法应用于提取月动力预报有用信息。利用中国气象局国家气候中心Ｔ６３Ｌ１６全球谱模式的５００百帕高度场月集合预报产品（集合成员数为８个，初始场的选取采用滞后方法（ＬＡＦ），即相邻两天的００００，０６００，１２００和１８００ＧＭＴ的初始化资料），就１９９７年１月至５月共１５次预报，分析了集合预报成员间的离散度与预报评分（距平相关系数和均方根误差）的关系，研究了用集合各成员预报离散度作为各个成员逐日预报的权重对月预报效果的影响。结果表明集合预报成员的离散度与预报评分有显著的相关，是有效预报长度Ｎ的一个很好估计；用离散度作为权重平均的月预报高度距平相关系数明显高于算术平均和线性权重，此外个例分析表明月平均环流及其异常的预报得到明显的提高。