Hava rotalamasında tahmin belirsizliği: ensemble yayılımı ve karar penceresi

Aday rotalarda tahmin üyesi kimliğini koruyun, ilgili olayı sayın; sonraki güncellemenin anlamlı karar sınırından önce yetişip yetişmeyeceğini belirleyin.

Bu sayfada

Hava rotalaması kararı, yalnız tek noktadaki tahmine veya en düzgün ortalama alana değil, hareketli geminin yolu boyunca olası koşullara bağlıdır. Özgün on üyeli örnek ensemble yayılımını, rota boyunca eşik taramasını ve yeni tahmini kullanmak için kalan süreyi ayırır. Hesaplar kanıtı karşılaştırır; gerçek rota seçmez.

Ensemble okumadan önce kararı tanımlayın

Aday rotaları, beklenen geçiş zamanlarını ve kararla ilgili gemiye özgü büyüklükleri belirtin. Anlamlı dalga yüksekliği, dalga periyodu, yön, rüzgâr ve gemi tepkisi farklı değişkenlerdir. ECMWF dalga diyagramı rehberi anlamlı dalga yüksekliğini en büyük tekil dalga yüksekliğinden ayırır. Birine bağlı eşik diğerine sessizce taşınmamalıdır. Değişken tanımı, olasılık yorumundan önce gelir.

Buradaki eğitim olayı, belirtilen iki rota noktasından en az birinde Hs > 4,0 m olmasıdır. Eşik aritmetik için uydurulmuştur; işletme sınırı değildir. Kontrol noktaları her rotanın seçilmiş konum–zaman örnekleridir; aradaki her noktayı kapsamaz. Bu nedenle doğru sayım bile bütün hava maruziyetini veya gemi hasarını değil, tanımlanmış olayı tarar. Örneklemenin kapsamı sonuçta korunmalıdır.

Yayılımın neyi temsil ettiğini anlayın

ECMWF ensemble yaklaşımını olası tahmin gelişmelerindeki belirsizliği temsil etme yolu olarak açıklar. Üyeler bağımsız gözlenmiş gelecekler değil, belirsiz başlangıç koşulları ve modellemeye ilişkin varsayımları yansıtır. Yayılım, modelin temsil ettiği seçenekler hakkında bilgi verir; bütün tahmin hatası kaynaklarını veya yerel değişkenliği kendiliğinden kapsamaz. Üyelerin benzerliği tam bilgiyle eş tutulmamalıdır.

Sonlu eğitim kümesinde ortalama μ = ΣHi/N, yayılım s = √[Σ(Hi − μ)²/N] ve H’nin birimi metredir. Payda N’dir; dolayısıyla bu, listelenen üyelerin topluluk yayılımıdır. Küçük yayılım üyelerin birbirine yakın olduğunu söyler; gerçekle uyuştuklarını kanıtlamaz. Doğrulama ayrı kanıt gerektirir. Yayılımın hangi örnek ve değişken üzerinde hesaplandığı açık tutulmalıdır.

Tam eşleştirilmiş senaryo kümesi kullanın

Tablo eşit ağırlıklı on kurgusal üyeyi verir. A1 ve A2, A rotasının iki kontrol noktasıdır; B1 ve B2, B rotasının kendi varsayılan zamanlamasındaki karşılık gelen iki noktadır. Her satır dört değer boyunca aynı üye kimliğini korur. Bunlar indirilmiş hava verisi veya gerçek deniz alanı iddiası değil, analiz için üretilmiş değerlerdir. Satır eşleşmesi hesabın temel girdisidir.

Eşleştirme, ortalama almadan önce “bu üye rotanın örneklenen herhangi bir yerinde eşiği aşıyor mu?” sorusunu yanıtlatır. Her sütunu ayrı sıralamak veya her noktada ayrı uygun üye seçmek, hiçbir üyenin öngörmediği yapay yolculuk yaratır. İki seçenek de kendi beklenen geçiş zamanında değerlendirilmelidir; farklı rotalar aynı hava sistemiyle farklı evrelerde karşılaşabilir. Zamanlama, mekânsal seçimin parçasıdır.

ÜyeA1A2B1B2
12,83,13,13,2
23,03,43,23,3
33,23,73,33,4
43,44,13,43,5
53,64,43,53,6
63,84,63,63,7
74,23,83,73,8
84,43,53,83,9
93,53,33,94,0
103,73,64,14,2

