Sub-Surface Carbon Stocks in Northern Taiga Landscapes Exposed in the Batagay Megaslump, Yana Upland, Yakutia

: The most massive and fast-eroding thaw slump of the Northern Hemisphere located in the Yana Uplands of Northern Yakutia was investigated to assess in detail the cryogenic inventory and carbon pools of two distinctive Ice Complex stratigraphic units and the uppermost cover deposits. Di ﬀ erentiating into modern and Holocene near-surface layers (active layer and shielding layer), highest total carbon contents were found in the active layer (18.72 kg m − 2 ), while the shielding layer yielded a much lower carbon content of 1.81 kg m − 2 . The late Pleistocene upper Ice Complex contained 10.34 kg m − 2 total carbon, and the mid-Pleistocene lower Ice Complex 17.66 kg m − 2 . The proportion of organic carbon from total carbon content is well above 70% in all studied units with 94% in the active layer, 73% in the shielding layer, 83% in the upper Ice Complex and 79% in the lower Ice Complex. Inorganic carbon is low in the overall structure of the deposits.


Introduction
Pioneering fieldwork in the Batagay megaslump to study its cryogenic inventory and origin was undertaken in 2011 [1]. The prominent ice-rich Ice Complex units and ice-poor sand units as well as enclosed ground ice exposed in the Batagay megaslump accumulated during the Pleistocene under the continental conditions of the non-glaciated region of eastern Siberia belonging to the western part of Beringia [1]. Currently, the Batagay megaslump and its exposed permafrost are increasingly studied by various disciplines to assess its modern morphodynamics such as erosion rates and budgets [2,3] or the behavior of living microorganisms inside the frozen zone [4]. Most studies, however, deal with the Quaternary frozen deposits, the fossil inventory and ground ice.
Based on current cryolithological research and various dating approaches the permafrost exposed in the Batagay megaslump differentiates into the following seven units [1,5,6] starting from the slump Land 2020, 9,305 3 of 16 30 to 50 kg m −2 of organic carbon are stored in the uppermost 3 m of permafrost on a circum-arctic scale [19].
The carbon component in the landscapes of the circumpolar region is one of the criteria for observing the influence of anthropogenic and natural factors on permafrost. Changes occurred in the violation of pristine properties of the soil cover entail the inevitable loss of carbon from the active layer and permafrost as CO 2 , CH 4 and physical removal by erosion into riverine and marine realms. The permafrost exposed in the Batagay megaslump is an important example of how natural permafrost feature experiences overall degradation, changing as a result of human activities, the geological structure of the territory and due to its geographical location. However, according to Masyagina and Menyailo [20], the permafrost of Siberia will remain stable soon, but it is still a potential source of a vast amount of greenhouse gases, and ongoing release into the atmosphere will result in a further temperature increase, which is the permafrost carbon feedback.
At the present time, the permafrost vulnerability is recorded by researchers in various geographic regions of the Earth. Its most severe degradation occurs in the Arctic, caused by climate change [21,22]. In Alaska, increased precipitation over the past 5 years threatens once stable frozen soil, while more precipitation leads to deeper thawing of soil [23].
Batagay megaslump is the one and only large thermokarst formation in the Northern Hemisphere of Eurasia. Its frozen rocks have not been sufficiently studied in terms of biogenic elements; therefore, we have made attempts to study one of the main elements of the ecosystem such as carbon and its fund in the distinctive elements of the megaslump. It is important to note that the aim of the article is not to compare the results with other objects of the cryolithozone, since such objects do not currently exist. Moreover, we have to show the current state of the Batagay megaslump, as well as partially and indirectly alleged reasons for its formation. Undoubtedly, this is a small part of the preliminary studies which will be continued in the future. We provide a literary analysis of previously conducted studies on Batagay megaslump from different backgrounds and the results of international teams of scientists.

