# Impact of Isothermal Layering on the Stability of Saline Soil Subgrade in Cold Regions

Authors: Yu Zhang; Xuerui Chen; Meisi Zou; Runze Tian; Yunlong Hou; Bingbing Han
Journal: Geofluids
 (2022)

---

## Abstract

In order to determine the influence of heat on the stability of the saline soil subgrade slope in cold regions. Firstly, four typical representative saline soils were selected as subgrade fillers, and geotextiles were added at the bottom of the subgrade as a study subject. A two-dimensional numerical model for the temperature field of subgrade soil was established based on COMSOL Multiphysics. Secondly, combined with the main strength parameters of subgrade soil at different temperatures, an isothermal stratification-strength parameter function calculation model was proposed. Then, the SLOPE/W module in GeoStudio is used to stratify the established subgrade model according to isotherms and calculate the stability safety factor of saline soil subgrades in cold regions under different months, slopes, and subgrade heights, respectively. Finally, the influence of submergence at the foot of the slope on the stability of the subgrade slope is discussed. The results show: After the cold season coming, with the decreasing of temperature, the stability of HC-1 and HC-2 saline soil subgrade slope increased to different extents when the ground temperature dropped to -10°C. When the local temperature dropped below -10°C, the stability has slight decline; CS-1 and CS-2 saline soil subgrade in the ground temperature dropped to 10°C (HC-1 and HC-2 denote two groups of saline soil with high chloride-salt content (≥15%) and low sulfate salt content (≤0.5%), whereas CS-1 and CS-2 denote two groups of saline soil with low chloride-salt (≤5%)), the slope stability was significantly improved, and the stability of CS-1 increases slightly with further decrease in temperature.


---

## Body

## 1. Introduction

Highway distress in China, such as salt swelling, muddying, wet (solution) subsidence, and corrosion, occurs frequently in saline soil areas, which seriously affects the performance of subgrade [1]. Research from McRobert and Foley [2] showed that soil at areas of high salinity and high groundwater levels could produce salt crystals with the change of external temperature. The result could be the destruction of the internal structure of the soil, which will inevitably lead to the decline of the overall strength and stability of the subgrade.In cold regions, different kinds of salt in saline soil could produce corresponding salt crystals or mirabilite at different ambient temperatures, which will change the volume and strength of the soil. Therefore, the internal structure of the soil will be changed accordingly and the strength of the soil will be destroyed. Swenne [3] pointed out that salt will crystallize and NaCl•2H2O will form out of the chloride-salt and unfrozen water when the temperature is lower than 0°C. Liao and Archer [4, 5] mentioned that Na2SO4 • 10H2O will be produced out of sulfate and unfrozen water when the temperature reaches 33° C. Different kinds of salt in the soil affects its freezing temperature. Bing and Ma [6] found that the freezing temperature will drop when salinity increases, and the temperature rises when water content increases. The freezing temperature of the soil varies in accordance with the type of salt contained. The effects of common anions in the soil on freezing temperature can be listed as Cl-, CO32-, and SO42-. The influence of common cations in the soil on the freezing temperature can be ranked in descending order of magnitude as K+, Na +, and Ca2+. Banin and Anderson [7] proposed an equation to calculate the freezing temperature through the soluble salt content in porous medium. Predecessors explored the relationship between different temperature and different salinity crystallization. However, the different salt crystals generated by the high salt content saline soil in the cold season under non-freezing conditions directly affect the soil strength. Therefore, it is meaningful to discuss the effects of temperature and soil strength parameters on the stability of saline soil subgrade.Stability theories of subgrade slopes in cold region have previously been done. Weeks [8] deduced the “ice blocking seepage” model, in which they discovered that the excess pore water pressure was formed, for the lower layer of molten soil prevented the outward discharge of pore water. Another research done by Hutchinson [9] used the total stress analysis method. He found that the reduction of undrained shear strength of soil, which was caused by the increase of the water content when the soil melted, could trigger the instability of the slope. In addition, the increase of water content could also cause the accumulation of fragment ice on the frozen frontal surface. McRoberts and Morgenstern [10] pointed out that the melting or consolidation process of frozen soil on slopes increases the super-pore water pressure inside the soil, which eventually leads to slope instability. Although research on the mechanism of slope instability has been done, reports on the slope stability of saline soil subgrade and the change of the strength of the soil at its different depths caused by the temperature field change are seldom seen.In summary, the value of soil strength in most of the studies on slope stability of road base usually takes soil type into consideration only, instead of combining the climate condition at the construction site and soil characteristics in cold region, because under the effect of heat, the change amplitude of soil strength parameter is large, which will cause the change of slope stability of road base. Due to the large variation amplitude of soil strength parameters under thermal action, the stability of subgrade slope changes. Therefore, it is particularly important to study the stability of saline soil subgrade slope with different salinity and different salt types under thermal action.


## 2. Project Background

The highway from Qarham to Golmud, which is located at the area where its altitude varies greatly, crosses subtropical and cold temperate zones. It is dry here and lack of rainfall with large annual evaporation. The place has four distinctive seasons. The temperature can be as low as -33°C in cold season and reaches 35°C in warm season. The temperature difference between day and night is great [10–12].According to the actual engineering geology, a meteorological station was established at K664+360 of the Qarham-Golmud highway to monitor the weather conditions and the surface and underground temperature in Qarham, where temperature probes were placed at the surface, 40 cm underground, 80 cm underground, 160 cm underground, and 320cm underground, respectively, and collected once or twice a day. Finally, the data were aggregated to obtain the temperature and ground temperature changes in different months of the year.As Figure1 shows, the monthly average temperature of this region in July was the highest, reaching 23°C. The monthly average temperature in June, July, and August was above 15°C; the lowest temperature in December was as low as -9.2°C; 11, 12, 1, 2,. The monthly average temperature in November, December, January, February, and March was below 0°C. The temperature difference between the highest temperature and the lowest temperature was 32.2°C. The temperature changes significantly in four seasons, and conforms to the law of sinusoidal function in general. The temperature difference between day and night in winter and summer is large. The land surface temperature is close to the air temperature here. There is a significant difference between the ground temperature at 40 cm below the surface and the air temperature. The deeper the depth, the greater the temperature difference. The ground temperature at 320 cm below the surface is quite different from the air temperature. The ground temperature generally conforms to the trend of a sinusoidal function. The deeper the depth, the smaller the temperature change, and the ground temperature changes slightly around 10°C at 320 cm below the surface.Figure 1
Average temperature and ground temperature from 2016 to 2020.The Qarham-Golmud highway passes Qarham Salt Lake, Garsu Station, Yushui River, with a total length of 80.052 kilometers. The expressway was constructed with bi-directional and four-lane. The design speed is 100 kilometers per hour, the subgrade is 26 meters wide, and the pavement structure is bituminous concrete. The saline soil along the highway covers a large area, whose salt composition and content have regional differences [13]. The high salinity soil sampling points along the Qarham-Golmud highway are shown in Figure 2.Figure 2
High salinity soil sampling points along Qarham-Golmud highway.
(a)
Visible salt crystals (white)(b)
Drainage along the highwayThe temperature in the cold and warm seasons in this area is quite different, and the saline soils in the test section of soluble salts are mainly chloride-saline soils, whose physical and mechanical properties are susceptible to temperature changes. In this case, a large amount of salt crystals are generated at the foot of the slope and the drainage ditch, which will change the internal structure of the soil, its strength could be reduced, and the stability of the subgrade slope will change accordingly. According to the chemical analysis of soluble salts towards salinity of the soil along the high by Zhang et al. [14], the soil along the Qarham-Golmud highway was mainly chloride-salinized soil with a small amount of sulfate salty soil. The salinity is mainly between 3% and 50%. Most of the soils here belong to over-salted soil and strong saline soil. According to the freezing temperature test of some saline soils in this area by Zhang et al. [15], three types of saline soils (HC-1, HC-2, and CS-1) selected from the article will not freeze at -20° C, so they will not freeze under the ground temperature in the cold season in this area.


## 3. Computational Model

COMSOL Multiphysics is based on the finite element method, which can be used to describe the temperature and moisture fields of saline soil subgrade using two partial differential equations [16].


### 3.1. Temperature Field Governing Equation

According to the basic theory of heat transfer, using the latent heat of phase transition as the internal heat source, the differential equation of heat conduction can be expressed by equation (1):
(1)ρCθ∂T∂t=∇⋅λθ∇T+LI⋅ρI∂θI∂t+Lc⋅ρc∂θc∂t.In the equation, the latent heat of phase change is used as the heat source, and the latent heat of ice is extended to the latent heat of ice and salt crystals. Salt crystals only consider salt crystals (NaCl⋅2H2O) and mirabilite (NaSO4⋅10H2O). ρ, ρI, and ρc are the density of the soil, the density of ice, and the density of salt crystalskg/m3;Cθ is the specific heat J/kg⋅°C)); T is the transient temperature of the soil; t is the time s; λθ is the thermal conductivity of the soil W/m⋅°C)); ∇ is the Hamiltonian; LI is the latent heat of phase transition of water; Lc is the latent heat of phase transition of salt; θ, θI, and θc are the volumetric water content, the pore ice volume content, and the salt crystal volume content. For parameter values, refer to reference [17].


