Ice-bearing Permafrost Research Articles

In situ bottom sediment temperatures in the Siberian Arctic seas: Current state of subsea permafrost in the Kara sea vs laptev and East Siberian seas

Features of sediment temperature on the shelf and continental slope areas of the Russian Arctic seas and its physical properties are important for understanding the current state of subsea permafrost and the gas hydrates stability zone. New data are reported for the Kara Sea region where the bottom sediment temperatures are influenced by warming effects from great Siberian rivers and the Atlantic currents. The data collected during marine expeditions in 2019–2022 are combined with results of earlier marine studies, drilling operations, and geophysical surveys in the Laptev and East Siberian seas, in order to identify major trends of in situ temperature and properties distribution of bottom sediments in the Russian Arctic region.Most (85%) of bottom sediments in the Kara Sea shelf, as well as in the Laptev and East Siberian shelves, consist of water-saturated silty clay and silt with rather uniform particle size distribution. The obtained thermal conductivity and heat capacity values for the Kara Sea sediments agree with the values of 1.0 W/(m·K) and 2900 kJ/m3, respectively, obtained previously from other Arctic seas. Thermal conductivity becomes up to 40% higher depthward from 0 to 2 m subbottom depth, possibly, because of lower moisture content and porosity in more lithified sediments.The bottom sediment temperatures in the Arctic seas are distributed unevenly, especially in the Kara Sea shelf (from +5.0 °C in the west to −1.4 °C in the east), where the high sediment temperatures in the western and central parts of the Kara Sea being due to the effect of warm water inputs. The distribution of bottom sediment temperatures correlates well with distribution of relic subsea permafrost. Ice-bearing permafrost in the Siberian Arctic shelf extends from the shoreline till sea depths of 80–100 m, within the respective offshore distances of ∼800–1000 km in some areas, but permafrost remnants may exist locally at sea depths within 120 m. Buried 100–600 m thick continuous subsea permafrost may occur in the Kara, Laptev, and East Siberian shelves under unfrozen (cryotic) saline shallow sediments. However, subsea permafrost is discontinuous and sporadic at sea depths ∼70 m and more. Thus, the bottom sediment temperature features in the Arctic seas can be used as a proxy of subsea permafrost extent contenting intrapermafrost and subpermafrost gas and gas hydrate accumulations.

Geologic controls on the genesis of the Arctic permafrost and sub-permafrost methane hydrate-bearing system in the Beaufort–Mackenzie Delta

The Canadian Mackenzie Delta exhibits a high volume of proven sub-permafrost gas hydrates that naturally trap a significant amount of deep-sourced thermogenic methane (CH4) at the Mallik site. The present study aims to validate the proposed Arctic sub-permafrost gas hydrate formation mechanism, implying that CH4-rich fluids were vertically transported from deep overpressurized zones via geologic fault systems and formed the present-day observed GH deposit since the Late Pleistocene. Given this hypothesis, the coastal permafrost began to form since the early Pleistocene sea-level retreat, steadily increasing in thickness until 1 Million years (Ma) ago. Data from well logs and 2D seismic profiles were digitized to establish the first field-scale static geologic 3D model of the Mallik site, and to comprehensively study the genesis of the permafrost and its associated GH system. The implemented 3D model considers the spatially heterogeneously distributed hydraulic properties of the individual lithologies at the Mallik site. Simulations using a proven thermo-hydro-chemical numerical framework were employed to gain insights into the hydrogeologic role of the regional fault systems in view of the CH4-rich fluid migration and the geologic controls on the spatial extent of the sub-permafrost GH accumulations during the past 1 Ma. For >87% of the Mallik well sections, the predicted permafrost thickness deviates from the observations by less than 0.8%, which validates the general model implementation. The simulated ice-bearing permafrost and GH interval thicknesses as well as sub-permafrost temperature profiles are consistent with the respective field observations, confirming our introduced hypothesis. The spatial distribution of GHs is a result of the comprehensive interaction between various processes, including the source-gas generation rate, subsurface temperature, and the hydraulic properties of the structural geologic features. Overall, the good agreement between simulations and observations demonstrates that the present study provides a valid representation of the geologic controls driving the complex permafrost-GH system. The model’s applicability for the prediction of GH deposits in permafrost settings in terms of their thicknesses and saturations can provide relevant contributions to future GH exploration and exploitation.

