Na obszarze badawczym historyczne zanieczyszczenie CH było najwyższe w obrębie centralnej grupy studni (MW-25-MW-30) oraz w pobliżu Sherman Road (Rys. 5). W 1983 roku z terenu składowiska odpadów (na północ od obszaru badawczego) usunięto znaczne ilości zanieczyszczeń, a w 2001 roku przeprowadzono dodatkowe wykopy. Stężenia CH spadły po usunięciu źródła, szczególnie w pobliżu dawnych wykopów (Sherman Road), jednak w regionie centralnej grupy studni nadal występuje trwały pióropusz zanieczyszczeń. Znane jest, że sezonowe opady deszczu przejściowo zwiększają stężenia CH, a pozostałe zanieczyszczenia desorbują z gleb27. Gleby w tym obszarze to głównie dawne piaski pogłębiarskie. Możliwa interferencja z opisaną metodą mogłaby wystąpić w przypadku obecności starodawnych skał węglanowych i bardzo niskiego pH wód gruntowych (<~5). Mogłoby to prowadzić do rozpuszczania węglanów i pojawienia się sygnału starodawnego w generowanym CO2. W okolicy nie znane są znaczące pokłady CaCO3, niemniej jednak zmierzono zawartość kationów oraz pH, a następnie poddano je analizie regresji i analizie głównych składowych (PCA). Główną obawą było to, że niskie pH może sprzyjać rozpuszczaniu węglanu wapnia (CaCO3), co mogłoby zaburzyć analizę radiowęglową (starodawne skały węglanowe mogłyby dostarczyć starodawnego CO2 w przypadku rozpuszczenia przez wody kwasowe). Zawartość Na+ była nieco wyższa po południowej stronie terenu (najbliżej oceanu), jednak żadne wartości nie znajdowały się w zakresie wskazującym na znaczną intruzję wody morskiej. Stężenia jonów wapnia wahały się od 8.0 do 58 mg L-1. Relacja między stężeniem jonów wapnia a pH nie wskazywała na rozpuszczanie węglanów (r2 < 0.3). Wykresy bi-plot PCA nie wykazały silnych obciążeń dla żadnej ze zmiennych. Różnice między studniami również nie wskazywały na rozpuszczanie węglanów (Rys. 6). Ta analiza konformacyjna powinna być uznana za krytyczną przy adaptacji metodologii do nowych obszarów — szczególnie tych, których geologia regionalna wskazuje na znaczące formacje skalne węglanowe.
Szybkość produkcji CO2 wynosiła od 0 do 34 mg CO2 d-1. Produkcja CO2 była najniższa w centralnej grupie otworów w obszarze, w którym historyczne zanieczyszczenie było najwyższe (Rys. 5). Produkcja CO2 w otworze MW-01 (otwór tła — niepokazany, ale znajdujący się ok. 500 metrów na północny zachód od głównej grupy otworów) była bardzo wysoka i wynosiła 31 mg CO2 d-1). Błędy standardowe dla powtórnych analiz respiracji wynosiły od 0,03 do 6% CO2, a średnio mniej niż 1% (0,98). Dwa 2-tygodniowe okresy pomiarów w sezonie suchym zostały uśrednione do dalszych obliczeń. Pomiary respiracji nie różniły się znacząco pomiędzy poszczególnymi 2-tygodniowymi okresami. Błąd standardowy respiracji między okresami wynosił od <1 do 51%, ale średnio 13% (Tabela 1). Uśrednienie respiracji pozwoliło na obliczenie pojedynczej objętości usuniętego CH w okresie jednego miesiąca. Otwór tła (MW-01) miał wiek radiowęglowy 1280 lat przed teraźniejszością (ybp) lub 85 procent węgla nowoczesnego (pMC) — co mieści się w typowym zakresie dla starej organicznej materii glebowej26. Wartość z tego otworu posłużyła jako tło dla modelu mieszania izotopowego. Ponownie, ze względu na to, że pobieranie próbek ograniczono do łącznego okresu jednego miesiąca, wykorzystano dwa następujące po sobie okresy w tym samym sezonie, aby „reprezentować” sezon suchy — który zazwyczaj uznaje się za okres najbardziej stagnacyjnych warunków, a tym samym konserwatywny dla ekstrapolowanych szacunków. Podobnie jak w przypadku szybkości produkcji DIC, pomiary radiowęglowe były zbliżone pomiędzy poszczególnymi 2-tygodniowymi okresami. Błąd standardowy między okresami wynosił od 0,25 do 18%, a średnio 6%. Wiek radiowęglowy CO2 wynosił od ~34 do 85 pMC lub od ~1340 do 8700 ybp (Tabela 1). Próbki z MW-27 i MW-32, w przypadku których podejrzewano nieszczelność pompy, wykazały nowoczesne wartości radiowęglowe i tym samym potwierdzono ich zanieczyszczenie. Próbki te nie zostały włączone do dalszych analiz.
Poprzednie raporty zostały wykorzystane do określenia właściwości hydraulicznych wód podziemnych oraz właściwości rozpuszczonego CO2 w celu opracowania modelu ZOI26,27 (Tabela 2). Dane pogodowe (z lat 2007, 2011 i 2012) ze stacji CIMIS w San Diego (ID stacji 184) zostały użyte do oszacowania szybkości zasilania warstwy wodonośnej. Dane pływowe z tego samego okresu ze stacji NOAA w San Diego (ID stacji: 9410170) zostały wykorzystane do zdefiniowania warunków brzegowych. Kalibracja modelu zakładała stały gradient hydrauliczny oraz stałe szybkości zbierania CO2. W parametryzacji modelu pomogły dodatkowe symulacje z różnicującymi się średnimi szybkościami zbierania CO2 i początkowym poziomem tła CO2, połączone ze zwiększeniem gradientu hydraulicznego o 10%. Dodatkowa symulacja z wykorzystaniem średniej szybkości zbierania CO2 wykazała około 46-procentowy wzrost szacowanego tła CO2 (tj. wzrost z 6,5 do 9,5 g m-3), jeżeli szybkość zbierania zmieniła się z 0,00530 (+10%) do 0,00434 g hr-1 (-10%) w ciągu 2-tygodniowego okresu pobierania próbek (Tabela 3). Założenia modelu ZOI obejmowały pomijalną produkcję CO2 wynikającą z degradacji CH w okresie zbierania oraz jednorodny początkowy rozkład CO2 w celu opracowania końcowej symulacji (Rys. 7). Szybkość reakcji CO2 dla badanego obszaru może być niedoszacowana.
Wykorzystując tempo produkcji CO2, ilość CO2 przypisywaną degradacji CH oraz szacunki z modelu ZOI, obliczono masowy usuw proces CH w każdym otworze w jednostce czasu. Dane z Table 1 zostały wykorzystane w dwuskładnikowym modelu mieszania (równ. (1)) w celu wyznaczenia fpet dla każdego otworu. Ponieważ w obrębie terenu stwierdzono jedynie zanieczyszczenie CH i nie znaleziono żadnego innego źródła CO2 na terenie lub w pobliżu terenu, założono, że degradacja CH jest głównym źródłem CO2. Wartości fpet na terenie obiektu wahały się od 1 do 60% (Table 4). Udział ten przeliczono na zawartość węgla i pomnożono przez tempo produkcji CO2, aby obliczyć tempo degradacji CH (Table 4). Wykorzystując objętość ZOI (Table 3), określono tempo degradacji zanieczyszczeń w jednostce czasu i objętości (Table 4). Wartość ta mieściła się w przedziale od 0 do 32 mg C m-3 d-1 (Table 4). Degradacja CH była najniższa w regionach o najwyższym historycznym zanieczyszczeniu CH (MW-25 - MW-30). Najwyższą degradację CH zmierzono w otworach położonych w pobliżu obrzeży terenu (w okolicy Sherman Road). Produkcja CO2 była w tych obszarach wyższa, podczas gdy fpet wskazywała na znaczną wymianę CH (Fig. 8).