### 3.2. Governing Equation of Water Field

The movement of unfrozen water in the subgrade soil follows Darcy’s law. Water migration in unsaturated soils can be represented by Richard’s equation [14] using equation (2):
(2)∂θu∂t+ρIρw⋅∂θI∂t=∇⋅Dθu∇θu+kgθu.In the equation,∂θu/∂t represents the change of unfrozen water content with time, ρI⋅∂θI/∂t represents the change of ice content with time, and ∇⋅Dθu∇θu+kgθu represents the migration of water. kg is the permeability coefficient of unsaturated soil in the direction of gravity acceleration; Dθu is the diffusivity of unfrozen water in frozen soil [18]. The calculation equation is shown in equation (3):
(3)Dθu=kθucθu⋅I.Among them,kθu is the permeability (m/s) of the unsaturated soil; cθu is the specific water capacity (1/m); I is the impedance factor, which indicates that the pore ice blocks the unfrozen water migration, I=1010θI.


### 3.3. Simultaneous Equations

Based on experimental data, Xu [19] gave an empirical expression for the unfrozen water content in frozen soil, as shown in equation (4):
(4)w0wu=TTfB,T<Tf.In the equation,w0 is the initial moisture content of the soil (%); wu is the moisture content of the unfrozen water (%) when the negative temperature is T; Tf is the freezing temperature of the soil (°C), which can be determined through experiments. According to the physical characteristics analysis of the soil sample, the soil sample is silty clay, B is a constant, approximately 0.56.The “solid-liquid ratio” [20] represents the ratio of the volume of pore ice to unfrozen water in frozen soil, and is written as. According to equation (5), there are:
(5)BI=θIθu=1.1TTfB−1.1T<Tf0T≥Tf.In the equation, the coefficient 1.1 is the ratio of the density of water to ice.From equation (5), we know that the solid-liquid ratio BI is a single value function of temperature. Therefore, the relational equation of pore ice, unfrozen water, and temperature in frozen soil can be expressed by equation (6) as:
(6)θI=BIT⋅θu


### 3.4. Calculation Parameters and Boundary Conditions

#### 3.4.1. Thermal Characteristics of Soil

The thermal properties of soil mainly include thermal conductivity and volume-specific heat capacity. By consulting the literature [20, 21], the specific values are shown in Table 1.Table 1
Soil thermal parameters.
Nameλf/(W·m-1·°C-1)λu/(W·m-1·°C-1)Cf/(J·m-3·°C-1)Cu/(J·m-3·°C-1)Lc/(J·m-3)Silty clay1.3511.1252.357 × 1061.879 × 10660.3 × 106


#### 3.4.2. Boundary Conditions

The boundary conditions of the subgrade temperature field include the first type of boundary conditions, the second type of boundary conditions, and the third type of boundary conditions. In the subgrade temperature field hydrothermal coupling model, the first type of boundary conditions is adopted for the pavement, shady slope, and sun slope, and the second boundary condition of constant geothermal flow is adopted for the lower boundary. The thermal boundary condition at the bottom is taken as the temperature gradient, 0.03°C/m [22–24].The first type of boundary conditions [14] used on pavements, shades, and sun slopes can be fitted by equation (11):
(7)T=T0+Gt+Asinωt+φIn the equation,T0 on the right side of the equation represents the average temperature, and A represents the amplitude, which is the maximum temperature fluctuation in a cycle. The road, shady slope, and sun slope correspond to different amplitudes. ω is to determine the period of temperature change, generally 365 × 24 h, and here the unit is calculated in seconds, 365×24×3600s. φ usually the physical meaning in the sine function is the primary phase, and the value is generally determined according to the moment of starting the simulation. Gt is a function of time that increases the temperature of the Tibetan Plateau over the next 50 years by 2.2~2.6° C [25, 26]. The maximum temperature rise in this model in the next 50 years is 2.6° C, and its expression can be expressed by equation (12):
(8)Gt=C×t.In the equation,C is 0.052/8760×3600 on the right side of the equation, and t represents time in seconds.According to the attached surface theory [15], the least square method is used for fitting, and the values of T0, A, and φ are obtained, as shown in equations (13) to (16).
(9)Roadsurface:Tt=10.71+16.9sin2π365×24×3600t+1.43+0.052t365×24×3600(10)Sunslope:Tt=8.85+16.32sin2π365×24×3600t+1.41+0.052t365×24×3600(11)Shadyslope:Tt=7.65+16.04sin2π365×24×3600t+1.42+0.052t365×24×3600(12)Naturalsurface:Tt=8.6+15.89sin2π365×24×3600t+1.42+0.052t365×24×3600


### 3.5. Geometric Model

The subjects of the study are subgrade and partial foundation. Combined with the actual engineering situation of the Qarham-Golmud highway test section, a two-dimensional axisymmetric subgrade model is established, with a road width of 26 m, a subgrade height of 4 m, and a slope of 1 : 1.5. The bottom of the subgrade is paved with geotextiles at a distance of 0.5 meters on the ground.Subgrade filling is compacted silty clay, and the foundation soil is silty clay (Figure3).Figure 3
Subgrade schematic.


### 3.6. Validation of Ground Temperature Simulation

Due to the low water content of the subgrade soil with scarce precipitation in the area, the groundwater source was blocked by the addition of geotextiles to the subgrade. Therefore, the problem of salt and water migration due to temperature gradients was not considered.Comparing the theoretical calculation results of the ground temperature below the natural ground and at the mid-line of the subgrade with the actual ground-bore ground temperature test data in the cold season (Figure4), it can be seen that:Figure 4
High salinity soil sampling points along Qarham-Golmud highway.
(a)
natural ground surface(b)
subgrade central line(a) Natural ground surface(b) Subgrade central line(1)
The simulation results of soil temperature at different depths are consistent with the field measurement results. In the shallow layer, the simulation result and the actual measurement are quite different; in the deep layer, the difference between the simulation result and the actual measurement is small. The main reason is that the external environment changes greatly, and the surface layer is greatly affected by the accumulation effect of the external environment. As the depth increases, this effect gradually decreases(2)
The measured results of the subgrade temperature in the four months of the cold season (December, January, February, and March) are in good agreement with the simulation results. From the subgrade depth of 0 cm to the subgrade depth of 60 cm, the measured values are different from the simulation results, but the differences are not significant. When the depth reaches 160 cm to 320 cm, the measured value is very close to the simulation result(3)
On the whole, the numerical simulation results are very close to the measured data of boreholes. The change trends are consistent, which verifies the validity of the model


## 3.1. Temperature Field Governing Equation

According to the basic theory of heat transfer, using the latent heat of phase transition as the internal heat source, the differential equation of heat conduction can be expressed by equation (1):
(1)ρCθ∂T∂t=∇⋅λθ∇T+LI⋅ρI∂θI∂t+Lc⋅ρc∂θc∂t.In the equation, the latent heat of phase change is used as the heat source, and the latent heat of ice is extended to the latent heat of ice and salt crystals. Salt crystals only consider salt crystals (NaCl⋅2H2O) and mirabilite (NaSO4⋅10H2O). ρ, ρI, and ρc are the density of the soil, the density of ice, and the density of salt crystalskg/m3;Cθ is the specific heat J/kg⋅°C)); T is the transient temperature of the soil; t is the time s; λθ is the thermal conductivity of the soil W/m⋅°C)); ∇ is the Hamiltonian; LI is the latent heat of phase transition of water; Lc is the latent heat of phase transition of salt; θ, θI, and θc are the volumetric water content, the pore ice volume content, and the salt crystal volume content. For parameter values, refer to reference [17].


## 3.2. Governing Equation of Water Field

The movement of unfrozen water in the subgrade soil follows Darcy’s law. Water migration in unsaturated soils can be represented by Richard’s equation [14] using equation (2):
(2)∂θu∂t+ρIρw⋅∂θI∂t=∇⋅Dθu∇θu+kgθu.In the equation,∂θu/∂t represents the change of unfrozen water content with time, ρI⋅∂θI/∂t represents the change of ice content with time, and ∇⋅Dθu∇θu+kgθu represents the migration of water. kg is the permeability coefficient of unsaturated soil in the direction of gravity acceleration; Dθu is the diffusivity of unfrozen water in frozen soil [18]. The calculation equation is shown in equation (3):
(3)Dθu=kθucθu⋅I.Among them,kθu is the permeability (m/s) of the unsaturated soil; cθu is the specific water capacity (1/m); I is the impedance factor, which indicates that the pore ice blocks the unfrozen water migration, I=1010θI.


## 3.3. Simultaneous Equations

