The Origin of the Natural Water Chemical Composition in the Permafrost Region of the Eastern Slope of the Polar Urals

This article presents the results of the study of water chemical composition and formation processes in the eastern slope of the Polar Ural where permafrost is widely spread. To date, this region has not been studied in detail. However, it is very important to have information on the natural waters in this region because they play a significant role in all geochemical processes, including climate formation. For this study, 107 water samples were collected from lakes, rivers and active layer waters. The studied waters are ultrafresh; total dissolved solids vary from 14 to 438 mg/L. pH ranges from 3.5 to 9.0. The chemical type of the studied waters is mostly HCO3–Ca–Mg. The geochemical environment of the studied waters contributes to the accumulation of such trace elements as Fe, Mn, Al, Ni, Au, Co, Li, Sc, Ti, Cr, Sr, Nb, Mo, Cs, La, Eu, Lu, Hg, and Se. The chemical composition peculiarities of active layer waters, rivers, lakes and waters of stone pits and exploration trenches and the description of their chemical composition formation processes are given in the present study. The water–rock interaction is a dominant process of water chemical composition formation in the studied region. The obtained results have both potential theoretical and practical applications.


Introduction
The Polar Urals is a region of northern Russia on which great hopes are pinned for the development and renewal of the mineral resource base of the Russian Federation [1][2][3].This remote region has not been studied well, particularly with respect to hydrogeochemical aspects.The specific feature of this region is widespread permafrost, as in other northern regions including Canada and the USA.
The formation of the natural water composition in the northern regions, where permafrost is widely spread, occurs under specific conditions.The natural waters of this region are in close contact with the atmosphere; they undergo seasonal phase changes and they are in solid phase most of the year.These and other factors determine the peculiarities of water composition formation processes in the permafrost region.The problems of water geochemistry in regions with a cold climate and permafrost have been well studied.
The geochemical processes occurring in the temperate region also operate in the permafrost region [4].Cationic denudation in the permafrost regions is substantially higher than the world average [5] indicating the intensity of chemical weathering in such an environment [6][7][8][9].Evidently the low average temperatures do not inhibit chemical weathering reactions.
According to [10][11][12][13] the main geochemical processes that control both surface water and groundwater chemistry in the studied regions are silicate weathering: the dissolution of freshly weathered, easily soluble rocks (mainly gypsum, anhydrite, calcite, dolomite, halite and limestone) sometimes combined with sulfide oxidation [14].Moreover, the weathering rate is determined by the permafrost degradation factor [15,16].Soluble cations in permafrost and the active layer of static cryosols at recently disturbed sites were two orders of magnitude higher than in the active layer at undisturbed sites [16].
The detailed analysis [17] of Sr isotopic composition and river chemical composition shows that silicate weathering reactions are not the only mechanisms that control solute concentrations.An atmospheric component, constituted by the dissolution of evaporate and carbonate aerosols, is necessary to explain the dispersion of chemical ratios such as Ca/Na, Mg/Na, Sr/Na and Cl/Na.These aerosols are probably of local origin.
Different aspects of hydrogeochemistry of permafrost regions are described in these works including specificity of water-rock (sulphate, carbonate, aluminosilicates) interaction and water chemical composition formation during phase transformations etc.The study of water chemical composition formation based on the research of water-rock interaction is up-to-date.Not enough attention has been paid to the role of dissolved organic matter in the processes of water chemical composition formation in the permafrost regions.Data on water enrichment in dissolved organic matter in the arctic regions are given in some articles [18][19][20][21][22][23][24][25].Most of them focus on studying river water composition, while less attention is paid to groundwater research, and lake water chemical composition formation has not been studied at all.The lakes of the arctic zone are unique objects worth studying and there are some publications on lake water balance characteristics and their influence on permafrost [26,27].
Little work has been done on the geochemistry of water of the Polar Urals.Previous studies on lake waters circulating in this area have been carried out [28] and aimed at investigating the distribution and size of lakes on the eastern slope of the Polar Urals, but there is no information on the origin of lake water chemical composition.Nevertheless, the lake, river and groundwater geochemistry of the arctic regions has not been well studied yet.
The purpose of this work is to study the geochemical characteristics of water in order to identify the main hydrogeochemical processes controlling the chemical composition based on major ion chemistry.Another aim of the study is to assess the concentration of trace elements in lakes, rivers and groundwater on the eastern slope of the Polar Urals in order to evaluate the regional geochemical background.This study would be useful for conducting hydrogeochemical prospecting for ore deposits on the eastern slope of the Polar Urals, and for investigating sources of drinking-water.

Study Area
The study site covers an area of approximately 170 km 2 on the eastern slope of the Polar Urals, between latitudes 67 • 36 29 N and 67 • 46 48 N and longitudes 68 • 08 50 E and 68 • 25 58 E (Figure 1).The morphology status of the region is a denudation plain; the elevation range is from 55 m to 202 m.The studied area covers the watershed of the rivers Enzoryakha and Yunyakha.The climate of the studied area is subarctic.The mean annual air temperature is about −5 • C. The annual precipitation is 400-600 mm.In the region, the continuous permafrost thickness varies between 100 and 320 m [29].The thickness of the active layer, the layer of soil above the permafrost thawing out during the summer, reaches a maximum 2.5 m.From the geological point of view, the studied area consists of two structure levels such as the Paleozoic and the Mesozoic (Figure 2).The Paleozoic structural level is represented by three structural stages such as the Lodlow-lower Eifelian, the upper Eifelian-Givetian and the Givetianmid-Carboniferous. The Lodlow-lower Eifelian stage is composed of the rhyolite-andesit-basaltic unit intruded by gabbroid.The upper Eifelian-Givetian stage lies unconformably and is represented by bituminous limestones intruded by gabbro-diorite unit.The Givetian-mid-Carboniferous stage is represented by terrigenous rocks.The Mesozoic stratified rocks are represented by the Carnian-Norian stage of the Triassic period [31].