Rycina 1. Uszczelnianie i przygotowanie pomp recyrkulacyjnych. Uszczelnianie pomp recyrkulacyjnych do zastosowań w terenie.

Rycina 2. Pułapki z NaOH przygotowane do zastosowania w terenie. Butelki po surowicy o pojemności 120 ml z dodaną pułapką z NaOH i zamknięte poprzez zakładkę (crimp seal).

Rysunek 3. Układ terenowy. Przewody prowadzone do wyposażonych studni (lewo), pułapka rozmieszczona przy studni (prawy górny róg) oraz system dystrybucji energii słonecznej (prawy dolny róg). Studnie są wyposażone w terenie w systemy poboru (w tym okablowanie, dystrybucję zasilania oraz pompy/pułapki).

Rycina 4. Zmodyfikowane pokrywy otworów ukazujące linie recyrkulacji gazu. Na tej rycinie przedstawiono pokrywy otworów zmodyfikowane o linie doprowadzania i powrotu gazu.

Rycina 5. Historyczne zanieczyszczenie chlorowanymi węglowodorami (µg L-1). Na rycinie przedstawiono historyczne zanieczyszczenie chlorowanymi węglowodorami na terenie badawczym.

Rysunek 6. Biplot PCA wykazujący brak współkorelacji między kationami rozpuszczonymi a pH. Na tym rysunku przedstawiono biplot wyników i ładunków PCA utworzony na podstawie danych hydrogeologicznych (pH i kationów) dla terenu badań.

Rycina 7. Skalibrowany model ZOI dla średniej szybkości gromadzenia CO2 (0,0048 g m-3). Skalibrowane tło stężenia CO2 wynosiło 6,5 g m-3, a progowe stężenie ZOI wynosiło 6,18 g m-3 (ciągła czarna linia). Średnice podłużna i poprzeczna ZOI wynosiły odpowiednio 2,28 m i 0,72 m. Głębokość ZOI wynosiła 0,12 m. Zmodyfikowano na podstawie 18. Rycina ta przedstawia graficzną reprezentację modelu ZOI w 3 wymiarach.

