Assessments of Organic Carbon Stabilization Using the Spectroscopic Characteristics of Humic Acids Separated from Soils of the Lena River Delta

: In the Arctic zone, where up to 1024 × 1013 kg of organic matter is stored in permafrost-affected soils, soil organic matter consists of about 50% humic substances. Based on the analysis of the molecular composition of humic acids, we assessed the processes of accumulation of the key structural fragments, their transformations and the stabilization rates of carbon pools in soils in general. The landscape of the Lena River delta is the largest storage of stabilized organic matter in the Arctic. There is active accumulation and deposition of a signiﬁcant amount of soil organic carbon from terrestrial ecosystems in a permafrost state. Under ongoing climate change, carbon emission ﬂuxes into the atmosphere are estimated to be higher than the sequestration and storing of carbon compounds. Thus, investigation of soil organic matter stabilization mechanisms and rates is quite an urgent topic regarding polar soils. For study of molecular elemental composition, humic acids were separated from the soils of the Lena River delta. Key structural fragments of humic matter were identiﬁed and quantiﬁed by CP/MAS 13 C NMR spectroscopy: carboxyl (–COOR); carbonyl (–C=O); CH 3– ; CH 2– ; CH-aliphatic; –C-OR alcohols, esters and carbohydrates; and the phenolic (Ar-OH), quinone (Ar = O) and aromatic (Ar–) groups as benchmark Cryosols of the Lena delta river terrestrial ecosystem. Under the conditions of thermodynamic evolutionary selection, during the change between the dry and wet seasons, up to 41% of aromatic and carboxyl fragments accumulated in humic acids. Data obtained showed that three main groups of carbon played the most important role in soil organic matter stabilization, namely C, H-alkyls ((CH 2 )n/CH/C and CH 3 ), aromatic compounds (C-C/C-H, C-O) and an OCH group (OCH/OCq). The variations of these carbon species’ content in separated humics, with special reference to soil–permafrost organic proﬁles’ recalcitrance in the current environment, is discussed.


