Abstract

In the context of climate warming and the frequent wind-sand hazards in the Qinghai-Tibet Engineering Corridor (QTEC), the construction of the embankment will affect the thermal regime of permafrost underground. The influence of embankment construction on the variation of the permafrost table beneath it is different, especially for the regime with different mean annual ground temperatures (MAGTs). In this study, the effects of the unilateral sand particles deposition on the spatial distribution and temporal evolution pattern of temperature beneath the embankment are investigated through the numerical simulations, in which the heat transfer is considered. The model is validated by the field observed data of soil temperatures around an experimental zone built at the sand hazard area in Honglianghe, the interior of Qinghai-Tibet Plateau (QTP). The simulated results indicate that the temperature field beneath the embankment is asymmetrically distributed under the condition of unilateral sand particles deposition. This asymmetry gradually weakened with the increase of operation time and the gradual adjustment of the permafrost temperature field. By comparing the permafrost table beneath the natural surface, the sand deposition center, and the middle of the embankment center, it could be found that the unilateral sand particles deposition has less effect on the degradation of the permafrost table in the center of the embankment. However, for the center of the sand deposition, the change of the permafrost table is larger with the increase of time and the corresponding rate of permafrost table degradation is higher than that without sand particles deposition, especially for the high-temperature permafrost. In addition, with different sand thickness and width conditions, the effect of “narrow-thick” form sand particles deposition on the temperature field beneath embankment is greater than that of “wide-thin” form sand deposition. Hence, in order to reduce its impact on the long-term thermal condition beneath the embankment, it is necessary to clean the thicker deposition sand particles at the toe of the embankment.

1. Introduction

The Qinghai-Tibet Engineering Corridor (QTEC) is the most important link between the Tibet Autonomous Region and inner China, the width of which is several kilometers at wider sections but only hundreds of meters at narrow section [1]. Within the corridor, there are several major linear infrastructures including the Qinghai-Tibet Highway, Qinghai-Tibet Railway, Qinghai-Tibet Power Transmission Line, and Golmud to Lhasa Oil Pipeline and some optical communication cables. In recent years, human activities and increasing air temperature have caused significant grassland degradation and desertification in the QTEC [25]. Following the impacts, particularly the construction of the highway and the railway embankments, the energy of near-surface wind flow within the corridor has been disturbed, as well as the original relatively stable sand movement [6]. Hence, the wind-sand hazards were very serious around the highway and railway embankments, particularly in some river valleys and basins that the highway and railway traverse (Figure 1).

In the permafrost, the construction of the embankment will disturb the heat balance of the original surface. The strong heat absorption of asphalt pavement will lead to more rapid permafrost degradation beneath the highway [79]. In addition, the change of the thermal-mechanical properties of permafrost will cause embankment and pavement damage, including longitudinal cracks and uneven settlement [1012]. Hence, in order to reduce the impact of the human activities, a lot of effective engineering structural measures have been considered and used in the construction of highway and railway embankment including the crushed-rock and duct-ventilated embankments. Due to the excellence of the crushed-rock and duct-ventilated embankments, they have been widely used to ensure the thermal stability of highways and railways in permafrost regions of the QTEC [1319]. However, in the context of climate warming and the expanding scope of human activities, the wind-sand hazards around the embankment in the QTEC have become more and more serious in the last few years [2023]. The sand particles deposition around the embankment will not only change the surface boundary but also affect the cooling effect of the construction of crushed-rock or duct-ventilated embankments [2426]. Research studies have shown that with the increasing thickness of sand in the rock layer, the critical temperature difference between the sand-free layer increases, and the Ra number decreases, and the natural convection intensity weakens gradually [2730]. However, in these research studies, several assumptions about the sand deposition around the embankment are proposed. The distribution of sand particles deposition is generally assumed to be the same at the two slopes of the embankment. The prevailing wind direction within the corridor is west and northwest, and the corridor is from northeast to southwest [1]. Then, the two slopes of embankments of both the railway and highway are generally windward and leeward slopes, respectively. Thus, there should be a considerable difference in the sand deposition on the two slopes. However, the difference was not considered within these previous research studies, which will undermine the accuracy of evaluation on thermal impacts of sand particles deposition on the embankments installed with air-cooled structures.

