Depozycja atmosferyczna radionuklidów: interpolacja, modele transportu i monitoring

Streszczenie

Stężenie radionuklidu w powietrzu i depozycja na gruncie to dwie różne wielkości. Chmura może przez krótki czas dawać podwyższone Bq/m3, ale dopiero osadzenie na powierzchni tworzy mapę Bq/m2, która wpływa na moc dawki od gruntu, resuspensję, żywność, wodę i monitoring długoterminowy.

Ten artykuł porządkuje drogę od plumy do depozycji: suchą i mokrą depozycję, jednostki, próbkowanie, interpolację wyników pomiarowych, kriging, kokriging z opadem, modele transportu chemicznego oraz ograniczenia map. Celem jest metrologiczne zrozumienie, skąd biorą się mapy skażenia i dlaczego nie wolno traktować ich jak fotografii rzeczywistości.^1,2

Rozszerzenie tematu

Dwie warstwy tego samego zdarzenia

Radionuklid w atmosferze można opisać przynajmniej na dwa sposoby. Pierwszy to stężenie w powietrzu, zwykle w Bq/m3. Drugi to depozycja, czyli aktywność osadzona na jednostce powierzchni, zwykle w Bq/m2.

Te wielkości są powiązane, ale nie są wymienne. Bq/m3 mówi, co w danej chwili może zostać wdychane. Bq/m2 mówi, co pozostało na powierzchni po przejściu chmury albo po opadzie. Dawkę inhalacyjną liczy się z jednej ścieżki, a moc dawki od gruntu i przejście do żywności z drugiej.

To rozróżnienie jest podstawą wszystkich kalkulatorów środowiskowych serwisu. Pluma radiologiczna zaczyna od stężenia w powietrzu i dopiero potem liczy depozycję. Dawka od depozycji gruntu startuje już od Bq/m2. Food-chain dose pyta o przejście z powierzchni lub środowiska do żywności.³

Sucha depozycja

Sucha depozycja zachodzi bez deszczu. Cząstki aerozolu, pyły, krople, jod w różnych postaciach chemicznych albo inne nośniki radionuklidów osiadają na powierzchniach wskutek grawitacji, turbulencji, dyfuzji, impakcji i oddziaływań z podłożem.

W prostym modelu zapisuje się:

D_sucha = C * v_d * t,

gdzie C jest stężeniem w powietrzu, v_d efektywną prędkością depozycji, a t czasem ekspozycji. Taki zapis jest wygodny, bo od razu pokazuje sens jednostek:

Bq/m3 * m/s * s = Bq/m2.

Jedna prędkość depozycji jest jednak skrótem. W rzeczywistości zależy od rozmiaru cząstek, wilgotności, rodzaju powierzchni, roślinności, chropowatości, stabilności warstwy przyziemnej i postaci chemicznej radionuklidu. Dla dydaktyki to dobre pierwsze przybliżenie; dla mapy terenowej to tylko jeden z elementów modelu.

Mokra depozycja

Mokra depozycja zachodzi przez opad. Materiał może zostać wymyty z chmury przez krople deszczu albo śnieg, a następnie trafić na powierzchnię. Dlatego ten sam epizod atmosferyczny może dać niską depozycję w suchym miejscu i wyraźny hot spot tam, gdzie akurat padało.

W kalkulatorze pluma radiologiczna deszcz jest potraktowany jako efektywna prędkość depozycji mokrej. To jest świadomie prosty model: pokazuje kierunek i rząd wielkości, ale nie rozwiązuje chmury, pionowego profilu opadu, zmian w czasie ani rozkładu wielkości kropel.³

W praktyce mapowania mokra depozycja wymaga danych opadowych. Wytyczne do oceny depozycji atmosferycznej w Polsce zalecają wykorzystywanie danych o wysokości opadu jako zmiennej pomocniczej, ponieważ depozycja mokra silnie zależy od ilości opadu docierającego do danego obszaru.²

Jednostki i pułapki jednostek

W jednym raporcie mogą wystąpić różne jednostki:

• Bq/m3 - stężenie aktywności w powietrzu,
• Bq/m2 - depozycja powierzchniowa,
• Bq/kg - aktywność właściwa gleby, żywności albo materiału,
• Bq/L - aktywność w wodzie lub opadzie,
• uSv/h - moc dawki od pola promieniowania,
• mSv - dawka po czasie i po zastosowaniu modelu narażenia.

W chemii atmosfery spotyka się także kg/ha, g/ha lub mg/m2 dla masy zdeponowanego składnika. Wytyczne depozycji chemicznej podają przykład łączenia wysokości opadu i stężenia składnika w opadzie do mapy depozycji w kg/ha.² Dla radionuklidów analogiczna logika jednostkowa działa z aktywnością zamiast masy.

Pułapka polega na tym, że sama aktywność nie mówi jeszcze o dawce. 10 000 Bq/m2 Cs-137, I-131 i Sr-90 nie oznacza tego samego dozymetrycznie. Różne są półokresy, promieniowanie, biokinetyka, przejście do żywności i znaczenie drogi inhalacyjnej lub pokarmowej.⁴

Co mierzy monitoring

Monitoring depozycji może obejmować kilka typów próbek:

• filtry aerozolowe z powietrza,
• próbki opadu atmosferycznego,
• opad całkowity zbierany na powierzchniach depozycyjnych,
• glebę z określonej warstwy,
• roślinność, mleko, żywność albo wodę,
• pomiary mocy dawki w terenie,
• spektrometrię gamma próbek i pomiary radiochemiczne dla nuklidów słabo widocznych w gamma.