Based on experimental data, Xu [19] gave an empirical expression for the unfrozen water content in frozen soil, as shown in equation (4):
(4)w0wu=TTfB,T<Tf.In the equation,w0 is the initial moisture content of the soil (%); wu is the moisture content of the unfrozen water (%) when the negative temperature is T; Tf is the freezing temperature of the soil (°C), which can be determined through experiments. According to the physical characteristics analysis of the soil sample, the soil sample is silty clay, B is a constant, approximately 0.56.The “solid-liquid ratio” [20] represents the ratio of the volume of pore ice to unfrozen water in frozen soil, and is written as. According to equation (5), there are:
(5)BI=θIθu=1.1TTfB−1.1T<Tf0T≥Tf.In the equation, the coefficient 1.1 is the ratio of the density of water to ice.From equation (5), we know that the solid-liquid ratio BI is a single value function of temperature. Therefore, the relational equation of pore ice, unfrozen water, and temperature in frozen soil can be expressed by equation (6) as:
(6)θI=BIT⋅θu


## 3.4. Calculation Parameters and Boundary Conditions

### 3.4.1. Thermal Characteristics of Soil

The thermal properties of soil mainly include thermal conductivity and volume-specific heat capacity. By consulting the literature [20, 21], the specific values are shown in Table 1.Table 1
Soil thermal parameters.
Nameλf/(W·m-1·°C-1)λu/(W·m-1·°C-1)Cf/(J·m-3·°C-1)Cu/(J·m-3·°C-1)Lc/(J·m-3)Silty clay1.3511.1252.357 × 1061.879 × 10660.3 × 106


### 3.4.2. Boundary Conditions

The boundary conditions of the subgrade temperature field include the first type of boundary conditions, the second type of boundary conditions, and the third type of boundary conditions. In the subgrade temperature field hydrothermal coupling model, the first type of boundary conditions is adopted for the pavement, shady slope, and sun slope, and the second boundary condition of constant geothermal flow is adopted for the lower boundary. The thermal boundary condition at the bottom is taken as the temperature gradient, 0.03°C/m [22–24].The first type of boundary conditions [14] used on pavements, shades, and sun slopes can be fitted by equation (11):
(7)T=T0+Gt+Asinωt+φIn the equation,T0 on the right side of the equation represents the average temperature, and A represents the amplitude, which is the maximum temperature fluctuation in a cycle. The road, shady slope, and sun slope correspond to different amplitudes. ω is to determine the period of temperature change, generally 365 × 24 h, and here the unit is calculated in seconds, 365×24×3600s. φ usually the physical meaning in the sine function is the primary phase, and the value is generally determined according to the moment of starting the simulation. Gt is a function of time that increases the temperature of the Tibetan Plateau over the next 50 years by 2.2~2.6° C [25, 26]. The maximum temperature rise in this model in the next 50 years is 2.6° C, and its expression can be expressed by equation (12):
(8)Gt=C×t.In the equation,C is 0.052/8760×3600 on the right side of the equation, and t represents time in seconds.According to the attached surface theory [15], the least square method is used for fitting, and the values of T0, A, and φ are obtained, as shown in equations (13) to (16).
(9)Roadsurface:Tt=10.71+16.9sin2π365×24×3600t+1.43+0.052t365×24×3600(10)Sunslope:Tt=8.85+16.32sin2π365×24×3600t+1.41+0.052t365×24×3600(11)Shadyslope:Tt=7.65+16.04sin2π365×24×3600t+1.42+0.052t365×24×3600(12)Naturalsurface:Tt=8.6+15.89sin2π365×24×3600t+1.42+0.052t365×24×3600


## 3.4.1. Thermal Characteristics of Soil

The thermal properties of soil mainly include thermal conductivity and volume-specific heat capacity. By consulting the literature [20, 21], the specific values are shown in Table 1.Table 1
Soil thermal parameters.
Nameλf/(W·m-1·°C-1)λu/(W·m-1·°C-1)Cf/(J·m-3·°C-1)Cu/(J·m-3·°C-1)Lc/(J·m-3)Silty clay1.3511.1252.357 × 1061.879 × 10660.3 × 106


## 3.4.2. Boundary Conditions

The boundary conditions of the subgrade temperature field include the first type of boundary conditions, the second type of boundary conditions, and the third type of boundary conditions. In the subgrade temperature field hydrothermal coupling model, the first type of boundary conditions is adopted for the pavement, shady slope, and sun slope, and the second boundary condition of constant geothermal flow is adopted for the lower boundary. The thermal boundary condition at the bottom is taken as the temperature gradient, 0.03°C/m [22–24].The first type of boundary conditions [14] used on pavements, shades, and sun slopes can be fitted by equation (11):
(7)T=T0+Gt+Asinωt+φIn the equation,T0 on the right side of the equation represents the average temperature, and A represents the amplitude, which is the maximum temperature fluctuation in a cycle. The road, shady slope, and sun slope correspond to different amplitudes. ω is to determine the period of temperature change, generally 365 × 24 h, and here the unit is calculated in seconds, 365×24×3600s. φ usually the physical meaning in the sine function is the primary phase, and the value is generally determined according to the moment of starting the simulation. Gt is a function of time that increases the temperature of the Tibetan Plateau over the next 50 years by 2.2~2.6° C [25, 26]. The maximum temperature rise in this model in the next 50 years is 2.6° C, and its expression can be expressed by equation (12):
(8)Gt=C×t.In the equation,C is 0.052/8760×3600 on the right side of the equation, and t represents time in seconds.According to the attached surface theory [15], the least square method is used for fitting, and the values of T0, A, and φ are obtained, as shown in equations (13) to (16).
(9)Roadsurface:Tt=10.71+16.9sin2π365×24×3600t+1.43+0.052t365×24×3600(10)Sunslope:Tt=8.85+16.32sin2π365×24×3600t+1.41+0.052t365×24×3600(11)Shadyslope:Tt=7.65+16.04sin2π365×24×3600t+1.42+0.052t365×24×3600(12)Naturalsurface:Tt=8.6+15.89sin2π365×24×3600t+1.42+0.052t365×24×3600


## 3.5. Geometric Model

The subjects of the study are subgrade and partial foundation. Combined with the actual engineering situation of the Qarham-Golmud highway test section, a two-dimensional axisymmetric subgrade model is established, with a road width of 26 m, a subgrade height of 4 m, and a slope of 1 : 1.5. The bottom of the subgrade is paved with geotextiles at a distance of 0.5 meters on the ground.Subgrade filling is compacted silty clay, and the foundation soil is silty clay (Figure3).Figure 3
Subgrade schematic.


## 3.6. Validation of Ground Temperature Simulation

Due to the low water content of the subgrade soil with scarce precipitation in the area, the groundwater source was blocked by the addition of geotextiles to the subgrade. Therefore, the problem of salt and water migration due to temperature gradients was not considered.Comparing the theoretical calculation results of the ground temperature below the natural ground and at the mid-line of the subgrade with the actual ground-bore ground temperature test data in the cold season (Figure4), it can be seen that:Figure 4
High salinity soil sampling points along Qarham-Golmud highway.
(a)
natural ground surface(b)
subgrade central line(a) Natural ground surface(b) Subgrade central line(1)
The simulation results of soil temperature at different depths are consistent with the field measurement results. In the shallow layer, the simulation result and the actual measurement are quite different; in the deep layer, the difference between the simulation result and the actual measurement is small. The main reason is that the external environment changes greatly, and the surface layer is greatly affected by the accumulation effect of the external environment. As the depth increases, this effect gradually decreases(2)
The measured results of the subgrade temperature in the four months of the cold season (December, January, February, and March) are in good agreement with the simulation results. From the subgrade depth of 0 cm to the subgrade depth of 60 cm, the measured values are different from the simulation results, but the differences are not significant. When the depth reaches 160 cm to 320 cm, the measured value is very close to the simulation result(3)
On the whole, the numerical simulation results are very close to the measured data of boreholes. The change trends are consistent, which verifies the validity of the model


## 4. Test and Value of Strength Parameters of Subgrade Soil

### 4.1. Physical and Chemical Test Parameters

In order to study the situation of saline soil along the line comprehensively and accurately, and to investigate and analyze the distribution of saline soil along the way, four soil samples with different salinity were selected at the following four mileage piles. Results from test by Zhang et al. [14] about physical properties are as shown in Table 2 and its salinity are as shown in Table 3. In the table, HC-1 and HC-2 represent samples of high-chlorine saline soils, and CS-1 and CS-2 represent samples of relatively high-sulfate chloride-saline soils.Table 2
Saline soil properties parameter table.
Sampling positionSoil sampleProportionLiquid limit (%)Plastic limit (%)Soil sample categoryOptimal moisture content (%)Maximum dry density (g/cm3)Saturation (%)K603+260HC-12.4226.317.4Low-plasticity clay6.61.9675.6K602+750HC-22.5828.518.2Low-plasticity clay8.71.8358.6K636+000CS-12.6538.725.4Low-plasticity clay16.91.7795.5K642+200CS-22.5542.027.3Silt17.81.6689.5Table 3
Saline soil salt content.
Soil sampleContent of selected salts (weight %)Total salt content (%)CO32-HCO3-Cl-SO42-Ca2+Mg2+TotalHC-10.030.1317.930.480.600.9820.1436.1HC-20.000.0314.730.240.601.3416.9525.4CS-10.000.013.371.490.700.245.8111.9CS-20.000.160.530.720.300.241.962.6


