04 · A vizsgálat részletes bemutatása

Az adatoktól a következtetésig

Rövid válasz: a 2012–2025 közötti meteorológiai adatok alapján nem találtunk olyan megbízható statisztikai mintázatot, amely alátámasztaná, hogy a JÉGER 2018-as országos bevezetése fokozta volna Magyarország aszályosságát. A működési szezonban kapott becslések az aszályosító hipotézissel ellentétes irányba mutattak, de a módszertani ellenőrzések miatt ezeket nem tekintjük oksági hatásnak.

Miért kezdtük el vizsgálni?

Az aszály Magyarországon nem elvont kérdés. Közvetlenül érinti a termést, a talaj vízkészletét, a természetes élőhelyeket és a vízgazdálkodást. A különösen száraz években érthetően felerősödik az igény arra, hogy az emberek okokat találjanak. A JÉGER országos jégkármérséklő rendszer körül ezért visszatérő társadalmi vita alakult ki. Ennek egyik állítása szerint a rendszerben alkalmazott ezüst-jodid „elveszi az esőt”, vagyis hozzájárul a szárazsághoz.

Ezt az aggodalmat nem érdemes sem gúnnyal, sem előre gyártott válasszal kezelni. A hipotézis vizsgálható részeit adatokkal kell összevetni. Elemzésünk nem laboratóriumi kísérlet, és nem követi a kibocsátott ezüst-jodidot. Azt kérdezi meg, látható-e a 2018-as országos bevezetés után Magyarország meteorológiai mutatóiban a kontrollterületekhez viszonyítva olyan változás, amely összhangban lenne egy erős, országos aszályosító hatással.

Pontosan mit kérdeztünk?

2018 után, hasonló európai kontrollterületekhez képest Magyarországon csökkent-e a csapadék, hosszabbá vált-e a leghosszabb összefüggő száraz időszak, illetve csökkent-e a felső talajréteg nedvessége? Elsősorban az áprilistól szeptemberig tartó aktív szezont vizsgáltuk. Ugyanezt elvégeztük az október–márciusi téli ellenőrző időszakra, amikor nem várható JÉGER-hatás (placebo-időszak) adataival is.

Az előzetes várakozást előre rögzíthetjük: ha a JÉGER országos bevezetése erős aszályosító hatással járt volna, akkor az aktív szezonban a kontrollokhoz képest kevesebb csapadékot, hosszabb száraz periódust és alacsonyabb talajnedvességet várnánk. Ez irányhipotézis, nem bizonyíték.

Milyen adatokat használtunk?

Az elemzés az Európai Unió Copernicus programjához tartozó ERA5-Land reanalízisre épül. A reanalízis mérési információkat és időjárási modelleket egyesít, így térben és időben egységes meteorológiai adatsort ad. A vizsgálat 2012-től 2025 végéig tart, 0,1 fokos térbeli rácson. Az órás ERA5-Land állományokból napi értékek készültek; ezekből számoltuk a szezonális mutatókat.

Az audit szerint mind a 168 vizsgált hónap nyers és napi feldolgozott állománya rendelkezésre állt és érvényes volt. Magyarország rácscelláit Ausztria, Szlovákia, Románia, Horvátország, Szerbia és Szlovénia potenciális kontrollcelláihoz hasonlítottuk.

Miért kell kontrollterület?

Nem elég összevetni Magyarország 2018 előtti és utáni átlagát. Képzeljünk el két, hasonló földet: az egyiken történik egy beavatkozás, a másikon nem, miközben mindkettőt ugyanaz a hőhullám vagy regionális csapadékhiány éri. A második föld segít elkülöníteni a közös időjárási változást a beavatkozással egy időben megjelenő eltéréstől.

Az igazi kérdés így hangzik: mi történt volna Magyarországgal 2018 után JÉGER nélkül? Ezt az alternatív történetet nem figyelhetjük meg. Hasonló kontrollterületekkel közelítjük, de a közelítés minősége döntő. Ha a kontroll már 2018 előtt sem mozog együtt Magyarországgal, az utána mért különbséget nem lehet biztonságosan a bevezetéshez rendelni.

