Pasquill-Gifford i modele plumy: od pogody do dawki

Streszczenie

Wyobraź sobie punktowe uwolnienie znacznika do powietrza: ten sam materiał może po kilku minutach tworzyć szeroką, rozcieńczoną smugę albo wąski pas podwyższonego stężenia przy ziemi. O różnicy decyduje nie tylko ilość uwolnionej aktywności, lecz także wiatr, pionowe mieszanie atmosfery, wysokość uwolnienia, deszcz i droga narażenia.

Ten artykuł wyjaśnia model plumy Pasquilla-Gifforda jako narzędzie pierwszego przybliżenia: od klasy stabilności A-F i szerokości sigma_y, sigma_z, przez stężenie w powietrzu, depozycję suchą i mokrą, aż do dawki inhalacyjnej i mocy dawki od skażonego gruntu. Nie jest to instrukcja oceny realnego zdarzenia radiologicznego, tylko opis dydaktycznego modelu, który pomaga zrozumieć, dlaczego wyniki kalkulatora tak silnie zależą od meteorologii.^1,2

Rozszerzenie tematu

Co naprawdę modelujemy

Najprostsze pytanie brzmi: ile radionuklidu znajdzie się w powietrzu w punkcie obserwacji? To pytanie szybko rozpada się na kilka mniejszych. Trzeba wiedzieć, ile materiału weszło do chmury, z jaką prędkością wieje wiatr, jak mocno atmosfera miesza chmurę w bok i w pionie, czy źródło jest przy ziemi czy wyżej, oraz czy cząstki osiadają na powierzchni.

W praktyce dydaktycznej nie zaczyna się od pełnego modelu numerycznej pogody. Zaczyna się od płaskiego, stacjonarnego przybliżenia: wiatr wieje w jednym kierunku, źródło jest punktowe, teren jest jednorodny, a rozkład stężenia poprzecznie do wiatru i pionowo ma kształt gaussowski. Taki model nie opisuje zabudowy, dolin, frontów atmosferycznych ani zmiennego deszczu. Daje jednak bardzo użyteczny język: x oznacza odległość z wiatrem, y odchylenie boczne, u prędkość wiatru, a sigma_y i sigma_z mówią, jak szeroka jest chmura.

To jest sens modelu Pasquilla-Gifforda w serwisie. Nie odpowiada on na pytanie "jaka jest prawdziwa mapa skażenia w terenie", lecz na pytanie "jak zmienia się skala wyniku, gdy atmosfera jest chwiejna, neutralna albo stabilna". Taki poziom jest wystarczający do nauki zależności i do interpretacji prostych kalkulatorów środowiskowych.^1,3

Pluma, nie opad po wybuchu

Słowo "pluma" oznacza tu pióropusz zanieczyszczenia niesiony wiatrem. Może to być epizodowe uwolnienie z instalacji, aerozol w ćwiczeniu modelowym, chmura po awarii technicznej albo abstrakcyjny znacznik używany do nauki dyspersji.

Nie należy automatycznie utożsamiać tego z lokalnym falloutem po naziemnej detonacji jądrowej. Fallout po wybuchu zależy od wyniesienia gleby, rozkładu rozmiarów cząstek, wysokości chmury, frakcjonowania produktów rozszczepienia, opadania grawitacyjnego i meteorologii na wielu wysokościach. Model plumy Gaussa może pomóc zrozumieć część języka dyspersji, ale nie zastępuje modelu chmury powybuchowej.⁴

To rozróżnienie jest ważne także z punktu widzenia bezpieczeństwa treści. Artykuł uczy, jak czytać proste stężenie Bq/m3, depozycję Bq/m2 i dawkę, a nie jak planować realne uwolnienie lub przewidywać skutki operacyjne.

Układ współrzędnych

W modelu przyjmuje się zwykle trzy osie:

• x - odległość z wiatrem od źródła,
• y - odległość w bok od osi plumy,
• z - wysokość nad gruntem.

Jeżeli obserwator stoi na osi plumy, ma y = 0. Jeżeli stoi z boku, czynnik gaussowski exp(-y^2 / (2 sigma_y^2)) obniża stężenie. Dla y = sigma_y stężenie boczne jest już zauważalnie mniejsze, a dla kilku sigma_y spada do małej części wartości osiowej.

W pionie działa podobna logika. sigma_z mówi, jak szybko chmura miesza się w górę i w dół. Małe sigma_z oznacza wąską warstwę; duże sigma_z oznacza intensywne pionowe mieszanie. Przy źródle wyniesionym nad ziemię pionowe rozmycie decyduje, kiedy istotna część chmury dotrze do poziomu gruntu.

Równanie plumy Gaussa

W klasycznym, ciągłym zapisie dla źródła o strumieniu emisji q w Bq/s stężenie na poziomie gruntu można zapisać schematycznie:

C(x,y,0) = q / (2 pi u sigma_y sigma_z) * exp(-y^2 / (2 sigma_y^2)) * F_H.

C ma jednostkę Bq/m3, u jest prędkością wiatru, a F_H jest czynnikiem wysokości uwolnienia i odbicia od gruntu. Wzór pokazuje trzy intuicyjne zależności:

• większy strumień emisji zwiększa stężenie,
• szybszy wiatr i większe sigma rozcieńczają chmurę,
• oddalenie boczne od osi plumy obniża stężenie wykładniczo.

Kalkulator pluma radiologiczna używa tej samej struktury zależności, ale wprowadza aktywność epizodu i czas ekspozycji jako parametry dydaktyczne. Dlatego wyniki należy czytać porównawczo: jako odpowiedź modelu na zmianę klasy stabilności, wiatru, wysokości, opadu i prędkości depozycji, a nie jako gotową ekspertyzę terenową.^1,2