In this study, a two-dimensional numerical model of the heat transfer for highway embankment is established. In the numerical model, three different mean annual ground temperatures (MAGTs) of the permafrost, in which the embankment is located, are considered, as well as the unilateral sand particles deposition at the toe of the embankment. Field observed data of soil temperatures around an experimental zone are used to validate the numerical model. Through numerical simulations, the effects of unilateral sand deposition on the spatial distribution and temporal evolution pattern of temperature beneath the embankment are investigated. The results of this study could provide informative references for highways constructed in permafrost zones, in which the wind-sand hazards frequently occur around the embankment.

2. Methods

2.1. Governing Equations

According to the experiment, in freeze-thaw soil layers, the ratio of heat conduction is far larger than that of the heat convention [31]. In this article, by ignoring the heat convection, the heat transfer progress considering the heat conduction and phase change in freeze-thaw soil layers can be described as follows [3234]:where is the equivalent volume heat capacity and is the equivalent thermal conductivity. According to the method of sensible heat capacity, the equivalent volume heat capacity and the equivalent thermal conductivity in freeze-thaw soils can be written as follows [3234]:where is the temperature range of phase change; and are the volume heat capacity and thermal conductivity of unfrozen soil; and are the volume heat capacity and thermal conductivity of frozen soil; L is the latent heat of phase change per unit volume.

2.2. Physical Model

A physical model of a highway without sand particles deposition near the toe of the embankment (NSE) is shown in Figure 2(a), and a highway with sand particles deposition near the toe of the embankment (SE) in Figure 2(b). In the model, part I is embankment fill, part II (0∼3 m) is gravel and clayey layer, part III (3∼8 m) is silty clay, and part IV(8∼30 m) is weathered mudstone. The height of the embankment is 4 m and the width of its paving is 10 m. The gradient of the slope is 1 : 1.5. The widths of the sand deposition are determined as 7.1, 10, and 16.7 m. The thicknesses of sand particles deposition are determined as 0.3, 0.5, and 0.7 m. The computational domain is extended by 30 m wide from the outside slop toe of the embankment in the horizontal direction and 30 m height beneath the natural surface in the vertical direction. The thermal parameters of soil layers are given in Table 1 [33].

2.3. Boundary

According to the IPCC report, the air temperature in Qinghai-Tibetan Plateau (QTP) will be warmed up by 2.6°C in the future 50 years because of climate change [35]. Based on the adhered layer theory [12, 36], the thermal boundary conditions of the computational domain are expressed as follows:

The temperature at natural surfaces of the NSE model (AB and EF) and the SE model (AM and EF) is as follows:

The temperature at the side slopes of the NSE model (BC and DE) and the SE model (OC and DE) is as follows:

The temperature at the asphalt pavement surface of the NSE model and SE model (CD) is as follows:

The temperature at the sand particles deposition surface of the SE model (NO and NM) is as follows:where is the time; Ta is the mean annual air temperature at which the embankment is located, being determined as −3.0, −3.5, and −4.0°C, respectively. The geothermal heat flux of 0.03 W/m2 is applied to the bottom boundary (IJ) in both models. The lateral boundaries (ALKI and FGHI) are assumed to be adiabatic. With the governing equations and boundary conditions above, the problem is solved numerically using the commercial software of Fluent 14.0.

The temperature boundary of the natural ground surface without consideration of climate warming is used to calculate the initial temperature fields beneath the embankment (parts II, III, and IV). The obtained stable temperature fields on July 15 are taken as the initial temperature condition of these parts, as shown in Figure 3.

2.4. Model Validation

