共查询到16条相似文献，搜索用时 109 毫秒
1.
线性化物理过程能够改善四维变分同化中极小化收敛的稳定性和增加极小化过程中对大气物理过程和动力更加精确的描述,它是四维变分同化中非常重要的一部分。通过在GRAPES全球模式中研究线性化物理过程,尤其是两个湿线性化物理过程,改善切线性模式预报精度,来提高GRAPES全球四维变分同化的分析和预报效果。线性化物理过程的开发首先需要简化原非线性化物理过程中的强非线性项,然后对线性化物理过程进行规约化,以抑制切线性扰动的异常增长。目前GRAEPS全球模式中的线性化物理过程主要包括次网格尺度地形参数化、垂直扩散、积云深对流和大尺度凝结。线性化物理过程预报精度的检验方法是通过选择合适大小的初始扰动（同化分析增量）,来比较非线性模式和切线性模式中的扰动演化的纬向平均误差。然后以绝热版本的切线性模式为基础,通过冬、夏两个个例试验来分别检验4个线性化物理过程的12 h预报效果。试验结果表明,通过添加次网格地形参数化和垂直扩散两个干线性化物理过程方案,可以有效抑制住绝热版本切线性模式低层扰动的异常增长,大幅度改善切线性模式预报效果。通过添加积云深对流和大尺度凝结两个湿线性化物理过程,可以在热带区域和中、高纬度地区提高切线性模式中湿变量和温度变量的近似精度,提高GRAPES全球四维变分同化的分析和预报效果。 相似文献
2.
持续发展和优化切线性模式的线性化物理过程,保持与非线性模式一致是改善四维变分同化(4DVar)分析和预报效果的有效方法之一。目前业务系统的CMA-GFS模式采用基于Charney-Phillips(C-P)跳点的边界层参数化方案,而CMA-GFS 4DVar系统中采用基于Lorenz跳点的边界层线性化方案。为改善CMA-GFS 4DVar系统的边界层分析和预报效果,基于C-P跳点的边界层参数化方案研发了新边界层线性化方案,并通过对方案中地表热量通量和水汽通量扰动、自由大气的理查逊系数扰动、边界层的热量和动量交换系数扰动等进行更加精细地规约化约束,在确保CMA-GFS切线性和伴随模式稳定运行的情况下,减少线性化过程对切线性模式预报精度的影响。切线性近似试验检验表明:相较于原方案,新边界层线性化方案可以减少边界层位温和比湿的相对误差,最大可减少10%。批量4DVar循环同化试验表明:新边界层线性化方案可以有效改善切线性模式对低层位温、风场和比湿扰动的预报精度,减少4DVar内外循环目标泛函的相对差异,并提高700 hPa位势高度的可预报时效。 相似文献
3.
基于总能量模的奇异向量扰动常用于构造集合预报的初始条件。以建立GRAPES（Global and Regional Assimilation PrEdiction System）全球集合预报系统为目的,基于前期研发的GRAPES全球模式奇异向量方法,在GRAPES全球切线性模式和伴随模式2.0版的框架下,开展了引入线性化边界层方案来改善奇异向量结构,并提高奇异向量计算效率的研究。通过连续试验,从奇异向量的扰动能量结构、扰动能量谱及扰动空间分布等方面,综合分析改进GRAPES全球奇异向量的结构及演变特征。试验结果表明,改进后的GRAPES奇异向量方法有效抑制了之前扰动能量在近地面层不合理的快速增长,同时,奇异向量最优扰动的结构更客观地体现了中高纬度区域大气初始条件中的斜压不稳定扰动及其演变,如在初始时刻奇异向量扰动能量主要位于对流层中层,并呈现出随高度向西倾斜的大气斜压特征;经过线性化演变,扰动能量向较大水平尺度转移,并在垂直结构上表现出向对流层高层上传及向对流层低层下传的特征等。针对GRAPES奇异向量迭代求解中伴随模式计算耗时为主的情况,改进伴随模式中广义共轭余差方案的调用方式,并采用大内存存储法来提高其计算效率,进而将奇异向量总计算时间缩短了25%。总之,改进后的GRAPES奇异向量方法,可应用于构建面向业务应用的GRAPES全球集合预报系统。 相似文献
4.
伴随技术是四维变分同化（4DVar）系统中计算代价函数梯度的最佳办法,切线性和伴随模式的效果和效率直接影响着4DVar系统的发展。基于GRAPES（Global and Regional Assimilation PrEdiction System）全球切线性和伴随模式1.0版本,利用GRAPES全球模式2.0版本在并行框架和性能等方面的改善,重新优化和设计了GRAPES全球切线性伴随模式2.0版本,提高了GRAPES全球切线性和伴随模式的效果和效率,优化了切线性模式程序结构,使其计算时间最优可控制在非线性模式的1.2倍以内;采用在切线性模式中保存基态的方法,重构了伴随模式的程序结构,使其计算时间最优控制在非线性模式的1.5倍以内;在GRAPES全球切线性物理过程的设计中,将线性物理过程的轨迹基态计算和切线性扰动计算解耦,提高了GRAPES全球切线性和伴随模式的计算效果和效率。 相似文献
5.
GRAPES模式中地面通量在一次台风数值模拟中的敏感性试验研究 总被引：1，自引：0，他引：1
在中尺度天气数值预报模式（GRAPES）的边界层方案（MRF）中,地面通量的计算与模式分层有直接关系,可能引起计算结果的不准确。在模式中引入Beljaars方案,通过比较分析模拟的台风地表特征量及形势场和风场表明,引入的方案对台风强度和路径的预报效果都起到改善作用。对台风“巨爵”个例的模拟表明,引入Beljaars方案之后,其地表的潜热、感热通量均明显增强,摩擦速度最大值增大,说明地表向上传输的热通量、水汽通量、动量通量均增强,这都有利于台风的发展和加强,因此模拟的台风路径和强度也更接近实况;从台风结构来看,方案二模拟的台风倾斜垂直结构没有方案一的明显,暖心结构有所加强;从风场分析来看,台风眼区的风速低值中心可延伸到500 hPa的高度上,明显高于方案一。 相似文献
6.
以发展基于奇异向量技术为初值扰动的GRAPES全球集合预报系统为目的,在GRAPES模式及其干动力框架下的切线性、伴随模式基础上开展了以总能量模为权重算子的奇异向量计算技术研究,建立奇异向量的计算求解模块,并通过奇异向量检验方法和切线性近似方法验证了奇异向量求解的正确性.通过对中高纬度的GRAPES奇异向量水平结构的线性演变分析,证实了在最优时间间隔内GRAPES奇异向量能够快速增长,并能描述中高纬度大气的斜压不稳定特征.分析在初始时刻和最优化时间间隔时刻的GRAPES奇异向量总能量及其分量(动能和势能)的垂直分布特征,发现在中高纬度区域,GRAPES奇异向量能够描述对流层不同层次的斜压不稳定增长特征. 相似文献
7.
GRAPES全球四维变分资料同化系统需要积云深对流参数化方案的线性化与伴随方案,直接采用原始复杂参数化方案进行线性化并不可行,需要发展简化光滑方案来减缓非线性与非连续性特征。GRAPES全球模式采用NSAS积云对流参数化方案,积云深对流对环境的反馈主要通过补偿下沉来实现,研究突出补偿下沉作用,忽略降水蒸发、动量反馈等贡献,形成简化方案。采用输入温、湿度廓线加入不同幅度小扰动方法,评估参数化方案计算的温度、比湿时间倾向对输入扰动的敏感性,检验非线性与非连续特征。提出避免或减缓非连续“开关”的方法,在简化方案的基础上发展了简化光滑方案。简化光滑方案与原始积云深对流方案相比,在对流触发上一致,在对流的位温与比湿倾向、降水的时序模拟等方面相似,而在减缓非线性、避免非连续性方面显著优于原始方案。基于简化光滑方案发展的线性化方案表明,对小于2倍分析增量幅度的扰动,线性化方案可以较好地模拟非线性方案的扰动发展。发展的简化光滑方案具有合理性和实用性。 相似文献
8.
设计了适用于四维变分同化系统的扰动预报模式GRAPES_PF。根据GRAPES的地形追随坐标非静力原始方程组,采用小扰动分离方法推导微分形式的线性扰动预报方程组,并利用与GRAPES非线性模式相似的数值求解方案求解线性扰动微分方程组。在设计扰动预报模式时采用了两个时间层半隐式半拉格朗日方案对动量方程、热力学方程、水汽方程和连续方程进行时间差分,空间差分方案的变量分布水平方向采用Arakawa C跳点网格,垂直方向采用Charney/Phillips跳层。利用代数消元法进一步推导得到只包含未来时刻扰动Exner气压的亥姆霍兹方程,进而通过广义共轭余差法（GCR）求解,在此基础上得到未来时刻扰动量的预报值。基于所开发的扰动模式开展了数值试验。首先在非线性模式中施加一个中尺度初始扰动高压,得到初始扰动在非线性模式中的后续演变,然后将相同的初始扰动作为扰动模式的初值进行时间积分,将扰动模式预报的结果与非线性模式的结果做了对比。结果表明,所开发的扰动模式GRAPES_PF较好地模拟了惯性重力内波的传播过程:初始高压扰动激发了一个迅速向外传播的惯性重力内波,在气压场向风场适应的过程中,水平风场、垂直运动、位温和湿度等变量均出现了扰动增量,与非线性模式得到的结果相当接近。GRAPES_PF作为GRAPES非线性模式的合理线性模式为建立基于线性扰动预报的区域四维变分同化系统奠定了科学基础。 相似文献
9.
在非均匀分层下,目前GRAPES（Global/Regional Assimilation and Prediction System）模式中使用的垂直差分方案只能达到一阶精度。本文设计了一种适用于非均匀分层的二阶精度垂直差分方案,并将它应用于改进GRAPES模式动力框架的垂直离散化过程。一维廓线理想试验结果表明:二阶精度方案可以减少差分计算误差,而这种改进的幅度相对于差分计算本身引起的误差来说仍然是比较小的。通过密度流试验对修改后的模式动力框架进行测试,结果表明二阶方案可以保持模式动力框架的准确性和稳定性。进一步利用实际资料开展批量测试,发现二阶方案可以降低模式高空要素场的预报误差,而且这种改进随着预报时间的延长变得更为明显。最后选择一次典型的华南暴雨过程进行模拟,同样发现二阶精度方案对于48小时之后的降水会有一定程度的改进。 相似文献
10.
GRAPES区域集合预报系统模式不确定性的随机扰动技术研究 总被引：5，自引：2，他引：3
为进一步描述GRAPES(Global/Regional Assimilation and Prediction System)区域集合预报系统(GRAPES Me—soscale Ensemble Prediction System,GRAPES—MEPS)中GRAPES—Meso模式的不确定性特征,本研究在GRAPES—MEPS系统中引人了模式物理参数化倾向随机扰动方案(Stochastically Perturbed Parameterization Tendencies,SPPT),随机扰动型的产生是基于-阶马尔科夫链,其具有时间相关性特征,并服从正态分布,另外经过谱展开随机场具有空间结构特征,在水平结构上较平滑和连续。本文开展了基于SPPT方案的GRAPES—MEPS集合预报试验,针对SPPT方案中随机场的扰动幅度和时间相关尺度参数开展了一系列敏感性试验,并对试验结果进行了较全面的集合预报客观检验,此外,针对一次强降水过程,分析了SPPT方案对降水预报的影响。试验结果表明,引入SPPT方案能在一定程度上提高GRAPES—MEPS系统的预报技巧,降低系统的漏报率;且能显著改进预报后期大雨量级降水的预报技巧。通过敏感性试验发现,对于GRAPES—MEPS系统,SPPT方案的效果与随机扰动场幅度的范围,及扰动场的时间相关尺度选择相关,需经过敏感性试验确定出较适合的参数。 相似文献
11.
GRAPES伴随模式是其四维变分同化系统的核心组成部分。由于其计算过程复杂,临时数据较多,实现中采用断点存储策略可以有效减少伴随模式的计算时间和存储空间。极限断点存储策略是在单积分步内以全存储策略实现为基础,将其中部分基态以计算代替的一种类断点存储策略。在该策略的支持下,需要一种新的数据管理结构,来保证程序的正确运行。文章提出了在已有栈基础上优化的新数据存储管理方式——嵌套多链栈,这种结构可有效满足使用极限断点存储技术实现GRAPES伴随模式的初态管理需求。试验表明:相比断点存储技术,在总内存增加不超过30%的情况下可使GRAPES的运行效率提高1倍。 相似文献
12.
GRAPES非静力数值预报模式的三维变分资料同化系统的发展 总被引：21，自引：3，他引：18
为了减少分析变量与模式状态变量之间的插值误差,改善业务预报模式的初值质量,在GRAPES等压面三维变分资料同化系统的基础上,研究发展了针对GRAPES区域模式的非静力模式变量三维变分资料同化系统(GRAPES m3DVAR).该资料同化系统的垂直坐标及其分析变量的水平分布格式、垂直跳点方案与GRAPES预报模式保持完全一致.由于垂直坐标的变化和非静力关系,m3DVAR分析系统中设计了求解动力学约束方程的新方案.通过有效的高精度数学方案,避免了地形追随坐标下平衡方程的非线性项造成的复杂计算,有效解决了非静力平衡条件下求解平衡方程中非线性项的切线性方程和伴随方程引起的困难.重新构造各种观测算子,并考虑了质量场和风场之间的平衡约束关系、背景误差协方差结构,实现对探空、地面资料、船舶报等常规观测的同化.理想单点试验和实际资料的多变量资料同化分析结果表明,非静力模式变量三维变分资料同化系统能够正确地描写多变量之间的相互作用以及物理约束关系,分析结果合理,能够有效减少原等压面三维变分资料同化系统的分析与模式变量之间需要相互插值、变换产生的误差,在一定程度上提高了分析场质量,对预报模式的初值具有一定改善. 相似文献
13.
设计了一个以MM5模式为基础的遗传算法同化系统,并对一次暴雨过程进行了实际降水的模拟,通过对比遗传同化和伴随同化的降水预报效果,对遗传算法同化系统的同化性能进行验证。试验结果表明,遗传算法与四维变分相结合的同化系统能有效地改善模式的初始场,使MM5模式要素预报和降水预报的准确率得到提高。 相似文献
14.
15.
Four-dimensional variational data assimilation for mesoscale and storm-scale applications 总被引：1，自引：0，他引：1
Summary ?The status and progress of the four-dimensional variational data assimilation (4DVAR) are briefly reviewed focusing on application
to prediction of mesoscale/storm-scale atmospheric phenomena. Theoretical background is provided for each important component
of the 4DVAR system – forecast and adjoint models, observations, background, cost function, preconditioning, and minimization.
An overview of practical issues specific for mesoscale/storm-scale 4DVAR is then presented in terms of high-resolution observations,
nonlinearity and discontinuity problem, model error, errors from lateral boundary condition, and precipitation assimilation.
Practical strategies for efficient and simplified 4DVAR are also introduced, e.g., incremental 4DVAR, poor man’s 4DVAR, and
inverse 3DVAR. A new concept on hybrid approach is proposed to combine an efficient 4DVAR scheme and the standard 4DVAR scheme
aiming at reducing computational demand required by the standard 4DVAR while improving the accuracy of the simplified 4DVAR.
Applications to both hydrostatic and nonhydrostatic models are illustrated and our vision on opportunities and directions
for future research is provided.
Received March 12, 2001; revised July 24, 2001; accepted September 5, 2001 相似文献
16.
The variational assimilation theory is generally based on unbiased observations. In practice, however, almost all observations
suffer from biases arising from observational instruments, radiative transfer operator, precondition of data, and so on. Therefore,
a bias correction scheme is indispensable. The current scheme for radiance bias correction in the GRAPES 3DVar system is an
offline scheme. It is actually a static correction for the radiance bias before the process of cost function minimization.
In consideration of its effects on forecast results, this kind of scheme has some shortcomings. Thus, this study provides
a variational bias correction (VarBC) scheme for the GRAPES 3DVar system following Dee’s idea. In the VarBC scheme, the observation
operator is modified and a new control variable is defined by taking the predictor coefficients as the control parameters.
According to the feature of the GRAPES-3DVAR, an incremental formulation is applied and the original bias correction scheme
is maintained in the actual process of observations. The VarBC is designed to co-exist with the original scheme, because it
is a dynamic revision to the observational operator on the basis of the old method, i.e., it adjusts the model state vector
along with the control parameters to an unbiased state in the process of minimization and the assimilation system remains
consistent with available information automatically. Preliminary experimental results show that the mean departures of background-minus-observation
and analysis-minus-observation are reduced as expected. In a case study of the heavy rainfall that happened in South China
on 11–13 June 2008, the 500-hPa geopotential height is better simulated using the analyzed field from the VarBC as the initial
condition. 相似文献