Sampling and Analysis
The waters were sampled for chemical characterization during the summer period with 107 samples collected from active layer waters (32 samples), rivers (29 samples), lakes (34 samples) and waters of stone pits and exploration trenches (12 samples).The locations of the water sampling points are shown in Figure 1.At each point, two water samples were collected.The samples for ionic analyses were stored in 500 ml clean polyethylene bottles which were rinsed twice by the sampled water.The samples for trace elements analyses were stored in 50 mL clean polyethylene bottles.Each sample was filtered in situ through a 0.45 μm acetate cellulose membrane filter of 30 mm in diameter using a pre-cleaned syringe.First, 5 mL of the filtered water was eliminated at the beginning of filtration.The samples were immediately acidified with concentrated HNO3 in order to prevent metal precipitation.
The field measurements including the water temperature, pH and redox potential (Eh) were carried out using a portable device Water Test (Hanna Instruments, Cluj-Napoca, Romania).The pH measurements were performed after calibration with pH 4.01 (HI 7004) and 7.01 (HI 7007) buffer solutions.The accuracy of the instrument is ±0.2 pH and ±5 mV.The concentration of HCO3 was determined by titration with 0.01 N solution of HCl against methyl orange indicator.The titration method was used for the determination of the fulvic and humic acid concentration.Cations and anions, including SO4 2− ,Cl − , Ca 2+ , Mg 2+ , Na + , K + were analyzed by Dionex 1000 and 2000 ion The groundwater stratification is not so much determined by geological structure and water bearing rocks as by permafrost thickness and framework and its temperature regime [30].According to the groundwater localization in relation to permafrost in the studied area, the following types of groundwater can be distinguished: suprapermafrost (water saturated part of the active layer), intra-and subpermafrost waters.In the present study only suprapermafrost waters are described, which refer to the active layer waters in the present research.
Several large rivers flow in the studied area, such as the Nganotaykha, Nerutsyakha, Khohoreytanyo and Neyashor; all of these are within the Ob basin.These rivers are fed by snowmelt waters, precipitation and permafrost thawing.Besides these rivers, about 140 lakes are situated here.They are located in clusters having different sizes.The water surface of the largest lake is 0.87 km 2 , other lakes have smaller water surfaces which are less than 0.25 km 2 .The depth of the lakes is not more than 2-3 m, while the depth of a few shallow lakes is less than 1 m.Only few lakes are connected to rivers.
From the geological point of view, the studied area consists of two structure levels such as the Paleozoic and the Mesozoic (Figure 2).The Paleozoic structural level is represented by three structural stages such as the Lodlow-lower Eifelian, the upper Eifelian-Givetian and the Givetian-mid-Carboniferous.The Lodlow-lower Eifelian stage is composed of the rhyolite-andesit-basaltic unit intruded by gabbroid.The upper Eifelian-Givetian stage lies unconformably and is represented by bituminous limestones intruded by gabbro-diorite unit.The Givetian-mid-Carboniferous stage is represented by terrigenous rocks.The Mesozoic stratified rocks are represented by the Carnian-Norian stage of the Triassic period [31].
chromatographs.Dissolved trace elements were determined by inductively coupled mass spectrometry (ICP-MS, ELAN-DRC-e, Perkin Elmer, Shelton, CT, USA).All the standard solutions for it were prepared with ultra-pure deionized water and Perkin Elmer Multi Element Standard Solutions.The physico-chemical modelling program HydroGeo [32] was used for calculation of the ion activities in the studied waters.To show the water saturation state with respect to key mineral phases, the diagrams of some secondary and primary mineral stability fields were used.Statistica, version 8.0, was used to perform statistical analyses.Correlation was estimated using the Pearson correlation coefficient (R-coefficient) with statistical significance p = 0.05.The mean content of chemical elements was determined according to their distribution laws with the use of excess and asymmetry coefficients.

General Characteristics of Waters
The results of the chemical composition analysis and field measurement are given in Appendix A. Total dissolved solids (TDS) range from 14 to 438 mg/L.Most waters (at 94 points) of the studied area are fresh and have TDS of less than 200 mg/L.pH values are from 3.5 to 9.0 (Figure 3).Most waters in the region have pH 6.5-7.5.The lowest pH values 3.5-4.9were obtained in a wetland area.Lake waters are mostly alkaline with pH ranging from 8.5 to 9.2.