Każdy typ próbki odpowiada na inne pytanie. Filtr powietrza mówi o tym, co przechodziło przez atmosferę w czasie poboru. Próbka opadu mówi o mokrej ścieżce usuwania zanieczyszczenia. Gleba integruje historię depozycji, migracji i wcześniejszych zdarzeń. Żywność dodaje biologię i sezon.

Dlatego nie ma jednego "prawdziwego pomiaru skażenia". Jest zestaw wielkości, które trzeba opisać metadanymi: gdzie, kiedy, jaka powierzchnia, jaki czas integracji, jaka głębokość próbki, jaka wydajność detektora, jakie tło, jaka granica wykrywalności i jaki radionuklid.

Od punktów pomiarowych do mapy

Pomiary terenowe są punktowe. Mapa jest powierzchnią ciągłą. Między nimi zawsze stoi model interpolacji.

Najprostsze podejście to przypisanie wartości przez odległość, na przykład IDW - inverse distance weighting. Punkt blisko stacji dostaje większą wagę niż punkt odległy. Taka metoda jest intuicyjna, ale słabo opisuje strukturę błędu i nie wykorzystuje jawnie korelacji przestrzennej.

Kriging idzie dalej. Traktuje mierzoną wielkość jako pole przestrzenne o korelacji zależnej od odległości. Wariogram opisuje, jak różnice między punktami rosną wraz z odległością. Dzięki temu kriging daje nie tylko przewidywaną wartość, ale też informację o niepewności interpolacji.

Wytyczne depozycji atmosferycznej wskazują, że dla danych jednorodnych można stosować kriging zwykły bez transformacji i ze sferyczną funkcją wariogramu. Gdy dane są dodatnio skośne albo zawierają wartości odstające, zalecają rozważenie transformacji logarytmicznej i porównanie wariantów przez statystyki walidacyjne.²

Wartości odstające

Depozycja bardzo często ma wartości odstające. Lokalny deszcz, pył, błąd poboru, rzeczywisty hot spot albo problem analityczny mogą dać wynik dużo większy niż reszta.

Nie wolno automatycznie usuwać takiego punktu, bo może być prawdziwy. Nie wolno też bezkrytycznie interpolować go na duży obszar, bo można stworzyć fałszywą plamę na mapie. Najpierw trzeba rozstrzygnąć, czy wartość jest:

• błędem pomiarowym,
• artefaktem przygotowania próbki,
• realnym lokalnym zjawiskiem,
• skutkiem opadu,
• skutkiem geometrii terenu lub powierzchni.

Wytyczne podkreślają analizę jednorodności danych i możliwość transformacji dla rozkładów skośnych.² W wersji radionuklidowej trzeba dodać jeszcze kontrolę spektrometryczną, tło, interferencje linii gamma i zgodność z innymi nuklidami.

Walidacja interpolacji

Mapa bez walidacji jest rysunkiem. Mapa z walidacją staje się wynikiem pomiarowo-modelowym.

Wytyczne wskazują zestaw parametrów używanych do wyboru wariantu interpolacji, m.in. średni błąd ME, pierwiastek średniego błędu kwadratowego RMSE, średni błąd standardowy ASE, średni standaryzowany błąd MSE i pierwiastek średniego błędu standaryzowanego RMSSE.²

Ich sens można streścić tak:

• ME pokazuje, czy model systematycznie zawyża albo zaniża,
• RMSE karze duże błędy,
• ASE mówi, jaką niepewność przewiduje sam model,
• MSE i RMSSE pomagają sprawdzić, czy niepewność jest realistycznie skalibrowana.

Student powinien zapamiętać, że ładna mapa kolorów nie wystarcza. Trzeba pokazać, jak mapa zachowuje się przy krzyżowej walidacji i gdzie niepewność jest największa.

Siatka mapowania

Wytyczne depozycji atmosferycznej rekomendują mapowanie w siatce 10 x 10 km dla ocen rocznych i miesięcznych.² Taka siatka jest kompromisem. Jest wystarczająco szczegółowa, by pokazać regionalne różnice, ale nie udaje precyzji pojedynczej ulicy.

Dla radionuklidów rozdzielczość trzeba dobierać do skali zjawiska. Globalny fallout historyczny, długookresowa depozycja Cs-137, epizod lokalnego deszczu i uwolnienie z punktowego źródła wymagają innych oczek siatki. Im drobniejsza siatka, tym większe ryzyko, że rysujemy szczegóły nieobecne w danych.

Dlatego mapa depozycji powinna mieć obok siebie co najmniej trzy warstwy: wartości interpolowane, punkty pomiarowe i mapę niepewności. Bez punktów pomiarowych odbiorca nie widzi, gdzie mapa jest podparta danymi, a gdzie głównie modelem.

Kokriging i opad

Kokriging wykorzystuje dodatkową zmienną, która jest skorelowana z wielkością główną i ma lepsze pokrycie przestrzenne. W depozycji mokrej naturalną zmienną pomocniczą jest wysokość opadu.

W wytycznych opisano podejście, w którym użyto danych o stężeniach składników opadu z sieci chemizmu oraz znacznie liczniejszych danych opadowych. Warstwy stężenia i wysokości opadu można połączyć, aby uzyskać bardziej zróżnicowaną mapę mokrej depozycji niż z samych punktów chemicznych.²

