Differentiation of Trace Metal Contamination Level between Different Urban Functional Zones in Permafrost Affected Soils (the Example of Several Cities in the Yamal Region, Russian Arctic)

: Dynamically developing urbanization causes a number of environmental effects, including those related to the chemical transformation of soils. Relatively less information about the urban areas of the Arctic and Subarctic zones, constructed mostly on permafrost and intensively populated areas can be found. By the example of the analysis of basic soil properties and concentrations of trace metals in the soils of the cities of Salekhard, Urengoy, Nadym, Novy Urengoy and Gaz Sale (the Yamalo-Nenets Autonomous District), as well as various functional zones within the cities, the relationship between the age of the cities, the level of anthropogenic pressure and the type of parent materials and the character of accumulation of metals in the soil proﬁle of urban soils have been described. The direct correlation was found between the content of Pb, Cr, Ni, As and soil sorption characteristics. In young cities built on sandy sediments, there is less accumulation of heavy metals in the topsoil horizons. Relatively higher concentrations of Cu and Cd were noted in soils of industrialized cities, regardless of functional zones. The higher content of Cr, Ni, Cu, Zn, As and Pb has been registered in older zones also frequently used for residential purposes. The calculated values of the PI index for some functional zones of young cities show the medium and high content of heavy metals. The analysis of Igeo and PLI indices shows a large diversity both in relation to individual cities and their functional zones. Soil quality, in spite of the high level of anthropogenic load, was assessed as mostly satisfactory.


Introduction
Urbanization and urban sprawl has to a large extent a demographic backgroundmodern cities accept a constantly growing number of inhabitants-constituting already 55.3%, and in 2050, possibly 68.4% of the total world population [1]. Strong urbanization processes will be noted in previously sparsely populated and urbanized areas. It should be expected that they will proceed very intensively in areas of key importance to the economies of individual countries. Such areas are strongly related to the presence of raw materials and dynamically developing industry.
Climatic conditions, the isolation of settlements from each other, and transport inaccessibility cause a high degree of urbanization of the Arctic regions. According to various estimates, more than 85% of the population of the Russian Arctic lives in cities. This is the highest rate of urbanization in all of Russia [2]. Large settlements in the Arctic regions were formed as a result of the industrial development of the north. The effective mining of natural resources was impossible without the logistics and transport chain, the basis of which are large settlements. In the Yamal-Nenets Autonomous District, the largest settlements are inextricably linked to industrial production [3]. The first attempts to apply scientifically based methods of complex development and management of the Arctic territories began during the period of mass industrial activity in the oil and gas fields of Western Siberia. As a result, it was considered optimal to create cities of different categories: large, medium or small, depending on the specific needs of the region. As a result, in the period from 1960 to 1990 there were created 6 cities and 14 small settlements, which were focused on the search and extraction of hydrocarbon raw materials. At present about 75% of the population of the Yamal region lives in these towns [4,5].
Most of the Arctic cities were formed relatively recently. The oldest city in which our study was conducted is Salekhard, it has an ancient history and was founded in 1595 and received city status in 1938. The ages of the rest of the cities we studied do not exceed 45 years, for example, Gaz Sale was formed as a municipality only in 2004 [4].
The understanding of the conceptual role of the Arctic territories and the elaboration of projects for the development of the North as a special space for state development are widely associated with the works of academicians Agranat [6] and Slavin [7]. Currently, the development of the Arctic regions is being raised more and more frequently. The Russian Federation has adopted at the governmental level the program "Social-economic development of the Arctic zone of the Russian Federation for the period until 2025," which provides for the development of the northern regions in various directions, from mining to agriculture [8].
It is reasonable to assume that such global projects will necessarily lead to population growth and density in the Arctic regions and, as a consequence, to the increase and development of cities in the Arctic and Subarctic zones. However, intensive urban growth carries with it adverse environmental effects on the fragile Arctic ecosystem.
The environmental problems of cities and their impact on the global ecosystem are well studied and unquestioned [9,10]. Urban development, and the resulting point-in-time extreme anthropogenic impact on ecosystems, has been well studied in the examples of the largest cities and their agglomerations in Europe, China, Latin America and India [11][12][13][14]. The development of urban infrastructure cardinally transforms landscapes and changes the temperature regime due to the concentration of industrial enterprises and vehicles. This leads to the formation of so-called "heat islands" which are one of the causes of climate change at the global and local levels [9,10].
Urban areas located in the Arctic and Subarctic zones require special management in terms of their impact on the environment. Built mostly on permafrost, large populated settlements may become one of the causes of the acceleration of permafrost melting, which has enormous impacts. One of the effects of thawing permafrost is a massive increase in the emission of greenhouse gases such as CO 2 or methane into the atmosphere [15,16]. Since more than 60% of Russia's territory is located in the permafrost zone, questions of controlling the anthropogenic load on Arctic territories will always remain relevant [17,18].
Dynamically developing urbanization causes a number of environmental effects, including those related to the transformation of soils. The pace of land development is proportional to the scale of these transformations [19,20]. The soil cover of the northern regions is highly vulnerable to anthropogenic impact. Soil formation in the Arctic is complicated due to extreme climatic conditions and the prevalence of low temperatures. However, it is the analysis of the soil cover that may be the key to assessing the degree of anthropogenic impact on Arctic ecosystems. Additionally, all areas with poorly developed soil cover (poorly developed soils) are susceptible to their devastation as a result of investment processes [21]. The mechanical transformations of these soils lead to the shortening of the soil profiles and the exposure of the bedrock level [22]. This enables other manifestations of soil degradation, including chemical ones, with a direct threat to the lower parts of the soil profile and geological layers [22,23]. Accumulation of trace metals in urban soils is an unavoidable effect of excessive anthropogenic pressure on ecosystems [24,25].
Since the adverse environmental impact of the urban ecosystem of the Yamal region is associated not only with traditional factors such as emissions from industrial plants and motor vehicles, but also active extraction of natural resources (gas, oil, chromite ores), it is extremely important to understand the degree of anthropogenic load which falls not only on natural ecosystems, but also on urbanized areas [26][27][28]. Thus, the goal of our work is to assess the content of trace metals in the soils of Arctic cities and their functional zones (Salekhard, Nadym, Urengoy, Gaz Sale and Novy Urengoy) with different histories of anthropogenic load and different soil-forming bedrock.