Sampling and Analysis
The waters were sampled for chemical characterization during the summer period with 107 samples collected from active layer waters (32 samples), rivers (29 samples), lakes (34 samples) and waters of stone pits and exploration trenches (12 samples).The locations of the water sampling points are shown in Figure 1.At each point, two water samples were collected.The samples for ionic analyses were stored in 500 mL clean polyethylene bottles which were rinsed twice by the sampled water.The samples for trace elements analyses were stored in 50 mL clean polyethylene bottles.Each sample was filtered in situ through a 0.45 µm acetate cellulose membrane filter of 30 mm in diameter using a pre-cleaned syringe.First, 5 mL of the filtered water was eliminated at the beginning of filtration.The samples were immediately acidified with concentrated HNO 3 in order to prevent metal precipitation.
The field measurements including the water temperature, pH and redox potential (Eh) were carried out using a portable device Water Test (Hanna Instruments, Cluj-Napoca, Romania).The pH measurements were performed after calibration with pH 4.01 (HI 7004) and 7.01 (HI 7007) buffer solutions.The accuracy of the instrument is ±0.2 pH and ±5 mV.The concentration of HCO 3 was determined by titration with 0.01 N solution of HCl against methyl orange indicator.The titration method was used for the determination of the fulvic and humic acid concentration.Cations and anions, including SO 4 2− ,Cl − , Ca 2+ , Mg 2+ , Na + , K + were analyzed by Dionex 1000 and 2000 ion chromatographs.Dissolved trace elements were determined by inductively coupled mass spectrometry (ICP-MS, ELAN-DRC-e, Perkin Elmer, Shelton, CT, USA).All the standard solutions for it were prepared with ultra-pure deionized water and Perkin Elmer Multi Element Standard Solutions.
The physico-chemical modelling program HydroGeo [32] was used for calculation of the ion activities in the studied waters.To show the water saturation state with respect to key mineral phases, the diagrams of some secondary and primary mineral stability fields were used.Statistica, version 8.0, was used to perform statistical analyses.Correlation was estimated using the Pearson correlation coefficient (R-coefficient) with statistical significance p = 0.05.The mean content of chemical elements was determined according to their distribution laws with the use of excess and asymmetry coefficients.

