Analysis of the Impact of Coupled Gravity Waves on the Low-altitude Catastrophic Gust at Kunming Changshui International Airport

  • Wendong HU , 1, 2 ,
  • Yu PENG 1, 2 ,
  • Tiangui XIAO 1, 2 ,
  • Wenjie ZHOU 1 ,
  • Yixue DENG 3 ,
  • Jingyi TAO 4
Expand
  • 1. College of Atmospheric Sciences/ Sichuan Key Laboratory of Plateau Atmosphere and Environment/ Sichuan Engineering Laboratory for Meteorological Disaster Prediction and Warning,Chengdu University of Information Technology,Chengdu 610225,Sichuan,China
  • 2. Yunnan Research and Development Center for Natural Disaster Prevention Technology,Chengdu University of Information Technology,Kunming 650034,Yunnan,China
  • 3. Yunnan Provincial Meteorological Observatory,Kunming 650034,Yunnan,China
  • 4. College of Earth Sciences,Yunnan University,Kunming 650091,Yunnan,China

Received date: 2025-02-20

  Revised date: 2025-07-29

  Online published: 2026-05-22

Copyright

© Editorial Department of Plateau Meteorology (CC BY-NC-ND)

Abstract

On March 27, 2024, a disastrous gust damaged the terminal building at Kunming Changshui International Airport in Yunnan, China.The role of coupled atmospheric gravity waves in this event is examined using conventional surface meteorological observations, high-resolution mesoscale surface network, upper-air sounding data, ERA5 reanalysis data and FY-2G geostationary satellite imagery.The results show that: (1) After solar heating in the morning, an atmospheric gravity wave originating over the Hengduan Mountains propagated southeastward under the influence of westerly flow.Meanwhile, another gravity wave generated over the Ailao Mountains propagated northeastward.These two gravity waves met near 100.56°E, 25°N about 09:53(Beijing Time), producing a significant coupling effect.(2) Based on the force analysis within a statically stable atmospheric layer and the initial perturbations induced by orographic lifting of the background flow, the analytical mathematics formula of the gravity wave was derived.Using satellite observations and infrared cloud imagery, the key parameters of each gravity waves including significant amplitude, wavelength, intrinsic frequency, and phase speed-were quantitatively obtained.(3) By harmonic analysis, the coupled wave with its zonal and meridional components were determined.The coupled wave induced substantial vertical motions above the airport: the zonal component generated a downward motion of 18 m·s-1, contributing positively to the surface winds, while the meridional component produced a weaker upward motion of 1 m·s-1, exerting a negative contribution.(4) The comprehensive effect of the coupled gravity waves and the DCAPE releasing generated a strong subsiding airflow, which enhanced the downward transport of upper-level jet momentum.This momentum transition was combined with the near-surface horizontal wind field, ultimately resulting in a catastrophic surface gust reaching 86.5 m·s-1, causing severe structural damage to the main terminal building of Kunming Airport.Quantitative assessment indicates that gravity waves contributed approximately 19.65% to the total wind speed.In contrast, the presence of either gravity wave A or B alone without coupling effect would account for a maximum wind speed contribution of only about 10 m·s-1.The role of gravity wave coupling is of critical in intensifying subsidence and promoting the extreme low-level wind hazards.

Cite this article

Wendong HU , Yu PENG , Tiangui XIAO , Wenjie ZHOU , Yixue DENG , Jingyi TAO . Analysis of the Impact of Coupled Gravity Waves on the Low-altitude Catastrophic Gust at Kunming Changshui International Airport[J]. Plateau Meteorology, 2026 , 45(3) : 678 -690 . DOI: 10.7522/j.issn.1000-0534.2025.00091

1 引言

低空经济是以航空载运与作业装备技术为主要工具、 以低空空域为主要活动场域、 以低空飞行活动为最终产出形式的系列经济活动领域, 正展现出巨大的发展潜力和广阔的市场前景(覃睿, 2023)。
尽管低空经济前景光明, 但却面临着诸多挑战, 其中航空气象保障的欠缺尤为突出。在大气科学发展的历史进程中, 人类一直关注地表的气象条件, 如1.5 m气温、 10 m风向风速等等, 同时积累了高空如500 hPa、 700 hPa、 850 hPa的相关预报分析经验。但对低空, 特别是地面之上300 m左右的低空开放区域, 却并无成熟可靠的预报能力, 低空气象保障存在很大不足。因此, 加强航空气象保障研究, 提升低空飞行的安全性和效率, 是当前亟待解决的重要课题(胡壮等, 2022)。近年来对大风的研究不断深入, 如杨显玉等(2023)张珊等(2023)李叶晴等(2023)李晨轩和韦志刚(2024)付文卓等(2024)马文倩等(2025)徐景峰等(2025)杨坤琳等(2025)董泽新等(2025)分别对大风的时空特征、 气候变化、 复杂情况的风场模拟进行了分析研究。特别是在航空气象方面, 鲍艳松等(2024)杨坤琳等(2025)以及董泽新等(2025)在大风相关保障领域取得了新进展。
在众多影响低空飞行的气象因素中, 重力波作为一种广泛存在的大气现象往往重视不够, 但其对飞行安全的影响却不容忽视。重力波是层结稳定大气在受到垂直方向扰动后产生的一种周期性振荡运动, 这种运动不仅在大气中传播, 还通过动量的垂直输送促进不同大气层之间物质和能量的交换, 对大气中的动量和能量平衡起着重要影响(魏栋等, 2016)。重力波可以传播很长的距离, 将大量的动量和能量在很大的范围内转移, 从而影响广大区域(Houghton, 1978Lindzen, 1981Holton et al, 1995Richter et al, 2010Limpasuvan et al, 2012Butchart, 2014)。
重力波可由多种机制产生(Fritts and Alexander, 2003), 对流、 地形、 与锋面和急流相关的流场不平衡、 切变不稳定以及日间的海陆温差等, 均可产生不同特性和机制的重力波(孙淑清, 1983Durran, 1990陆汉城等, 1992Lane and Clark, 2002林永辉等, 2003Plougonven and Zhang, 2014Du and Rotunno, 2015Lane, 2015Wang et al, 2018Fang and Du, 2022Yang et al, 20232024)。尤其在复杂地形下, 如气流经过山脉而激发出重力波的情况十分常见(Malkus and Stern, 1953; Tjernstrom and Mauritsen, 2009; 张云等, 2011谢家旭和李国平, 2021)。
重力波在对流等灾害性天气中发挥着重要作用。在特定条件下, 重力波能触发积云对流, 进而引发强对流天气, 产生暴雨、 冰雹以及大风等极端天气现象(李麦村, 1978Clark et al, 1986Adams-Selin and Johnson, 2013覃卫坚等, 2013崔新艳等, 2021)。目前, 关于重力波对极端天气的影响分析主要集中在暴雨和冰雹方面, 如郭虎等(2006)发现重力波会影响暴雨的发生和时空分布。孙继松等(2009)研究表明, 地形重力波能够造成冰雹直径增加, 从而增强冰雹的破坏力。王文等(2011)研究表明, 重力波与对流相互作用加强暴雨的强度, 影响暴雨的空间分布。吴迪等(2016)发现重力波是中尺度强降水的发生机制之一。Liu et al(2018)的研究表明, 重力波与降水之间存在正反馈机制。汪科均等(2024)提出包含重力波在内的中尺度波动在夜间山地暴雨过程中, 各物理量场表现出特定特征。
近年来关于重力波对大风天气的影响研究取得了一些成果, 如孙艳辉等(20122015)对暴雪、 大风过程的分析, 汤浩等(2019)对新疆大风过程的分析,但总体来看成果相对较少。
昆明长水国际机场位于中国西南高原地区, 周边山地众多, 地形对大气动力过程的作用显著。已有研究表明, 复杂地形条件下, 重力波是影响低空大风的重要因素(Smith, 1979; 卢冰等, 2014)。目前重力波如何影响低空大风的研究相对较为薄弱, 特别是重力波耦合效应在高原机场环境下, 这一领域存在诸多未知之处, 值得进一步探索。因此, 深入分析重力波对大风灾害的影响机制, 有助于提升航空气象保障能力, 为低空经济的安全发展提供有力支撑。