Exploring the capabilities of electrical resistivity tomography to study subsea permafrost

Abstract. Sea level rise and coastal erosion have inundated large areas of Arctic permafrost. Submergence by warm and saline waters increases the rate of inundated permafrost thaw compared to sub-aerial thawing on land. Studying the contact between the unfrozen and frozen sediments below the seabed, also known as the ice-bearing permafrost table (IBPT), provides valuable information to understand the evolution of sub-aquatic permafrost, which is key to improving and understanding coastal erosion prediction models and potential greenhouse gas emissions. In this study, we use data from 2D electrical resistivity tomography (ERT) collected in the nearshore coastal zone of two Arctic regions that differ in their environmental conditions (e.g., seawater depth and resistivity) to image and study the subsea permafrost. The inversion of 2D ERT data sets is commonly performed using deterministic approaches that favor smoothed solutions, which are typically interpreted using a user-specified resistivity threshold to identify the IBPT position. In contrast, to target the IBPT position directly during inversion, we use a layer-based model parameterization and a global optimization approach to invert our ERT data. This approach results in ensembles of layered 2D model solutions, which we use to identify the IBPT and estimate the resistivity of the unfrozen and frozen sediments, including estimates of uncertainties. Additionally, we globally invert 1D synthetic resistivity data and perform sensitivity analyses to study, in a simpler way, the correlations and influences of our model parameters. The set of methods provided in this study may help to further exploit ERT data collected in such permafrost environments as well as for the design of future field experiments.

Numerical Simulation of Coastal Sub-Permafrost Gas Hydrate Formation in the Mackenzie Delta, Canadian Arctic

The Mackenzie Delta (MD) is a permafrost-bearing region along the coasts of the Canadian Arctic which exhibits high sub-permafrost gas hydrate (GH) reserves. The GH occurring at the Mallik site in the MD is dominated by thermogenic methane (CH4), which migrated from deep conventional hydrocarbon reservoirs, very likely through the present fault systems. Therefore, it is assumed that fluid flow transports dissolved CH4 upward and out of the deeper overpressurized reservoirs via the existing polygonal fault system and then forms the GH accumulations in the Kugmallit–Mackenzie Bay Sequences. We investigate the feasibility of this mechanism with a thermo–hydraulic–chemical numerical model, representing a cross section of the Mallik site. We present the first simulations that consider permafrost formation and thawing, as well as the formation of GH accumulations sourced from the upward migrating CH4-rich formation fluid. The simulation results show that temperature distribution, as well as the thickness and base of the ice-bearing permafrost are consistent with corresponding field observations. The primary driver for the spatial GH distribution is the permeability of the host sediments. Thus, the hypothesis on GH formation by dissolved CH4 originating from deeper geological reservoirs is successfully validated. Furthermore, our results demonstrate that the permafrost has been substantially heated to 0.8–1.3 °C, triggered by the global temperature increase of about 0.44 °C and further enhanced by the Arctic Amplification effect at the Mallik site from the early 1970s to the mid-2000s.

Thermokarst Lagoons: A Core-Based Assessment of Depositional Characteristics and an Estimate of Carbon Pools on the Bykovsky Peninsula