To validate the reliability of the computational model and parameters, numerically simulated results of soil temperature beneath the sand particles deposition surface are compared with the field measured data on April 15 and July 15 [37], as shown in Figure 4. It could be seen that the agreement between field measured data and numerically simulated results is good. However, there are some slight discrepancies between the simulated results and field measured data, especially on July 15. The simulated results of soil temperature are lower than those of the measured data within the active layer. The discrepancies may be attributed to the simplifications of boundary and soil strata in the computational model. Overall, the comparison shows that the computational model and the parameters can be used for simulating the spatial distribution and temporal evolution pattern of temperature beneath the embankment in permafrost zones.

3. Results and Analysis

3.1. Variation of Temperature beneath the Embankment without Sand Deposition

For the NSE, without consideration of the sunny-shady slope effect, Figure 5 shows soil temperature distribution beneath the NSE, in which the embankment is located with different MAGTs, in October 15 in the 5th, 20th, and 50th year after the embankment construction. The soil temperature distribution is given from the centerline of the embankment to 20 m away from the slop toe (Parts 1–3 in Figure 2(a)).

In the 5th year after the embankment construction (Figure 5(a)), the distribution pattern of the ground temperature field under the NSE in different MAGT permafrost zones is basically the same. The permafrost table beneath the centerline of the embankment (CE) is higher than that beneath the natural surface (NS). As for the permafrost zone with different MAGTs of −0.5, −1.0, and −1.5°C, the increase of the permafrost table is 0.3, 1.6, and 2.0 m, respectively. In addition, due to the increase in the thermal area of the asphalt pavement of the embankment and the structure, the internal heat absorption of the embankment will increase sharply. The trend of the geothermal line changes in different permafrost zones. In the permafrost zone with the MAGT of −0.5°C, it can be seen that the depth of the geothermal line of −0.4°C beneath the NS is the same as that beneath the CE. In contrast, in the permafrost zone with two MAGTs of −1.0 and −1.5°C, the depth of the geothermal line of −0.4°C beneath the NS is lower than that beneath CE (Figure 5(a)). It reveals that the thermal disturbance to different MAGT permafrost zones is different in the early stage of embankment construction.

With the operational time increasing, the permafrost table and the permafrost warming beneath the NS and the embankment decline, especially for the CE. The internal heat absorption of the embankment will lead to a significant downward trend of the geothermal line beneath the CE. In the 20th year after construction of the embankment, the permafrost table beneath the NS in the permafrost zone with the MAGT of −0.5, −1.0, and −1.5°C is 2.7, 2.2, and 2.0 m, respectively. Meanwhile, the permafrost table beneath the CE in the permafrost zone with the MAGT of −0.5, −1.0, and −1.5°C is 3.8, 2.6, and 1.2 m, respectively (Figure 2(b)). In the 50th year after the construction of the embankment, due to the “heat gathering” effect of the asphalt pavement and the structure of the embankment, the permafrost table beneath the CE in permafrost zones with three different MAGTs is lower than that beneath the NS. The difference in the permafrost table is 5.5, 2.1, and 1.7 m, respectively (Figure 5(c)).

Hence, considering the warming of the permafrost beneath NS and the CE in the context of climate warming and human engineering activities, the change rate of the permafrost can be used to evaluate the degradation of the permafrost. The variation of the permafrost table beneath the CE and the NS is listed in Table 2. As for the three different MAGT permafrost zones, it can be seen that the higher the MAGT is, the greater the rate of degradation of the permafrost beneath the CE and NS will be. In addition, with the operation of the embankment, the degradation of permafrost beneath the CE is higher than that of beneath the NS, which shows the thermal influence of the embankment.

3.2. Variation of Temperature beneath the Embankment with Unilateral Sand Deposition

As an obstacle structure, the construction of the embankment will change the progress of the initial flow field, including the wind flow and wind-sand flow. The redistribution of the wind-sand flow around the embankment will influence the capacity of the flow field, resulting in the sand particles deposition around the embankment. In this part, considering the influence of the ambient wind speeds, the difference of sand particles deposition around the embankment is significant, which may affect the temperature beneath the embankment. Hence, the influence of the unilateral sand particles deposition (USPD) at the toe of the embankment with the height of 4 m is investigated in this part. The width and thickness of the sand layer are 10 m and 0.5 m, respectively. The variations of the temperature beneath the embankment over the simulation period of 50 years for three initial MAGTs of −0.5, −1.0, and −1.5°C are shown in Figures 68.