Dla radionuklidów ten pomysł jest szczególnie ważny po epizodach atmosferycznych. Jeśli chmura przechodzi nad regionem, ale deszcz występuje tylko miejscami, mapa depozycji może być bardziej podobna do mapy opadu niż do gładkiej plumy. Bez danych opadowych łatwo wygładzić hot spot albo umieścić go w złym miejscu.

Modele transportu chemicznego

Interpolacja zaczyna od pomiarów. Model transportu chemicznego zaczyna od emisji, meteorologii, transportu, przemian i depozycji. Takie modele mogą dawać pole przestrzenne nawet tam, gdzie nie ma stacji, ale wymagają wielu założeń.

W wytycznych jako model rekomendowany do przyszłego wykorzystania wskazano EMEP. Dokument zwraca uwagę na wersje open-source, europejskie zastosowanie i możliwość pracy z emisjami oraz meteorologią w ustalonej siatce.²

Równocześnie model nie zastępuje pomiarów. Wytyczne omawiają testy porównania z danymi monitoringowymi i podkreślają, że wykorzystanie danych modelowych do interpolacji wymaga dalszych testów na nowych danych.² To właściwa postawa metrologiczna: model jest źródłem informacji pomocniczej, a nie automatycznym arbitrem.

Model źródło-do-punktu i mapa z pomiarów

W artykule o Pasquillu-Giffordzie zaczynaliśmy od źródła. Przyjmowaliśmy aktywność, wiatr, stabilność atmosfery i odległość, a potem liczyliśmy stężenie oraz depozycję w punkcie.

W monitoringu depozycji problem jest często odwrotny. Mamy punkty pomiarowe i chcemy dowiedzieć się, jak wygląda pole między nimi. Źródło może być nieznane, rozproszone albo historyczne. Zamiast jednego pióropusza mamy sumę transportu, opadu, migracji, tła i wcześniejszych zdarzeń.

Te dwa podejścia powinny się kontrolować. Model plumy pomaga zrozumieć, czy wzór na mapie jest fizycznie wiarygodny. Pomiary pomagają sprawdzić, czy model nie pominął deszczu, zmiany wiatru albo lokalnego osadzania.

Depozycja a łańcuch pokarmowy

Depozycja na gruncie lub roślinach jest początkiem kolejnych ścieżek. Dla I-131 ważna jest świeża depozycja na paszy i mleko. Dla Cs-137 znaczenie ma dłuższe zachowanie w glebie, roślinach, produktach leśnych i mięsie. Dla Sr-90 istotna jest chemiczna analogia do wapnia i przejście do układu kostnego.

Prosty wynik Bq/m2 nie odpowiada sam z siebie na pytanie o żywność. Potrzebne są współczynniki transferu, czas od depozycji do zbioru, pora roku, typ gleby, opad, praktyki rolnicze i biokinetyka radionuklidu. Dlatego kalkulator food-chain dose pokazuje scenariusze niskiego, centralnego i wysokiego transferu, zamiast jednej pozornie dokładnej liczby.³

To także powód, dla którego mapy depozycji powinny być traktowane jako wejście do modeli środowiskowych, a nie jako końcowa ocena ryzyka. Ta sama depozycja może mieć inne znaczenie zimą, latem, na torfie, na glebie mineralnej, w lesie i na polu uprawnym.

Depozycja a dawka od gruntu

Dla radionuklidów gamma-emisyjnych depozycja powierzchniowa może zostać przeliczona na moc dawki od gruntu. W najprostszym modelu:

Hdot = A_s * k,

gdzie A_s jest depozycją w kBq/m2, a k współczynnikiem mocy dawki dla danego radionuklidu. Potem dawka po czasie to całka z mocy dawki, z poprawkami na rozpad, migrację i ewentualne usuwanie materiału.

To jest nadal uproszczenie. Realne pole gamma zależy od rozkładu głębokościowego w glebie, energii fotonów, geometrii, ekranowania przez budynki, chropowatości terenu i czasu. Dla edukacji ważne jest jednak rozdzielenie: Bq/m2 jest wielkością środowiskową, uSv/h wielkością pola dawki, a mSv wynikiem po czasie i scenariuszu ekspozycji.

Spektrometria gamma i radionuklidy trudniejsze

Wiele map depozycji radionuklidów opiera się na spektrometrii gamma próbek gleby, filtrów lub opadu. Cs-137, Co-60 czy część produktów rozszczepienia ma linie gamma możliwe do identyfikacji. Wtedy problemem stają się kalibracja wydajności, geometria próbki, tło, samopochłanianie i granica wykrywalności.

Nie wszystkie ważne radionuklidy są tak wygodne. Sr-90 jest beta-emiterem i często wymaga radiochemii. Tryt wymaga innych metod. Alfa-emitery wymagają cienkich źródeł i kontroli odzysku. Dlatego mapa "skażenia gamma" nie jest automatycznie mapą wszystkich radionuklidów.

To jest częsty błąd w czytaniu monitoringu. Jeśli raport pokazuje tylko Cs-137, nie mówi bezpośrednio o Sr-90. Jeśli pokazuje moc dawki, może nie rozdzielać radionuklidów. Jeśli pokazuje aktywność w glebie, nie musi mówić o aktywności w żywności.

Tło historyczne i nowe zdarzenie

