Pochodzenie wody na Ziemi

Skąd pochodzi woda w ziemskich oceanach i czy mogły ją dostarczyć komety?

Porównanie stosunku deuteru do wodoru (D/H) w wodzie kometarnej i ziemskiej pozwala przetestować tę hipotezę. Jako pozostałości po procesie powstawania planet komety przechowują pierwotną materię, której badanie może pomóc nam poznać zarówno pochodzenie ziemskiej wody, jak i przebieg formowania się Układu Słonecznego oraz innych układów planetarnych. Zawarta w nich woda niesie izotopowy zapis warunków fizyko-chemicznych, w których powstała i ulegała dalszemu przetworzeniu. Ziemskim punktem odniesienia jest wiedeński wzorzec średniej wody oceanicznej (Vienna Standard Mean Ocean Water, VSMOW), dla którego D/H = 1,5576 × 10⁻⁴. Wyniki zgodne z tą wartością wspierają hipotezę o udziale komet w dostarczeniu wody do ziemskich oceanów. Pomimo dużego znaczenia naukowego takie porównania nadal napotykają trudności obserwacyjne. Dostępna obecnie próbka obejmuje zaledwie 13 komet Układu Słonecznego, w tym dwie zbadane in situ przez sondy Giotto i Rosetta Europejskiej Agencji Kosmicznej (ESA). Ponieważ związki zawierające deuter wytwarzają niezwykle słabe sygnały spektroskopowe, udane pomiary zdalne ograniczały się głównie do obiektów obserwowanych w wyjątkowo sprzyjających warunkach, z których wiele było widocznych gołym okiem. Nawet przy wykorzystaniu najpotężniejszych dostępnych instrumentów, dotychczasowe znaczne niepewności i różnice metodologiczne utrudniają interpretację wyników.

Rysunek 1 podsumowuje postępy w pomiarach zawartości deuteru w wodzie kometarnej. Badania te rozpoczęły się wraz z przelotem sondy Giotto w pobliżu komety 1P/Halley w 1986 roku. Dwa instrumenty — spektrometr mas jonów (Ion Mass Spectrometer, IMS) i spektrometr mas cząsteczek neutralnych (Neutral Mass Spectrometer, NMS) — dostarczyły niezależnych oszacowań stosunku D/H. Pierwotne analizy dały zgodne ze sobą wartości wynoszące około 3 × 10⁻⁴, czyli mniej więcej dwukrotnie więcej niż w ziemskich oceanach (Balsiger et al. 1995; Eberhardt et al. 1995). Późniejsza ponowna analiza danych NMS przeprowadzona przez Brown et al. (2012) doprowadziła do skorygowania wyniku do (2,1 ± 0,3) × 10⁻⁴, bliższego VSMOW. Następnie wyjątkowo jasne pojawienia się komet C/1996 B2 (Hyakutake) i C/1995 O1 (Hale-Bopp) umożliwiły pierwsze naziemne detekcje kometarnego HDO — cząsteczki wody, w której jeden atom wodoru jest zastąpiony deuterem. Obserwacje jej przejścia przy częstotliwości 464,925 GHz, wykonane odpowiednio za pomocą Caltech Submillimeter Observatory (CSO) i James Clerk Maxwell Telescope (JCMT), pozwoliły wyznaczyć stosunki D/H zasadniczo zgodne z pierwotnymi wynikami Giotto (Bockelée-Morvan et al. 1998; Meier et al. 1998). Pomiary te wymagały jednak niezależnego oszacowania tempa produkcji wody na podstawie obserwacji produktów jej dysocjacji, co wiązało się z różnicami w czasie obserwacji i badanych obszarach komy.

Dwie inne komety, C/2002 T7 (LINEAR) i C/2001 Q4 (NEAT), zbadano za pomocą produktów powstających podczas rozpadu cząsteczek wody pod wpływem promieniowania słonecznego. Zwykła woda rozpada się na OH i H, natomiast HDO może wytwarzać albo OD i H, albo OH i D. Za pomocą spektrografu Ultraviolet and Visual Echelle Spectrograph (UVES) przy Very Large Telescope (VLT) Hutsemékers et al. (2008) obserwowali pasma ultrafioletowe OH i OD w pobliżu 310 nm w komecie C/2002 T7. Nie wykryto żadnej pojedynczej linii OD, jednak połączenie najsilniejszych oczekiwanych linii umożliwiło detekcję na poziomie około 3σ. Zmierzony stosunek OD/OH wynoszący (3,5 ± 1,0) × 10⁻⁴ przeliczono, z uwzględnieniem względnych prawdopodobieństw poszczególnych kanałów fotodysocjacji, na stosunek D/H w wodzie równy (2,5 ± 0,7) × 10⁻⁴. W przypadku C/2001 Q4 obserwacje wykonane w 2004 roku za pomocą instrumentu Space Telescope Imaging Spectrograph (STIS) na Kosmicznym Teleskopie Hubble’a były natomiast ukierunkowane na emisję Lyman-α atomowego wodoru i deuteru. Weaver et al. (2008) wyznaczyli D/H = (4,6 ± 1,4) × 10⁻⁴, a więc wartość znacznie wyższą od VSMOW. Kolejną podwyższoną wartość, (4,09 ± 1,45) × 10⁻⁴, uzyskano dla komety typu Halleya 8P/Tuttle na podstawie obserwacji wykonanych w 2008 roku za pomocą instrumentu Cryogenic Infrared Echelle Spectrograph (CRIRES) przy VLT, obejmujących przejścia HDO w pobliżu 3,7 μm (Villanueva et al. 2009). Wyniki te poszerzyły obserwowany zakres kometarnych stosunków D/H, choć ich niepewności pozostawały znaczne.