2 资料与方法

2.1 数据资料

分析采用的资料包括常规地面气象观测资料、 加密中尺度地面气象观测资料、 高空气象探测资料、 ERA5(欧洲中期天气预报中心第五代)0.25°× 0.25°逐时大气再分析资料, FY-2G静止气象卫星可见光(波长0.55~1.05 µm)和红外波段(波长10.3~11.5 µm)云图资料。

2.2 研究方法

运用天气学分析、 天气诊断分析等方法, 对环流形势和影响天气系统进行分析。运用卫星气象学方法, 对可见光与红外云图进行分析。
通过空气质点的受力分析, 结合初始条件推导出重力波运动方程。利用可见光云图的条带状重力波云系空间分布, 测量重力波的波长; 利用重力波云系的整体传播速度, 计算出重力波的频率、 周期; 通过红外云图的云顶亮温, 结合大气垂直结构与ERA5再分析资料, 分别计算重力波云系的云顶与云底高度, 得到重力波的显著振幅及全部参数。利用谐波分析方法, 研究两个重力波耦合后的波动及其在东西和南北两个方向的分量, 从而计算出重力波做功对大风造成的影响。

2.3 研究区域

研究区域为昆明长水机场及其上游横断山脉和哀牢山脉[图1, 此图及文中涉及的研究区域地图是基于国家地理信息公共服务平台网站审图号为GS(2025)1508号的标准地图制作, 底图无修改]。横断山脉位于云南西北地区, 平均海拔超4000 m, 为南北走向。哀牢山脉位于云南中部地区, 平均海拔2000 m以上, 主峰海拔3166 m, 为西北至东南走向。昆明长水国际机场位于102.93°E, 25.1°N, 海拔2102 m。机场处于两个山岭之间的狭长平原,即“坝子”区域, 航站楼在跑道西南端靠近昆明市区一侧, 是本次过程的机场前端迎风处。
图1 昆明长水国际机场与云南山脉地形

A: 横断山系, B: 哀牢山系, C: 昆明机场

Fig.1 Kunming Changshui International Airport and the terrain of Yunnan mountainsA: Hengduan Mountain, B: Ailao Mountain, C: Kunming Airport

3 灾害实况、 环流形势和影响系统

3.1 灾害实况

2024年3月27日11:50(北京时, 下同)左右, 昆明长水国际机场受罕见瞬时大风袭击。气象观测到地面10 m高度处阵风速度27 m·s-1, 风向为230°。大风导致昆明长水机场T1航站楼的金属屋顶局部被大风掀起, 部分区域受损(图2)。
图2 长水机场受损实况图

黑色方框表示受损区域

Fig.2 Damage to Changshui International Airport.The black rectangle represents the damaged area

3.2 环流形势和影响系统

2024年3月27日08:00, 500 hPa[图3(a)]亚欧大陆的环流为两脊一槽, 我国大范围地区受到一个显著的高空大槽影响。昆明位于槽前, 槽前有西南急流, 最大风速为28 m·s-1, 槽后为较强冷平流。08:00 -12:00, 500 hPa高空槽不断东移并持续发展。至12:00[图3(b)], 高空槽移至长水机场上空, 为此处天气系统发展提供了有利的条件。
图3 2024年3月27日500 hPa环流形势

(a)08:00, (b)12:00, 黑色三角表示长水机场

Fig.3 500 hPa circulation at 08:00 (a) and 12:00 (b) on 27 March 2024.The black triangles represent Kunming Changshui International Airport

700 hPa[图4(a)]有槽逼近长水机场, 槽前为暖干西南急流, 最大风速为24 m·s-1。槽后则伴随着强烈的冷平流。这种冷暖平流的交汇, 加剧了大气的不稳定性。根据槽的位置可知, 12:00长水机场上空500 hPa高度槽超前于700 hPa[图4(b)]高度槽, 是前倾槽空间配置, 即高空的冷空气已经移到长水机场附近, 而低空仍处于低空槽前暖空气之中, 这种上冷下暖前倾槽的热力配置非常有利于对流的发展。昆明上空500 hPa大气非常干燥, 700 hPa露点温度差为7 ℃, 大气为中等湿度状态。受到水汽条件限制, 该区域不会有深厚湿对流发展, 不会产生强烈的降水。
图4 2024年3月27日700 hPa环流形势

(a)08:00, (b)12:00, 黑色三角表示长水机场

Fig.4 700 hPa circulation at 08:00 (a) and 12:00 (b) on 27 March 2024.The black triangles represent Kunming Changshui International Airport