General Characteristics of Waters
The results of the chemical composition analysis and field measurement are given in Appendix A. Total dissolved solids (TDS) range from 14 to 438 mg/L.Most waters (at 94 points) of the studied area are fresh and have TDS of less than 200 mg/L.pH values are from 3.5 to 9.0 (Figure 3).Most waters in the region have pH 6.5-7.5.The lowest pH values 3.5-4.9were obtained in a wetland area.Lake waters are mostly alkaline with pH ranging from 8.5 to 9.2.The redox conditions are characterized by values Eh −68 to 347 mV.There is a weakly positive relationship between Eh and pH (R 2 = 0.0032) (Figure 3) in the studied waters.This pH and Eh relationship is not typical for the natural water; it is caused by the wide variety of the chemical composition of the studied waters.The most common waters in the studied area have the values of redox potential of 150-200 mV.More oxygenated environments (where Eh is greater than 200 mV) are formed mainly in acidic and neutral waters in watershed areas and upstreams, where water exchange is more active.More reducing conditions (where Eh is below 200 mV) are formed in water with low water exchange intensity.The enrichment processes of organic substances in water play a significant role in the formation of the reduction environment [33].
The concentrations of fulvic and humic acids at a number of points are shown in Table 1.The fulvic acid concentrations in waters vary from 1.3 to 19.4 mg/L and humic acid concentrations vary from 0.5 to 2.8 mg/L.The highest fulvic acid concentrations are found in waters with a reducing environment where Eh value is less than 200 mV (Figure 4).The maximum concentrations of fulvic acid are noted in the most acidic waters with a pH less than 5.5.The fulvic acid concentrations decrease with an increase in pH.The behavior of soluble humic acid is practically independent of geochemical parameters of the aquatic environment [33,34].The redox conditions are characterized by values Eh −68 to 347 mV.There is a weakly positive relationship between Eh and pH (R 2 = 0.0032) (Figure 3) in the studied waters.This pH and Eh relationship is not typical for the natural water; it is caused by the wide variety of the chemical composition of the studied waters.The most common waters in the studied area have the values of redox potential of 150-200 mV.More oxygenated environments (where Eh is greater than 200 mV) are formed mainly in acidic and neutral waters in watershed areas and upstreams, where water exchange is more active.More reducing conditions (where Eh is below 200 mV) are formed in water with low water exchange intensity.The enrichment processes of organic substances in water play a significant role in the formation of the reduction environment [33].
The concentrations of fulvic and humic acids at a number of points are shown in Table 1.The fulvic acid concentrations in waters vary from 1.3 to 19.4 mg/L and humic acid concentrations vary from 0.5 to 2.8 mg/L.The highest fulvic acid concentrations are found in waters with a reducing environment where Eh value is less than 200 mV (Figure 4).The maximum concentrations of fulvic acid are noted in the most acidic waters with a pH less than 5.5.The fulvic acid concentrations decrease with an increase in pH.The behavior of soluble humic acid is practically independent of geochemical parameters of the aquatic environment [33,34].Among inorganic anions bicarbonate is predominant, the concentrations of which range from 2.6 to 325 mg/L.Waters with HCO3 concentrations less than 150 mg/L are widely spread.The second most abundant anion is SO4, ranging from 0.5 to 50 mg/L.The Cl content is low (less than 0.3-0.5 mg/L).Ca and Mg are the most abundant cations with concentrations of 1.6-91.7 mg/L and 0.7-21.2mg/L, respectively (Appendix A).The concentrations of Na range from 0.4 to 15.9 mg/L.K is the least abundant major cation with a content of 0.01-7.1 mg/L.
Hydrogeochemical types of the studied waters are shown in Figure 5. Thus, the HCO3 contents are mostly more than 70% of the total anion, except for a few data points (50%-70%).At these points, the SO4 concentration (30%-50%) increases and these waters belong to HCO3-SO4 type.The waters with a higher concentration of SO4 (10%-50%) are widely spread in the north-western part of the studied area.The Cl content increase in waters in the north-western part of the study are.Among inorganic anions bicarbonate is predominant, the concentrations of which range from 2.6 to 325 mg/L.Waters with HCO 3 concentrations less than 150 mg/L are widely spread.The second most abundant anion is SO 4 , ranging from 0.5 to 50 mg/L.The Cl content is low (less than 0.3-0.5 mg/L).Ca and Mg are the most abundant cations with concentrations of 1.6-91.7 mg/L and 0.7-21.2mg/L, respectively (Appendix A).The concentrations of Na range from 0.4 to 15.9 mg/L.K is the least abundant major cation with a content of 0.01-7.1 mg/L.
Hydrogeochemical types of the studied waters are shown in Figure 5. Thus, the HCO 3 contents are mostly more than 70% of the total anion, except for a few data points (50%-70%).At these points, the SO 4 concentration (30%-50%) increases and these waters belong to HCO 3 -SO 4 type.The waters with a higher concentration of SO 4 (10%-50%) are widely spread in the north-western part of the studied area.The Cl content increase in waters in the north-western part of the study are.
The cation composition mostly varies between Ca, Ca-Mg, Mg-Ca and Ca types.A significant feature of cation composition of the studied waters is largely a function of rock composition.Mostly magnesium waters (with Mg of more than 50%) are located in the area with widely spread granite and gabbro.Calcium waters (with Ca of more than 50%) are widely spread in the area with carbonates.Na concentration is generally less than 10% in the studied waters.The waters with Na concentration of more than 10% have low TDS and they are located in the northern part of the studied area.
abundant major cation with a content of 0.01-7.1 mg/L.
Hydrogeochemical types of the studied waters are shown in Figure 5. Thus, the HCO3 contents are mostly more than 70% of the total anion, except for a few data points (50%-70%).At these points, the SO4 concentration (30%-50%) increases and these waters belong to HCO3-SO4 type.The waters with a higher concentration of SO4 (10%-50%) are widely spread in the north-western part of the studied area.The Cl content increase in waters in the north-western part of the study are.The silicon concentrations vary from 0.16 to 28.75 mg/L.A significant feature of these waters is hydrolysable cation accumulation such as Al, Fe, and Mn.Their mean concentrations are high for waters with such low TDS and reach 0.03, 0.25, and 0.01 mg/L respectively.
The average concentrations of trace elements in the studied waters correspond to both the average of world rivers and mean contents in groundwaters of tundra landscapes [33].The average concentrations of Sc, Eu, Fe, Ho, Dy, Tb, Gd, Pr, Er, Yb, Sm, Se, Nd, Lu, Nb, Tm, Ce, Mn, La, and Ni in the studied waters exceed the average of world rivers by 1.5-10 times (Table 2).
The trace element concentrations in the studied waters are lower than those in waters of tundra landscape; exceptions are Sc, Fe, Mn, and Ni concentrations.The Cu, V, As, Sr, I, and Ba concentrations in the studied waters are close to those in waters of tundra landscape.Moreover, new information about Bi, Cd, Pt, B, Ge, Se, Br, Y, Rh, Pd, Te, Cs, TR, Hf, Ta, W, Re, Ir, Hg, Tl, and Th distribution in the natural waters in the eastern slope of the Polar Urals was found (see Table 2).
Figure 6a,b show the concentration of the main ions plotted against TDS.On the whole, the data show a strong positive linear correlation pointing to the significant dependence of HCO 3 and Ca on TDS (R = 0.97 and R = 0.90, p = 0.05 respectively).concentrations in the studied waters are close to those in waters of tundra landscape.Moreover, new information about Bi, Cd, Pt, B, Ge, Se, Br, Y, Rh, Pd, Te, Cs, TR, Hf, Ta, W, Re, Ir, Hg, Tl, and Th distribution in the natural waters in the eastern slope of the Polar Urals was found (see Table 2).Figures 6a,b show the concentration of the main ions plotted against TDS.On the whole, the data show a strong positive linear correlation pointing to the significant dependence of HCO3 and Ca on TDS (R = 0.97 and R = 0.90, p = 0.05 respectively).The dependence of the calcium concentration on the bicarbonate concentration, with a positive connection (R = 0.87, p = 0.05), is shown in Figure 6c.However, water enrichment in these ions is not determined only by simple calcite dissolution but it is controlled by the processes of aluminosilicates hydrolyses because this type of rock is widely spread in the studied area (Figure 6d).The exception is waters in stone pits and exploration trenches.Calcium entering in these waters occurs primarily due to the calcium dissolution processes.The dependence of the calcium concentration on the bicarbonate concentration, with a positive connection (R = 0.87, p = 0.05), is shown in Figure 6c.However, water enrichment in these ions is not determined only by simple calcite dissolution but it is controlled by the processes of aluminosilicates hydrolyses because this type of rock is widely spread in the studied area (Figure 6d).The exception is waters in stone pits and exploration trenches.Calcium entering in these waters occurs primarily due to the calcium dissolution processes.
No strong dependence of other ions on TDS was found (Table 3).However, there are positive correlations between K and Al, K and Si, demonstrating likely their common source.In Table 3, a positive correlation between Al and Si (R = 0.85, p = 0.05) is observed in the study waters: concentrations of Al increase with the growth of Si.

