Vela 1979 jako problem metrologii, a nie tylko historii

Streszczenie

Incydent Vela z 22 września 1979 jest zwykle opowiadany jako zagadka polityczna: czy nad południowym Atlantykiem doszło do tajnego testu jądrowego, a jeśli tak, kto go przeprowadził. Dla metrologii ciekawsze jest inne pytanie: co zrobić z sygnałem bardzo podobnym do znanej sygnatury jądrowej, jeśli nie zgadzają się wszystkie kanały pomocnicze i część danych pozostaje niejawna.^1,2

Ten artykuł traktuje Vela 1979 jako studium niejednoznacznego alarmu. Nie rozstrzyga, czy był to test jądrowy; pokazuje raczej, jak porównywać hipotezy, ważyć kanał optyczny, oceniać brak potwierdzeń, rozumieć fałszywe alarmy i odróżniać "nieudowodnione" od "obalone".

Rozszerzenie tematu

Statystyka testowania hipotez w kontekście monitoringu jądrowego

Metrologiczny problem Vela można sformalizować przez pryzmat teorii testowania hipotez Neymana-Pearsona. Mamy dwie hipotezy:

• H₀ (null): zdarzenie jest artefaktem (meteoroid, refleksja, błąd instrumentu)
• H₁ (alternative): zdarzenie jest detonacją jądrową

Błąd I rodzaju (α): Odrzucamy H₀ gdy jest prawdziwa → fałszywy alarm (false positive). System zaalarmuje rząd, gdy nie ma wybuchu. Konsekwencje: dyplomatyczne napięcia, niepotrzebne działania.

Błąd II rodzaju (β): Nie odrzucamy H₀ gdy jest fałszywa → przeoczenie (false negative/miss). System nie wykryje prawdziwego wybuchu. Konsekwencje: proliferacja bez konsekwencji.

Strategia projektowania systemu monitoringu: Przy incydencie Vela panel prezydencki wybrał niskie α (konserwatywnie nie alarmował rządów i nie wydał publicznego werdyktu pro-jądrowego) kosztem wyższego β. To jest racjonalne, gdy:

• skutki fałszywego alarmu są wysokie politycznie
• środki detekcji są niepełne
• brak redundancji potwierdzającej

ROC (Receiver Operating Characteristic): Dla każdego progu sygnału bhangmetru, który uznajemy za „wybuchy", mamy parę punktów (wskaźnik fałszywych alarmów α, moc testu 1−β). Krzywa ROC dla systemu z 41/41 poprawnych detekcji wcześniejszych sugeruje wysoką czułość i swoistość. Ale jeden punkt na historycznej krzywej (Vela 1979) z odchyleniem od typowego stosunku kanałów może oznaczać:

• błąd (poruszamy się wzdłuż krzywej ROC)
• inną krzywą ROC (nowe zjawisko fizyczne lub uszkodzenie instrumentu)

Próg decyzji a informacja a priori: W klasycznym podejściu Neymana-Pearsona próg sygnału alarm/nie-alarm jest wyznaczany dla ustalonego α — bez uwzględniania prawdopodobieństwa a priori. W Bayesowskim — próg zależy od kosztu błędów i prawdopodobieństwa a priori. Vela jest przypadkiem, gdzie różne grupy analityczne używały różnych priorów (wojsko vs cywilni naukowcy) i różnych funkcji kosztów (skutki dyplomatyczne vs skutki nieproliferacyjne).

System IMS CTBTO — architektura i możliwości

Incydent Vela pokazał, że sieć monitoringu 1979 roku była niewystarczająca. Dało to impuls do budowy IMS (International Monitoring System) w ramach CTBT:

Cztery technologie IMS:

1. Sejsmika: 170 stacji (50 głównych, 120 pomocniczych). Wyznacza lokalizację zdarzenia z dokładnością ~10 km i magnitudo do 0,5 M_b (odpowiada ~0,1 kt dla testu podziemnego).
2. Infradźwięki: 60 stacji, wykrywają eksplozje atmosferyczne o mocy ≥ 1 kt w odległości < 5000 km. Kluczowe dla wybuchu Vela-type.
3. Hydroakustyka: 11 stacji oceanicznych (w tym stacje na Azorach, Crozet, Diego Garcia, Iwo Jima), wykrywają wybuchy podwodne i nadbrzeżne do setek km od pobrzeża.
4. Radionuklidy: 80 stacji aerozoli (wysokorozdzielczy HPGe) + 40 stacji ksenonu (SAUNA/ARIX). Pokrycie atmosferyczne globalne, wrażliwość na uwolnienia z testów i incydentów reaktorowych.

IDC (International Data Centre, Wiedeń): Centrum obliczeniowe CTBTO przetwarza ~15 GB danych/dobę, produkuje automatyczne listy zdarzeń (SEL — Standard Event List), które analitycy manualnie weryfikują. Zdarzenia klasyfikowane jako "NDC" (National Data Centre) — każde państwo strona CTBT może pobierać pełne dane do własnej analizy.

Próg detekcji IMS vs 1979:

• Sejsmika 1979: próg ~4,5 M_b (odpowiada ~5 kt podziemny), pokrycie regionalne
• Sejsmika IMS 2024: próg ~3,5 M_b w pokryciu globalnym (poza biegunami)
• Infradźwięki 1979: brak globalnej sieci
• Infradźwięki IMS 2024: 60 stacji, wykrycie wybuchu Vela-type byłoby niemal pewne przy mocy > 2 kt

Symulacja IMS dla Vela 1979: Gdyby 22 września 1979 zdarzenie powtórzyło się z siecią IMS 2024:

