Ice Thickness Measurements Research Articles

Subglacial Topography of an Icefall Inferred From Repeated Terrestrial Laser Scanning

Information about ice thickness and subglacial topography is a key element for ice dynamics modeling and, consequently, for better understanding of current response of glaciers to climatic forcing. While there has been a substantial progress in the measurements of glacial surface, ice thickness measurements still suffer from limited coverage and strong interpolation. Furthermore, there are some glaciated areas where typical remote sensing techniques provide only limited data due to unfavorable climatic conditions and steep topography. A perfect example of such a feature is icefalls dominating landscapes of the maritime Antarctic. In order to close that gap, this letter presents a simple method of inferring the subglacial topography based on repeated terrestrial laser scanning of the ice surface and inverse shallow ice approximation modeling. Emerald Icefalls, King George Island, were surveyed twice within eight-day period, allowing to derive ice surface velocity field by feature tracking analysis and to estimate ice thickness and the subglacial topography. Calculated ice thickness is low, with the mean value of 72±18 m. Therefore, in line with former studies, the overall ice flux of Emerald Icefalls is small despite their relatively high surface ice flow velocities.

Circum-Greenland, ice-thickness measurements collected during PROMICE airborne surveys in 2007, 2011 and 2015

The Greenland ice sheet has experienced an average mass loss of 142 ± 49 Gt/yr from 1992 to 2011 (Shepherd et al. 2012), making it a significant contributor to sea-level rise. Part of the ice- sheet mass loss is the result of increased dynamic response of outlet glaciers (Rignot et al. 2011). The ice discharge from outlet glaciers can be quantified by coincident measurements of ice velocity and ice thickness (Thomas et al. 2000; van den Broeke et al. 2016). As part of the Programme for monitoring of the Greenland Ice Sheet (PROMICE; Ahlstrøm et al. 2008), three airborne surveys were carried out in 2007, 2011 and 2015, with the aim of measuring the changes in Greenland ice-sheet thicknesses. The purpose of the airborne surveys was to collect data to assess the dynamic mass loss of the Greenland ice sheet (Andersen et al. 2015). Here, we present these datasets of observations from ice-penetrating radar and airborne laser scanning, which, in combination, make us able to determine the ice thickness precisely. Surface-elevation changes between surveys are also presented, although we do not provide an in-depth scientific interpretation of these.

The Greenland ice sheet – snowline elevations at the end of the melt seasons from 2000 to 2017

The Greenland ice sheet has experienced an average mass loss of 142 ± 49 Gt/yr from 1992 to 2011 (Shepherd et al. 2012), making it a significant contributor to sea-level rise. Part of the ice- sheet mass loss is the result of increased dynamic response of outlet glaciers (Rignot et al. 2011). The ice discharge from outlet glaciers can be quantified by coincident measurements of ice velocity and ice thickness (Thomas et al. 2000; van den Broeke et al. 2016). As part of the Programme for monitoring of the Greenland Ice Sheet (PROMICE; Ahlstrøm et al. 2008), three airborne surveys were carried out in 2007, 2011 and 2015, with the aim of measuring the changes in Greenland ice-sheet thicknesses. The purpose of the airborne surveys was to collect data to assess the dynamic mass loss of the Greenland ice sheet (Andersen et al. 2015). Here, we present these datasets of observations from ice-penetrating radar and airborne laser scanning, which, in combination, make us able to determine the ice thickness precisely. Surface-elevation changes between surveys are also presented, although we do not provide an in-depth scientific interpretation of these.

Estimation of Glacier Thickness From Surface Mass Balance and Ice Flow Velocities: A Case Study on Argentière Glacier, France

