Yedoma Permafrost Genesis: Over 150 Years of Mystery and Controversy

Since the discovery of frozen megafauna carcasses in Northern Siberia and Alaska in the early 1800s, the Yedoma phenomenon has attracted many Arctic explorers and scientists. Exposed along coastal and riverbank bluffs, Yedoma often appears as large masses of ice with some inclusions of sediment. The ground ice particularly mystified geologists and geographers, and they considered sediment within Yedoma exposures to be a secondary and unimportant component. Numerous scientists around the world tried to explain the origin of Yedoma for decades, even though some of them had never seen Yedoma in the field. The origin of massive ice in Yedoma has been attributed to buried surface ice (glaciers, snow, lake ice, and icings), intrusive ice (open system pingo), and finally to ice wedges. Proponents of the last hypothesis found it difficult to explain a vertical extent of ice wedges, which in some cases exceeds 40 m. It took over 150 years of intense debates to understand the process of ice-wedge formation occurring simultaneously (syngenetically) with soil deposition and permafrost aggregation. This understanding was based on observations of the contemporary formation of syngenetic permafrost with ice wedges on the floodplains of Arctic rivers. It initially was concluded that Yedoma was a floodplain deposit, and it took several decades of debates to understand that Yedoma is of polygenetic origin. In this paper, we discuss the history of Yedoma studies from the early 19th century until the 1980s—the period when the main hypotheses of Yedoma origin were debated and developed.

Microbial and Geochemical Evidence of Permafrost Formation at Mamontova Gora and Syrdakh, Central Yakutia

Biotracers marking the geologic history and permafrost evolution in Central Yakutia, including Yedoma Ice Complex (IC) deposits, were identified in a multiproxy analysis of water chemistry, isotopic signatures, and microbial datasets. The key study sections were the Mamontova Gora and Syrdakh exposures, well covered in the literature. In the Mamontova Gora section, two distinct IC strata with massive ice wedges were described and sampled, the upper and lower IC strata, while previously published studies focused only on the lower IC horizon. Our results suggest that these two IC horizons differ in water origin of wedge ice and in their cryogenic evolution, evidenced by the differences in their chemistry, water isotopic signatures and the microbial community compositions. Microbial community similarity between ground ice and host deposits is shown to be a proxy for syngenetic deposition and freezing. High community similarity indicates syngenetic formation of ice wedges and host deposits of the lower IC horizon at the Mamontova Gora exposure. The upper IC horizon in this exposure has much lower similarity metrics between ice wedge and host sediments, and we suggest epigenetic ice wedge development in this stratum. We found a certain correspondence between the water origin and the degree of evaporative transformation in ice wedges and the microbial community composition, notably, the presence of Chloroflexia bacteria, represented by Gitt-GS-136 and KD4-96 classes. These bacteria are absent at the ice wedges of lower IC stratum at Mamontova Gora originating from snowmelt, but are abundant in the Syrdakh ice wedges, where the meltwater underwent evaporative isotopical fractionation. Minor evaporative transformation of water in the upper IC horizon of Mamontova Gora, whose ice wedges formed by meltwater that was additionally fractionated corresponds with moderate abundance of these classes in its bacterial community.

Paleontological monument of nature “Mamontova Gora” exposure (Yakutia, Eastern Siberia, Russia)

The Mamontova Gora (“Mammoth Mountain”) exposure is the Neogene – Pleistocene key section of Siberia. This outcrop is located in the lower reaches of the Aldan River, 325 km above its mouth and extends for almost 12 km. It consists of an 80-meter structural plateau (80-meter terrace), 50- and 30-meter alluvial terraces. Sediments from the Middle Miocene (16-10 Ma) to the Upper Pleistocene are exposed on the 80-meter terrace. The basement of the 50-meter terrace is composed of Middle Miocene sediments, overlain by Pleistocene sediments. On a younger 30 m terrace, the deposits are dated from the Upper Pliocene to the Upper Pleistocene. The Mamontova Gora outcrop is one of the richest localities of the Neogene flora of Eurasia. There are numerous finds of remains of Miocene evergreen and thermophilic plants (tree stumps, leaf imprints, cones, nuts, seeds). More than 250 genera of fossil plants have been found on Mamontova Gora. This outcrop is also well known to paleontologists due to the abundance of bone remains of mammals of the Middle Pleistocene (early type mammoth, eastern horse, broad-fronted moose, long-horned bison) and Late Pleistocene (representatives of the mammoth fauna: woolly mammoth, woolly rhinoceros, Lena horse, reindeer, saiga-antelope, steppe bison, Arctic fox, wolverine, cave lion, etc.). It was revealed that the ancient frozen sediments on the Mamontova Gora outcrop abound with viable microorganisms and traces of their vital activity. A strain of microbe Bacillus sp. was isolated from ~ 2 - 3 Ma permafrost layers of this outcrop. A large group of microorganisms including fungi was isolated from the ancient ice wedge. Pleistocene permafrost deposits contain invertase, urease, catalase and dehydrogenase enzimes. Mamontova Gora is a unique geological object in Russia. By the decree of the Council of Ministers of the Yakutia Autonomous Soviet Socialist Republic of 18.02.1987 No. 56 Mamontova Gora was given the status of a “natural monument” and a specially protected natural area of regional significance. The article presents the main results of studies of ancient flora and faunas of Mamontova Gora.