• IS-02 (Juan Fernández), IS-07 (Tristan da Cunha, Atlantyk Południowy): oba w zasięgu 5 000 km, wykrycie 90% pewne dla ≥1 kt
• RN-53 (Argentyna), RN-49 (Nowa Zelandia): próbki radionuklidów po 5–8 dniach transportu
• Ksenon SAUNA: wykrycie Xe-133 przy stężeniu > 0,1 mBq/m³ po 3–7 dniach

Wniosek: Nowoczesna sieć IMS z 90% prawdopodobieństwem rozstrzygnęłaby incydent Vela w ciągu 7–14 dni od zdarzenia. Brak sygnału infradźwiękowego byłby sam w sobie bardzo silnym argumentem przeciw hipotezie atmosferycznego wybuchu jądrowego.

Leczka historyczna: Żaden z publicznych raportów z 1979 roku nie zawierał systematycznej analizy możliwości detekcji przy ówczesnej sieci (false negative rate). Eksperci dyskutowali o tym, co „powinno być widoczne", nie o tym, jakie było rzeczywiste prawdopodobieństwo detekcji dla danej mocy wybuchu i lokalizacji. To fundamentalny błąd metrologiczny — ocena negatywnego wyniku bez znajomości czułości systemu detekcji jest niekompletna. W dzisiejszych analizach IDC CTBTO oblicza detection probability explicite dla każdego scenariusza przed wydaniem orzeczenia.

Technika bhangmetrów satelitarnych — szczegóły metrologiczne

Bhangmetry (bang-meters, oficjalnie: Vela Uniform Detection Sensors, BUDS) na satelitach Vela były instrumentami optycznymi projektowanymi do rejestracji charakterystycznego przebiegu jasności atmosferycznej detonacji jądrowej. Oto ich metrologiczne szczegóły:

Dwa kanały detekcji: Każdy bhangmetr zawierał dwa fotodetektory o różnych zakresach dynamicznych i czasach odpowiedzi:

• Kanał szybki (fast channel): odpowiedź czasowa ~1 ms; rejestruje pierwszy krótki impuls (pierwszą fazę jaśnienia kuli ognistej)
• Kanał wolny (slow channel): odpowiedź ~10 ms; rejestruje drugi, dłuższy i słabszy impuls

Oczekiwany przebieg jądrowy: pierwszy impuls o czasie wzrostu ~1 ms, następnie spadek do minimum (breakaway — faza gdy fronton uderzeniowy wyprzedza kule ognistą i przesłania ją), potem wolniejsze narastanie drugiego impulsu trwające setki milisekund.

Stosunek amplitud: Dla znanych detonacji jądrowych istnieje charakterystyczny stosunek A₂/A₁ (amplituda drugiego impulsu do pierwszego). Ten stosunek zależy od mocy wybuchu, wysokości, warunków atmosferycznych i geometrii obserwacji, ale mieści się w charakterystycznym zakresie kalibracyjnym.

Kalibracja historyczna: Do 1979 roku Vela 6911 (wystrzelony 1970, lata na orbicie geosynchronicznej ~90 000 km) zarejestrował 41 wcześniejszych podwójnych błysków — wszystkie z późniejszym potwierdzeniem detonacji jądrowych (testy atmosferyczne w latach 1963–1979, głównie Francja i Chiny). Ten historyczny rekord był podstawą zaufania do instrumentu.

Ograniczenia metrologiczne bhangmetrów:

• Brak pozycjonowania przestrzennego (nieobrazujące — tylko natężenie, nie kierunek)
• Geometria obserwacji zależy od pozycji Słońca (reflektancja od chmur, oceanu)
• Starzenie fotomultiplikatorów po 9 latach w przestrzeni kosmicznej — zmiana czułości i proporcji kanałów
• Brak możliwości weryfikacji in situ kalibracji po wystrzeleniu

Anomalia kanałów w 1979: Zmierzony stosunek amplitud odbiegał od zakresu kalibracji historycznej. To było źródłem kontrowersji — czy to wskazywało na artefakt, czy na „niestandardowy" wybuch (np. inną wysokość, konfigurację chmur, mniejsza moc, inną geometrię)?

Atmosferyczna detekcja radionuklidów — co można było wykryć po wybuchu

Po powierzchniowym lub atmosferycznym wybuchu jądrowym o mocy rzędu 1–10 kt (typowa moc dla niedużego testu) uwalniane są produkty rozszczepienia i aktywacji. Prognozowanie tego, co byłoby możliwe do wykrycia w 1979 roku:

Produkty rozszczepienia kluczowe w krótkim czasie (<2 tygodnie):

• I-131 (T₁/₂ = 8,02 d), I-132 (T₁/₂ = 2,28 h), I-133 (T₁/₂ = 20,8 h): lotne izotopy jodu, łatwo detekcja filtrem węglowym
• Ba-140/La-140 (T₁/₂ 12,75 d / 1,68 d): silne emitenty gamma
• Zr-95/Nb-95 (T₁/₂ 64 d / 35 d): długo mierzalne, charakterystyczne linie 756/766 keV
• Ce-141 (T₁/₂ 32,5 d): wskaźnik fissji
• Sr-89 (T₁/₂ 50,5 d): emiter beta, ważny dla bilansu strontu
• Sr-90/Y-90 (T₁/₂ 28,8 lat): długotrwały depozyt, trudny do przypisania do konkretnego zdarzenia

Charakterystyczny stosunek izotopowy (signature):
Stosunek Zr-95/Nb-95, Ce-141/Ce-144, Ba-140/Ce-141 zależy od czasu po wybuchu (aktywność Nb-95 narasta z Zr-95). Zmierzony stosunek w próbce pozwala wstecznie obliczyć czas od detonacji (age dating opadu) — pod warunkiem pobrania próbki.