Site Locations and Their Anthropogenic Load
The study focused on topsoil horizons for differentiation of basic soil properties and content of trace elements between several large settlements in the Arctic region ( Figure 1). A comparison was also made between the different functional zones of the northern cities. Soil samples were taken in five major towns of the Yamal-Nenets Autonomous District, namely: -Salekhard-settlement foundation~1595, town status 1938, total area 84.5 km 2 , population 51 thous. (2020), capital city of the Yamalo-Nenets Autonomous Okrug, chromite deposits, oil and gas industry; -Nadym-settlement foundation~1597, town status 1972, total area 99.8 km 2 , population 45 thous. (2020), oil and gas industry; -Novy Urengoy-construction of a research station 1966, town status 1980, total area 111 km 2 , population 118 thous. (2020), oil and gas industry; -Gaz Sale-settlement, from 2025 will be included in Tazovskiy municipal area, total area 0.03 km 2 , population 1721 (2019), oil and gas industry; -Urengoy-construction of the geological station 1966, town status 1979, total area 2.8 km 2 , population 10 thous. (2020), oil and gas industry.
Anthropogenic load of cities in the Yamal region is associated not only with traditional factors for urbanized areas, such as emissions from industrial plants and vehicles, burning of fossil fuels, degradation of green spaces and many others. The world's largest chromite deposit is located near the towns of Labytnangi and Salekhard. It is mined using the open-pit method and the ore is transported to processing plants [26,27]. Additionally on the Yamal Peninsula is one of Russia's largest natural gas and oil production and transportation facilities. Of course, all of these processes lead to the inflow of trace metals into the environment, surface waters, and soils [28,29]. Young and small cities such as Gaz Sale and Nadym are vulnerable to emissions from local sources of emissions (e.g., combined heat and power plants) and vehicles [30]. Minerals 2021, 11, x FOR PEER REVIEW 4 of 21

