共查询到18条相似文献，搜索用时 178 毫秒
1.
基于总能量模的奇异向量扰动常用于构造集合预报的初始条件。以建立GRAPES（Global and Regional Assimilation PrEdiction System）全球集合预报系统为目的,基于前期研发的GRAPES全球模式奇异向量方法,在GRAPES全球切线性模式和伴随模式2.0版的框架下,开展了引入线性化边界层方案来改善奇异向量结构,并提高奇异向量计算效率的研究。通过连续试验,从奇异向量的扰动能量结构、扰动能量谱及扰动空间分布等方面,综合分析改进GRAPES全球奇异向量的结构及演变特征。试验结果表明,改进后的GRAPES奇异向量方法有效抑制了之前扰动能量在近地面层不合理的快速增长,同时,奇异向量最优扰动的结构更客观地体现了中高纬度区域大气初始条件中的斜压不稳定扰动及其演变,如在初始时刻奇异向量扰动能量主要位于对流层中层,并呈现出随高度向西倾斜的大气斜压特征;经过线性化演变,扰动能量向较大水平尺度转移,并在垂直结构上表现出向对流层高层上传及向对流层低层下传的特征等。针对GRAPES奇异向量迭代求解中伴随模式计算耗时为主的情况,改进伴随模式中广义共轭余差方案的调用方式,并采用大内存存储法来提高其计算效率,进而将奇异向量总计算时间缩短了25%。总之,改进后的GRAPES奇异向量方法,可应用于构建面向业务应用的GRAPES全球集合预报系统。 相似文献
2.
集合数值天气预报的关键问题就是如何生成有效的初始扰动。奇异向量反映了初始扰动在大气系统相空间中演变发展的最不稳定方向,基于奇异向量产生的集合样本是模拟概率密度函数的最合理方法。以非静力、半隐式半拉格朗日GRAPES-Meso中尺度数值预报模式为基础,采用Lanczos迭代算法,利用GRAPES-Meso的切线伴随模式计算GRAPES奇异向量。为了检验求得的奇异向量的正确性,提出了两种检验奇异向量正确性的方法:一是比较计算的奇异值的一致性;二是依据特征向量在矩阵变换中的方向不变性意义,验证GRAPES奇异向量空间结构的正确性。最后研究了不同的时间间隔对GRAPES奇异向量的影响,结果表明GRAPES奇异向量在36小时的最优时间间隔误差增长速度最快,这表明在非静力、半隐式半拉格朗日格点模式中利用切线伴随技术计算奇异向量是可行的。 相似文献
3.
目前国际上采用的奇异向量集合预报初值扰动法对于初值不确定性的描述存在一定的不足,为了更有效地反映初始误差的时空多尺度特性,基于GRAPES全球奇异向量计算技术,计算了不同空间分辨率及不同最优时间间隔的多个尺度的奇异向量,并采用基于高斯分布的线性组合法来构造多尺度奇异向量的扰动初值,以代表在相空间中增长最快的多尺度初值误差模态。通过2019年1月19日的初值扰动集合预报试验,对比分析了单一尺度奇异向量初值扰动法与多尺度初值扰动法的扰动特征及集合预报效果。结果表明,多尺度奇异向量初值扰动法为区域集合预报提供的初始扰动场是合理的,扰动的大小随时间增长,且在空间分布上较好地反映了当前大气的斜压不稳定特征。此外,多尺度奇异向量扰动可以描述一定的大尺度以及中小尺度运动误差特征,较单一尺度奇异向量扰动能反映出更多初始场的不确定性信息。检验分析表明,GRAPES多尺度奇异向量集合预报在集合一致性、连续等级概率评分、离群值等方面有一定的优势,相比于单一尺度奇异向量法有较好的预报技巧。因此,基于GRAPES的多尺度奇异向量初值扰动法对于集合预报的预报效果有一定的提高,能为构建一套完善的GRAPES区域奇异向量集合预报系统提供一定的科学依据和应用基础。 相似文献
4.
采用线性化物理过程方案的GRAPES全球模式奇异向量在进行非线性模式积分时会有部分奇异向量出现崩溃问题,这说明奇异向量结构可能存在扰动变量之间不协调之处,需要对奇异向量扰动的计算方法优化,进而改进基于奇异向量的集合预报初值扰动,提高GRAPES全球集合预报效果。基于原有的GRAEPS全球奇异向量计算方法,在求解奇异向量时,对气压扰动的处理进行改进,将初始时刻的气压扰动分量通过位温扰动根据静力平衡关系导出获得,其他保持一致,发展了静力平衡奇异向量改进方法。基于有两个台风过程的个例（2019年8月8日12时（世界时））,分别采用原奇异向量方法和静力平衡奇异向量改进方法进行热带气旋目标区奇异向量的计算求解,并进行相应奇异向量的非线性模式积分,对比分析奇异向量非线性积分的稳定性。进而,对比分析奇异向量求解方法改进前、后热带气旋奇异向量的结构特征和初值扰动特征,开展了集合预报试验,评估改进后的奇异向量求解方法对GRAPES全球集合预报系统预报性能的影响。试验结果表明,静力平衡奇异向量改进方法通过产生协调的气压扰动和位温扰动场,解决了奇异向量非线性积分崩溃的问题,消除了原来不利于积分稳定性的气压扰动过于局地化的小尺度结构。静力平衡奇异向量改进方法对奇异向量中位温扰动分量和纬向风扰动分量结构影响较小,使得气压扰动分量的大值区位于台风附近,更好地描述热带气旋初值不确定性,与位温扰动分量的分布更加协调。采用静力平衡奇异向量改进方法,可以提高GRAPES全球集合预报在北半球和南半球等压面要素集合预报技巧和中国地区24 h累计降水概率预报技巧,增大台风路径集合离散度。 相似文献
5.
伴随技术是四维变分同化（4DVar）系统中计算代价函数梯度的最佳办法,切线性和伴随模式的效果和效率直接影响着4DVar系统的发展。基于GRAPES（Global and Regional Assimilation PrEdiction System）全球切线性和伴随模式1.0版本,利用GRAPES全球模式2.0版本在并行框架和性能等方面的改善,重新优化和设计了GRAPES全球切线性伴随模式2.0版本,提高了GRAPES全球切线性和伴随模式的效果和效率,优化了切线性模式程序结构,使其计算时间最优可控制在非线性模式的1.2倍以内;采用在切线性模式中保存基态的方法,重构了伴随模式的程序结构,使其计算时间最优控制在非线性模式的1.5倍以内;在GRAPES全球切线性物理过程的设计中,将线性物理过程的轨迹基态计算和切线性扰动计算解耦,提高了GRAPES全球切线性和伴随模式的计算效果和效率。 相似文献
6.
针对GRAPES四维变分同化系统的升级, 研究了GRAPES模式垂直扩散方案线性化问题。通过2005年8月7—27日21个个例的批量试验, 发现在GRAPES模式垂直扩散方案源代码的基础上逐句线性化得到的切线性垂直扩散方案即使能通过正确性测试试验, 在少数情况下也会存在很大误差。切线性模式计算的扰动气压场和扰动风场可能出现明显异常, 这种异常与垂直扩散方案中地表动量通量的强非线性有关。如果在切线性垂直扩散方案中忽略地表动量通量扰动, 既可以避免异常的出现, 又不影响其他正常时刻的计算精度。修改后的切线性垂直扩散方案能够在所有变量上一致地提高切线性模式的计算精度。 相似文献
7.
线性化物理过程能够改善四维变分同化中极小化收敛的稳定性和增加极小化过程中对大气物理过程和动力更加精确的描述,它是四维变分同化中非常重要的一部分。通过在GRAPES全球模式中研究线性化物理过程,尤其是两个湿线性化物理过程,改善切线性模式预报精度,来提高GRAPES全球四维变分同化的分析和预报效果。线性化物理过程的开发首先需要简化原非线性化物理过程中的强非线性项,然后对线性化物理过程进行规约化,以抑制切线性扰动的异常增长。目前GRAEPS全球模式中的线性化物理过程主要包括次网格尺度地形参数化、垂直扩散、积云深对流和大尺度凝结。线性化物理过程预报精度的检验方法是通过选择合适大小的初始扰动（同化分析增量）,来比较非线性模式和切线性模式中的扰动演化的纬向平均误差。然后以绝热版本的切线性模式为基础,通过冬、夏两个个例试验来分别检验4个线性化物理过程的12 h预报效果。试验结果表明,通过添加次网格地形参数化和垂直扩散两个干线性化物理过程方案,可以有效抑制住绝热版本切线性模式低层扰动的异常增长,大幅度改善切线性模式预报效果。通过添加积云深对流和大尺度凝结两个湿线性化物理过程,可以在热带区域和中、高纬度地区提高切线性模式中湿变量和温度变量的近似精度,提高GRAPES全球四维变分同化的分析和预报效果。 相似文献
8.
为描述GRAPES全球模式初始条件的不确定性,基于适合集合预报应用的GRAPES全球奇异向量技术,依据大气初始误差符合正态分布的特征,采用高斯取样奇异向量来构造全球集合预报初始扰动,在此基础上建立了GRAPES全球集合预报系统(GRAPES-GEPS)。利用GRAPES全球同化分析场,对采用初始扰动的GRAPES-GEPS连续试验预报结果进行检验和分析。结果表明:GRAPES-GEPS中高度场、风场及温度场预报的集合离散度能有效快速增加,集合平均均方根误差与集合离散度的关系合理;相对控制预报的均方根误差,集合平均的预报优势在预报中期非常显著。为进一步体现GRAPES-GEPS中模式物理过程的不确定性,发展了模式物理过程倾向随机扰动技术(SPPT)。试验结果表明:SPPT方案的应用有效提高了GRAPES-GEPS在南、北半球和热带地区等压面要素预报的集合离散度,同时一定程度减小了集合平均误差,进而改进了集合平均误差与集合离散度的关系,其中SPPT方案在热带地区的改进最为显著。本文发展的基于奇异向量的初始扰动方法和模式扰动SPPT方案在中国气象局2018年12月业务化运行的GRAPES-GEPS中得到了应用。 相似文献
9.
一种求解条件非线性最优扰动的快速算法及其在台风目标观测中的初步检验 总被引：2，自引：0，他引：2
条件非线性最优扰动(CNOP)是Mu等2003年提出的一个新的理论方法,它是线性奇异向量在非线性情形的推广,克服了线性奇异向量不能代表非线性系统最快发展扰动的缺陷,成为非线性系统可预报性和敏感性等研究新的有效工具.然而,由于以往CNOP的求解需要采用伴随技术,计算量相当巨大,限制了该方法的推广应用.为了克服这一困难,本文基于经验正交分解(EOF),提出了一种求解CNOP的快速算法,利用GRAPES区域业务预报模式实现了CNOP快速计算,并在台风"麦莎"的目标观测研究中得到初步检验,通过观测系统模拟实验(OSSE)检验了该方法确定敏感性区域(瞄准区)的有效性和可行性.试验结果表明,用快速算法求解的CNOP,其净能量随时间快速地发展,而且发展呈非线性.在台风"麦莎"个例的目标观测试验中,用快速算法得到的预报时间为24 h的CNOP可以有效地识别瞄准区,并通过瞄准区内初值的改善,可明显减少目标区域(检验区)内24 h累计降水预报误差.尤其,累计降水预报的这种改进效果能够延伸到更长时间(如72 h),尽管检验时间是设在第24小时.进一步分析发现,24 h累计降水预报误差的减少是通过利用瞄准区内改善的初值改进初始时刻台风暖心结构、高空相对涡度以及水汽条件等而得以实现的. 相似文献
10.
基于副热带奇异向量的初值扰动方法已应用于GRAPES （Global and Regional Assimilation PrEdiction System）全球集合预报系统,但存在热带气旋预报路径离散度不足的问题。通过分析发现,热带气旋附近区域初值扰动结构不合理导致预报集合不能较好地估计热带气旋预报的不确定性,是路径集合离散度不足的可能原因之一。通过建立热带气旋奇异向量求解方案,将热带气旋奇异向量和副热带奇异向量共同线性组合生成初值扰动,以弥补热带气旋区域初值扰动结构不合理这一缺陷,进而改进热带气旋集合预报效果。利用GRAPES全球奇异向量计算方案,以台风中心10个经纬度区域为目标区构建热带气旋奇异向量求解方案,针对台风“榕树”个例进行集合预报试验,并开展批量试验,利用中国中央气象台最优台风路径和中国国家气象信息中心的降水观测资料进行检验,对比分析热带气旋奇异向量结构特征和初值扰动特征,评估热带气旋奇异向量对热带气旋路径集合预报和中国区域24 h累计降水概率预报技巧的影响。结果表明,热带气旋奇异向量具有局地化特征,使用热带气旋奇异向量之后,热带气旋路径离散度增加,路径集合平均预报误差和离散度的关系得到改善,路径集合平均预报误差有所减小,集合成员更好地描述了热带气旋路径的预报不确定性;中国台风降水的小雨、中雨、大雨、暴雨各量级24 h累计降水概率预报技巧均有一定提高。总之,当在初值扰动的生成中考虑热带气旋奇异向量后,可改进热带气旋初值扰动结果,并有助于改善热带气旋路径集合预报效果。 相似文献
11.
In this paper,a nonlinear optimization method is used to explore the finite-time instability of the atmospheric circulation with a three-level quasigeostrophic model under the framework of the conditional nonlinear optimal perturbation (CNOP).As a natural generalization of linear singular vector (SV),CNOP is defined as an initial perturbation that makes the cost function the maximum at a prescribed forecast time under certain physical constraint conditions.Special attentions are paid to the different structures and energy evolutions of the optimal perturbations.The results show that the most instable region of the global atmospheric circulation lies in the midlatitude Eurasian continent.More specially,SV and CNOP in the total energy norm with an optimization time of 2 days both present localness:they are mainly located in the midlatitude Asian continent and its east coast.With extension of the optimization time,SVs are more upstream and less localized in the zonal direction,and CNOPs differ essentially from SVs with broader zonal and meridional coverages; as a result,CNOPs acquire larger kinetic and available potential energy amplifications than SVs in the nonlinear model at the corresponding optimization time.For the climatological wintertime flow,it is seen that the baroclinic terms remain small over the entire time evolution,and the energy production comes essentially from the eddy kinetic energy,which is induced by the horizontal shear of the basic flow.In addition,the effects of SVs and CNOPs on the Eurasian atmospheric circulation are explored.The results show that the weather systems over the Eurasian continent in the perturbed fields by CNOPs are stronger than those by SVs at the optimization time.This reveals that the CNOP method is better in evaluating the instability of the atmospheric circulation while the SV method underestimates the possibility of extreme weather events. 相似文献
12.
The singular vector (SV) initial perturbation method can capture the fastest-growing initial perturbation in a tangent linear model (TLM). Based on the global tangent linear and adjoint model of GRAPES-GEPS (Global/Regional Assimilation and Prediction System—Global Ensemble Prediction System), some experiments were carried out to analyze the structure of the moist SVs from the perspectives of the energy norm, energy spectrum, and vertical structure. The conclusions are as follows: The evolution of the SVs is synchronous with that of the atmospheric circulation, which is flow-dependent. The moist and dry SVs are located in unstable regions at mid-to-high latitudes, but the moist SVs are wider, can contain more small- and medium-scale information, and have more energy than the dry SVs. From the energy spectrum analysis, the energy growth caused by the moist SVs is reflected in the relatively small-scale weather system. In addition, moist SVs can generate perturbations associated with large-scale condensation and precipitation, which is not true for dry SVs. For the ensemble forecasts, the average anomaly correlation coefficient of large-scale circulation is better for the forecast based on moist SVs in the Northern Hemisphere, and the low-level variables forecasted by the moist SVs are also improved, especially in the first 72 h. In addition, the moist SVs respond better to short-term precipitation according to statistical precipitation scores based on 10 cases. The inclusion of the large-scale condensation process in the calculation of SVs can improve the short-term weather prediction effectively. 相似文献
13.
Previous studies have shown that Regional Climate Models (RCM) internal variability (IV) fluctuates in time depending on synoptic events. This study focuses on the physical understanding of episodes with rapid growth of IV. An ensemble of 21 simulations, differing only in their initial conditions, was run over North America using version 5 of the Canadian RCM (CRCM). The IV is quantified in terms of energy of CRCM perturbations with respect to a reference simulation. The working hypothesis is that IV is arising through rapidly growing perturbations developed in dynamically unstable regions. If indeed IV is triggered by the growth of unstable perturbations, a large proportion of the CRCM perturbations must project onto the most unstable singular vectors (SVs). A set of ten SVs was computed to identify the orthogonal set of perturbations that provide the maximum growth with respect to the dry total-energy norm during the course of the CRCM ensemble of simulations. CRCM perturbations were then projected onto the subspace of SVs. The analysis of one episode of rapid growth of IV is presented in detail. It is shown that a large part of the IV growth is explained by initially small-amplitude unstable perturbations represented by the ten leading SVs, the SV subspace accounting for over 70% of the CRCM IV growth in 36?h. The projection on the leading SV at final time is greater than the projection on the remaining SVs and there is a high similarity between the CRCM perturbations and the leading SV after 24–36?h tangent-linear model integration. The vertical structure of perturbations revealed that the baroclinic conversion is the dominant process in IV growth for this particular episode. 相似文献
14.
大气对外强迫低频遥响应的数值模拟 II:对欧亚中高纬“寒潮”异常的响应 总被引：6，自引：2，他引：4
本文通过用IAP GCM所作的数值模拟研究了欧亚大陆中高纬度地区的外强迫在全球大气中激发的响应.结果清楚表明,同赤道地区的热源强迫一样,中高纬度地区的外强迫也可以在全球大气中产生低频遥响应;通过低频波列EAP,欧亚大陆中高纬度地区的寒潮异常可以对赤道中西太平洋地区的大气运动及全球大气环流有重要影响.对30—60天振荡的强迫激发来讲,地球大气的气候基本态是极为重要的,热带大气对于全球大气的低频振荡活动具有尤为突出的作用. 相似文献
15.
The effect of barotropic shear in the basic flow on baroclinic instability is investigated using a linear multilevel quasi-geostrophic β-plane channel model and a nonlinear spherical primitive equation model. Barotropic shear has a profound effect on baroclinic instability. It reduces the growth rates of normal modes by severely restricting their structure, confirming earlier results with a two-layer model. Dissipation, in the form of Ekman pumping and Newtonian cooling, does not change the main characteristics of the effect of the shear on normal mode instability.Barotropic shear in the basic state, characterized by large shear vorticity with small horizontal curvature, also effects the nonlinear development of baroclinic waves. The shear limits the energy conversion from the zonal available potential energy to eddy energy, reducing the maximum eddy kinetic energy level reached by baroclinic waves. Barotropic shear, which controls the level of eddy activity, is a major factor which should be considered when parameterizing the eddy temperature and momentum fluxes induced by baroclinic waves in a climate model. 相似文献
16.
全球变暖减缓期陆地地表气温变化特征和CMIP5多模式的未来情景预估 总被引：3，自引：1，他引：2
2000年后全球气温的增温率显著下降,全球进入变暖减缓期.本文基于CRU(Climatic Research Unit) 观测资料,分析讨论了2000年后全球及欧亚中高纬度地区全球变暖的减缓特征,评估了CMIP5(Coupled Model Intercomparison Project Phase 5)试验多模式对全球变暖减缓的模拟及未来气温变化预估.结果表明,2000年后全球陆地平均地面气温的增温率大幅下降至0.14℃ (10 a)-1,仅为1976~1999年加速期增温率的一半.全球陆地13个区域中有9个地区的增温率小于2000年前,4个地区甚至出现了降温.其中以欧亚中高纬地区最为特殊.加速期(1976~1999年)增温率达到0.50℃ (10 a)-1,为全球陆地最大,2000年后陡降至-0.17℃ (10 a)-1,为全球最强降温区,为全球变暖的减缓贡献了49.13%.并且具有显著的季节依赖,减缓期冬季增温率下降了-2.68℃ (10 a)-1,而秋季升高了0.86℃ (10 a)-1,呈现反位相变化特征.CMIP5多模式计划中仅BCC-CSM1.1在RCP2.6情景下和MRI-ESM1模式在RCP8.5下的模拟较好地预估了全球及欧亚中高纬地区在2000年后增温率的下降以及欧亚中高纬秋、冬温度的反位相变化特征.BCC-CSM1.1在RCP2.6情景下预估欧亚中高纬地区2012年后温度距平保持在1.2℃左右,2020年后跃至2℃附近振荡.而MRI-ESM1在RCP8.5情景下预估的欧亚中高纬度温度在2030年前一直维持几乎为零的增温率,之后迅速升高. 相似文献
17.
In this study,singular vectors related to a heavy rainfall case over the Korean Peninsula were calculated using the fifth-generation Pennsylvania State University-National Center for Atmospheric Research Mesoscale Model (MM5) adjoint modeling system.Tangent linear and adjoint models include moist physical processes,and a moist basic state and a moist total energy norm were used for the singular-vector calculations.The characteristics and nonlinear growth of the first singular vector were analyzed,focusing on the relationship between the basic state and the singular vector.The horizontal distribution of the initial singular vector was closely related to the baroclinicity index and the moisture availability of the basic state.The temperature-component energy at a lower level was dominant at the initial time,and the kinetic energy at upper levels became dominant at the final time in the energy profile of the singular vector.The nonlinear growth of the singular vector appropriately reflects the temporal variations in the basic state.The moisture-component energy at lower levels was dominant at earlier times,indicating continuous moisture transport in the basic state.There were a large amount of precipitation and corresponding latent heat release after that period because the continuous moisture transport created favorable conditions for both convective and nonconvective precipitation.The vertical propagation of the singular-vector energy was caused by precipitation and the corresponding latent heating in the basic state. 相似文献
18.
Sei-Young Park Hyun Mee Kim Tae-Young Lee Michael C. Morgan 《Meteorology and Atmospheric Physics》2013,120(3-4):107-122
To investigate the statistical sensitivity distributions of tropical cyclone (TC) forecasts over the Korean Peninsula, total energy (TE) singular vectors (SVs) were calculated and evaluated over a 10-year period. TESVs were calculated using the fifth-generation Pennsylvania State University/National Center for Atmospheric Research Mesoscale Model (MM5) and its tangent linear and adjoint models with a Lanczos algorithm over a 48-h period. Chosen cases were 21 TCs that affected the Korean Peninsula among 230 TCs that were generated in the western North Pacific from 2001 to 2010. Sensitive regions indicated by TESVs were mainly located near mid-latitude troughs and TC centers but varied depending on TC track and environmental conditions such as subtropical high and mid-latitude trough. The cases were classified into three groups by clustering TC tracks based on the finite mixture model. The two groups that passed through the western and southern sea of the Korean Peninsula had maximally sensitive regions in the mid-latitude trough and largely sensitive regions around the TC center, while the other group that passed straight through the eastern sea of the Korean Peninsula had maximally sensitive regions near the northeastern region of the TC center. Vertically, the former two clustered groups had the westerly tilted TESVs and potential vorticity structures under the mid-latitude troughs at the initial time, indicating the TCs were in a baroclinic environment. Conversely, the straight-moving TCs were not in a baroclinic environment. Based on the results in the present study, the TCs moving toward a fixed verification region over the Korean Peninsula have different sensitivity regions and structures according to their moving tracks and characteristic environmental conditions, which may provide guidance for targeted observations of TCs affecting the Korean Peninsula. 相似文献