Glacier thickness distribution is a prerequisite to simulate the future of glaciers. Inaccurate thicknesses may lead to significant uncertainties in the timing of future changes to glaciers and their consequences for water resources or sea level rise. Unfortunately, glacier thickness distribution is rarely measured and consequently has to be estimated. In this study, we present an approach developed on the well documented Argentiere Glacier (French Alps) that uses surface mass balance (SMB) together with surface flow velocity data to quantify glacier thickness distribution over the entire surface of the glacier. We compare the results of our approach to those obtained applying Farinotti et al. (2009) approach. Our results show that glacier thickness distribution are significantly biased when the glacier SMB profile used to quantify the ice fluxes is not constrained with in situ measurements. We also show that even with SMB measurements available on the studied glacier, ice flux estimates can be inaccurate. This inability to correctly estimate ice fluxes from the apparent SMB may be due to the steady state assumption that is not respected from the available glacier surface topography data. Therefore, ice thickness measurements on few cross sections (four are used in this study) are required to constrain the ice flux estimates and lead to an overall agreement between the ice thickness estimations and measurements. Using our approach, the ice thicknesses only differ by 10% from observations in average, but can differ by up to 150 m (or 30%) locally. We also show that approaches that use the glacier surface slope can lead to large uncertainties given that the quantification of the slope is highly uncertain. The approach presented here does not pretend to be applied globally but rather as a tool to quantify ice thickness distribution over the entire surface of glaciers for which a few in situ surface mass balance and thickness data are available together with surface flow velocities that can be obtained for example from remote sensing.

A method for correcting seal-borne oceanographic data and application to the estimation of regional sea ice thickness

Abstract The high-latitude oceans surrounding Antarctica are substantially undersampled compared to lower latitudes. Mammal based instruments such as Conductivity-Temperature-Depth Satellite Relay Data Loggers (CTD-SRDLs) present one possible solution. Unfortunately, these are subject to instrument-dependent offsets in absolute salinity. This study investigates a set of satellite-transmitted data collected by CTD-SRDLs mounted on Weddell seals (Leptonychotes weddellii) in the South-western Ross Sea in 2011. The uncorrected salinity offset between devices was found to be up to 1.4 g kg−1, making the data unsuitable for some oceanographic studies without correction. Here, a correction method was developed that uses profiles from pairs of CTD-SRDLs that are considered to be co-located and to sample the same body of water if they occur within defined time and space windows. Using least squares, a best-fit solution to the matrix of offsets in co-located pairs was found that reduces salinity offsets between the CTD-SRDLs. These offsets are smaller than the original offsets by a factor of 10. A calibrated reference instrument, that was co-located with some of the devices, provided further improvement in the absolute accuracy of all the CTD-SRDLs. Using the corrected CTD-SRDL data we estimate the rejection of salt into the water column by sea ice formation, and derived the time evolution of sea ice thickness in the South-western Ross Sea. Our estimates of regional sea ice thickness are in agreement with direct sea ice thickness measurements taken over a limited area in November 2011, providing further affirmation of our method.

Marginal Ice Zone Thickness and Extent due to Wave Radiation Stress

AbstractOcean surface wave radiation stress represents the flux of momentum due to the waves. When waves are dissipated or reflected by sea ice, that momentum is absorbed or reflected, resulting in a horizontal forcing that frequently compresses the ice. In this work, wave radiation stress is used to estimate the compressive force applied by waves to the marginal ice zone (MIZ). It is balanced by an ice internal compressive stress based on Mohr–Coulomb granular materials theory. The ice internal stress can be related to ice thickness, allowing this force balance to be used as a model for the estimation of MIZ ice thickness. The model was validated and tested using data collected during two field campaigns in the St. Lawrence estuary in 2016 and 2017. Modeled ice thickness was found to be consistent with the mean measured ice thickness over the conditions available. The range of validity of the model is discussed, and a definition of MIZ extent, based on the relative strength of wind and wave forcing, is proposed.

Routine determination of ice thickness for cryo-EM grids.