地面天气图显示, 08:00[图5(a)]冷空气聚集在青藏高原一带, 冷高压前缘形成一条狭长的冷锋, 在云南与贵州的交界地带存在一个低压中心, 昆明长水机场正处于这个低压系统的控制之下。到13:00[图5(b)], 冷高压中心有所加强, 冷锋向昆明靠近, 由此判断大风天气发生时刻, 有冷锋逼近长水机场, 5 h内冷锋移动了约200 km, 移速约为40 km·h-1, 移动相对较快。
图5 2024年3月27日海平面气压场

(a)08:00, (b)13:00, 黑色三角表示长水机场

Fig.5 Sea level pressure at 08:00 (a) and 13:00 (b) on 27 March 2024.The black triangles represent Kunming Changshui International Airport

4 卫星云图上的重力波云系

4.1 两个重力波云系向机场移动

从3月27日09:00和10:00可见光卫星云图上[图6(a), (b)] 可以清晰地观察到, 滇西北横断山脉地区, 即昆明以西存在重力波云系A, 且云系呈现出向东南延伸的态势, 表明重力波正向东南方向传播。11:00和12:00[图6(c), (d)]在昆明西南方向, 即哀牢山脉地区出现另一片重力波云系B, 并向东北方向移动。两组重力波均向昆明传播, 并相互逼近。09:53两者在哀牢山区附近100.56°E, 25°N处相遇, 耦合为一个新的重力波C, 继续向昆明移动发展, 并对长水机场大风造成影响。
图6 2024年3月27日可见光云图

(a)09:00, (b)10:00, (c)11:00, (d)12:00, 三角形表示长水机场, 五角星表示两个重力波相遇的点, 箭头表示重力波及其传播方向

Fig.6 Visible image at 09:00(a), 10:00 (b), 11:00(c), 12:00 (d) on 27 March 2024.Triangle is Changshui International Airport, arrow is gravity wave propagation direction, pentagram is the point where two waves met

4.2 云系的云顶和云底高度

10:00云图上可见昆明上空晴朗无云, 亮温为21 ℃, 此为卫星所测得的昆明市区地面温度。根据加密中尺度地面气象观测资料, 昆明此时的实际地面温度为21.1 ℃, 与云图测量结果存在0.1 ℃误差。红外云图表明重力波A云系云顶亮温约为266 K, 对误差进行订正后, 由云顶实际亮温结合探空资料, 计算得出云顶高度为5100 m。根据ERA5资料得出重力波造成的云底高度为2600 m。同样可得出重力波B云系的云顶高度和云底高度, 发现两片云系的云顶和云底大致处于同一高度。
A区云系是西风遇横断山脉时强迫产生, 由于上午阳光照射, 造成山脉阳坡有更强的上升气流, 重力波也更强, 11:00较09:00重力波波长更短, 频率更高能量更强, 更加有利于昆明机场大风天气。

5 重力波动力学方程和相关参数分析

5.1 受力分析与重力波表达式分析

在平衡高度EL附近的一小段垂直距离上, 通过受力分析研究气块垂直方向的运动方程(图7)。因气块垂直方向所受浮力与气块和环境大气之间的温差成正比。由图7可知, 该温差与气块到EL的距离成正比, 该气块所受的力f可用以下公式表示:
f = m a = m h 2 t 2 = - k h
式中: k为常数; h为重力波中空气块偏离平衡位置的垂直位移。由于位移和受力方向相反, 所以式(1)中有负号。T ve表示环境空气虚温; T vp表示气块虚温。
图7 气块垂直受力分析示意

Fig.7 Schematic diagram of vertical force analysis on an air parcel

求解该二阶偏微分方程可知, 此方程为简谐振动即大气重力波, 该方程的解是正弦或余弦函数。
因空气的水平运动, 即风吹向山脉区域时受地形影响产生强迫上升(图8)。由初始条件可知重力波的波动为正弦波, 初相位为0, 且波动沿着风向下游传播, 其波动表达式为:
图8 重力波传播示意图

箭头为风向, 黑色三角为山脉, 竖长划线为长水机场位置

Fig.8 Schematic diagram of gravity wave propagation.The arrows represent wind direction, the black triangle represents the mountain range, the long dash line represents the location of Changshui International Airport

h = A s i n 2 π λ ( x - c t )
式中: A为振幅; λ为重力波波长; c为波速。

5.2 重力波相关参数分析

从卫星云图和探空资料上, 读取上述重力波云系的云顶和云底高度作为重力波的波峰和波谷高度, 即重力波波峰高度为5100 m, 波谷高度为2600 m, 可计算出两个重力波的振幅(5100-2600)/2=1250 m。而重力波实际振幅应该比云层展现的该数值更大, 因为在云顶之上一般仍会存在重力波的上升气流, 但若大气十分干燥, 其中的水汽已经全部凝结, 则更高处的重力波不会有成云的水汽, 所以无法通过云体反映重力波的上限。同样, 其下限存在水汽受热蒸发、 云体消散的情况, 此时尽管存在重力波的下沉运动, 但更低处其真实位置无法通过云体表现。因此本文讨论的重力波振幅, 是由云图表现的较小值, 即直接可测振幅, 真实的振动相对更强, 振幅更大。
根据10:00 -11:00云图得出重力波A向东南传播, 在1 h内向东传播了约64 km, 向南传播了约11 km。在笛卡尔坐标系中分别取向东和向北为正, 由此可计算出重力波速度Vax =18 m·s-1Vay =-3 m·s-1。测量云图上波状云之间的距离, 得出波长λa=12165 m, 得x方向上的重力波波长λax 约为12000 m, y方向上的重力波波长λay 约为2000 m。则x方向上重力波A的波动方程为:
h a x = 1250 s i n 2 π 12000 ( x - 18 t )
同理可得, y方向上的重力波A的波动方程为:
h a y = 1250 s i n 2 π 2000 ( - y - 3 t )
根据波动方程可计算出重力波周期:
T a x = λ a x V a x = 12   k m 18   m · s - 1 667 ( s ) 11.1 m i n
即在这一周期内, 会有一次重力波的强迫上升和下沉的振动过程。此外, 还可计算出重力波A的角速度ω和重力波频率F
ω a x = 2 π T = 2 π 667 s 0.009 r a d s - 1
F a x = 1 T = 1 667 s 1.5 × 10 - 3 H z
重力波B向东北偏东方向传播, Vbx =15 m·s-1Vby =5 m·s-1, 测量得出波长λb =15000 m, 将其分解为沿x方向(即向东方向)与y方向(即向北方向)的波动, 则x方向波长为14230 m, y方向波动为4743 m。得到其波动方程:
h b x = 1250 s i n 2 π 14230 ( x - 15 t )
h b y = 1250 s i n 2 π 4743 ( y - 5 t )
该重力波波动周期Tb =949 s≈15.8 min, 角速度为0.0067 rad·s-1, 频率为0.001 Hz。