Klasy stabilności Pasquilla-Gifforda

Klasa stabilności jest skrótem opisującym, jak silnie atmosfera miesza zanieczyszczenie. W materiałach dydaktycznych spotyka się klasy od A do F: od bardzo chwiejnej do stabilnej.⁵

Klasy można czytać praktycznie:

• A - atmosfera bardzo chwiejna; silne mieszanie konwekcyjne, szeroka pluma, mniejsze stężenie osiowe w pobliżu ziemi,
• B i C - warunki chwiejne lub lekko chwiejne,
• D - warunki neutralne, często używane jako punkt odniesienia,
• E i F - atmosfera stabilna; słabe mieszanie, wąska pluma, większa wrażliwość wyniku na położenie osi i wysokość uwolnienia.

W kodzie kalkulatora serwisu klasy mają przypisane współczynniki a_y, a_z, c_z i e_z, z których liczone są sigma_y i sigma_z.¹ Nie jest to pełna meteorologia. To tabela korelacji pozwalająca studentowi zobaczyć, że zmiana z klasy A na F nie jest kosmetyką, tylko zmianą szerokości chmury.

Jak liczone są sigma_y i sigma_z

W wersji używanej w kalkulatorze poziome rozmycie ma postać:

sigma_y = a_y x / sqrt(1 + 0,0001 x).

Pionowe rozmycie zaczyna się od:

sigma_z = a_z x,

a dla części klas jest dodatkowo mnożone przez:

(1 + c_z x)^e_z.

Dla klasy A współczynniki są duże, więc sigma_y i sigma_z rosną szybko. Dla klasy F współczynniki są małe, więc chmura pozostaje wąska. W obu przypadkach odległość x ma znaczenie: im dalej od źródła, tym więcej czasu na rozmycie, ale zależność nie musi być liniowa w całym zakresie.

To wyjaśnia typowy paradoks dydaktyczny. Stabilna atmosfera może dać większe stężenie na osi, ale bardzo mały błąd położenia bocznego może wtedy mocno zmienić wynik. Chwiejna atmosfera daje mniejsze maksimum, lecz szerszy obszar niskich stężeń.

Prędkość wiatru

W równaniu Gaussa prędkość wiatru u znajduje się w mianowniku. Jeżeli pozostałe parametry są stałe, większa prędkość wiatru obniża stężenie w danym punkcie, bo ta sama emisja jest przenoszona przez większy przepływ powietrza.

Nie znaczy to, że silny wiatr zawsze "zmniejsza problem". Silniejszy wiatr może szybciej przenieść chmurę dalej, zmienić czas dotarcia do punktu i wpłynąć na depozycję. W prostym modelu jednopunktowym widać jednak najważniejszą zależność: przy stałej geometrii i tej samej klasie stabilności C maleje w przybliżeniu odwrotnie proporcjonalnie do u.

W kalkulatorze warianty niepewności celowo zmieniają prędkość wiatru wraz z uwolnioną aktywnością i depozycją. To nie jest probabilistyczna prognoza pogody, lecz demonstracja, że wynik radiologiczny nie jest liczbą oderwaną od meteorologii.²

Wysokość uwolnienia i odbicie od gruntu

Jeżeli źródło jest przy ziemi, pluma jest od razu blisko poziomu oddychania i powierzchni. W prostym modelu uwzględnia się także odbicie od gruntu przez metodę źródła lustrzanego: chmura nie może dyfundować swobodnie pod powierzchnię ziemi, więc pionowy rozkład jest korygowany tak, jakby istniało drugie, lustrzane źródło.

W kodzie kalkulatora czynnik przy ziemi ma postać:

F_H = 2 * exp(-0,5 * (H / sigma_z)^2).

Dla H = 0 czynnik dąży do 2. Dla wysokiego uwolnienia i małego sigma_z stężenie przy ziemi jest małe, bo chmura nie zdążyła jeszcze wymieszać się pionowo do poziomu gruntu. Dopiero dalej z wiatrem, gdy sigma_z rośnie, wkład przy ziemi może stać się istotny.

To jest jedna z najważniejszych intuicji dla emisji kominowych i wyniesionych. Maksimum przy gruncie nie musi pojawić się tuż przy źródle. Może wystąpić dalej, gdy pionowe rozmycie "sprowadzi" część plumy do poziomu obserwatora.

Depozycja sucha

Depozycja sucha oznacza przechodzenie materiału z powietrza na powierzchnie bez udziału opadu. W prostym zapisie:

D_sucha = C * v_d * t,

gdzie D_sucha jest depozycją w Bq/m2, C stężeniem w powietrzu, v_d efektywną prędkością depozycji, a t czasem ekspozycji.

v_d streszcza wiele procesów: rozmiar cząstek, ich gęstość, chropowatość podłoża, turbulencję przy powierzchni i chemiczną postać radionuklidu. Jedna liczba nie opisuje pełnej fizyki aerozolu. W kalkulatorze jest jednak użyteczna, bo pozwala oddzielić pytanie o stężenie w powietrzu od pytania o osadzanie na gruncie.²

Praktyczny wniosek jest prosty: dwa scenariusze o podobnym Bq/m3 mogą dać różną depozycję Bq/m2, jeśli cząstki mają inną prędkość osiadania lub inny czas kontaktu z powierzchnią.

Depozycja mokra