Az analógia határa is fontos. Két mezőgazdasági tábla esetében könnyebb lehet hasonló párt találni, mint egy egész ország sok száz rácscellájához. A domborzat, a tengerektől való távolság és az időjárási rendszerek útvonala olyan különbségeket okozhat, amelyeket néhány átlagos klímamutató nem fog meg. A kontrollterületek kiválasztása (matching) ezért nem garancia, hanem ellenőrizendő kiindulópont.

Miért 2018 a választóvonal?

2018 a JÉGER országos intézményi indulásának éve. A modellben ez a „kezelés”. Ez azonban nem azonos a tényleges ezüst-jodid-expozícióval. Nem tudjuk belőle, hogy egy adott napon melyik generátor mennyi ideig működött, mennyi anyagot használt, merre fújt a szél, vagy egy zivatarfelhő milyen koncentrációval találkozott. A vizsgálat ezért a kontrollterületekhez viszonyított változásvizsgálat, nem közvetlen fizikai hatásmérés.

Három mutató, amit különösen figyeltünk

1. Szezonális csapadék: az aktív szezon vagy a JÉGER-hatás nélküli téli ellenőrző időszak napi csapadékainak összege milliméterben.
2. CDD: a leghosszabb egymást követő száraz napok száma a szezonban. A nagyobb érték hosszabb megszakítás nélküli száraz periódust jelent.
3. Felső talajréteg nedvessége: a talaj legfelső rétegének átlagos térfogati víztartalma. Ez nem azonos a teljes gyökérzóna vízkészletével, de hasznos felszínközeli mutató.

Mit találtunk?

A magyarországi és a kontrollterületi változások különbségét mérő párosított modell (matched DID) az aktív szezonban a következő, kontrollterületekhez viszonyított eltéréseket becsülte 2018 után:

• 17,1 mm-rel több szezonális csapadék (95%-os konfidenciaintervallum: 14,3–19,9 mm);
• 0,71 nappal rövidebb leghosszabb száraz periódus (a becslés előjele negatív; 95%-os intervallum: -0,90–-0,52 nap);
• 0,00356-tal magasabb felső talajnedvesség (95%-os intervallum: 0,00292–0,00421).

Mindhárom előjel az aszályosító hipotézissel ellentétes. Ha erős országos aszályosítás lenne a meghatározó jel, első közelítésben nem ezt a mintát várnánk. Ez fontos leíró eredmény, de a „nem ezt látjuk” még nem azonos azzal, hogy a JÉGER okozta volna a kedvezőbb értékeket.

A konfidenciaintervallumok a statisztikai bizonytalanságot mutatják az adott modell feltételei mellett. Nem tartalmazzák automatikusan a rossz kontrollválasztásból, a hiányzó expozíciós adatokból vagy az eltérő előtrendből eredő bizonytalanságot. A szűk intervallum és a nagyon kis p-érték ezért sem teszi okságivá a becslést.

Akkor a JÉGER több esőt okoz?

Nem ezt mutattuk ki. A magyarországi és a kontrollterületi változások különbsége (DID-becslés) egy feltételekhez kötött statisztikai együttmozgást számszerűsít: hogyan változott Magyarország és a párosított kontrollok különbsége a JÉGER 2018-as országos indulása körül. Oksági állításhoz többek között hihető párhuzamos előtrend, megfelelő kontrollilleszkedés és olyan téli ellenőrző eredmény kellene, amely nem jelez hasonló „hatást” ott, ahol annak a feltételezett mechanizmus szerint nem kellene jelentkeznie.

A pozitív csapadékbecslés tehát nem bizonyítja, hogy a JÉGER esőt hoz létre. Lehet regionális meteorológiai különbség, nem megfelelő kontroll, eltérő trend vagy más, a modellel nem kezelt folyamat eredménye. A helyes állítás szűkebb: az aktív szezon relatív eltérései nem az országos aszályosító hipotézis által várt irányba mutattak.

Három szintet következetesen különválasztunk. Az adat éves és térbeli meteorológiai értékeket mutat. A modell ezek feltételes, kontrollokhoz viszonyított különbségét becsüli. Az oksági következtetés csak akkor lenne indokolt, ha az összehasonlítás azonosítási feltételei is hihetőek. Itt az első két szintről tudunk beszélni; a harmadikhoz a diagnosztika nem ad elegendő alapot.

Megpróbáltuk megcáfolni a saját eredményünket