In the 5th year after the sand particles deposition at the toe of the embankment, the USPD will induce the asymmetrical distribution of the temperature beneath the embankment. As for three different MAGTs, the lower the MAGT of the permafrost zone is, the more significant the asymmetry along the depth direction will be. From Figure 6, it can be seen that due to the USPD at the toe of the embankment, the maximum depth of zero annual amplitude of ground temperature is offset from the CE. The offsets are 5, 7.5, and 13.0 m for three different MAGTs, respectively (Figure 6). Otherwise, the permafrost table beneath the centerline of the sand particles deposition (CSD) has moved 1.0, 1.3, and 1.4 m above the NS for three different MAGTs, respectively. It reveals that the sand particles deposition can delay the degradation of the permafrost in the early stage of the sand particles deposition.

With sand particles deposition time increasing at the toe of the embankment, the influence of the USPD for the temperature beneath the embankment is different in the three different MAGT permafrost zones. In the 20th year after the sand particles deposition at the toe of the embankment, in the permafrost zone with the MAGTs of −0.5°C, the temperature field beneath the embankment is symmetrically distributed with the maximum depth of zero annual amplitude of ground temperature locating at the CE. The influence of the USPD on the temperature beneath the embankment gradually weakens. In contrast, for the permafrost zones of the MAGT of −1.0 and −1.5°C, the asymmetrical distribution of the temperature beneath the embankment is still significant. The offsets of the centerline are 3.0 and 5.0 m for two different MAGTs, respectively (Figure 7). With the sand particles deposition for 20 years, the permafrost table beneath the CSD has moved 0.4, 1.1, and 1.3 m upward above the NS for three different MAGTs of −0.5, −1.0, and −1.5°C, respectively. It reveals that with the time of USPD increasing, the sensitivity of the permafrost to sand particles deposition gradually weakens and the protective effect of sand particles deposition on permafrost gradually disappears (Figure 8). In addition, in the 50th year after the sand particles deposition at the toe of the embankment, the influence of the USPD at the toe of the embankment could be ignored for three different MAGT permafrost zones. The distribution of the temperature beneath the embankment is symmetrical. However, with the sand particles deposition lasting for 50 years, the permafrost table beneath the CSD is lower than that beneath the NS in the permafrost zone with MAGT of −0.5°C. Meanwhile, the permafrost table beneath the centerline of the sand particles deposition is basically the same as the NS for the two MAGTs of −1.0 and −1.5°C, respectively (Figure 9). Thus, by comparing the variation of the permafrost beneath the sand particles deposition and the NS in three different MAGTs, it could be found that with the time of sand particles at the toe of the embankment increasing, the asymmetrical distribution of the temperature beneath the embankment is gradually weakening. What is more, the influence of the sand particles deposition can mitigate the degradation of the permafrost beneath it, especially for the permafrost zone with the MAGRs of −1.0 and −1.5°C. Within 50 years of the sand particles deposition, the permafrost table beneath the CSD is higher than that beneath the NS (Table 3).

3.3. Effect of Unilateral Sand Particles Deposition on Temperature at Different Locations around the Embankment

The variations of the temperature beneath the CE and the toe of the embankment (TES), at which the sand particles are deposited, are shown in Figures 9 and 10. It could be seen that with the increase of time, the variations of the temperature beneath different surface boundaries in three different MAGT permafrost zones vary significantly. In Figures 9 and 10, it can be seen that 50 years after the sand particles deposition at the toe of the embankment, the maximum depth of zero annual amplitude of ground temperature beneath the CE and TES is about 16 and 14 m, respectively. Compared with the embankment without sand particles deposition at the toe, the sand particles deposition has little influence on the value of the maximum depth of zero annual amplitude of ground temperature.