Pierwsze wyznaczenie stosunku D/H w komecie z rodziny Jowisza przyniosło wyraźnie odmienny wynik. Obserwacje 103P/Hartley 2 za pomocą instrumentu Heterodyne Instrument for the Far-Infrared (HIFI) na pokładzie należącego do ESA Kosmicznego Obserwatorium Herschela dały D/H = (1,61 ± 0,24) × 10⁻⁴, zgodne ze składem wody w ziemskich oceanach (Hartogh et al. 2011). W pomiarze wykorzystano przejście HDO przy częstotliwości 509,292 GHz. Zastosowanie podobnej metody do komety C/2009 P1 (Garradd), pochodzącej z Obłoku Oorta, dało D/H = (2,06 ± 0,22) × 10⁻⁴ — wartość pośrednią między VSMOW a około dwukrotnym wzbogaceniem stwierdzonym w najwcześniejszych pomiarach (Bockelée-Morvan et al. 2012). Wyniki te początkowo sugerowały możliwą różnicę izotopową między kometami z rodziny Jowisza, napływającymi z obszaru transneptunowego, a kometami z Obłoku Oorta. Późniejsze pomiary pokazały jednak, że prosty podział między tymi populacjami nie wyjaśnia dostępnych danych. Herschel znacznie poprawił spójność pomiarów, obserwując izotopologi wody niemal jednocześnie i przy zbliżonych rozmiarach wiązek teleskopu. Wyznaczenie referencyjnego tempa produkcji wody nadal wymagało jednak albo modelowania optycznie grubej linii H₂O, albo obserwacji rzadszego izotopologu H₂¹⁸O i przyjęcia stosunku izotopowego tlenu w celu określenia całkowitej obfitości wody.

Po zakończeniu pracy Herschela w 2013 roku pomiary zależały od pozostałych dostępnych instrumentów i sprzyjających pojawień się komet. Kometa C/2014 Q2 (Lovejoy) unaoczniła trudności w porównywaniu wyników uzyskanych różnymi technikami. Paganini et al. (2017) zmierzyli HDO i H₂O jednocześnie w podczerwieni za pomocą 10-metrowego teleskopu Keck II, uzyskując D/H = (3,02 ± 0,87) × 10⁻⁴, czyli 1,94 ± 0,56 wartości VSMOW. Biver et al. (2016) połączyli natomiast obserwacje linii HDO przy częstotliwości 241,561 GHz z 30-metrowego radioteleskopu IRAM z późniejszymi obserwacjami H₂O i H₂¹⁸O przy częstotliwościach 556,936 i 547,676 GHz z satelity Odin. Uzyskany przez nich wynik, około (1,4 ± 0,4) × 10⁻⁴, był zgodny ze składem ziemskiej wody. Oba oszacowania różniły się mniej więcej dwukrotnie i były zgodne dopiero na poziomie około 2σ. W tej samej pracy Biver et al. podali D/H = (6,5 ± 1,6) × 10⁻⁴ dla C/2012 F6 (Lemmon), na podstawie detekcji HDO na poziomie około 4σ. Wartość ta, przekraczająca VSMOW ponad czterokrotnie, silnie kontrastowała z ich wynikiem dla Lovejoy, mimo że obie komety pochodziły z Obłoku Oorta i zostały zbadane podobnymi metodami. Takie porównania mogą odzwierciedlać rzeczywiste zróżnicowanie, ale ujawniają również istotne ograniczenia obserwacyjne: pomiary wykonane w różnych terminach i przy różnych rozmiarach wiązek mogą obejmować odmienne fazy zmienności związanej z obrotem jądra oraz różny udział wody uwalnianej z sublimujących ziaren lodowych. Rozbieżność wyników dla Lovejoy nie stanowi zatem sama w sobie dowodu na zmienność składu izotopowego jej jądra w czasie.

