飛機雲預測

Contrails API 使用兩種不同的模型預測飛機雲:一種是以機器學習 (ML) 為基礎的模型,可預測可能出現飛機雲的區域;另一種是以物理為基礎的模型,稱為飛機雲卷雲預測 (CoCiP),可預測暖化。

機器學習模型

飛機雲可能區域 (CLZ) 預測模型採用機器學習技術,可預測飛機雲形成的機率。這個模型是深層類神經網路,會以天氣特徵做為輸入內容,並根據衛星偵測到的飛機雲預測 CLZ (Geraedts et al. 2023)。

輸入內容主要包含 HRES 天氣特徵:比濕、溫度、風的 u 分量、風的 v 分量、垂直速度、相對渦度、雲量比例、比雲冰水含量、比雪水含量和散度。相對濕度是根據比濕和溫度計算得出。模型的輸入內容也包括當地太陽時、一年中的日期,以及飛行路線控點的高度。

在某些地理區域 (例如美國),模型經過微調,會使用經度、緯度和海拔高度做為輸入特徵,而模型的全球備用方案則只會使用海拔高度 (不使用經度或緯度)。

根據觀測到的飛機雲資料評估,該模型可達到頂尖效能。

CoCiP 模型

飛機雲卷雲預測模型會預測飛機雲的能量強迫,這是衡量飛機雲對氣候衝擊的指標。

能量強迫的定義如下:

\[ EF [J] = \int_{0}^{t} RF'(t) \times L(t) \times W(t)dt \]

也就是飛機雲在生命週期內的瞬時輻射強迫積分 (Teoh et al. 2020)。能量強迫值會依飛行距離正規化,單位為 (J/m)。

CoCiP 是以物理學為基礎的模型,可利用大氣條件、飛機類型、飛行路徑和其他特徵,模擬飛機雲的形成、演變和影響 (Schumann 2012;Schumann et al. 2012)。模型會使用 ECMWF 高解析度預測集合 (HRES ENS) 的 10 個集合成員做為輸入內容,用於平流傳輸已形成飛機雲的飛行路線控點 (Hersbach et al. 2020)。CoCiP 模型也會使用雲微物理理論,判斷哪些雲凝結尾會持續存在,並考量初始下降氣流、下降和昇華。CoCiP 會根據模擬的飛機雲演變,計算飛機雲屬性和周圍天氣條件的能量強迫。

除了 CoCiP 的能量強迫估算值外,系統也會計算能量強迫的氣候估算值,方法是將一年的 CoCiP 輸出結果平均,並依時段、季節和緯度分組。

最終的能量強迫量是 CoCiP 集合成員的能量強迫平均值 (EF 非零),以及氣候平均值 (一律非零)。將氣候學納入平均值,可確保即使 CoCiP 未使用任何天氣集合成員預測飛機雲的形成,也能估算飛機雲的影響。

預期有效能量強迫

接著,系統會將機器學習模型計算出的飛機雲形成機率,乘以 CoCiP 模型計算出的飛機雲有效能量強迫,得出預期有效能量強迫。

飛機雲強迫指數值

接著,系統會將計算出的預期有效能量輻射強迫值對應至 0 到 4 的量表,產生凝結尾輻射強迫指數值。凝結尾強迫指數值是根據亂流預測值而定。

Contrails API 會透過裁剪和線性比例,將 expected_effective_energy_forcing 轉換為 contrails 嚴重程度指數:

ds["contrails"] = ds["expected_effective_energy_forcing"].clip(min=2e7, max=2e8)
ds["contrails"] = ((ds["contrails"] - 2e7) / (2e8 - 2e7)) * 4

如要進一步瞭解這項對應關係的背景資訊,請參閱「能量強迫解讀」。

此外,您也可以參閱 Contrails API 版本資訊,並訂閱公告,以便接收這項量表的異動通知。

二氧化碳當量

預期有效能量強迫也可解讀為飛機雲對氣候的影響,與二氧化碳相當,即二氧化碳當量。

如要計算預測飛機雲暖化的 $\text{CO}_2$ 等效值,請使用下列轉換公式:

\[ \text{CO}_2\text{eq [kg]} = EEEF \text{ [J]} \times \text{GWP Factor} \]

其中:

  • EEEF 是指預期有效能量強迫
  • GWP 係數是全球暖化潛勢係數,也就是 1 公斤 $\text{CO}_2$ 在特定時間範圍 (H,以年為單位) 內產生的總整合能量強迫的倒數。可以是 GWP 20、50 或 100。

GWP 係數產品

根據 20、50 和 100 年的時間範圍,提供的 GWP 係數如下:

時間範圍 GWP 係數 (kg-$\text{CO}_2$/J)
GWP20 2.497e-9
GWP50 1.173e-9
GWP100 6.779e-10

GWP 係數計算

這些 GWP 係數是根據 $\text{CO}_2$ 的絕對全球暖化潛勢 (AGWP),使用 IPCC AR5 係數和下列計算方式得出:

\[ \text{GWP Factor} = \frac{1}{AGWP_{\text{CO}_2, H} \times A_{\text{earth}} \times S_{\text{year}}} \]

其中:

  • $AGWP_{\text{CO}_2, H}$:所選時間範圍內 $\text{CO}_2$ 的絕對全球暖化潛勢,資料來源為 IPCC AR5
  • $A_{\text{earth}}$:地球表面積 ($5.101 \times 10^{14} \text{ m}^2$)
  • $S_{\text{year}}$:一年的秒數 ($31,536,000 \text{ s}$)

授權

Forecast API 提供的資料依據 CC BY 4.0 授權。

後續步驟

參考資料

  • Geraedts, Scott, Erica Brand, Thomas R. Dean、Sebastian Eastham、Carl Elkin、Zebediah Engberg、Ulrike Hager 等人,2023 年。"A Scalable System to Measure Contrail Formation on a per-Flight Basis." arXiv [physics.ao-Ph]. arXiv. http://arxiv.org/abs/2308.02707.
  • Hersbach, Hans, Bill Bell, Paul Berrisford, Shoji Hirahara, András Horányi, Joaquín Muñoz-Sabater, Julien Nicolas, et al. 2020. "The ERA5 Global Reanalysis." Quarterly Journal of the Royal Meteorological Society 146 (730): 1999-2049.
  • Schumann, U. 2012 年。「飛機雲卷雲預測模型。」Geoscientific Model Development 5 (3): 543-80.
  • Schumann, U.、B. Mayer, K. Graf 與 H. Mannstein。2012 年。「A Parametric Radiative Forcing Model for Contrail Cirrus」(凝結尾卷雲的參數化輻射強迫模型)。Journal of Applied Meteorology and Climatology 51 (7): 1391-1406.
  • Shapiro, Marc、Zeb Engberg、Roger Teoh、Marc Stettler 和 Tom Dean。2023. Pycontrails:用於模擬航空氣候衝擊的 Python 程式庫。 https://doi.org/10.5281/zenodo.825291
  • Teoh, Roger, Ulrich Schumann, Arnab Majumdar, and Marc E. J. Stettler。2020. 〈Mitigating the Climate Forcing of Aircraft Contrails by Small-Scale Diversions and Technology Adoption〉。Environmental Science & Technology 54 (5): 2941-50.
  • Teoh, Roger, Zebediah Engberg, Ulrich Schumann, Christiane Voigt, Marc Shapiro, Susanne Rohs, and Marc E. J. Stettler。2024。 「2019 年至 2021 年全球航空飛機雲對氣候的影響」。 Atmospheric Chemistry and Physics 24: 6071–6093. https://www.researchgate.net/publication/384141259_Global_aviation_contrail_climate_effects_from_2019_to_2021.