Na terenie badawczym historyczne zanieczyszczenie CH było najwyższe w obrębie centralnej grupy piezometrów (MW-25-MW-30) oraz w pobliżu Sherman Road (Rys. 5). W 1983 roku ze składowiska odpadów (na północ od terenu badawczego) usunięto znaczne ilości zanieczyszczeń, a dodatkowe wykopaliska przeprowadzono w 2001 roku. 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 piezometrów 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 przede wszystkim dawne piaski pogłębiarskie. Możliwa interferencja z opisaną metodą mogłaby wystąpić w przypadku obecności starożytnych skał węglanowych przy bardzo niskim pH wód gruntowych (<~5). Mogłoby to prowadzić do rozpuszczania węglanów i pojawienia się starożytnego sygnału w generowanym CO2. W tym obszarze nie są znane żadne znaczące pokłady CaCO3, niemniej jednak zmierzono stężenia 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ą (starożytne skały węglanowe mogłyby dostarczyć starożytnego CO2 w przypadku rozpuszczenia przez wody kwasowe). Zawartość Na+ była nieznacznie wyższa po południowej stronie terenu (najbliżej oceanu), ale żadne wartości nie mieściły się w zakresie wskazującym na znaczną intruzję wód morskich. Stężenia jonów wapnia wahały się od 8,0 do 58 mg L-1. Brak było wskazań na rozpuszczanie węglanów przy zestawieniu stężenia jonów wapnia z pH (r2 < 0,3). Wykresy bi-plot PCA nie wykazały silnych obciążeń dla żadnej zmiennej. Różnice między piezometrami również nie wskazywały na rozpuszczanie węglanów (Rys. 6). Taka analiza potwierdzająca powinna być uznana za krytyczną przy adaptacji metodologii do nowych terenów — szczególnie tych, w których geologia regionalna wskazuje na obecność znacznych formacji skał węglanowych.
Szybkość produkcji CO2 wynosiła od 0 do 34 mg CO2 d-1. Produkcja CO2 była najniższa w środkowej grupie otworów w obszarze, gdzie historyczne zanieczyszczenie było najwyższe (Rys. 5). Produkcja CO2 w otworze MW-01 (otwór tła — niepokazany, ale 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 powtórzonych analiz respiracji wynosiły od 0,03 do 6% CO2, a średnio mniej niż 1% (0,98). W celu dalszych obliczeń wyciągnięto średnią z pomiarów przeprowadzonych w dwóch 2-tygodniowych okresach pory suchej. 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 CH usuniętej w ciągu jednego miesiąca. Otwór tła (MW-01) miał wiek radiowęglowy 1280 lat przed teraźniejszością (ybp) lub 85 procent współczesnego węgla (pMC) — co mieści się w typowym zakresie dla starej materii organicznej w glebie26. Wartość z tego otworu została przyjęta jako tło dla izotopowego modelu mieszania. Ponownie, ze względu na to, że pobieranie próbek ograniczono do łącznie jednego miesiąca, dwa następujące po sobie okresy w tej samej porze roku zostały wykorzystane do „reprezentowania” pory suchej — uznawanej zazwyczaj za czas najbardziej stagnacyjnych warunków, a zatem konserwatywną 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). Otwory MW-27 i MW-32, w których podejrzewano nieszczelność pompy, wykazały współczesne wartości radiowęglowe, co potwierdziło ich zanieczyszczenie. Próbki te nie zostały uwzględnione w dalszych analizach.
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 San Diego (ID stacji 184) posłużyły do oszacowania tempa zasilania warstwy wodonośnej. Dane pływowe z tego samego okresu ze stacji NOAA 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 tempo poboru CO2. W parametryzacji modelu pomogły dodatkowe symulacje z zastosowaniem zmiennych średnich temp poboru CO2 i początkowego tła CO2, powiązane z 10% wzrostem gradientu hydraulicznego. Dodatkowa symulacja z wykorzystaniem średniego tempa poboru CO2 wykazała około 46% wzrost oszacowanego tła CO2 (i.e., wzrost z 6,5 do 9,5 g m-3), jeżeli tempo poboru zmieniło się z 0,00530 (+10%) do 0,00434 g hr-1 (-10%) w ciągu 2-tygodniowego okresu poboru (Tabela 3). Założenia modelu ZOI obejmowały pomijalną produkcję CO2 wynikającą z degradacji CH podczas okresu poboru oraz jednorodny początkowy rozkład CO2 w celu opracowania końcowej symulacji (Rys. 7). Szybkość reakcji CO2 dla obszaru badań mogła zostać niedoszacowana.
Wykorzystując tempo produkcji CO2, ilość CO2 przypisaną degradacji CH oraz szacunki z modelu ZOI, obliczono masę usuwanego CH w każdym otworze w jednostce czasu. Dane z Tabeli 1 wykorzystano w dwuskładnikowym modelu mieszania (równanie (1)) w celu wyznaczenia fpet dla każdego otworu. Ponieważ teren ten jest znany jedynie z zanieczyszczenia CH i nie stwierdzono żadnego innego źródła CO2 na terenie obiektu lub w jego pobliżu, przyjęto, że degradacja CH jest głównym źródłem CO2. Wartości fpet w obrębie terenu wahały się od 1 do 60% (Tabela 4). Proporcję tę przeliczono na bazę węglową i pomnożono przez tempo produkcji CO2, aby obliczyć tempo degradacji CH (Tabela 4). Wykorzystując objętość ZOI (Tabela 3), określono tempo degradacji zanieczyszczeń na jednostkę czasu i objętości (Tabela 4). Wartość ta mieściła się w zakresie od 0 do 32 mg C m-3 d-1 (Tabela 4). Degradacja CH była najniższa w obszarach o najwyższym historycznym zanieczyszczeniu CH (MW-25 - MW-30). Najwyższą degradację CH zmierzono w otworach znajdujących się w pobliżu obrzeży terenu (w okolicy Sherman Road). Produkcja CO2 była w tych obszarach wyższa, podczas gdy fpet wskazywał na znaczną wymianę CH (Rys. 8).

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

