Relating Paramagnetic Properties to Molecular Parameters of Humic Acids Isolated from Permafrost Peatlands in the European Arctic

Free radicals (FRs) are intermediate participants in the transformation process of soil organic matter, and free radical activity is a fundamental property of humic substances. The aim of this work was to conduct a comparative study of the paramagnetic properties of humic acids (HAs) isolated from Histosols by electron paramagnetic resonance (EPR) spectroscopy. The studied Histosols are found in permafrost peatlands in four natural geographic subzones of the European Arctic (from forest tundra to northern tundra). The results obtained showed that in anaerobic conditions on the peatlands in the tundra zone, the formation of semiquinone-type radicals occurs through the reduction of quinone fragments of HAs and leads to an increase in the concentration of paramagnetic centres within HAs. PCA analysis allowed us to reveal relationships between the properties of the initial raw peat samples, the molecular composition of the isolated HAs, and their paramagnetic parameters. It was found that FR localization occurs predominantly on aromatic fragments of lignin nature, which are confined to the low molecular weight fraction of HAs. The g-factor values of the EPR spectra of HAs indicate the presence of carbon- and oxygen-centred FRs in the HA structure, with a predominance of the latter.


Introduction
Many studies on soil organic matter confirm the important role of semiquinone-type organic free radicals (FRs) in biochemical processes [1,2].FRs are atoms or molecular groupings with an unpaired electron, which gives them paramagnetic properties and allows them to be detected by magnetic resonance spectroscopy.In inorganic systems, SRs are generally short-lived, but in complex organic compounds, they become stable, especially in the presence of aromatic nuclei.Semiquinone radicals are a classic example of such radicals.Modern concepts indicate that these radicals are intermediates in the formation of humic substances in soils [3,4], and free-radical activity is a fundamental property of soil organic matter (SOM) [5,6].
The use of electron paramagnetic resonance (EPR) to study organic matter in soils provides a variety of information about the structure of radical-containing substances [7,8].The paramagnetic properties of humic acids (HAs) depend on zonal bioclimatic factors and the level of anthropogenic impact [9][10][11].Previous studies have shown a positive correlation between the degree of humification and FR concentration [5].In addition, the content of paramagnetic centres in HA is determined by the degree of condensation and the colour coefficient [12].As a result of experiments, it was assumed that the formation of HAs occurs during the condensation of polyphenols, amino acid residues, etc., carried out through the stage of free-radical structures.The EPR method has high sensitivity and allows for obtaining original information on the structure of radical-containing substances.The study of the free-radical structure of biopolymers, including humic substances, makes it possible to estimate their ability to undergo polymerisation reactions [13].
Arctic wetlands contain significant carbon stocks stored in SOM [14].Permafrost degradation in this region leads to the release of melted SOM and increased carbon cycling [15].This is also facilitated by geomorphological processes such as thermokarst and erosion, which expose deep layers of organic matter [16].In addition, these areas are traditional reindeer grazing areas, which affects the stability and biodegradability of SOM as reindeer trample lichens and mosses, which are the main drivers of plant communities in the Arctic zone [17].
Studies using modern physicochemical methods to analyse the molecular structure of HAs in mineral soils have now been published [5,6,9].However, studies of HAs in European Arctic peatlands are rare and lack systematicity.Peatlands are specific natural soillike organic formations where the nature of humic compounds has been least studied [3].Analysis of study results showed that both natural factors and agricultural use of soils have a significant influence on the concentration of HAs in humic compounds [1,10,18,19].In this context, the aim of our work was to investigate the paramagnetic properties of HA samples isolated from peat hummocky bogs in four natural geographic subzones of the European Arctic and to reveal the relationship between paramagnetic properties and molecular parameters of HAs.