Sampling Strategy
The sampling strategy was based on searching for different functional zones within several cities of the Yamal region and collecting samples from topsoil (0-20) horizons. Soil samples were taken in the industrial, residential, and park zones of Salekhard, Novy Uregnoi, Gaz Sale, Nadym, and Urengoi. Background (mature stand) samples were also collected from outside the territories of the above-mentioned cities ( Table 1). Within each functional zone, samples were taken from individual sites typical for a given form of use, in three replicates. A total of 60 soil samples were taken.
In the studied areas, soil pits were also made, exposing the soil profile to the permafrost level or to the subcutaneous waters. The described limiting elements of the soil profiles from below were located not deeper than 100 cm below the ground level ( Figure 2). According to the International Soil Classification System [31], the studied soils belong to the taxonomic sections of Podzols, Technosols, and Cryosols.

Sampling Strategy
The sampling strategy was based on searching for different functional zones within several cities of the Yamal region and collecting samples from topsoil (0-20) horizons. Soil samples were taken in the industrial, residential, and park zones of Salekhard, Novy Uregnoi, Gaz Sale, Nadym, and Urengoi. Background (mature stand) samples were also collected from outside the territories of the above-mentioned cities ( Table 1). Within each functional zone, samples were taken from individual sites typical for a given form of use, in three replicates. A total of 60 soil samples were taken.
In the studied areas, soil pits were also made, exposing the soil profile to the permafrost level or to the subcutaneous waters. The described limiting elements of the soil profiles from below were located not deeper than 100 cm below the ground level ( Figure 2). According to the International Soil Classification System [31], the studied soils belong to the taxonomic sections of Podzols, Technosols, and Cryosols.  . Soil diversity in the study cities and functional zones (a1-a5-mature; b1-b5-residential; c1-c5-park; d1-d5industrial).

Chemical Analyses and Data Processing
Soil samples were transported to the laboratory of the Department of Applied Ecology, St. Petersburg University. The samples were air dried, grounded and passed through a 2 mm sieve. The pH values of soil solution were measured by the using a pH-metermillivoltmeter pH-150MA (Belarus). Soil solution was prepared in the ratio of 1:2.5 with water or 1M CaCl2 (for mineral soils the optimal soil weight for solution preparation is 8 g) [32]. Basal soil respiration was evaluated by measuring CO2 in Sodium Hydroxide. Incubation of CO2 was conducted for 10 days in plastic sealed containers [33].
The total carbon and nitrogen contents were evaluated using a CHN analyzer (Leco CHN-628; Leco Corporation, St. Joseph, MI).

Chemical Analyses and Data Processing
Soil samples were transported to the laboratory of the Department of Applied Ecology, St. Petersburg University. The samples were air dried, grounded and passed through a 2 mm sieve. The pH values of soil solution were measured by the using a pH-metermillivoltmeter pH-150MA (Belarus). Soil solution was prepared in the ratio of 1:2.5 with water or 1M CaCl 2 (for mineral soils the optimal soil weight for solution preparation is 8 g) [32]. Basal soil respiration was evaluated by measuring CO 2 in Sodium Hydroxide. Incubation of CO 2 was conducted for 10 days in plastic sealed containers [33].
The total carbon and nitrogen contents were evaluated using a CHN analyzer (Leco CHN-628; Leco Corporation, St. Joseph, MI, USA).
The Geoaccumulation Index (Igeo) allows us to classify seven levels of soil contamination, from Practically unpolluted (Igeo ≤ 0) to Extremely polluted (Igeo > 5) [35,36]. The generic calculation formula is as follows: where, C n -the measured concentration of the element in soil; B n -the geochemical background value. One more index that assists in evaluating what heavy metal contributes the most to soil contamination is the Single Pollution Index (PI) [36]. The basic formula is as follows: where, C n -the content of heavy metal in soil; GB-values of the geochemical background. Single Pollution Index is used to calculate complex indices, which allow you to estimate the total degree of pollution. One of such indices is Pollution Load Index (PLI), it is calculated as a geometric average of PI values [36,37]. The basic formula for calculating this complex index is as follows: where, n-the number of analyzed metals; PI-calculated values for the Single Pollution Index.
All applied indices have their own evaluation scales (Table 2), which allows a qualitative assessment of soil contamination by trace metals.