Ksenon radioaktywny (Xe-131m, Xe-133, Xe-133m, Xe-135):
Ksenon radioaktywny jest lotny, nie adsorbuje na ziemi, przemieszcza się jako gaz szlachetny w troposferze. CTBTO IMS używa detektorów SAUNA (Swedish Automatic Unit for Noble gas Acquisition) do pomiaru izotopów ksenonu — to najcenniejszy wskaźnik testu jądrowego w atmosferze.
W 1979 roku żaden odpowiednik SAUNA nie istniał. Sieć była słaba i skupiona na produktach partykularnych (filtry do aerozoli).

Transport atmosferyczny:
Produkty rozszczepienia w 1 kt wybuchu atmosferycznym nad południowym Atlantykiem (ok. 47°S) rozpraszają się przez zachodniopacyficzne strumienie powietrza. Czas dotarcia do stacji zlokalizowanych w USA, Europie lub Australii: 5–15 dni w zależności od sezonu i wzorców cyrkulacji.

Modelowanie wsteczne trajektory (backward trajectory) z potencjalnych lokalizacji testu wskazuje, że próby zbierania danych musiałyby uwzględniać kilkudniowe okno czasowe i właściwe regiony geograficzne. Program lotów AFTAC (US Air Force Technical Applications Center) wykonywał właśnie takie zbiórki.

Sejsmika, infradźwięki i hydroakustyka wybuchu atmosferycznego

Sejsmika podziemna vs atmosferyczna: Wybuchy podziemne generują sejsmiczne fale P i S o charakterystycznym stosunku p/s, kierunku i głębokości. Wybuchy atmosferyczne sprzęgają energię do gruntu przez ciśnienie akustyczne — to tzw. airburst coupling. Sprzężenie jest 100–10 000 razy słabsze niż bezpośrednio podziemne.

Dla 5 kt wybuchu atmosferycznego w odległości 5000 km: magnitudzie ~3,5–4,0 M_b (zamiast ~5,5–6,0 dla podziemnego). To poniżej progu stałego monitoringu w 1979 roku.

Infradźwięki: Eksplozje atmosferyczne generują fale infradźwiękowe (< 20 Hz), które rozchodzą się w kanale dźwiękowym troposfery (SOFAR channel). Zasięg detekcji dla 1–10 kt: tysiące km. IMS CTBTO posiada 60 stacji infradźwiękowych — żadna nie istniała w 1979 roku. Pojedyncze stacje naukowe mogły zarejestrować słaby sygnał.

Hydroakustyka: Wybuchy nadwodne lub w płytkim oceanie generują fale hydroakustyczne T. Globalna sieć SOSUS (US Navy) wykrywała okręty podwodne w całym oceanie — mogła też wykryć wybuchy nad/w oceanie. Dane SOSUS z 1979 są w dużej mierze wciąż niejawne.

Analiza porównawcza:

Kanał	Wrażliwość na wybuch atm.	Status w 1979	Status w 2024 (IMS)
Bhangmetr satelitarny	Wysoka (bezpośrednia)	Aktywny (Vela)	GPS/SBIRS (niejawne)
Sejsmika	Niska (pośrednia)	Ograniczona	170 stacji IMS
Infradźwięki	Wysoka (bezpośrednia)	Brak stacji	60 stacji IMS
Hydroakustyka	Wysoka (nad oceanem)	Niejawna SOSUS	11 stacji IMS
Radionuklidy (aerozole)	Wysoka (ale opóźniona)	Słaba sieć	80 stacji IMS
Radionuklidy (ksenon)	Bardzo wysoka	Brak techniki	40 stacji SAUNA/IMS

Rekonstrukcja scenariuszy interpretacyjnych

Metrologicznie poprawne jest rozważenie każdej hipotezy z pytaniem: „jakie dane byłyby spójne z tą hipotezą i jakie dane ją falsyfikują?"

Scenariusz A (Wybuch jądrowy, mała moc ~2–5 kt, wysokość ~0–3 km):

• Spójne z: podwójnym błyskiem, brakiem silnej sejsmiki (mała moc + atmosfera), słabym opadem rozcieńczonym przez ocean
• Problemowe: stosunek kanałów bhangmetru, brak pewnych potwierdzeń w danych AFTAC
• Falsyfikuje: gdyby satelita meteorologiczny w tym samym czasie pokazał jasny błysk z tym samym CG koordynatem

Scenariusz B (Meteoroid uderzający w satelitę):

• Spójne z: brakiem radionuklidów, brakiem sejsmiki, brakiem hydroakustyki
• Problemowe: odtworzenie dokładnej krzywej czasowej (szybki wzrost, minimum, wolny wzrost) ze zwykłego odbicia

Eksperci podkreślali: aby meteoroid wytworzył podwójny błysk o kształcie podobnym do jądrowego, musiałby spełnić bardzo specyficzne warunki geometryczne — mało prawdopodobne, ale nie niemożliwe.

Scenariusz C (Wybuch niejądrowy — konwencjonalny, duży):

• Moc wybuchu konwencjonalnego zdolna wytworzyć podobny sygnał: > 100 kt ekwiwalentu TNT — niezgodne z militarnymi możliwościami 1979 dla testu tajnego
• Wykluczone przez rozmiar dostępnych ładunków konwencjonalnych

Polska perspektywa — monitoring i CTBT

Polska ratyfikowała CTBT (Comprehensive Nuclear-Test-Ban Treaty) w 1999 roku. Na terytorium Polski w siatce IMS działają:

• Stacja sejsmiczna: PS-26 Suwałki (stacja główna IMS dla Polski)
• Stacja infradźwiękowa: IS-26 Warszawa (oddana ok. 2010)
• Stacja hydroakustyczna: brak na terytorium Polski
• Stacja radionuklidów: RN-56 Suwałki

PAA koordynuje przekazywanie danych z polskich stacji IMS do Centrum Danych CTBTO w Wiedniu (IDC — International Data Centre). Dane są analizowane automatycznie i manualnie przez analityków IDC.

W kontekście Vela 1979: gdyby to samo zdarzenie miało miejsce dziś nad Antarktydą lub południowym Atlantykiem, polskie stacje sejsmiczne (PS-26) i infradźwiękowe (IS-26) byłyby nieczułe na ten region (zbyt daleko), ale globalna sieć IMS rozlokowana na stacjach antarktycznych (Palmer, Mawson, Syowa) i południowoatlantyckich mogłaby dostarczyć rozstrzygające dane.

Trzy przykłady rachunkowe

Przykład 1 — Estymacja zasięgu detekcji infradźwiękowej (2024 IMS vs 1979)

Stacja infradźwiękowa IS-26 Warszawa (53°N, 21°E). Potencjalny wybuch Vela: ~47°S, 40°E.
Odległość: Δφ = 100°, Δλ = 19°. Odległość kołowa ≈ R_Ziemi × (Δkąt w radianach).
d ≈ 6371 km × √((100° × π/180)² + (19° × π/180)²) = 6371 × √(3,046 + 0,110) = 6371 × 1,776 ≈ 11 313 km.

Prędkość propagacji infradźwięków w troposferze: v ≈ 340 m/s (nominalnie), ale z kanałem troposferycznym v_eff ≈ 310–350 m/s.

Czas dotarcia: t = 11 313 km / 0,340 km/s ≈ 33 273 s ≈ 9,24 godzin.

Amplituda sygnału przy stacji: dla 5 kt wybuchu, amplituda spada ~1/d. Przy d = 11 313 km i progu detekcji stacji IS-26 ≈ 0,1 Pa, niezbędna amplitude przy źródle: ~1000 Pa. Dla 5 kt wybuchu: typ. amplituda ~10 Pa w odległości 1000 km → 10 × (1000/11313)^1.5 ≈ 0,03 Pa → poniżej progu IS-26. Stacja polonijna byłaby czuła na ten wybuch tylko przy mocy > ~50 kt.

Wniosek: stacja infradźwiękowa w Polsce nie byłaby w stanie samodzielnie wykryć małego (<10 kt) wybuchu nad południowym Atlantykiem. Konieczna jest stacja Antarctica lub Południowoatlantycka (np. IS-02 Juan Fernández, IS-27 Antarktyka).

Przykład 2 — Transport atmosferyczny izotopów z Antarktydy

Wybuch 5 kt w lokalizacji 47°S, 40°E. Dominujący wiatr w tej szerokości: Wiatry westerlies, 15–25 m/s ze wschodu na zachód (a więc niosące z E na W), prędkość horyzontalnej advection ~20 m/s.

Czas dotarcia do stacji w Australii Zachodniej (~120°E, 34°S): dystans ~7000 km, czas advection = 7000 km / (20 m/s × 86,4) ≈ 4,05 dnia.

W tym czasie Ba-140 (T₁/₂ = 12,75 d): f = exp(−ln2/12,75 × 4,05) = exp(−0,220) = 0,802. Aktywność 80% startowej po transporcie.
I-131 (T₁/₂ = 8,02 d): f = exp(−ln2/8,02 × 4,05) = exp(−0,350) = 0,705. Aktywność 70,5% startowej.

Stężenie przy stacji: dla 5 kt wybuchu, aktywność I-131 ~10¹⁵ Bq uwolnionego. Rozmywanie w kolumnie powietrza (stożek dyspersji): przy 7000 km i prędkości 20 m/s, objętość efektywna ~10¹⁸ m³. Stężenie: ~10¹⁵/10¹⁸ = 10⁻³ Bq/m³ = 1 mBq/m³.

Typowa czułość stacji IMS dla I-131 na filtrze: ~0,1–1 mBq/m³ przy 24h próbkowaniu. Wniosek: 5 kt wybuch mógłby być ledwo wykrywalny przy sprzyjającym transporcie.

Przykład 3 — Bayesowski update po obserwacji braku radionuklidów

Zakładamy a priori prawdopodobieństwo: P(H_jądrowy) = 0,60 (silny sygnał bhangmetru daje wysoki prior).
Obserwacja: brak detekcji radionuklidów na 3 stacjach AFTAC.

Prawdopodobieństwo braku detekcji przy H_jądrowym: P(brak detect | H_jądrowy) = f_transport × f_threshold = 0,30 × 0,60 = 0,18.
(transport 30% szansy dostarczenia próbki do odpowiedniej stacji × 60% szansy przekroczenia progu detekcji 1979)

Prawdopodobieństwo braku detekcji przy H_artefakt: P(brak detect | H_artefakt) = 1,00 (artefakt = brak źródła jądrowego = zawsze brak radionuklidów).

Posterior według Bayesa:
P(H_jądrowy | brak) ∝ P(brak | H_jądrowy) × P(H_jądrowy) = 0,18 × 0,60 = 0,108
P(H_artefakt | brak) ∝ P(brak | H_artefakt) × P(H_artefakt) = 1,00 × 0,40 = 0,400

Normalizacja: 0,108 + 0,400 = 0,508.
P(H_jądrowy | brak) = 0,108 / 0,508 = 0,213.
P(H_artefakt | brak) = 0,400 / 0,508 = 0,787.