Shared-Energy Prediction Model for Ship-Ice Interactions

Low- and non-ice-class ship-ice interactions are modelled with a shared-energy approach, which typically models the internal mechanics with nonlinear finite element methods. For applications like the preliminary design phase and quick operational assessments of the ship’s structural capabilities, a finite element shared-energy approach can be time consuming and information intensive, therefore, an analytical share-energy algorithm is proposed. The proposed algorithm applies the upper bound energy methodology by equating the external collision energy, determined with the Popov collision model (Popov, et al., 1967), to the sum of the internal ice and structural response energies. The distribution of the internal energy, between the ice and the structure, is determined by iterating through possible shared contact forces until the sum of the internal response energies equals the external energy introduced into the system. The ice-crushing energy is modelled with Daley’s (1999) energy based ice collision force models, and the internal structural strain energy is modelled through a combination of classical beam theory and design of experiments methodology. The proposed model is benchmarked against a finite element ice wedge-ship grillage structure interaction.

SMALL-AMPLITUDE DISCONTINUOUS DISTURBANCES IN PLEISTOCENE DEPOSITS OF THE VOLYN-PODILSKA UPLAND

The paper presents the results of the study of the small-amplitude discon¬tinuous disturbances of the possibly cryogenic (thermokarst) origin. The dislocations were found in the outcrops of Middle and Upper Pleistocene sediments of the Volyn-Podilska Upland, accumulated in periglacial or sub-periglacial conditions. The distur¬bances are represented mostly by the micro-normal faults and also by sheared fractures and are very similar to tectonic (seismogenic) discontinuities. The tectonotypic fractures in the near-surface deposits of the Pleistocene terraces of Western Bug and Styr (five sections within Volhynian Upland, four of them – in the valley of Bug), as well as in the cover of the Late Pleistocene sediments on the slope of the valley of Dniester (Galician Prydnisterya) are subjected to consideration, analysis and interpretation. In the last location the ruptures are represented mostly by the dis¬turbances identified as sheared fractures. In all others there are small-amplitude normal faults. One reverse fault, timed to an ice-wedge cast, was also revealed. Typical micro-normal faults of all sections are steep and have a number of other common features, which testifies to the same or almost identical mechanism of their formation. These features, in particular, are as follows: 1) insignificant (usually up to 2–2.5 m) length in cross-section and small (several centimeters) amplitude of displacement along the rupture plane; 2) gradual attenuation of the fractures up and down the section. All micro-normal faults are confined to sediments (thicknesses) that are partially or completely composed of sand. The formation of the micro-normal faults and other examined ruptures can be ex¬plained by the uneven compaction and the gravitational subsidence of the rocks, and in the section on the slope of the Dniester valley – also by their displacement down along the slope. It is probable that these processes occurred due to: 1) the degradation of the permafrost; 2) the dehydration of the sand deposits during a significant decrease in the groundwater levels; 3) the melting of the buried layers and lenses of snow, which were accumulated during the winter season in the thickness of sandy the niveo-aeolian deposits. In the outcrops of this terrace, they occur no less frequently than the confidently identified ice wedge pseudomorphs. Key words: small-amplitude disturbances; microfaults; thermokarst; Volyn-Podilska Upland.

New insights into the drainage of inundated ice-wedge polygons using fundamental hydrologic principles