Brief Natural Conditions of the Study Area
The climate of the Yana Upland region is severe and sharply continental. It is affected by the geographical location between the arctic and subarctic climatic zones as well as by isolation between mountain ranges. A detailed climate description of the area is presented in publication [11], so in our article, we will not focus on this.
The vegetation of the northern taiga is characterized by low species diversity, in mountain climates with predominantly insufficient atmospheric moisture. In the study area, forests are dominated by larch (Larix cajanderi) with undergrowth of birch, cedar dwarf pine, lingonberry and moss-lichen communities [24]. Landscapes are formed on permafrost in continuous distribution with mean annual ground temperatures from −7 to −4 • C. In the Verkhoyan-Kolyma region, the maximum permafrost thickness reaches 350 m in depth [25]. The presence of ice-rich permafrost in combination with widespread cryogenic processes such as thermokarst, thermo-erosion, solifluction and frost cracking, determine the appearance of periglacial landscapes and their development.
The research area is located in the tectonic Verkhoyansk anticlinal zone. The regional structure of the Verkhoyan-Kolyma folded region, including the Verkhoyansk meganticlinorium, bordered from the west by the Priverkhoyansk trough, and from the east by the structures of the Yano-Indigir synclinal zone. The territory is strongly dissected by a mountainous relief developing on the Mesozoic folded base [26].
Currently, the Verkhoyansk-Kolyma Mesozoic folding is distinguished by its outer part, which consists of the Verkhoyansk fold-thrust belt adjacent to the Siberian platform. Internally, it has a complex structure and includes a large number of intense folded linear zones and individual exotic blocks (terranes). The border of the inner and outer zones is the Adycha-Taryn thrust, which is located 50 and 130 km away from Batagay on a section of the Srednyaya river and the Verkhnyaya Adycha river. It is also a boundary general fault of the Chersky seismic zone [27]. The loose deposits in the area are underlain by a bedrock of the Triassic age of the Ladinsky tier such as sandstones, siltstones and conglomerates. Quaternary formations are represented by undivided frozen eluvium and diluvium consisting of sand, silt and clay including ground ice and organic matter.