Variance Analysis and Data Visualization
Statistical analysis of the data was performed using multivariate analysis of variance (MANOVA) in the Statistica v12.0 software package (Statsoft, Tulsa, OK, USA). Data visualization was performed using QGIS 3.16 and GraphPad Prizm 9.0.0.

Bedrock Material
The investigated urban areas are located on modern alluvial floodplain and lakealluvial sediments. The city of Salekhard and its vicinity are characterized by bedrock of predominantly clay and loamy structure. Urengoy, Nadym, Novy Urengoy and Gaz Sale are located on Quaternary sediments of predominantly sandy composition [40,41] ( Figure 2). The particle size distribution of soils in the Yamal region is highly differentiated. In the vicinity of Salekhard, the clay content reaches 28% [42]. In the pristine soils of Belyi Island the content of the clay fraction does not reach 5%, the largest percentage is the sand fraction, its content can reach 78% [43]. In the topsoil horizons in the cities of Gaz-Sale and Urengoy, the content of clay particles does not exceed 12%, the largest proportion in the fine soil is taken by sand particles, and their content can reach 87%.
The bedrock of the soil, located close to the ground surface, largely determines the possibility of implementing municipal and industrial investments [21]. It also influences the behavior of chemical substances from atmospheric emissions [22,44]. Compact clay material with the ability to swell and shrink is defective from the point of view of construction activities as some green areas as well [45]. This favors its removal from the investment area or mixing with skeletal materials, e.g., gravel and construction rubble. Sandy soils are usually treated differently, being easy to build. Nevertheless, they are characterized by high water permeability. As a result, they show a much deeper deposition of pollutants in the soil profile compared to the finer materials [46]. Regardless of texture, soils are affected to varying degrees by the deposition of anthropogenic materials, including construction and municipal waste [47][48][49].

Main Soil Properties
The analysis of the reaction of the soil solution showed that all of the studied soils are characterized as acidic (pH < 7). There are strong differences between the cities and between the functional zones within the cities (Figure 3). predominantly clay and loamy structure. Urengoy, Nadym, Novy Urengoy and Gaz Sale are located on Quaternary sediments of predominantly sandy composition [40,41] ( Figure  2). The particle size distribution of soils in the Yamal region is highly differentiated. In the vicinity of Salekhard, the clay content reaches 28% [42]. In the pristine soils of Belyi Island the content of the clay fraction does not reach 5%, the largest percentage is the sand fraction, its content can reach 78% [43]. In the topsoil horizons in the cities of Gaz-Sale and Urengoy, the content of clay particles does not exceed 12%, the largest proportion in the fine soil is taken by sand particles, and their content can reach 87%. The bedrock of the soil, located close to the ground surface, largely determines the possibility of implementing municipal and industrial investments [21]. It also influences the behavior of chemical substances from atmospheric emissions [22,44]. Compact clay material with the ability to swell and shrink is defective from the point of view of construction activities as some green areas as well [45]. This favors its removal from the investment area or mixing with skeletal materials, e.g., gravel and construction rubble. Sandy soils are usually treated differently, being easy to build. Nevertheless, they are characterized by high water permeability. As a result, they show a much deeper deposition of pollutants in the soil profile compared to the finer materials [46]. Regardless of texture, soils are affected to varying degrees by the deposition of anthropogenic materials, including construction and municipal waste [47][48][49].