Wniosek: Brak radionuklidów obniżył prawdopodobieństwo hipotezy jądrowej z 60% do 21%, przy założeniu słabych możliwości transportu 1979. Gdyby transport był lepszy (szansa 80%), update byłby mniejszy: P(H_jądrowy|brak) = 0,60×(1−0,8×0,6) / [(1−0,48)×0,6 + 0,40] = ... ≈ 44%. Kluczowe: odpowiedź bayesowska silnie zależy od założeń o możliwościach transportu i detekcji.

Pytania otwarte

1. Bhangmetr Vela 6911 miał 9 lat w przestrzeni kosmicznej w 1979 roku. Jakie procesy starzenia (radiacyjne uszkodzenia CCD/PMT, zmiany transmisji optycznej, zmiana wzmocnienia) mogły wpłynąć na stosunek amplitud kanałów? Jak można było zrekonstruować charakterystykę po 9 latach bez kalibracji in situ?

2. Sieć IMS CTBTO posiada 40 stacji monitoringu ksenonu radioaktywnego (SAUNA/ARIX). Jakie izotopy ksenonu są najbardziej diagnostyczne dla wybuchu jądrowego vs wyciek reaktorowy vs naturalna emanacja z gleby? Jak stosunek Xe-131m/Xe-133 pomaga to rozróżnić?

3. Modelowanie transportu atmosferycznego (HYSPLIT, FLEXPART) używane przez CTBTO IDC oblicza wsteczną trajektorię dla wykrytych radionuklidów. Jakie są główne źródła niepewności w takim modelu dla Southern Hemisphere (Antarktyda, Południowy Atlantyk), szczególnie dla 1979 roku bez reanaliz ERA5?

4. Bayesowska analiza incydentu Vela zależy od prawdopodobieństwa a priori hipotezy jądrowej. Jakie nieradiometryczne dane (wywiadowcze, geopolityczne, komercyjne satelitarne zdjęcia) mogą wpływać na ten prior? Czy takie dane powinny być uwzględniane w naukowej analizie technicznej czy tylko w ocenie wywiadowczej?

5. CTBTO IDC automatycznie klasyfikuje zdarzenia sejsmiczne i oblicza parametry źródłowe. Dla wybuchu atmosferycznego o mocy 5 kt w odległości 5000 km od najbliższej stacji sejsmicznej, jakie parametry (głębokość, magnituda, stosunek Mo) wskazałyby na wybuch powietrzny zamiast trzęsienia ziemi?

6. Panel prezydencki Carter (1979) stwierdził, że hipoteza meteoroidu jest „najbardziej prawdopodobna", podczas gdy laboratoria (LLNL, LANL) skłaniały się ku hipotezie jądrowej. Jakie procedury zarządzania ekspertyzami powinny obowiązywać, żeby takie rozbieżności były formalizowane i archiwizowane, a nie zanikały w poufnych raportach?

7. Tryt z detonacji termojądrowej (boosted) trafia do atmosfery i powoli miesza się z globalnymi zasobami. W 1979 roku tło trytu było wysokie po dekadzie testów atmosferycznych. Czy pomiar trytu w opadach mógłby byłby dodatkowym wskaźnikiem wybuchu w 1979 roku, czy byłby przysłonięty przez historyczne tło?

8. Incydent Vela jest przypadkiem, gdzie polityczna niejawność części danych utrudnia naukowe rozstrzygnięcie. Jakie są argumenty za i przeciw odtajnieniu wszystkich danych z 1979 roku po 45 latach? Jakie standardy archiwizacji danych detekcyjnych w systemach monitoringu zbrojeń powinny obowiązywać dla potomnych?

Podsumowanie dydaktyczne

1. Bhangmetry Vela były instrumentami o wysokiej czułości, ale bez pozycjonowania przestrzennego. Dwa kanały optyczne z różnymi czasami odpowiedzi dawały kształt krzywej jasności, nie mapę. Sama kształt krzywej nie jest dowodem nieodwracalnym — wymaga potwierdzenia przez inne kanały.

2. Incydent Vela jest klasycznym studium przypadku problemu identyfikowalności: ten sam zbiór obserwacji jest spójny z więcej niż jedną hipotezą. Rozróżnienie wymaga niezależnych kanałów, co było ograniczone przez infrastrukturę 1979 roku.

3. Brak detekcji radionuklidów zmniejsza prawdopodobieństwo hipotezy jądrowej, ale jej nie falsyfikuje. Wartość dowodowa braku jest odwrotnie proporcjonalna do możliwości detekcji — niska czułość sieci, słaby transport, czy krótki czas reakcji zmniejszają wagę negatywnego wyniku.

4. Analiza Bayesowska jest właściwym formalizmem dla przypadków z niepewnym prior i niejednoznacznymi obserwacjami. Kluczowe: każda hipoteza musi mieć przypisane prawdopodobieństwo generowania każdej obserwacji. Niespełnienie tego wymogu prowadzi do pseudonauki.

5. Dzisiejsza sieć IMS (sejsmika, infradźwięki, hydroakustyka, radionuklidy) daje gęstą macierz detekcji, która w 1979 roku nie istniała. To zmienia możliwości rozstrzygnięcia podobnych incydentów — nie eliminuje niepewności, ale dramatycznie zmniejsza przestrzeń niejednoznaczności.

6. Rola polityki i niejawności w naukowo-technicznej analizie jest nieuchronna w kwestiach proliferacji. Ale dobre standardy archiwizacji i transparentności danych metrologicznych powinny być oddzielone od interpretacji wywiadowczej.