Chemical Composition of Water Types
Chemical composition variety of the studied waters is determined by various existence factors and processes in the study area.The chemical compositions of rivers, lakes and groundwaters (waters of active layer) form in natural conditions of landscape.At the same time, chemical composition of waters in stone pits and exploration trenches form in other hydrodynamic and physico-chemical conditions, which are mostly the result of anthropogenic activity.Meanwhile, the specific feature of the study area is the absence of a great anthropogenic influence.Therefore, it is necessary to describe each type of the studied waters separately.
Waters of the active layer are presented by small swampy land sometimes with water table, but more often by water-bearing peat soils.These waters have a color similar to "strong tea", which is conditioned by high concentrations of organic matter.Active layer waters are mostly weak-acid or neutral, pH values range from 3.5 to 7.5.Their TDS vary from 24.1 to 273.0 mg/L, average TDS is 96.2 mg/L (Appendix A).On average, the chemical composition of active layer waters is close to that in river waters.The values of TDS were estimated using only main inorganic ions, but for this type of water, organic matter plays a significant part in chemical composition.In this manner, values of TDS can be increased more than 10%.The studied waters are characterized by HCO 3 type with various ratios of Ca and Mg.
The high concentrations of such elements as Si, Ag, Bi, Cu, Cd, Pb, Be, Sc, Rb, Y, Nb, Te, Cs, Eu, Lu, Hf, Ta, W, Th, V, Ga, Ge, As, I, and Ba are observed in this type of water (Table 2).
It should be noted that the concentrations of Mn, Co, Zn, Pd, La, Tl, Sn, Fe, and Ir are significantly higher in active layer waters than those in rivers and lakes.
The lake waters are ultrafresh; the maximum value of TDS is less than 200 mg/L.Geochemical environment of lake waters is characterized by pH from 5.3 to 9.2 with an average value of 7.3, and Eh from 47 to 341 mV with an average value of 175 mV.The values of these parameters are very close to those in rivers in the studied area.Calcium and magnesium are the dominant cations in lake waters; together they account for 72%-92% of the cations measured (Figure 5).Anion chemistry is dominated by HCO 3 constituting from 80% to 99% of the major anions.Higher concentrations are characteristic of some trace elements in lake waters among which are lithophilous elements-Li, Be, B, Ti, Cs, La, Eu, I, Ba and chalkophilous elements-Pb, Rb, Pd, Te, Hg, Tl, and As (Table 2).For siderophilous elements, only Ni reveals a higher content.The high concentration of Ni in the waters is a one of the specific features of the studied area [35].
The river waters are characterized by values of Eh from 65 to 225 mV.Measured pH values range from 4.8 to 8.2.The average value of TDS is 80 mg/L, which is equal that of active layer waters.The maximum value of TDS is 319.2 mg/L.River waters belong to HCO 3 -Ca hydrogeochemical type.Geochemical conditions in river waters are appropriate for enrichment in lithophilous elements-Si, Al, Mn, Sc, Cr, Y, Zr, Nb, Hf, Ta, W, and U, chalkophilous elements-Au, Zn, Cu, Cd, Sb, V, and Se, and siderophilous elements-Co, Pt, Mo, Ru, Rh, Ir, Re, and Ge (Table 2).
Waters in stone pits and exploration trenches are characterized by different chemical composition in all studied water types.The measured pH values are neutral to alkaline, ranging from 7.4 to 8.9.The water exhibit redox potential characteristic of a more oxygenated environment; measured values of Eh vary from 121 to 218 mV.On the whole, TDS (45 < TDS < 438 mg/L; average 218 mg/L) levels are greater than those of all studied water types.As in the previous water type, bicarbonate is the most abundant major anion; its values range from 32.9 to 330.4 mg/L.However, the highest contents of sulfate are observed in waters in stone pits and exploring trenches for all studied water types; their content ranges from 1.0 to 50.0 mg/L.The dominant cations are calcium and magnesium.Therefore, these waters belong to the HCO 3 -SO 4 -Ca-Mg type.The highest concentrations of Cl and K (up to 6.3 and up to 5.7 mg/L respectively) were found in this type of water.Also among all the studied water types, the maximum concentrations of such elements as Al, Fe, Au, Pt, Ti, Ni, Cr, Sr, Mo, Ru, Rh, Sb, Re, Ir, and U were detected in the waters of stone pits and exploring trenches.

Water-Rock Interaction
According to the dependence of TDS on Na/(Ca + Na) [36] (Figure 7), study waters are characterized by low values of Na/(Na + Ca) at less than 0.5, with moderate TDS values of less than 250 mg/L, typical of rock-dominated waters.Several regularities can be identified in Figure 5. On one hand, values of Na/(Na + Ca) decrease with TDS growth in active layer waters.On the other hand, a Na/(Na + Ca) value decrease does not lead to TDS increase in river and lake waters.It may be conditioned by dilution waters with precipitation.Moreover, waters in stone pits and exploration trenches with minimum values of Na/(Na + Ca) are characterized by maximum values of TDS.Therefore, in the studied area, water chemical composition formation is mostly dominated by rock weathering processes.The point positions in Figure 7 confirm this fact.
studied water types; their content ranges from 1.0 to 50.0 mg/L.The dominant cations are calcium and magnesium.Therefore, these waters belong to the HCO3-SO4-Ca-Mg type.The highest concentrations of Cl and K (up to 6.3 and up to 5.7 mg/L respectively) were found in this type of water.Also among all the studied water types, the maximum concentrations of such elements as Al, Fe, Au, Pt, Ti, Ni, Cr, Sr, Mo, Ru, Rh, Sb, Re, Ir, and U were detected in the waters of stone pits and exploring trenches.