5.3 谐波分析

重力波A和重力波B相遇耦合为一个新的波动, 在x方向上其表达式如下:
h x = 1250 s i n 2 π 12000 x - 18 t + s i n 2 π 14230 x - 15 t
写为波动形式, 则上式为:
h x = 2500 c o s 2 π 223 34152000 x - 7614 34152000 t · s i n 2 π 2623 34152000 x - 43614 34152000 t
该波动如图9所示, 其振幅本身是一个余弦波动, 振幅表达式为:
A ' = 2500 c o s 2 π 223 34152000 x - 7614 34152000 t
图9 重力波耦合后x方向波列图

x轴上零点为两个重力波相遇的点, y轴上零点为3850 m的高空, 竖线为机场位置, 239000 m为重力波耦合点到机场的东西方向距离, 自A点(三角, 波峰, 振幅2403 m)到B点(圆点, 机场, 振幅645 m)为影响机场大风的重力波做功区域

Fig.9 The wave train diagram in the x-direction of the coupled gravity waves.The zero point on the x-axis indicates the intersection of the two gravity waves, the zero point on the y-axis represents an altitude of 3850 m, the vertical line denotes the airport location, 239000 m is the east-west distance from the gravity wave coupling point to the airport, The zone from point A (triangle, wave crest, amplitude 2403 m) to point B (circle, airport, amplitude 645 m) is the work region of the gravity wave affecting the strong winds at the airport

该耦合波动的最大振幅为2500 m, 波动平均波长为11074 m, 波群速 C g = λ b V a - λ a V b λ b - λ a = 34 m s - 1
y方向上耦合后的波动方程为:
h y = 1250 s i n 2 π 2000 ( - y - 3 t ) + s i n 2 π 4743 ( y - 5 t )
写成波动形式上式为:
h y = 2500 c o s 2 π - 6743 18972000 y - 4299 18972000 t · s i n 2 π - 2743 18972000 y - 24229 18972000 t
此波动的波列如图10所示, 该耦合波动的最大振幅为2500 m, 平均波长为2001 m, 波群速Cg =1.5 m·s-1, 向北传播。
图10 耦合重力波y方向波列图

x轴上零点为两个重力波相遇的点, y轴上零点为3850 m的高空, 竖线为机场位置, 10500 m为重力波耦合点到机场南北方向的距离, 自A点(三角, 波谷, 振幅-39 m)到B点(圆点, 机场, 振幅-32 m)为影响机场大风的重力波做功区域

Fig.10 The wave train diagram in the y-direction of the coupled gravity waves.The zero point on the x-axis indicates the intersection of the two gravity waves, the zero point on the y-axis represents an altitude of 3850 meters, the vertical line denotes the airport location, 10500 m is the north-south distance from the gravity wave coupling point to the airport, the zone from point A (triangle, wave trough, amplitude -39 m) to point B (circle, airport, amplitude -32 m) is the work region of the gravity wave affecting the strong winds at the airport

6 重力波对大风的影响

6.1 重力波对大风的直接和间接影响

由卫星云图可知, 09:53两个重力波在100.56°E, 25°N相遇, 耦合为一个新的重力波向前传播, 重力波相遇的点在x方向上距离机场约239 km, y方向上距离机场10.5 km, 其能量以波动群速度传播。根据波动的群速度, 此耦合的重力波在11:50到达机场上空。
x方向上的重力波到达机场时, 波动振幅为645 m(图9中B点), 在到达机场前约74 s, 即11:48:46, 波动在机场西侧2510 m, 位于波峰处, 振幅为2403 m(图9中A点), 从A点到B点过程中, 有向下的负浮力做功, 产生向下的垂直运动。做功的表达式为:
W x * = 0 h f d h = - 0 h k h d h = - 1 2 k h 2
计算出系数k式(16)]和波峰到机场时下沉强迫做功Wx *式(17)]的数值:
k = - f h = - g T v e 2 - T v p 2 h T v e 2 = 6.2 × 10 - 5   ( N · m - 1 )
W x * = 1 2 × 6.2 × 10 - 5 × 2403 2 - 1 2 × 6.2 × 10 - 5 × 645 2 = 166.1   ( J · k g - 1 )
理想状态下, 如果这些能量全部转化为垂直运动, 可造成的最大下沉运动速度为W max
W m a x = 2 × W x * = 2 × 166.1 18   ( m · s - 1 )
x方向上耦合的重力波可贡献约18 m·s-1的下沉速度, 对大风形成起促进作用。
y方向上的重力波到达机场时, 波动振幅为-32 m(图10中B点), 在到达机场前23 s, 即11:49:37, 波动位于波谷附近, 振幅为-39 m(图10中A点), 从A点到B点有向上的正浮力做功产生上升运动, 计算表明可造成约1 m·s-1的垂直上升运动, 对大风形成有一定抑制作用。综合x方向和y方向, 耦合重力波的总体作用可产生约17 m·s-1的垂直下沉运动, 不考虑其他因素, 在地表刚性下垫面上转化为水平大风。
从11:00与12:00云图[图6(c), (d)]上可见在昆明上空一带存在槽前西南急流云系, 由雷达观测(图略)发现, 其径向速度为27.5 m·s-1。重力波可将中低空环境西南急流向地面传输, 造成高效的动量下传。

6.2 低空区域航站楼受损点的风力分析