Main Soil Properties
The analysis of the reaction of the soil solution showed that all of the studied soils are characterized as acidic (pH < 7). There are strong differences between the cities and between the functional zones within the cities (Figure 3).  The soils of the Yamal region are predominantly acidic (pH < 7) [50][51][52][53][54]. For natural Al-Fe humus soils, high acidity is determined by the processes of formation and transposition of humus, aluminum and iron compounds [55]. According to the data given above (Figure 3), it can be seen that in all the studied cities the lowest values of pH are in the functional zone with the lowest anthropogenic load (mature) pH values mainly pH H 2 O < 5, pH CaCl 2 < 4. Within cities, pH values are somewhat higher, especially in residential and industrial areas. Particularly highlighted values of pH in the anthropogenic functional zones of the city of Gaz Sale, here pH is characterized as slightly acidic, may indicate the presence of anthropogenic alkalization of the soil, which is typical for cities with strong industrial development and high traffic load [56,57].
The level of soil basal respiration is highly variable between cities and their functional zones (Figure 4). The soils of the Yamal region are predominantly acidic (pH < 7) [50][51][52][53][54]. For natural Al-Fe humus soils, high acidity is determined by the processes of formation and transposition of humus, aluminum and iron compounds [55]. According to the data given above (Figure 3), it can be seen that in all the studied cities the lowest values of pH are in the functional zone with the lowest anthropogenic load (mature) pH values mainly pH H2O < 5, pH CaCl2 < 4. Within cities, pH values are somewhat higher, especially in residential and industrial areas. Particularly highlighted values of pH in the anthropogenic functional zones of the city of Gaz Sale, here pH is characterized as slightly acidic, may indicate the presence of anthropogenic alkalization of the soil, which is typical for cities with strong industrial development and high traffic load [56,57].
The level of soil basal respiration is highly variable between cities and their functional zones (Figure 4). In most cases, the maximum level of biological activity is recorded in park and mature areas. Since outside the city and in its green areas the processes of accumulation of soil organic matter occur most actively, it causes a high level of carbon dioxide emissions by soil microorganisms. Nevertheless, the maximum level of carbon dioxide emissions was recorded in the residential area of Salekhard. The minimum values are recorded in the city of Gaz Sale for all its functional areas.
To understand the degree of degradation of soil organic matter and the rate of its mineralization, it is essential to know the ratio of total organic carbon to total nitrogen (C/N). With a permanent source of organic matter, high values indicate good preservation of organic material [58].
At a high rate of mineralization of organic matter for the zone of permafrost spreading, the C/N ratio is usually less than 12, the moderate rate of mineralization corresponds to values of 13-25, at higher values the degree of mineralization is characterized as low [59,60]. For urban and anthropogenically transformed soils, where the processes of fresh In most cases, the maximum level of biological activity is recorded in park and mature areas. Since outside the city and in its green areas the processes of accumulation of soil organic matter occur most actively, it causes a high level of carbon dioxide emissions by soil microorganisms. Nevertheless, the maximum level of carbon dioxide emissions was recorded in the residential area of Salekhard. The minimum values are recorded in the city of Gaz Sale for all its functional areas.
To understand the degree of degradation of soil organic matter and the rate of its mineralization, it is essential to know the ratio of total organic carbon to total nitrogen (C/N). With a permanent source of organic matter, high values indicate good preservation of organic material [58].
At a high rate of mineralization of organic matter for the zone of permafrost spreading, the C/N ratio is usually less than 12, the moderate rate of mineralization corresponds to values of 13-25, at higher values the degree of mineralization is characterized as low [59,60]. For urban and anthropogenically transformed soils, where the processes of fresh organic matter input are highly limited and the soil-plant chain is disrupted, the C/N ratio can fluctuate greatly [61]. Additionally, high values of the C/N ratio may be associated with the presence of aromatic hydrocarbons and petroleum products in anthropogenically disturbed soils [62].
According to our results ( Figure 5), we can see that the level of carbon, nitrogen and their ratios are highly variable. For the city of Salekhard, as the largest of the studied cities, the high level of carbon content is associated rather with the accumulation in the surface soil horizons of various oil products for residential and industrial zones. organic matter input are highly limited and the soil-plant chain is disrupted, the C/N ratio can fluctuate greatly [61]. Additionally, high values of the C/N ratio may be associated with the presence of aromatic hydrocarbons and petroleum products in anthropogenically disturbed soils [62].
According to our results ( Figure 5), we can see that the level of carbon, nitrogen and their ratios are highly variable. For the city of Salekhard, as the largest of the studied cities, the high level of carbon content is associated rather with the accumulation in the surface soil horizons of various oil products for residential and industrial zones. In the park functional zone of the city, the C/N ratio is 17, which indicates a low rate of mineralization of organic matter and a lack of nitrogen in the soil profile (as leaf litter is removed and it is not retained on the soil surface). In the mature zone, on the contrary, the C/N ratio is 7, there is no deficiency of nitrogen, as the processes of accumulation of organic matter is not disturbed by anthropogenic influences.