Depozycja mokra zachodzi, gdy opad usuwa materiał z powietrza i przenosi go na powierzchnię. Deszcz może więc lokalnie zwiększać depozycję, nawet jeśli jednocześnie zmniejsza stężenie w powietrzu.

Pełny model mokrej depozycji wymagałby współczynników wymywania, rozkładu opadu w przestrzeni i czasie, wysokości chmury, rozmiarów kropli oraz rodzaju aerozolu. Kalkulator serwisu stosuje prostszy zamiennik: intensywność deszczu jest przeliczana na efektywną prędkość depozycji mokrej, z ograniczeniem górnym zabezpieczającym przed nierealistycznymi wartościami.^1,2

To przybliżenie dobrze pokazuje kierunek efektu, ale nie powinno być traktowane jak model hydrometeorologiczny. W rzeczywistej ocenie depozycji opad jest jednym z najważniejszych pól pomocniczych. Polskie wytyczne do oceny depozycji atmosferycznej wskazują, że uwzględnienie danych o wysokości opadu może poprawiać mapowanie mokrej depozycji, a do interpolacji rekomendują m.in. kriging zwykły i kokriging z opadem jako zmienną pomocniczą.⁶

Od powietrza do dawki inhalacyjnej

Stężenie w powietrzu nie jest jeszcze dawką. Dla inhalacji potrzebny jest co najmniej czas ekspozycji, tempo oddychania i współczynnik dawki dla danego radionuklidu oraz drogi wniknięcia.

Schematycznie:

I = C * BR * t,

gdzie I jest aktywnością wprowadzoną przez oddychanie, BR jest wentylacją oddechową, a t czasem. Potem:

E = I * e_inh,

gdzie e_inh jest współczynnikiem dawki inhalacyjnej. Współczynniki Sv/Bq lub ich przekształcone wersje pochodzą z modeli dozymetrii wewnętrznej, takich jak kompendia ICRP.⁷

Dlatego w wyniku kalkulatora nie wolno mieszać etapów. Bq/m3 opisuje powietrze. Bq/m2 opisuje powierzchnię. mSv opisuje dawkę po zastosowaniu modelu drogi narażenia. Każdy etap ma osobne założenia i osobne źródła niepewności.

Dawka od zdeponowanego gruntu

Depozycja na gruncie może dawać zewnętrzną moc dawki, szczególnie dla radionuklidów emitujących promieniowanie gamma. Kalkulator zapisuje ten etap jako przeliczenie Bq/m2 na uSv/h przez współczynnik zależny od nuklidu.^1,2

To jest inna ścieżka niż inhalacja. Dla inhalacji ważne jest stężenie w powietrzu w czasie oddychania. Dla gruntu ważne jest to, co zostało zdeponowane, jak długo pozostaje na powierzchni, jaki ma półokres rozpadu, czy jest wymywane, migruje w glebie, przykrywane lub usuwane.

Prosty model jednopunktowy liczy tylko początkowe przejście: pluma -> depozycja -> moc dawki od powierzchni. Długoterminowe narażenie środowiskowe wymaga osobnych modeli, takich jak dawka od depozycji gruntu, food-chain dose albo skażenie środowiskowe.

Dlaczego wynik ma tak dużą niepewność

W modelu plumy wrażliwe są niemal wszystkie wejścia:

• aktywność i frakcja uwolniona do aerozolu,
• klasa stabilności atmosfery,
• prędkość i kierunek wiatru,
• odległość z wiatrem i odchylenie boczne,
• wysokość uwolnienia,
• rozmiar i postać chemiczna cząstek,
• prędkość depozycji suchej,
• opad i mokra depozycja,
• tempo oddychania i czas ekspozycji,
• współczynniki dawki dla radionuklidu.

Niepewność nie jest tu dodatkiem kosmetycznym. Jeżeli wiatr zmieni kierunek, punkt obserwacji może przestać leżeć na osi plumy. Jeżeli atmosfera z neutralnej stanie się stabilna, wzrośnie stężenie osiowe. Jeżeli pojawi się lokalny deszcz, depozycja może zmienić się bardziej niż samo stężenie w powietrzu.

Dlatego kalkulator pokazuje warianty niskiej, centralnej i wysokiej depozycji. To jest lepsze dydaktycznie niż jedna liczba z wieloma cyframi znaczącymi. W prawdziwym modelowaniu środowiskowym potrzebne byłyby pola meteorologiczne, siatka przestrzenna, dane opadowe, walidacja pomiarami i analiza niepewności.

Jak to łączy się z monitoringiem depozycji

Model plumy idzie od źródła do punktu: przyjmujemy źródło, pogodę i geometrię, a potem liczymy stężenie oraz depozycję. Monitoring depozycji często pracuje odwrotnie: mamy pomiary z kilku lub kilkudziesięciu stacji i chcemy zbudować mapę przestrzenną.

Wytyczne do oceny depozycji atmosferycznej w Polsce omawiają właśnie ten drugi problem. Zwracają uwagę na interpolację statystyczną, liczbę punktów pomiarowych, wartości odstające, rozkład opadu i możliwość użycia modeli transportu chemicznego, w szczególności modelu EMEP.⁶

To ważne połączenie dla serwisu. Prosty model Pasquilla-Gifforda uczy fizycznego języka plumy. Kriging, kokriging i modele transportu chemicznego uczą, jak przejść od pomiarów do mapy. Obie warstwy są potrzebne, ale odpowiadają na inne pytania.

Typowe błędy interpretacji

Pierwszy błąd to traktowanie klasy stabilności jako mało ważnej etykiety. W rzeczywistości klasa steruje sigma_y i sigma_z, a więc szerokością chmury i stężeniem osiowym.