受损建筑区域距地表71 m, 与地面风的10 m观测高度相比, 位置高出60 m以上。在近地层大气中, 由于地表复杂地形的影响, 摩擦作用显著, 且在风速较大时摩擦阻力作用更为强烈。所以在建筑受损区域的70 m以上位置, 实际风速应远大于地面10 m风杆测得的气象记录风速。长水国际机场的航站楼作为国家大型公共建筑, 其设计标准和施工水平均有极高的国家规范要求。在防风抗风方面, 该航站楼的设计构造应足以抵御超强台风的侵袭。但这次大风造成该机场屋顶受损, 表明当时在该航站楼顶部的风速应该超过了强台风(中心附近风速41.5~50.9 m·s-1)标准。
在理想状态下, 重力波直接贡献(风速17 m·s-1)与高空急流(风速27.5 m·s-1)的动量下传相结合, 不再考虑其他因素, 即可产生44.5 m·s-1的强大水平风速。受损区域处于机场主体建筑群正面的迎风点, 在如此强烈大风的冲击下必然受到极大的考验。
本次过程中, 机场记录到的最大风速仅为27 m·s-1, 远远小于本文计算所得风速。如果航站楼受损区域仅为上述实测风速, 则不可能造成破坏。由此表明低空对风速的分析预报应慎重参考地面风速记录, 更应综合考虑各种影响因素对大风的贡献。

6.3 重力波对低空大风贡献占比分析

6.3.1  DCAPE对大风的影响

下沉对流有效位能(Downward Convective Available Potential Energy, DCAPE)是空气块从自由下沉起始高度经干绝热过程下降至地面时, 气块因较环境大气温度更低、 密度更高而产生负浮力做功所获得的最大动能, 是一个衡量雷暴大风潜势的重要参数。
大风发生时刻, 机场DCAPE值为610 J·kg-1。不考虑其他因素, 假设下沉空气块起始高度处垂直速度为零, 其接触下垫面时冷性出流速度可通过以下公式计算:
E = D C A P E = 1 2 m W 2
式中: 能量E为DCAPE; m为下沉空气质量, 取单位质量则为1; W为垂直下沉运动速度。当下沉空气接地时, DCAPE被完全释放, 下沉运动速度达到最大:
W = 2 × D C A P E = 2 × 610 35 m · s - 1
即理想状态下, 若DCAPE完全释放可产生35 m·s-1下沉气流, 并在地表转化为水平风速。

6.3.2 重力波对大风贡献占比

耦合重力波和DCAPE协同作用, 可产生强烈的下沉气流, 理想状态在地表全部转化为水平风, 同时将27.5 m·s-1的高空急流引导至地面。与地面原有的7 m·s-1风速叠加时, 地表风速可激增至7+17+35+27.5=86.5 m·s-1。大风最高等级是17级, 对应风速为56.1~61.2 m·s-1。而这里计算出的数值是最高等级风速上限的142%, 这一速度极高, 必然会造成灾害。在此情况下分析各项占比可知, DCAPE对最终风速的贡献约为40.46%, 中低空急流贡献了31.79%, 而地面原有的水平风速占8.09%。耦合重力波则贡献了约19.65%的风速。

6.4 重力波耦合效应对低空大风的贡献

如果仅有重力波A或B并未形成两者耦合, 其振幅为1250 m, 则最大可造成10 m·s-1的下沉气流, 约可贡献12.58%的风速。这说明耦合重力波通过相互作用, 加强了下沉气流, 进一步增强了大风风速, 导致风速增加了7 m·s-1, 为大风额外贡献了8.8%的风速。

7 讨论

本次过程中耦合重力波的影响体现在, 其下沉波动一方面与DCAPE造成的下沉气流共同作用, 将中低空环境大气中的西南急流向地面传输, 造成高效的动量下传; 另一方面耦合重力波产生的下沉运动能量在低空转化为水平运动, 两个效应均有助于形成大风, 且耦合重力波可显著增强下沉气流促成灾害。
本文在计算中, 仅以卫星云图可视的云顶和大气垂直剖面资料为依据, 对重力波的振幅进行了相对保守的估计, 所以推算出重力波的贡献偏低, 实际情况重力波更强且对大风的影响更大。本文对DCAPE能量转化为地面水平风这一过程进行了理想化假设, 但因近地面摩擦力增大、 夹卷以及湍流等作用, 其转化率不可能达到100%(Cohen et al, 2007王秀明等, 2023), 实际上DCAPE的贡献应该更小, 则重力波的作用相对更大。耦合重力波也会因近地面的摩擦力而减小, 而摩擦力的特点是速度越大受到的摩擦阻力越大, 效果越明显, 速度为0 m·s-1时摩擦力消失。在DCAPE造成的下沉运动强于耦合重力波形成的下沉运动时, DCAPE产生的下沉运动受摩擦阻力作用急剧衰减, 而耦合重力波的下沉运动受摩擦阻力影响相对较小衰减有限, 实际上DCAPE贡献的风速更小, 则耦合重力波作用应较上文分析更大。
机场只有地面10 m的风速记录而航站楼受损处无实测资料, 所以这里真实风速仍然无法精准确定。若有垂直梯度观测资料, 特别是破损处的实际资料, 则可以进行更加细致的低空推算, 从而为深入认识耦合重力波, 加强低空气象保障提供更为可靠的支撑。

8 结论

本文通过对2024年3月27日昆明长水国际机场遭受的罕见瞬时大风事件进行分析, 揭示了两支重力波耦合后在此次大风灾害中的重要作用。
(1) 源于横断山区的大气重力波A和哀牢山区的重力波B向昆明传播。两者在长水机场上游100.56°E, 25°N附近相遇产生耦合效应, 促成机场低空灾害性大风。
(2) 根据稳定气层空气受力情况, 分析出重力波的数学表达式。通过红外云图的云顶亮温结合大气垂直结构, 分别得到重力波云系的云顶与云底高度及显著振幅2500 m。利用卫星云图确定波长分别为12165 m和15000 m。通过重力波云系的整体传播速度, 得到其周期分别为11.1 min和15.8 min。从而获得了重力波的全部参数。利用谐波分析方法得到了耦合重力波数学表达式。
(3) 11:48:46东西方向的耦合重力波在机场西侧2510 m开始下沉, 引发机场18 m·s-1的下沉运动促进大风形成。11:49:37南北方向的耦合重力波在机场南侧开始上升, 产生1 m·s-1的上升运动。耦合重力波净效果是产生17 m·s-1的下沉气流。
(4) 耦合重力波与DCAPE能量释放共同作用产生动量下传, 将27.5 m·s-1的西南急流带至低空, 并与地面7 m·s-1的水平风速结合, 形成低空灾害性大风, 破坏了航站楼主体建筑。重力波对总体风速产生了大约19.65%贡献。
Adams-Selin R D Johnson R H2013.Examination of gravity waves associated with the 13 March 2003 bow echo[J].Monthly Weather Review141(11): 3735-3756.DOI: 10.1175/mwr-d-12-00343.1 .