Trace metals content
The soil cover, due to its sorption capacity, accumulates various chemical compounds from artifacts of anthropogenic impact [63,64]. Soils may accumulate various chemical compounds over time and act as geosphere-anthroposphere indicators of changes in environmental properties [65]. Soil absorption capacity depends on many factors, such as local features of soil-forming processes, the enrichment of the surface horizons with organic matter, as well as their particle size distribution [66][67][68][69]. It is known that soils of heavy (with a predominance of clay particles in the fine-grained soil) particle size distribution have a higher sorption capacity [70][71][72].
The accumulation of trace metals in the soil profile is one of many negative environmental factors of anthropogenic impact on ecosystems [24,73,74]. The harmful effects of trace metals on human health and their effects on plant productivity are well understood, In the park functional zone of the city, the C/N ratio is 17, which indicates a low rate of mineralization of organic matter and a lack of nitrogen in the soil profile (as leaf litter is removed and it is not retained on the soil surface). In the mature zone, on the contrary, the C/N ratio is 7, there is no deficiency of nitrogen, as the processes of accumulation of organic matter is not disturbed by anthropogenic influences.

Trace Metals Content
The soil cover, due to its sorption capacity, accumulates various chemical compounds from artifacts of anthropogenic impact [63,64]. Soils may accumulate various chemical compounds over time and act as geosphere-anthroposphere indicators of changes in environmental properties [65]. Soil absorption capacity depends on many factors, such as local features of soil-forming processes, the enrichment of the surface horizons with organic matter, as well as their particle size distribution [66][67][68][69]. It is known that soils of heavy (with a predominance of clay particles in the fine-grained soil) particle size distribution have a higher sorption capacity [70][71][72].
The accumulation of trace metals in the soil profile is one of many negative environmental factors of anthropogenic impact on ecosystems [24,73,74]. The harmful effects of trace metals on human health and their effects on plant productivity are well understood, and some heavy metals also have carcinogenic and mutagenic effects [75,76]. Trace metals are satellites of industrial production processes, combustion of fossil fuels in the engines of motor vehicles and from various fuel and energy enterprises [65,77,78].
The level of soil contamination by trace metals in various urban functional zones is an indicator of excessive anthropogenic impact both in the past and at the present time [79]. The results obtained in the process of our investigation reflect not only the results of the anthropogenic pressure stretched over time on urban ecosystems, but also confirm the working hypothesis about the differences in the sorption capacities of soils of fundamentally different particle size distribution. Under conditions of climate homogeneity in the studied area, the processes of soil formation (and as a consequence, the sorption capacity of soils to accumulate trace metals) are mainly determined by lithogenic conditions, landscape features [80], and the time factor of anthropogenic pressure. Urban conditions in any part of the world can hardly be called homogeneous in terms of micro-landscape diversity. Differences are visible even at the level of functional zones, if the park area represents landscapes close to natural (but still modified in the process of urban development), the residential and industrial areas represent a completely transformed landscape, subjected to constant anthropogenic load.
As can be seen in Figure 6, among all the cities we investigated, the content of trace metals in the soils of Salekhard stands out, and in all the studied functional zones. Compared to other cities, the content of Cr, Cu, Ni, Zn and As is higher in Salekhard soils. Differentiation of metal concentrations in the functional zones of Salekhard is marked by high concentrations of Cu in the residential zone, Cd in the park zone, a very high (>120 mg × kg −1 ) concentration of Zn in the residential zone and an increased concentration of As in the industrial zone. For other cities, the highest levels of metal content were recorded mainly in industrial and residential areas. The highest Cd concentrations were found in the industrial areas of Nadym, Urengoy and Gaz Sale. A similar situation was observed for Cu concentrations, with the highest values observed in industrial and residential zones.
Particular attention should be paid to the high values of metal concentrations in the mature areas. In Gaz Sale, the maximum concentrations of Cr, Ni, Cu, Zn, As and Pb are registered exactly in the mature zones, such a situation may be due to both the past raw anthropogenic impact, and the peculiarities of the selected location for sampling. Samples were taken close to a major road, and metals could get into the soil with runoff from the roadway.
To understand the variability in metal content within each functional zone (not including cities), violin diagrams of the content of each studied element within each studied zone were plotted (Figure 7). It is well seen that the concentration of trace metals mainly ranges from 0.1 to 10 mg × kg −1 (except for Cd). The peak values occur in the soils of the functional zones of Salekhard ( Figure 6). Thus, the variation in metal concentrations between the functional zones of the young cities (Urengoy, Nadym, Gaz Sale, Novy Urengoy) is not so pronounced. Particular attention should be paid to the high values of metal concentrations in the mature areas. In Gaz Sale, the maximum concentrations of Cr, Ni, Cu, Zn, As and Pb are registered exactly in the mature zones, such a situation may be due to both the past raw anthropogenic impact, and the peculiarities of the selected location for sampling. Samples were taken close to a major road, and metals could get into the soil with runoff from the roadway.
To understand the variability in metal content within each functional zone (not including cities), violin diagrams of the content of each studied element within each studied zone were plotted (Figure 7). It is well seen that the concentration of trace metals mainly ranges from 0.1 to 10 mg × kg −1 (except for Cd). The peak values occur in the soils of the functional zones of Salekhard ( Figure 6). Thus, the variation in metal concentrations between the functional zones of the young cities (Urengoy, Nadym, Gaz Sale, Novy Urengoy) is not so pronounced. We analyzed the literature for studies on the content of metals in natural and anthropogenically loaded soils in the Yamal region. Analysis of trace metal content in soils according to literature sources (Table 3) shows a high differentiation of concentrations both within individual cities and between urbanized and pristine areas. The largest differences for undisturbed areas are reflected in the work [50,81], concentrations of Pb, Cr, As and Zn vary greatly within a single zone of the Gydan Peninsula [82] recorded strong differences in Pb and Cu content in the Kharsaim city. In most other sources, metal concentrations do not widely differ within the same study area.
Compared with the data obtained in our study, the greatest differences are noted in the concentrations of As, especially with the results obtained for the city of Salekhard. For the soils of Novy Urengoy, our data are very different from those previously published, our recorded Zn content in the industrial zone more than two times higher than the published values for urban soils. Cr and Ni concentrations also strongly differ by more than twice. Combining the literature and our data we can say that our recorded concentrations of trace metals in urban soils are predominantly higher compared to those previously published, especially differing concentrations of their industrial and residential functional zones we studied. We analyzed the literature for studies on the content of metals in natural and anthropogenically loaded soils in the Yamal region. Analysis of trace metal content in soils according to literature sources (Table 3) shows a high differentiation of concentrations both within individual cities and between urbanized and pristine areas. The largest differences for undisturbed areas are reflected in the work [50,81], concentrations of Pb, Cr, As and Zn vary greatly within a single zone of the Gydan Peninsula [82] recorded strong differences in Pb and Cu content in the Kharsaim city. In most other sources, metal concentrations do not widely differ within the same study area.
Compared with the data obtained in our study, the greatest differences are noted in the concentrations of As, especially with the results obtained for the city of Salekhard. For the soils of Novy Urengoy, our data are very different from those previously published, our recorded Zn content in the industrial zone more than two times higher than the published values for urban soils. Cr and Ni concentrations also strongly differ by more than twice. Combining the literature and our data we can say that our recorded concentrations of trace metals in urban soils are predominantly higher compared to those previously published, especially differing concentrations of their industrial and residential functional zones we studied.