Depozycja mierzona w glebie jest sumą historii. W Europie wciąż można wykrywać wkłady globalnych testów atmosferycznych i awarii w Czarnobylu; lokalnie mogą dochodzić inne źródła. Nowe zdarzenie trzeba więc odróżnić od tła historycznego.

Pomagają w tym:

• stosunki izotopowe,
• radionuklidy krótkotrwałe,
• porównanie z wcześniejszymi pomiarami,
• datowanie przez półokresy,
• rozkład przestrzenny zgodny z meteorologią,
• obecność lub brak nuklidów wskaźnikowych.

W praktyce to zadanie forensyczne i metrologiczne. Sama mapa wartości nie wystarcza; potrzebny jest profil radionuklidów i czas poboru.

Jak opisywać mapę depozycji

Dobra mapa depozycji powinna mieć kartę metadanych:

• radionuklid lub wskaźnik,
• jednostkę,
• okres poboru lub datę próbki,
• typ próbki i głębokość,
• metodę analityczną,
• granicę wykrywalności,
• liczbę punktów pomiarowych,
• metodę interpolacji,
• rozdzielczość siatki,
• statystyki walidacji,
• mapę lub opis niepewności,
• informację, czy użyto opadu, modelu transportu albo danych pomocniczych.

Bez tego mapa jest trudna do audytu. Może wyglądać profesjonalnie, ale nie wiadomo, czy powstała z gęstej sieci pomiarowej, z kilku punktów, z modelu, czy z połączenia tych źródeł.

Najczęstsze błędy

Pierwszy błąd to mylenie stężenia w powietrzu z depozycją. Krótkotrwały pik Bq/m3 nie musi oznaczać dużego Bq/m2, jeśli cząstki nie osiadły lub epizod był krótki.

Drugi błąd to ignorowanie opadu. W mokrej depozycji mapa deszczu może być równie ważna jak mapa samej chmury.

Trzeci błąd to nadinterpretacja interpolacji. Kriging nie tworzy danych tam, gdzie ich nie ma. Daje najlepsze oszacowanie według przyjętego modelu korelacji przestrzennej.

Czwarty błąd to mieszanie radionuklidów. Mapa Cs-137 nie jest mapą I-131, a moc dawki nie jest składem izotopowym.

Piąty błąd to pomijanie czasu. I-131 po tygodniach znika zupełnie inaczej niż Cs-137; świeża depozycja i stara depozycja mają inne znaczenie.

Podsumowanie

Depozycja atmosferyczna jest pomostem między chmurą w powietrzu a długoterminowym problemem środowiskowym. Fizyka plumy mówi, skąd może wziąć się osadzenie. Monitoring mówi, co rzeczywiście zmierzono. Interpolacja i modele transportu próbują wypełnić przestrzeń między punktami.

Najważniejsza zasada brzmi: mapa depozycji jest wynikiem pomiarowo-modelowym. Trzeba czytać ją razem z jednostką, metodą próbkowania, radionuklidem, opadem, siatką, walidacją i niepewnością.

Sieć monitoringu radionuklidów w Polsce i Europie

Polska posiada rozbudowaną sieć monitoringu skażeń promieniotwórczych powietrza i środowiska, koordynowaną przez kilka instytucji o komplementarnych rolach.

Centralne Laboratorium Ochrony Radiologicznej (CLOR) w Warszawie jest główną instytucją koordynującą radiologiczny monitoring środowiska w Polsce, w tym sieć stacji aerozolowych i opadowych. CLOR prowadzi ciągłe pomiary Cs-137, Sr-90, trytu i innych radionuklidów w powietrzu, opadzie, glebie i żywności. Wyniki są raportowane do IAEA, UNSCEAR i Komisji Europejskiej.

Sieć stacji Wczesnego Wykrywania Skażeń Promieniotwórczych (WCKZ) obejmuje kilkadziesiąt stacji automatycznych, mierzących moc dawki i stężenie radionuklidów w czasie rzeczywistym. System umożliwia szybkie wykrycie epizodów atmosferycznych i generuje alerty dla Państwowej Agencji Atomistyki (PAA).

Wojewódzkie Inspektoraty Ochrony Środowiska (WIOŚ) prowadzą lokalne sieci monitoringu powietrza i pobierają próbki środowiskowe zgodnie z centralnymi wytycznymi metodycznymi.

Europejska sieć EURDEP (EURopean Radiological Data Exchange Platform) łączy krajowe systemy monitoringu z państw Unii Europejskiej i pozwala na wizualizację aktualnej mapy mocy dawki w czasie rzeczywistym. Dane z Polski, Francji, Niemiec, Finlandii i innych krajów trafiają do wspólnego systemu, skąd są dostępne zarówno dla organów regulacyjnych, jak i publicznie.

Program IAEA ALMERA (Analytical Laboratories for the Measurement of Environmental Radioactivity) ocenia zdolności analityczne laboratoriów krajowych poprzez testy biegłości z certyfikowanymi materiałami referencyjnymi. Uczestnictwo w ALMERA jest standardem jakości dla laboratoriów mierzących radionuklidy środowiskowe.

CTBTO IMS Radionuclide Network — Polska hostuje jedną ze stacji IMS sieci CTBTO (System Monitoringu Kompleksowego Układu Zakazu Prób Jądrowych), która pobiera próbki pyłu atmosferycznego i oznacza produkty rozszczepienia i aktywacji. Stacje te mają podstawowe znaczenie dla globalnej weryfikacji traktatu CTBT.