Water-Rock Interaction
According to the dependence of TDS on Na/(Ca + Na) [36] (Figure 7), study waters are characterized by low values of Na/(Na + Ca) at less than 0.5, with moderate TDS values of less than 250 mg/L, typical of rock-dominated waters.Several regularities can be identified in Figure 5. On one hand, values of Na/(Na + Ca) decrease with TDS growth in active layer waters.On the other hand, a Na/(Na + Ca) value decrease does not lead to TDS increase in river and lake waters.It may be conditioned by dilution waters with precipitation.Moreover, waters in stone pits and exploration trenches with minimum values of Na/(Na + Ca) are characterized by maximum values of TDS.Therefore, in the studied area, water chemical composition formation is mostly dominated by rock weathering processes.The point positions in Figure 7 confirm this fact.In the studied area, aluminosilicate rocks are mostly predominant, whereas carbonate rocks are less common.Thus, the main process of water chemical composition formation is hydrolysis of aluminosilicate minerals and less the dissolution of carbonates.These processes lead to the release of chemical elements in solution and secondary mineral precipitation.For the study of these processes, In the studied area, aluminosilicate rocks are mostly predominant, whereas carbonate rocks are less common.Thus, the main process of water chemical composition formation is hydrolysis of aluminosilicate and less the dissolution of carbonates.These processes lead to the release of chemical elements in solution and secondary mineral precipitation.For the study of these processes, the major component activities of the water samples were plotted on diagrams of some secondary and primary mineral stability fields (Figure 8).
The point positions outside the anorthite, albite, Mg-chlorite, analcim, and muscovite stability fields (Figure 8) indicate that all the analyzed waters are not in equilibrium with these minerals of primary rocks.Therefore, the studied waters are likely to promote dissolution of primary aluminosilicate minerals and constantly dissolve them.Weathering of aluminosilicate minerals is a slow process [37], releasing chemical elements in solution from primary rocks and also producing secondary minerals; removing chemical elements from solution.As previously reported [33], the processes of water-rock interaction have a stage character, one of the main factors in which is water-rock interaction time, determined by water exchange intensity.
The chemical composition formation of the studied waters occurs under different conditions of water-rock interaction.Active layer waters, rivers, lakes and waters of stone pits and exploration trenches are characterized by different time of interaction with rocks and therefore they reach different water-rock interaction stages.In Figure 8, the point positions in the gibbsite, kaolinite, montmorillonites, calcite stability fields demonstrate that the studied waters are in equilibrium with the products of primary rock weathering.Each point position in the secondary mineral stability fields shows that the studied waters reach four stages of water-rock interaction [38].
the major component activities of the water samples were plotted on diagrams of some secondary and primary mineral stability fields (Figure 8).The point positions outside the anorthite, albite, Mg-chlorite, analcim, and muscovite stability fields (Figure 8) indicate that all the analyzed waters are not in equilibrium with these minerals of primary rocks.Therefore, the studied waters are likely to promote dissolution of primary The relatively short time of interaction between waters and rocks defines the development of water-rock system until the first stage, namely water equilibrium with mineral groups of oxides, hydroxides.At this stage of water-rock interaction, the aluminosilicates hydrolyses processes can be described by following reactions: In the first stage of water-rock interaction, the primary aluminosilicates of various compositions actively dissolve by water.This leads to the release of a great quantity of chemical elements in solution.The accumulation of Al, Si and other chemical elements in waters occurs until saturation with respect to the least soluble secondary mineral-gibbsite, ferrihydrite, geothite and other minerals of the oxides, hydroxides group.According to the diagram in Figure 8, the large part of samples is located in the gibbsite stability field.This is characteristic, especially in acid lake waters and rivers and partly for active layer waters.In addition, studied waters are in equilibrium with goethite; the average value of SI is 32.
In the state of saturation with respect to gibbsite, the Al content growth slows down because of its removal from solution with the formation of secondary phase but it does continue.Accumulation of Si in solution occurs during the continued aluminosilicate hydrolyses.When the appropriate ion activity product is gained, waters reach the state of saturation with respect to kaolinite, the second stage of water-rock interaction.It can be shown by following reactions: 2Al(OH) 3  This stage is mostly reached by active layer waters and rivers at the floodplain areas.
With the Si accumulation in solution, the degree of water saturation with respect to clay minerals grows and the water-rock system can reach the next stage-saturation with respect to montmorillonites.However, the waters are oversaturated with respect to montmorillonites only in single points among all the studied water samples (Figure 8).
The behavioral peculiarity of Al and Si during the process of interaction between studied waters and rocks is plotted on Figure 9.The expected rise of Al concentration in solution occurs with the Si content growth and therefore with the increase of interaction time between water and aluminosilicate minerals (Figure 9).The intensive growth of the Al content with the Si content growth is observed for the active layer waters, the water chemistry of which is mostly determined by the water-rock interaction processes.Meanwhile, the staging of Al removal from solution during the formation of secondary mineral phase is shown in Figure 9.In the first stage of water-rock interaction there are at least three limitations for Al accumulation in the studied waters: water saturation with respect to gibbsite, kaolinite and montmorillonites.In Figure 9a, the points of the Al content or accumulation rate decreases can be noted at the Si concentrations of 2.5 and 6-15 mg/L, which correspond to the points of saturation with respect to kaolinite and montmorillonites [37].This dependence is characteristic, mostly for active layer waters.However, it should be noted that the Al accumulation rate decreases with Si content increase in river waters.In addition, there is no dependence of Al content on Si concentration in waters of stone pits and exploring trenches, because their water chemistry is dominated by mostly carbonate rock dissolution processes.The expected Si decrease in solution during secondary mineral precipitation in all mentioned stages of water-rock interaction does not occur; this is because silica concentration in primary minerals is higher than that in clay minerals [39].
mostly for active layer waters.However, it should be noted that the Al accumulation rate decreases with Si content increase in river waters.In addition, there is no dependence of Al content on Si concentration in waters of stone pits and exploring trenches, because their water chemistry is dominated by mostly carbonate rock dissolution processes.The expected Si decrease in solution during secondary mineral precipitation in all mentioned stages of water-rock interaction does not occur; this is because silica concentration in primary minerals is higher than that in clay minerals [39].The point position in Figure 8 shows that the studied waters are mostly undersaturated with respect to calcite, but the saturation is indicated in the single points, corresponding mainly to waters of stone pits and exploring trenches.In addition, alkaline lake waters are oversaturated with respect to calcite, or in a state close to saturation.The state of water saturation with respect to calcite is not typical for natural waters in the permafrost regions [33].The state of studied waters is conditioned by the influence of local factors such as interaction between water and carbonates and weak water exchange.Weak water exchange is caused by flat land surface and presence permafrost layer, which do not assist to precipitation infiltration.Active layer waters and rivers are far from saturation with respect to calcite.
The presented results of the thermodynamic calculations of water-rock interaction are proven by field observations of precipitation of the secondary minerals in the studied permafrost region or under similar conditions, presented in publications.Thus, under similar conditions in the Bolshezemelskaya tundra, the presence of highly dispersed clay minerals, aluminium and iron colloids in large amounts were found in the active layer rocks [40,41].These published data, our thermodynamic calculations and field observations show the ability of mentioned secondary mineral formation in permafrost regions.