Rycina 8. Szybkość degradacji zanieczyszczeń na jednostkę czasu i jednostkę powierzchni. Zmodyfikowano na podstawie 18. Przedstawiono interpolowaną szybkość degradacji CH na obszarze badawczym w okresie pobierania próbek.

Wideo 1. Rozwój ZOI przy użyciu symulacji MT3DMS23 - MODFLOW (kliknij prawym przyciskiem myszy, aby pobrać). Pobierz, zainstaluj, zainicjuj i stwórz symulację dla ZOI.
| Otwór | δ13C
(‰VPDB) | Δ14C
(‰) | Wiek konwencjonalny
(ybp) | Procent węgla współczesnego
(pMC) |
| MW-01 | -34 | -147 | 1280 | 85 |
| MW-21 | -28 | -663 | 8730 | 34 |
| MW-25 | -23 | -153 | 1340 | 85 |
| MW-26 | -25 | -298 | 2845 | 70 |
| MW-27 | -18 | N.D.* | N.D.* | N.D.* |
| MW-28 | -25 | -190 | 1695 | 81 |
| MW-30 | -35 | -254 | 2365 | 75 |
| MW-32 | -20 | N.D.* | N.D.* | N.D.* |
| MW-34 | -32 | -283 | 2670 | 72 |
| MW-35 | -25 | -598 | 7320 | 40 |
| MW-38 | -32 | -354 | 3515 | 65 |
| MW-41 | -28 | -232 | 2125 | 77 |
| MW-42 | -23 | -482 | 5280 | 52 |
| *N.D. Brak danych – nieszczelna pompa |
Tabela 1. Pomiary i przeliczniki izotopów CO2. Pomiary stabilnych izotopów i radiowęgla CO2 oraz przeliczniki na jednostki użyte w manuskrypcie.
| Parametr | Jednostki | Wartość |
| Hydrologia | | |
| Przewodność hydrauliczna | ml hr-1 | 0,44 (warstwa wodonośna) |
| 10 (otwór studzienny) |
| Porowatość (warstwa wodonośna) | | 0,48 (warstwa wodonośna) |
| 0,99 (otwór studzienny) |
| Gęstość objętościowa | g cm-3 | 1,4 |
| Wydajność właściwa | cm3 cm-3 | 0,2 |
| Gradient hydrauliczny | m m-1 | 0,015 |
| Transport substancji rozpuszczonej CO2 | | |
| Współczynnik dyfuzji | m2 hr-1 | 5,77 x 10-5 |
| Wzdłużna | m | 6,1 |
| Dyspersyjność |
| Poprzeczna pozioma | m | 0,61 |
| Dyspersyjność |
| Poprzeczna pionowa | m | 0,061 |
| Dyspersyjność |
| CO2 w gazie glebowym | % | 0,56 |
Tabela 2. Parametry modelu ZOI. Parametry wykorzystane w modelu ZOI i symulacjach.
| Poziom szybkości zbierania | Szybkość zbierania | Stężenie tła | Wielkość ZOI |
| Wzdłużna | Poprzeczna | Głębokość | Objętość |
| (g/hr) | (g/m3) | (m) | | | (m3) |
| Maksymalny | 0.0131 | 17.6 | 2.47 | 0.77 | 0.13 | 0.193 |
| Średni | 0.0048 | 6.5 | 2.28 | 0.72 | 0.12 | 0.176 |
| Minimalny | 0.0003 | 4 | 2.16 | 0.68 | 0.11 | 0.149 |
Tabela 3. Wyniki modelu ZOI. Wyniki modelu dla ZOI. Tabela ta opisuje trójwymiarową objętość dla ZOI.
| Studnia | fpozytonowa tomografia emisyjna
(%) | Szybkość degradacji zanieczyszczeń
(mg C d-1 ±10%) | Degradacja zanieczyszczeń na jednostkę czasu i objętości
(mg C m-3 d-1 ±15%) |
| MW-01 | 0 | N.D. | N.D. |
| MW-21 | 60 | 5.6 | 32 |
| MW-25¥ | 1 | 0 | 0 |
| MW-26 | 18 | 0.18 | 1 |
| MW-28 | 5 | 0.017 | 0.098 |
| MW-30 | 12 | 0.34 | 1.9 |
| MW-34 | 16 | 0.1 | 0.58 |
| MW-35 | 53 | 3.6 | 20 |
| MW-38 | 24 | 1.4 | 8.1 |
| MW-41 | 10 | 0.44 | 2.5 |
| MW-42 | 39 | 1.7 | 9.8 |
| N.D. Nie dotyczy – MW-01 użyto jako tło (np. brak kontaminacji) |
| ¥Zakładając, że proces ten jest napędzany wyłącznie przez stan równowagi (np., brak oddychania) |
Tabela 4. Skalowane szacunki degradacji zanieczyszczeń. Szacunki degradacji zanieczyszczeń na jednostkę czasu i jednostkę objętości dla pobranych próbek z otworów wiertniczych.