7. Polska jako strona CTBT dostarcza dane z PS-26 i IS-26 do IDC w Wiedniu. Zrozumienie możliwości i ograniczeń tych stacji jest ważną częścią świadomości roli Polski w globalnym systemie nieproliferacji.

8. Najważniejsza lekcja Vela dla dydaktyki: „niepewność" i „brak dowodu" to nie to samo co „dowód braku". Uczciwa analiza metrologiczna formalizuje niepewność i przedstawia alternatywne hipotezy z oceną ich wagi — nie wydaje werdyktu bez podstaw.

Co wiemy jako punkt wyjścia

W publicznie dostępnej literaturze stałe są przynajmniej trzy elementy. Po pierwsze, satelita Vela 6911 zarejestrował sygnał optyczny typu podwójnego błysku. Po drugie, zdarzenie przypisano do bardzo szerokiego obszaru obserwacji obejmującego południowy Atlantyk, Ocean Indyjski, południową Afrykę i część Antarktyki. Po trzecie, późniejsza interpretacja była sporna: część analiz wspierała hipotezę jądrową, a panel prezydencki uznał, że sygnał prawdopodobnie nie pochodził od detonacji jądrowej.^1,2

To nie jest sytuacja "brak danych". To jest sytuacja "dane nie układają się w jeden publicznie bezsporny obraz". W metrologii to dużo ciekawsze. Jednoznaczne przypadki uczą mniej niż przypadki graniczne.

Najbardziej podstawowa lekcja brzmi: detektor może zarejestrować prawdziwy sygnał, a mimo to interpretacja może pozostać niepewna. Przyczyną może być geometria, tło, stan instrumentu, brak danych z innych kanałów albo polityczna niejawność części materiału.

Sygnał optyczny

Sygnał optyczny był ważny dlatego, że bhangmetry Vela miały rozpoznawać podwójny błysk atmosferycznej detonacji jądrowej. Taki przebieg wynika z fizyki kuli ognistej: szybki pierwszy impuls, chwilowe osłabienie, a potem drugi, wolniejszy wzrost jasności po przejściu etapu breakaway.^3,4

National Security Archive przytacza, że w poprzednich 41 przypadkach, w których satelita Vela zarejestrował taki podwójny błysk, późniejsze dane potwierdzały detonację jądrową.¹ To wzmacnia wagę sygnału, ale jej nie absolutyzuje. Statystyka działania instrumentu nie eliminuje możliwości rzadkiego artefaktu.

W przypadku 22 września 1979 pojawił się problem proporcji między dwoma kanałami bhangmetru. Nie oczekiwano identycznych wartości liczbowych, bo kanały nie musiały mieć tej samej czułości, ale oczekiwano podobnego stosunku między nimi jak w poprzednich jednoznacznych detonacjach. Ten stosunek miał odbiegać od doświadczenia.^1,2

Co oznacza niezgodność kanałów

Niezgodność kanałów nie mówi automatycznie: "to nie był wybuch". Mówi: "model pomiaru jest niepełny". Trzeba wtedy rozważyć kilka klas hipotez:

• realny wybuch, ale nietypowa geometria obserwacji,
• realny wybuch, ale częściowe przesłonięcie lub rozproszenie przez chmury,
• błąd albo starzenie jednego kanału,
• artefakt wywołany cząstką, mikrometeoroidem lub odbiciem światła,
• zjawisko naturalne o rzadkim przebiegu,
• złożenie kilku niezależnych zdarzeń.

Metrologicznie to jest problem identyfikowalności. Ten sam wynik z dwóch kanałów może mieć więcej niż jedną przyczynę. Aby rozdzielić hipotezy, potrzebne są kanały niezależne od optyki.

Hipoteza meteoroidu

Panel prezydencki rozważał scenariusz, w którym meteoroid uderza w satelitę, wyrzuca drobne cząstki, a odbite światło słoneczne tworzy sygnał podobny do podwójnego błysku.¹ Taka hipoteza jest atrakcyjna, bo tłumaczy, dlaczego sygnał pochodziłby z samej geometrii satelita-Słońce-detektor, a nie z powierzchni Ziemi.

Jej słabość polega na tym, że musi odtworzyć nie tylko "jakiś błysk", lecz konkretny przebieg podobny do znanej krzywej jądrowej. National Security Archive przytacza dokumenty, w których część analiz uznawała takie scenariusze za mało prawdopodobne, a niektóre raporty wskazywały, że hipoteza jądrowa lepiej pasowała do danych.¹

Właśnie tu widać różnicę między możliwością a prawdopodobieństwem. Rzadki artefakt może istnieć fizycznie. Pytanie brzmi, czy jest bardziej prawdopodobny niż zdarzenie, które instrument miał wykrywać.

Brak radionuklidów

Jednym z najsilniejszych argumentów przeciw rozstrzygnięciu jądrowemu był brak publicznie jednoznacznego potwierdzenia radiochemicznego. Po potencjalnym atmosferycznym lub powierzchniowym wybuchu można oczekiwać produktów rozszczepienia w powietrzu, opadzie albo próbkach środowiskowych. W przypadku Vela publiczny obraz pozostaje niejednoznaczny.^1,2

Brak detekcji nie jest jednak automatycznie dowodem braku zdarzenia. Może znaczyć:

• zdarzenia nie było,
• próbki pobrano w złym miejscu lub czasie,
• uwolnienie było małe,
• transport atmosferyczny ominął stacje i trasy poboru,
• radionuklidy rozcieńczyły się poniżej granicy wykrywalności,
• dane istnieją, ale nie są publiczne,
• sygnały były słabe i sporne.

