Heat and Mass Transfer by Vapour in Freezing Soils †

: Vapour mass transfer is often underestimated when designing the bases for structures in frost susceptible soils. Intensive and long-term vapour transport may lead to excessive frost heaving and associated issues. A vapour transport model and the algorithm of its calculation is presented in this study based on the results of experimental freeze–thaw cycles of nine soil samples with varied density. The temperature ﬁeld distribution, air voids volume and the energy comprising latent heat for the phase transition and heat extracted during the temperature drop are the main parameters for determining the vapour velocity and the amount of ice formed. According to the results, the average speed of vapour transport in frozen soils was about 0.4 m/h. The amount of ice built in 1 h during uniaxial freezing due to the saturated vapour pressure difference was 1.64 × 10 − 5 –3.6 × 10 − 5 g/h in loose samples and 1.41 × 10 − 6 g/h to 5.61 × 10 − 7 g/h in dense samples of 10 cm diameter and 10 cm high sections. The results show that vapour mass transfer can increase the risk of ice growth and related problems.


Introduction
Frost heave poses huge geotechnical challenges for foundations in cold areas with both seasonally freezing and seasonally thawing soils in arctic regions [1]. Highways and dams are the most vulnerable structures prone to frost heaving [2]. While soil ice formation due to the migration of liquid water has been extensively studied, the role of vapour mass transfer is less well understood.
The focus for the research presented here is Nur-Sultan, Kazakhstan. This is the second coldest capital in the world and is classified as a climate zone III with a minimum day temperature of −41 • C in winter and a record low temperature of −52 • C [3]. It is also one of the most resource-intensive cities in terms of annual investment in the repair of highways and other structures with shallow foundations.
Soils in the city are represented by alluvial mid-Quaternary modern deposits with a thickness of 0.9 to 10.0 m, consisting of sandy clay soils with interlayering loamy sands and sandy clay loams, and clays and silty soils, which can be considered as a highly compressible and frost susceptible base. Eluvial formations of the weathering crust occur as sandy clay loams located beneath the sand and gravel alluvial formations at a depth of 6.0-10.0 m. Gravelly, gruss and rubble soils of eluvial formations of the weathering crust with sufficient bearing capacity have been found at a depth of 7.0 to 23.0 m. [4]. The filtration coefficient for alluvial medium-upper Quaternary modern soils ranges from 0.2-0.6 m/day for alluvial sandy clay soils, 3.8 ÷ 15.7 for medium size sands and 17.32-30.29 m/day for gravel and gravel soils. The filtration coefficient for eluvial formations underlying Quaternary sediments ranges from 0.001-0.007 m/day for sandy clay soils, 2.4 ÷ 20.0 m/day for grit-crushed stone soils and 0.37 ÷ 11.5 m/day for fractured sandstones. The sites are periodically flooded with rain and melt water and the groundwater level is relatively high and prone to frost heaving in the winter and weakening during thaws in the spring. The increased density and thermal conductivity of pavement materials and structures creates low-temperature fields near structures and enhances increased moisture transfer due to cryosuction forces that exacerbate frost heaving.
There are few widely recognized concepts of frost heave, notably the segregation potential perspective introduced by Konrad and Morgenstern, which implies water migration [5], and the discrete ice lens theory modified by Gilpin [6]. Both imply that moisture mass transfer occurs in a liquid state, i.e., a water-ice phase transfer. The water-ice energy balance was also considered in the coupled heat-moisture transfer, the application of which was later used in the FROSTB model [7]. Gorelik and Kolunin V.S. demonstrated that in an ice body with low hydraulic permeability all solutes and inclusions concentrate in air bubbles, which slowly move towards the warm side, where the surface tension has a reduced value [8]. This approach explains the formation of uniform blanket ice body formation containing soil particles. Arenson and Sego, in their conceptual model for ice growth in the coarse soils, observed that ice segregation in the coarse soils starts from the grain skeleton, as the thermal conductivity of the sand particles surpasses the water conductivity [9]. Moreover, the growth of ice needles starts from the sand particles' surface in a perpendicular direction and the concentrated brine was steadily moved to the pore centres [9].
Water-ice transformation modelling in fully saturated soils and moisture mass transfer in a liquid state have been widely studied and explained by Ming and Li, among others [10,11]. Henry, in her review, has summarized that frost heave is the conversion from the liquid part to the solid, which is dependent on: a removal of heat, a means of transporting the ice away from the pores and a water supply [12]. Arenson et al., discovered that vertical ice vein growth precedes the horizontal ice lens formation [13] and noted that vertical veins do not grow in thickness over time, unlike horizontal lenses. They did not identify in what phase the moisture was transported, but did note concerns about the suction required to drive the hydraulic conductivity in atmospheric pressure by determining that the negative pressure should be no less than 900 kPa to draw up the water [13]. Concerning the crack formation concept, Azmatch et al., recognised that it is still unclear whether the horizontal cracking of the soil during ice lens formation occurs due to desiccation shrinkage or thermal strain [14].
In this work, the moisture transfer in a gas state was considered. The pressure reduction in the cold temperatures on the top of the sample and the higher pressure in the warm layers induced saturated vapour transport from the high pressure towards the freezing fringe zone. Moreover, the vapour does not have a capillary surface tension; hence, its transportation is more reliable in terms of hydrodynamics and tension. This article is an extended version of Sarsembayeva et al. [15] with detailed examples of the calculation presented by the height of the samples with various density.