Czarnobyl jako historyczny przypadek mapowania depozycji

Awaria elektrowni jądrowej w Czarnobylu (26 kwietnia 1986 r.) była pierwszym przypadkiem, w którym mapy depozycji atmosferycznej radionuklidów na skalę kontynentalną miały bezpośrednie znaczenie dla decyzji zdrowotnych i rolniczych. Analiza tego zdarzenia jest standardem dydaktycznym w edukacji o depozycji atmosferycznej.

W ciągu pierwszych dni po awarii transportowany materiał obejmował szerokie spektrum radionuklidów: jodki (I-131), cezy (Cs-137, Cs-134), telur (Te-132), ruten (Ru-103, Ru-106), bar (Ba-140) i inne. Każdy z tych nuklidów miał inny półokres, inną prędkość suchej depozycji, inne właściwości mokrej depozycji i inną biologiczną ścieżkę absorpcji. Mapa „Czarnobyla" zależy więc od radionuklidu, daty próbki i metody pomiaru.

Mokra depozycja odegrała decydującą rolę w kształtowaniu map: regiony, nad którymi padało w dniach 27 kwietnia – 4 maja 1986 r., otrzymały wielokrotnie wyższe ładunki Cs-137 niż obszary, które pozostały suche podczas przejścia chmury. To właśnie dlatego mapa depozycji Cs-137 w Europie ma mozaikowy charakter — silne skażenie Szwecji, Finlandii, Bawarii i zachodniej Ukrainy wynikało z miejscowego deszczu, a nie z bliższego położenia względem źródła.

Mapy z 1986 roku były tworzone w oparciu o nierównomierne i niepełne dane pomiarowe — w warunkach kryzysowych sieć stacji była zbyt rzadka, metadane często niekompletne, a pomiary prowadzone różnymi metodami w różnych krajach. Późniejsze retrospektywne analizy (m.in. projekt ETEX, ECDC, UNSCEAR 2008) pokazały rozbieżności między mapami z różnych źródeł, podkreślając jak ważna jest harmonizacja metodologii i dokumentacja metadanych.

Dziś mapy depozycji czarnobylskiej służą jako benchmark dla systemów modelowania atmosferycznego: czy model potrafi retrospektywnie odtworzyć obserwowane rozkłady Cs-137, gdy podane są rzeczywiste dane meteorologiczne i szacunkowe emisje? Modele HYSPLIT, NAME i ARGOS były wielokrotnie walidowane właśnie na danych czarnobylskich.

Modele dyspersji atmosferycznej stosowane w ocenie depozycji

Poza prostym modelem gaussowskim Pasquilla-Gifforda istnieje kilka bardziej zaawansowanych narzędzi modelowania transportu i depozycji radionuklidów:

HYSPLIT (Hybrid Single-Particle Lagrangian Integrated Trajectory): Model opracowany przez NOAA, szeroko stosowany do analizy trajektorii powietrza i dyspersji zanieczyszczeń. W wersji eulero-lagranżowej oblicza zarówno trajektorie powietrza wsteczne (skąd przyszła chmura), jak i prognozowane rozkłady stężeń i depozycji. Dostępny online przez NOAA ARL Ready, bezpłatny, standardowo stosowany w studiach przypadku środowiskowych.

NAME (Numerical Atmospheric-dispersion Modelling Environment): Model lagranżowski opracowany przez UK Met Office. Standardowo stosowany przez Wielką Brytanię i CTBTO do analiz incydentów radiologicznych. Model NAME był użyty do analizy epizodu zanieczyszczeń Ru-106 z 2017 roku, identyfikując prawdopodobne lokalizacje uwolnienia na Uralu.

ARGOS (Accident Reporting and Guidance Operational System): Cywilny system decyzyjny stosowany przez wiele krajów europejskich, łączący modele dispersionowe z bazami danych reaktorów i modułami oceny dawek. Integruje w czasie rzeczywistym dane meteorologiczne, emisje i wyniki monitoringu.

RODOS (Real-time Online Decision Support System): Europejski system wspomagania decyzji w sytuacjach radiologicznych, rozwijany w ramach projektów EU. Integruje modele atmosferyczne, dane z sieci monitoringu, modele transferu do żywności i łańcucha pokarmowego, obliczenia dawek i moduły zarządzania kryzysowego.

Wspólną cechą wszystkich tych modeli jest zależność od jakości danych meteorologicznych (ECMWF, GFS), charakterystyki emisji (czas, masa, skład izotopowy) i danych weryfikacyjnych z monitoringu. Niepewność w charakterystyce emisji jest często dominującym źródłem błędu w prognozach depozycji.

Asymilacja danych: fuzja modelu z pomiarami

Sama interpolacja punktów pomiarowych i sam model transportu są każde z osobna niepełne. Interpolacja nie ma fizyki — nie wie nic o tym, skąd przyszła chmura ani jak wiatr ją ukształtował. Model nie ma rzeczywistych pomiarów — pracuje na szacunkach emisji i modelu meteorologicznym, które mogą mieć błędy.

Asymilacja danych łączy oba podejścia: model dostarcza pierwotne pole przestrzenne (tło modelu), a pomiary kalibrują to pole i korygują błędy systematyczne. Algorytmy takie jak Kalman Filter, 3D-Var i 4D-Var (znane głównie z numerycznej prognozy pogody) mogą być stosowane do oceny pola depozycji przy jednoczesnym użyciu danych z sieci monitoringu i wyników modelu transportu.