Długotrwałe badania komety 67P/Churiumow–Gerasimenko prowadzone przez Rosettę stworzyły bezprecedensową możliwość pomiaru kometarnych obfitości izotopowych in situ. Korzystając ze spektrometru mas o podwójnym ogniskowaniu ROSINA DFMS, Müller et al. (2022) wyznaczyli stosunek D/H w wodzie na (5,01 ± 0,41) × 10⁻⁴, czyli około 3,2 razy więcej niż VSMOW, przy niewielkiej zmienności w analizowanych pomiarach. Rosetta umożliwiła również badania wzbogacenia w deuter innych cząsteczek, w tym metanu, etanu, propanu, butanu, metanolu, siarkowodoru i amoniaku. Jednak ponowna analiza ponad 4000 pomiarów stosunku D/H w wodzie, przeprowadzona przez Mandt et al. (2024), zakwestionowała interpretację wyników dotyczących wody. Autorzy zaproponowali, że preferencyjna adsorpcja HDO na ziarnach pyłu, a następnie uwalnianie wody wzbogaconej w deuter z ich lodowych otoczek mogą podwyższać stosunek D/H mierzony lokalnie w pobliżu jądra. Na podstawie pomiarów wykonanych po przejściu przez peryhelium, w odległości ponad 120 km od jądra, wyznaczyli D/H = (2,59 ± 0,36) × 10⁻⁴, czyli około 1,7 wartości VSMOW, co uznali za wynik bardziej reprezentatywny dla jądra. Interpretacja ta pozostaje przedmiotem dyskusji: późniejsza analiza Müller et al. (2026), uwzględniająca udział jonu hydroniowego oraz różnice w kalibracji skali mas, ponownie dała wyższą średnią wartość dla całej misji, wynoszącą (4,67 ± 0,59) × 10⁻⁴. Rozbieżność ta podkreśla znaczenie odróżnienia składu izotopowego całkowitego zasobu lodu w jądrze od lokalnych efektów zachodzących w komie oraz niepewności analizy nakładających się sygnałów w widmach masowych.

Bliskie zbliżenie komety z rodziny Jowisza 46P/Wirtanen w 2018 roku stworzyło kolejną możliwość porównania wody kometarnej z wodą ziemskich oceanów. Obserwacje za pomocą Stratospheric Observatory for Infrared Astronomy (SOFIA) obejmowały przejścia H₂¹⁸O i HDO w pobliżu 547 i 509 GHz. Uzyskany stosunek D/H = (1,61 ± 0,65) × 10⁻⁴ był zgodny z VSMOW, choć jego wyznaczenie wymagało przyjęcia stosunku izotopowego tlenu, podobnie jak w odpowiadających im pomiarach Herschela. Na podstawie tego wyniku Lis et al. (2019) zaproponowali możliwy związek między wartością D/H a procentem aktywnej sublimacyjnie powierzchni komety, potencjalnie powiązany z udziałem sublimujących ziaren lodowych. Kolejną wartość zbliżoną do ziemskiej podali Cordiner et al. (2025) dla komety typu Halleya 12P/Pons–Brooks: obserwacje ALMA wykonane w 2024 roku dały D/H = (1,71 ± 0,44) × 10⁻⁴. Chociaż ALMA dostarczyła map zarówno HDO, jak i H₂O, obie cząsteczki obserwowano w różnych terminach. W analizie wykorzystano zatem równoczesne obserwacje metanolu do powiązania tych zbiorów danych, a zasadność tego podejścia wsparto niezależną weryfikacją stosunku obfitości metanolu do wody. Obserwacje linii wody przy częstotliwości 183 GHz wymagały również starannej korekcji absorpcji w ziemskiej atmosferze. Inne poszukiwania, obejmujące między innymi komety 153P/Ikeya–Zhang, C/2007 N3 (Lulin), C/2012 S1 (ISON) i 45P/Honda–Mrkos–Pajdušáková, pozwoliły wyznaczyć jedynie górne granice stosunku D/H, bez uzyskania detekcji.

Udział komet w dostarczeniu wody na Ziemię pozostaje zatem nierozstrzygnięty. Pomiary zgodne z VSMOW pokazują, że niektóre komety zawierają wodę o składzie izotopowym zgodnym ze składem naszych oceanów, jednak niewielka próba i różnorodność metod obserwacyjnych uniemożliwiają wiarygodną ocenę całej populacji komet. Niskie stosunki sygnału do szumu, obserwacje prowadzone w różnych terminach, różnice w badanych obszarach komy, efekty grubości optycznej oraz założenia dotyczące referencyjnych obfitości izotopowych dodatkowo utrudniają porównania. Przezwyciężenie tych ograniczeń wymaga znacznie większej próby jednorodnych, wewnętrznie spójnych pomiarów o wysokim stosunku sygnału do szumu. Dedykowany instrument umożliwiający jednoczesny pomiar wodoru i deuteru w tym samym obszarze komy, i z odpowiednio dużą czułością, pozwoliłby dokonać zasadniczego postępu. Oprócz przetestowania udziału komet w powstaniu ziemskich oceanów taki zbiór danych dostarczyłby ograniczeń dla modeli chemii wody we wczesnym Układzie Słonecznym oraz powstawania i przemieszczania się planetozymali, z których wyłaniają się układy planetarne.