Materials and Methods
The laboratory testing was implemented by freeze-thaw cycles of 9 soil columns chilled from the top downwards and supplied with deionized water from the base (open system) [16]. The environmental chamber was designed based on the ASTM D 5918-06 standard [17], although the cooling rate and the sample preparation has been modified. Each soil column of 50 cm length and 10 cm diameter was compacted with the dry density varying from very loose to very dense soils, with a range of 1177 to 1800 kg/m 3 (Table 1). A mould for packing a soil sample was assembled from plastic rings 10 cm high. Gradual compaction of the samples was carried out with alternating filling of a portion of soil to ensure uniform density along the height. The samples with loose density received a smaller number of blows and energy effort per 5 cm soil layer. Each subsequent column had a higher density and accordingly more mechanical energy, i.e., hammer blows, applied during packing. After packing, each column was securely wrapped with cling film and installed with water supplied to the base (Figure 1). A thermocouple was inserted into the middle of each section, creating a 10 cm height interval. The remoulded samples were prepared using a similar technique to Sarsembayeva and Zhussupbekov [18], with the only difference being the height of the sample and the absence of chemical reagents. The sandy clay samples were modelled in grade and grain size to reproduce natural soils in Nur-Sultan city with similar plasticity limits and engineering properties. The initial physical characteristics of the soil samples are presented in Table 2.
Energies 2022, 15, x FOR PEER REVIEW 3 of 17 had a higher density and accordingly more mechanical energy, i.e., hammer blows, applied during packing. After packing, each column was securely wrapped with cling film and installed with water supplied to the base ( Figure 1). A thermocouple was inserted into the middle of each section, creating a 10 cm height interval. The remoulded samples were prepared using a similar technique to Sarsembayeva and Zhussupbekov [18], with the only difference being the height of the sample and the absence of chemical reagents. The sandy clay samples were modelled in grade and grain size to reproduce natural soils in Nur-Sultan city with similar plasticity limits and engineering properties. The initial physical characteristics of the soil samples are presented in Table 2.  Two freeze-thaw cycles were performed with a slow freeze technique to provide sufficient time for moisture transport and phase transfer and thereby ensure the most favourable conditions for frost heaving [18,19]. The freezing rate was 2 • C/day; the setting temperature dropped up to −24 • C in the chilling collars. However, the minimum temperature in the topsoil layer was only −16 • C to −18 • C. Temperature changes in the sample sections were recorded by 96 thermocouples and recorded via Pica loggers over the entire period of the freeze-thaw cycles. Additional 3.5 kPa surcharges were placed on the soil surfaces, representing the weight of the pavement and evenly distributing the heat extraction over the surface of the sample. Solid plastic rods screwed to surcharges served to measure the vertical linear movements during the freeze-thaw cycles.
After the second freezing cycle, all the samples were immediately removed and moisture was determined in each section by oven drying at 105 • C. The redistribution of dry density was determined by weighing each section and determining the moisture content in each section separately.