### 4.2. Freezing Temperature Experiment

The freezing temperature of the soil is affected by salinity. CS-1 and CS-2 soil samples with low total salt content were selected for the freezing temperature experiment. At the beginning of the experiment, a starting temperature was preset in order to reduce the data collection effort while achieving the experimental purpose. After the temperature was reached, the configured soil samples were successively put into the prepared bags, inserted into the temperature collection probes, sealed, and put into the cold liquid of the low-temperature thermostat. Then, the low-temperature thermostat was adjusted again, and the preset temperature was set to -10°C and -20°C, respectively. The computer was started, the software was opened, and the temperature collection interval was set to 10s, so that the temperature changes inside the sample could be collected on the computer. When the display temperature of the low-temperature thermostat reaches -10°C and -20°C, the acquisition is stopped and the instrument is turned off. Figure5 shows the cooling curve of the soil sample.Figure 5
Cooling time history curve of the CS-1 and CS-2.The freezing temperature of free water at standard atmospheric pressure is 0 °C, while the pore water in the soil body makes its freezing temperature lower than 0°C due to its interaction with the surface of soil particles [27]. The freezing temperature of the soil is the temperature of pore water when it freezes stably. It is verified according to the fitting equation for freezing temperature of saline soils given by Chen [21].
(13)−Tf=2.8752W−1.0509Cl−0.8427SO42−0.3575.Calculations based on the major negative ion content of the soil samples, brought into Equation 22, yielded a theoretical freezing temperature of -2.8°C for CS-2 and -23°C for CS-1. For CS-2 soil, the freezing temperature is -2.7°C, which is about 4% lower than the value predicted by the equation. However, the cooling curve of the CS-1 soil sample did not show any indication of freezing or latent heat release when the temperature reached -20°C, which indicates that the freezing temperature of CS-1 was at least -20°C lower, which is also consistent with the prediction of the equation. In fact, the field survey data of the road route showed that its salt content was generally greater than 10%, which implies that most saline soils along the route do not freeze in winter.


### 4.3. Strength Parameters

The relationship between the temperature measured by Zhang et al. [14] and the internal friction angle of cohesion according to the temperature interval is as follows:The piece wise function relationship of HC-1 cohesion as a function of temperature and internal friction angle as a function of temperature is shown in equation (14):
(14)C=58.2−1.86TT∈10,2050.1−1.05TT∈0,1050.1−0.81TT∈−10,069+1.08TT∈−20,−10,φ=37.1−0.03TT∈10,2036.9−0.01TT∈0,1036.9−0.02TT∈−10,037.6+0.05TT∈−20,−10.The piece wise function relationship of HC-2 cohesion as a function of temperature and internal friction angle as a function of temperature is shown in equation (15):
(15)C=57−1.38TT∈10,2050.4−0.72TT∈0,1050.4−1.28TT∈−10,078.1+1.49TT∈−20,−10,φ=34.1−0.02TT∈10,2033.9T∈0,1033.9−0.04TT∈−10,034.4+0.01TT∈−20,−10.The piecewise function relationship of CS-1 cohesion with temperature and internal friction angle with temperature (16) is shown as:
(16)C=69.6−1.747TT∈10,2063−1.08TT∈0,1063−0.9TT∈−10,071.4−0.06TT∈−20,−10,φ=34.9−0.03TT∈10,2034.7−0.01TT∈0,1034.7−0.09TT∈−10,034.6−0.1TT∈−20,−10.The piece wise function relationship of CS-2 cohesion with temperature and internal friction angle with temperature (17) is shown as:
(17)C=90.6−1.32TT∈10,2084−0.66TT∈0,1084−2.16TT∈−10,0100.8−0.48TT∈−20,−10,φ=38.3−0.06TT∈10,2037.8−0.01TT∈0,1037.8−0.13TT∈−10,039−0.01TT∈−20,−10.Using Zhang et al. [12] to select natural foundation soil samples, direct shear experiments at different temperatures have obtained the relationship between the cohesion and internal friction angle of the foundation soil samples at different temperatures as shown in (18):
(18)C=14−0.33TT∈10,2012.1−0.14TT∈0,1012.1−0.75TT∈−10,027.7+0.81TT∈−20,−10,φ=21.6−0.33TT∈10,2019.4−0.11TT∈0,1019.4−0.13TT∈−10,021.3+0.06TT∈−20,−10.


### 4.4. Values Based on Temperature Field Strength Parameters

When the temperature of the subgrade soil at a certain depth is determined, its strength parameter will be determined by the function of temperature and strength parameter. Applying the method above to theoretical calculations, as shown in Figure6, it is assumed that the sliding surface is an arc sliding surface, i represents the ordinal number of soil strips, j represents the isotherm ordinal number passing through the bottom edge of the soil strip, Cij represents the cohesion of the soil strip i when the j-th isotherm passes through the bottom edge of the soil strip i. Similarly, φij is the magnitude of the internal friction angle of soil strip i when the j-th isotherm passes through the bottom edge of soil strip i. The value of Cij and φij can be determined by equations (14)~(18) by interpolation.Figure 6
Determine cohesion based on isotherms.


## 4.1. Physical and Chemical Test Parameters

In order to study the situation of saline soil along the line comprehensively and accurately, and to investigate and analyze the distribution of saline soil along the way, four soil samples with different salinity were selected at the following four mileage piles. Results from test by Zhang et al. [14] about physical properties are as shown in Table 2 and its salinity are as shown in Table 3. In the table, HC-1 and HC-2 represent samples of high-chlorine saline soils, and CS-1 and CS-2 represent samples of relatively high-sulfate chloride-saline soils.Table 2
Saline soil properties parameter table.
Sampling positionSoil sampleProportionLiquid limit (%)Plastic limit (%)Soil sample categoryOptimal moisture content (%)Maximum dry density (g/cm3)Saturation (%)K603+260HC-12.4226.317.4Low-plasticity clay6.61.9675.6K602+750HC-22.5828.518.2Low-plasticity clay8.71.8358.6K636+000CS-12.6538.725.4Low-plasticity clay16.91.7795.5K642+200CS-22.5542.027.3Silt17.81.6689.5Table 3
Saline soil salt content.
Soil sampleContent of selected salts (weight %)Total salt content (%)CO32-HCO3-Cl-SO42-Ca2+Mg2+TotalHC-10.030.1317.930.480.600.9820.1436.1HC-20.000.0314.730.240.601.3416.9525.4CS-10.000.013.371.490.700.245.8111.9CS-20.000.160.530.720.300.241.962.6


## 4.2. Freezing Temperature Experiment

The freezing temperature of the soil is affected by salinity. CS-1 and CS-2 soil samples with low total salt content were selected for the freezing temperature experiment. At the beginning of the experiment, a starting temperature was preset in order to reduce the data collection effort while achieving the experimental purpose. After the temperature was reached, the configured soil samples were successively put into the prepared bags, inserted into the temperature collection probes, sealed, and put into the cold liquid of the low-temperature thermostat. Then, the low-temperature thermostat was adjusted again, and the preset temperature was set to -10°C and -20°C, respectively. The computer was started, the software was opened, and the temperature collection interval was set to 10s, so that the temperature changes inside the sample could be collected on the computer. When the display temperature of the low-temperature thermostat reaches -10°C and -20°C, the acquisition is stopped and the instrument is turned off. Figure5 shows the cooling curve of the soil sample.Figure 5
Cooling time history curve of the CS-1 and CS-2.The freezing temperature of free water at standard atmospheric pressure is 0 °C, while the pore water in the soil body makes its freezing temperature lower than 0°C due to its interaction with the surface of soil particles [27]. The freezing temperature of the soil is the temperature of pore water when it freezes stably. It is verified according to the fitting equation for freezing temperature of saline soils given by Chen [21].
(13)−Tf=2.8752W−1.0509Cl−0.8427SO42−0.3575.Calculations based on the major negative ion content of the soil samples, brought into Equation 22, yielded a theoretical freezing temperature of -2.8°C for CS-2 and -23°C for CS-1. For CS-2 soil, the freezing temperature is -2.7°C, which is about 4% lower than the value predicted by the equation. However, the cooling curve of the CS-1 soil sample did not show any indication of freezing or latent heat release when the temperature reached -20°C, which indicates that the freezing temperature of CS-1 was at least -20°C lower, which is also consistent with the prediction of the equation. In fact, the field survey data of the road route showed that its salt content was generally greater than 10%, which implies that most saline soils along the route do not freeze in winter.


## 4.3. Strength Parameters