Mechanisms of Water Chemical Composition Formation
According to our investigation, the main process of water chemical composition formation is the interaction between water and rocks.
As it was shown in the present study, the main anion in the studied water is HCO3.Water enrichment in HCO3 results the neutralization of alkalinity, forming as a result of aluminosilicates hydrolysis, by the products of organic matter mineralization [42].The point position in Figure 8 shows that the studied waters are mostly undersaturated with respect to calcite, but the saturation is indicated in the single points, corresponding mainly to waters of stone pits and exploring trenches.In addition, alkaline lake waters are oversaturated with respect to calcite, or in a state close to saturation.The state of water saturation with respect to calcite is not typical for natural waters in the permafrost regions [33].The state of studied waters is conditioned by the influence of local factors such as interaction between water and carbonates and weak water exchange.Weak water exchange is caused by flat land surface and presence permafrost layer, which do not assist to precipitation infiltration.Active layer waters and rivers are far from saturation with respect to calcite.
The presented results of the thermodynamic calculations of water-rock interaction are proven by field observations of precipitation of the secondary minerals in the studied permafrost region or under similar conditions, presented in publications.Thus, under similar conditions in the Bolshezemelskaya tundra, the presence of highly dispersed clay minerals, aluminium and iron colloids in large amounts were found in the active layer rocks [40,41].These published data, our thermodynamic calculations and field observations show the ability of mentioned secondary mineral formation in permafrost regions.

Mechanisms of Water Chemical Composition Formation
According to our investigation, the main process of water chemical composition formation is the interaction between water and rocks.
As it was shown in the present study, the main anion in the studied water is HCO 3 .Water enrichment in HCO 3 results the neutralization of alkalinity, forming as a result of aluminosilicates hydrolysis, by the products of organic matter mineralization [42].
The neutralization reaction may also occur in the presence of organic acids existing in essential amounts in waters.
The main source of SO 4 in the studied waters is the oxidation processes.It is the result of oxidation processes of sulfide ores, which are in contact with the atmosphere, especially in stone pits and exploring trenches.
The formation of active layer water chemical composition occurs under the influence of different factors and processes, with differing results.On one hand there is weak water exchange because of flat ground surface on the studied area.On the other hand, there is a process of seasonal freezing-melting of water-bearing rocks.This determines the water existence in the solid state for the most part of the year.For that reason, water-rock interaction time decreases and water enrichment in chemical elements is low.
However, in the waters of active layer there are favorable conditions for chemical element accumulation, such as acid conditions and high concentrations of organic matter.Acidic environment of the studied waters defines the water aggressivity on water-bearing rocks and organic acids contribute to formation of organic species of chemical elements.
Various chemical compositions of the studied waters predefine the peculiarity of water-rock interaction.Those active layer waters are oversaturated with respect to gibbsite and kaolinite.The degree of saturation with respect to these minerals increases with a water exchange intensity decrease-from watershed to river floodplain.Waters in watershed areas are characterized by low TDS and pH values, because the water-rock interaction time is very short.Waters on river floodplain areas have higher pH values, which is defined by water enrichment in OH − and chemical elements, forming as a result of water-rock interaction.
The lake water chemistry is mostly dominated by dissolution of carbonates.This process defines the geochemical peculiarities of lake waters such as the highest values of pH and the lowest values of TDS among all the studied waters.Moreover, more alkaline lakes with pH 7.3-9.2are located at the upper main river parts where carbonates are widely spread.These lake waters are very close to saturation with respect to calcite.This is one of the main peculiarities of lake waters.Among alkaline lakes, there are lakes with low pH (4.6-7.0).These lakes are located at the floodplain of main rivers, where lake waters enrich in organic matter, resulting in the decrease of pH.At the same time, the lowest values of TDS may be conditioned by dilution waters with precipitation.Thus, the main factors and processes of lake water chemical composition formation are carbonate dissolution, morphology position, scant quantity of organic matter, and short time of water-rock interaction.
As for waters in stone pits and exploring trenches, the main processes and factors of water chemical composition formation are weak water exchange, scant quantity of organic matter and active oxidation processes.According to the lithology of the studied area, marble rocks with sulfide inclusions are widely spread.
Oxidation of sulfide minerals causes water enrichment with the sulfate ion.The result of this process is a pH decrease [43], whereas the studied waters are characterized by high pH values.This is due to the neutralization processes accompanying water interaction with carbonated rocks according to following reactions.