Water Mass Transfer in Freezing Soils
The method presented in this paper, although outwardly similar to the previous work of Konrad and Morgenstern, on the concept of the ice segregation and water migration based on hygroscopic water transport [5,6], is fundamentally different in the presence of the gaseous component in the transport of water in frozen soil. The further stages of frost heave have been considered for unsaturated soils without a chemical component: In Figure 2a, an unfrozen soil sample with relatively low moisture content is presented. The degree of saturation is S 0 < 1, therefore the sample is unsaturated. The soil is considered as a closed system with a phase equilibrium condition. The moisture here is considered in two phases: the liquid part presented with gravitational, hydroscopic and capillary water and the gas phase in the form of saturated vapour. The equilibrium between the water and vapour phases in a porous soil media is established when the temperature T and pressure P are stabilized. ii.
In Figure 2b, the sample is considered at the beginning of the freezing period, when the ice lenses' cores have just started to segregate. The soil structure in the freezing fringe includes the moisture in three phases: the solid part-ice lenses; the liquid phase-hydroscopic; and the capillary water and the gas phase-saturated and unsaturated vapour. iii.
In Figure 2c, the moisture mass transfer in terms of vapour migration continues, progressively filling all the pore volume with ice, which starts with the pores nearest to the vertical channel or the so-called ice veins. iv.
In Figure 2d, all pores of the top layer are filled with ice and the suction in the vertical vein is still inducing the vapour mass transfer; the uneven filling of the soil layer induces a horizontal crack. This crack continues to fill with ice and lifts the soil layer in the only possible direction-upwards. As a result, the formed frost heaving is equal to ∆ (Figure 2d). The thickness of the horizontal crack increases until the overburden pressure equilibrates with cryosuction, after which the vapour mass continues to fill the lower layer pores ( Figure 3). of the gaseous component in the transport of water in frozen soil. The further stages of frost heave have been considered for unsaturated soils without a chemical component: i.
In Figure 2a, an unfrozen soil sample with relatively low moisture content is presented. The degree of saturation is S₀ < 1, therefore the sample is unsaturated. The soil is considered as a closed system with a phase equilibrium condition. The moisture here is considered in two phases: the liquid part presented with gravitational, hydroscopic and capillary water and the gas phase in the form of saturated vapour. The equilibrium between the water and vapour phases in a porous soil media is established when the temperature T and pressure P are stabilized.
ii. In Figure 2b, the sample is considered at the beginning of the freezing period, when the ice lenses' cores have just started to segregate. The soil structure in the freezing fringe includes the moisture in three phases: the solid part-ice lenses; the liquid phase-hydroscopic; and the capillary water and the gas phase-saturated and unsaturated vapour. iii. In Figure 2c, the moisture mass transfer in terms of vapour migration continues, progressively filling all the pore volume with ice, which starts with the pores nearest to the vertical channel or the so-called ice veins. iv. In Figure 2d, all pores of the top layer are filled with ice and the suction in the vertical vein is still inducing the vapour mass transfer; the uneven filling of the soil layer induces a horizontal crack. This crack continues to fill with ice and lifts the soil layer in the only possible direction-upwards. As a result, the formed frost heaving is equal to ∆ (Figure 2d). The thickness of the horizontal crack increases until the overburden pressure equilibrates with cryosuction, after which the vapour mass continues to fill the lower layer pores ( Figure 3).

Vapour Mass Tranfer in the Gaseous State
The calculation of the moisture transfer in a gaseous state was based on the temperature data and a volumetric ratio of gaseous liquid and solid medium in timeline. Porosity and the volume of air are the crucial characteristics for vapour mass transfer in freezing soils. The reason for this is that the mass transfer of vapour is based on convection of a vapour flow from warmer to colder parts due to the difference in pore pressures during uniaxial freezing and accompanying phase transitions, i.e., condensation and freezing of excess gas on cold surfaces, as well as evaporation in warmer lower soil layers. A greater volume of air means more convection capacity in the soil. The amount of vapour in the pores is also highly temperature dependent. Evaporation/sublimation occurs in warm parts and condensation/deposition in the cold due to the pressure of saturated vapour over ice or water, therefore the transfer of moisture in the gaseous state explains the cryosuction forces. It is notable that the phase transfer from the gas to the liquid state happens with the high compression of molecules and volume release. The latent heat for condensation to 1 kg of water = 2.3 × 10 J is more than six times higher than that for segregation to ice = 335 × 10 J. Since there is no capillary surface tension during the vapour transport, it can pass through the smallest pores. An algorithm for determining mass transfer in a gaseous state is shown in Figure 4. Therefore, by knowing the volume and mass of the soil, and , as well as the humidity and temperatures T at the boundary points along the height, it is possible to determine the bulk and dry density of soil and the corresponding volumetric weight characteristics of each constituent: water , solid particles , air voids and porosity e.