Recent advances in instrumentation and automation have made cryo-EM a popular method for producing near-atomic resolution structures of a variety of proteins and complexes. Sample preparation is still a limiting factor in collecting high quality data. Thickness of the vitreous ice in which the particles are embedded is one of the many variables that need to be optimized for collection of the highest quality data. Here we present two methods, using either an energy filter or scattering outside the objective aperture, to measure ice thickness for potentially every image collected. Unlike geometrical or tomographic methods, these can be implemented directly in the single particle collection workflow without interrupting or significantly slowing down data collection. We describe the methods as implemented into the Leginon/Appion data collection workflow, along with some examples from test cases. Routine monitoring of ice thickness should prove helpful for optimizing sample preparation, data collection, and data processing.

High-Resolution Ice-Sounding Radar Measurements of Ice Thickness Over East Antarctic Ice Sheet as a Part of Chinese National Antarctic Research Expedition

This paper presents the ice thickness, fine resolution internal reflecting horizons (IRHs), and distinct bottom topography measurements of Chinese Kunlun Station and Grove Mountains, Antarctica, derived from sounding these glaciers with a high-resolution radar. To enable the development of next-generation ice-sheet models, we need information on IRHs, bottom topography, and basal conditions. To this end, we performed measurements with the progressively improved ice-sounding radar system, currently known as the high-resolution ice-sounding radar developed by the Key Laboratory of Electromagnetic Radiation and Sensing Technology of Institute of Electronics, Chinese Academy of Sciences, Beijing, China. We processed the collected data using focused synthetic aperture radar (SAR) algorithm named the modified range migration algorithm using curvelets and the modified nonlinear chirp scaling algorithm to improve radar sensitivity and reduce along-track surface clutter. Representative results from selected transects indicate that we successfully sounded 3-km-thick ice with a fine resolution of 0.75 m. In this paper, we provide a brief description of the radar system, discuss the focused SAR processing algorithms, and provide sample results to demonstrate the successful sounding of the ice sheet in Antarctica.

Double-frequency ground penetrating radar for measurement of ice thickness and water depth in rivers and canals: Development, verification and application

In situ measurement of the ice thickness and water depth distributions in a river, canal, or lake is one of the most important means of measuring and responding to an ice flood. In view of the shortcomings of the use of single-frequency ground penetrating radar (GPR) for simultaneous measurement of ice thickness and water depth, a double-frequency GPR with 100 and 1500 MHz antennas was developed in the present study. To avoid interference between the two antennas with different frequencies, there is a 15-ms time-sharing operation. Similarity analysis of neighboring GPR signals was used to obtain the accuracy value of two-way travel time. The proposed GPR system also enables the real-time determination of the longitude and latitude coordinates of the measurement site using a real-time kinematic (RTK) system. The system measured the ice thickness and water depth at 97 and 99 sites, respectively, in the Mohe section of the Heilongjiang River, and the Togtoh, Inner Mongolia section of the Yellow River; both in China. Comparison of the measurements to those obtained by traditional drill-hole methods, yielded a mean value of absolute errors of 0.12 m and 0.04 m for water depth and ice thickness, respectively. The system was further used to measure the distribution of the ice thickness and water depth over a distance of 0.5 km in the Longdao Wharf section of the Heilongjiang River, and over a distance of 67 km in the Mohe section of the river,between the entrance of the Luogu River and Beiji village. The measurement results indicate that the developed double-frequency GPR is suitable for in situ measurement of ice thickness and water depth in rivers during ice periods.

Ice discharge error estimates using different cross-sectional area approaches: a case study for the Canadian High Arctic, 2016/17