Permafrost region subsurface organic carbon (OC) pools are a major component of the terrestrial carbon cycle and vulnerable to a warming climate. Thermokarst lagoons are an important transition stage with complex depositional histories during which permafrost and lacustrine carbon pools are transformed along eroding Arctic coasts. The effects of temperature and salinity changes during thermokarst lake to lagoon transitions on thaw history and lagoon deposits are understudied. We analyzed two 30-m-long sediment cores from two thermokarst lagoons on the Bykovsky Peninsula, Northeast Siberia, using sedimentological, geochronological, hydrochemical, and biogeochemical techniques. Using remote sensing we distinguished between a semi-closed and a nearly closed lagoon. We (1) characterized the depositional history, (2) studied the impact of marine inundation on ice-bearing permafrost and taliks, and (3) quantified the OC pools for different stages of thermokarst lagoons. Fluvial and former Yedoma deposits were found at depth between 30 and 8.5 m, while lake and lagoon deposits formed the upper layers. The electrical conductivity of the pore water indicated hypersaline conditions for the semi-closed lagoon (max: 108 mS/cm), while fresh to brackish conditions were observed beneath a 5 m-thick surface saline layer at the nearly closed lagoon. The deposits had a mean OC content of 15 ± 2 kg/m3, with higher values in the semi-closed lagoon. Based on the cores we estimated a total OC pool of 5.7 Mt-C for the first 30 m of sediment below five mapped lagoons on the Bykovsky Peninsula. Our results suggest that paleo river branches shaped the middle Pleistocene landscape followed by late Pleistocene Yedoma permafrost accumulation and early Holocene lake development. Afterward, lake drainage, marine flooding, and bedfast ice formation caused the saline enrichment of pore water, which led to cryotic talik development. We find that the OC-pool of Arctic lagoons may comprise a substantial inventory of partially thawed and partially refrozen OC, which is available for microbial degradation processes at the Arctic terrestrial-marine interface. Climate change in the Arctic leading to sea level rise, permafrost thaw, coastal erosion, and sea ice loss may increase the rate of thermokarst lagoon formation and thus increase the importance of lagoons as biogeochemical processors of former permafrost OC.

Spatial distribution and controls of permafrost development in non-glacial Arctic catchment over the Holocene, Fuglebekken, SW Spitsbergen

This article presents the distribution and properties of the permafrost based on electrical resistivity tomography (ERT) and multichannel analysis of surface waves (MASW) data collected at the Fuglebekken coastal catchment area in SW Spitsbergen. This work summarizes the development of permafrost in this area during the Holocene, from the mountain environment through to the system of elevated marine terraces found around the coast. The ERT models were analysed taking into consideration the non-unique nature of the data inversion process and the physical limitations of this method. Comparing the ERT and the MASW results allows a zonal characterization of the occurring ice-bearing permafrost and its correlation with the evolution history of the catchment area. Maritime transgression as well as intensive watercourses during past degradation episodes have altered the permafrost presence and ice-accumulating abilities of different sediment zones. Permafrost development depends greatly on the presence of surface watercourses in talus slopes. The youngest elevated uplifted marine terrace did not develop an ice-rich permafrost, but the presence of permafrost in a cryotic form is possible. The significant range of the fjord water infiltration found within the sedimentary cover have influenced the development of the coastal permafrost. The current structure of ice-bearing permafrost found in the research area seems to be very sensitive to the climatic changes. Based on these results, we propose a model for the formation of the current permafrost in the studied area.

Heat and Salt Flow in Subsea Permafrost Modeled with CryoGRID2.