A téli ellenőrző időszak, amikor nem várható JÉGER-hatás (placebo-időszak), azért fontos. Ha a feltételezett kapcsolat a JÉGER működési szezonjához kötődik, akkor az október–márciusi időszakban megjelenő nagy „hatás” gyanús. Arra utalhat, hogy a modell által mért eltérést nem maga a kezelés, hanem a területek közti más különbség okozza.

A téli ellenőrző időszak becslései jelentősek voltak: csapadékra -32,1 mm, a leghosszabb száraz periódusra (CDD) +0,96 nap, talajnedvességre -0,00403. Ezek éppen aszályosabb irányba mutatnak, noha nem ez a fő működési szezon. Ez nem „bizonyítja”, hogy a téli ellenőrző eredmény valódi hatás. Ellenkezőleg: figyelmeztet, hogy a modell szezonális összehasonlítása olyan eltéréseket is kezelésként érzékelhet, amelyeknek más oka van.

Az előtrend-probléma

A magyarországi és a kontrollterületi változások különbségét mérő modell (DID) egyik kulcsfeltétele, hogy Magyarország és a kontrollterületek már 2018 előtt is hasonló pályán mozogjanak. Az évenkénti változásokat bemutató elemzés (event study) ezt évenként vizsgálja egy 2018 előtti referenciaévhez képest. Az eredmény szerint legalább egy, valójában több 2018 előtti csapadékkoefficiens szignifikánsan eltér nullától. Az összehasonlított területek kiinduló hasonlóságát mérő táblában (matching-egyensúly) a kilenc illesztési változóból 8 abszolút standardizált különbsége nagyobb 0,1-nél. Ez gyenge egyensúlyra utal.

Ez nem apró technikai lábjegyzet. Azt jelenti, hogy a kontroll nem épít fel kellően meggyőző alternatív történetet Magyarország számára. Emiatt az aktív szezonban számított a két területcsoport eltérő változását számszerűsítő együtthatókat nem értelmezzük oksági hatásként.

Az évenkénti változásokat bemutató elemzés 2018 utáni koefficiensei is erősen ingadoznak. Ez összhangban van az egyes évek nagy természetes meteorológiai változékonyságával. Egy átlagos post-koefficiens ezt összesűríti, de nem oldja meg az előidőszakban már jelen lévő különbségeket. Ezért az ábra a teljes előidőszakot és a 95%-os intervallumokat is megmutatja.

Adatokból összeállított mesterséges kontrollterület (synthetic control)

Külön módszerrel donorországok súlyozott kombinációjából „szintetikus Magyarországot” is készítettünk. A vizsgált időszakban mért eltérés (post-gap) átlagosan -180 mm volt, vagyis Magyarország értéke alacsonyabb volt a szintetikus kontrollnál. Ezt önmagában könnyű lenne erős ellenérvként bemutatni, de az összevetés már 2018 előtt is rosszul illeszkedett: a vizsgálat előtti illeszkedés hibája (pre-RMSPE) 187 mm. A nagy előilleszkedési hiba miatt a vizsgált időszakban mért eltérést nem tekintjük megbízható bizonyítéknak.

Ebben a futásban a hat donorország egyenlő súlyt kapott. Ez technikailag érvényes eredmény, de a magas előhiba mellett nem jelent jó kontrollt. Ha a módszer már a bevezetés előtt rendszeresen jóval csapadékosabb „szintetikus Magyarországot” rajzol, a bevezetés utáni negatív rés nem választható el ettől a korábbi szintkülönbségtől.

Mit mondhatunk ezek után?

Nem találtunk megbízható bizonyítékot arra, hogy a JÉGER fokozná az aszályt. A működési szezon becslései az aszályosító hipotézissel ellentétes irányúak, de a téli ellenőrző időszak, a 2018 előtti eltérő változások és az adatokból összeállított mesterséges kontrollterület (synthetic control) gyenge illeszkedése megakadályozzák az oksági értelmezést.