ABSTRACTWe analyse the various error sources in the estimation of ice discharge through flux gates, distinguishing the cases with ice-thickness data available for glacier cross-sections or only along the centreline. For the latter, we analyse the performance of three U-shaped cross-sectional approaches. We apply this methodology to glaciers of the Canadian High Arctic. The velocity field is the main error source for small and medium-size glaciers (discharge <100 Mt a−1) with low velocities (<100 m a−1), while for large glaciers (discharge >100 Mt a−1) with high velocities (>100 m a−1) the error in cross-sectional area dominates. Thinning/thickening between ice-thickness and velocity measurements should be considered, as it implies systematic errors up to 8% in our study. The U-shaped parabolic approach, which allows for an adjusted estimation when the ice-thickness measurement point is displaced from the glacier centreline, performs best, with small bias and admissible standard error. We observe an increase of ice discharge from the main glaciers (Trinity and Wykeham) of the Prince of Wales Icefield from 2015 to 2016, by 5 and 20%, respectively, followed by a decrease in 2017, by 10 and 15%, respectively. Belcher Glacier, of the Devon Ice Cap, maintains similar discharges during 2015–17.

Ice thickness measurements in Lake Erie during the winter of 2010–2011

Abstract Time series measurements of ice thicknesses were made at either 1 or 2 Hz at 6 locations in the western part of Lake Erie's central basin during the winter of 2010–2011. Ice was observed over approximately 80 days beginning in late December and continuing through mid-March. Deformation and ridging of ice occurred frequently and produced ice thicknesses of up to 10 m, and over 6 m at all stations. The measurements show considerable variability (up to several meters) between stations, even when the distance between them is

Antarctic ice sheet thickness estimation using the horizontal-to-vertical spectral ratio method with single-station seismic ambient noise

Abstract. We report on a successful application of the horizontal-to-vertical spectral ratio (H / V) method, generally used to investigate the subsurface velocity structures of the shallow crust, to estimate the Antarctic ice sheet thickness for the first time. Using three-component, five-day long, seismic ambient noise records gathered from more than 60 temporary seismic stations located on the Antarctic ice sheet, the ice thickness measured at each station has comparable accuracy to the Bedmap2 database. Preliminary analysis revealed that 60 out of 65 seismic stations on the ice sheet obtained clear peak frequencies (f0) related to the ice sheet thickness in the H / V spectrum. Thus, assuming that the isotropic ice layer lies atop a high velocity half-space bedrock, the ice sheet thickness can be calculated by a simple approximation formula. About half of the calculated ice sheet thicknesses were consistent with the Bedmap2 ice thickness values. To further improve the reliability of ice thickness measurements, two-type models were built to fit the observed H / V spectrum through non-linear inversion. The two-type models represent the isotropic structures of single- and two-layer ice sheets, and the latter depicts the non-uniform, layered characteristics of the ice sheet widely distributed in Antarctica. The inversion results suggest that the ice thicknesses derived from the two-layer ice models were in good concurrence with the Bedmap2 ice thickness database, and that ice thickness differences between the two were within 300 m at almost all stations. Our results support previous finding that the Antarctic ice sheet is stratified. Extensive data processing indicates that the time length of seismic ambient noise records can be shortened to two hours for reliable ice sheet thickness estimation using the H / V method. This study extends the application fields of the H / V method and provides an effective and independent way to measure ice sheet thickness in Antarctica.

Assessing Sea Ice Trafficability in a Changing Arctic