Terms Used in the Article
There are some terms from Geocryology. The object of study of Geocryology is frozen rocks, including underground and ground accumulations of ice and snow [28]. The terms are widely used by Russian researchers and a number of foreign scientific teams conducting scientific experiments on the territory of the Russian cryolithozone. For understanding the terms in our work, their main definitions are given with the reference: Active layer is the uppermost layer of soil alternately freezing and thawing in the presence of permafrost at its base [29].
Ataxitic cryostructure is inherent only for finely dispersed biogenic soils, in which a complex combination of ice inclusions and mineral aggregates suspended in ice is observed [29].
Cryogenic processes are set of processes of physical, chemical and mineralogical changes and transformation of soils and rocks of the weathering crust, as well as hydrosphere at a negative temperature [29].
Cryolithozone or permafrost zone is part of the cryosphere, representing the upper part of the Earth's crust, characterized by negative temperatures of soils and rocks, also by the presence of ice or supercooled water in them [29].
Ice wedge is the ice, filling frost fissure, as well as other cracks in the field of permafrost. During thawing out of ice wedges, ditches up to 5-10 m at the depth and up to 10-15 m in width are formed, which are filled with soil mass [29].
Loess rocks (sediments) or loess are the loose, silty loamy rock of a fawn or gray-yellow color. The particle size distribution is dominated by the coarse dust fraction of 0.05-0.01 mm. [31].
Permafrost is a phenomenon of long-term freezing of the upper part of lithosphere, as well as the thickness and range of permafrost [29].
Shielding layer is a surface layer, which thaws during favorable climatic conditions, joins the active layer, but during unfavorable conditions it does thaw and is the upper layer of permafrost [32].
Ice schlieren is the ice inclusions in frozen ground, represented by veins, interlayers, and lenses of different orientations [29].
Thermokarst is the process of thawing of frozen grounds and ground ice, accompanied by the formation of subsidence forms of relief [29].
The Ice Complex or yedoma is loess-like syngenetic frozen soils of the Pleistocene age with high thickness, forming the northern coastal lowlands, the plains of Central Yakutia and the North-East. The Ice Complex is characterized by high ice content of sediments up to 70%-90% and large syngenetic ice lodes, with heights of up to 60-80 m [29].
Frosty weathering is a process of physical destruction of rocks that occurs under conditions of frequent temperature fluctuations above and below 0 • C, i.e., phase transitions of water into ice and back. Water freezing in cracks, increasing in volume, acts as a wedge separating pieces of rock from its main mass [29].
Land 2020, 9, 305 5 of 16 In this work, the use of the terms such as permafrost, soil and sediments should be identified with the synonym for the word "soil", since this group of terms is a multi-component system undergoing a complex weathering process. During fieldwork in March-April 2019, the thaw slump walls were sampled at four locations ( Figure 1). The sampled stratigraphic units exposed in the wall of the slump were marked as follows: A1 (lower Ice Complex), B1 (upper Ice Complex) close to the headwall of the slump, D1 (ravine) in the south-facing exposure and C1 (upper Ice Complex) in the north-facing exposure. Firstly, the walls of the sections were carefully cleaned to remove the outermost material in order to avoid contamination. Secondly, the sampling profiles were cryolithologically described in detail. Finally, we took samples of the stratigraphical horizons for carbon analysis and soil density calculations. Soil samples were taken between ice wedges, where there were no ice inclusions. This was done in order to eliminate ice from the carbon samples. It should be noted that the outcrop of icy rocks is characteristic of southern exposure, marker A1. For the rest of megaslump, we did not observe the same situation. The Batagay natural complex provides a detailed study of the system of typical key areas, divided into the above markers; in geography, these are called facies. Key areas boil down to the fact that a number of components of the local geographic system and their interconnection with each other, in our case, involves individual elementary layers with carbon reserves. This includes the features of the complex and the current processes taking place in it. Marking the position in space and in time allows for checking the observation material at another selected time, which ensures the receiving of mass and discrete data. So, the main differences between the markers are: A1 is the lower Ice Complex (deep permafrost layer) of the southern exposure of the megaslump at the maximum possible depth from the surface of 59.6 m, composed of sandy loam and silt particles. B1 is the upper Ice Complex, an annually increasing ravine passes through it and includes all elementary layers such as active, protective, and permafrost. It is a homogeneous genetic profile. C1 is the upper Ice Complex of the northern exposure of the megaslump, which includes only the active layer. D1 is a ravine (active layer), formed and expanding by temporary streams in the spring-summer period, and water flows down into the Batagay River.
The soil density (g cm −3 ) in natural composition was determined by taking samples every 10 cm into steel cylinders with a defined volume of 100 cm −3 . Knowing the mass of the cylinder with soil and the mass of the empty cylinder, we determined the difference in soil mass at certain moisture content. By determining the moisture percentage, we calculated the dry mass of the material. The density of the soil of undisturbed composition was obtained by dividing the mass of dry soil by cylinder volume.
The soil samples were processed in the laboratory at room temperature. The sediment was laid out on paper and dried to an air-dry state. After that, the samples were oven-dried for 10 h at a temperature of 105 • C to exclude excess moisture from the samples.
We determined the organic carbon (OC), inorganic carbon (IC) and total carbon (TC) content using a certified LECO RC612 multiphase carbon analyzer, according to DIN19539. The applied temperature programming allows dividing the various carbon forms into organic and inorganic forms. The maximum temperature for burning samples was 1100 • C in an oxidizing atmosphere.
Carbon stocks (kg m −2 ) were calculated for each selected layer. We calculated the stocks for the active and shielding layers and underlying permafrost and summed them to get the total carbon stock for the entire thickness. The calculation of carbon stocks was carried out, according to the Equation [33]: where S-carbon stock (kg m −2 ), H-soil layer depth, p-soil density (g cm −3 ) and X-average carbon content (%).
Statistical processing of the carbon data was carried out in the program StatSoft STATISTICA for Windows 13.3. To exclude distorted indicators in the sample set, a typical sample from the general set of observations was used. The data are presented as arithmetic means with the standard error of the mean values.
In this work, the use of the terms such as permafrost, soil and sediments should be identified with the synonym for the word "soil", since this group of terms is a multi-component system undergoing a complex weathering process. During fieldwork in March-April 2019, the thaw slump walls were sampled at four locations ( Figure 1). The sampled stratigraphic units exposed in the wall