The EPR Data of HAs
The EPR spectra of the HA samples (Figure 1) showed a broad line with a g-factor close to that of a free electron (g = 2.0023), indicating the presence of a strongly delocalized molecular orbital in the HA structure.The g-factor of the EPR spectra is a useful parameter for identifying the nature of FRs.Radicals with g-factors in the range of 2.0020 to 2.0050 are organic radicals [20].Studies have shown that g-factors less than 2.0030 are carbon-centred radicals, whereas the g-factor of oxygen-centred radicals is greater than 2.0040 [21].
out through the stage of free-radical structures.The EPR method has high sensitivity and allows for obtaining original information on the structure of radical-containing substances.The study of the free-radical structure of biopolymers, including humic substances, makes it possible to estimate their ability to undergo polymerisation reactions [13].
Arctic wetlands contain significant carbon stocks stored in SOM [14].Permafrost degradation in this region leads to the release of melted SOM and increased carbon cycling [15].This is also facilitated by geomorphological processes such as thermokarst and erosion, which expose deep layers of organic matter [16].In addition, these areas are traditional reindeer grazing areas, which affects the stability and biodegradability of SOM as reindeer trample lichens and mosses, which are the main drivers of plant communities in the Arctic zone [17].
Studies using modern physicochemical methods to analyse the molecular structure of HAs in mineral soils have now been published [5,6,9].However, studies of HAs in European Arctic peatlands are rare and lack systematicity.Peatlands are specific natural soillike organic formations where the nature of humic compounds has been least studied [3].Analysis of study results showed that both natural factors and agricultural use of soils have a significant influence on the concentration of HAs in humic compounds [1,10,18,19].In this context, the aim of our work was to investigate the paramagnetic properties of HA samples isolated from peat hummocky bogs in four natural geographic subzones of the European Arctic and to reveal the relationship between paramagnetic properties and molecular parameters of HAs.

The EPR Data of HAs
The EPR spectra of the HA samples (Figure 1) showed a broad line with a g-factor close to that of a free electron (g = 2.0023), indicating the presence of a strongly delocalized molecular orbital in the HA structure.The g-factor of the EPR spectra is a useful parameter for identifying the nature of FRs.Radicals with g-factors in the range of 2.0020 to 2.0050 are organic radicals [20].Studies have shown that g-factors less than 2.0030 are carboncentred radicals, whereas the g-factor of oxygen-centred radicals is greater than 2.0040 [21].The EPR signal has the shape of a Lorentzian curve, which transitions to a Gaussian distribution on the wings.This indicates that, in addition to spin and spin-lattice interactions, the unpaired electron has sufficient freedom for spin-spin contacts [22,23].A characteristic feature of HAs is the redistribution of electron density in their molecular π-orbitals, which accounts for the paramagnetic properties of HAs [3,24].
Calculations of the integral intensity of the absorption line allowed us to estimate the concentrations of unpaired electrons in the samples (Table 1).The results showed that the concentration of paramagnetic centres in HA preparations is higher than in raw peat samples [6], due to the presence of a large number of aromatic and other polyconjugated structures in HA molecules, which may contain unpaired electrons. 1CPC-concentration of paramagnetic centres. 2 -the absolute error for determination of the g-factor is 0.00020.
The CPC in the investigated HA samples usually has a bimodal character, with maxima on the surface and permafrost horizons (Table 1).This is due to solar irradiation of the upper layer and a high degree of peat decomposition in the permafrost horizons, resulting in a higher content of aromatic and other polyconjugated structures in the peat HA composition, on which unpaired electrons can delocalize.Below 60 cm, there is an increase in the paramagnetic activity of HAs caused by the change in the botanical composition of the original peat and an increase in the proportion of wood residues and sedges enriched with syringyl-and guaiacyl-propane units of lignin fragments, which contribute to the stabilization of FRs and the formation of paramagnetic properties of HAs.
The g-factor values of the EPR spectra of the HAs from the peatlands studied vary between 2.00350 and 2.00372, indicating a mixture of carbon-and oxygen-centred radicals in the HA structure, with the latter predominating.The general trend of g-factor values observed for peatland HAs (Table 1) is not as pronounced, probably due to the different nature of the sediments.There is only a tendency for g-factor values to decrease with depth.The decrease in index values is explained by the delocalization of the unpaired electron from O,N-substituted structures (mainly semiquinones surrounded by O-containing fragments) to gradually increasing condensed aromatic fragments.Such a phenomenon has been observed in the natural carbonization processes of humic substances from compost, peats, and various types of coal [18,25,26].Experimental data for living peat-forming plants show relatively high values of g-factor: 2.0043 for Carex sp. and 2.0044 for Sphagnum sp.[27].Such high g-factor values have not been observed under the regenerative conditions of Arctic peatlands, even for the upper peat layers.