Arctic sea ice has undergone rapid changes during the last few decades, with negative implications for over-ice travel and on-ice operations, which benefit from services provided by the sea ice. A Parameter-based Trafficability Hierarchy (PATH) is presented here as a framework for developing quantitative assessment strategies that can guide planning and execution of operations on or near sea ice and quantify the impacts of recent changes on ice use. A PATH assessment has been completed for three case studies in Arctic Alaska. These cases, which correspond to a range of different icescapes and ice uses, identify and quantify different parameters linked to trafficability and safe operations. For ice road applications, PATH was used to determine an ice thickness compensation factor, a factor increasing the minimal thickness threshold for operations, to help translate sporadic auger ice thickness measurements along the Kotzebue–Kiana community ice road into an envelope for safe operations. A compensation factor as high as 1.5 was found to be necessary to ensure safety because of the high local thickness variability that is currently a concern for ice road operators. A PATH assessment of ice roughness for ice trail routing at Utqiaġvik draws on satellite remote sensing and is relevant for over-ice travel in general, including escape, evacuation, and rescue. We compared the routing of local snowmobile trails with Synthetic Aperture Radar (SAR) data products to identify specific ranges of ice conditions, roughness, and topography favored for ice trail construction. The same combination of data sources was used to identify potentially beneficial trail routes. Finally, an ice stability and safety assessment was completed for ice road construction and maintenance by industry near the Northstar Island oil production facility. We evaluated small-scale ice displacement data obtained from SAR interferometry to infer internal ice strain and stress and used these data in assessing the potential for fractures to reduce load-bearing capacity.

Spatially distributed ice-thickness modelling for Chhota Shigri Glacier in western Himalayas, India

ABSTRACTEstimation of glacier ice-thickness distribution is important for many glacio-hydrological applications such as runoff projections, glacial lake outburst flood (GLOF) predictions, future evolution of glaciers. Varieties of modelling approaches are available for estimating ice-thickness distribution depending upon the data availability. In the present study, we estimated the ice-thickness distribution and total ice volume of Chhota Shigri Glacier using an optimally parameterized Glacier Bed Topography version 2 (GlabTop2) model, a shallow ice approximation (SIA)-based spatially distributed approach. Among the model input parameters, the shape factor (f), a non-measurable factor compensating for unaccounted effects such as valley shape needs to be calibrated using field measurements of ice-thickness. However, the lack of direct measurements of ice-thickness over many glaciers worldwide restricts model calibration and effective implementation. Therefore, to overcome this limitation, in this study, a novel approach using a relationship between shape factor, glacier cross-sectional width, and ice-thickness at the centre of a cross-section has been proposed and also tested to estimate optimal shape factor of the study glacier. Additionally, a detailed analysis of the effect of Digital Elevation Model (DEM) resolution and shape factor parameterization on the modelled ice-thickness estimates indicate that improving either the DEM resolution or calibrating the shape factor individually will not lead to improved ice-thickness estimates. In fact, both are necessary for better estimation of ice-thickness distribution. The high resolution DEM used in this study is TerraSAR-X add-on for Digital Elevation Measurement (TanDEM-X) DEM of 10 × 10 m grid size. Finally, a comparison of the results from a previous study where they used Shuttle Radar Topography Mission (SRTM) 90 m DEM indicates that the improved parameterisation of GlabTop2 model has led to a reduction in the error bounds of the estimated ice-thickness including the total glacier stored ice volume for the year 2013, which is estimated to be 1.74 ± 0.25 km3. Furthermore, based on the obtained results, it can be said that the GlabTop2 model combined with the proposed parameterization approach is having enormous potential to be applied over the wide range of data scarce Himalayan glaciers to quantify reliable ice-thickness estimates.

Comparison of Freeboard Retrieval and Ice Thickness Calculation From ALS, ASIRAS, and CryoSat‐2 in the Norwegian Arctic to Field Measurements Made During the N‐ICE2015 Expedition

AbstractWe present freeboard measurements from airborne laser scanner (ALS), the Airborne Synthetic Aperture and Interferometric Radar Altimeter System (ASIRAS), and CryoSat‐2 SIRAL radar altimeter; ice thickness measurements from both helicopter‐borne and ground‐based electromagnetic‐sounding; and point measurements of ice properties. This case study was carried out in April 2015 during the N‐ICE2015 expedition in the area of the Arctic Ocean north of Svalbard. The region is represented by deep snow up to 1.12 m and a widespread presence of negative freeboards. The main scattering surfaces from both CryoSat‐2 and ASIRAS are shown to be closer to the snow freeboard obtained by ALS than to the ice freeboard measured in situ. This case study documents the complexity of freeboard retrievals from radar altimetry. We show that even under cold (below −15°C) conditions the radar freeboard can be close to the snow freeboard on a regional scale of tens of kilometers. We derived a modal sea‐ice thickness for the study region from CryoSat‐2 of 3.9 m compared to measured total thickness 1.7 m, resulting in an overestimation of sea‐ice thickness on the order of a factor 2. Our results also highlight the importance of year‐to‐year regional scale information about the depth and density of the snowpack, as this influences the sea‐ice freeboard, the radar penetration, and is a key component of the hydrostatic balance equations used to convert radar freeboard to sea‐ice thickness.