Thawing of subsea permafrost can impact offshore infrastructure, affect coastal erosion, and release permafrost organic matter. Thawing is usually modeled as the result of heat transfer, although salt diffusion may play an important role in marine settings. To better quantify nearshore subsea permafrost thawing, we applied the CryoGRID2 heat diffusion model and coupled it to a salt diffusion model. We simulated coastline retreat and subsea permafrost evolution as it develops through successive stages of a thawing sequence at the Bykovsky Peninsula, Siberia. Sensitivity analyses for seawater salinity were performed to compare the results for the Bykovsky Peninsula with those of typical Arctic seawater. For the Bykovsky Peninsula, the modeled ice‐bearing permafrost table (IBPT) for ice‐rich sand and an erosion rate of 0.25 m/year was 16.7 m below the seabed 350 m offshore. The model outputs were compared to the IBPT depth estimated from coastline retreat and electrical resistivity surveys perpendicular to and crossing the shoreline of the Bykovsky Peninsula. The interpreted geoelectric data suggest that the IBPT dipped to 15–20 m below the seabed at 350 m offshore. Both results suggest that cold saline water forms beneath grounded ice and floating sea ice in shallow water, causing cryotic benthic temperatures. The freezing point depression produced by salt diffusion can delay or prevent ice formation in the sediment and enhance the IBPT degradation rate. Therefore, salt diffusion may facilitate the release of greenhouse gasses to the atmosphere and considerably affect the design of offshore and coastal infrastructure in subsea permafrost areas.

Remote sensing evidence of lava–ground ice interactions associated with the Lost Jim Lava Flow, Seward Peninsula, Alaska

Thermokarst terrains develop when ice-bearing permafrost melts and causes the overlying surface to subside or collapse. This process occurs widely throughout Arctic regions due to environmental and climatological factors, but can also be induced by localized melting of ground ice by active lava flows. The Lost Jim Lava Flow (LJLF) on the Seward Peninsula of Alaska provides evidence of former lava–ground ice interactions. Associated geomorphic features, on the scale of meters to tens of meters, were identified using satellite orthoimages and stereo-derived digital terrain models. The flow exhibits positive- and mixed-relief features, including tumuli (N = 26) and shatter rings (N = 4), as well as negative-relief features, such as lava tube skylights (N = 100) and irregularly shaped topographic depressions (N = 1188) that are interpreted to include lava-rise pits and lava-induced thermokarst terrain. Along the margins of the flow, there are also clusters of small peripheral pits that may be the products of meltwater or steam escape. On Mars, we observed morphologically similar pits near lava flow margins in northeastern Elysium Planitia, which suggests a common formation mechanism. Investigating the LJLF may therefore help to elucidate processes of lava–ground ice interaction on both Earth and Mars.

Distribution and origin of ground ice in University Valley, McMurdo Dry Valleys, Antarctica

AbstractGround ice is one of the most important and dynamic geologic components of permafrost; however, few studies have investigated the distribution and origin of ground ice in the McMurdo Dry Valleys of Antarctica. In this study, ice-bearing permafrost cores were collected from 18 sites in University Valley, a small hanging glacial valley in the Quartermain Mountains. Ground ice was found to be ubiquitous in the upper 2 m of permafrost soils, with excess ice contents reaching 93%, but ground ice conditions were not homogeneous. Ground ice content was variable within polygons and along the valley floor, decreasing in the centres of polygons and increasing in the shoulders of polygons towards the mouth of the valley. Ground ice also had different origins: vapour deposition, freezing of partially evaporated snow meltwater and buried glacier ice. The variability in the distribution and origin of ground ice can be attributed to ground surface temperature and moisture conditions, which separate the valley into distinct zones. Ground ice of vapour-deposition origin was predominantly situated in perennially cryotic zones, whereas ground ice formed by the freezing of evaporated snow meltwater was predominantly found in seasonally non-cryotic zones.

Subsea ice‐bearing permafrost on the U.S. Beaufort Margin: 1. Minimum seaward extent defined from multichannel seismic reflection data

AbstractSubsea ice‐bearing permafrost (IBPF) and associated gas hydrate in the Arctic have been subject to a warming climate and saline intrusion since the last transgression at the end of the Pleistocene. The consequent degradation of IBPF is potentially associated with significant degassing of dissociating gas hydrate deposits. Previous studies interpreted the distribution of subsea permafrost on the U.S. Beaufort continental shelf based on geographically sparse data sets and modeling of expected thermal history. The most cited work projects subsea permafrost to the shelf edge (∼100 m isobath). This study uses a compilation of stacking velocity analyses from ∼100,000 line‐km of industry‐collected multichannel seismic reflection data acquired over 57,000 km2 of the U.S. Beaufort shelf to delineate continuous subsea IBPF. Gridded average velocities of the uppermost 750 ms two‐way travel time range from 1475 to 3110 m s−1. The monotonic, cross‐shore pattern in velocity distribution suggests that the seaward extent of continuous IBPF is within 37 km of the modern shoreline at water depths < 25 m. These interpretations corroborate recent Beaufort seismic refraction studies and provide the best, margin‐scale evidence that continuous subsea IBPF does not currently extend to the northern limits of the continental shelf.