Principal Component Analysis
The paramagnetic properties of peat HAs are determined by a number of factors.Therefore, PCA was performed to determine the influence of the molecular composition of HAs isolated from peat samples with different chemical properties on the free radical concentration and g-factor of HAs.
The PCA results explained 63.21% of the total variability of the HA properties.The dimensionality of the 19 input variables was reduced by PCA to two principal components with eigenvalues greater than two: the first axis (PC1) explained 46.34% and the second axis (PC2) explained 16.87% of the total variability (Table 2).
The results indicate that the paramagnetic properties of HAs depend on both the molecular composition of the HAs themselves and the properties of the original raw peat samples (Figure 2).The results indicate that the paramagnetic properties of HAs depend on both molecular composition of the HAs themselves and the properties of the original raw p samples (Figure 2).  1 and 2.
The PCA results show that PC1 is positively related to CPC, low molecular we fraction (x(LMF)), C atom content of methoxyl (O-CH3) and aromatic fragments (O arom and C,H-arom), degree of peat decomposition (R), acidity (pH) and negativ related to g-factor, medium (x(MMF)) and high molecular weight (x(HMF)) fracti acidity (pH) and is negatively related to g-factor, medium (x(MMF)) and high molec weight   1 and 2.
The PCA results show that PC1 is positively related to CPC, low molecular weight fraction (x(LMF)), C atom content of methoxyl (O-CH 3 ) and aromatic fragments (O,N-arom and C,H-arom), degree of peat decomposition (R), acidity (pH) and negatively related to g-factor, medium (x(MMF)) and high molecular weight (x(HMF)) fractions, acidity (pH) and is negatively related to g-factor, medium (x(MMF)) and high molecular weight (x(HMF)) fractions, number average molecular weight (Mn), content of alkyl fragments (C,H-alkyl and O,N-alkyl), H/C element ratios and peat redox potential (Eh).
The proximity of the CPC, the content of C atoms of methoxyl and aromatic fragments in the HA structure, and the x(LMF) of HAs in the graph may indicate the localization of unpaired electrons predominantly on "lignin" fragments in the HA composition.This is also confirmed by the significant correlation between the concentration of paramagnetic centres and the content of O-CH 3 groups (r = 0.40, n = 43, r cr = 0.30), O,N-arom (r = 0.32, n = 43, r cr = 0.30) and C,H-arom fragments (r = 0.33) as well as the x(LMF) (r = 0.46, n = 43, r cr = 0.30).The findings are consistent with results indicating that the HMF of HAs consists mainly of C atoms of aliphatic structures, whereas the LMF of HA tends to contain aromatic fragments [28].Studies by Shi Y. et al. [29] also showed that a high proportion of oxygen-centred FRs is characteristic of low molecular weight HA samples (<3.5, <7, and <14 kDa).
A negative relationship between CPC and Eh values of native peat is observed, due to a shift in the equilibrium towards the formation of semiquinone radicals from quinone fragments as the redox potential decreases: Molecules 2023, 28, x FOR PEER REVIEW 6 of 14 that a high proportion of oxygen-centred FRs is characteristic of low molecular weight HA samples (<3.5, <7, and <14 kDa).
A negative relationship between CPC and Eh values of native peat is observed, due to a shift in the equilibrium towards the formation of semiquinone radicals from quinone fragments as the redox potential decreases:

Semiquinone
The peat layers of Arctic bogs are predominantly moderately and weakly reducing due to high humidity and the presence of anaerobic processes.Under these conditions, the humification of plant litter decreases, which increases the FR content, since in such regimes, paramagnetic centres of the semiquinone type are formed by the reduction of quinones.At the same time, semiquinone in the form of radical ions is most stable at higher pH values [5,30].This may be one of the reasons for the positive relationship between CPC and pH in peat samples (Figure 2).
The close relationship between the g-factor and the content of O,N-substituted alkyl fragments of HA is confirmed by a significant correlation coefficient between these parameters (r = 0.45, n = 43, rcr = 0.30).A statistically significant inverse relationship was found between the g-factor and the C-atom content of aromatic fragments.These results are in agreement with those obtained for HAs from the Baltic-type raised bog [31].The effect of a decrease in the g-factor of semiquinone radicals observed for HAs with increasing aromaticity [32] was confirmed and theoretically explained on the basis of quantum mechanical calculations of the g-factor for model semiquinone radicals [33].
On the second PC2 axis, a significant negative relationship was observed for the degree of peat decomposition, O/C ratio, and the content of carboxyl groups, and quinone fragments (Figure 2).This is due to the fact that an increase in the degree of decomposition of peat-forming plants, which are GC precursors, is accompanied by an increase in the proportion of O-containing structures.

Study Area
Four key wetland ecosystem polygons in the Arctic Belt were selected for this study (Figure 3), located in the Vorkuta District of the Komi Republic (KR) and the Nenets Autonomous Okrug (NAO).
These polygons represent key zonal biomes and are located in different areas of the Bolshezemelskaya tundra, including the forest tundra (Plot-1, Usa River basin, KR), the southern tundra (Plot-2, Khasyrei-Ty-Vis River basin, KR, and NAO border), the transitional ecotone of the northern and southern tundra (Plot-3, watershed of the Padimei-Ty-Vis and Korotaikha rivers, NAO), and the northern tundra (Plot-4, eastern part of the Pechora Bay of the Barents Sea, NAO).The studies were conducted in an area with massive islands and continuous permafrost [34].The peat layers of Arctic bogs are predominantly moderately and weakly reducing due to high humidity and the presence of anaerobic processes.Under these conditions, the humification of plant litter decreases, which increases the FR content, since in such regimes, paramagnetic centres of the semiquinone type are formed by the reduction of quinones.At the same time, semiquinone in the form of radical ions is most stable at higher pH values [5,30].This may be one of the reasons for the positive relationship between CPC and pH in peat samples (Figure 2).
The close relationship between the g-factor and the content of O,N-substituted alkyl fragments of HA is confirmed by a significant correlation coefficient between these parameters (r = 0.45, n = 43, rcr = 0.30).A statistically significant inverse relationship was found between the g-factor and the C-atom content of aromatic fragments.These results are in agreement with those obtained for HAs from the Baltic-type raised bog [31].The effect of a decrease in the g-factor of semiquinone radicals observed for HAs with increasing aromaticity [32] was confirmed and theoretically explained on the basis of quantum mechanical calculations of the g-factor for model semiquinone radicals [33].
On the second PC2 axis, a significant negative relationship was observed for the degree of peat decomposition, O/C ratio, and the content of carboxyl groups, and quinone fragments (Figure 2).This is due to the fact that an increase in the degree of decomposition of peat-forming plants, which are GC precursors, is accompanied by an increase in the proportion of O-containing structures.