Near-margin ice thickness and subglacial water routing, Leverett Glacier, Greenland

ABSTRACT Ice thickness measurements near the margin of the Greenland Ice Sheet (GrIS) are relatively sparse, presenting issues for modeling ice-flow dynamics, ice-sheet change, and subglacial hydrology. We acquired near-margin ice thickness data at Leverett Glacier, west Greenland, using a highly portable, low power, ground-penetrating radar operating at 10–80 MHz. Ice-thickness measurements, to a maximum of 270 m, were incorporated into the BedMachine model of ice thickness, created using mass conservation methods. The new data significantly modified the modeled ice thickness, and hence bed elevation and routing of subglacial water, in both the Leverett and adjacent Russell Glacier. Although the revised modeled basal topography and subglacial hydrology are consistent with observations, our new data unrealistically reduced the overall size of the Leverett Glacier hydrological catchment. Additional ice-thickness measurements are therefore required to realistically constrain subglacial topography and subglacial hydrological routing in this area. Our work improves understanding of the basal topography and the subglacial hydrology of Leverett Glacier, with implications for glacier dynamics and assessments of water piracy between catchments in the marginal zone of the GrIS, and for the interpolation of ice-thickness grids using mass conservation methods.

Linking Regional Winter Sea Ice Thickness and Surface Roughness to Spring Melt Pond Fraction on Landfast Arctic Sea Ice

The Arctic sea ice cover has decreased strongly in extent, thickness, volume and age in recent decades. The melt season presents a significant challenge for sea ice forecasting due to uncertainty associated with the role of surface melt ponds in ice decay at regional scales. This study quantifies the relationships of spring melt pond fraction (fp) with both winter sea ice roughness and thickness, for landfast first-year sea ice (FYI) and multiyear sea ice (MYI). In 2015, airborne measurements of winter sea ice thickness and roughness, as well as high-resolution optical data of melt pond covered sea ice, were collected along two ~5.2 km long profiles over FYI- and MYI-dominated regions in the Canadian Arctic. Statistics of winter sea ice thickness and roughness were compared to spring fp using three data aggregation approaches, termed object and hybrid-object (based on image segments), and regularly spaced grid-cells. The hybrid-based aggregation approach showed strongest associations because it considers the morphology of the ice as well as footprints of the sensors used to measure winter sea ice thickness and roughness. Using the hybrid-based data aggregation approach it was found that winter sea ice thickness and roughness are related to spring fp. A stronger negative correlation was observed between FYI thickness and fp (Spearman rs = −0.85) compared to FYI roughness and fp (rs = −0.52). The association between MYI thickness and fp was also negative (rs = −0.56), whereas there was no association between MYI roughness and fp. 47% of spring fp variation for FYI and MYI can be explained by mean thickness. Thin sea ice is characterized by low surface roughness allowing for widespread ponding in the spring (high fp) whereas thick sea ice has undergone dynamic thickening and roughening with topographic features constraining melt water into deeper channels (low fp). This work provides an important contribution towards the parameterizations of fp in seasonal and long-term prediction models by quantifying linkages between winter sea ice thickness and roughness, and spring fp.

Record low sea-ice concentration in the central Arctic during summer 2010

