Modele cyrkulacji zwykle potrafią wiarygodnie przepowiadać pogodę na co najwyżej 5 do 6 dni naprzód. Dlatego prognozy sezonowe w większości są nie lepsze niż wyniki rzucania monetą. Są jednak wyjątki, z czego najbardziej znanym jest przepowiadanie znaku zimowego indeksu NAO z wyprzedzeniem wielu miesięcy. Sam to robiłem z powodzeniem wiele razy zanim kosa trafiła na kamień tej zimy, stosując metody z opublikowanych prac naukowych (ostatnim razem zabrakło mi nieco wiary, że możemy mieć kolejną zimę z dodatnim NAO). A ponieważ indeks NAO bardzo silnie wpływa na to jaką będziemy mieli zimę, nie jest to coś be znaczenia praktycznego. Szczególnie gdy właśnie ukazał się artykuł, twierdzący, że da się zimowe NAO przewidzieć… na wiele lat naprzód.
Met Office już kilka lat temu twierdził w prasie naukowej [1], że jest w stanie zimowe NAO (czyli zasadniczo i temperaturę) przewidywać na rok na przód, ale… tego nie robił. Coś jak “A dziad wiedział, nie powiedział, a to było tak”. Ale coś wiemy o tym “tak” z tego artykułu. A mianowicie wiadomo jakie parametry używali w tej czysto statystycznej prognozie (i jak się one korelują z zimowym NAO). A są to:
- ENSO (indeks cyklu El Nino/La Nina),
- siła “wiru polarnego” (ale w listopadzie co mi niezbyt pasuje do rocznej prognozy),
- atlantycki układ temperatur powierzchni morza (słynna “trójbiegunowka” Met Office’u),
- lód morski na Morzu Kara.
Praca ta była kulminacją ich wieloletnich wysiłków, dających stopniowo coraz lepsze prognozy zimowego NAO, najpierw za pomocą modeli statystycznych, czyli opartych o obserwowane korelacje między zimowym NAO, a różnymi zjawiskami (jak wyżej) wiele miesięcy naprzód, a później także modelami cyrkulacji asymilującymi wartości tych parametrów. W pracy [1] twierdzili już, że są w stanie przewidzieć NAO na rok naprzód.
Teraz poszli jeszcze dalej i twierdzą, że da się znaku zimowego NAO, a zatem temperaturę i opady np. w Europie przewidzieć w skali dekadowej przy pomocy wiązki modeli. Tak, do 10 lat naprzód. Jest o tym nowy artykuł datowany dzisiaj [2]. Większość autorów obu prac jest ta sama.
W pracy użyli wiązki 169 realizacji 9 modeli z projektu Coupled Model Intercomparison Project (CMIP), w którym rozmaite modele klimatyczne stosują te same wymuszenia w celu ich porównania, a także przygotowywania projekcji klimatycznych do raportów IPCC. W tej pracy użyto zarówno CMIP5 (do poprzedniego raportu z 2013 roku) jak i CMIP6 do następnego, który się niedługo ukaże. Rysunek 1 pokazuje po lewej wyniki (średnia to czerwona linia i zasięg wartości na różowo) prognozowania zimowego NAO (grudzień-marzec) z wyprzedzeniem 2 do 9 lat dla lat 1960-2005,w porównaniu do jego obserwowanych wartości (czarna linia). W sumie porażka, nieprawdaż? I parę poprzednich prac tak podobne wyniki oceniło. Smith i koledzy z Met Office zauważyli, że modele przewidują za małą zmienność zimowego NAO (to jedna z konkluzji pracy) postanowili przeskalować wyniki wiązki do wartości zmienności wyników obserwacji. Czyli po prostu powiększyli kilka razy wartość średniej prognozy z wiązki przy pomocą sprytnego triku z powiększeniem ilości prognoz (patrz niżej). I oto jak królik z cylindra ukazuje się całkiem udana prognoza zimowego NAO, przypominam na 2-9 lat naprzód (Rys 1 po prawej).
Wynik niezły ale autorzy poszli jeszcze dalej. Jeśli się sztucznie zwiększy zmienność prognozy to „pompuje” się nie tylko wpływ czynników sterujących wartością NAO ale także i pozostałych, np. rosnącego wymuszenia gazów cieplarnianych. Jak tego uniknąć? Rozwiązali to wybierając podzbiór kilku modeli najlepiej przewidujących NAO (w praktyce tych, które najbliższe były do średniej prognozy NAO dla wiązki). Istotnie polepszyło to jakość prognozy temperatur (Rys. 2), a także pola ciśnień i co ciekawe opadów (te najtrudniej przewidzieć).
Czyli wynika z tego, że NAO jest przewidywalne na wiele lat naprzód. Przyczyną musi być wolna zmienność temperatur wód północnego Atlantyku. Co ciekawe ta ich pod-wiązka dobrana pod względem wierności przebiegu indeksu NAO lepiej przepowiada temperatury powierzchniowe Atlantyku niż cała wiązka, ale na odwrót (pod-wiązka dobrana pod względem wiernego prognozowania temperatur) już trochę gorzej co autorzy interpretują jako znak, że zimowe NAO zależy nie tylko od temperatur Atlantyku (nic nowego: sami powinni to wiedzieć ze swoich poprzednich prac, np. [1]).
Konkluzją artykułu jest to, ze modele mają za małą zmienność indeksów klimatycznych, nie tylko NAO z czego ma płynąć wniosek, ze nie tylko NAO da się przewidzieć na wiele lat naprzód. Ale nie uwierzę póki nie zobaczę wyników.
Mam tez kilka wątpliwości (prawdziwy naukowiec zawsze musi coś wyłapać bo inaczej ie nadawałby się na recenzenta): – Nie rozumiem do końca sposobu w jaki zwiększono zmienność i poprawiono prognozę. Z tekstu zrozumiałem że jest to odrzucenie prognoz na 2 i 3 lata, a zostawienie na 4 do 9 które w jakiś sposób polepsza wyniki. Czy nie powinno być odwrotnie? Z korespondencji z pierwszym autorem wiem, że jest to “użycie lat 2-9 z czterech ostatnich prognoz” ale też nie powiem abym to do końca rozumiał. I dlaczego nie pokazano jak jakość prognozy zależy od liczby lat od jej startu. – Dlaczego prognozy kończą się na roku 2005? Po części pewnie dlatego, ze CMIP5 nie był startowany później ale CMIP6 był. I dlaczego wynik obcięty jest przed rokiem 2010 z rekordowo ujemnym NAO, skoro 2005 + 9 = 2014? – I czy rzeczywiście da się to powtórzyć dla jakiegoś innego indeksu klimatycznego lub parametru meteorologicznego? Ale tak czy siak prognoza NAO na kilka lat naprzód ze skillem 0,79 to szok.
PS. Poprawiłem nieco opis tej “opóźnionej wiązki”, która powiększa zmienność i polepsza prognozę dzięki odpowiedzi jakiej na moje pytanie udzielił Doug Smith, pierwszy autor (7 sierpnia 2020 r.).
[1] Dunstone, N., Smith, D., Scaife, A. et al. Skilful predictions of the winter North Atlantic Oscillation one year ahead. Nature Geosci 9, 809–814 (2016). https://doi.org/10.1038/ngeo2824
[2] Smith, D.M., Scaife, A.A., Eade, R. et al. North Atlantic climate far more predictable than models imply. Nature 583, 796–800 (2020). https://doi.org/10.1038/s41586-020-2525-0
Niniejszy wpis ukazuje się równolegle na blogu Meteomodel
Hits: 917
Bardzo się ucieszyłem, zauważając niespodziewanie nowy wpis w czytniku RSS! Świetny artykuł!
Mam mocne postanowienie pisania tu znów regularnie.
Pisałem już na blogu meteomodel, powtórzę i tutaj, może nieco innymi słowami. Jestem mocno sceptyczny, bo wygląda to bardziej na “zabawę” z wiązkami modeli post factum tak, by odpowiednio dopasować je do zaobserwowanych już danych rzeczywistych. Plus to, co piszesz o podejrzanym zakresie czasowym, obcięciu wyników przed 2010…Na podobnej zasadzie być może dałoby się dobrać i “przeskalować” wiązki GFS, ECMWF itd i pokazać jak to wspaniale przewidziały temperaturę za 10-14 dni.
A przede wszystkim: dlaczego w takim razie państwo z Met Office nie opublikowali prognozy NAO na 10 lat np do 2030, byśmy mogli ją w najbliższych latach weryfikować? Zamiast tego pokazują nam “prognozy” na okres ponad 10 lat wstecz…
Dobrze wyczułeś, ze mi się to też do końca nie podoba. Powody wymieniłem pod koniec wpisu. A co do braku przepowiedni na najbliższą zimę to ono tak robią w każdym artykule. Ten Dunstone i inni (2017) cytowany we wpisie był tego jeszcze lepszym przykładem bo celowo nie pokazali prognozy na następną zimę chociaż wszystkie dane do niej mieli przed złożeniem manuskryptu. Dlatego podejrzewam, ze mają zakaz tego typu prognoz, zapewne po klęsce nieprzewidzenia ujemnego NAO zima 2010 roku (przewidzieli neutralne, a było rekordowo ujemne).
O ile dobrze rozumiem artykuł, nie przewiduje się w nim jednorocznych wartości NAO ani innych indeksów z wyprzedzeniem 2 do 9 lat (ani 4-9 lat), tylko NAO uśrednione w okresie 2-9 lat. Przewidywana jest więc ruchoma średnia 8-letnia, co wyjaśnia dlaczego “wynik obcięty jest przed rokiem 2010”, i dlaczego tego “rekordowo niskiego” NAO z 2010 nie widać w prognozie. Bardziej więc przewidują klimat, niż pogodę.
Wydaje mi się też, że wyniki są polepszone nie po “odrzuceniu prognoz na 2 i 3 lata”, ale po zwiększeniu liczby członków wiązki z 169 do 676 poprzez stworzenie “lagged ensemble”, co zwiększa SNR. Wyjście poza rok 2005 i ograniczenie się tylko do symulacji z CMIP6 oznaczałoby redukcję rozmiaru wiązki o 49, więc można podejrzewać że byłaby ona za mała by dało się pokazać interesujące wyniki.
Prognoz do 2030 roku (zakładając, że chodzi o średnią wartość prognozy na lata 2023-2030) nie można wyprodukować metodą opisaną w artykule z tego prostego powodu, że wymagałoby to wykonania kilkuset symulacji inicjalizowanych obserwacjami z roku 2021. Biorąc pod uwagę to, jak czasochłonne jest zaprojektowanie, skoordynowanie, przygotowanie, przeprowadzenie i opublikowanie wyników takich symulacji, nieuniknione jest, w najlepszym wypadku, kilka lat opóźnienia, więc dopóki ktoś nie wyłoży pieniędzy żeby takie duże wiązki produkować regularnie, najlepsze co można osiągnąć to pokazanie że dla danych z przeszłości ta metoda działa. Z naukowego punktu widzenia przyszłość jest mniej użyteczna, bo to tylko kilka dodatkowych lat potencjalnych danych, a hindcasty można robić dla lat kilkudziesięciu.
No i o ile rozumiem artykuł Dunstone’a, to pokazuje on prognozę na następną zimę, tj. 2016-2017, jest to ostatni punkt na rys. 1b, inicjalizowaną danymi z listopada 2015.
Witam na odnowionym blogu. I cieszę się z bardzo trafionego pytania. Będę musiał chyba dopisać małe sprostowanie bo z korespondencji z pierwszym autorem dowiedziałem się, że chyba nie do końca zrozumiałem tę część z opóźnionymi prognozami. Zapytałem Douga Smita i to i odpowiedział:
To ma więcej sensu niż tak jak to zrozumiałem z tekstu artykułu, dokładniej z dość ciężko napisanego opisu w podrozdziale “Lagged ensemble” w Metodach.
No tak, i oczywiście przyłączam się do głosów przedmówców – dzięki za reaktywację bloga!
A wymowę omawianego artykułu mam za bardziej pesymistyczną niż wynikałoby z tego co mówią autorzy. Tak, fajnie że wymuszony sygnał NAO jest bardziej przewidywalny niż sugerują modele, ale co z tego jeśli trzeba te modele poprawić, albo wyciągać sygnał z absurdalnie dużych wiązek. W obu przypadkach oznacza to, że dekadowe prognozy klimatu są co najmniej o dekadę od praktycznego wdrożenia… I to wciąż przy założeniu, że ten wymuszony sygnał faktycznie tam jest.
Chyba się zgadzamy, że ta praca ma większe znaczenie teoretyczne niż praktyczne. Pokazuje, ze przynajmniej w przypadku NAO jest gdzieś w systemie klimatu informacja o jego przyszłych wartościach. Z tym że jeśli, jak uważam od lat, NAO i AMO to dwie strony tego samego medalu (ale z przesunięciem fazowym) to ta informacja jest po prostu w temperaturach oceanu. Ale co innego to podejrzewać a co innego pokazać, ze nawet modele klimatyczne, które przecież są tylko naszą próbą przedstawienia procesów klimatycznych za pomocą równań matematycznych, tę informacje w jakiejś formie posiadają.
Dodam jeszcze, ze metodologia tej pracy nie spadła z kosmosu. Inspirowana jest, jak wynika z samej jej treści dwoma wcześniejszymi pracami, pochodzącymi z innych zespołów, które trzebaby niestety przeczytać:
Dobrynin, M., Domeisen, D. I. V., Müller, W. A., Bell, L., Brune, S., Bunzel, F., et al. (2018). Improved teleconnection-based dynamical seasonal predictions of boreal winter. Geophysical Research Letters, 45, 3605–3614. https://doi.org/10.1002/2018GL077209
Simpson, I.R., Yeager, S.G., McKinnon, K.A. et al. (2019). Decadal predictability of late winter precipitation in western Europe through an ocean–jet stream connection. Nat. Geosci. 12, 613–619. https://doi.org/10.1038/s41561-019-0391-x