The relationship between the temperature measured by Zhang et al. [14] and the internal friction angle of cohesion according to the temperature interval is as follows:The piece wise function relationship of HC-1 cohesion as a function of temperature and internal friction angle as a function of temperature is shown in equation (14):
(14)C=58.2−1.86TT∈10,2050.1−1.05TT∈0,1050.1−0.81TT∈−10,069+1.08TT∈−20,−10,φ=37.1−0.03TT∈10,2036.9−0.01TT∈0,1036.9−0.02TT∈−10,037.6+0.05TT∈−20,−10.The piece wise function relationship of HC-2 cohesion as a function of temperature and internal friction angle as a function of temperature is shown in equation (15):
(15)C=57−1.38TT∈10,2050.4−0.72TT∈0,1050.4−1.28TT∈−10,078.1+1.49TT∈−20,−10,φ=34.1−0.02TT∈10,2033.9T∈0,1033.9−0.04TT∈−10,034.4+0.01TT∈−20,−10.The piecewise function relationship of CS-1 cohesion with temperature and internal friction angle with temperature (16) is shown as:
(16)C=69.6−1.747TT∈10,2063−1.08TT∈0,1063−0.9TT∈−10,071.4−0.06TT∈−20,−10,φ=34.9−0.03TT∈10,2034.7−0.01TT∈0,1034.7−0.09TT∈−10,034.6−0.1TT∈−20,−10.The piece wise function relationship of CS-2 cohesion with temperature and internal friction angle with temperature (17) is shown as:
(17)C=90.6−1.32TT∈10,2084−0.66TT∈0,1084−2.16TT∈−10,0100.8−0.48TT∈−20,−10,φ=38.3−0.06TT∈10,2037.8−0.01TT∈0,1037.8−0.13TT∈−10,039−0.01TT∈−20,−10.Using Zhang et al. [12] to select natural foundation soil samples, direct shear experiments at different temperatures have obtained the relationship between the cohesion and internal friction angle of the foundation soil samples at different temperatures as shown in (18):
(18)C=14−0.33TT∈10,2012.1−0.14TT∈0,1012.1−0.75TT∈−10,027.7+0.81TT∈−20,−10,φ=21.6−0.33TT∈10,2019.4−0.11TT∈0,1019.4−0.13TT∈−10,021.3+0.06TT∈−20,−10.


## 4.4. Values Based on Temperature Field Strength Parameters

When the temperature of the subgrade soil at a certain depth is determined, its strength parameter will be determined by the function of temperature and strength parameter. Applying the method above to theoretical calculations, as shown in Figure6, it is assumed that the sliding surface is an arc sliding surface, i represents the ordinal number of soil strips, j represents the isotherm ordinal number passing through the bottom edge of the soil strip, Cij represents the cohesion of the soil strip i when the j-th isotherm passes through the bottom edge of the soil strip i. Similarly, φij is the magnitude of the internal friction angle of soil strip i when the j-th isotherm passes through the bottom edge of soil strip i. The value of Cij and φij can be determined by equations (14)~(18) by interpolation.Figure 6
Determine cohesion based on isotherms.


## 5. Calculation and Analysis of Slope Stability

### 5.1. Calculation of Safety Factor

#### 5.1.1. Calculation Model of Isothermal Stratification—Strength Parameters

Combined with the function of the soil strength parameters and temperature measured in the laboratory, the isothermal layering-strength parameter calculation model is to layer the subgrade laterally in accordance with the soil isotherm of the subgrade, so as to calculating the safety factor of the subgrade slope and determining its slope stability. The value of strength to a certain depth is based on the temperature. The dashed line in Figure7 indicates the isotherm. The strength parameters of each layer of soil were given according to the type of subgrade soil in the previous research. In the isothermal layering-strength parameter model, however, the subgrade is subdivided on the basis of previous achievements in accordance with the isotherm. It takes both the strength parameters of different soils and the changes in outside air temperature that affect the strength parameters of the soil into account. Obviously, the changes of soil strength parameters can be more accurate.Figure 7
Schematic diagram of isothermal stratification-strength parameter model calculation.


#### 5.1.2. Solution Steps of SLOPE/W Module in GeoStudio

(1)
Model rendering(2)
Material area generating. Divide the soil layer into different areas according to the isotherm and draw them separately. Number of isothermal lines = material area+1(3)
Material properties defining. Select the Mohr-Coulomb criterion and enter parameters such as cohesion, internal friction angle, bulk density, and elastic modulus(4)
Material area assignment. Add the material properties defined in step ③ to the corresponding material area defined in step ②(5)
Draw the water pressure line(6)
Draw the entrance and exit ranges of the sliding surface(7)
Model inspection. After the parameters of all models are defined, check whether the models are defined correctly and ensure that the subsequent calculations are performed normally(8)
Calculation and solutionThe model calculation diagram is shown in Figure8.Figure 8
The calculation of the safety factor for January in the subgrade (4 m in height and 1 : 1.5 in slope) built on HC-1 soil.


### 5.2. Stability Analysis of Subgrade Slope in Cold Season

#### 5.2.1. Stability Analysis of Subgrade Slope under Heat

(1) Monthly Changes in Slope Stability. The subgrade fills are HC-1, HC-2, CS-1, and CS-2. When the subgrade width is 26 m, the subgrade height is 4 m, and the slope is 1: 1.5, the subgrade slope safety factor is shown in Figure 9.
(a)
It can be seen from Figure9 that, in general, the slope safety factors of the four groups of subgrades filled with soil samples change during these months, and its trend conforms to the variation law of sinusoidal function curve. In February, the slope safety coefficients reached their maximum values, and in August, it reached their minimum values, but the minimum values were greater than 2. In each month of the year, the subgrade slopes were stable(b)
In Figure8, the safety coefficients of subgrade slopes in descending order throughout the year are CS-2, CS-1, HC-2, and HC-1. HC-1 soil-filled subgrade slope safety factor reached a maximum of about 3.1 in February and August reached a minimum value of about 2.0, which decreased by 35.5% in August compared with February; HC-2 soil-filled subgrade slope safety factor reached a maximum of about 3.2 in February, a minimum of about 2.1 in August, and a decrease of 34.4% in August compared with February; the safety factor of CS-1 soil subgrade slope reached a maximum of about 3.4 in February and a minimum of August. The value was about 2.2, which decreased by 35.3% in August compared with February; the safety factor of CS-2 soil-filled subgrade slope reached a maximum of about 3.6 in February, reached a minimum of about 2.6 in August, and decreased by 27.8 in August from February %. Among them, HC-1, HC-2, and CS-1 have similar declines, while CS-2 has a smaller decline(c)
Above all, the stability of the subgrade slope changes during the months in cold season. The stability in February is the best, and in August is relatively poor, but the minimum value of the safety factor is greater than 2. Comparing the information of air temperature and ground temperature, it can be concluded that the temperature below the surface reaches the minimum value in February, but reaches the maximum value in October. The temperature has a certain influence on the safety factor of the slope of saline soil subgrade, and the stability of the slope of saline soil subgrade is higher in the cold season than in the summer, but the slope of the subgrade is stable throughout the year, but the subgrade slope is stable throughout the year. Under the same construction conditions, the stability of subgrade slopes filled with saline soil with low chloride ion content is higher than that of saline soil with high chloride ion contentFigure 9
Monthly safety factor of subgrade slope.(2) Influence of Subgrade Height on Subgrade Slope Stability. Subgrade fillings are HC-1, HC-2, CS-1, and CS-2; the width of pavement is 26 m; slope of 1: 1.5; subgrade height is 3 m, 3.5 m, 4 m, 4.5 m, and 5 m, respectively; the safety factors of the subgrade slope in December, January, February, and March are shown in (a), (b), (c), and (d) in Figure 10.Figure 10
Different subgrade filling, slope 1 : 1.5, slope safety coefficient of different subgrade height in cold season.
(a)
Subgrade filling HC-1(b)
Subgrade filling HC-2(c)
Subgrade filling CS-1(d)
Subgrade filling CS-2(a) Subgrade filling HC-1(b) Subgrade filling HC-2(c) Subgrade filling CS-1(d) Subgrade filling CS-2(a)
It can be seen from Figure10 that in the cold season, when the height of the subgrade increases from 3 m to 3.5 m, 4 m, 4.5 m, and 5 m, the safety factor gradually decreases(b)
As can be seen from Figure10(a), when the height of the HC-1 soil-filled subgrade increased from 3 m to 3.5 m, 4 m, 4.5 m, and 5 m, the safety factors in December were the lowest of the four months in the cold season, which dropped by 29.4% from about 3.4 to around 2.4; the safety factor in February was the highest of the four months in the cold season, which increased by 28.2% from about 3.9 to about 2.8(c)
As can be seen from Figure10(b), when the height of the HC-2 soil-filled subgrade increased from 3 m to 3.5 m, 4 m, 4.5 m, and 5 m, the safety factor in December was the lowest of the four months in the cold season, which decrease by 29.4% from about 3.4 reduced to about 2.4; the safety factor in February was the highest of the four months in the cold season, which decrease by 30% from about 4.0 to about 2.8(d)
As can be seen from Figure10(c), when the height of the CS-1 soil-filled subgrade increased from 3 m to 3.5 m, 4 m, 4.5 m, and 5 m, the safety factor in December was the lowest of the four months in the cold season; it was reduced from about 3.6 to about 2.6, a decrease of 27.8%; the safety factor in February was the highest in the four months of the cold season, from about 4.2 to about 3, a decrease of 28.6%.(e)
As can be seen from Figure10(d), when the height of the subgrade filled with CS-2 soil is increased from 3 m to 3.5 m, 4 m, 4.5 m, and 5 m, the safety factor in December is the lowest of the four months in the cold season, from about 3.8 to about 3.0, a decrease of 21.1%; the safety factor in February was the highest of the four months in the cold season, from about 4.5 to about 3.3, a decrease of 26.7%.(f)
As the height of the subgrade increases (3~5 m), the safety factor decreases. In December, HC-1 and HC-2 soil-filled subgrade slopes experienced the fastest decline in safety factor, and in February, the HC-2 soil-filled subgrade slopes experienced the fastest decline in safety factor. The stability of the slope also decreases with the increase of the height of the subgrade, but the minimum value of all safety factors is greater than 2. When the height of the subgrade is 3~5 m, the slope is stable(3) Influence of Subgrade Slope on Its Slope Stability. The subgrade fillings are HC-1, HC-2, CS-1, and CS-2, the width of the road surface is 26 m, and the height of the subgrade is 4 m. When the slopes are selected to be 1 : 1.5, 1 : 1.75, and 1 : 2, the safety factors of the subgrade slopes in December, January, and March are shown as (a), (b), (c), and (d) in Figure 11.Figure 11
Different subgrade filling, slope safety factors with different gradients in the cold season.
(a)
Subgrade filling HC-1(b)
Subgrade filling HC-2(c)
Subgrade filling CS-1(d)
Subgrade filling CS-2(a) Subgrade filling HC-1(b) Subgrade filling HC-2(c) Subgrade filling CS-1(d) Subgrade filling CS-2(a)
As can be seen from Figure11, in the cold season, if the slope reduced from 1 : 5 to 1 : 1.75 and 1 : 2, the slope safety factor generally increases(b)
As can be seen from Figure11(a), if the slope of the subgrade filled with HC-1 soil decreased from 1 : 5 to 1 : 1.75 and 1 : 2, the slope safety factor in February is higher than the other months in the cold season, and the minimum value It is about 3.1, the maximum value is about 3.3, which is an increase of 6%; the slope safety factor in December is lower than other months in the cold season, the minimum value is about 2.7, the maximum value is about 2.9, and an increase of 7.4%(c)
As can be seen from Figure11(b), when the slope of the subgrade filled with HC-2 soil reduced from 1 : 5 to 1 : 1.75 and 1 : 2, the slope safety factor in February was higher than that in other months in the cold season, and the minimum value is about 3.1, the maximum value is about 3.3, which increased by 6%; the slope safety factor in December is lower than other months in the cold season, the minimum value is about 2.7, the maximum value is about 3.1, and an increase of 14.8%(d)
As can be seen from Figure11(c), when the slope of the subgrade filled with CS-1 soil slowed down from 1 : 5 to 1 : 1.75 and 1 : 2, the slope safety factor in February was higher than that in other months in the cold season, and the minimum value is about 3.4, the maximum value is about 3.6, which increased by 5.9%; the slope safety factor in December is lower than other months in the cold season, with the minimum value of about 3.0, the maximum value of about 3.2, and an increase of 6.7%(e)
As can be seen from Figure11(d), when the slope of the subgrade filled with CS-2 soil slowed down from 1 : 5 to 1 : 1.75 and 1 : 2, the slope safety factor in February was higher than that in other months in the cold season, and the minimum value is about 3.6, the maximum value is about 3.7, an increase of 2.8%; the slope safety factor in December is lower than other months in the cold season, the minimum value is about 3.2, the maximum value is about 3.3, and an increase of 3.1%(f)
As the slope of the subgrade slows down (1 : 1.5 to 1 : 2), the safety factor increases. The fastest increase in the safety factor in February was the subgrade slopes filled with HC-1 and HC-2 soil, and the fastest increase in the safety factor in December was in the subgrade slopes filled with HC-2 soil. The stability of the slope also increases as the slope of the subgrade slows down. All safety factors are greater than 2.7. When the slope of the subgrade is 1 : 1.5~1 : 2, its slope is stable(g)
Measures to improve the stability of subgrade slopes. During the construction process, the salt content of the saline soil subgrade fill is strictly controlled; insulate the slope of the subgrade to reduce the amplitude of the change of temperature gradient and prevent the influence of large temperature difference on the stability of the slope of the subgrade; geotextiles or geomembranes are placed between the different layers of the roadbed slope to block the migration of water and salts and prevent the accumulation of salts that can lead to the reduction of the stability of the subgrade slope