Study Area
Four key wetland ecosystem polygons in the Arctic Belt were selected for this study (Figure 3), located in the Vorkuta District of the Komi Republic (KR) and the Nenets Autonomous Okrug (NAO).
These polygons represent key zonal biomes and are located in different areas of the Bolshezemelskaya tundra, including the forest tundra (Plot-1, Usa River basin, KR), the southern tundra (Plot-2, Khasyrei-Ty-Vis River basin, KR, and NAO border), the transitional ecotone of the northern and southern tundra (Plot-3, watershed of the Padimei-Ty-Vis and Korotaikha rivers, NAO), and the northern tundra (Plot-4, eastern part of the Pechora Bay of the Barents Sea, NAO).The studies were conducted in an area with massive islands and continuous permafrost [34].
The peat layers of Arctic bogs are predominantly moderately and weakly reducing due to high humidity and the presence of anaerobic processes.Under these conditions, the humification of plant litter decreases, which increases the FR content, since in such regimes, paramagnetic centres of the semiquinone type are formed by the reduction of quinones.At the same time, semiquinone in the form of radical ions is most stable at higher pH values [5,30].This may be one of the reasons for the positive relationship between CPC and pH in peat samples (Figure 2).
The close relationship between the g-factor and the content of O,N-substituted alkyl fragments of HA is confirmed by a significant correlation coefficient between these parameters (r = 0.45, n = 43, r cr = 0.30).A statistically significant inverse relationship was found between the g-factor and the C-atom content of aromatic fragments.These results are in agreement with those obtained for HAs from the Baltic-type raised bog [31].The effect of a decrease in the g-factor of semiquinone radicals observed for HAs with increasing aromaticity [32] was confirmed and theoretically explained on the basis of quantum mechanical calculations of the g-factor for model semiquinone radicals [33].
On the second PC2 axis, a significant negative relationship was observed for the degree of peat decomposition, O/C ratio, and the content of carboxyl groups, and quinone fragments (Figure 2).This is due to the fact that an increase in the degree of decomposition of peat-forming plants, which are GC precursors, is accompanied by an increase in the proportion of O-containing structures.

Study Area
Four key wetland ecosystem polygons in the Arctic Belt were selected for this study (Figure 3), located in the Vorkuta District of the Komi Republic (KR) and the Nenets Autonomous Okrug (NAO).
These polygons represent key zonal biomes and are located in different areas of the Bolshezemelskaya tundra, including the forest tundra (Plot-1, Usa River basin, KR), the southern tundra (Plot-2, Khasyrei-Ty-Vis River basin, KR, and NAO border), the transitional ecotone of the northern and southern tundra (Plot-3, watershed of the Padimei-Ty-Vis and Korotaikha rivers, NAO), and the northern tundra (Plot-4, eastern part of the Pechora Bay of the Barents Sea, NAO).The studies were conducted in an area with massive islands and continuous permafrost [34].According to climatic zoning, the study area belongs to the eastern subarctic subzone [35].The climate here is dominated by Arctic air masses, with continentality increasing inland.The annual air temperature amplitude reaches 31-35 °C.The average annual air temperature in this area varies from minus 4 to minus 9 °C.The climatic conditions in this area are characterized by harsh conditions.Winter lasts about 250 days, and the frost period (when the average daily temperature is below minus 20 °C) lasts 80-95 days.The snow cover lasts from 230 to 240 days.The period of positive temperatures lasts from 100 to 115 days.
The moisture regime in the area is characterized by a high precipitation rate, moderate evaporation, and high air and soil moisture.Annual rainfall decreases from south to north and ranges from 500 to 360 mm.Maximum rainfall occurs during the warm season, usually in August, and accounts for about 70 percent of the total annual rainfall (200 to 300 mm).Depending on the type of vegetation cover, the depth of the roof of perennially frozen rocks in loamy soils varies and can range from 0.7-0.8 to 2.0-3.0 m or more [36].Moss-shrub and coarse-shrub (bird and willow) tundras are widespread on plateaus.In areas with poor drainage, sedges and flat-bumped bogs are found.In the Arctic, vegetation cover is affected by the recurring seasonal processes of thawing and freezing of soil water, fluidity of over-watered clay soils, soil frost heaving, crack formation, snow corrasion of the soil surface, turf breaks, and spotting [37,38].Reindeer grazing and caterpillar transport have a major impact on vegetation cover [17].