General Cryostratigraphy of Permafrost Exposed in the Batagay Megaslump
The Batagay megaslump is the largest one in the Northern Hemisphere occupying an area of 81 hectares, and with a length of more than 1000 m and a width of 800 m in 2019 (see Figure 1). The slump expands at high rates along the perimeter of up to 0.026 km 2 year −1 (between 1991 and 2018) [3]. The exposed depth along the slump perimeter varies depending on the ice content of the permafrost deposits and the activity of thermo-erosion. It reaches up to about a 60 m depth in the central part of the headwall. The slope inclination is pre-defined by exogenous geological processes that formed the slopes of the adjacent Kirgilyakh and Khatyngnakh mounts, subsequent long-term aggradation of partly very ice-rich permafrost, i.e., the Ice Complex and thermo-erosion along with valley structures on the slopes by ground ice melt. The Batagay megaslump morphodynamics is characterized by active thermo-erosion at the walls and the fluvial erosion by temporary streams on the slump bottom. The latter transports eroded fine material with meltwater along the central slump channel into the floodplain of the Batagay River. However, substantial amounts of the eroded material, as well as permafrost remnants, remain on the slump bottom, forming a high relief of ridges and hills with different heights.
The The Batagay megaslump is the largest one in the Northern Hemisphere occupying an area of 81 hectares, and with a length of more than 1000 m and a width of 800 m in 2019 (see Figure 1). The slump expands at high rates along the perimeter of up to 0.026 km 2 year −1 (between 1991 and 2018) [3]. The exposed depth along the slump perimeter varies depending on the ice content of the permafrost deposits and the activity of thermo-erosion. It reaches up to about a 60 m depth in the central part of the headwall. The slope inclination is pre-defined by exogenous geological processes that formed the slopes of the adjacent Kirgilyakh and Khatyngnakh mounts, subsequent long-term aggradation of partly very ice-rich permafrost, i.e., the Ice Complex and thermo-erosion along with valley structures on the slopes by ground ice melt. The Batagay megaslump morphodynamics is characterized by active thermo-erosion at the walls and the fluvial erosion by temporary streams on the slump bottom. The latter transports eroded fine material with meltwater along the central slump channel into the floodplain of the Batagay River. However, substantial amounts of the eroded material, as well as permafrost remnants, remain on the slump bottom, forming a high relief of ridges and hills with different heights.
The   The wedge ice of the upper Ice Complex was pillar-shaped, resembling a solid wall of ice from a distance. Table 1 summarizes descriptions of the cover deposits (samples in profiles B1, C1 and D1) of the upper Ice Complex, mainly composed of sandy and dusty particles with a localized separate newly formed structure (ironing) and individual mechanical inclusions of pebbles and gravel, which formed one soil pedon for the Batagay megaslump. The unit sampled at A1 represents a separate feature of the lower Ice Complex. The formation of the lower Ice Complex and its preservation remains uncertain. Perhaps, the lower Ice Complex survived at least partly due to the thick and rapidly accumulated lower sand unit that is up to about 25 m thick and is relatively ice-poor. But in this case, what contributed to such an active accumulation of sedimentary cover over the lower Ice Complex? The lower Ice Complex is surrounded by a heterogeneous composition of materials such as sandy loam and loess. The origin and existence of loess to a depth of 58 m suggests the There were ice schlieres with angled directions from the sides of the ice wedges. Above the ice wedges, there was a layer consisting of icy sandy loam of dark gray color which was penetrated by roots, and below was a layer of silty sandy loam of brown color alternating with crushed stones, inclusions of charcoal and narrow strips of peat.
The wedge ice of the upper Ice Complex was pillar-shaped, resembling a solid wall of ice from a distance. Table 1 summarizes descriptions of the cover deposits (samples in profiles B1, C1 and D1) of the upper Ice Complex, mainly composed of sandy and dusty particles with a localized separate newly formed structure (ironing) and individual mechanical inclusions of pebbles and gravel, which formed one soil pedon for the Batagay megaslump. The unit sampled at A1 represents a separate feature of the lower Ice Complex. The formation of the lower Ice Complex and its preservation remains uncertain. Perhaps, the lower Ice Complex survived at least partly due to the thick and rapidly accumulated lower sand unit that is up to about 25 m thick and is relatively ice-poor. But in this case, what contributed to such an active accumulation of sedimentary cover over the lower Ice Complex? The lower Ice Complex is surrounded by a heterogeneous composition of materials such as sandy loam and loess. The origin and existence of loess to a depth of 58 m suggests the accumulation of loess formations occurred as a result of frost weathering and solifluction processes due to periglacial phenomena in the Batagay area. Table 1. Cryolithological characterization of the studied stratigraphic features exposed in the Batagay megaslump.