W kontekście depozycji radionuklidów asymilacja danych umożliwia:

• korektę biasu modelu na podstawie pomiarów (jeśli model systematycznie zawyża, pomiary go obniżają),
• ekstrapolację w obszarach bez stacji monitoringowych przy jednoczesnym zachowaniu fizycznie wiarygodnego pola,
• retrospektywną analizę zdarzeń przeszłych, gdy pierwotny monitoring był niewystarczający.

Europejski system ENSEMBLE CTBTO stosuje kilka modeli transportu równolegle i łączy ich prognozy, redukując błędy systematyczne wynikające z wyboru jednego modelu.

Trzy przykłady obliczeniowe

Przykład 1: Obliczenie suchej depozycji Cs-137 z chmury

Chmura powietrza zawiera stężenie Cs-137: C = 120 Bq/m³. Prędkość depozycji dla cząstek aerozolu (Cs-137 jest zwykle związany z cząstkami submikronowymi): v_d = 0,001 m/s. Czas przejścia chmury nad punktem: t = 6 h = 21 600 s.

Sucha depozycja: D_s = C × v_d × t = 120 Bq/m³ × 0,001 m/s × 21 600 s = 2 592 Bq/m².

Jeżeli jednocześnie padał deszcz 15 mm z intensywnością scavenging coefficient Λ = 3 × 10⁻⁵ s⁻¹ (typowa wartość dla małych cząstek), to mokra depozycja w czasie t dodaje:

D_m = C × (Λ/λ) × (1 − e^{−λt}) × ... (uproszczony model wykładniczy);

W bardzo prostym przybliżeniu przy stężeniu opadu 500 Bq/L (analogia dla scavenging z chmury) i 15 mm deszczu: D_m = 500 Bq/L × 0,015 m = 7 500 Bq/m².

Łączna depozycja z opadem wielokrotnie przewyższa suchą (7 500 + 2 592 ≈ 10 000 Bq/m²), co ilustruje decydującą rolę opadów w kształtowaniu mapy skażeń.

Przykład 2: Przeliczenie depozycji Cs-137 na moc dawki zewnętrznej

Depozycja Cs-137 wynosi A_s = 10 kBq/m² = 10 000 Bq/m². Współczynnik konwersji dawki zewnętrznej dla Cs-137 (gamma 662 keV, zaktualizowany dla geometrii terenowej): k ≈ 2,5 × 10⁻⁷ Sv/h na Bq/m² (wartość przybliżona dla równomiernej depozycji na gruncie, osoba stojąca).

Moc dawki: Ḣ = A_s × k = 10 000 Bq/m² × 2,5 × 10⁻⁷ Sv/h·(Bq/m²)⁻¹ = 2,5 × 10⁻³ Sv/h = 2,5 mSv/h.

Dawka roczna przy stałym pobycie na zewnątrz: D_rok = 2,5 mSv/h × 8760 h = 21 900 mSv/rok = 21,9 Sv/rok.

Oczywiście w rzeczywistości: (1) moc dawki zmniejsza się z czasem przez rozpad Cs-137 (T₁/₂ = 30,17 lat, więc zmiany są powolne w skali roku), (2) ludzie spędzają dużą część czasu w budynkach (czynnik ochrony K_shelter ≈ 0,2–0,5), (3) depozycja migruje w głąb gleby redukując emisję powierzchniową. Dla realistycznej oceny dawki te czynniki muszą być uwzględnione.

Przykład 3: Kriging — ocena niepewności interpolacji

Sieć 15 stacji pomiarowych mierzyła depozycję Cs-137. Variogram sferyczny dopasowany do danych ma parametry: nugget c₀ = 50 (Bq/m²)², sill c₁ = 500 (Bq/m²)², zasięg a = 120 km.

Dla punktu predykcji odległego od najbliższej stacji o d = 80 km, kriging zwykły daje niepewność interpolacji σ_K² ≈ c₀ + c₁ × h(d/a) = 50 + 500 × (1,5 × 80/120 − 0,5 × (80/120)³) = 50 + 500 × (1 − 0,630 + 0,074) = 50 + 500 × 0,856 = 50 + 428 = 478 (Bq/m²)².

Odchylenie standardowe: σ_K = √478 ≈ 21,9 Bq/m².

Dla punktu tuż obok stacji (d = 2 km): σ_K² ≈ c₀ + c₁ × h(2/120) ≈ 50 + 500 × 0,025 = 50 + 12,5 ≈ 62,5, σ_K ≈ 7,9 Bq/m².

Wniosek: niepewność interpolacji wzrasta z odległością od stacji — jest 3 razy większa dla punktów oddanych (80 km) niż blisko stacji (2 km). Mapa niepewności kriggingu powinna towarzyszyć mapie wartości interpolowanych jako obowiązkowy element dokumentacji.

Pytania otwarte

1. W jaki sposób sezonowość (zima, lato, jesień) wpływa na prędkość suchej depozycji Cs-137? Jak zmienia to mapę depozycji z jednego epizodu atmosferycznego w różnych porach roku?

2. Jak poprawnie rozróżnić historyczne tło depozycji (testy jądrowe z lat 50.–60. i Czarnobyl 1986) od ewentualnego nowego zdarzenia mierzonego w glebie? Jakie stosunki izotopowe lub dodatkowe radionuklidy mogą w tym pomóc?