The pathways and timing of drainage from the inundated centers of ice-wedge polygons in a warming climate have important implications for carbon flushing, advective heat transport, and transitions from methane to carbon dioxide dominated emissions. Here, we expand on previous research using a recently developed analytical model of drainage from a low-centered polygon. Specifically, we perform (1) a calibration to field data identifying necessary model refinements and (2) a rigorous model sensitivity analysis that expands on previously published indications of polygon drainage characteristics. This research provides intuition on inundated polygon drainage by presenting the first in-depth analysis of drainage within a polygon based on hydrogeological first principles. We verify a recently developed analytical solution of polygon drainage through a calibration to a season of field measurements. Due to the parsimony of the model, providing the potential that it could fail, we identify the minimum necessary refinements that allow the model to match water levels measured in a low-centered polygon. We find that (1) the measured precipitation must be increased by a factor of around 2.2, and (2) the vertical soil hydraulic conductivity must decrease with increasing thaw depth. Model refinement (1) accounts for runoff from rims into the ice-wedge polygon pond during precipitation events and possible rain gauge undercatch, while refinement (2) accounts for the decreasing permeability of deeper soil layers. The calibration to field measurements supports the validity of the model, indicating that it is able to represent ice-wedge polygon drainage dynamics. We then use the analytical solution in non-dimensional form to provide a baseline for the effects of polygon aspect ratios (radius to thaw depth) and coefficient of hydraulic conductivity anisotropy (horizontal to vertical hydraulic conductivity) on drainage pathways and temporal depletion of ponded water from inundated ice-wedge polygon centers. By varying the polygon aspect ratio, we evaluate the relative effect of polygon size (width), inter-annual increases in active-layer thickness, and seasonal increases in thaw depth on drainage. The results of our sensitivity analysis rigorously confirm a previous analysis indicating that most drainage through the active layer occurs along an annular region of the polygon center near the rims. This has important implications for transport of nutrients (such as dissolved organic carbon) and advection of heat towards ice-wedge tops. We also provide a comprehensive investigation of the effect of polygon aspect ratio and anisotropy on drainage timing and patterns, expanding on previously published research. Our results indicate that polygons with large aspect ratios and high anisotropy will have the most distributed drainage, while polygons with large aspect ratios and low anisotropy will have their drainage most focused near their periphery and will drain most slowly. Polygons with small aspect ratios and high anisotropy will drain most quickly. These results, based on parametric investigation of idealized scenarios, provide a baseline for further research considering the geometric and hydraulic complexities of ice-wedge polygons.

COUNTING ICE-WEDGE POLYGONS FROM SPACE: USE OF COMMERCIAL SATELLITE IMAGERY TO MONITOR CHANGING ARCTIC POLYGONAL TUNDRA

The microtopography associated with ice wedge polygons (IWPs) governs the Arctic ecosystem from local to regional scales due to the impacts on the flow and storage of water and therefore, vegetation and carbon. Increasing subsurface temperatures in Arctic permafrost landscapes cause differential ground settlements followed by a series of adverse microtopographic transitions at sub decadal scale. The entire Arctic has been imaged at 0.5 m or finer resolution by commercial satellite sensors. Dramatic microtopographic transformation of low-centered into high-centered IWPs can be identified using sub-meter resolution commercial satellite imagery. In this exploratory study, we have employed a Deep Learning (DL)-based object detection and semantic segmentation method named the Mask R-CNN to automatically map IWPs from commercial satellite imagery. Different tundra vegetation types have distinct spectral, spatial, textural characteristics, which in turn decide the semantics of overlying IWPs. Landscape complexity translates to the image complexity, affecting DL model performances. Scarcity of labelled training images, inadequate training samples for some types of tundra and class imbalance stand as other key challenges in this study. We implemented image augmentation methods to introduce variety in the training data and trained models separately for tundra types. Augmentation methods show promising results but the models with separate tundra types seem to suffer from the lack of annotated data.

A Quantitative Graph-Based Approach to Monitoring Ice-Wedge Trough Dynamics in Polygonal Permafrost Landscapes

In response to increasing Arctic temperatures, ice-rich permafrost landscapes are undergoing rapid changes. In permafrost lowlands, polygonal ice wedges are especially prone to degradation. Melting of ice wedges results in deepening troughs and the transition from low-centered to high-centered ice-wedge polygons. This process has important implications for surface hydrology, as the connectivity of such troughs determines the rate of drainage for these lowland landscapes. In this study, we present a comprehensive, modular, and highly automated workflow to extract, to represent, and to analyze remotely sensed ice-wedge polygonal trough networks as a graph (i.e., network structure). With computer vision methods, we efficiently extract the trough locations as well as their geomorphometric information on trough depth and width from high-resolution digital elevation models and link these data within the graph. Further, we present and discuss the benefits of graph analysis algorithms for characterizing the erosional development of such thaw-affected landscapes. Based on our graph analysis, we show how thaw subsidence has progressed between 2009 and 2019 following burning at the Anaktuvuk River fire scar in northern Alaska, USA. We observed a considerable increase in the number of discernible troughs within the study area, while simultaneously the number of disconnected networks decreased from 54 small networks in 2009 to only six considerably larger disconnected networks in 2019. On average, the width of the troughs has increased by 13.86%, while the average depth has slightly decreased by 10.31%. Overall, our new automated approach allows for monitoring ice-wedge dynamics in unprecedented spatial detail, while simultaneously reducing the data to quantifiable geometric measures and spatial relationships.