To nie jest wygodna lista, ale jest uczciwa metrologicznie. Detekcja radionuklidów ma ogromną wartość dowodową, gdy występuje. Jej brak trzeba interpretować przez model transportu, czułość, czas i geometrię poboru.

Brak mocnych danych sejsmicznych i hydroakustycznych

Wybuch atmosferyczny lub nadwodny może sprzęgać się z oceanem, atmosferą i ziemią słabiej niż podziemna próba. Dlatego nie należy oczekiwać tak prostego sygnału jak przy teście podziemnym DPRK. Mimo to niezależne dane sejsmiczne, hydroakustyczne lub infradźwiękowe byłyby bardzo cenne.

National Security Archive opisuje poszukiwania danych w systemach takich jak SOSUS, MILS, satelity meteorologiczne i inne kanały obserwacji. Część raportów nie znajdowała potwierdzenia, część wskazywała sygnały słabe lub zaszumione, a część kwestionowała wnioski panelu prezydenckiego.¹

To jest przykład problemu korelacji. Jeśli kanał optyczny mówi "zdarzenie mogło być tu i teraz", to hydroakustyka lub sejsmika muszą nie tylko coś wykryć, lecz wykryć coś zgodnego czasowo, geometrycznie i energetycznie z hipotezą.

Ionosfera, pogoda i pobór próbek

Po Vela analizowano także dane cywilne i naukowe: obserwacje jonosferyczne, satelity meteorologiczne, ślady możliwego opadu i materiały środowiskowe. National Security Archive wskazuje m.in. anomalną falę jonosferyczną z Arecibo i intensywne próby poboru materiału przez samoloty AFTAC oraz inne działania terenowe.¹

To pokazuje, jak szeroki staje się problem po niejednoznacznym alarmie. Nie wystarczy odczytać jeden instrument. Trzeba zadać pytania o atmosferę, ocean, jonosferę, promieniowanie, ścieżki transportu i lokalne tło.

Z punktu widzenia studenta ważny jest porządek dowodowy. Dane pomocnicze nie mają tej samej wagi. Anomalia jonosferyczna może być interesująca, ale nie jest tym samym co widmo produktów rozszczepienia. Brak potwierdzenia z jednego satelity nie jest tym samym co brak zdarzenia, jeśli geometria była niekorzystna.

Dlaczego oceny ekspertów się rozeszły

Różne zespoły mogą dojść do różnych wniosków bez złej wiary. Przyczyną może być:

• dostęp do różnych danych,
• różne założenia o niezawodności instrumentu,
• inne priorytety: minimalizacja fałszywego alarmu albo minimalizacja przeoczenia,
• różne modele tła,
• różne rozumienie prawdopodobieństwa hipotezy a priori,
• niejawność części danych,
• presja polityczna na ostrożność publiczną.

Panel prezydencki miał powód, aby bardzo ostrożnie traktować sygnał, którego nie dało się publicznie domknąć. Z drugiej strony laboratoria i analitycy znający historię bhangmetrów mogli uważać, że optyczny ślad jest zbyt charakterystyczny, by łatwo przyjąć hipotezę artefaktu.¹

To jest dobry przykład, dlaczego w nauce i kontroli zbrojeń warto oddzielać trzy zdania:

• "zdarzenie było jądrowe",
• "zdarzenie prawdopodobnie było jądrowe",
• "publiczne dane nie wystarczają do rozstrzygnięcia".

Podejście bayesowskie

Najczystszy sposób myślenia o Vela jest bayesowski, nawet jeśli nie znamy liczb. Mamy hipotezy:

• H_N: doszło do detonacji jądrowej,
• H_A: doszło do artefaktu lub zjawiska naturalnego,
• H_M: doszło do innego zdarzenia technicznego lub mieszanego.

Następnie pytamy, jak prawdopodobne są obserwacje przy każdej hipotezie:

• podwójny błysk,
• nietypowy stosunek kanałów,
• brak mocnego publicznego potwierdzenia radionuklidowego,
• brak jednoznacznego potwierdzenia falowego,
• geopolityczny kontekst RPA i Izraela,
• możliwe dane jonosferyczne lub środowiskowe,
• historia wcześniejszych detekcji Vela.

Nie chodzi o podstawienie liczb z powietrza. Chodzi o dyscyplinę: każda obserwacja powinna zmieniać ocenę hipotez, ale żadna pojedyncza obserwacja nie powinna sama udawać pełnego dowodu.

Macierz dowodowa

Dobrą formą notatki analitycznej byłaby macierz:

Obserwacja	Wspiera H_N	Wspiera H_A	Uwagi metrologiczne
Podwójny błysk	silnie	słabo/średnio	zależy od zaufania do bhangmetru
Nietypowy stosunek kanałów	osłabia	wzmacnia	może wynikać też z geometrii
Brak radionuklidów	osłabia	wzmacnia	zależy od transportu i poboru
Brak jednoznacznej sejsmiki	osłabia	wzmacnia	zależy od typu i miejsca zdarzenia
Geopolityczny kontekst	wzmacnia a priori	neutralnie	nie jest pomiarem fizycznym
Hipoteza meteoroidu	konkuruje	wzmacnia	wymaga odtworzenia krzywej czasowej

Taka tabela nie rozwiązuje sporu, ale wymusza uczciwość. Nie pozwala ukryć słabego punktu we frazie "wszyscy wiedzą" ani "nic nie udowodniono".

Czym Vela różni się od przypadków DPRK

Testy DPRK są dobrym kontrastem. Tam sejsmika bardzo szybko lokalizowała zdarzenia podziemne, a w części przypadków późniejsze radionuklidy wzmacniały interpretację jądrową.⁵ Kanały były zgodne z typem testu: podziemny wybuch -> sejsmika; późniejsze produkty rozszczepienia -> radionuklidy.

