Gas Hydrate Formation and Dissipation Histories in the Northern Margin of Canada: Beaufort-Mackenzie and the Sverdrup Basins

Abstract

Gas hydrates (GHs) are a prominent subsurface feature on the Canadian Arctic continental margin. They occur both onshore and offshore, although they formed generally terrestrially, during the last glacial sea level low-stand, both in a region that was persistently glaciated (Queen Elizabeth Islands Group, Canadian Arctic Archipelago (QEIG)), and in a region that was not persistently glaciated (Mackenzie Delta-Beaufort Sea (MD-BS)). Parts of both regions were transgressed in the Holocene. We study the dynamic permafrost and GH history in both regions using a numerical model to illustrate how changes in setting and environment, especially periodic glacial ice cover, affected GH stability. MD-BS models represent the Mallik wellsite and these models successfully match current permafrost and GH bases observed in the well-studied Mallik wells. The MD-BS models show clearly that GHs have persisted through interglacial episodes. Lower surface temperatures in the more northerly QEIG result in an earlier appearance of GH stability that persists through glacial-interglacial intervals, although the base of GH base stability varies up to 0.2 km during the 100 ka cycles. Because of the persistent glacial ice cover QEIG models illustrate pressure effects attributed to regional ice sheet loading on the bases of both permafrost and GHs since 0.9 MYBP. QEIG model permafrost and GH depths are 572 m and 1072 m, respectively, which is like that observed commonly on well logs in the QEIG. In order to match the observed GH bases in the QEIG it is necessary to introduce ice buildup and thaw gradually during the glacials and interglacials. QEIG sea level rose 100–120 m about 10 ka ago following the most recent glaciation. Shorelines have risen subsequently due to isostatic glacial unloading. Detailed recent history modeling in QEIG coastal regions, where surface temperatures have changed from near zero in the offshore to −20°C in the onshore setting results in a model GH stability base, that is, <0.5 km. These coastal model results are significantly shallower than the inferred average GH base about 1 km in wells, Smith and Judge (1993). QEIG interisland channels are generally shallow and much of the previous shoreline inundated by the Holocene transgression was above the glacial sea level low-stand during the last ice age, resulting in a QEIG setting somewhat analogous to the relict terrestrial GH now transgressed by the shallow Beaufort Sea. It is also possible that the marine conditions were present at emergent shorelines for a shorter time or that the pretransgression subsurface temperatures persisted or were influenced by coastal settings, especially where lateral effects may not be well represented by 1D models.