Drugi błąd to porównywanie samych aktywności źródła. Dwa uwolnienia o tej samej aktywności mogą dać różne dawki, jeśli jedno zachodzi przy silnym mieszaniu i braku opadu, a drugie przy stabilnej atmosferze i deszczu.

Trzeci błąd to mylenie stężenia z depozycją. Bq/m3 w powietrzu i Bq/m2 na gruncie są powiązane, ale nie są tą samą wielkością. Łączy je czas i prędkość depozycji.

Czwarty błąd to automatyczne przeliczanie Bq/m2 na ryzyko zdrowotne. Dawka zależy od drogi narażenia, czasu, nuklidu, wieku osoby referencyjnej, zachowania radionuklidu w środowisku i przyjętych współczynników.

Co powinien umieć student po tym artykule

Po przejściu przez model student powinien umieć wytłumaczyć, dlaczego stabilna nocna atmosfera może być gorsza dla osiowej plumy niż chwiejny dzień, dlaczego zwiększenie prędkości wiatru obniża stężenie w prostym równaniu, dlaczego wysokość uwolnienia przesuwa maksimum przy gruncie i dlaczego deszcz jest jednocześnie procesem oczyszczania powietrza oraz wzmacniania lokalnej depozycji.

Powinien też umieć oddzielić trzy poziomy rachunku:

• dyspersję w powietrzu,
• depozycję na powierzchni,
• dozymetrię inhalacyjną i zewnętrzną.

Najkrótsze podsumowanie jest takie: Pasquill-Gifford nie jest "prawdą o pogodzie", lecz zwięzłym językiem pierwszego przybliżenia. Daje prostą odpowiedź na pytanie, jak meteorologia zmienia stężenie, depozycję i dawkę w modelu edukacyjnym.

Modele alternatywne: od Gaussa do HYSPLIT

Model Pasquilla-Gifforda jest klasycznym modelem analitycznym dla bliskiego zasięgu (do ~10–20 km) w homogenicznej atmosferze. Dla dalszego zasięgu lub bardziej złożonych warunków stosuje się modele numeryczne:

HYSPLIT (Hybrid Single-Particle Lagrangian Integrated Trajectory): model trajektorii i dyspersji NOAA (USA). Oblicza trajektorie cząstek powietrza lub cząstek aerozolu przez trójwymiarowe pola meteorologiczne z GFS, ECMWF lub innych modeli numerycznych. Stosowany do: analizy źródeł przy wykryciu radionuklidu w stacji IMS, wstecznej rekonstrukcji trasy plumy, prognozowania przy awariach reaktorowych.

NAME (Numerical Atmospheric-dispersion Modelling Environment): model UK Met Office używany przez brityjskie IPEM i NRPB. Podobna filozofia do HYSPLIT, ale z inną parametryzacją turbulencji i różnymi danymi meteorologicznymi (UK Met Office UM model).

ARGOS/RIMPUFF: duński model mesoskalowy do ~100 km. Stosowany przez NATO i kraje skandynawskie. Puff-based (model obłoków, nie pióropusza), co lepiej opisuje epizodyczne uwolnienia i zmienne warunki wiatrowe.

MELCOR i MACCS2 (USA NRC): modele integralne dla oceny ryzyka elektrowni jądrowych (PRA — Probabilistic Risk Assessment). Łączą modele termiczno-hydrauliczne reaktora z modelem dyspersji i modelem dawki dla populacji. Stosowane w licencjonowaniu reaktorów.

Kluczowa różnica od Gaussa:

Cecha	Gauss P-G	HYSPLIT/NAME
Zasięg	do ~10–20 km	do globalnego
Meteorologia	stała, jeden wektor	3D pole w czasie
Terrain	jednorodny	topografia, użytkowanie ziemi
Turbulencja	sigma_y, sigma_z empiryczne	parametryzacja PBL
Zastosowanie	edukacja, ocena wstępna	operacyjne, śledztwa

Wpływ topografii: doliny i budynki

Model Gaussa w standardowej formie nie uwzględnia topografii. W rzeczywistości:

Doliny: lokalizacja w dolinie może prowadzić do "kanalizowania" plumy wzdłuż osi doliny zamiast swobodnego rozmycia w kierunku poprzecznym. Nocna inwersja w dnie doliny może "zamknąć" piumę bardzo blisko gruntu przy bardzo stabilnej atmosferze. Elektrownie zlokalizowane w dolinach górskich (np. typowe lokalizacje szwajcarskie i austriackie) wymagają specjalnych analiz meteorologicznych.

Budynki i "wake" (zawir): wiatr za budynkiem tworzy zawirowanie, które może lokalne zwiększyć stężenie. "Building wake effect" jest istotny dla emisji niskiego źródła (np. komin niższy od budynku) i wymaga modelowania turbulencji CFD lub specjalnych tablic korekcji.

Wzgórza (terrain lifting): chmura wiatrowa, trafiając na wzgórze, częściowo omija przeszkodę, częściowo podnosi się — "Froude number terrain flow" — efekty trudne do opisania prostym Gaussem.

Dla celów dydaktycznych ważne jest rozumienie, że liczba z kalkulatora na jednolitym terenie może nie oddawać rzeczywistości w skomplikowanej topografii. Planowanie awaryjne elektrowni zawsze wymaga site-specific meteorological assessment.

Polskie plany awaryjne i strefy ochrony

Polskie plany awaryjne dla obiektów jądrowych (projekt elektrowni w Choczewie) bazują na wytycznych IAEA i europejskich standardach. Dyrektywa Rady 2013/59/EURATOM i polskie Prawo atomowe definiują:

Strefa planowania awaryjnego: dla elektrowni jądrowych dwie strefy:

• PAZ (Precautionary Action Zone): typowo 3–5 km. Natychmiastowe działania ochronne (schronienie, jodowanie) bez oczekiwania na wyniki pomiarów.
• UPZ (Urgent Protective Action Planning Zone): typowo 10–30 km. Działania ochronne oparte na prognozowaniu i pomiarach.

Obliczenia do wyznaczania stref są oparte na modelach dyspersji dla "referencyjnego zdarzenia projektowego" z założoną aktywnością uwolnienia, meteorologią niekorzystną (klasa stabilności F, prędkość wiatru 2 m/s) i dawką referencyjną (np. 100 mSv efektywna w pierwszych 7 dniach jako kryterium ewakuacji).

Dla Polski projekt elektrowni jądrowej w Choczewie (Morze Bałtyckie, Puck): trwają analizy środowiskowe i meteorologiczne dla oceny zasięgu stref awaryjnych. Dane meteorologiczne z IMGW-PIB są wejściem do modeli dyspersji.

Tabela wartości sigma dla klas stabilności

Poniżej zestawienie typowych wartości sigma_y i sigma_z (m) dla klas stabilności A–F w funkcji odległości od źródła, na podstawie korelacji Pasquilla-Gifforda (przybliżone wartości dla terenów otwartych):

Odległość x [m]	σ_y (A)	σ_y (D)	σ_y (F)	σ_z (A)	σ_z (D)	σ_z (F)
100	22	9	5	20	6	2
500	90	38	20	150	18	5
1000	160	70	36	350	30	8
2000	280	130	65	700	50	12
5000	600	300	145	>1000	100	22
10000	1000	550	265	>2000	170	38

Interpretacja: dla klasy A (bardzo chwiejna, dzień słoneczny, słaby wiatr) sigma_z = 350 m przy x = 1000 m oznacza, że chmura sięga prawie do granicy warstwy mieszania atmosferycznej. Dla klasy F (stabilna, nocna inwersja) sigma_z = 8 m przy tej samej odległości oznacza cienką, wąską wstęgę unoszącą się nad ziemią — maksymalne stężenia osiowe mogą być o 2–3 rzędy wielkości wyższe niż w klasie A.

Ta tabela liczbowo ilustruje, dlaczego klasa stabilności jest najważniejszym parametrem kalkulatora plumy: zmiana z A na F może zmienić wynik o 1000-krotność przy tym samym źródle i geometrii.

Równanie washout: mokra depozycja z koeficjentem wymywania

W uproszczonych modelach mokrą depozycję opisuje się przez koeficjent wymywania Λ (lambda, 1/s):

dC/dt = −Λ × C

czyli stężenie w powietrzu zmniejsza się wykładniczo z powodu opadu. Po czasie przejścia plumy nad punktem obserwacji t_p, frakcja pozostała w powietrzu: exp(−Λ × t_p).

Wartości Λ zależą od intensywności opadu i rozmiaru cząstek:

Opad [mm/h]	Λ (gazy i drobne aerozole) [s⁻¹]	Λ (grubsze cząstki) [s⁻¹]
1	0,5 × 10⁻⁴	1–2 × 10⁻⁴
5	2 × 10⁻⁴	5 × 10⁻⁴
10	4 × 10⁻⁴	1 × 10⁻³
25	8 × 10⁻⁴	2 × 10⁻³

Deponowana aktywność na jednostkę powierzchni w wyniku mokrej depozycji:
D_wet = Λ × ∫C(t) dt × H_cloud

gdzie H_cloud jest wysokością warstwy, przez którą pada deszcz (typowo 1000–3000 m).

Dla celów dydaktycznych: przy intensywnym deszczu (10 mm/h) i czasie przejścia plumy ~30 min (t_p = 1800 s), frakcja wymyta: 1 − exp(−4×10⁻⁴ × 1800) = 1 − exp(−0,72) ≈ 51%. Połowa aktywności z powietrza opada na ziemię.

Systemy wspomagania decyzji: RODOS i ARGOS

W sytuacjach awaryjnych decydenci nie korzystają z ręcznych kalkulatorów Gaussowskich — używają zintegrowanych systemów wspomagania decyzji radiologicznych:

RODOS (Real-time Online Decision Support System): europejski system wspomagania decyzji radiologicznych opracowany przez KIT (Karlsruhe Institut fur Technologie) i Forschungszentrum Jülich. Stosowany przez wiele krajowych organizacji zarządzania kryzysowego, w tym polskie centrum zarządzania kryzysowego. Zawiera: modele dyspersji atmosferycznej (ATSTEP i DERMA), modele depozycji, modele dawki dla populacji, zarządzanie danymi pomiarowymi z sieci stacji, assimilację danych i interfejs map. W Polsce RODOS jest wdrożony w systemie monitoringu PAA.

ARGOS (Accident Reporting and Guiding Operational System): duński system DST (Decision Support Tool), stosowany przez NATO, Demark DCRS i kilkanaście krajów europejskich. Łączy modele dyspersji (RIMPUFF — Risø i MEPS-meteorology) z narzędziami planowania ochrony ludności (strefy ewakuacji, jodowanie jodem stałym).

Różnica od modelu Gaussowskiego:

• Modele w RODOS/ARGOS są "mesoscale" (skala mezometeorologiczna, do 100 km) lub "microscale" (lokalna, teren złożony)
• Pola meteorologiczne pobierane w czasie rzeczywistym z ECMWF, DWD lub IMGW (Polska)
• Możliwość asymilacji danych pomiarowych z sieci DOZ
• Modelowanie depozycji z uwzględnieniem topografii i pokrycia terenu