Vapour Mass Tranfer in the Gaseous State
The calculation of the moisture transfer in a gaseous state was based on the temperature data and a volumetric ratio of gaseous liquid and solid medium in timeline. Porosity and the volume of air are the crucial characteristics for vapour mass transfer in freezing soils. The reason for this is that the mass transfer of vapour is based on convection of a vapour flow from warmer to colder parts due to the difference in pore pressures during uniaxial freezing and accompanying phase transitions, i.e., condensation and freezing of excess gas on cold surfaces, as well as evaporation in warmer lower soil layers. A greater volume of air means more convection capacity in the soil. The amount of vapour in the pores is also highly temperature dependent. Evaporation/sublimation occurs in warm parts and condensation/deposition in the cold due to the pressure of saturated vapour over ice or water, therefore the transfer of moisture in the gaseous state explains the cryosuction forces. It is notable that the phase transfer from the gas to the liquid state happens with the high compression of molecules and volume release. The latent heat for condensation to 1 kg of water C cond = 2.3 × 10 6 J is more than six times higher than that for segregation to ice C = 335 × 10 3 J. Since there is no capillary surface tension during the vapour transport, it can pass through the smallest pores.
An algorithm for determining mass transfer in a gaseous state is shown in Figure 4. Therefore, by knowing the volume and mass of the soil, V soil and m soil , as well as the humidity W and temperatures T at the boundary points along the height, it is possible to determine the bulk and dry density of soil ρ dry and the corresponding volumetric weight characteristics of each constituent: water V water , solid particles V solid , air voids V air and porosity e.
Another important parameter to determine the mass transfer is time t i . It is presumed that the concentration of vapour in an enclosed space tends to full saturation. To determine the mass of the vapour m vapour,t i at the initial or final moment of the period, it is necessary to know the pressure of saturated vapour over water or ice corresponding to the phase state and temperature: where µ-molar mass of vapour and equal to 18.01528 (33) g/mol; P sw , P si -saturated vapour pressure over water, ice surface corresponding to the temperature, Pa; R-Universal gas constant, 8.3144598 J/(mol·K), which corresponds to 8.3144598 m 3 ·Pa K·mol , T-temperature in Kelvin.
Further, by the mass of vapour, its density is determined for a certain time and temperature ( Figure 4). However, to find the volume of vapour passing through the soil section over time the cross-section of the air channel must be determined.  Another important parameter to determine the mass transfer is time . It is presumed that the concentration of vapour in an enclosed space tends to full saturation. To determine the mass of the vapour , at the initial or final moment of the period, it is necessary to know the pressure of saturated vapour over water or ice corresponding to the phase state and temperature: where μ-molar mass of vapour and equal to 18.01528 (33) g/mol; , -saturated vapour pressure over water, ice surface corresponding to the temperature, Pa; R-Universal gas constant, 8.3144598 J/(mol·K), which corresponds to 8.3144598 Further, by the mass of vapour, its density is determined for a certain time and temperature ( Figure 4). However, to find the volume of vapour passing through the soil section over time the cross-section of the air channel must be determined.
Based on experimental studies [1,12,16,18], porosity and volume of air are the most important characteristics for moisture transfer in freezing soils. The classical understanding of the volumetric relationship of a soil sample is given on the left side of Figure 5. In this study, the volumes are regrouped as separate cylinders ( Figure 5) as this makes it possible to calculate the cross-sectional area of each of them. The calculation of vapour volume passed through the section over time is found via the speed of vapour passing through the air voids' cross-section over the time t: Heat transfer in time period:  Based on experimental studies [1,12,16,18], porosity and volume of air are the most important characteristics for moisture transfer in freezing soils. The classical understanding of the volumetric relationship of a soil sample is given on the left side of Figure 5. In this study, the volumes are regrouped as separate cylinders (Figure 5) as this makes it possible to calculate the cross-sectional area of each of them. The calculation of vapour volume passed through the section over time is found via the speed of vapour passing through the air voids' cross-section A air over the time t: where v-average speed of vapour, cm/h; and A air -cumulative cross-section of the air voids A air = π·d air 2 4 , cm 2 , corresponding to the porosity coefficient and moisture content.
where -average speed of vapour, cm/h; and -cumulative cross-section of the air voids = • , cm², corresponding to the porosity coefficient and moisture content. Further, returning to the algorithm for calculating the vapour mass transfer in Figure  4, the heat transfer for the gaseous constitutive Q during the calculated period t of the test in section was found as total energy spent on the temperature drop ΔΤ of the gaseous phase Q₁ and the latent heat of the vapour part transitioning from the gas phase to the solid ice phase Q₂, called the deposition.
Expressing the mass with volume and density and substituting Equation (1) for the vapour volume it became possible to calculate the speed of vapour transfer. The vapour speed is found at the beginning and at the end of the time period: The mass of ice built from the vapour passing over time with speed over the cumulative air channel cross-section is calculated: where -mass of built ice in grams; -is taken as an average density value of the vapour densities at the start and the end time point, g/cm³.
The volume ratio between the phases of a liquid and a solid medium changes with time and, accordingly, so does the pore space. Consequently, it is necessary to consider the formed ice lenses and cracks, and the negative pressure present due to the formation of lenses during frost heaving. These volumetric changes need to be tracked over time and updated in computer simulations of ice lens growth.
The main benefit of the proposed calculation method is the determination of vapour mass transfer based on the known physical laws. The calculation considers the data obtained through the ongoing monitoring of volumetric and temperature data without applying empirical formulas or coefficients.