Carbon Stocks of the Lower and the Upper Ice Complex Units, and the Uppermost Cover
It is necessary to take into account that microorganisms begin to multiply actively in thawing permafrost [20]. When organic matter is oxidized, they remove carbon from the bound state into the atmosphere. The degree of decomposition of soil organic matter depends on the properties of the soil as well as the depth of the profile and community of bacteria, archaea and eukaryotes, which are mostly found in frozen soils. This leads to the intensification of microbial metabolic activity and the possible creation of positive feedback to the conditions of anticipating thawing of the permafrost [34][35][36][37]. The main fraction of OC was concentrated in the active layer topping the upper Ice Complex due to the productivity of this layer, both in terms of greater seasonal thawing and processing of soil organic matter by microbes (Table 2). The resulting determinants are soil warming up to positive temperatures and activation of microbiological processes; as a result, destruction and intensive mineralization of animals and plant residues occur. Variations in the OC content in the active layer varied widely from a minimum value of 0.3% to a maximum of 12.5%. Higher content of OC was prevalent on the slopes of the southern exposure. Consequently, for the northern exposure slope (C1), the average OC value was 2.1%, while for the southern exposure (B1) and the ravine (D1) it was higher with 3.5% and 2.8%, respectively. The deep penetration of heat on the southern exposure slopes causes permafrost thaw and collapse in the warm period of the year compared with the northern exposure slopes, as is noticeable in the picture (see Figure 1). This increases the physical runoff of soil material into the Batagay River, as well as its subsidence on the bottom of the slump.
The properties of the shielding layer are closely related to the active (seasonally thawed) layer. The susceptibility of the seasonally thawed layer to the constant cycle of freezing and thawing determines the presence or absence of the frozen shielding layer at the uppermost part of permafrost. In case of destruction of the vegetation cover, the active layer degrades as a result of soil temperature Land 2020, 9, 305 11 of 16 increase, and the shielding layer loses its usual boundaries and merges with seasonally thawed soil, which leads to a lack of shielding properties leading to permafrost degradation and thaw subsidence. The average OC content in the shielding layer was low if compared to the active layer, reaching only 0.3%. This low OC content is most likely due to the buffer properties of the shielding layer, its low thickness, low moisture content and the paretic rate of chemical reactions. In the permafrost layer, the average OC content for the lower Ice Complex (A1) was 0.65% and for the upper Ice Complex (B1) is 1.04%. There was a concentration increase of 58.5% and 74.0% compared to the shielding layer. Perhaps this is the result of the genetic features of lower cryogenic horizon formation and the water presence in it with ongoing metabolic processes even at temperatures below 0 • C. It is worth noting that in the past, organic residues deposited in combination with heterotrophic microbial activity were capable of carbon remineralizing, despite the negative temperatures of the deposits.
The deposits of the Batagay megaslump had a small amount of carbonate in their composition. This is apparent in low inorganic carbon values, and there were no significant variations over depth. The inorganic carbon distribution was homogeneous both in the active layer and permafrost with systematic low values over the entire profile of sediments. Consequently, the oxide group of minerals took part in the formation of the deposits, for example, quaternary formations of the slump mainly consisted of quartz minerals in the form of sand and they are represented by the loess fraction.
TC content is formed from the two previous carbon fractions (OC and IC); therefore, it basically repeats the nature of their distribution in differentiated (elementary) layers. Table 3 shows the carbon stocks of the studied stratigraphic features of the Batagay megaslump. The enormous amount of organic carbon was stored in the active layer for sections B1 and D1, where 60% of organic carbon was concentrated. At the same thickness of the active layer, the OC in C1 amounted to 9.52 kg m −2 , while the OC of D1 reached 18.03 kg m −2 , which was almost twice as much. Even with a changing of slope exposure, the stocks for these sections in the active layer increased from the northern to the southern exposure by 42%-47% of organic carbon.
The lowest carbon stocks of all studied stratigraphic features were recorded in the shielding layer, reaching 1.3 kg OC m −2 and 0.5 kg IC m −2 . IC stocks in the active layer of section B1 did not differ from the other two profiles (C1 and D1) and ranged from 2 to 2.2 kg m −2 . Carbonates enriched the first meter from the surface in non-permafrost soils. Under permafrost conditions, rather high inorganic carbon data were observed, 1.8 kg m −2 close to the level of the active layer. Section A1 had a higher IC stock of 3.6 kg m −2 . By moving to the lower Ice Complex, the IC reserves were in an an inactive state due to the influence of negative soil temperature, with substances trying to penetrate with soil solutions to the Yedoma horizons, and most likely the reserves remained unchanged for an extended period. Low estimates of inorganic carbon in the deposits of the Batagay slump are the result of low natural carbonate, as well as frosty weathering of rocks. Carbonate formation occurs in soil pores and frost cracks, indicating their predominant dissolution and leaching in cold climate conditions. The upper Ice Complex (B1) shows the distribution of the carbon in the elementary layers and the conjugation of carbon from the thickness and the amount of carbon available for mobilization in them. The maximum reserves of total carbon are contained in the active layer and amount to 18.72 kg m −2 , the minimum in the shielding layer is 1.81 kg m −2 , and in permafrost rocks, it is 10.34 kg m −2 . In total, 30.87 kg m −2 was deposited in a two-meter thickness. The ratio of organic and inorganic carbon in the section varies in scale for the active layer-15, the shielding layer-3 and permafrost deposits-5. The ratio indicators show that organic carbon is the dominant feature in the soil horizon, layer or sediment. The higher this value, the more organic part is contained in a particular element of the Ice Complex. Also, this indicates the ability of organic carbon to be a nutrient medium for microorganisms and a source of CO 2 . Inorganic carbon indicates the level of carbonate in a given horizon. In the upper two meters of the Ice Complex, carbon stocks were 30.87 kg m −2 . In the study area, the dominating forest species is Larix cajanderi. Its phytomass reserves barely exceeded 6.00 kg m −2 , which is formed mainly due to green, lignified aerial parts, roots and moss-lichen complexes. Stocks of dead organic matter are 3.00 kg m −2 , and litter is 2.50 kg m −2 . In the northern taiga zone, they are considered to be the forests with the lowest production, with about 11.50 kg m −2 per year [38]. As a result, carbon accumulation rates are low in the sediments of the upper Ice Complex. The main part of the organic matter is stored on the surface as a humus horizon of the soil, which undergoes degradation and weak decomposition, and the mineralization rate in such natural conditions remains low.
There is no doubt that the continental [1,6] genesis of the formation of the lower Ice Complex (A1) and its archaic past, compared with the upper Ice Complex (B1), affected the loss, impossibility of renewal and accumulation of plant organic matter. The cycling of the biogenic elements was interrupted by the rapid accumulation of sandy loam layers and sediment deposits [5] and most likely by the removal of organic layer beyond the meadow steppes [39] due to aeolian processes. We estimated the amount of ancient buried carbon at 17.66 kg m −2 for a thickness of 1.2 m.
Permafrost degradation caused by global warming is not always a catastrophic consequence for the ecumene. To justify this point of view, we give a few well-founded examples. Siberian larch forests occupy 40% of all coniferous wood in Russia, which develop on permafrost. The amount of accessible nutrients for trees increases with permafrost degradation, which causes an increase in their biomass [40]. An increase in the average annual temperature and rainfall contributes to the advancement of natural zones deep into the north, expanding the agricultural area and changing the structure of the soil cover with an increasing share of organic matter due to an increase in the productivity of vegetation [41].
The permafrost zone has changed repeatedly, giving way to periods of warming and cooling, which favoured the formation and accumulation of permafrost in the past. With a certain degree of probability, it can be assumed that the multi-tiered features of the Batagay megaslump and their heterogeneous origin are the results of the participation of several geological factors: (a) prolonged permafrost conditions of the territory with cryogenic and exogenous processes occurring in it; (b) the transfer of aeolian material and its accumulation; (c) the impact of the Quaternary climate and paleoclimate; (d) the movement and deformation as a result of the tectonic influence of the territory. That means the Adycha-Taryn general snap has or had an impact on the formation of the Batagay megaslump in conjunction with the other listed predictors. This is because the formation of such a large-scale current geocryological object in the North of Yakutia could not have taken place without the participation of intense seismotectonic conditions of the territory.
The anthropogenic activity was also one of the significant reasons for the formation of the megaslump. The economic and industrial development of the northern territories has never left an integral ecosystem behind, especially in such vulnerable natural conditions as boreal forests or tundra landscapes, which undoubtedly provoked the degradation of the Batagay Ice Complex through the destruction of vegetation and soil cover. Cutting down forests for timber for consumer needs, including the destabilizing pyrogenic situation at the site, has also had an impact.

