一次强飑线过程极端大风的中尺度分析

利用常规气象观测资料、区域自动气象站加密观测资料和GFS 0.25°×0.25°逐6 h的分析场数据以及多普勒雷达、风廓线等资料,通过背景形势场分析、物理量诊断分析和中尺度分析,对2018年3月4日发生在华东地区的强飑线天气过程进行了诊断分析。结果表明,这次过程具有发生时间(季节)早、移速快、范围广、致灾强等特点,是一次比较少见的早春(冬末)十分强烈的飑线天气过程,是在高空急流辐散区、低空西南急流轴前端、低涡南侧的暖区中发展起来的。飑线过程的地面要素变化十分剧烈,地面有强冷池,与飑线前暖空气之间构成了强的水平温度梯度,致使飑线强度更强;飑线经过时气压涌升所形成的雷暴高压、强气压梯度以及飑线的快速移动均有利于地面极端大风的出现。飑线发展过程中观测到弓形回波、超级单体等强天气系统。中高层动量下传和光滑湖面、喇叭口、狭管效应等特殊地形对于大风的增强效应比较显著,这些因素也加剧了地面极端大风的形成。      

位场曲化平的插值-迭代法

将起伏曲面B上的位场向下延拓至曲面最低点的平面A的插值-迭代法步骤是：1）将曲面B上的场值放置在水平面A上具有相同水平坐标的点上,作为A上的初值;2）用若干水平面切割B,从A的初值,用快速傅里叶变换法(FFT)向上延拓出这些平面的场值,用插值的方法从这些平面的场值计算曲面B的场值;3）根据B上的实测值与计算值的差值,对A上的值进行加权改正;4）重复步骤2）和3）,直到B上的差值小到可以忽略.这种插值-迭代法具有高的计算速度,比通常的FFT法延拓得更深,可以超过10倍点距.文中给出计算实例.     

用平面波延拓方程进行地震数据的叠前速度反演

本文讨论地震勘探数据的叠前速度反演方法及其在海洋地震勘探数据上的反演试验.反演主要的计算步骤是:1.采用Fourier-Hankel变换把球面波分解为平面谐波;2.用平面谐波的延拓方程将上行波与下行波同时向下延拓,并计算每一层底部的反射系数和下一层的波阻抗;3.用最小二乘法从波阻抗中确定该层的声波速度.重复第2步与第3步,直到某一预定深度时结束.通过反演试验,对地震振幅比例的改变,子波变形,以及第1层速度和密度的误差对反演方法的稳定性及其精度的影响进行了分析.还通过实际海洋地震勘探数据的反演试验,对这一方法在地震勘探中的应用前景作了论述.     

卫星重力梯度向下延拓的频域最小二乘配置法

本文深入研究频域最小二乘配置法的基本原理及其在求解卫星重力梯度向下延拓问题中的应用。与一般空域最小二乘配置法相比，该方法具有高效稳定的特点，特别适合于大规模重力场数据处理，并为利用卫星重力梯度数据精化局部重力场提供了可供参考的方法。     

高低空急流在闽西北大暴雨过程中的作用及数值模拟

本文对2002年6月14~18日福建西北部地区的连续暴雨~大暴雨过程的诊断及模拟分析表明,在有利的环流背景条件下,高空急流的加强和动量下传,促使低空急流的建立和维持;低空急流的加强则提供了有利暴雨产生的高温高湿的大气条件;高低空急流的适宜配置,产生了动力场的耦合作用,为大暴雨的发生、发展提供了动力条件。     

A new technique to determine geoid and orthometric heights from satellite positioning and geopotential numbers

This paper takes advantage of space-technique-derived positions on the Earth’s surface and the known normal gravity field to determine the height anomaly from geopotential numbers. A new method is also presented to downward-continue the height anomaly to the geoid height. The orthometric height is determined as the difference between the geodetic (ellipsoidal) height derived by space-geodetic techniques and the geoid height. It is shown that, due to the very high correlation between the geodetic height and the computed geoid height, the error of the orthometric height determined by this method is usually much smaller than that provided by standard GPS/levelling. Also included is a practical formula to correct the Helmert orthometric height by adding two correction terms: a topographic roughness term and a correction term for lateral topographic mass–density variations.     