Introduction
Soil organic matter (SOM) is a product that accumulates in soil after the partial decomposition of various types of materials derived from microorganisms and plant residues. It constitutes a key element of the global carbon cycle in the atmosphere, vegetation, soil, rivers and oceans [1][2][3][4]. Soil organic matter supports key functions of soil and ecosystem services, as it is critical for stabilizing soil structure, retaining and releasing plant nutrients and allowing water to penetrate and be stored in the soil. Soils provide a wide range of ecosystem services (benefits) for human beings. These include water purification and pollutant reduction, climate regulation, nutrient cycling, providing a habitat for various living organisms, being a source of medicinal substances and building materials as well as a basis for infrastructure, and also participating in the sequestration of organic carbon in soil [1]. The most important ecosystem service of the soils of the Arctic belt is the sequestration and deposition of SOM within the permafrost. Nowadays, up to 496 × 10 13 kg of SOM are buried in soil in a permafrost state [5,6].
Under the process of climate change, the emission of carbon into the atmosphere can be greater than its sequestration in the soil and permafrost. This could affect the global carbon cycle and accelerate climate change on the planet. Ecosystems with the highest SOM content (meadows, swamps, and permafrost) can make the greatest contributions to climate change. Therefore, to assess the interactions of the biogeochemical carbon cycle with climate, it is important to characterize the areas affected by permafrost as well as the soils' organic carbon stocks and the rate of the soils' transformation. An important direction in the study of the transformation of SOM is the mechanism of stabilization of organic matter; its resistance to biodegradation by soil microorganisms, plants and fungi; and its chemical decomposition [7]. The permafrost zone covers an area of more than 8.6 million km 2 , which is about 27% of the total land area above 50 • N. It accumulates a huge amount of organic carbon, and therefore it is considered one of the most important elements of the cryosphere. There, carbon accumulates in soils in huge amounts due to low temperatures, which lead to slow degradation of SOM [7][8][9][10][11]. During the transformation of SOM, greenhouse gases (carbon dioxide and methane) are emitted into the atmosphere. If this happens too quickly, soils can contribute to the climate crisis on our planet. On the other hand, many soils have the potential to increase SOM stocks, thus mitigating climate change by reducing CO 2 emissions [10,12,13].
To study the molecular composition of SOM, the analysis of humic acids (HAs), which are formed in soils during the humification of organic residues, has been used. HAs are heterogeneous systems of different (high and low) molecular weights; they are condensed compounds formed as a result of the decomposition of plant and animal remains in terrestrial and aquatic ecosystems [14,15].The question of the structure of HAs still relevant based on the work carried out in the gel chromatography by Piccolo et al. [16]; it was concluded that the HAs were a supramolecular system; they were a set of molecules with relatively low molecular weight (up to several hundred D) combined into a kind of micelle-like structure of hydrogen and hydrophobic bonds. However, this technique allowed us the right to doubt the correctness of the conclusions. According to the work of Perminova [17], the use of a strong mineral acid set to pH 2 in a HA solution changed its ionic strength by only 0.026 M; when using a weak organic acid (acetic acid), the content of substances in the system increased by 4.3 M. Moreover, data have been given on the existence of macromolecules, which are transformed biomolecules (derivatives of lipids, lignin, carbohydrates and proteins), in the composition of HAs, which does not contradict the "macromolecular model". Thus, modern data on the structure of HAs indicate that it is rather a "molecular ensemble" of macromolecules in which molecules with different molecular weights are present [14,15].
Climatic parameters, precursors of humification and local position in the landscape determine the diversity of the composition and properties of HAs [9,[18][19][20][21][22]. Various methods are used to study the HAs in soil organic matter. There are also many methods for determining the composition of HAs. Fourier transform infrared spectroscopy, ultraviolet and visible spectroscopy, cross polarization, magic angle rotation, molecular fluorescence spectroscopy and electron paramagnetic resonance spectroscopy are all useful and often used in SOM studies [23,24].
The method of electron paramagnetic resonance (EPR) is used to determine the structure of SOM, which consists of free radicals. The free radicals participate in the process of polymerization humus formation, and free radical activity is a fundamental property. In a number of works, a correlation between the content of free radicals and the degree of humification was noted. Thus, the thesis was put forward that the condensation of HA molecules led to the formation of polycondensed systems that were biochemically resistant to biodegradation. This method was highly sensitive and did not destroy the studied HA material [24,25]. Vibrational (IR) spectroscopy is a method used to determine the structural fragments of HAs for the analysis of the structure of molecules. The method is based on weakening the intensity of light transmitted through the test substance. This method can be used to determine aliphatic and aromatic structures. Nevertheless, this method, in comparison with 13 C-NMR, has lower accuracy due to the fact that when determining any group of structural fragments in the HA molecule, the group is not isolated from the rest of the molecule, which leads to a change in the intensity of the transmitted light. Taking into account the high heterogeneity of HAs molecules, absorption bands in the IR spectra can form, which cannot be attributed to vibrations related to any structural fragments [26].
Densitometric fractionation of SOM is a fairly simple and informative way to isolate chemically and biologically heterogeneous structures of organic matter. The method is based on the structural features of SOM; when using this method of analysis, organic matter can be divided into three groups (free, occluded and organo-mineral). Their differences are associated with their composition and resistance to biodegradation [27,28].
One-dimensional (1D) solid-state 13 C NMR spectra provide some structural information about HAs, including the quantification of various types of chemical groups [9,18,20,29,30]. The advantage of using the method of one-dimensional NMR spectroscopy is the ability to quantitatively determine the content of groups of individual and structural fragments in HA molecules. This method is also used to assess changes in SOM during decomposition and humification. Nowadays, studies of the quality of SOM from the Arctic belt have revealed the generalized, slightly degraded nature of organic molecules that retain most of the chemical nature of their precursor humification material (plants and soil living organisms) due to a weak humification process under cryogenic processes [19]. However, using a 1D structure, the differences among the atoms cannot be characterized in a prominent way due to wide lines and the diversity of the structures of functional groups [31][32][33]. Using (2D) NMR techniques, such peak overlap can be eliminated, and the chemical structure corresponding to a given chemical shifts can be investigated more specifically and provide additional structural information about dipolar interaction among the unprotonated carbon atoms [31].
For further study of the fundamental processes of humification and structural fragments in polar regions, it is necessary to use modern, high-precision instrumental methods. Analysis of molecular composition is necessary to understand the fundamental processes of SOM transformation and the structures of the natural, high-molecular compounds of humic acids in soils affected by permafrost [9,20,30,34].
In this way, the aim of this work was to use the composition of SOM to deduce trends in the decomposition rates of HAs in Arctic sediments. Thus, the following tasks were set: -Investigate the elemental composition of HAs separated from soils of the Lena River delta; -Determine the molecular composition of HAs separated from soils of the Lena River delta; -Estimate the rate of stabilization of organic matter in the studied soils and buried organic remains.

The Study Site
The Lena River delta is the largest northern delta in the world; it is located in the Arctic zone and has an area of about 30,000 km 2 . The delta was formed as a result of river activity: sediment removal, erosion, abrasion under the influence of sea level fluctuations and movement of the earth's crust [35] (Figure 1). The Lena River delta can be defined by three different geomorphological terraces with active flooded levels. The active flooded level is the first terrace, which is the youngest part of the river delta. The first terrace was formed during the Middle Holocene and is mainly represented in the eastern part of the Lena delta. The second terrace was formed between the Late Pleistocene and the Early Holocene. The polygonal relief is less pronounced, and a large number of thermokarst lakes are typical of this area. The third terrace is the oldest terrace in the Lena delta. It consists of fine-grained sand that is rich in organic matter; it is icy material that accumulated in front of the Chekanovsky and Kharulakh Ridges.
The Lena River delta is located in a zone with an arctic continental climate. The average annual air temperature is −13 °C; the average January temperature is −32 °C and the average July temperature is −6.5 °C. Annual precipitation is 190 mm. Most of the land is characterized by the presence of permafrost at a depth of about 1 m. The depth of the seasonally thawed layer is different; in depressions it can reach 30 cm at the end of August, and on the sand materials, up to 1 m.

Sampling Procedure
The sampling and classification procedure of soils and soil horizons was carried out according to the standard procedure [36,37]. In the course of the work, samples were taken from the upper organo-mineral as well as the lower horizons in which suprapermafrost accumulation of organic matter occurred.
The key areas of research were soils from the Kurungnah Islands (third terrace), soils from Jipyries Island (second terrace) and soils from the Chekanovsky Ridge (Table 1). The Lena River delta can be defined by three different geomorphological terraces with active flooded levels. The active flooded level is the first terrace, which is the youngest part of the river delta. The first terrace was formed during the Middle Holocene and is mainly represented in the eastern part of the Lena delta. The second terrace was formed between the Late Pleistocene and the Early Holocene. The polygonal relief is less pronounced, and a large number of thermokarst lakes are typical of this area. The third terrace is the oldest terrace in the Lena delta. It consists of fine-grained sand that is rich in organic matter; it is icy material that accumulated in front of the Chekanovsky and Kharulakh Ridges.
The Lena River delta is located in a zone with an arctic continental climate. The average annual air temperature is −13 • C; the average January temperature is −32 • C and the average July temperature is −6.5 • C. Annual precipitation is 190 mm. Most of the land is characterized by the presence of permafrost at a depth of about 1 m. The depth of the seasonally thawed layer is different; in depressions it can reach 30 cm at the end of August, and on the sand materials, up to 1 m.

Sampling Procedure
The sampling and classification procedure of soils and soil horizons was carried out according to the standard procedure [36,37]. In the course of the work, samples were taken from the upper organo-mineral as well as the lower horizons in which suprapermafrost accumulation of organic matter occurred.
The key areas of research were soils from the Kurungnah Islands (third terrace), soils from Jipyries Island (second terrace) and soils from the Chekanovsky Ridge (Table 1). Soil samples were taken during the 2019 field season. Samples were taken from the superficial soil horizons as well as from the lower horizons located on the permafrost table in order to assess the degree of stabilization of organic matter and study the evolutionary selection of HA molecules within the soil profile. A sample from the Chekanovsky Ridge was taken from a wind shelter inside rocks ( Figure 2).

Laboratory Analysis
Soil samples were air-dried (24 h, 20 °C), ground and passed through a 2 mm sieve. Chemical analyses were performed using classical methods: C and N content were determined using an element analyzer (EA3028-HT EuroVector, Pavia PV, Italy), their pH in water (soil-dissolvent ratios were 1:2.5 in the case of mineral horizons and 1:25 in the case of organo-mineral horizons) suspensions using a pH meter (pH-150M, Teplopribor, Moscow, Russia). The particle-size distribution analysis was carried out according to the Kachinsky "wet sedimentation" method, which is the Russian analogue of analysis by Bowman and Hutka [38].

Laboratory Analysis
Soil samples were air-dried (24 h, 20 • C), ground and passed through a 2 mm sieve. Chemical analyses were performed using classical methods: C and N content were determined using an element analyzer (EA3028-HT EuroVector, Pavia PV, Italy), their pH in water (soil-dissolvent ratios were 1:2.5 in the case of mineral horizons and 1:25 in the case of organo-mineral horizons) suspensions using a pH meter (pH-150M, Teplopribor, Moscow, Russia). The particle-size distribution analysis was carried out according to the Kachinsky "wet sedimentation" method, which is the Russian analogue of analysis by Bowman and Hutka [38].

Protocol for the Extraction of Humic Acids from Soils
Humic acids were extracted according to a published IHSS protocol [39] in modification by Chukov [15] that involved: -Soil sieving: We manually broke up a portion of the dried soil sample into powder or small portions using a mortar and pestle and sieved the ground sample through 2 mm mesh; After centrifugation and washing, the HA gel was placed in bags made of dialysis cellophane and placed in large containers with distilled water. During the first 3 days, the water was changed every day, and then it was change after 2 days. Usually, it took 7-10 days to completely remove excess sodium sulfate. The completeness of removal was controlled by the qualitative reaction to sulfates in the water flowing down from the bag; -Drying of HA preparations: HAs preparations from dialysis bags were transferred into Petri dishes or small crystallizers, after which they were dried in a vacuum oven over containers with dry CaCl 2 . Further, the HA preparations were transferred into weighing bottles and, if necessary, dried in a desiccator over P 2 O 5 . The ash content of the HA preparations obtained by this method did not exceed 5% and significantly depended on the speed and duration of the first centrifugation.

Elemental Composition of HAs
The elemental composition of an HA is the percentage of the elements C, H, N and O. For graphical analysis of the elemental composition, we used the van Krevelen diagram [40], using the H/C and O/C ratios to determine the direction of the transformation processes of various organic compounds into natural conditions. The elemental composition was corrected for weight moisture and ash content. The oxygen content was calculated from the differences in the masses of the whole samples and the gravimetric concentrations of C, N, H and ash. The determination was carried out in the same way on an elemental analyzer (EA3028-HT EuroVector, Pavia PV, Italy). Ethylenediaminetetraacetic acid disodium salt dehydrate (EDTA) was used as a reference standard in the analysis (CAS Number: 6381-92-6).

Molecular Composition of HAs
Solid-state spectra of humic acids were determined by CP/MAS 13 C-NMR spectroscopy on a NMR spectrometer (Bruker Avance 500, Billerica, MA, USA) in a 3.2 mm ZrO 2 rotor. The magic angle rotation speed was 20 kHz, and the cross-polarization nutation frequency u1/2p 1/4 was 62.5 kHz. The repetition delay was 3 s. The number of scans was 6500-35,000. Temperature was 294.5 • C, according to a 4 mm MAS BB/1H H13925/0005 probe.
Molecular weight was calculated using the minimum molecular weight; its value corresponds to the simplest formula of humic acid that is calculated from the elemental composition and determines the lower limit of possible molecular weights [41,42].

Statistical Survey
The statistical analysis was performed in Paleontological Statistics (PAST) software. Principal component analysis The method is based on projecting a multidimensional dataset of up to several dimensions in such a way as to preserve as much variation as possible and to facilitate data visualization. In addition, the axes of maximum variance (principal components) can be identified and interpreted.
Spearman's correlation coefficient The method is based on measuring the linear relationship between random variables. Hierarchical clustering: Ward's method. This analysis uses a method of variance to estimate the distances between clusters. The method minimizes the sum of squares (SS) for any two (hypothetical) clusters that could be formed at each step.

The Physico-Chemical Characteristics of Study Soils
The studied soils were characterized by slightly acidic and neutral pH reactions ( Table 2). This was due to the erosion of carbonate rocks in the middle course of the Lena River; the presence of carbonates was also noted in the ridges adjacent to the delta [43]. In the studied soils, relatively little SOM accumulated, with the highest content noted in the wind shelter site, which was associated with more favorable conditions for humification [44]. The enrichment of organic matter with nitrogen is low, according to Orlov [45]. The studied soils were dominated by the sandy fraction, which was associated with the activity of the river and the features of the formation of the delta complex. The low content of the clay fraction was associated with the relatively high speed of the river; in such conditions, mainly large particles (sand, dust) settle [46].

Elemental Composition
The elemental composition of an HA is the most important indicator that determines the progress of humification and oxidation as well as the HA's degree of condensation. The characteristic features of HAs formed under arctic conditions, and especially in soils affected by permafrost, are relatively high H and low O concentrations compared to boreal and subboreal soils (Table 3) [30]. The studied soils accumulated up to 52% of carbon in the form of HAs, and the highest content was noted in the organo-mineral horizon in the soils of a drained lake (Y 1-3). This site, a drained lake, was quite different from the typical conditions prevailing on this island (zonal tundra). In the studied soils, there was a thick (up to 25 cm) organo-mineral horizon, which formed under the vascular plants. The processes of peat accumulation and active cryogenesis dominated in the soils located outside the drained lake. A characteristic feature of the soils of a drained lake is the presence of buried, medium-decomposed organic remains that accumulate on the permafrost table. Different dynamics of carbon accumulation in the composition of HAs were revealed within the three studied pingos (intra-permafrost ice-cored hill) in the area of a drained lake. In the Y-1 profile, the carbon content decreased with depth and the oxygen content increased. That dynamic could apparently be associated with an increase in oxygen-containing fragments (O-alkyl and anomers) due to an increase in soil hydromorphism (the result of temporary or permanent waterlogging conditions of the soil profile or soil horizon, wherein the amount of moisture exceeded 70-80% of the total soil moisture capacity) in this pingo. This was also reflected in the closer occurrence of the permafrost table and the height of the pingo (1.5 m). In pingo Y-2, the same carbon content was observed in the upper and lower organo-mineral horizons, with a slight increase in the oxygen content of the buried organic material (from 41% to 43%). In pingo Y-3, there was an increase in carbon content (from 47% to 50%) and a decrease in oxygen content. The increase in the carbon content in the buried organic horizon may have been associated with the processes of humification and humus accumulation at the boundary with the permafrost. The different dynamics of carbon in the studied organic residues are explained by the high variability of soil formation conditions in this region, the active development of cryogenic landforms (pingos), humus accumulation and soil hydromorphism [47][48][49][50][51].
The H/Cmod ratio is an indicator of the stability of HAs in soils [45]. An increase in the proportion of hydrogen in the HA structure indicates the significant development of side chains with C-H, CH 2 and CH 3 groups (aliphatic compounds) [47,52]. The lower this indicator, the greater is the process of condensation of monomers into high-molecular substances. In general, this shows a low variation, which is associated with the homogeneous conditions of the formation of organic matter in the Lena River delta. The smallest indicators were noted in the samples of the drained lake (Soil ID numbers1 and 6), which indicated the process of condensation in HA molecules. The greatest value was noted in sample number 7 on the island of Jipyries, the second terrace. The most active cryogenic processes in soils were noted here. This was noted in the morphological characteristics of the soil profile (cryoturbation and suprapermafrost accumulation of organic matter) [53,54].
For a graphical representation of the elemental composition of HAs from soils, the method by Kleinhempel and van Krevelen was used (Figure 3). The method was based on the construction of H/Cmod and O/C diagrams and served as a method to demonstrate the contributions of oxidation and condensation to changes in the elemental composition of HAs.
indicated the process of condensation in HA molecules. The greatest value was noted in sample number 7 on the island of Jipyries, the second terrace. The most active cryogenic processes in soils were noted here. This was noted in the morphological characteristics of the soil profile (cryoturbation and suprapermafrost accumulation of organic matter) [53,54].
For a graphical representation of the elemental composition of HAs from soils, the method by Kleinhempel and van Krevelen was used (Figure 3). The method was based on the construction of H/Cmod and O/C diagrams and served as a method to demonstrate the contributions of oxidation and condensation to changes in the elemental composition of HAs.  Table 1. From the data obtained, the integral index of H/Cmod and O/C in most of the studied HAs was relatively low, which indicated a low content of oxygen-containing fragments in HAs. As mentioned above, high values of H/Cmod and O/C were observed in the sample from the second terrace (Jipyries Island) and indicated that the sample was formed in an environment with a low level of microbiological activity, which contributed to better preservation of the carbohydrate and amino acid fragments of HAs. A decrease in H/Cmod indicated the accumulation of fragments more resistant to biodegradation in the composition of HAs [55]. The most condensed molecules were formed in buried organic remains, which might have been related to the evolutionary selection of stable HA molecules [56,57].

Molecular Composition of HAs Separated from the Soils of the Lena River Delta
Numerous fragments of molecules were identified by CP/MAS 13 C NMR spectroscopy: carboxyl (-COOR); carbonyl (-C=O); CH3-, CH2-, CH-AL; -C-OR alcohols, esters and carbohydrates; and phenolic (Ar-OH), quinone (Ar=O) and aromatic (Ar-) groups, which indicated the great complexity of the structure of HA [20]. Signals from non-polar  Table 1. From the data obtained, the integral index of H/Cmod and O/C in most of the studied HAs was relatively low, which indicated a low content of oxygen-containing fragments in HAs. As mentioned above, high values of H/Cmod and O/C were observed in the sample from the second terrace (Jipyries Island) and indicated that the sample was formed in an environment with a low level of microbiological activity, which contributed to better preservation of the carbohydrate and amino acid fragments of HAs. A decrease in H/Cmod indicated the accumulation of fragments more resistant to biodegradation in the composition of HAs [55]. The most condensed molecules were formed in buried organic remains, which might have been related to the evolutionary selection of stable HA molecules [56,57].
The obtained spectra of the studied HAs separated from the soil of the Lena River delta are shown in Figure 4.   (Table 4).   (Table 4). The studied HAs accumulated up to 71% of aliphatic fragments. The highest content was noted in the organic matter accumulating at the permafrost table in the profile from the Jipiries Island (second terrace). The predominance of aliphatic fragments indicated a deficiency of vascular plants or the low «maturation» of humic substances in the soil. The predominance of aliphatic structures was characteristic of humic substances formed under reduction conditions including aqueous humic substances [18,19,45,48].
In the studied HAs isolated from pingos, weak dynamics in the molecular composition could be noted. In Y-1, there was a slight increase in aromatic fragments (aromatic C-O) in the buried remains. This was apparently related to the evolutionary selection of stable macromolecules over time. This was also noted in Y-2. Normally, buried substances have a more stable structure, which is noted in a number of other investigated works. The soils of the Lena River delta are formed in soft climatic conditions in contrast to the soils that form in continental parts at the same latitude.
In the samples from the island of Jipyries, negative dynamics were observed, with a decrease in aromatic compounds with depth [58]. This was apparently due to the low content of lignin-containing fragments in the precursors of humification [15,20]. Due to the hydrogenation process, the C-C chains could be broken, and the aromatic structures in the composition of humic acids could degrade [34,45,51].
The highest contents of aromatic structures were observed in the sample of buried material in Y-2 as well as in the sample from the wind shelter (number 9). The formation of such sites apparently led to the accumulation of high-molecular compounds in the structures of HAs. Under the conditions of wind shelters, the formation of soil horizons was less influenced by climatic features and the microclimate was more stable, which had a positive effect regarding an increase in the content of aromatic structures as well as resistance to biodegradation.

The Stabilization of Organic Matter Based on 13 C NMR Spectroscopy
The studied samples accumulated up to 40% aromatic structures (aromatic C and -COOR fragments). This content was relatively high in relation to the soils in the continental part of the Arctic belt. Soil cryogenesis, climatic features and humification precursors have a rather strong effect on the molecular composition of HAs. Low microbiological activity, a relatively high level of hydromorphism in Cryosols and the predominance of moss and lichens in the plant composition lead to the accumulation of aliphatic fragments. In the studied soils, the vegetation cover was dominated by vascular plants, in which up to 30% of lignin-containing fragments were accumulated [59]. Under conditions of demitylation, during the decomposition of lignin, the HA structure accumulates a carboxyl group, which leads to the stabilization of organic matter [8,9,21,23,30,32,34,52,58].
Under the conditions of the transformation of mosses and lichens, carbohydrates (mono-and oligosaccharides, cellulose, fats, alcohols) come into the soil [26,45,56]. During their transformation, aliphatic acids with a content of fatty acids and paraffins are predominantly formed. Plant residues are the sources of various chemical structures, while the formation of aliphatic and aromatic structures is largely associated with the activity of the soil microbiota. The increase in the proportion of aliphatic fragments is associated with a high content of carbohydrates in plant residues as well as the possibility of producing these fragments by the soil microbiota itself. These substances are a larger part of the biochemical cycle and are consumed faster by living organisms [18,49,60,61].
Lignin makes the greatest contribution to the aromaticity of molecules [59]. This component is relatively resistant to biodegradation. It is an irregular polymer with branched macromolecules built from residues substituted by phenol alcohols. This leads to the possibility of the appearance of a large number of various products of its decomposition. Depending on the composition of plant residues, methoxyl groups of compounds as well as various alcohols can form. In organogenic and peat horizons, various resins and waxes, which consist of esters, mono and dihydric alcohols, accumulate to a greater extent [7,18,21,42,45,49,60,61].
In our study, a change in the molecular composition of soil profiles was noted. Under the conditions of transformation of various compounds (carbohydrates, arenas), the condensation of molecules and the formation of various compounds occurred. The longer this process, the more stable molecules accumulated in the HAs, which led to the stabilization of the organic material. HAs are compounds of variable composition. Therefore, their transformation/condensation over time is associated with the evolutionary selection of fragments that are more resistant to biodegradation.
To standardize the quantitative characteristics of HA molecules, the following parameters were used: the degree of decomposition of organic matter (C-alkyl/O-alkyl) and the integral index of HA hydrophobicity (AL h, r + AR h, r) ( Figure 5). process, the more stable molecules accumulated in the HAs, which led to the stabilization of the organic material. HAs are compounds of variable composition. Therefore, their transformation/condensation over time is associated with the evolutionary selection of fragments that are more resistant to biodegradation.
To standardize the quantitative characteristics of HA molecules, the following parameters were used: the degree of decomposition of organic matter (C-alkyl/O-alkyl) and the integral index of HA hydrophobicity (AL h, r + AR h, r) ( Figure 5). From the data obtained, we noted that the most humified remains were formed in buried organic remains (numbers 6,2,8). The increased proportion of aromatic substances was associated with local position in the relief, the duration of the humification process, wind conditions, hydromorphism of the territory and the precursors of humification (vegetation) [44]. According to degradation hypotheses, stable plant and microbial biopolymers, during transformation, form aromatic benzene rings and labile aliphatic fragments are almost completely mineralized. Soil contains a wide range of various enzymes; various chemical reactions take place here at the same time. The most rapidly mineralized residues are carbohydrates and amino acids [24,25,31,32,34,62,63] and then oligo-and polysaccharides and lignin. The most stable compounds are organic substances of a "secondary" nature, which are products that are formed when the initial substances change (condensed macromolecules). The high stability of such compounds is explained by the fact that they have already passed the stage of biochemical transformations and are now capable of only certain limited reactions, for example, with molecular oxygen. On the other hand, these substances have acquired structural features that different from typical bioorganic molecules. As a result, they gradually become less susceptible to the effects of enzymes, narrowly specialized for strictly defined substrates. Thus, these substances represent the first forms of thermodynamically stable organic carbon compounds outside living organisms under the conditions of the soil cover. HAs are able to persist in soils for quite From the data obtained, we noted that the most humified remains were formed in buried organic remains (numbers 6,2,8). The increased proportion of aromatic substances was associated with local position in the relief, the duration of the humification process, wind conditions, hydromorphism of the territory and the precursors of humification (vegetation) [44]. According to degradation hypotheses, stable plant and microbial biopolymers, during transformation, form aromatic benzene rings and labile aliphatic fragments are almost completely mineralized. Soil contains a wide range of various enzymes; various chemical reactions take place here at the same time. The most rapidly mineralized residues are carbohydrates and amino acids [24,25,31,32,34,62,63] and then oligo-and polysaccharides and lignin. The most stable compounds are organic substances of a "secondary" nature, which are products that are formed when the initial substances change (condensed macromolecules). The high stability of such compounds is explained by the fact that they have already passed the stage of biochemical transformations and are now capable of only certain limited reactions, for example, with molecular oxygen. On the other hand, these substances have acquired structural features that different from typical bioorganic molecules. As a result, they gradually become less susceptible to the effects of enzymes, narrowly specialized for strictly defined substrates. Thus, these substances represent the first forms of thermodynamically stable organic carbon compounds outside living organisms under the conditions of the soil cover. HAs are able to persist in soils for quite a long time; the age of such compounds can reach up to 35 thousand years. They give soil stability and buffering capacity as well as a certain biochemical background [30,45,60,61].
Long-term processes of humification under the influence of cryogenesis and low microbiological activity contribute to the formation of thermodynamically stable compounds and lead to the storage of organic carbon in the permafrost state. An increase in the degree of hydrophobicity in HA molecules is also associated with the degradation of organic matter. At a high level of hydrophobicity of molecules, hydrogen does not attach to the double and triple bonds of carbon and prevents the breaking of these bonds and the degradation of organic compounds in the structure of HAs. Buried organic matter serves as a reliable indicator of the processes taking place in soils because it is able to preserve the molecular structure [45,55,64].
Statistical analysis using Spearman's correlation and analysis of principal components in the composition of HAs was performed (Figures 6 and 7). crobiological activity contribute to the formation of thermodynamically stable compounds and lead to the storage of organic carbon in the permafrost state. An increase in the degree of hydrophobicity in HA molecules is also associated with the degradation of organic matter. At a high level of hydrophobicity of molecules, hydrogen does not attach to the double and triple bonds of carbon and prevents the breaking of these bonds and the degradation of organic compounds in the structure of HAs. Buried organic matter serves as a reliable indicator of the processes taking place in soils because it is able to preserve the molecular structure [45,55,64].
Statistical analysis using Spearman's correlation and analysis of principal components in the composition of HAs was performed (Figures 6 and 7). High correlation was observed in the nitrogen with carboxyl group (r = 0.9), hydrogen with N-alkyl (r = 1) and hydrogen with quinone (r = 0.89). As we said earlier, the composition of HAs from cryoconites varies greatly depending on biological (biochemical) processes, which is reflected in the elemental and qualitative composition of HAs.   Based on the statistical analysis, it could be concluded that there was a high correlation between the studied samples associated with aliphatic fragments (C,H-Alkyl PC1 (89.5%)) and aromatic fragments PC2 (7.3%). Thus, the formation of long aliphatic chains presented in lipids (fatty acids, paraffins), which were the result of the decomposition of moss/lichen residues and the activity of microorganisms. The accumulation of aromatic compounds was associated with the long-term transformation of lignin-containing components in the composition of vascular plants. In general, C,H-Alkyl and carboxylic acids corresponded to the composition of HAs separated from the studied soils. To identify dif- High correlation was observed in the nitrogen with carboxyl group (r = 0.9), hydrogen with N-alkyl (r = 1) and hydrogen with quinone (r = 0.89). As we said earlier, the composition of HAs from cryoconites varies greatly depending on biological (biochemical) processes, which is reflected in the elemental and qualitative composition of HAs.
Based on the statistical analysis, it could be concluded that there was a high correlation between the studied samples associated with aliphatic fragments (C,H-Alkyl PC1 (89.5%)) and aromatic fragments PC2 (7.3%). Thus, the formation of long aliphatic chains presented in lipids (fatty acids, paraffins), which were the result of the decomposition of moss/lichen residues and the activity of microorganisms. The accumulation of aromatic compounds was associated with the long-term transformation of lignin-containing components in the composition of vascular plants. In general, C,H-Alkyl and carboxylic acids corresponded to the composition of HAs separated from the studied soils. To identify differences between the studied samples, cluster analysis (Ward's method) was carried out (Figure 8). Based on the data obtained, we could distinguish two groups that were closest in terms of values (Soil ID 4, 9 and Soil ID 1-3, 5-8). The greatest differences were observed between the studied groups.
Our data corresponded with previously published materials from scientists working in the Arctic belt [9,45,47,52,55,64,65]. The increase of aliphatic HA fragments was associated with the specific composition of the vegetation cover and the microbiological composition and activity of the soil as well as climatic conditions [30]. The tundra and taiga zones are characterized by the predominance of moss and lichen vegetation, which is a source of carbohydrates and various lipids. Depending on the climatic zone (tundra, taiga, etc.), the dynamics of the contents of aliphatic and aromatic fragments were noted; in more southern areas, an increase in aromatic fragments occurred, which was associated with a change in vegetation cover and in the climate of the territory [34,48,51,52]. Thus, the annual changes in climatic parameters, cryogenic processes and precursors of humification determined the composition of HAs in the study area. The predominance of moss/lichen communities contributed to the formation of variable aliphatic chains in HA molecules. The replacement of plant communities with vascular plants and the alternation of wet and dry seasons promoted the condensation and demitylation of aromatic and carboxyl fragments of HAs, which were associated with the resistance of organic material to biodegra- Based on the data obtained, we could distinguish two groups that were closest in terms of values (Soil ID 4, 9 and Soil ID 1-3, 5-8). The greatest differences were observed between the studied groups.
Our data corresponded with previously published materials from scientists working in the Arctic belt [9,45,47,52,55,64,65]. The increase of aliphatic HA fragments was associated with the specific composition of the vegetation cover and the microbiological composition and activity of the soil as well as climatic conditions [30]. The tundra and taiga zones are characterized by the predominance of moss and lichen vegetation, which is a source of carbohydrates and various lipids. Depending on the climatic zone (tundra, taiga, etc.), the dynamics of the contents of aliphatic and aromatic fragments were noted; in more southern areas, an increase in aromatic fragments occurred, which was associated with a change in vegetation cover and in the climate of the territory [34,48,51,52]. Thus, the annual changes in climatic parameters, cryogenic processes and precursors of humification determined the composition of HAs in the study area. The predominance of moss/lichen communities contributed to the formation of variable aliphatic chains in HA molecules. The replacement of plant communities with vascular plants and the alternation of wet and dry seasons promoted the condensation and demitylation of aromatic and carboxyl fragments of HAs, which were associated with the resistance of organic material to biodegradation [18,20,23,24,31,34,47,52,60].

Conclusions
In the soils of the Lena River delta, there is an active accumulation of aliphatic fragments in the composition of Has of up to 70%. Within the profile dynamics, the processes of transformation of aliphatic compounds into aromatic ones were traced. Long-term formation ("maturation") of HAs led to an increase in the content of aromatic fragments in the buried organic matter. Evolutionary thermodynamic selection of condensed molecules led to the stabilization of organic matter and its deposition in the permafrost. It was noted that in soils formed under relatively homogeneous conditions (wind shelters), a relatively high content of aromatic fragments in the composition of HAs accumulated. Climatic features, microbiological activity and precursors of humification as well as time played important roles in the formation of stable ("mature") organic matter.
Author Contributions: V.P. performed the expedition with fieldwork and soil sampling; V.P. and E.A. wrote the paper. All authors have read and agreed to the published version of the manuscript.
Funding: This work was supported by a grant from the Russian Foundation for Basic Research, No. 19-05-50107.