Butchart N2014.The Brewer-Dobson circulation[J].Reviews of Geophysics52(2): 157-184.DOI: 10.1002/ 2013RG000448 .

Clark T L Hauf T Kuettner J P1986.Convectively forced internal gravity waves: Results from two-dimensional numerical experiments[J].Quarterly Journal of the Royal Meteorological Society112(474): 899-925.DOI: 0.1002/qj.49711247402 .

Cohen A E Coniglio M C Corfidi S F, et al, 2007.Discrimination of mesoscale convective system environments using sounding observations[J].Weather and Forecasting22(5): 1045-1062.DOI: 10.1175/waf1040.1 .

Du Y Rotunno R2015.Thermally driven diurnally periodic wind signals off the east coast of China[J].Journal of the Atmospheric Sciences72(7): 2806-2821.DOI: 10.1175/jas-d-14-0339.1 .

Durran D R1990.Mountain waves and downslope winds[M].Atmospheric Processes over Complex Terrain.Boston, MA: American Meteorological Society: 59- 81.DOI: 10.1007/978-1-935704-25-6_4 .

Fang J Y Du Y2022.Aglobal survey of diurnal offshore propagation of rainfall[J].Nature Communications, 13: 7437.DOI: 10. 1038/s41467-022-34842-0 .

Fritts D C Alexander M J2003.Gravity wave dynamics and effects in the middle atmosphere[J].Reviews of Geophysics41(1): 1003.DOI: 10.1029/2001RG000106 .

Holton J R Haynes P H Mcintyre M E, et al, 1995.Stratosphere-troposphere exchange[J].Reviews of Geophysics33(4): 403-439.DOI: 10.1029/95RG02097 .

Houghton J T1978.The stratosphere and mesosphere[J].Quarterly Journal of the Royal Meteorological Society104(439): 1-29.DOI: 10.1256/smsqj.43901 .

Lane T P2015.GRAVITY WAVES| Convectively generated gravity waves[M].Encyclopedia of Atmospheric Sciences 2nd Edition, 3: 171-179.DOI: 10.1016/b978-0-12-382225-3.00489-8 .

Lane T P Clark T L2002.Gravity waves generated by the dry convective boundary layer: Two-dimensional scale selection and boundary-layer feedback[J].Quarterly Journal of the Royal Meteorological Society128(583): 1543-1570.

Limpasuvan V Richter J H Orsolini Y J, et al, 2012.The roles of planetary and gravity waves during a major stratospheric sudden warming as characterized in WACCM[J].Journal of Atmospheric and Solar-Terrestrial Physics78-79: 84-98.DOI: 10.1016/j.jastp.2011.03.004 .

Lindzen R S1981.Turbulence and stress owing to gravity wave and tidal breakdown[J].Journal of Geophysical Research Oceans86(C10): 9707-9714.DOI: 10.1029/JC086iC10p09707 .

Liu L Ran L K Gao S T2018.Analysis of the characteristics of inertia-gravity waves during an orographic precipitation event[J].Advances in Atmospheric Sciences35(5), 604-620.DOI: 10. 1007/s00376-017-7159-2 .

Malkus J S Stern M E1953.The flow of a stable atmosphere over a heated island, part 1[J].Journal of Meteorology10(1): 30-41.DOI: 10.1175/1520-0469(1953)010<0030: tfoasa>2.0.co; 2 .

Plougonven R Zhang F Q2014.Internal gravity waves from atmospheric jets and fronts[J].Reviews of Geophysics52(1): 33-76.DOI: 10.1002/2012RG000419 .

Richter J H Sassi F Garcia R R2010.Toward a physically based gravity wave source parameterization in a general circulation model[J].Journal of the Atmospheric Sciences67(1): 136-156.DOI: 10.1175 /2009 jas3112.1 .

Smith R B2010.The influence of mountains on the atmosphere[J].Advances in Geophysics, 21: 87-230.DOI: 10.1016/S0065-2687(08)60262-9 .

Tjernström M Mauritsen T2009.Mesoscale variability in the summer Arctic boundary layer[J].Boundary-Layer Meteorology130(3): 383-406.DOI: 10.1007/s10546-009-9354-x .

Wang Y Zhang L F Peng J, et al, 2018.Mesoscale gravity waves in the mei-yu front system[J].Journal of the Atmospheric Sciences75(2): 587-609.DOI: 10.1175/jas-d-17-0012.1 .

Yang H P Du Y Chen Z J, et al, 2024.Could developing frontal rainfall influence warm-sector rainfall [J].Geophysical Research Letters51(15): e2024GL110430.DOI: 10.1029/2024gl110430 .

Yang H P Du Y Wei J H2023.Generation of multiple gravity wave couplets from convection[J].Journal of the Atmospheric Sciences80(9): 2323-2343.DOI: 10.1175/jas-d-22-0212.1 .

鲍艳松, 季凌潇, 李欢, 等, 2024.面向空投的青藏高原风场大涡模拟研究[J].高原气象43(2): 293-302.DOI: 10.7522/j.issn.1000-0534.2023.00051.Bao Y S

Ji L X Li H, et al, 2024.Application of large eddy simulation on Qinghai-Xizang plateau wind for airdrop[J].Plateau Meteorology43(2): 293-302.DOI: 10.7522/j.issn.1000-0534.2023.00051 .

崔新艳, 陈明轩, 秦睿, 等, 2021.对流初生机理的研究进展[J].气象47(11): 1297-1318.DOI: 10.7519/ j.issn.1000-0526.2021.11.001.Cui X Y

Chen M X Qin R, et al, 2021.Research advances in the convective initiation mechanisms[J].Meteorological Monthly47(11): 1297-1318.DOI: 0.7519/j.issn.1000-0526.2021.11.001 .

董泽新, 吴硕岩, 叶芳, 等, 2025.基于WRF模式的CFD与LSTM技术对低空风切变数值模拟研究[J].高原气象44(2): 546-562.DOI: 10.7522/j.issn.1000-0534.2024.00119.Dong Z X