Aynı ortalamayı farklı yayılımla karşılaştırın

İlk kontrol noktasında iki rotanın ortalaması 3,56 m’dir. Altı ondalığa yuvarlandığında A yayılımı 0,473709 m, B yayılımı 0,303974 m olur. İlk nokta aşılma oranları 0,2 ve 0,1’dir. Dolayısıyla aynı ortalama, özellikle küçük seçilmiş bu senaryo kümesinde bile eşiğe göre farklı dağılımı gizler. Merkez değeri tek başına kuyruk davranışını açıklamaz.

İkinci noktada ortalamalar A için 3,75 m, B için 3,66 m’dir. Dört nokta ortalaması da 4,0 m’nin altındayken tekil üyeler eşiği aşar. Yalnız ortalamayı taramak bu üyeleri kaçırır. Ortalama sonrasında eşik uygulamak, üye düzeyindeki eşik sonuçlarını ortalamaktan farklı işlemdir; iki işlem genellikle yer değiştiremez. Hesap sırası tanımlanan olayı etkiler.

Rota olayını her üyenin içinde sayın

max(Hm,r,1, Hm,r,2) > 4,0 m ise Im,r = 1, değilse sıfır tanımlayın. Rota oranı pr = ΣIm,r/10 olur. A’da 4, 5 ve 6 numaralı üyeler ikinci noktada; 7 ve 8 birinci noktada aşar. Böylece beş ayrı üye olayı tetikler; örneklenen rota için 0,5 elde edilir. Üye kimliği bu toplamın doğru kurulmasını sağlar.

B’de yalnız 10 numaralı üye iki noktada da aşar; aynı üye iki kez sayılmadan 0,1 bulunur. 9 numaralı üyenin 4,0 m değeri katı büyüktür eşiğini aşmaz. Eşitsizliği belirtmek olay tanımının parçasıdır; büyük veya eşite geçmek sonucu değiştirir. Yerel şekil toplam oranlardan önce üye düzeyindeki mantığı gösterir. Sayımın hangi satırlardan geldiği görünürdür.

On eşli kurgusal üye A ve B rotalarının iki noktasında eşik göstergeleriyle sunulur. A’da 4–8 üyeleri en az bir noktada 4 metre anlamlı dalga yüksekliğini aşar; oran 0,5’tir. B’de yalnız 10 numaralı üye aşar; oran 0,1’dir. 14:00 ayrımı, 45 dakika eylem hazırlığı ve 30 dakika inceleme için son yararlı veri zamanı 12:45 UTC’dir.
Her rota için yalnız iki noktada örneklenmiş özgün kurgusal eşit ağırlıklı üyeler. Katı Hs > 4,0 m olayı eğitim eşiğidir; gemi sınırı veya en büyük tekil dalga yüksekliği değildir. Oranlar kalibre edilmemiş senaryo sayımlarıdır. Ayrı zaman kutusu varsayılan inceleme ve eylem sürelerini kullanır.

Dayanaksız bağımsızlık kısayolundan kaçının

ECMWF olasılık rehberi birleşik zaman veya alan olaylarının neden üye düzeyindeki bağımlılığı gerektirdiğini açıklar. A’da nokta oranları 0,2 ve 0,3, bağımsızlık altında 1 − (1 − 0,2)(1 − 0,3) = 0,44 verir; gerçek eşli sayım 0,5’tir. Aşımlar ayrık üyelerde olur; bağımsızlık kısayolu başka ortak modeli anlatır. Aynı marjinaller birleşik olayı tek başına belirlemez.

B’de iki nokta oranı da 0,1’dir; kısayol 0,19 verirken bütünüyle aynı üyedeki olay 0,1 verir. Zıt yönlü hatalar bağımsızlığın genel muhafazakâr ikame olmadığını gösterir. Zaman ve mekân boyunca korelasyon tahmin yörüngesinin parçasıdır; nokta tahminini rota olayına çevirirken korunması gerekir. Bağımlılığın atılması sonucu yalnız büyütmez veya yalnız küçültmez.

Kalibrasyonu ve sonlu üye sınırını denetleyin

ECMWF güvenilirlik diyagramları tahmin olasılıklarını çok sayıda durumdaki gözlenmiş olay sıklıklarıyla karşılaştırır. Karşılaştırma, ilgili ortamda olasılığın görünen değeriyle kullanılıp kullanılamayacağını sınar. Dar ve güven veren ensemble sistematik yanlı olabilir. Kalibrasyon ayrıca ilgisiz toplu puan yerine temsil edici değişken, bölge, tahmin ufku ve olay tanımı gerektirir. Uygunluk bağlamla birlikte doğrulanmalıdır.