## 5.1. Calculation of Safety Factor

### 5.1.1. Calculation Model of Isothermal Stratification—Strength Parameters

Combined with the function of the soil strength parameters and temperature measured in the laboratory, the isothermal layering-strength parameter calculation model is to layer the subgrade laterally in accordance with the soil isotherm of the subgrade, so as to calculating the safety factor of the subgrade slope and determining its slope stability. The value of strength to a certain depth is based on the temperature. The dashed line in Figure7 indicates the isotherm. The strength parameters of each layer of soil were given according to the type of subgrade soil in the previous research. In the isothermal layering-strength parameter model, however, the subgrade is subdivided on the basis of previous achievements in accordance with the isotherm. It takes both the strength parameters of different soils and the changes in outside air temperature that affect the strength parameters of the soil into account. Obviously, the changes of soil strength parameters can be more accurate.Figure 7
Schematic diagram of isothermal stratification-strength parameter model calculation.


### 5.1.2. Solution Steps of SLOPE/W Module in GeoStudio

(1)
Model rendering(2)
Material area generating. Divide the soil layer into different areas according to the isotherm and draw them separately. Number of isothermal lines = material area+1(3)
Material properties defining. Select the Mohr-Coulomb criterion and enter parameters such as cohesion, internal friction angle, bulk density, and elastic modulus(4)
Material area assignment. Add the material properties defined in step ③ to the corresponding material area defined in step ②(5)
Draw the water pressure line(6)
Draw the entrance and exit ranges of the sliding surface(7)
Model inspection. After the parameters of all models are defined, check whether the models are defined correctly and ensure that the subsequent calculations are performed normally(8)
Calculation and solutionThe model calculation diagram is shown in Figure8.Figure 8
The calculation of the safety factor for January in the subgrade (4 m in height and 1 : 1.5 in slope) built on HC-1 soil.


## 5.1.1. Calculation Model of Isothermal Stratification—Strength Parameters

Combined with the function of the soil strength parameters and temperature measured in the laboratory, the isothermal layering-strength parameter calculation model is to layer the subgrade laterally in accordance with the soil isotherm of the subgrade, so as to calculating the safety factor of the subgrade slope and determining its slope stability. The value of strength to a certain depth is based on the temperature. The dashed line in Figure7 indicates the isotherm. The strength parameters of each layer of soil were given according to the type of subgrade soil in the previous research. In the isothermal layering-strength parameter model, however, the subgrade is subdivided on the basis of previous achievements in accordance with the isotherm. It takes both the strength parameters of different soils and the changes in outside air temperature that affect the strength parameters of the soil into account. Obviously, the changes of soil strength parameters can be more accurate.Figure 7
Schematic diagram of isothermal stratification-strength parameter model calculation.


## 5.1.2. Solution Steps of SLOPE/W Module in GeoStudio

(1)
Model rendering(2)
Material area generating. Divide the soil layer into different areas according to the isotherm and draw them separately. Number of isothermal lines = material area+1(3)
Material properties defining. Select the Mohr-Coulomb criterion and enter parameters such as cohesion, internal friction angle, bulk density, and elastic modulus(4)
Material area assignment. Add the material properties defined in step ③ to the corresponding material area defined in step ②(5)
Draw the water pressure line(6)
Draw the entrance and exit ranges of the sliding surface(7)
Model inspection. After the parameters of all models are defined, check whether the models are defined correctly and ensure that the subsequent calculations are performed normally(8)
Calculation and solutionThe model calculation diagram is shown in Figure8.Figure 8
The calculation of the safety factor for January in the subgrade (4 m in height and 1 : 1.5 in slope) built on HC-1 soil.


## 5.2. Stability Analysis of Subgrade Slope in Cold Season

### 5.2.1. Stability Analysis of Subgrade Slope under Heat