Rycina 2. Pułapki NaOH przygotowane do rozmieszczenia w terenie. Butelki na surowicę o pojemności 120 ml z dodaną pułapką NaOH, zamknięte za pomocą zakrętek zaciskowych.

Rycina 3. Układ terenowy. Przewody poprowadzone do wyposażonych studni (lewo), pułapka rozmieszczona przy studni (prawy górny róg) oraz system dystrybucji zasilania solarnego (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 nasadki otworów pokazujące linie recyrkulacji gazu. Na tej rycinie przedstawiono nasadki otworów zmodyfikowane poprzez dodanie linii wlotowych i powrotnych dla gazu.

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

Rysunek 6. Wykres biplot PCA wykazujący brak korelacji między kationami rozpuszczonymi a pH. Rysunek ten przedstawia biplot wyników (scores) i ładunków (loadings) PCA utworzonych na podstawie danych hydrogeologicznych (pH i kationów) dla terenu badań.

Rysunek 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 stężenie progowe 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. Rysunek ten przedstawia graficzną reprezentację modelu ZOI w 3 wymiarach.

Rycina 8. Szybkość degradacji zanieczyszczeń na jednostkę czasu i jednostkę powierzchni. Zmodyfikowano na podstawie 18. Jest to interpolowana szybkość degradacji CH na obszarze badań w pobranym okresie czasowym.

Wideo 1. Rozwój ZOI przy użyciu symulacji MT3DMS23 - MODFLOW (kliknij prawym przyciskiem myszy, aby pobrać). Pobierz, zainstaluj, zainicjuj i utwórz symulację dla ZOI.
| Otwór | δ13C
(‰VPDB) | Δ14C
(‰) | Wiek konwencjonalny
(ybp) | Procent nowoczesnego C
(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 konwersje izotopów CO2. Pomiary stabilnych izotopów i radiowęgla CO2 oraz konwersje 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 |
| Użyteczność 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 |
| Podłużna | m | 6,1 |
| Dyspersyjność |
| Pozioma poprzeczna | m | 0,61 |
| Dyspersyjność |
| Pionowa poprzeczna | m | 0,061 |
| Dyspersyjność |
| CO2 w gazach glebowych | % | 0,56 |
Tabela 2. Parametry modelu ZOI. Parametry wykorzystane w modelu ZOI i symulacjach.
| Poziom szybkości poboru | Szybkość pobierania | Stężenie tła | Wielkość strefy zahamowania wzrostu |
| Podłużne | Poprzeczny | Głębokość | Objętość |
| (g/h) | (g/m3) | (m) | | | (m3) |
| Maksymalny | 0.0131 | 17.6 | 2.47 | 0.77 | 0.13 | 0.193 |
| Średnia | 0.0048 | 6.5 | 2.28 | 0.72 | 0.12 | 0.176 |
| Minimum | 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łada się, że proces ten jest napędzany wyłącznie przez równowagę (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.