On üyeli alıştırmada bir üyenin olay durumu değişince oran 0,1 basamakla değişir. Bu kabalık belirsizlik aralığı değildir; sıfır aşım gerçek dünyada sıfır şans kanıtlamaz. Listelenen oranlar kalibrasyon verisi eklenmemiş açık senaryo özetleridir. Ek basamaklar bunları doğrulanmış gemi hasarı veya kabul edilebilir geçiş olasılığına dönüştüremez. Sayısal kesinlik ile tahmin doğruluğu ayrılmalıdır.

Rota karşılaştırmasını çok boyutlu tutun

Seçilen taramada B’de daha az aşan üye vardır. Ayrı olarak sabit 12 kn hızda 24 deniz mili eklediğini varsayın; nominal ek süre 2,0 h olur. Aritmetik yakıt tüketimini, akıntıları, hız kaybını, varış kısıtlarını veya gemi tepkisini içermez. B’yi en iyi ilan etmek yerine ek girdilere ihtiyaç duyan karşılaştırmayı açığa çıkarır. Tek ölçüt bütün amacı temsil etmez.

Gerçek rotalama modelinde hız ve varış zamanı tahmin üyesine ve seçilen kontrol politikasına bağlı olabilir. O zaman her yörünge kendi oluşan zamanlarında örneklenmelidir. Gemi yavaşladığında veya fırtınanın evresi kaydığında sabit zaman karşılaştırması değişebilir. Eğitim tablosu, bu ek geri beslemeyi sonuçlara gizlice katmamak için zamanlamaları bilerek sabit tutar. Modelde bulunmayan etki hesaplanmış gibi sunulmaz.

Son yararlı güncelleme zamanını hesaplayın

Kurgusal aynı günde son anlamlı rota ayrımına 14:00 UTC’de ulaşıldığını varsayın. Eylemin hazırlık süresine 45 min, yeni veriyi değerlendirmeye 30 min ayırın. Son yararlı veri hazır zamanı 14:00 − 45 min − 30 min = 12:45 olur. Sınır meteorolojik tahmin geçerlilik zamanı değil, karar varsayımıdır. Kararın uygulanma süresi hesaba açıkça katılır.

Kullanılabilir hâle 13:00’te gelen yeni veri, nominal model çevrimi çok daha önce başlamış olsa bile pencereyi 15 min kaçırır. Model başlangıcı kadar indirme, işleme ve inceleme süresi de önemlidir. Güncellemeyi beklemek ancak yorumlandıktan sonra karar uygulanabilir kaldığında yararlıdır; sonraki tahmin kaybedilmiş seçeneği kendi başına geri getiremez. Beklemenin değeri kalan seçeneklerle birlikte değerlendirilir.

Koşullu karar kaydını koruyun

Tahmin kaynağını ve çevrimini, değişken tanımlarını, rota zaman varsayımlarını, olay eşiğini, üyelerin ele alınışını ve kalibrasyon sınırlarını kaydedin. Hangi gözlemin veya yeni tahminin karşılaştırmayı değiştireceğini ve bilginin ne zamana kadar kullanılabileceğini belirtin. Böylece her yeni hava haritasını planı tersine çevirmek için bağımsız neden saymak yerine sınırları belirli yeniden değerlendirme oluşur.

Özgün örnek iki noktalı olayında 0,5 ve 0,1 rota oranlarını, aynı ilk nokta ortalamasıyla farklı yayılımları ve 12:45 UTC son veri hazır zamanını verir. Bunlar verilen alıştırmanın yeniden üretilebilir özellikleridir. Gerçek rota kararları güncel yetkili hava bilgisini, belirli geminin sınır ve tepkisini ve alıştırmanın bilerek modellemediği daha geniş seyir kısıtlarını gerektirir. Sonuç kendi olay tanımının sınırlarında kalır.

Kaynaklar

  1. ECMWF — Quantifying forecast uncertainty.
  2. ECMWF Forecast User Guide — Waves and wavegrams.
  3. ECMWF Forecast User Guide — Consideration of Probability.
  4. ECMWF Forecast User Guide — Statistical Concepts: Probabilistic Data.