Creep behaviors of methane hydrate coexisting with ice

A series of creep tests were carried out in order to clarify the creep behaviors of methane hydrate coexisting with ice, which may exist in the ice-bearing permafrost. Measurements on synthetic methane hydrate-ice specimens indicate that (1) higher external load will increase the initial strain and the axial strain during the whole creep test. And the creep strain rate will increase with the increase of external load; (2) the initial strain, axial strains and creep strain rate of the specimens under certain external load increase with the decrease of confining pressure; (3) temperature drop leads to the decrease of initial strain and axial strains of the specimens under certain external load. And the creep strain rate firstly increases and then decreases with the increase of temperature, due to the refreezing of dissociated water from methane hydrate; (4) the specimens run to the secondary creep stage earlier when they endure lower external load, higher confining pressure and lower temperature; (5) the specimens under higher confining pressure and lower temperature show higher failure strength.

Inferred gas hydrate and permafrost stability history models linked to climate change in the Beaufort-Mackenzie Basin, Arctic Canada

Abstract. Atmospheric methane from episodic gas hydrate (GH) destabilization, the "clathrate gun" hypothesis, is proposed to affect past climates, possibly since the Phanerozoic began or earlier. In the terrestrial Beaufort-Mackenzie Basin (BMB), GHs occur commonly below thick ice-bearing permafrost (IBP), but they are rare within it. Two end-member GH models, where gas is either trapped conventionally (Case 1) or where it is trapped dynamically by GH formation (Case 2), were simulated using profile (1-D) models and a 14 Myr ground surface temperature (GST) history based on marine isotopic data, adjusted to the study setting, constrained by deep heat flow, sedimentary succession conductivity, and observed IBP and Type I GH contacts in Mallik wells. Models consider latent heat effects throughout the IBP and GH intervals. Case 1 GHs formed at ~0.9 km depth only ~1 Myr ago by in situ transformation of conventionally trapped natural gas. Case 2 GHs begin to form at ~290–300 m ~6 Myr ago in the absence of lithological migration barriers. During glacial intervals Case 2 GH layers expand both downward and upward as the permafrost grows downward through and intercalated with GHs. The distinctive model results suggest that most BMB GHs resemble Case 1 models, based on the observed distinct and separate occurrences of GHs and IBP and the lack of observed GH intercalations in IBP. Case 2 GHs formed >255 m, below a persistent ice-filled permafrost layer that is as effective a seal to upward methane migration as are Case 1 lithological seals. All models respond to GST variations, but in a delayed and muted manner such that GH layers continue to grow even as the GST begins to increase. The models show that the GH stability zone history is buffered strongly by IBP during the interglacials. Thick IBP and GHs could have persisted since ~1.0 Myr ago and ~4.0 Myr ago for Cases 1 and 2, respectively. Offshore BMB IBP and GHs formed terrestrially during Pleistocene sea level low stands. Where IBP is sufficiently thick, both IBP and GHs persist even where inundated by a Holocene sea level rise and both are also expected to persist into the next glacial even if atmospheric CO2 doubles. We do not address the "clathrate gun" hypothesis directly, but our models show that sub-IBP GHs respond to, rather than cause GST changes, due to both how GST changes propagates with depth and latent heat effects. Models show that many thick GH accumulations are prevented from contributing methane to the atmosphere, because they are almost certainly trapped below either ice-filled IBP or lithological barriers. Where permafrost is sufficiently thick, combinations of geological structure, thermal processes and material properties make sub-IBP GHs unlikely sources for significant atmospheric methane fluxes. Our sub-IBP GH model histories suggest that similar models applied to other GH settings could improve the understanding of GHs and their potential to affect climate.