Model Gaussa jest prostszym, transparentnym narzędziem do edukacji i pierwszego oszacowania. Systemy wspomagania decyzji są operacyjnymi narzędziami dla profesjonalistów.

Trzy przykłady obliczeniowe

Przykład 1: Stężenie w punkcie na osi plumy

Dane: strumień emisji q = 10⁸ Bq/s (Cs-137). Prędkość wiatru u = 4 m/s. Odległość x = 2000 m. Klasa stabilności D (neutralna). Uwolnienie przy ziemi (H = 0).

Z tabeli: σ_y(D, 2000 m) = 130 m, σ_z(D, 2000 m) = 50 m.

C(2000, 0, 0) = q / (π × u × σ_y × σ_z) × exp(0) × F_H

Dla H = 0: F_H = 2

C = 10⁸ / (3,1416 × 4 × 130 × 50) × 2 = 10⁸ / 81,68 × 10³ × 2 = 10⁸ / 81681 × 2 = 2449 Bq/m³.

Dla H = 30 m (komin): F_H = 2 × exp(−0,5 × (30/50)²) = 2 × exp(−0,18) = 2 × 0,835 = 1,670.
C = 10⁸ / (3,14 × 4 × 130 × 50) × 1,670 = 10⁸ / 81681 × 1,670 = 2042 Bq/m³.

Wynik: uwolnienie naziemne daje ~2450 Bq/m³, a kominem z H = 30 m — ~2040 Bq/m³ na osi plumy w 2000 m. Dla tej odległości różnica jest niewielka. Dla mniejszych odległości (x = 200 m) komin daje dużo mniejsze stężenie przy ziemi.

Przykład 2: Dawka inhalacyjna z 1 godz. ekspozycji

Stężenie Cs-137 w powietrzu: C = 500 Bq/m³ (przez 1 godzinę). Wentylacja oddechowa dorosłego: BR = 1,2 m³/h.

Aktywność wdychana: I = C × BR × t = 500 × 1,2 × 1 = 600 Bq.

Współczynnik dawki inhalacyjnej ICRP 119 dla Cs-137, dorosły: e_inh = 8,6 × 10⁻⁹ Sv/Bq.

Dawka efektywna: E = I × e_inh = 600 × 8,6×10⁻⁹ = 5,16 × 10⁻⁶ Sv = 5,2 µSv.

Porównanie: 5 µSv jest porównywalne z dawką od jednego prześwietlenia rentgenowskiego klatki piersiowej (~20 µSv) lub lotem samolotem Warszawa-Londyn (~20 µSv).

Przykład 3: Moc dawki od skażonego gruntu

Depozycja Cs-137 na gruncie: D = 10 000 Bq/m².
Współczynnik przeliczenia Cs-137 depozycja na moc dawki gamma (dla geometrii nieskończonej płaszczyzny, na wysokości 1 m): ~1,5 × 10⁻⁷ mSv/h per Bq/m².

Moc dawki gamma od gruntu: H_rate = 10000 × 1,5×10⁻⁷ = 1,5 × 10⁻³ mSv/h = 1,5 µSv/h.

Kontekst: normalne tło naturalne w Polsce wynosi 0,07–0,15 µSv/h. Depozycja 10 000 Bq/m² Cs-137 (typowa dla strefy 30 km od Czarnobyla) daje 10-krotność naturalnego tła. Dla depozycji 1 mln Bq/m² (strefa wykluczona, <10 km od Czarnobyla) moc dawki od gruntu: ~150 µSv/h = 1,3 mSv/rok (tło) + 1300 mSv/rok od gruntu — co uzasadnia ewakuację.

Pytania otwarte

1. Dlaczego model Gaussa zakłada stały wiatr i jednorodną atmosferę, i jak te założenia są naruszane przez rzeczywiste zjawiska meteorologiczne (fronty, inwersje termiczne, topografia)? W jakich warunkach model Gaussa jest najgorzej przybliżony do rzeczywistości?

2. Jak zmienia się depozycja mokra na jednostkę powierzchni, gdy plumy przelatują nad obszarami o różnej intensywności opadu? Dlaczego mapy depozycji Cs-137 po Czarnobylu mają nieregularny, "mozaikowy" kształt zamiast regularnych elips zgodnych z kierunkiem wiatru?

3. W jakich sytuacjach model Pasquilla-Gifforda jest zdecydowanie zbyt uproszczony i wymaga zastąpienia modelami mesoskalowymi (np. RODOS, HYSPLIT)? Jakie minimalne dane meteorologiczne są potrzebne do uruchomienia modelu mesoskalowego?

4. Jak zachowanie chemiczne różnych radionuklidów wpływa na prędkość depozycji suchej v_d? Dlaczego cesowy v_d jest inny niż jodowy, a gazowe formy jodu (I₂, CH₃I) mają v_d bliskie zeru nawet na trawę?

5. Jak definiuje się "strefę planowania awaryjnego" wokół elektrowni jądrowej i jakie obliczenia modelu plumy są podstawą tej definicji? Jakie aktywności uwolnienia i klasy stabilności przyjmuje się jako "planistyczny scenariusz referencyjny"?

6. Jak koincydencja stabilnej atmosfery nocnej z kierunkiem wiatru w stronę zaludnionego obszaru wpływa na planowanie systemów monitoringu i wczesnego ostrzegania wokół instalacji nuklearnych? Jakie "najgorsze warunki meteorologiczne" należy uwzględnić w planie?