Additionally, Tables 4 and 5 show the variation of the permafrost table beneath the CE and CES (the center of the embankment without sand particles deposition), TES and the toe without sand particles deposition (TE) in October of the 5th, 20th, and 50th years with three different MAGTs. Comparing the permafrost table beneath the CE and CES with different MAGTs and times, it can be found that the USPD at the toe of the embankment has little influence on the variation of the permafrost table beneath the CE (Table 4). However, as for the toe of the embankment (Table 5), 5 years after the sand particles deposition, the permafrost table beneath the toe of TES in three different MAGTs is 1.52, 1.04, and 0.64 m, the value of which is 0.67, 0.80, and 0.91 m higher than that beneath the TE, respectively. With the sand particles deposition increasing to 20 years, the permafrost table beneath the TES is 0.29, 0.68, and 0.81 m higher than that beneath the TE, respectively. Moreover, after 50 years of the sand particles deposition, the value of the increase is −0.21, 0.00, and −0.15 m, respectively. Hence, from the result of the comparison between the permafrost table beneath the TES and TE, it could be found that as for the three different MAGTs permafrost zones, the USPD at the toe of the embankment has a significant influence on the permafrost table beneath the toe of embankment. From the 5th year of sand particles deposition to 20th year, the sand particles deposition could improve the permafrost table beneath the toe of the embankment, especially for the permafrost zone with the MAGT of −1.5°C. However, with the time of sand particles deposition increasing, the influence of the sand particles deposition may induce the temperature of permafrost to increase and accelerate the degradation of the permafrost. Hence, as for the permafrost zone with different MAGTs, the adoption of sand control measures is beneficial to protect the permafrost by considering the time of sand particles deposition around the embankment.

3.4. Effect of Different Forms of Sand Particles Deposition on the Temperature of Permafrost

Based on the study in the previous section, it can be concluded that within 20 years of the sand particles deposition, the sand particles deposition has a significant protective effect on the permafrost beneath it. Hence, in order to study the thermal effects of different forms of sand particles deposition (thick and thin of the sand layer) on the permafrost, under the condition of the MAGT −1.0°C, the thermal effects of sand particles deposition on the permafrost ground temperature with different thin and thick of the sand layer are given in Figures 11 and 12.

In Figure 11, with the thickness of the sand layer being 0.3 m, the isotherms beneath the toe of the embankment with sand particles deposition vary significantly, inducing the asymmetrical distribution of the geothermal field beneath the embankment. The depths of −1.0°C isotherm beneath the NS are compared, which are 5 m away from the two toes of the embankment. The depth of the −1.0°C isotherm beneath the toe without sand particles deposition has moved 6.5, 4.0, 3.3, and 1.7 m upward compared with the toe with 0.3 m thick sand layer after the sand particles deposition of 5, 10, 15, and 20 years, respectively. In addition, in Figure 12, with the thickness of sand layer being 0.7 m, it can be found that the distribution of the ground temperature field is basically the same as that of the 0.3 m condition. As for the −1.0°C isotherm, the depth of it beneath the toe without sand particles deposition has moved 4.8, 2.8, 1.4, and 1.3 m upward compared with the toe with 0.7 m thick sand layer after the sand particles deposition of 5, 10, 15, and 20 years, respectively. This means that with the time of the sand particles deposition increasing, the “wide and thin” type of sand particles deposition has a greater influence on the depth of the lower ground temperature field than that of the “narrow and thick” type of sand particles deposition.

To investigate permafrost warming beneath different thin and thick sand layers in the context of climate warming, the variations of the permafrost table beneath the CE, TE, TES, and CSD are shown in Tables 6 and 7 within 20 years of the sand particles deposition. It could be found that the two different sand particles deposition forms have less effect on the variation of permafrost table beneath CE and TE. However, there is a significant effect on the permafrost table beneath the CSD and TES. In Table 6, with the thickness of the sand layer being 0.3 m, the degradation rate of the permafrost beneath the CSD, TES, and TE is 0.03, 0.05, and 0.04 m·a−1, respectively. Additionally, in Table 7, with the thickness of the sand layer being 0.7 m, the degradation rate of the permafrost beneath the CSD, TES, and TE is 0.03, 0.04, and 0.04 m·a−1, respectively. This means that, within 20 years of sand particles deposition, no matter which form of the sand layer is, the deposition of the sand particles could weaken the degradation rate of the permafrost. However, as for the TES, the more the thickness of the sand layer, the smaller the degradation rate of permafrost.