Results
Temperature distribution along the x axis can be conditionally divided into: first, the consolidation zone; second, the freezing and thawing cycles ( Figure 6). Temperature distribution over the sample length showed a cooling deceleration in the transition zone through zero degrees due to the latent heat of freezing. This is especially evident in large Further, returning to the algorithm for calculating the vapour mass transfer in Figure 4, the heat transfer for the gaseous constitutive Q during the calculated period t of the test in section was found as total energy spent on the temperature drop ∆T of the gaseous phase Q 1 and the latent heat of the vapour part transitioning from the gas phase to the solid ice phase Q 2 , called the deposition.
Expressing the mass with volume and density and substituting Equation (1) for the vapour volume it became possible to calculate the speed of vapour transfer. The vapour speed is found at the beginning and at the end of the time period: The mass of ice built from the vapour passing over time t with speed v over the cumulative air channel cross-section A air is calculated: where m ice -mass of built ice in grams; ρ vapour -is taken as an average density value of the vapour densities at the start and the end time point, g/cm 3 . The volume ratio between the phases of a liquid and a solid medium changes with time and, accordingly, so does the pore space. Consequently, it is necessary to consider the formed ice lenses and cracks, and the negative pressure present due to the formation of lenses during frost heaving. These volumetric changes need to be tracked over time and updated in computer simulations of ice lens growth.
The main benefit of the proposed calculation method is the determination of vapour mass transfer based on the known physical laws. The calculation considers the data obtained through the ongoing monitoring of volumetric and temperature data without applying empirical formulas or coefficients.