(1) Monthly Changes in Slope Stability. The subgrade fills are HC-1, HC-2, CS-1, and CS-2. When the subgrade width is 26 m, the subgrade height is 4 m, and the slope is 1: 1.5, the subgrade slope safety factor is shown in Figure 9.
(a)
It can be seen from Figure9 that, in general, the slope safety factors of the four groups of subgrades filled with soil samples change during these months, and its trend conforms to the variation law of sinusoidal function curve. In February, the slope safety coefficients reached their maximum values, and in August, it reached their minimum values, but the minimum values were greater than 2. In each month of the year, the subgrade slopes were stable(b)
In Figure8, the safety coefficients of subgrade slopes in descending order throughout the year are CS-2, CS-1, HC-2, and HC-1. HC-1 soil-filled subgrade slope safety factor reached a maximum of about 3.1 in February and August reached a minimum value of about 2.0, which decreased by 35.5% in August compared with February; HC-2 soil-filled subgrade slope safety factor reached a maximum of about 3.2 in February, a minimum of about 2.1 in August, and a decrease of 34.4% in August compared with February; the safety factor of CS-1 soil subgrade slope reached a maximum of about 3.4 in February and a minimum of August. The value was about 2.2, which decreased by 35.3% in August compared with February; the safety factor of CS-2 soil-filled subgrade slope reached a maximum of about 3.6 in February, reached a minimum of about 2.6 in August, and decreased by 27.8 in August from February %. Among them, HC-1, HC-2, and CS-1 have similar declines, while CS-2 has a smaller decline(c)
Above all, the stability of the subgrade slope changes during the months in cold season. The stability in February is the best, and in August is relatively poor, but the minimum value of the safety factor is greater than 2. Comparing the information of air temperature and ground temperature, it can be concluded that the temperature below the surface reaches the minimum value in February, but reaches the maximum value in October. The temperature has a certain influence on the safety factor of the slope of saline soil subgrade, and the stability of the slope of saline soil subgrade is higher in the cold season than in the summer, but the slope of the subgrade is stable throughout the year, but the subgrade slope is stable throughout the year. Under the same construction conditions, the stability of subgrade slopes filled with saline soil with low chloride ion content is higher than that of saline soil with high chloride ion contentFigure 9
Monthly safety factor of subgrade slope.(2) Influence of Subgrade Height on Subgrade Slope Stability. Subgrade fillings are HC-1, HC-2, CS-1, and CS-2; the width of pavement is 26 m; slope of 1: 1.5; subgrade height is 3 m, 3.5 m, 4 m, 4.5 m, and 5 m, respectively; the safety factors of the subgrade slope in December, January, February, and March are shown in (a), (b), (c), and (d) in Figure 10.Figure 10
Different subgrade filling, slope 1 : 1.5, slope safety coefficient of different subgrade height in cold season.
(a)
Subgrade filling HC-1(b)
Subgrade filling HC-2(c)
Subgrade filling CS-1(d)
Subgrade filling CS-2(a) Subgrade filling HC-1(b) Subgrade filling HC-2(c) Subgrade filling CS-1(d) Subgrade filling CS-2(a)
It can be seen from Figure10 that in the cold season, when the height of the subgrade increases from 3 m to 3.5 m, 4 m, 4.5 m, and 5 m, the safety factor gradually decreases(b)
As can be seen from Figure10(a), when the height of the HC-1 soil-filled subgrade increased from 3 m to 3.5 m, 4 m, 4.5 m, and 5 m, the safety factors in December were the lowest of the four months in the cold season, which dropped by 29.4% from about 3.4 to around 2.4; the safety factor in February was the highest of the four months in the cold season, which increased by 28.2% from about 3.9 to about 2.8(c)
As can be seen from Figure10(b), when the height of the HC-2 soil-filled subgrade increased from 3 m to 3.5 m, 4 m, 4.5 m, and 5 m, the safety factor in December was the lowest of the four months in the cold season, which decrease by 29.4% from about 3.4 reduced to about 2.4; the safety factor in February was the highest of the four months in the cold season, which decrease by 30% from about 4.0 to about 2.8(d)
As can be seen from Figure10(c), when the height of the CS-1 soil-filled subgrade increased from 3 m to 3.5 m, 4 m, 4.5 m, and 5 m, the safety factor in December was the lowest of the four months in the cold season; it was reduced from about 3.6 to about 2.6, a decrease of 27.8%; the safety factor in February was the highest in the four months of the cold season, from about 4.2 to about 3, a decrease of 28.6%.(e)
As can be seen from Figure10(d), when the height of the subgrade filled with CS-2 soil is increased from 3 m to 3.5 m, 4 m, 4.5 m, and 5 m, the safety factor in December is the lowest of the four months in the cold season, from about 3.8 to about 3.0, a decrease of 21.1%; the safety factor in February was the highest of the four months in the cold season, from about 4.5 to about 3.3, a decrease of 26.7%.(f)
As the height of the subgrade increases (3~5 m), the safety factor decreases. In December, HC-1 and HC-2 soil-filled subgrade slopes experienced the fastest decline in safety factor, and in February, the HC-2 soil-filled subgrade slopes experienced the fastest decline in safety factor. The stability of the slope also decreases with the increase of the height of the subgrade, but the minimum value of all safety factors is greater than 2. When the height of the subgrade is 3~5 m, the slope is stable(3) Influence of Subgrade Slope on Its Slope Stability. The subgrade fillings are HC-1, HC-2, CS-1, and CS-2, the width of the road surface is 26 m, and the height of the subgrade is 4 m. When the slopes are selected to be 1 : 1.5, 1 : 1.75, and 1 : 2, the safety factors of the subgrade slopes in December, January, and March are shown as (a), (b), (c), and (d) in Figure 11.Figure 11
Different subgrade filling, slope safety factors with different gradients in the cold season.
(a)
Subgrade filling HC-1(b)
Subgrade filling HC-2(c)
Subgrade filling CS-1(d)
Subgrade filling CS-2(a) Subgrade filling HC-1(b) Subgrade filling HC-2(c) Subgrade filling CS-1(d) Subgrade filling CS-2(a)
As can be seen from Figure11, in the cold season, if the slope reduced from 1 : 5 to 1 : 1.75 and 1 : 2, the slope safety factor generally increases(b)
As can be seen from Figure11(a), if the slope of the subgrade filled with HC-1 soil decreased from 1 : 5 to 1 : 1.75 and 1 : 2, the slope safety factor in February is higher than the other months in the cold season, and the minimum value It is about 3.1, the maximum value is about 3.3, which is an increase of 6%; the slope safety factor in December is lower than other months in the cold season, the minimum value is about 2.7, the maximum value is about 2.9, and an increase of 7.4%(c)
As can be seen from Figure11(b), when the slope of the subgrade filled with HC-2 soil reduced from 1 : 5 to 1 : 1.75 and 1 : 2, the slope safety factor in February was higher than that in other months in the cold season, and the minimum value is about 3.1, the maximum value is about 3.3, which increased by 6%; the slope safety factor in December is lower than other months in the cold season, the minimum value is about 2.7, the maximum value is about 3.1, and an increase of 14.8%(d)
As can be seen from Figure11(c), when the slope of the subgrade filled with CS-1 soil slowed down from 1 : 5 to 1 : 1.75 and 1 : 2, the slope safety factor in February was higher than that in other months in the cold season, and the minimum value is about 3.4, the maximum value is about 3.6, which increased by 5.9%; the slope safety factor in December is lower than other months in the cold season, with the minimum value of about 3.0, the maximum value of about 3.2, and an increase of 6.7%(e)
As can be seen from Figure11(d), when the slope of the subgrade filled with CS-2 soil slowed down from 1 : 5 to 1 : 1.75 and 1 : 2, the slope safety factor in February was higher than that in other months in the cold season, and the minimum value is about 3.6, the maximum value is about 3.7, an increase of 2.8%; the slope safety factor in December is lower than other months in the cold season, the minimum value is about 3.2, the maximum value is about 3.3, and an increase of 3.1%(f)
As the slope of the subgrade slows down (1 : 1.5 to 1 : 2), the safety factor increases. The fastest increase in the safety factor in February was the subgrade slopes filled with HC-1 and HC-2 soil, and the fastest increase in the safety factor in December was in the subgrade slopes filled with HC-2 soil. The stability of the slope also increases as the slope of the subgrade slows down. All safety factors are greater than 2.7. When the slope of the subgrade is 1 : 1.5~1 : 2, its slope is stable(g)
Measures to improve the stability of subgrade slopes. During the construction process, the salt content of the saline soil subgrade fill is strictly controlled; insulate the slope of the subgrade to reduce the amplitude of the change of temperature gradient and prevent the influence of large temperature difference on the stability of the slope of the subgrade; geotextiles or geomembranes are placed between the different layers of the roadbed slope to block the migration of water and salts and prevent the accumulation of salts that can lead to the reduction of the stability of the subgrade slope


## 5.2.1. Stability Analysis of Subgrade Slope under Heat