Forest Tundra
In the forest tundra, the studies were conducted within the flat-bumpy-mossy According to climatic zoning, the study area belongs to the eastern subarctic subzone [35].The climate here is dominated by Arctic air masses, with continentality increasing inland.The annual air temperature amplitude reaches 31-35 • C. The average annual air temperature in this area varies from minus 4 to minus 9 • C. The climatic conditions in this area are characterized by harsh conditions.Winter lasts about 250 days, and the frost period (when the average daily temperature is below minus 20 • C) lasts 80-95 days.The snow cover lasts from 230 to 240 days.The period of positive temperatures lasts from 100 to 115 days.
The moisture regime in the area is characterized by a high precipitation rate, moderate evaporation, and high air and soil moisture.Annual rainfall decreases from south to north and ranges from 500 to 360 mm.Maximum rainfall occurs during the warm season, usually in August, and accounts for about 70 percent of the total annual rainfall (200 to 300 mm).Depending on the type of vegetation cover, the depth of the roof of perennially frozen rocks in loamy soils varies and can range from 0.7-0.8 to 2.0-3.0 m or more [36].Moss-shrub and coarse-shrub (bird and willow) tundras are widespread on plateaus.In areas with poor drainage, sedges and flat-bumped bogs are found.In the Arctic, vegetation cover is affected by the recurring seasonal processes of thawing and freezing of soil water, fluidity of over-watered clay soils, soil frost heaving, crack formation, snow corrasion of the soil surface, turf breaks, and spotting [37,38].Reindeer grazing and caterpillar transport have a major impact on vegetation cover [17].

Forest Tundra
In the forest tundra, the studies were conducted within the flat-bumpy-mossy complex on Hemic Folic Cryic Histosol-Plot-1.In the following, the names of the soils are given in brackets according to the international WRB classification [39].
The bog massif under consideration occupies a high terrace above the floodplain in the left-bank part of the Sedyakha River valley (a right-bank tributary of the Usa River).The absolute height of the terrace is about 90 m a.s.l., and the height difference at the watershed of the Sedyakha and Seyda Rivers is 100-160 m a.s.l.The bog massif is a boggy marsh massif, which is a boggy marshland.The peat bog is a ridge-hollow complex with thermokarst lakes.Peat hummocks of various shapes prevail within the bog boundaries; the height of hummocks is from 1 to 3 m.The studied bog is characterized by a well-defined microrelief, largely caused by permafrost processes.Approximately 60% of the massif is occupied by rises or hills, the rest being occupied by hollows, sinks, watered runoff hollows, and secondary lakes.The hills are rounded or oval in shape and their cross-sectional dimensions vary from 10-15 to 20-45 m.The surface is flat or finely punctuated.As a result of permafrost degradation, saucer-shaped depressions with hygrophyte vegetation are formed on some hillocks.The thickness of the seasonally thawed layer varies from 40 to 50 cm.The peat is dark brown throughout its thickness, with a high and medium degree of decomposition, and is of the humic type.
Data on the composition and physico-chemical properties of peat soils and HA preparations isolated from them are presented in Table 3.A more detailed characterization and description of the molecular properties of HAs have been published earlier [38,40].