Multivariate Analysis of Variance and Correlation Analysis
To support the working hypothesis of the difference between the cities in the context of metal content, a statistical backtracking of the resulting concentration data set was performed. A multivariate analysis of variance (MANOVA) was conducted, and the following factors were taken as factors: city, functional zones within cities, and the crossinfluence of these factors. As a result (Table 4), the city factor was found to be significant (p < 0.0001), and consequently the influence of the particle size distribution and duration of anthropogenic pressure on the concentration of metals in the surface horizon of urban soils. A more detailed analysis for each metal (Table 5) revealed differences (p < 0.01) between concentrations by city in concentrations of Pb, Cr, Ni, and As.
The correlation analysis revealed the relationship of concentrations of some metals with the main soil properties. As can be seen in Table 6, the content of Pb, Cr, Ni and As is in direct correlation with the content of nitrogen and carbon in the soil profile; although the correlation coefficient in most cases does not exceed 0.6 this relationship is recognized as significant. Additionally, most trace metals are inversely related to pH values, but the highest and significant correlation coefficients with pH values were found for As.

Qualitative Assessment of Soil Contamination by Trace Metals
Qualitative assessment of soil contamination in the Arctic cities showed strong differences between the cities studied.
Among all the metals studied, Cd pollution was found to be the most significant for almost all cities and their functional zones according to Igeo index calculations ( Table 7). The soils of Salekhard were found to be the most contaminated with Pb, Cr, Cd and As. For Pb and As, the character of contamination of all functional zones is assessed as moderately polluted for Cr and Cd the character of pollution is assessed as moderately to highly polluted.
For other cities, point pollution in different functional zones was detected: Gaz Sale, mature zone-Cr-unpolluted to moderately polluted; Nadym, residential, park and industrial zones-Pb and As-unpolluted to moderately polluted. Additionally, in Novy Urengoy, the pollution by Zn and As in industrial and residential zones is estimated as unpolluted to moderately polluted. Application of the pollution index mainly showed results identical to the Igeo index, the soils of Salekhard city and its functional zones were recognized as the most polluted. The character of soil contamination in most cases is assessed as strong. The maximum values are revealed for Cr in all functional zones of the city of Salekhard, the nature of the occultation is assessed as very strong.
Application of the Pollution Load Index (PLI) showed ( Table 8) that soils in all functional zones of the city of Salekhard, mature zone of Gaz Sale and the industrial zone of Novy Urengoy are rated as deterioration of soil quality.