Mount Elbert Gas Hydrate Stratigraphic Test Well, Alaska North Slope: Coring operations, core sedimentology, and lithostratigraphy

In February 2007, BP Exploration (Alaska), the U.S. Department of Energy, and the U.S. Geological Survey completed the BPXA-DOE-USGS Mount Elbert Gas Hydrate Stratigraphic Test Well (Mount Elbert well) in the Milne Point Unit on the Alaska North Slope. The program achieved its primary goals of validating the pre-drill estimates of gas hydrate occurrence and thickness based on 3-D seismic interpretations and wireline log correlations and collecting a comprehensive suite of logging, coring, and pressure testing data. The upper section of the Mount Elbert well was drilled through the base of ice-bearing permafrost to a casing point of 594 m (1950 ft), approximately 15 m (50 ft) above the top of the targeted reservoir interval. The lower portion of the well was continuously cored from 606 m (1987 ft) to 760 m (2494 ft) and drilled to a total depth of 914 m. Ice-bearing permafrost extends to a depth of roughly 536 m and the base of gas hydrate stability is interpreted to extend to a depth of 870 m. Coring through the targeted gas hydrate bearing reservoirs was completed using a wireline-retrievable system. The coring program achieved 85% recovery of 7.6 cm (3 in) diameter core through 154 m (504 ft) of the hole. An onsite team processed the cores, collecting and preserving approximately 250 sub-samples for analyses of pore water geochemistry, microbiology, gas chemistry, petrophysical analysis, and thermal and physical properties. Eleven samples were immediately transferred to either methane-charged pressure vessels or liquid nitrogen for future study of the preserved gas hydrate. Additional offsite sampling, analyses, and detailed description of the cores were also conducted. Based on this work, one lithostratigraphic unit with eight subunits was identified across the cored interval. Subunits II and Va comprise the majority of the reservoir facies and are dominantly very fine to fine, moderately sorted, quartz, feldspar, and lithic fragment-bearing to -rich sands. Lithostratigraphic and palynologic data indicate that this section is most likely early Eocene to late Paleocene in age. The examined units contain evidence for both marine and non-marine lithofacies, and indications that the depositional environment for the reservoir facies may have been shallower marine than originally interpreted based on pre-drill wireline log interpretations. There is also evidence of reduced salinity marine conditions during deposition that may be related to the paleo-climate and depositional conditions during the early Eocene.

Pore fluid geochemistry from the Mount Elbert Gas Hydrate Stratigraphic Test Well, Alaska North Slope

The BPXA-DOE-USGS Mount Elbert Gas Hydrate Stratigraphic Test Well was drilled and cored from 606.5 to 760.1 m on the North Slope of Alaska, to evaluate the occurrence, distribution and formation of gas hydrate in sediments below the base of the ice-bearing permafrost. Both the dissolved chloride and the isotopic composition of the water co-vary in the gas hydrate-bearing zones, consistent with gas hydrate dissociation during core recovery, and they provide independent indicators to constrain the zone of gas hydrate occurrence. Analyses of chloride and water isotope data indicate that an observed increase in salinity towards the top of the cored section reflects the presence of residual fluids from ion exclusion during ice formation at the base of the permafrost layer. These salinity changes are the main factor controlling major and minor ion distributions in the Mount Elbert Well. The resulting background chloride can be simulated with a one-dimensional diffusion model, and the results suggest that the ion exclusion at the top of the cored section reflects deepening of the permafrost layer following the last glaciation (∼100 kyr), consistent with published thermal models. Gas hydrate saturation values estimated from dissolved chloride agree with estimates based on logging data when the gas hydrate occupies more than 20% of the pore space; the correlation is less robust at lower saturation values. The highest gas hydrate concentrations at the Mount Elbert Well are clearly associated with coarse-grained sedimentary sections, as expected from theoretical calculations and field observations in marine and other arctic sediment cores.

