[論文レビュー] A Generalized Richardson Number Diagnostic for Turbulence in the Free Atmosphere

Mohamed Foudad, Miguel A. C. Teixeira+2 more|arXiv (Cornell University)|Feb 25, 2026

Meteorological Phenomena and Simulations被引用数 0

ひとこと要約

本論文は水平風切れを含む一般化リチャードソン数Ri_newを提案し、ERA5データと247M件の航空機乱流報告を用いてRi_oldおよびTI1より優れた予測性能を示す。

ABSTRACT

A new Richardson number formulation, Ri_new, is introduced to improve the diagnosis of turbulence in the stratified free atmosphere, particularly near jet stream regions. The formulation is derived from the turbulent kinetic energy budget and accounts for both vertical wind shear and horizontal shear (deformation and divergence), weighted by the ratio of horizontal to vertical eddy viscosities (Kmh/Kmv). This extends the classical Richardson number Ri_old, which includes only vertical shear, and provides a physically based measure of the balance between stratification and three-dimensional shear production. The diagnostics Ri_new, Ri_old, and the widely used Turbulence Index 1 (TI1), computed from ERA5 reanalysis, are evaluated using more than 247 million automated turbulence reports from commercial aircraft (2017-2024). Across various turbulence intensity thresholds, Ri_new consistently outperforms the other diagnostics, resulting in higher AUC values and improved probability of detection at operationally relevant false-alarm rates. Sensitivity analyses show that the predictive skill of Ri_new is maximized for Kmh/Kmv values in the range 10^3-10^4, with peak performance near 5000 and weak dependence on horizontal resolution. Seasonal and regional evaluations indicate that the added value of Ri_new is largest where turbulence generation involves both vertical and horizontal shear, such as over the contiguous United States and during summer. Over oceans, performance remains high and Ri_new still provides the best overall discrimination skill. These results demonstrate that incorporating horizontal wind shear into the Richardson number yields a physically consistent and statistically robust improvement in turbulence diagnostics, with relevance for research and operational applications.

研究の動機と目的

• 成層圏自由大気、特にジェット気流付近での乱流診断を改善する動機付け。
• 水平変形と発散を含む全TKE予算からRi_newを導出。
• ERA5再分析と大規模MADIS ACARS乱流データを用いてRi_newをRi_oldおよびTI1と比較評価。
• Kh/Kv比が予測 skillへ与える影響を特定。
• 季節・地域・乱流強度の多様性に対する適用性を評価。

提案手法

• 水平勾配を保持し、一次的乱流閉じ込めでRi_newを導出することで、全TKE予算からRi_newを得る。
• Ri_newをRi_g/Pr_tとして表現し、Ri_gにはN^2と拡張された水平変形項を取り入れる。
• Ri_g = N^2 / [ Sv^2 + (Kh/ Kv)(Div^2 + D_ST^2 + D_SH^2) ] を定義。Svは鉛直せん断、Deformation項は水平せん断を表す。
• Kh/Kv = 0を設定してRi_newをRi_oldに結びつけ、古典的Ri_oldを回復。
• ACARS観測に合わせてNES均一補間したERA5場を用いてRi_new、Ri_old、TI1を比較。
• 247百万件のMADIS乱流報告（2017–2024）とROC/AUCを主要検証指標として使用。

Figure 1: Spatial distribution of ACARS turbulence reports from the NOAA MADIS archive for the period 1 January 2017 to 31 December 2024. Shading indicates the number of reports on a logarithmic scale, binned onto a 1° × 1° latitude–longitude grid. Only reports above 8 km flight altitude are include

実験結果

リサーチクエスチョン

• RQ1 Ri_newは乱流強度の異なる領域で古典的Ri_oldおよびTI1と比べて乱流識別性を改善するか。
• RQ2 予測技能を最大化するKh/Kv比はどれで、Ri_newは水平解像度にどの程度敏感か。
• RQ3 三次元的（垂直+水平）せん断が乱流発生を支配する地域・季節でRi_newはどう機能するか。
• RQ4 オーシャンと陸上での上部対流圏/下部成層圏乱流に対してRi_newは頑健か。
• RQ5 Ri_newベースの診断は大規模なグローバルデータセットのACARS乱流報告とどれだけ良く整合するか。

主な発見

• Ri_newは複数の乱流強度閾値でRi_oldおよびTI1を一貫して上回り、AUC値が高い。
• Ri_newの最適なKh/Kv比は約5×10^3で最良のROC性能を示し、10^3–10^4の範囲でもスキルが改善。
• Ri_newは解像度に頑健で、垂直および水平せん断の両方が乱流発生を駆動する領域（例：米国本土周辺・夏季）で最大の付加価値を示す。
• オーシャン領域でも性能は高く、Ri_newはそれらの地域で最良の全体識別性を提供。
• 本手法は物理的に一貫し統計的にも頑健な乱流診断の改善を提供し、航空予測および気候影響評価に潜在的な利益をもたらす。

Figure 2: ROC curves for TI1 index (blue), the classical Richardson number (Ri old , black), and the optimal Ri new formulation (red) for turbulence intensity thresholds of (a) EDR $=$ 0.10, (b) EDR $=$ 0.20, and (c) EDR $=$ 0.30 $\mathrm{m}^{2/3}\,\mathrm{s}^{-1}$ . Gray curves show Ri new for diff