3. Kriging zakłada stacjonarność przestrzenną (taka sama struktura korelacji w całym obszarze). Jak to założenie jest naruszone w przypadku depozycji mokrej zależnej od lokalnych opadów? Jakie metody mogą tę niejednorodność uwzględnić?

4. W jaki sposób pionowy profil głębokościowy Cs-137 w glebie (zależny od infiltracji wody i bioturbacji) wpływa na obliczoną moc dawki zewnętrznej? Jak zmieniałaby się moc dawki przy tej samej depozycji powierzchniowej po 30 latach migracji radionuklidów w głąb profilu glebowego?

5. Jak porównać wyniki dwóch map depozycji zrobionych różnymi metodami (kriging vs. model transportu) bez znajomości "prawdziwej" odpowiedzi? Jakie kryterium oceny zastosować?

6. Jakie są ograniczenia detektora w samolocie lub dronie do mapowania mocy dawki jako proxy dla depozycji Cs-137? Jak przeliczać pomiary lotnicze na depozycję powierzchniową?

7. Dlaczego depozycja Sr-90 jest trudniejsza do zmapowania niż depozycja Cs-137, mimo podobnej historii atmosferycznej emisji? Co to mówi o roli spektrometrii gamma w monitoringu skażeń?

8. Jak system RODOS lub ARGOS integruje w czasie rzeczywistym dane z sieci monitoringu i modelu transportu podczas sytuacji awaryjnej? Jakie są podstawowe ograniczenia takich systemów w pierwszych godzinach zdarzenia?

Podsumowanie dydaktyczne

1. Dwie skale opisu — stężenie w powietrzu (Bq/m³) i depozycja powierzchniowa (Bq/m²) to różne pytania fizyczne i różne ścieżki dawkowania. Pierwsze mówi o wdychaniu teraz, drugie o osiadaniu na powierzchni i długoterminowym problemie środowiskowym. Nie można ich bezpośrednio porównywać bez modelu fizycznego łączącego oba opisy.

2. Dominacja mokrej depozycji — opad deszczu lub śniegu może wielokrotnie zwiększyć lokalną depozycję w stosunku do depozycji suchej. Mapa depozycji radionuklidów z epizodów atmosferycznych jest więc w dużej mierze mapą opadu podczas przejścia chmury, a nie samej chmury. Dane opadowe są kluczową zmienną pomocniczą w kokrigingu.

3. Interpolacja jako model — kriging i inne metody interpolacji przestrzennej nie tworzą danych. Dają najlepsze przewidywanie według przyjętego modelu korelacji przestrzennej, wraz z ilościową oceną niepewności. Mapa niepewności jest równie ważna jak mapa wartości i powinna być zawsze prezentowana razem z nią.

4. Walidacja jako obowiązek — parametry statystyczne (ME, RMSE, ASE, MSE, RMSSE) są standardowymi miarami jakości interpolacji. Mapa bez tych informacji jest czysto ilustracyjna i nie może być podstawą decyzji ochronnych.

5. Radionuklidy nie są wymienne — mapa Cs-137 nie zastępuje mapy Sr-90, I-131 czy Pu-239. Każdy nuklear ma inną prędkość depozycji, inną mobilność w glebie, inną biologiczną ścieżkę i inne znaczenie dozymetryczne. Monitoring musi rozróżniać składniki; raport powinien zawsze podawać, który nuklid jest mierzony.

6. Tło historyczne — europejskie gleby nadal zawierają Cs-137 z testów atmosferycznych (1945–1963) i awarii w Czarnobylu (1986). Wykrycie Cs-137 w glebie nie jest samo w sobie dowodem na nowe zdarzenie. Odróżnienie tła od nowego sygnału wymaga analizy stosunków izotopowych, wartości bezwzględnych i profilu głębokościowego.

7. Model transportu a pomiary — modele HYSPLIT, NAME i ARGOS dostarczają fizycznie spójnego pola przestrzennego, ale mogą mieć błędy wynikające z niedokładności emisji lub meteorologii. Pomiary monitoringowe służą do walidacji i kalibracji modelu. Asymilacja danych łączy oba podejścia formalnie, minimalizując błędy systematyczne każdego z nich.

8. Polska infrastruktura i europejska sieć — CLOR, WCKZ, WIOŚ i sieć EURDEP tworzą wielopoziomowy system detekcji i mapowania radionuklidów atmosferycznych. Uczestnictwo w sieci CTBTO IMS i programie IAEA ALMERA gwarantuje porównywalność i jakość wyników. Rozumienie tej sieci jest niezbędne do właściwej interpretacji publicznie dostępnych danych o stężeniu radionuklidów w Polsce i Europie.

Rutenowy incydent 2017 jako współczesny przykład

Doskonałym współczesnym przykładem problemów z mapowaniem depozycji atmosferycznej jest epizod podwyższonego stężenia Ru-106 (ruten-106) w Europie, wykryty we wrześniu i październiku 2017 roku. Ru-106 ma T₁/₂ = 373,6 dnia i nie jest naturalnym składnikiem atmosfery — jego obecność wskazuje jednoznacznie na źródło sztuczne. Charakterystycznie, Ru-106 nie ma produktów rozszczepienia ze stabilnych izotopów, jest tworzony przez długotrwałą aktywację neutronową w reaktorze jądrowym.