Conclusions
In Yedoma deposits, carbon stocks are at a low level, which is typical of Northern Yakutia. The potential for the formation of organic matter in the local natural conditions of the Yana Upland is limited by climatic features and primary production of the forest. First of all, the soft parts of plants and carbon fractions that are readily accessible to microorganisms decompose to the final mineralization product, and the bulk of the organic matter reserves are stored on the soil surface, forest litter and in the first meter of the Ice Complex. As long as these reserves remain intact, the organic system is in a stable state, and there is a dynamic equilibrium between the release of carbon and its decrease in the atmosphere.
Otherwise, with the vast Ice Complex degradation on the Yansky plateau, it is possible that new large-scale sites like the Batagay megaslump will appear. This is not just about the study area, but also about the whole cryolithozone of the planet. As a result, greenhouse gas emissions will increase, and the reserves of conserved organic matter will be significantly depleted, which will lead to irreversible processes of the biological system and climate changes.
At the present time the Batagay megaslump of North-Eastern Siberia remains the only large geocryological object on Earth which is demonstrating the impact of anthropogenic and natural cause-and-effect connections on permafrost. It is also an obvious example of how permafrost is unprotected and prone to rapid degradation in a very short period of time if law enforcement measures to protect and prevent the melting of permafrost are not taken.