Refine
Has Fulltext
- no (4)
Document Type
- Article (4) (remove)
Language
- English (4) (remove)
Is part of the Bibliography
- yes (4)
Keywords
- salt diffusion (2)
- subsea permafrost (2)
- Bykovsky Peninsula (1)
- CryoGRID (1)
- Lena Delta (1)
- Siberia (1)
- electrical resistivity (1)
- electromagnetic (1)
- geophysics (1)
- lagoon (1)
- lake (1)
- permafrost (1)
- talik (1)
- thermokarst (1)
- thermokarst lake (1)
Institute
Thermokarst lakes are prevalent in Arctic coastal lowland regions and sublake permafrost degradation and talik development contributes to greenhouse gas emissions by tapping the large permafrost carbon pool. Whereas lateral thermokarst lake expansion is readily apparent through remote sensing and shoreline measurements, sublake thawed sediment conditions and talik growth are difficult to measure. Here we combine transient electromagnetic surveys with thermal modeling, backed up by measured permafrost properties and radiocarbon ages, to reveal closed-talik geometry associated with a thermokarst lake in continuous permafrost. To improve access to talik geometry data, we conducted surveys along three transient electromagnetic transects perpendicular to lakeshores with different decadal-scale expansion rates of 0.16, 0.38, and 0.58m/year. We modeled thermal development of the talik using boundary conditions based on field data from the lake, surrounding permafrost and a borehole, independent of the transient electromagnetics. A talik depth of 91m was determined from analysis of the transient electromagnetic surveys. Using a lake initiation age of 1400years before present and available subsurface properties the results from thermal modeling of the lake center arrived at a best estimate talk depth of 80m, which is on the same order of magnitude as the results from the transient electromagnetic survey. Our approach has provided a noninvasive estimate of talik geometry suitable for comparable settings throughout circum-Arctic coastal lowland regions.
As the Arctic coast erodes, it drains thermokarst lakes, transforming them into lagoons, and, eventually, integrates them into subsea permafrost. Lagoons represent the first stage of a thermokarst lake transition to a marine setting and possibly more saline and colder upper boundary conditions. In this research, borehole data, electrical resistivity surveying, and modeling of heat and salt diffusion were carried out at Polar Fox Lagoon on the Bykovsky Peninsula, Siberia. Polar Fox Lagoon is a seasonally isolated water body connected to Tiksi Bay through a channel, leading to hypersaline waters under the ice cover. The boreholes in the center of the lagoon revealed floating ice and a saline cryotic bed underlain by a saline cryotic talik, a thin ice-bearing permafrost layer, and unfrozen ground. The bathymetry showed that most of the lagoon had bedfast ice in spring. In bedfast ice areas, the electrical resistivity profiles suggested that an unfrozen saline layer was underlain by a thick layer of refrozen talik. The modeling showed that thermokarst lake taliks can refreeze when submerged in saltwater with mean annual bottom water temperatures below or slightly above 0 degrees C. This occurs, because the top-down chemical degradation of newly formed ice-bearing permafrost is slower than the refreezing of the talik. Hence, lagoons may precondition taliks with a layer of ice-bearing permafrost before encroachment by the sea, and this frozen layer may act as a cap on gas migration out of the underlying talik.
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.25m/year was 16.7 m below the seabed 350m 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-20m below the seabed at 350m 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.
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.