Natural Gas Hydrates in the Offshore Beaufort–Mackenzie Basin—Study of a Feasible Energy Source II

In the offshore part of Beaufort–Mackenzie Basin depth of methane hydrate stability reaches more than 1.5 km. However, there are areas in the western part of the basin where there are no conditions of methane hydrate stability. Construction of the first contour maps displaying thickness of hydrate stability zones as well as hydrate stability zone thicknesses below permafrost in the offshore area, shows that these zones can reach 1200 m and 900 m, respectively. Depth to the base of ice-bearing relict permafrost under the sea (depth of the −1°C isotherm-ice-bearing permafrost base) and regional variations of geothermal gradient are the main controlling factors. Hydrostatic pressures in the upper 1500 m are the rule. History of methane hydrate stability zone is related mainly to the history of permafrost and it reached maximum depth in early Holocene. More recently, the permafrost and hydrate zone is diminishing because of sea transgression. Reevaluation of the location of possible gas hydrate occurrences is done from the analysis of well logs and other indicators in conjunction with knowledge of the hydrate stability zone. In the offshore Beaufort–Mackenzie Basin, methane hydrate occurs in 21 wells. Nine of these locations coincides with underlying conventional hydrocarbon occurrences. Previous analyses place some of the hydrate occurrences at greater depths than proposed for the methane hydrate-stability zone described in this study. Interpretation of geological cross sections and maps of geological sequences reveals that hydrates are occurring in the Iperk–Kugmallit sequence. Hydrate–gas contact zones, however, are possible in numerous situations. As there are no significant geological seals in the deeper part of the offshore basin (all hydrates are within Iperk), it is suggested that overlying permafrost and hydrate stability zone acted as the only trap for upward migrating gas during the last tens of thousand of years (i.e., Sangamonian to Holocene).

Stability Zone of Natural Gas Hydrates in a Permafrost-Bearing Region of the Beaufort—Mackenzie Basin: Study of a Feasible Energy Source1 (Geological Survey of Canada Contribution No.1999275)

Analysis of geological and geophysical data from 150 wells in the Beaufort—Mackenzie region(study area between 68°30′–70°00′N and 131°–39°W) led to reinterpretation of the depth ofmethane hydrate stability and construction of the first contour maps displaying thickness of hydratestability zones as well as hydrate stability zone thicknesses below permafrost. Calculations werebased on construction of temperature-depth profiles incorporating regional heat-flow values, temperatureat the base of ice-bearing permafrost, and models relating thermal conductivity with depth.Data analysis indicates the presence and extent of the methane hydrate stability zone is relatedmainly to the history of permafrost development and less so by the relatively small regionalvariations of temperature gradients. Analysis of well logs and other indicators in conjunction withknowledge of the hydrate stability zone allows reevaluation of the location of possible gas hydrateoccurrences. Log analysis indicates that in the onshore and shallow sea area of theBeaufort—Mackenzie Basin, methane hydrate occurs in 27 wells. Fifteen of these locations coincides withunderlying conventional hydrocarbon occurrences. Previous analyses place some of the hydrateoccurrences at greater depths than proposed for the methane hydrate stability zone described inthis study. Interpretation of geological cross sections reveals that hydratesare related mainly to sandy deltaic and delta-plain deposits in Iperk, Kugmallit, and Reindeer sequences althoughadditional hydrate picks have been inferred in other sequences, such as Richards. Overlyingpermafrost may act as seal for hydrate accumulations; however, the thickness of permafrost andits related hydrate stability zone fluctuated during geological time. It is interpreted that only inthe last tens of thousand of years (i.e., Sangamonian to Holocene), conditions for hydrates changedfrom nonstable to stable. During Early and Late Wisconsinan and Holocene time, conditions werefavorable for generation and trapping of hydrates. However, previously during Sangamonian time,less favorable conditions existed for hydrate stability. Gas release from hydrates may have occurredduring times when hydrate stability was nonexistent because of permafrost melting episodes. It isinterpreted that entrapment of gas in hydrate molecular structures is related to the existence ofconventional structural traps as well as less permeable sediments such as the Mackenzie BayFormation, which act as seal.