Wu S Y Ye F, et al, 2025.Numerical simulation of low-level wind shear using CFD and LSTM technology based on the WRF model[J].Plateau Meteorology44(2): 546-562.DOI: 10.7522/j.issn.1000-0534.2024.00119 .

付文卓, 陈斌, 徐祥德, 2024.青藏高原春季区域性极端大风频次下降成因[J].高原气象43(5): 1087-1101.DOI: 10.7522/j.issn.1000-0534.2024.00010.Fu W Z

Chen B Xu X D2024.On the causes of decreased regional extreme wind frequency over the Qinghai-Xizang Plateau during the spring season[J].Plateau Meteorology43(5): 1087-1101.DOI: 10.7522/j.issn.1000-0534.2024.00010 .

郭虎, 季崇萍, 张琳娜, 等, 2006.北京地区2004年7月10日局地暴雨过程中的波动分析[J].大气科学30(4): 703-711.DOI: 10.3878/j.issn.1006-9895.2006.04.15.Guo H

Ji C P Zhang L N, et al, 2006.A case study of local rainstorm in Beijing on 10 July 2004: the analysis of the gravity wave[J].Chinese Journal of Atmospheric Sciences30(4): 703-711.DOI: 10.3878/j.issn.1006-9895.2006.04.15 .

胡壮, 茅成铮, 姜楠, 等, 2022.我国低空飞行气象服务产品的需求及现状[J].气候变化研究快报11(5): 631-642.DOI: 0.12677/CCRL.2022.115066.Hu Z

Mao C Z Jiang N, et al, 2022.Demand and current situation of low-altitude flight meteorological service products in China[J].Climate Change Research Letters11(5): 631-642.DOI: 10.12677/ccrl.2022.115066 .

李晨轩, 韦志刚, 2024.1979-2021年中国南海海表风速和大风事件的变化特征[J].高原气象43(03): 696-710.DOI: 10.7522/j.issn.1000-0534.2023.00079.Li C X

Wei Z G2024.Variation characteristics of surface wind speed and gale events in the South China Sea from 1979 to 2021.[J].Plateau Meteorology43(3): 696-710.DOI: 10.7522/j.issn.1000-0534.2023.00079 .

李麦村, 1978.重力波对特大暴雨的触发作用[J].大气科学2(3): 201-209.

Li M C1978.Possible triggering of the gravity wave over torrential rain[J].Chinese Journal of Atmospheric Sciences2(3): 201-209.

李叶晴, 师春香, 沈润平, 等, 2023.基于WRF-LES的门头沟地区近地面风场模拟与边界层方案敏感性研究[J].高原气象42(3): 758-770.DOI: 10.7522/j.issn.1000-0534.2022.00084.Li Y Q

Shi C X Shen R P, et al, 2023.Simulation of near-surface wind over Mentougou with WRF-LES and sensitivity study of planetary boundary layer schemes[J].Plateau Meteorology42(3): 758-770.DOI: 10.7522/j.issn.1000-0534.2022.00084 .

林永辉, 廖清海, 王鹏云, 2003.低空急流形成发展的一种可能机制——重力波的惯性不稳定[J].气象学报61(3): 374-378.

Lin Y H Liao Q H Wang P Y2003.A possible developing mechanism of low level jet—the inertial instability of gravity wave[J].Acta Meteorologica Sinica61(3): 374-378.

卢冰, 史永强, 王光辉, 等, 2014.新疆克拉玛依强下坡风暴的机理研究[J].气象学报72(6): 1218-1230.DOI: 10.11676/qxxb2014.068.Lu B

Shi Y Q Wang G H, et al, 2014.Numerical study of severe downslope winds at Kelamayi, Xinjiang[J].Acta Meteorologica Sinica72(6): 1218-1230.DOI: 10.11676/qxxb 2014.068 .

陆汉城, 吕梅, 何齐强, 1992.一次冷锋后飑线的大振幅重力波特性分析[J].应用气象学报3(2): 136-144.

Lu H C M He Q Q1992.Characteristics of large amplitude gravity wave of squall line behind the cold front[J].Journal of Applied Meteorological Science3(2): 136-144.

马文倩, 李华宏, 陈小华, 等, 2025.云南不同类型大风时空特征分析[J].高原气象44(3): 780-794.DOI: 10.7522/j.issn.1000-0534.2024.00091.Ma W Q

Li H H Chen X H, et al, 2025.Analysis on the spatio-temporal characteristics for different types of strong winds over Yunnan[J].Plateau Meteorology44(3): 780-794.DOI: 10.7522/j.issn.1000-0534.2024.00091 .

孙继松, 王华, 2009.重力波对一次雹暴天气过程演变的影响[J].高原气象28(1): 165-172.

Sun J S Wang H2009.The impact of gravity wave on a hailstorm weather process evolution[J].Plateau Meteorology28(1): 165-172.

孙淑清, 1983.低空急流对内重力波不稳定发展的作用[J].大气科学7(2): 136-144.

Sun S Q1983.The effect of low level jet on the stability of gravity waves[J].Chinese Journal of Atmospheric Sciences7(2): 136-144.

孙艳辉, 李泽椿, 寿绍文, 2012.2007年3月3

-5 日辽宁省暴雪和大风天气的中尺度分析[J].气象学报70(5): 936-948.DOI: 0577-6619/2012/70(5)-0936-48.Sun Y H Li Z C Shou S W2012.A mesoscale analysis of the snowstorm event of 3-5 March 2007 in Liaoning Province[J].Acta Meteorologica Sinica, 70(5): 936-948.DOI: 0577-6619/2012/70(5)-0936-48 .

孙艳辉, 李泽椿, 寿绍文, 2015.一次暴风雪过程中的中尺度重力波特征及其影响[J].气象学报73(4): 697-710.DOI: 10.11676/qxxb2015.052.Sun Y H

Li Z C Shou S W2015.An investigation into the features and effects of the mesoscale gravity waves as in a snowstorm event[J].Acta Meteorologica Sinica73(4): 697-710.DOI: 10.11676/qxxb2015.052 .

覃睿, 2023.再论低空经济: 概念定义与构成解析[J].中国民航大学学报41(6): 59-64.

Qin R2023.Further discussion on low altitude economy: concept definition and composition analysis[J].Journal of Civil Aviation University of China41(6): 59-64.

覃卫坚, 寿绍文, 王咏青, 2013.大气对流层重力波研究进展[J].气象科技41(5): 864-869.DOI: 10.19517/j.1671-6345.2013.05.014.Qin W J