Ez nem ugyanaz az állítás, mint hogy bebizonyítottuk: a JÉGER-nek vagy az ezüst-jodidnak semmilyen meteorológiai hatása nincs. A jelenlegi elemzés a JÉGER országos indulása körüli, kontrollterületekhez viszonyított változást vizsgál, és nem képes egyedi beavatkozásokat, zivatarokat vagy szélirány szerinti, a generátoroktól lefelé fekvő területek kitettségét (downwind expozíciót) azonosítani. Ugyanígy azt sem bizonyítja, hogy a JÉGER növeli a csapadékot.

A „nem találtunk megbízható bizonyítékot” nem kitérő, hanem a bizonyíték erejének pontos leírása. Különbözik attól, hogy „bizonyítottuk a hatás hiányát”. Az utóbbihoz olyan vizsgálati erő és expozíciós pontosság kellene, amely mellett egy érdemi hatás hiánya elkülöníthető a mérési és módszertani bizonytalanságtól.

Mit nem tudunk még?

Nem állt rendelkezésre a generátorok pontos koordinátája, napi üzemideje, az ezüst-jodid-felhasználás időpontja és mennyisége, a kezelt zivatarhelyzetek listája, a szélirány alapján számított, a generátoroktól lefelé fekvő területek kitettsége (downwind expozíció), illetve a tényleges légköri ezüst-jodid-koncentráció. Ezek nélkül az egyes meteorológiai események fizikai kitettsége nem rekonstruálható.

Mi kellene egy erősebb vizsgálathoz?

Egy következő elemzésnek össze kellene kapcsolnia a generátorok helyét és üzemnaplóját a radar- és meteorológiai eseményekkel. A szélmezők alapján eseményenként a generátoroktól szélirányban lefelé és felfelé fekvő (downwind és upwind) területeket lehetne kijelölni, majd az összehasonlítható zivatarok csapadékképződését vizsgálni. Ez közelebb vinné az elemzést a tényleges expozícióhoz, és lehetőséget adna dózis–válasz, távolság- és időzítési minták ellenőrzésére.

Az erősebb design több, egymást ellenőrző összehasonlításból állhatna: azonos meteorológiai helyzetű működő és nem működő napokból; szél alatti és szél feletti területekből; eltérő távolságokból; valamint negatív kontrollként olyan órákból vagy térségekből, amelyek nem lehettek kitéve. A radar segítene a folyamat időzítésében, felszíni állomások pedig független ellenőrzést adhatnának.

Hogyan érdemes olvasni az eredményt?

Az ábrákon szereplő nagy számok a magyarországi és kontrollterületi változások különbségét mérő fő modell (matched DID) becslései, nem közvetlenül megfigyelt „JÉGER-hatások”. A Magyarország–kontroll idősorok segítenek látni az éves ingadozást; az évenkénti változásokat bemutató elemzés (event study) azt ellenőrzi, hogy a különbségek mikor jelentek meg; a téli ellenőrzés (placebo) azt kérdezi, kapunk-e látszólagos eredményt olyan szezonban is, ahol a feltételezett mechanizmus alapján nem ezt várnánk. Egyetlen panel sem értelmezhető a többi nélkül.

A bizonyíték mérlegelésekor nem egyszerű szavazást végzünk a módszerek között. Nem mondjuk azt, hogy három aktív szezonos előjel „legyőz” egy téli ellenőrzést vagy egy rossz előtrendet. A diagnosztikai hibák azt érintik, hogy a fő becslés egyáltalán tekinthető-e a bevezetés következményének. Ezért nagyobb súllyal esnek latba, mint a fő koefficiensek statisztikai szignifikanciája.

Az eredmény ugyanakkor nem értéktelen. Egy erős, egyszerű és országosan egységes aszályosító jel esetén várnánk, hogy a fő aktív szezonos mutatók legalább következetesen az aszályosabb irányba mozduljanak a kontrollokhoz képest. Ez nem történt meg. A vizsgálat tehát nem támasztja alá ezt az egyszerű országos hipotézist, miközben nem tudja kizárni a kisebb, helyi, időjárási helyzethez kötött vagy másként definiált hatásokat.

Végül az ERA5-Land reanalízis erőssége az egységes térbeli lefedettség, de nem helyettesít minden helyi mérést. A rácsátlag kisimíthat lokális szélsőségeket, és a talajnedvesség modellalapú mennyiség. Egy későbbi eseményszintű vizsgálatban radar-, állomási és üzemeltetési adatok együttes használata adna erősebb ellenőrzést.