The Arctic sea-ice extent has shown a declining trend over the past 30 years. Ice coverage reached historic minima in 2007 and again in 2012. This trend has recently been assessed to be unique over at least the last 1450 years. In the summer of 2010, a very low sea-ice concentration (SIC) appeared at high Arctic latitudes—even lower than that of surrounding pack ice at lower latitudes. This striking low ice concentration—referred to here as a record low ice concentration in the central Arctic (CARLIC)—is unique in our analysis period of 2003–15, and has not been previously reported in the literature. The CARLIC was not the result of ice melt, because sea ice was still quite thick based on in-situ ice thickness measurements. Instead, divergent ice drift appears to have been responsible for the CARLIC. A high correlation between SIC and wind stress curl suggests that the sea ice drift during the summer of 2010 responded strongly to the regional wind forcing. The drift trajectories of ice buoys exhibited a transpolar drift in the Atlantic sector and an eastward drift in the Pacific sector, which appeared to benefit the CARLIC in 2010. Under these conditions, more solar energy can penetrate into the open water, increasing melt through increased heat flux to the ocean. We speculate that this divergence of sea ice could occur more often in the coming decades, and impact on hemispheric SIC and feed back to the climate.

Geometry and ice dynamics of the Darwin–Hatherton glacial system, Transantarctic Mountains

ABSTRACTThe Darwin–Hatherton Glacial system (DHGS) connects the East Antarctic Ice Sheet (EAIS) with the Ross Ice Shelf and is a key area for understanding past variations in ice thickness of surrounding ice masses. Here we present the first detailed measurements of ice thickness and grounding zone characteristics of the DHGS as well as new measurements of ice velocity. The results illustrate the changes that occur in glacier geometry and ice flux as ice flows from the polar plateau and into the Ross Ice Shelf. The ice discharge and the mean basal ice shelf melt for the first 8.5 km downstream of the grounding line amount to 0.24 ± 0.05 km3 a−1 and 0.3 ± 0.1 m a−1, respectively. As the ice begins to float, ice thickness decreases rapidly and basal terraces develop. Constructed maps of glacier geometry suggest that ice drainage from the EAIS into the Darwin Glacier occurs primarily through a deep subglacial canyon. By contrast, ice thins to <200 m at the head of the much slower flowing Hatherton Glacier. The glaciological field study establishes an improved basis for the interpretation of glacial drift sheets at the link between the EAIS and the Ross Ice Sheet.

Near‐shore talik development beneath shallow water in expanding thermokarst lakes, Old Crow Flats, Yukon

AbstractIt is generally assumed that permafrost is preserved beneath shallow lakes and ponds in the Western North American Arctic where water depth is less than about two thirds of the late‐winter lake ice thickness. Here we present field observations of talik development beneath water as shallow as 0.2 m despite a lake ice thickness of 1.5 m, in Old Crow Flats (OCF), YT. Conditions leading to the initiation and development of taliks beneath shallow water were investigated with field measurements of shore erosion rates, bathymetry, ice thickness, snow accumulation, and lake bottom temperature near the shores of two expanding lakes in OCF. The sensitivity of talik development to variations in lake bottom thermal regime was then investigated numerically. Where ice reached the lake bottom, talik development was controlled by the ratio of freezing degree days to thawing degree days at the lake bottom (FDDlb/TDDlb). In some cases, spatial variations in on‐ice snow depth had a minimal effect on annual mean lake bottom temperature (Tlb) but caused sufficient variations in FDDlb/TDDlb to influence talik development. Where Tlb was close to but greater than 0°C simulations indicated that the thermal offset allowed permafrost aggradation to occur under certain conditions, resulting in irregular near‐shore talik geometries. The results highlight the sensitivity of permafrost to small changes in lake bottom thermal conditions where the water column freezes through in early winter and indicate the occurrence of permafrost degradation beneath very shallow water in the near‐shore zone of Arctic ponds and lakes.