Results
Temperature distribution along the x axis can be conditionally divided into: first, the consolidation zone; second, the freezing and thawing cycles ( Figure 6). Temperature distribution over the sample length showed a cooling deceleration in the transition zone through zero degrees due to the latent heat of freezing. This is especially evident in large time intervals between 0 • C and −5 • C isotherms, which is explained by the specific heat of crystallization of water in the soil and requires much more energy than for its cooling. There is no clear pattern in the distribution of the temperature field inside samples with different densities. Slow freezing in this test gives enough time for phase transitions and uniform cooling and the density of the soil samples in this case does not affect the freezing rate. However, in loose density samples (Figure 6a-d), the isotherms in the first freezing cycle penetrate deeper than in the second, this is due to the large moisture content drawn upward in the second freezing cycle than in the initial state of the samples before testing. Due to the high heat capacity of water, a larger amount of heat removal is required to cool the samples.
Energies 2022, 15, x FOR PEER REVIEW 9 of 17 time intervals between 0 °C and −5 °C isotherms, which is explained by the specific heat of crystallization of water in the soil and requires much more energy than for its cooling.
There is no clear pattern in the distribution of the temperature field inside samples with different densities. Slow freezing in this test gives enough time for phase transitions and uniform cooling and the density of the soil samples in this case does not affect the freezing rate. However, in loose density samples (Figure 6a-d), the isotherms in the first freezing cycle penetrate deeper than in the second, this is due to the large moisture content drawn upward in the second freezing cycle than in the initial state of the samples before testing. Due to the high heat capacity of water, a larger amount of heat removal is required to cool the samples. Moisture redistribution after freeze-thaw cycles is presented in Figure 7. In the top 10 cm layer, ice lens formation caused a highly irregular distribution of moisture, which was confirmed by the centimetre-resolution sampling. Since the ground water supply table was quite shallow at 45 cm depth the moisture intake is high and this induces the capillary rise. For this reason, the moisture increase at a depth of 15 to 45 cm from the soil surface was very stable and depended just on the density of the soil samples.
The frost heave that occurred in the soil samples with various densities indicates a strong connection between the packing of soil particles and the amount of frost heaving and its progress (Figure 8). The maximum rates of frost heave were achieved in columns #7 and #8, where the dry density was close to the maximum value 1650-1790 kg/m³, while the loose soil samples with dry density 1180-1470 kg/m³ registered very weak heaving in the first cycle and consolidation or compression in the second cycle compared with the initial volume. Moisture redistribution after freeze-thaw cycles is presented in Figure 7. In the top 10 cm layer, ice lens formation caused a highly irregular distribution of moisture, which was confirmed by the centimetre-resolution sampling. Since the ground water supply table was quite shallow at 45 cm depth the moisture intake is high and this induces the capillary rise. For this reason, the moisture increase at a depth of 15 to 45 cm from the soil surface was very stable and depended just on the density of the soil samples.
The frost heave that occurred in the soil samples with various densities indicates a strong connection between the packing of soil particles and the amount of frost heaving and its progress (Figure 8). The maximum rates of frost heave were achieved in columns #7 and #8, where the dry density was close to the maximum value 1650-1790 kg/m 3 , while the loose soil samples with dry density 1180-1470 kg/m 3 registered very weak heaving in the first cycle and consolidation or compression in the second cycle compared with the initial volume.
The dry density change in soil samples during two freeze-thaw cycles are presented in Figure 9. The columns with loose densities were consolidated throughout the column length with the densest part in the bottom of the samples. In comparison, in the dense soil samples only the 10 cm layer on top was loosened, which also attracted a large amount of water to form ice lenses.  The dry density change in soil samples during two freeze-thaw cycles are presented in Figure 9. The columns with loose densities were consolidated throughout the column length with the densest part in the bottom of the samples. In comparison, in the dense soil samples only the 10 cm layer on top was loosened, which also attracted a large amount of water to form ice lenses.  The dry density change in soil samples during two freeze-thaw cycles are presented in Figure 9. The columns with loose densities were consolidated throughout the column length with the densest part in the bottom of the samples. In comparison, in the dense soil samples only the 10 cm layer on top was loosened, which also attracted a large amount of water to form ice lenses.

Vapour Mass Transfer in Soil Samples
Since the density, moisture and void ratio have the greatest influence on the amount of vapour transported, the general volumetric parameters of the tested soil columns are presented in Table 3. Table 3. Volumetric parameters of soil columns after two freeze-thaw cycles.

Column
Mass of Soil Based on the volumetric data, temperature change and the distribution of these parameters over the height of the samples, the vapour transfer was calculated using the example of columns #1, #7, #8 and #9 during the calculated period 494-518 h in the second freezing cycle (Tables 4 and A1-A3). Columns #7, #8 and #9 were chosen as the densest and preferred as soil bases for roads, buildings and structures. The calculation was carried out according to the algorithm presented in Figure 4. The void ratio was determined according to the vertical movement gauges' data from the top of the environmental chamber. According to the temperature data for 494 h and 518 h, the corresponding mass of saturated vapour and density were calculated for each 10 cm high section.