South Tundra
The study area is located in the southern tundra on the eastern macro slope of the Tarju and Vorkuta River watersheds, with a few streams and small rivers belonging to the Vorkuta River basin.
Plot-2: Fibric Folic Cryic Histosol, 67 • 45 N, 63 • 18 E.The plot is located on a flat peat plateau in the NAO, 28 km NW of Vorgashor, 3.5 km west of Lake Lek'yamboty, in the basin of the Lek'yambo-Ty-Vis River.The surface of the plateau has a tuberous microrelief; Ledum decumbens, Rubus chamaemorus dominate the vegetation cover, and the moss cover consists of Polytrichum commune, Pleurosium schreberi, species g.Dicranum, and g.Cladonia.
The study area is located in the Bolshezemelskaya tundra in the north-eastern part of the Pechora Lowland.Low hills and ridges, apparently of glacial and glacio-marine origin, dominate most of the area, alternating with wide gullies and lake basins.The hydrogeophysical network in the study area is rather dense, mainly due to small lakes.These are mainly small water bodies with a small surface area.The lakes are mainly located in the inter-block joints and therefore have pronounced depressions and mineralfree shores.Studies in the ecotone zone of the northern tundra and the southern tundra were carried out on the soil of the Hemic Folic Cryic Histosol-Plot-3.
Plot-3: 68 • 02 N, 62 • 43 E; the plot is located on a flat hillock at the edge of a flat bog with a lake complex.The hill has a complex in shape, 15-20 m wide, 20-25 m long, and 0.7-1.5 m high.Microprofile of the hill: diameter 30-80 cm, height 10-15 cm.Interhilly spaces are irregularly shaped, up to 30-50 cm in diameter.The shrub layer consists of Ledum decumbens, Betula nana, Empetrum nigrum, Vaccinium vitis-idaea, and Rubus chamaemorus.Shrubs mainly occupy the surface of the hillocks and the peripheral part of the peat hillock, including its sloping surface.The interhilly spaces are mainly occupied by moss-lichen vegetation with dicranial and polytrichous mosses and various lichen species.The depth of the upper limit of the permafrost is 33-44 cm, depending on the micro-relief.

North Tundra
In the northern tundra, studies were carried out on the Barents Sea coast within the flat-topped peat mound-hollow complex: Fibric Folic Cryic Histosol-Plot-4.
Plot-4: 68 • 35 N, 55 • 55 E; the plot is located in the centre of a ridge of the flat-bumpy bog.The shape of the hill is complex, polygonal, elongated from north to south, with sharp depressions of hollow and lake depressions.The hill has frost cracks crossing the ridges and bare areas without vegetation cover, d = 0.3-0.7 m.The bog is characterised by a well-defined micro-relief, largely caused by permafrost processes.The relief of the hills is shallow-to medium-hilly.The diameter of the hills varies from 50 to 100 cm.The shape of the hills varies from elongated oval to polygonal, bounded on the sides by lacustrine-hollow depressions.The diameter of the ridge at the sampling site is about 8 m, with a height of 0.5-0.7 m.Vegetation is represented by Ledum decumbens, Rubus chamaemorus, Andromeda, Empetrum nigrum, and Vaccinium vitis-idaea, as well as lichens and bryophyte mosses.The depth of the permafrost is 29-30 cm.

Methods
The method recommended by the International Humic Substances Society [41] was used for the extraction of HAs.HAs were extracted from 50 g of dry peat samples by double extraction with a 0.1 M NaOH solution at a ratio of 1:10 for complete extraction of HAs.The mixture was centrifuged for 1 h at 10,000 rpm on a SIGMA 2-16 KL ultracentrifuge (Sigma Laborzentrifugen GmbH, Osterode am Harz, Germany).The HA was then precipitated with a 6 M HCl solution to pH 1.0.HA was further purified from ash constituents using a mixture of 0.1 M HCl and 0.3 M HF.The ash content of HA preparations does not exceed 1%.To remove low molecular weight compounds, the HAs were dialyzed, then dried at 35 • C in a laboratory oven with forced convection, ground, and sieved through a 0.1 mm sieve.
Registration of EPR spectra of HA samples was carried out at the Research Centre "Diagnostics of structure and properties of nanomaterials"-Kuban State University on the spectrometer JES FA 300 (JEOL, Tokyo, Japan) in the X-band.Measurement conditions: microwave power 1 mW, amplitude of high-frequency modulation 0.06 mTl.Manganese (II) oxide, with a known radical content, was used as an external standard to determine the concentration of paramagnetic centres (CPC) and as a reference for calculating the g-factors of the investigated HA samples.CPC in HA samples was determined by comparing the relative signal intensities of the sample and the standard using JES-FA swESR software v. 3.0.0.1 (JEOL, Tokyo, Japan).The absolute error in the determination of CPC by the EPR method is up to 10% [9].Bivariate correlation analyses were performed using Pearson's product-moment correlation coefficient (r), and statistical significance was assessed using the Neyman-Pearson (normal distribution) approach.The technically observed value of the coefficient (based on n pairs) was compared with the critical value (r cr ) for a two-tailed test and a significance level of 0.05.
Principal component analysis (PCA) was performed to determine correlations between paramagnetic and molecular parameters of HAs using Statistica v. 12.1 (Dell, Round Rock, TX, USA).The number of factors extracted from the variables was determined by Kaiser's rule.The first two principal components with an eigenvalue greater than two were retained for this criterion.All statistical analyses were performed at the specified significance level of p ≤ 0.05.