7. Dlaczego przy ocenie depozycji na obszarach rolniczych i w konsekwencji dawki przez żywność sam model plumy jest niewystarczający? Jakie modele "food chain" (łańcuch żywnościowy) są stosowane i jakie dodatkowe parametry wymagają?

8. Jak zmienia się wynik kalkulatora plumy, gdy uwolnienie nie jest punktowe i ciągłe, ale złożone z serii krótkich epizodów (puff model)? W jakich sytuacjach model puff jest bardziej odpowiedni niż model continuous plume?

Podsumowanie dydaktyczne

1. Model Gaussa to przybliżenie, nie prognoza — równanie Pasquilla-Gifforda zakłada stałą meteorologię, jednorodny teren, punktowe źródło i izotropową turbulencję. Żadne z tych założeń nie jest spełnione idealnie w terenie. Model jest jednak niezwykle użyteczny do zrozumienia kierunku efektów i gruboszacunkowych obliczeń.

2. Klasa stabilności steruje rozrzutem o kilka rzędów wielkości — zmiana z klasy A na F może zmienić stężenie osiowe o 100–1000 razy. To jest największy pojedynczy parametr niepewności w modelowaniu plumy, ważniejszy niż niepewność samej aktywności uwolnienia w wielu scenariuszach awaryjnych.

3. Sigma rośnie z odległością, ale inaczej dla y i z — sigma_y (boczna) rośnie prawie liniowo z odległością. Sigma_z (pionowa) rośnie szybciej w klasach chwiejnych, ale może być ograniczona przez inw ersję temperaturową. To wyjaśnia, dlaczego "pułap inersji" jest ważnym parametrem dla dalekiego transportu.

4. Prędkość wiatru rozcieńcza, ale nie tylko — zwiększony wiatr rozcieńcza chmurę przez większy przepływ powietrza (mianownik 1/u), ale też szybciej przenosi ją na dalszy dystans. Dla krótkotrwałych uwolnień (puff) efekt czasu ekspozycji może być ważniejszy niż rozcieńczenie.

5. Depozycja sucha i mokra to różne procesy — sucha depozycja zależy od prędkości osiadania cząstek, turbulencji przygruntowej i chropowatości terenu. Mokra depozycja zależy od intensywności opadu, rozmiarów kropelek i koeficientu wymywania. Mogą zachodzić jednocześnie, ale z różnymi wydajnościami dla różnych form chemicznych radionuklidu.

6. Od stężenia w powietrzu do dawki potrzeba łańcucha modeli — pluma → stężenie w powietrzu → dawka inhalacyjna (model wentylacji i ICRP e(g)) → depozycja na gruncie → moc dawki gamma → skażenie roślin → dawka przez spożycie. Każde ogniwo ma własne niepewności i zależy od dodatkowych parametrów.

7. RODOS i ARGOS to narzędzia operacyjne, Gauss to narzędzie edukacyjne — systemy wspomagania decyzji używane przez służby obrony cywilnej integrują meteorologię rzeczywistą, modele dyspersji i dane z sieci monitoringu. Model Gaussa pozwala zrozumieć, jakie pytania te systemy starają się odpowiedzieć.

8. "Jaka dawka?" to pytanie wieloetapowe — odpowiedź wymaga: aktywności uwolnienia, meteorologii, odległości i geometrii, drogi narażenia (inhalacja/spożycie/promieniowanie zewnętrzne), czasu ekspozycji, nuklidu i współczynnika dawki, wieku i charakterystyk osoby referencyjnej. Pominięcie któregokolwiek z tych elementów daje odpowiedź niekompletną lub błędną.

Specyfika radionuklidów w modelu plumy: I-131, Cs-137, Kr-85

Nie wszystkie radionuklidy zachowują się tak samo w piumie:

Jod-131 (T₁/₂ = 8,02 dnia):

• Krótki okres połowicznego zaniku — po miesiącu transportu atmosferycznego zanika o czynnik 16.
• Dwie formy chemiczne: jod gazowy I₂ (v_d = 0,01–0,02 m/s na trawę) i metylojod CH₃I (v_d = 0,001–0,002 m/s). Gazowy jod łatwo deponuje się na pastwiskach — stąd "droga mleczna" (pastwisko → krowa → mleko → tarczyca niemowlęcia).
• Kluczowe dla dawki: nie stężenie w powietrzu, ale akumulacja w tarczycy przez mleko. Model plumy musi być uzupełniony modelem łańcucha pokarmowego.

Cez-137 (T₁/₂ = 30,1 roku):

• Długi okres zaniku — praktycznie nie zanika w transporcie atmosferycznym.
• Forma chemiczna: cząstki stałe (aerozol Cs⁺ zaadsorbowany na cząstkach pyłu). Prędkość suchej depozycji zależy od rozkładu rozmiarów cząstek (typowo v_d = 0,001–0,01 m/s).
• Główna droga dawki po awarii: zewnętrzne promieniowanie gamma od skażonego gruntu (długoterminowe).

Krypton-85 (T₁/₂ = 10,76 roku):

• Gaz szlachetny — nie deponuje się na gruncie (v_d ≈ 0). Nie uczestniczy w łańcuchu pokarmowym.
• Rozprasza się globalnie. Jego obecność w powietrzu jest markerem przemysłu jądrowego (przede wszystkim przeróbki paliwa — PUREX).
• Dla modelu plumy: praktycznie bez depozycji → stężenie w powietrzu jest jedyną ścieżką narażenia → dawka przez promieniowanie beta zewnętrzne.

Tabela porównawcza radionuklidów w kontekście modelu plumy:

Nuklid	T₁/₂	Forma w piumie	v_d	Główna droga dawki
I-131	8d	Gaz+cząstki	0,005–0,02 m/s	Tarczyca przez mleko
Cs-137	30a	Aerozol	0,001–0,01 m/s	Zewn. gamma od gruntu
Kr-85	11a	Gaz szlachetny	~0 m/s	Beta zewnętrzne
Pu-239	24ka	Aerozol	0,001–0,005 m/s	Płuca przez inhalację
Sr-90	29a	Aerozol	0,001–0,01 m/s	Kości przez spożycie

Ta tabela pokazuje, że "jeden model plumy" dla wielu radionuklidów wymaga różnych ścieżek dozy — każdy nuklid ma własny "critical pathway".

Walidacja modeli plumy: dane z Czarnobyla i Fukushimy

Model plumy jest matematyczną konstrukcją, ale jego wiarygodność można testować przez porównanie z rzeczywistymi danymi:

Czarnobyl 1986: CTBTO, IAEA i instytucje krajowe zebrały dane z setek stacji po awarii. Model plumy (w wersji HYSPLIT i ARGOS) był używany retrospektywnie do rekonstrukcji transportu atmosferycznego z Czarnobyla na Europę w dniach 26 kwietnia – 10 maja 1986. Wyniki pokazują: model dobrze odtwarza ogólną trajektorię i czas dotarcia plumy, ale słabo odtwarza lokalną depozycję Cs-137 — zróżnicowaną przez silnie zmienne opady deszczu. Bez danych o opadach z meteoradaru niemożliwa jest precyzyjna rekonstrukcja depozycji.

Fukushima 2011: IAEA, CTBTO i JMA (Japan Meteorological Agency) uruchomiły modele w czasie rzeczywistym w ciągu pierwszych godzin. CTBTO używał HYSPLIT i FLEXPART do prognozowania, gdzie radionuklidy z Fukushimy mogą dotrzeć do stacji IMS. Stacje w Japonii, Chinach, Rosji i USA wykryły Cs-137 i I-131 w kolejności i czasie dobrze zgodnej z prognozą modeli. Stanowi to jeden z największych testów walidacyjnych modeli dyspersji w historii.

Ru-106 2017 (Majak?): Nieidentyfikowane uwolnienie Ru-106 wykryte w Europie w październiku 2017 r. Modele ATM (wstecznie z HYSPLIT) wskazały rejon południowego Uralu jako możliwy obszar źródłowy — spójny z lokalizacją zakładu Majak. Oficjalnego potwierdzenia nie było (Rosja zaprzeczyła), ale przypadek jest ważnym przykładem jak ATM może prowadzić do atrybucji źródła.

Te przykłady pokazują, że modele plumy nie są tylko akademickimi narzędziami — są aktywnie używane w sytuacjach awaryjnych i dochodzeniach środowiskowych, a ich walidacja przez dane rzeczywiste jest ciągłym procesem. Dla polskiej sieci monitoringu DOZ, stacji IS26 IMS (Stefanów) i laboratoriów CLOR, rozumienie modeli transportu atmosferycznego jest konieczne do właściwej interpretacji danych: kiedy wzrost aktywności jest lokalny, kiedy wynika z transportu dalekiego zasięgu, i w jakim stopniu brak sygnału w stacji może wykluczać uwolnienie w danym regionie.

Dodatkowe materiały multimedialne

Ćwiczenia praktyczne

Pierwsze ćwiczenie: dla klasy D, x = 1000 m, a_y = 0,08, a_z = 0,06, c_z = 0,0015, e_z = -0,5 obliczyć sigma_y i sigma_z z uproszczonych wzorów używanych w kalkulatorze. Następnie powtórzyć rachunek dla klasy A i porównać szerokość plumy.

Drugie ćwiczenie: przyjmując ten sam punkt na osi plumy, porównać wpływ wiatru u = 2 m/s i u = 6 m/s na stężenie w równaniu Gaussa. Nie trzeba znać aktywności źródła; wystarczy wskazać, ile razy zmienia się czynnik 1/u.

Trzecie ćwiczenie: dla H = 0 m, H = 30 m i sigma_z = 15 m obliczyć czynnik F_H = 2 * exp(-0,5 * (H / sigma_z)^2). Wyjaśnić, dlaczego wysokie uwolnienie może zmniejszać stężenie przy gruncie blisko źródła.

Czwarte ćwiczenie: dla stężenia C = 100 Bq/m3, prędkości depozycji v_d = 0,002 m/s i czasu t = 1800 s obliczyć depozycję suchą w Bq/m2. Następnie dodać efektywną prędkość mokrej depozycji 0,001 m/s i porównać wynik.

Piąte ćwiczenie: dla C = 50 Bq/m3, wentylacji 1,2 m3/h i czasu 30 min obliczyć wdychaną aktywność. Następnie wskazać, jakie dodatkowe dane są potrzebne, aby zamienić ją na dawkę.

Szóste ćwiczenie: wejść do kalkulatora pluma radiologiczna i zmienić tylko klasę stabilności z A na F. Zanotować zmianę sigma_y, sigma_z, stężenia i depozycji, a potem opisać, która wielkość zmieniła się najbardziej.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły

• Depozycja atmosferyczna radionuklidów: interpolacja, modele transportu i monitoring
• Wielkości dozymetryczne: dawka pochłonięta, równoważna, efektywna i operacyjna
• ICRP 119 i współczynniki dawki: skąd biorą się liczby w kalkulatorach
• Opad promieniotwórczy: fallout i mechanizm powstawania
• Pluma skażenia radiologicznego