Subsurface temperatures and geothermal gradients on the north slope of Alaska

On the North Slope of Alaska, geothermal gradient data are available from high-resolution, equilibrated well-bore surveys and from estimates based on well-log identification of the base of ice-bearing permafrost. A total of 46 North Slope wells, considered to be in or near thermal equilibrium, have been surveyed with high-resolution temperatures devices and geothermal gradients can be interpreted directly from these recorded temperature profiles. To augment the limited North Slope temperature data base, a new method of evaluating local geothermal gradients has been developed. In this method, a series of well-log picks for the base of the ice-bearing permafrost from 102 wells have been used, along with regional temperature constants derived from the high-resolution stabilized well-bore temperature surveys, to project geothermal gradients. Geothermal gradients calculated from the high-resolution temperature surveys generally agree with those projected from known ice-bearing permafrost depths over most of the North Slope. Values in the ice-bearing permafrost range from ≈ 1.5° C 100 m in the Prudhoe Bay area to ≈ 4.5δ C 100 m in the east-central portion of the National Petroleum Reserve in Alaska. Geothermal gradients below the ice-bearing permafrost sequence range from ≈ 1.6° C 100 m to ≈ 5.2° C 100 m .

Surface-temperature history during the last 1000 years near Prudhoe Bay, Alaska: applying control theory to the inversion of borehole temperature profiles

The temperature-depth profile data by Lachenbruch et al. (1982) in ice-bearing permafrost near Prudhoe Bay, Alaska, are analyzed with control theory to retrodict the last 1000-yr history of the surface temperature. For this purpose, an inversion scheme for a nonlinear problem of heat diffusion in permafrost is developed. The following conclusions are made: (1) During the last 1000 years, the surface temperature near Prudhoe Bay, Alaska, oscillated with amplitude of 3 ± 1°C. (2) There was a minimum of the temperature in 1853 AD between two maxima at 1623 AD and 1937 AD. (3) A net warming by 4°C occurred from 1853 AD to 1973 AD.

Surface-temperature history during the last 1000 years near Prudhoe Bay, Alaska: applying control theory to the inversion of borehole temperature profiles

The temperature-depth profile data by Lachenbruch et al. (1982) in ice-bearing permafrost near Prudhoe Bay, Alaska, are analyzed with control theory to retrodict the last 1000-yr history of the surface temperature. For this purpose, an inversion scheme for a nonlinear problem of heat diffusion in permafrost is developed. The following conclusions are made: (1) During the last 1000 years, the surface temperature near Prudhoe Bay, Alaska, oscillated with amplitude of 3 ± 1°C. (2) There was a minimum of the temperature in 1853 AD between two maxima at 1623 AD and 1937 AD. (3) A net warming by 4°C occurred from 1853 AD to 1973 AD.

A constraint to the Wisconsinan glacial history, Canadian Arctic Archipelago

AbstractThick ice‐bearing permafrost is not observed today beneath the deeper channels of the central Queen Elizabeth Islands, Canadian Arctic Archipelago. Analysis of a precision temperature log recently obtained at an offshore well near Ellef Ringnes Island indicates that the thermal regime beneath the seabed is in equilibrium with today's marine environment. If thick permafrost similar to that observed on land today had existed in the Pleistocene in areas that are presently offshore, then such permafrost must have started melting no later than 25000 years ago in order to allow the present thermal regime to evolve. This suggests that the inter‐island channels must have been water‐filled at least by that date.