The topographic bias by analytical continuation in physical geodesy

This study emphasizes that the harmonic downward continuation of an external representation of the Earth’s gravity potential to sea level through the topographic masses implies a topographic bias. It is shown that the bias is only dependent on the topographic density along the geocentric radius at the computation point. The bias corresponds to the combined topographic geoid effect, i.e., the sum of the direct and indirect topographic effects. For a laterally variable topographic density function, the combined geoid effect is proportional to terms of powers two and three of the topographic height, while all higher order terms vanish. The result is useful in geoid determination by analytical continuation, e.g., from an Earth gravity model, Stokes’s formula or a combination thereof.     

The ellipsoidal corrections to the topographic geoid effects

The topographic effects by Stokes formula are typically considered for a spherical approximation of sea level. For more precise determination of the geoid, sea level is better approximated by an ellipsoid, which justifies the consideration of the ellipsoidal corrections of topographic effects for improved geoid solutions. The aim of this study is to estimate the ellipsoidal effects of the combined topographic correction (direct plus indirect topographic effects) and the downward continuation effect. It is concluded that the ellipsoidal correction to the combined topographic effect on the geoid height is far less than 1 mm. On the contrary, the ellipsoidal correction to the effect of downward continuation of gravity anomaly to sea level may be significant at the 1-cm level in mountainous regions. Nevertheless, if Stokes formula is modified and the integration of gravity anomalies is limited to a cap of a few degrees radius around the computation point, nor this effect is likely to be significant.AcknowledgementsThe author is grateful for constructive remarks by J Ågren and the three reviewers.     

雷暴大风过程中对流层中低层动量通量和动能通量输送特征研究

2014年5月31日北京发生一次雷暴大风过程。以雷达资料同化结果为初始场,对此次过程进行高分辨率数值模拟。采用非静力平衡和非地转平衡的经向动量方程和质量权重动能方程,利用模拟资料,对雷暴大风过程中经向动量和质量权重动能进行收支分析,以此来研究雷暴过程中对流层中低层动量通量和动能通量输送特征,讨论地面大风的可能成因。分析结果表明,在对流层中低层,经向动量通量散度是影响经向动量局地变化的主要强迫项。雷暴系统后部的入流把中低层的经向动量倾斜向下输送,系统前部对流云区中低层的下沉气流也向下输送经向动量。这两支下传动量通量先后与近地面经向动量的水平通量汇合,向系统前沿输送经向动量。在北京西北部地形阻挡作用下,经向动量通量在系统前端近地面辐合,促进那里的经向动量局地增长,有利于增强那里的南风。质量动能收支的特征与经向动量收支类似,在近地面层质量动能的局地变化主要是由质量动能通量散度引起的。系统后部入流把中层质量动能向下传输到近地面层,然后与近地面质量动能的水平通量汇合,向系统前沿输送质量动能。相对来说,近地面层经向动量和质量动能的水平通量比下传通量更重要,这主要与低层较强的东南急流有关。     

位场各阶垂向导数换算的新正则化方法

位场垂向导数大量应用于位场数据处理与解释中.当前广泛采用的位场各阶垂向导数换算方法为基于Laplace方程并结合波数域和空间域方法的具有递推特性的ISVD(integrated second vertical derivative)算法.本文在位场垂向导数换算的正则化方法和径向平均功率谱的基础上,提出一种位场各阶垂向导数换算的新正则化方法.新正则化方法仅需通过分析位场径向平均功率谱来确定一个截止波数,即可稳定换算位场各阶垂向导数.理论模型和实测数据实验结果表明:(1)新正则化方法物理意义明确、计算简单,且各阶垂向导数换算的稳定性和精度明显优于ISVD算法;(2)在用新正则化方法求得各阶垂向导数的基础上,利用泰勒级数法可以获得大深度、高精度的位场向下延拓结果.