Vapour Mass Transfer in Soil Samples
Since the density, moisture and void ratio have the greatest influence on the amount of vapour transported, the general volumetric parameters of the tested soil columns are presented in Table 3. Based on the volumetric data, temperature change and the distribution of these parameters over the height of the samples, the vapour transfer was calculated using the example of columns #1, #7, #8 and #9 during the calculated period 494-518 h in the second freezing cycle (Tables 4 and A1, Tables A2 and A3). Columns #7, #8 and #9 were chosen as the densest and preferred as soil bases for roads, buildings and structures. The calculation was carried out according to the algorithm presented in Figure 4. The void ratio was determined according to the vertical movement gauges' data from the top of the environmental chamber. According to the temperature data for 494 h and 518 h, the corresponding mass of saturated vapour and density were calculated for each 10 cm high section. The heat released from the vapour constitutive comprises the energy extracted for cooling down from the temperature at 494 h to the temperature at 518 h and the latent heat of the vapour part transitioning from the gas phase to the solid ice phase, called the deposition. Since the latent heat of deposition of the freezing part of the vapour is more than 100 times higher than the temperature difference in the heat capacity of the vapour during cooling in 24 h, this part of the energy has a greater influence.
As the temperature decreased approximately for 2 • C per 24 h in the top layers, the pressure of saturated vapour over the ice reduced to 11-13%, which caused vapour migration towards lower temperatures. At the same time, the similar part of the vapour deposited and froze on the cold surfaces of the ice crystals from above. In the lower sections, the temperature drop comprised 0.5-1.0 • C per 24 h, while the saturated vapour pressure over ice reduced to 3-7%.
The void ratio in the top section is two times that at 50 cm depth from the top in the dense samples; the loosening process is clearly observed in the top layers of the soil samples packed with 1.64-1.76 g/cm 3 dry density. In the soil loose density samples with 1.21-1.63 g/cm 3 the void ratio redistribution is not so unambiguous. The frost heaving and consolidation processes in the loose samples led to the restructuring of the skeleton of solid particles in the sample. The volume of air/vapour in the pores of one section was varying from 4 cm 3 in the top to 19 cm 3 in the dense samples, while the top layers' pores were filled with ice lenses. The loose packed samples had 170-215 cm 3 air volume and the intensive vapour mass transfer occurred there due to the increased air volume, as the cumulative cross-section is larger there.
The speed of vapour transfer is related to the temperature field distribution. The vapour rate was around 0.4 m/h in the heaved soils, while in the phase transition zone it accelerated due to the extensive energy of latent heat for deposition.
The obtained ice mass formation rate was in the range of 1.64 × 10 −5 -3.60 × 10 −5 g/h in the loose samples and 1.42 × 10 −6 g/h to 5.61 × 10 −7 g/h in the dense samples (Table 5). In the top layers with excessive ice lens formation the vapour transfer speed dropped to 6.40 × 10 −8 g/h (Table A1).