(1) Monthly Changes in Slope Stability. The subgrade fills are HC-1, HC-2, CS-1, and CS-2. When the subgrade width is 26 m, the subgrade height is 4 m, and the slope is 1: 1.5, the subgrade slope safety factor is shown in Figure 9.
(a)
It can be seen from Figure9 that, in general, the slope safety factors of the four groups of subgrades filled with soil samples change during these months, and its trend conforms to the variation law of sinusoidal function curve. In February, the slope safety coefficients reached their maximum values, and in August, it reached their minimum values, but the minimum values were greater than 2. In each month of the year, the subgrade slopes were stable(b)
In Figure8, the safety coefficients of subgrade slopes in descending order throughout the year are CS-2, CS-1, HC-2, and HC-1. HC-1 soil-filled subgrade slope safety factor reached a maximum of about 3.1 in February and August reached a minimum value of about 2.0, which decreased by 35.5% in August compared with February; HC-2 soil-filled subgrade slope safety factor reached a maximum of about 3.2 in February, a minimum of about 2.1 in August, and a decrease of 34.4% in August compared with February; the safety factor of CS-1 soil subgrade slope reached a maximum of about 3.4 in February and a minimum of August. The value was about 2.2, which decreased by 35.3% in August compared with February; the safety factor of CS-2 soil-filled subgrade slope reached a maximum of about 3.6 in February, reached a minimum of about 2.6 in August, and decreased by 27.8 in August from February %. Among them, HC-1, HC-2, and CS-1 have similar declines, while CS-2 has a smaller decline(c)
Above all, the stability of the subgrade slope changes during the months in cold season. The stability in February is the best, and in August is relatively poor, but the minimum value of the safety factor is greater than 2. Comparing the information of air temperature and ground temperature, it can be concluded that the temperature below the surface reaches the minimum value in February, but reaches the maximum value in October. The temperature has a certain influence on the safety factor of the slope of saline soil subgrade, and the stability of the slope of saline soil subgrade is higher in the cold season than in the summer, but the slope of the subgrade is stable throughout the year, but the subgrade slope is stable throughout the year. Under the same construction conditions, the stability of subgrade slopes filled with saline soil with low chloride ion content is higher than that of saline soil with high chloride ion contentFigure 9
Monthly safety factor of subgrade slope.(2) Influence of Subgrade Height on Subgrade Slope Stability. Subgrade fillings are HC-1, HC-2, CS-1, and CS-2; the width of pavement is 26 m; slope of 1: 1.5; subgrade height is 3 m, 3.5 m, 4 m, 4.5 m, and 5 m, respectively; the safety factors of the subgrade slope in December, January, February, and March are shown in (a), (b), (c), and (d) in Figure 10.Figure 10
Different subgrade filling, slope 1 : 1.5, slope safety coefficient of different subgrade height in cold season.
(a)
Subgrade filling HC-1(b)
Subgrade filling HC-2(c)
Subgrade filling CS-1(d)
Subgrade filling CS-2(a) Subgrade filling HC-1(b) Subgrade filling HC-2(c) Subgrade filling CS-1(d) Subgrade filling CS-2(a)
It can be seen from Figure10 that in the cold season, when the height of the subgrade increases from 3 m to 3.5 m, 4 m, 4.5 m, and 5 m, the safety factor gradually decreases(b)
As can be seen from Figure10(a), when the height of the HC-1 soil-filled subgrade increased from 3 m to 3.5 m, 4 m, 4.5 m, and 5 m, the safety factors in December were the lowest of the four months in the cold season, which dropped by 29.4% from about 3.4 to around 2.4; the safety factor in February was the highest of the four months in the cold season, which increased by 28.2% from about 3.9 to about 2.8(c)
As can be seen from Figure10(b), when the height of the HC-2 soil-filled subgrade increased from 3 m to 3.5 m, 4 m, 4.5 m, and 5 m, the safety factor in December was the lowest of the four months in the cold season, which decrease by 29.4% from about 3.4 reduced to about 2.4; the safety factor in February was the highest of the four months in the cold season, which decrease by 30% from about 4.0 to about 2.8(d)
As can be seen from Figure10(c), when the height of the CS-1 soil-filled subgrade increased from 3 m to 3.5 m, 4 m, 4.5 m, and 5 m, the safety factor in December was the lowest of the four months in the cold season; it was reduced from about 3.6 to about 2.6, a decrease of 27.8%; the safety factor in February was the highest in the four months of the cold season, from about 4.2 to about 3, a decrease of 28.6%.(e)
As can be seen from Figure10(d), when the height of the subgrade filled with CS-2 soil is increased from 3 m to 3.5 m, 4 m, 4.5 m, and 5 m, the safety factor in December is the lowest of the four months in the cold season, from about 3.8 to about 3.0, a decrease of 21.1%; the safety factor in February was the highest of the four months in the cold season, from about 4.5 to about 3.3, a decrease of 26.7%.(f)
As the height of the subgrade increases (3~5 m), the safety factor decreases. In December, HC-1 and HC-2 soil-filled subgrade slopes experienced the fastest decline in safety factor, and in February, the HC-2 soil-filled subgrade slopes experienced the fastest decline in safety factor. The stability of the slope also decreases with the increase of the height of the subgrade, but the minimum value of all safety factors is greater than 2. When the height of the subgrade is 3~5 m, the slope is stable(3) Influence of Subgrade Slope on Its Slope Stability. The subgrade fillings are HC-1, HC-2, CS-1, and CS-2, the width of the road surface is 26 m, and the height of the subgrade is 4 m. When the slopes are selected to be 1 : 1.5, 1 : 1.75, and 1 : 2, the safety factors of the subgrade slopes in December, January, and March are shown as (a), (b), (c), and (d) in Figure 11.Figure 11
Different subgrade filling, slope safety factors with different gradients in the cold season.
(a)
Subgrade filling HC-1(b)
Subgrade filling HC-2(c)
Subgrade filling CS-1(d)
Subgrade filling CS-2(a) Subgrade filling HC-1(b) Subgrade filling HC-2(c) Subgrade filling CS-1(d) Subgrade filling CS-2(a)
As can be seen from Figure11, in the cold season, if the slope reduced from 1 : 5 to 1 : 1.75 and 1 : 2, the slope safety factor generally increases(b)
As can be seen from Figure11(a), if the slope of the subgrade filled with HC-1 soil decreased from 1 : 5 to 1 : 1.75 and 1 : 2, the slope safety factor in February is higher than the other months in the cold season, and the minimum value It is about 3.1, the maximum value is about 3.3, which is an increase of 6%; the slope safety factor in December is lower than other months in the cold season, the minimum value is about 2.7, the maximum value is about 2.9, and an increase of 7.4%(c)
As can be seen from Figure11(b), when the slope of the subgrade filled with HC-2 soil reduced from 1 : 5 to 1 : 1.75 and 1 : 2, the slope safety factor in February was higher than that in other months in the cold season, and the minimum value is about 3.1, the maximum value is about 3.3, which increased by 6%; the slope safety factor in December is lower than other months in the cold season, the minimum value is about 2.7, the maximum value is about 3.1, and an increase of 14.8%(d)
As can be seen from Figure11(c), when the slope of the subgrade filled with CS-1 soil slowed down from 1 : 5 to 1 : 1.75 and 1 : 2, the slope safety factor in February was higher than that in other months in the cold season, and the minimum value is about 3.4, the maximum value is about 3.6, which increased by 5.9%; the slope safety factor in December is lower than other months in the cold season, with the minimum value of about 3.0, the maximum value of about 3.2, and an increase of 6.7%(e)
As can be seen from Figure11(d), when the slope of the subgrade filled with CS-2 soil slowed down from 1 : 5 to 1 : 1.75 and 1 : 2, the slope safety factor in February was higher than that in other months in the cold season, and the minimum value is about 3.6, the maximum value is about 3.7, an increase of 2.8%; the slope safety factor in December is lower than other months in the cold season, the minimum value is about 3.2, the maximum value is about 3.3, and an increase of 3.1%(f)
As the slope of the subgrade slows down (1 : 1.5 to 1 : 2), the safety factor increases. The fastest increase in the safety factor in February was the subgrade slopes filled with HC-1 and HC-2 soil, and the fastest increase in the safety factor in December was in the subgrade slopes filled with HC-2 soil. The stability of the slope also increases as the slope of the subgrade slows down. All safety factors are greater than 2.7. When the slope of the subgrade is 1 : 1.5~1 : 2, its slope is stable(g)
Measures to improve the stability of subgrade slopes. During the construction process, the salt content of the saline soil subgrade fill is strictly controlled; insulate the slope of the subgrade to reduce the amplitude of the change of temperature gradient and prevent the influence of large temperature difference on the stability of the slope of the subgrade; geotextiles or geomembranes are placed between the different layers of the roadbed slope to block the migration of water and salts and prevent the accumulation of salts that can lead to the reduction of the stability of the subgrade slope


## 6. Conclusion

The author used COMSOL simulation to simulate the temperature field of the subgrade in Qarham area. With that, an isothermal layering-strength parameter model is proposed. The author used the SLOPE/W module in GeoStudio to calculate the safety factor of the subgrade slope in the cold region. The main conclusions are:(1)
The numerical simulation results are consistent with the numerical change trend of the measured data of boreholes, and the validity of the temperature field simulation is verified(2)
As the temperature changes, the strength of the subgrade soil changes accordingly. Based on the establishment of isothermal layering-strength parameter model, the analysis towards stability process of the subgrade slope in the cold region is done(3)
The monthly variation curve of the safety factor of the subgrade slope is very close to the sinusoidal function curve. With the season changes, the safety factor also changes; the stability gradually decreases with the increase of the subgrade height in the cold season, of which HC-2 fill of subgrades fell the most in February, about 30%. As the slope of the subgrade reduced (1 : 1.5, 1 : 1.75, 1 : 2), the stability gradually increased. The subgrade filled with HC-2 had the largest increase in February, about 14.8%.



---

*Source: 3141964-2022-04-04.xml*