Vela jest przypadkiem odwrotnym. Najsilniejszy sygnał był optyczny, a reszta publicznych danych nie dała równie prostego domknięcia. To czyni go dużo trudniejszym metrologicznie.

Nie oznacza to, że Vela był słabszy jako podejrzenie. Oznacza, że był słabszy jako publicznie zamknięty przypadek dowodowy.

Jak wyglądałby taki przypadek dziś

Dzisiejszy świat ma więcej satelitów, lepsze dane meteorologiczne, globalny IMS, większe archiwa falowe i lepsze narzędzia transportu atmosferycznego. Podejrzany podwójny błysk byłby natychmiast zestawiany z:

• danymi innych satelitów optycznych i podczerwonych,
• sejsmiką,
• infradźwiękami,
• hydroakustyką,
• radionuklidami,
• modelami transportu atmosferycznego,
• danymi pogodowymi i chmurami,
• obserwacjami jonosferycznymi,
• kontekstem wywiadowczym.

To nie znaczy, że każdy współczesny przypadek byłby łatwy. Oznacza tylko, że macierz dowodowa byłaby gęstsza. Paradoksalnie więcej danych może też dać więcej sporów, jeśli poszczególne kanały są słabe i częściowo sprzeczne.

Lekcja o "braku dowodu"

W debatach o Vela często miesza się dwa zdania:

• "nie ma publicznego, rozstrzygającego dowodu testu",
• "wiadomo, że testu nie było".

To nie są te same zdania. Pierwsze jest ostrożne. Drugie jest mocniejsze i wymaga własnych dowodów. Analogicznie:

• "panel uznał hipotezę meteoroidu za najbardziej prawdopodobną" nie oznacza, że wszystkie inne analizy się z tym zgodziły,
• "część ekspertów uważała sygnał za jądrowy" nie oznacza, że sprawa została publicznie dowiedziona,
• "nie wykryto jednoznacznych radionuklidów" nie oznacza automatycznie, że nie było uwolnienia.

To są podstawowe reguły pracy z niepewnością. W kontroli zbrojeń są równie ważne jak same detektory.

Dlaczego to nie jest tylko historia

Vela 1979 uczy, jak projektować i czytać systemy monitoringu:

• trzeba znać historię fałszywych alarmów,
• trzeba mieć redundancję kanałów,
• trzeba zachować surowe dane,
• trzeba dokumentować kalibrację i stan instrumentu,
• trzeba rozdzielać dane fizyczne od hipotez politycznych,
• trzeba umieć powiedzieć "nie rozstrzygamy",
• trzeba budować procedury korelacji między sensorami.

W tym sensie Vela jest bardziej wartościowy dydaktycznie niż wiele jednoznacznych testów. Pokazuje nie triumf detektora, ale granice pewności.

Podsumowanie

Najkrótszy opis brzmi: Vela 6911 zobaczył sygnał optyczny, który bardzo poważnie wyglądał jak podwójny błysk atmosferycznej detonacji jądrowej. Jednocześnie proporcje kanałów, brak publicznie jednoznacznych potwierdzeń i możliwość rzadkiego artefaktu sprawiły, że przypadek nie został zamknięty w sposób powszechnie akceptowany.

Dla tego serwisu najważniejsza jest nie sensacja, lecz metoda. Vela pokazuje, że w świecie jądrowym dobry pomiar to nie tylko czuły czujnik. To także model błędu, alternatywne hipotezy, niezależne kanały, archiwizacja danych i uczciwa praca z niepewnością.

Dodatkowe materiały multimedialne

Ćwiczenia praktyczne

Pierwsze ćwiczenie: zbudować tabelę hipotez H_N, H_A, H_M i przypisać do każdej obserwacji z artykułu znak +, - albo 0. Następnie porównać, które obserwacje są naprawdę niezależne.

Drugie ćwiczenie: napisać dwa krótkie streszczenia tego samego przypadku. Pierwsze ma brzmieć: "dane są zgodne z hipotezą testu, ale nie rozstrzygają jej". Drugie: "dane są zgodne z artefaktem, ale nie obalają hipotezy testu". Porównać, które sformułowania są uczciwe, a które sugerują zbyt dużo.

Trzecie ćwiczenie: przeanalizować, dlaczego brak radionuklidów ma różną wagę dla wybuchu atmosferycznego, podziemnego, nadwodnego i bardzo dobrze odizolowanego. Nie obliczać dawki ani opadu; chodzi o logikę detekcji.

Czwarte ćwiczenie: dla syntetycznego sygnału optycznego z dwóch kanałów obliczyć stosunek amplitud drugiego piku do pierwszego w obu kanałach. Następnie opisać, co oznacza rozbieżność stosunków.

Piąte ćwiczenie: zaprojektować procedurę "co robimy po alarmie optycznym" bez politycznej interpretacji. Procedura ma obejmować dane satelitarne, sejsmikę, infradźwięki, hydroakustykę, radionuklidy, meteorologię i archiwizację surowych danych.

Szóste ćwiczenie: porównać Vela 1979 z jednym z testów DPRK opisanych przez CTBTO. Wskazać, który przypadek ma mocniejszą fuzję danych i dlaczego.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły

• Bhangmetry i podwójny błysk: optyczne wykrywanie detonacji jądrowych
• Międzynarodowy System Monitoringu CTBT: sejsmika, infradźwięki, hydroakustyka i radionuklidy
• Incydent Vela i systemy monitorowania eksplozji
• Kula ognista i Wilson Cloud
• Nuclear Forensics (Forensyka jądrowa)