Discussion
The modified freeze-thaw cycles with slow freezing rates and longer sample length allowed for approximating the freezing process of the soils close to natural conditions. According to Iushkov and Sergeev, the greatest values of the frost heave are achieved when the freezing rate in the unidirectional freezing test of the soil is compiled between 2 and 3 cm/day [19], while Arenson et al., noted that the steady state thermal conditions were achieved for the final ice lens growth at a cooling rate of 3 • C/day [20]. In the current study, the temperature decrease of the cold plate was set at 2 • C/day and was run independently from the freezing rate in the samples. The most favourable conditions for the frost heave formation were distinguished at the freezing rate of 2 cm/day, which agrees with the results obtained by Iushkov and Sergeev [19].
The conceptual model ice growth of Arenson et al., is correct for coarse grain soils with lower surfactant forces and a greater proportion of gravitational water [13]. However, it is not applicable to frost susceptible soils, where the surface tension and particle surfaces charge come to the fore, while the role of vapour transfer is undeservedly underestimated in the frost heaving processes.
Henry, while stating that the latent heat of freezing is part of the frost heave transport process, also proposed that there is a relationship between frost heave and the drying or evaporating from liquid water, yet went no further in defining this [12]. She emphasized how the initiation of ice lenses is associated with the phase transformation from liquid to solid state, accompanied by the replenishment of desiccated areas with water flow from the unfrozen soils. However, this understanding has been implied only regarding the liquid and solid states of the soil structure.
The amount of transported moisture is related to the withdrawn energy, which is distributed to the phase transfer energy and the cooling of each soil component, according to its heat conductivity. The supplementary requirement for the mass transfer is sufficient time. Incidentally, this is related to the travel distance of the vapour along the pore channels tortuosity.
Therefore, the freezing energy is equal to: where Q s , Q w , Q i , Q v -heat loss for cooling down the solids, water, ice and vapour components, respectively, and Q w−i -latent heat for the phase transfer from water to the ice state.
In Figure 10 the heat loss for the phase transfer from vapour to ice includes the heat required to transfer the vapour first from a gas to a liquid state, and subsequently into the solid state. The energy loss for the vapour to water phase transfer in the frozen fringe equilibrates with the water to vapour transition in the lower soil layers to stabilize the negative pressure. Eventually, the mass moisture transfer energy corresponds with the crystallization heat release from the liquid to solid state and the low energy-intensive cooling heat. Therefore, the vapour speed increases within the transition through zero degrees.
frost heave formation were distinguished at the freezing rate of 2 cm/day, which agrees with the results obtained by Iushkov and Sergeev [19].
The conceptual model ice growth of Arenson et al., is correct for coarse grain soils with lower surfactant forces and a greater proportion of gravitational water [13]. However, it is not applicable to frost susceptible soils, where the surface tension and particle surfaces charge come to the fore, while the role of vapour transfer is undeservedly underestimated in the frost heaving processes.
Henry, while stating that the latent heat of freezing is part of the frost heave transport process, also proposed that there is a relationship between frost heave and the drying or evaporating from liquid water, yet went no further in defining this [12]. She emphasized how the initiation of ice lenses is associated with the phase transformation from liquid to solid state, accompanied by the replenishment of desiccated areas with water flow from the unfrozen soils. However, this understanding has been implied only regarding the liquid and solid states of the soil structure.
The amount of transported moisture is related to the withdrawn energy, which is distributed to the phase transfer energy and the cooling of each soil component, according to its heat conductivity. The supplementary requirement for the mass transfer is sufficient time. Incidentally, this is related to the travel distance of the vapour along the pore channels tortuosity.
Therefore, the freezing energy is equal to: where , , , -heat loss for cooling down the solids, water, ice and vapour components, respectively, and -latent heat for the phase transfer from water to the ice state.
In Figure 10 the heat loss for the phase transfer from vapour to ice includes the heat required to transfer the vapour first from a gas to a liquid state, and subsequently into the solid state. The energy loss for the vapour to water phase transfer in the frozen fringe equilibrates with the water to vapour transition in the lower soil layers to stabilize the negative pressure. Eventually, the mass moisture transfer energy corresponds with the crystallization heat release from the liquid to solid state and the low energy-intensive cooling heat. Therefore, the vapour speed increases within the transition through zero degrees. The hypothesis formation of ice lenses presented in Figure 3 is perfectly combined with the theory of vapour transfer, when water in a gaseous state penetrates through any holes bypassing the surface tension of solid soil particles. The negative pressure occurring in the crack due to the volume expansion enhances the vapour transport to the crack.
Each 1 mm of crack formed in ice lens formation according to the presented examples caused an addition 78.5 cm³ of air volume corresponding to the existing 19 cm³ volume of gaseous state, which caused the reduction of saturated vapour pressure to −1115.5 Pa; not enough to draw in water, although this negative pressure is quite enough to advance the movement of vapour. The hypothesis formation of ice lenses presented in Figure 3 is perfectly combined with the theory of vapour transfer, when water in a gaseous state penetrates through any holes bypassing the surface tension of solid soil particles. The negative pressure occurring in the crack due to the volume expansion enhances the vapour transport to the crack.
Each 1 mm of crack formed in ice lens formation according to the presented examples caused an addition 78.5 cm 3 of air volume corresponding to the existing 19 cm 3 volume of gaseous state, which caused the reduction of saturated vapour pressure to −1115.5 Pa; not enough to draw in water, although this negative pressure is quite enough to advance the movement of vapour.
This research has implications for infrastructure development in areas that experience seasonal freezing. Ice formation cannot be entirely mitigated by simply draining soils to an unsaturated condition. Consideration of the potential impact of vapour-driven ice growth, and how this occurs, may reduce the risk of damage.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.