Shou S W Wang Y Q2013.Progresses in researches on tropospheric gravity waves[J].Meteorological Science and Technology41(5): 864-869.DOI: 10.19517/ j.1671-6345.2013.05.014 .

汤浩, 王旭, 储长江, 等, 2019.乌鲁木齐城区一次极端东南大风的形成机制——重力波与超低空急流耦合[J].干旱区地理42(6): 1229-1238.DOI: 10.12118/j.issn.1000-6060.2019.06.01.Tang H

Wang X Chu C J, et al, 2019.Formation mechanism of a southeast gale event in Urumqi urban area: The coupling of gravity wave and ultra-low level jet[J].Arid Land Geography42(6): 1229-1238.DOI: 10.12118/j.issn.1000-6060.2019.06.01 .

汪科均, 李国平, 谢家旭, 等, 2024.重力波在四川盆地东北部一次夜间山地暴雨中的作用初探[J].暴雨灾害43(2): 146-157.DOI: 10.12406/byzh.2023-115.Wang K J

Li G P Xie J X, et al, 2024.The role of gravity wave in a nocturnal mountain rainstorm in the northeastern Sichuan Basin[J].Torrential Rain and Disasters43(2): 146-157.DOI: 10.12406/byzh.2023-115 .

王文, 刘佳, 蔡晓军, 2011.重力波对青藏高原东侧一次暴雨过程的影响[J].大气科学学报34(6): 737-747.DOI: 10.13878/j.cnki.dqkxxb.2011.06.011.Wang W

Liu J Cai X J2011.Impact of mesoscale gravity waves on a heavy rainfall event in the east side of the Tibetan Plateau[J].Transactions of Atmospheric Sciences34(6): 737-747.DOI: 10.13878/j.cnki.dqkxxb.2011.06.011 .

王秀明, 俞小鼎, 费海燕, 等, 2023.下击暴流形成机理及监测预警研究进展[J].气象49(2): 129-145.DOI: 10.7519/j.issn.1000-0526.2022.100301.Wang X M

Yu X D Fei H Y, et al, 2023.A review of downburst genesis mechanism and warning[J].Meteorological Monthly49(2): 129-145.DOI: 10.7519/j.issn.1000-0526.2022.100301 .

魏栋, 田文寿, 陈泽宇, 等, 2016.青藏高原上空UTLS区域一次地形重力波过程中的物质上传[J].地球物理学报59(3): 791-802.DOI: 10.6038/cjg20160303.Wei D

Tian W S Chen Z Y, et al, 2016.Upward transport of air mass during a generation of orographic waves in the UTLS over the Tibetan Plateau[J].Chinese Journal of Geophysics59(3): 791-802.DOI: 10.6038/cjg20160303

吴迪, 王澄海, 何光碧, 2016.青藏高原地区夏季两次强降水过程中重力波特征分析[J].高原气象35(4): 854-864.DOI: 10.7522/j.issn.1000-0534.2015.00066.Wu D

Wang C H He G B2016.Gravity wave characteristics in two summer heavy rainfall in the Qinghai-Xizang Plateau[J].Plateau Meteorology35(4): 854-864.DOI: 10.7522/j.issn.1000-0534.2015.00066 .

谢家旭, 李国平, 2021.重力波与对流耦合作用在一次山地突发性暴雨触发中的机理分析[J].大气科学45(3): 617-632.DOI: 10.3878/j.issn.1006-9895.2009.20137.Xie J X

Li G P2021.Mechanism analysis of a sudden rainstorm triggered by the coupling of gravity wave and convection in mountainous area[J].Chinese Journal of Atmospheric Sciences45(3): 617-632.DOI: 10.3878/j.issn.1006-9895.2009.20137 .

徐景峰, 宋林烨, 陈婧, 等, 2025.北京冬奥会延庆复杂地形冬季和早春地面风场精细特征对比研究[J].高原气象44(2): 378-392.DOI: 10.7522/j.issn.1000-0534.2024.00071.Xu J F

Song L Y Chen J, et al, 2025.Comparative study on fine characteristics of surface wind field in winter and early spring over Yanqing complex terrain during Beijing Winter Olympics[J].Plateau Meteorology44(2): 378-392.DOI: 10.7522/ j.issn.1000-0534.2024.00071 .

杨坤琳, 高志伟, 高兵, 等, 2025.横断山区东南缘洱海东岸风场及风切变时空特征研究[J].高原气象44(2): 462-474.DOI: 10.7522/j.issn.1000-0534.2024.00078.Yang K L

Gao Z W Gao B, et al, 2025.Spatiotemporal distribution characteristics research of wind field and wind shear on the east coast of Erhai Lake in the southeast edge of Hengduan Mountain area[J].Plateau Meteorology44(2): 462-474.DOI: 10.7522/j.issn.1000-0534.2024.00078 .

杨显玉, 朱俊橙, 文军, 等, 2023.南疆大风气候特征分析及其对沙尘天气的影响[J].高原气象42(1): 186-196.DOI: 10.7522/j.issn.1000-0534.2021.00043.Yang X Y

Zhu J C Wen J, et al, 2023.Analysis on Characteristics of Gale Climate in South Xinjiang and Its Influence on Sandstrom[J].Plateau Meteorology42(1): 186-196.DOI: 10.7522/j.issn.1000-0534.2021.00043 .

张珊, 王宗敏, 黄刚, 等, 2023.基于WRF-LES的崇礼复杂地形局地风场模拟研究[J].高原气象42(1): 197-209.DOI: 10.7522/j.issn.1000-0534.2022.00011.Zhang S

Wang Z M Huang G, et al, 2023.Local wind simulation over complex terrain of Chongli using WRF-LES[J].Plateau Meteorology42(1): 197-209.DOI: 10.7522/j.issn.1000-0534.2022.00011 .

张云, 熊建刚, 万卫星, 2011.中层大气重力波的全球分布特征[J].地球物理学报54(7): 1711-1717.DOI: 10.3969/j.issn.0001-5733.2011.07.003.Zhang Y

Xiong J G Wan W X2011.Analysis on the global morphology of middle atmospheric gravity waves[J].Chinese Journal of Geophysics54(7): 1711-1717.DOI: 10.3969/j.issn.0001-5733.2011.07.003 .

Outlines

/