−
It should be noted that waters of stone pits and exploration trenches are oversaturated with respect to calcite.This is due to the weak intensity of water exchange, assisted to reach this stage of water-rock interaction.
Thus, chemical composition of waters in stone pits and exploration trenches is mostly conditioned by natural factors which were intensified by negligible exploration activity in the studied area.

Conclusions
The present investigation shows that the studied waters are ultrafresh with TDS from 14 to 438 mg/L.Geochemical environment is characterized as typical for swamp and groundwaters.The pH values range from 3.5 to 9.0, Eh values vary from −68 to 347 mV.The fulvic acid concentrations in waters vary from 1.3 to 19.4 mg/L and humic acids from 0.5 to 2.8 mg/L.The studied waters mostly belong to HCO 3 -Ca-Mg or HCO 3 -Mg-Ca chemical water types, rarely do waters have SO 4 -HCO 3 -Ca composition.
The studied water types are namely active layer waters; rivers, lakes and waters of stone pits and exploration trenches have their own chemical composition peculiarities caused by the different conditions of chemical composition formation.For example, active layer waters are enriched with organic matter more strongly than other water types in the studied area.However, lakes and waters of stone pits and exploration trenches are more alkaline.
The geochemical environment of studied waters assists in the accumulation of such trace elements as Fe, Mn, Al, Ni, Au, Co, Li, Sc, Ti, Cr, Sr, Nb, Mo, Cs, La, Eu, Lu, Hg, and Se, with concentrations much higher than the average for rivers around the world.
The chemical composition of the studied waters is the result of chemical element enrichment mainly during the processes of alunimosilicate hydrolysis and carbonate dissolution.The presented thermodynamic calculations of water equilibrium with rocks show that the studied waters are undersaturated with respect to primary minerals and oversaturated with respect to secondary minerals such as Al, Mn, Fe oxides, hydroxides, kaolinite, rarely montmorillonites and calcite.

Conflicts of Interest:
The author declares no conflict of interest.

Figure 1 .
Figure 1.Sketch map of the studied area showing the sampling sites.Inset shows location of the study area in the Russian Federation.

Figure 1 .
Figure 1.Sketch map of the studied area showing the sampling sites.Inset shows location of the study area in the Russian Federation.

Figure 2 .
Figure 2. Lithology sketch-map of the studied area.

Figure 2 .
Figure 2. Lithology sketch-map of the studied area.

Figure 3 .
Figure 3. pH vs. redox potential (Eh) plot for the studied area.

Figure 4 .
Figure 4. Relationship of organic acid content vs. pH (a); organic acid content vs. Eh (b) for same water points from the studied area

Figure 4 .
Figure 4. Relationship of organic acid content vs. pH (a); organic acid content vs. Eh (b) for same water points from the studied area

Figure 5 .
Figure 5. Piper diagram showing anion and cation compositions of the studied waters.1-waters of active layer; 2-lakes; 3-rivers; 4-waters of stone pits and exploration trenches.

Figure 5 .
Figure 5. Piper diagram showing anion and cation compositions of the studied waters.1-waters of active layer; 2-lakes; 3-rivers; 4-waters of stone pits and exploration trenches.

Figure 7 .
Figure 7. Gibbs plots indicating the mechanisms that determine major-ion composition of the studied waters.Symbols as in Figure 5.

Figure 7 .
Figure 7. Gibbs plots indicating the mechanisms that determine major-ion composition of the studied waters.Symbols as in Figure 5.

Figure 9 .
Figure 9. Relationships Al vs. Si for lakes, rivers, waters of stone pits and exploration trenches (a) and for waters of active layer (b).Symbols as in Figure 5.

Figure 9 .
Figure 9. Relationships Al vs. Si for lakes, rivers, waters of stone pits and exploration trenches (a) and for waters of active layer (b).Symbols as in Figure 5.

Table 1 .
The physico-chemical parameters and fulvic and humic acid contents in several water points in the studied area.

Table 1 .
The physico-chemical parameters and fulvic and humic acid contents in several water points in the studied area.

Table 2 .
Trace elements content in the studied waters, µg/L ( min-max average ).

Table A1 .
Analytical Data on Water Samples: Physico-Chemical Parameters and Major Components.Redox Potential (Eh, mV), Total Dissolved Solids (TDS) and Ionic Concentration in mg/L.