Conclusions
Although the traditional understanding of soil formation according to states of Dokuchaev and his followers is the equality of soil-forming factors, among which climate is a key determinant, in the Arctic under homogeneous humid climate conditions more attention should be paid to the lithogenic diversity and its influence on the soil cover. The investigation of soil urban ecosystems is largely complicated by a complex system of anthropogenic transformations of local landscapes in the territory of a single city. Thus, although the studied soils are in autonomous conditions at the level of the region, anthropogenic transformations of the city make it impossible to recognize them as truly autonomous. It is necessary to speak about heteronomous urban conditions of soil formation in humid climatic settings. Therefore, soil formation in urban ecosystems will be less diverse as compared to pristine areas. Taking into account similar climatic conditions and based mainly on the difference in the lithogenic structure of the bedrock and the time factor of anthropogenic impact, our study shows the relationship between these factors and the nature of accumulation of trace metals in the soils we investigated.
The city of Salekhard, whose territory is predominantly located on clay alluvial deposits, and which has a long history of anthropogenic pressure is logically characterized by a high content of trace metals in the topsoil. The absorption capacity of clay minerals backed up by a long-term anthropogenic load explains the high content of metals in the soils of Salekhard in all the studied functional urban zones.
Younger cities located on sandy textured parent materials are characterized by a lower degree of soil contamination with heavy metals. Evidently the lower sorption capacity of sandy soils can be offset by extreme levels of anthropogenic load, this is clearly seen in the content of metals in industrial areas. Local maximum content of Cu is observed in the industrial zone of Gaz Sale and Novy Urengoy, the highest Cd concentration is recorded in the soils of the industrial zone of Nadym, the concentration of Ni in the industrial zone of Novy Urengoy is also the highest. However, these local maximums in some cases are quantitatively lower compared to all studied functional zones of Salekhard.
The situation described above also found validation in the statistical analysis of the dataset. Significant statistical differences were found between the soils of different cities (p < 0.0001), but not the functional zones within the city. The cross-influence of the city and functional zone was also found to be not significant. Qualitative assessment of the pollution level is consistent with the situation described above; the application of various soil quality indices showed that the soils of Salekhard city are under significant anthropogenic pressure. While, for the other studied cities the level of soil contamination is characterized as minor (except for Cd contamination, it was found for almost all the functional zones of all the studied cities). Additionally, according to the calculated value of the index PI for some functional zones of the young cities' pollution of As is at medium and high levels. Particular attention should be paid to the vicinity of Gaz Sale city, here, according to PI index calculations, the character of Pb, As and Cr pollution is qualified as high, probably it is connected with the sampling place, which is 150 meters from the road, and there is a migration of these elements with emissions from motor vehicles.

Data Availability Statement:
The data that support the findings of this study are available from the corresponding author, upon reasonable request.