Pomiary stacji IMS CTBTO i krajowych sieci monitoringowych (Francja, Niemcy, Austria, Rumunia, Polska i inne kraje) zarejestrowały podwyższone stężenia Ru-106 w powietrzu na poziomie od kilku do kilkudziesięciu mBq/m³. Na podstawie analizy trajektorii wstecznych (HYSPLIT, NAME, ARGOS) zidentyfikowano prawdopodobny obszar emisji na Uralu w Rosji, a późniejsze doniesienia wskazywały na zakład przerobu paliwa jądrowego Majak.

Przypadek Ru-106 jest dydaktycznie wartościowy z kilku powodów:

• Wszystkie stacje IMS CTBTO zarejestrowały tylko Ru-106, bez towarzyszących nuklidów charakterystycznych dla awarii reaktora (brak Cs-137, Cs-134, I-131 na podwyższonym tle). To wskazywało na uwolnienie izolowanego Ru-106, zgodne ze scenariuszem przerobu chemicznego zużytego paliwa.
• Rekonstrukcja źródła z danych pomiarowych wielostacji jest przykładem odwrotnego zadania transportowego (inverse modelling): mając rozkłady stężeń w kilkudziesięciu punktach Europy, szacuje się metodami optymalizacyjnymi lokalizację i moc emisji.
• Brak oficjalnego przyznania się przez żadne państwo do emisji ilustruje trudność w politycznym zastosowaniu wyników monitoringu naukowego — nawet przy przekonujących dowodach analitycznych.

Incydent Ru-106 pokazał też siłę sieci monitoringowej: setki stacji w różnych krajach, każda działająca zgodnie z tymi samymi protokołami analitycznymi i metadanymi, stworzyły razem spójny obraz transnarodowego epizodu.

Przestrzenne korelacje i izotropowość pola depozycji

Jednym z kluczowych założeń kriggingu jest izotropowość: korelacja przestrzenna zależy tylko od odległości między punktami, nie od kierunku. W praktyce depozycja atmosferyczna jest anizotropowa — chmura jest transportowana w określonym kierunku przez dominujące wiatry, więc rozkład stężeń ma strukturę wydłużoną wzdłuż kierunku wiatru.

W przypadku silnej anizotropii można stosować kriging anizotropowy, gdzie variogram jest elliptyczny zamiast kolistego: różna zasięg korelacji w kierunku wiatru i prostopadle do wiatru. Bardziej zaawansowane podejście to sekwencyjna symulacja Gaussa lub symulacja sekwencyjna wskaźnikowa, które generują wiele równie prawdopodobnych realizacji pola depozycji — co pozwala na ocenę rozkładu niepewności (nie tylko jego wartości oczekiwanej, jak w kriggingu zwykłym).

Dla celów dydaktycznych ważna jest świadomość, że kriging zwykły daje jedno pole — „najlepszą" predykcję. Symulacja warunkowa daje wiele pól o tej samej wiarygodności statystycznej — co jest fundamentalną różnicą przy analizie scenariuszowej i ocenie ryzyka w obszarach bez pomiarów. W praktyce regulacyjnej ta różnica ma bezpośrednie konsekwencje: decyzje oparte na jednej mapie środkowej mogą być nieostrożne, gdy rozkład niepewności jest szeroki lub asymetryczny — jak to często bywa w warunkach rzadkiej sieci monitoringowej i nieregularnych opadów.

Dodatkowe materiały multimedialne

Ćwiczenia praktyczne

Pierwsze ćwiczenie: dla C = 80 Bq/m3, v_d = 0,0015 m/s i czasu 40 min obliczyć suchą depozycję w Bq/m2. Następnie podwoić czas i sprawdzić, czy wynik zmienia się liniowo.

Drugie ćwiczenie: przyjmij wysokość opadu 20 mm i stężenie radionuklidu w opadzie 5 Bq/L. Oblicz orientacyjną depozycję w Bq/m2, korzystając z faktu, że 1 mm opadu odpowiada 1 L/m2.

Trzecie ćwiczenie: przygotować tabelę pięciu punktów syntetycznych z depozycją Cs-137 i wskazać, który punkt jest potencjalnym odstającym. Zaproponować trzy testy kontrolne, zanim punkt zostanie użyty w interpolacji.

Czwarte ćwiczenie: dla trzech metod mapowania - IDW, kriging zwykły i kokriging z opadem - opisać po jednym założeniu, jednej zalecie i jednym ryzyku nadinterpretacji.

Piąte ćwiczenie: wziąć wynik 12 kBq/m2 dla Cs-137 i opisać, jakich dodatkowych danych potrzeba, aby oszacować zewnętrzną dawkę roczną. Nie wykonywać rachunku bez wskazania brakujących założeń.

Szóste ćwiczenie: porównać mapę depozycji i mapę punktów pomiarowych. Wskazać obszary, gdzie wynik jest dobrze podparty danymi, oraz obszary, gdzie dominuje model interpolacji.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły

• Pasquill-Gifford i modele plumy: od pogody do dawki
• Wielkości dozymetryczne: dawka pochłonięta, równoważna, efektywna i operacyjna
• Naturalna promieniotwórczość żywności i materiałów budowlanych: od widma do Bq/kg
• Opad promieniotwórczy: fallout i mechanizm powstawania
• Pluma skażenia radiologicznego
• Dawka od depozycji gruntu