Conclusions
In the present work, the paramagnetic properties of HA preparations isolated from Histosols of four polygons of Arctic Belt peat ecosystems were investigated.These polygons are key zonal biomes and are located in different areas of the Bolshezemelskaya tundra, including forest tundra, southern tundra, northern and southern tundra ecotone, and northern tundra.
The results showed that the redox conditions of peatlands as well as the molecular parameters of HAs determine their paramagnetic properties.The peat layers of Arctic peat bogs have predominantly moderate and weak redox conditions due to high humidity and the occurrence of anaerobic processes.Under these conditions, FRs of the semiquinoid type are formed by the reduction of quinone fragments of HAs, which leads to an increase in the content of paramagnetic centres in the composition of HAs.At the same time, their concentration in HAs usually has a bimodal character of distribution along the profile of the studied peatlands, with maxima on the surface and permafrost horizons.
PCA of the obtained data set showed that the localization of unpaired electrons occurs predominantly on "lignin" fragments, which are restricted to the low molecular weight fraction of peat HAs.The g-factor values of the EPR spectra of HAs indicate the presence of carbon-and oxygen-centred radicals in the HA structure, with a predominance of the latter.

Figure 2 .
Figure 2. Projection of HA parameters using PCA: abbreviations as in Tables1 and 2.
(x(HMF)) fractions, number average molecular weight (Mn), content of a fragments (C,H-alkyl and O,N-alkyl), H/C element ratios and peat redox potential (Eh The proximity of the CPC, the content of C atoms of methoxyl and arom fragments in the HA structure, and the x(LMF) of HAs in the graph may indicate localization of unpaired electrons predominantly on "lignin" fragments in the composition.This is also confirmed by the significant correlation between concentration of paramagnetic centres and the content of O-CH3 groups (r = 0.40, n = rcr = 0.30), O,N-arom (r = 0.32, n = 43, rcr = 0.30) and C,H-arom fragments (r = 0.33) as w as the x(LMF) (r = 0.46, n = 43, rcr = 0.30).The findings are consistent with results indica that the HMF of HAs consists mainly of C atoms of aliphatic structures, whereas the L of HA tends to contain aromatic fragments [28].Studies by Shi Y. et al. [29] also show

Figure 2 .
Figure 2. Projection of HA parameters using PCA: abbreviations as in Tables1 and 2.

←
Molecules 2023, 28, x FOR PEER REVIEW 6 of 14 that a high proportion of oxygen-centred FRs is characteristic of low molecular weight HA samples (<3.5, <7, and <14 kDa).A negative relationship between CPC and Eh values of native peat is observed, due to a shift in the equilibrium towards the formation of semiquinone radicals from quinone fragments as the redox potential decreases: Quinone ⎯⎯⎯ ⎯⎯⎯ Semiquinone

Table 1 .
Paramagnetic properties of HAs from Arctic Histosols.

Table 2 .
Analysis of principal components for HA parameters.

Table 3 .
Properties of raw peats and HAs.