As for the embankment in the permafrost zone, the main disease is thaw collapse, which is closely related to the variation of the permafrost table. After 5 years of the two forms of sand particles deposition (“wide and thin” and “narrow and thick”), the difference of the permafrost table, the location of which is beneath the two toes of the embankment, between the two toes of the embankment is 0.49 and 0.88 m, respectively (Tables 6 and 7). However, the degradation rate of the permafrost table with the sand layer thickness of 0.3 m is greater than that with the sand layer thickness of 0.7 m. With the sand particles deposition time increasing and the context of the climate warming, the difference of permafrost table between two toes of the embankment will induce uneven settlement significantly, especially for the embankment with the “narrow and thick” sand particles deposition at the toe of the embankment. Hence, in order to reduce its impact on the long-term thermal condition beneath the embankment, it was necessary to clean the thicker deposition sand particles at the toe of the embankment.

4. Conclusion

In permafrost areas, the engineering activities will have a considerable influence on the project and the underlying permafrost. As an engineering structure, the constructions of the road project will produce certain disturbance to the original relatively stable wind-sand flow along the route, leading to the redistribution of the flow field around the embankment and the emergence of sand particles deposition. The deposition and coverage of sand particles around the embankment can change the thermal conditions of the permafrost embankment. In this article, a mathematical model of heat transfer for freeze-thaw soil is constructed to investigate the long-term thermal effects of sand particles deposition at the toe of the embankment in the permafrost with different MAGTs. Based on the above numerical analyses and comparisons, the following conclusions can be drawn:(1)The thermal disturbance of the construction of the embankment in the permafrost zone with different MAGTs varies significantly. With time increasing, the degradation rate of the permafrost largely differs in different parts around the embankment. After 50 years of the construction of the embankment, the permafrost table beneath the CE in the permafrost zones with three different MAGTs is 11.5, 6.5, and 5.0 m. In contrast, the value of which beneath the NS is 6.0, 6.5, and 5.0 m. Therefore, the construction of the highway embankment will increase the heat absorption within the embankment and accelerates permafrost degradation.(2)The influence of the USPD will induce the distribution of the temperature beneath the embankment asymmetry significantly, especially for the embankment in the permafrost with different MAGTs. The lower the MAGT of the permafrost, the more significant the asymmetry of the temperature distribution. With the time of sand particles deposition increasing, the asymmetry of ground temperature field in three different MAGTs weakens gradually. Moreover, The influence of the sand particles deposition can mitigate the degradation of the permafrost beneath it, especially for the permafrost zone with the MAGRs of −1.0 and −1.5°C. Within 50 years of the sand particles deposition, the permafrost table beneath the CSD is higher than that beneath the NS.(3)From the variation of the temperature distribution at different parts of the embankment, the influence of the USPD on the permafrost table beneath the CE could be ignored. In contrast, for the TES, with the time of sand particles deposition increasing, the degradation rate of permafrost table at the TES is greater than that at the TE, especially for permafrost with the MAGT −0.5°C. After 50 years of the sand particles deposition at the toe of the embankment, the permafrost table has moved 0.21 m upward beneath the TE compared with that beneath the TES.(4)With different sand thickness and width conditions, the effect of “narrow and thick” form sand particles deposition on the temperature field beneath embankment was greater than that of “wide and thin” form sand deposition. Hence, in order to reduce its impact on the long-term thermal condition beneath the embankment, it was necessary to clean the thicker deposition sand particles at the toe of the embankment.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest related to this manuscript.

Acknowledgments

This research was supported by the Natural Science Foundation of China (41630636 and 41772325) and the Research Project of the State Key Laboratory of Frozen Soils Engineering (SKLFSE-ZY-17).