The comprehensive understanding of the variation law of soil thermal conductivity is the prerequisite of design and construction of engineering applications in permafrost regions. Compared with the unfrozen soil, the specimen preparation and experimental procedures of frozen soil thermal conductivity testing are more complex and challengeable. In this work, considering for essentially multiphase and porous structural characteristic information reflection of unfrozen soil thermal conductivity, prediction models of frozen soil thermal conductivity using nonlinear regression and Support Vector Regression (SVR) methods have been developed. Thermal conductivity of multiple types of soil samples which are sampled from the Qinghai-Tibet Engineering Corridor (QTEC) are tested by the transient plane source (TPS) method. Correlations of thermal conductivity between unfrozen and frozen soil has been analyzed and recognized. Based on the measurement data of unfrozen soil thermal conductivity, the prediction models of frozen soil thermal conductivity for 7 typical soils in the QTEC are proposed. To further facilitate engineering applications, the prediction models of two soil categories (coarse and fine-grained soil) have also been proposed. The results demonstrate that, compared with nonideal prediction accuracy of using water content and dry density as the fitting parameter, the ternary fitting model has a higher thermal conductivity prediction accuracy for 7 types of frozen soils (more than 98% of the soil specimens’ relative error are within 20%). The SVR model can further improve the frozen soil thermal conductivity prediction accuracy and more than 98% of the soil specimens’ relative error are within 15%. For coarse and fine-grained soil categories, the above two models still have reliable prediction accuracy and determine coefficient (R2) ranges from 0.8 to 0.91, which validates the applicability for small sample soils. This study provides feasible prediction models for frozen soil thermal conductivity and guidelines of the thermal design and freeze-thaw damage prevention for engineering structures in cold regions.
National Natural Science Foundation of China4150229251574037Applied Fundamental Research Project of China Communications Construction Co., Ltd.2018-ZJKJ-PTJS032016-ZJKJ-02China Postdoctoral Science Foundation2014M5607391. Introduction
With the increasing frequency of human activities and the trend of global warming, the temperature sensitive permafrost is changing significantly, most of which are in critical equilibrium state or even gradual degradation [1–3]. The widely distributed permafrost in the Qinghai-Tibet Plateau exhibits stronger temperature sensitivity and faster temperature rise rate than other permafrost regions with the same latitude [4–6]. The dynamic degradation process of frozen soil has a great impact on the stability of buildings in cold regions [7–10]. In the meanwhile, it will also lead to local environmental degradation, or further affect the global climate carbon cycle and climate change [11–14]. Therefore, many researchers had investigated the distribution, stability, and thermal regime of permafrost [15–17]. The thermal conductivity is one of the important thermophysical parameters for reconstructing the past and predicting the temperature status of permafrost under climate change conditions and also determining a parameter for engineering design in cold regions [18, 19]. Therefore, the comprehensive understanding of variation law of soil thermal conductivity is one of the important tasks in permafrost research.
Thermal conductivity is an inherent parameter to characterize the heat transfer performance of soil, which can usually be obtained by means of experimental test and prediction models [20–22]. The test of thermal conductivity can generally be divided into the steady-state method and transient method [23, 24]. Specifically, it includes steady-state comparison method, steady-state heat flow meter method, transient hot wire method, transient heat pulse method, and TPS method [25]. However, the disadvantages of the experimental test are usually strict technical requirements, long time consumption, and high cost, and sometimes limited test results cannot satisfy the needs of practical applications. Therefore, as reported [26, 27], many soil thermal conductivity prediction models had been proposed. Empirical models are usually based on statistical analysis of thermal conductivity experimental results and use water content, dry density, porosity, and other related soil properties as fitting parameters [28–31]. Kersten [28] analyzed the test results of 19 soil types and established the empirical formula of thermal conductivity using water content and dry density. Johansen [29] proposed the concept of normalized thermal conductivity and given the interpolation calculation model of soil thermal conductivity based on the relationship between normalized coefficient and soil saturation. Lu et al. [30] performed a series of thermo-TDR tests on twelve natural soils and proposed a linear prediction model across a wide range of soil moisture condition. Yan et al. [31] developed a generalized effective soil thermal conductivity model for soils of various textures from dry to saturation. Then, many scholars improved the model for wider applicability and higher prediction accuracy [32–35]. At the same time, many scholars have established theoretical predictive models of soil thermal conductivity based on specific physical theory and models [36–38]. Farouki [39] and Xu et al. [40], respectively, established the weighted calculation model and geometric average calculation formula based on the volume ratio and thermal conductivity of each component of the frozen soil on the basis of predecessors. Considering for the interactions among soil particle, water, and air of a soil unit cell, Haigh [41] derived a theoretical thermal conductivity model for sand soil. In recent years, with the maturity of machine learning and its advantages in dealing with complex nonlinear problems, many researchers had adopted it to the prediction of soil thermal conductivity [42, 43]. For instance, Bang et al. [44] investigated the application effect of linear regression and various machine learning methods in the prediction of thermal conductivity of compacted bentonite and verified the feasibility and superiority of the machine learning method in the prediction of the soil thermal conductivity.
Compared to unfrozen soil, the specimen preparation and experimental procedures of frozen soil thermal conductivity testing are more complex and challengeable. And it can be found that the prediction accuracy of the frozen soil thermal conductivity predictive model is lower than that of unfrozen soil [40]. Meanwhile, as a complex multiphase composition, previous research indicated that the thermal conductivity of frozen soil was associated with many factors [26, 27, 45].The existing frozen soil thermal conductivity empirical model (usually using dry density, water content, porosity, etc., as the fitting parameters) cannot comprehensively consider the influences of complex multiphase and porous structural characteristics of soil. Therefore, in the present work, prediction models of frozen soil thermal conductivity using nonlinear regression and SVR methods have been developed considering for essentially multiphase and porous structural characteristic information reflection of unfrozen soil thermal conductivity. Thermal conductivity of various types of soil samples which are sampled from the QTEC are tested by the TPS method. Correlation analysis of thermal conductivity of unfrozen and frozen soil has been adopted. Based on the measurement data of unfrozen soil thermal conductivity, the ternary fitting and SVR prediction models of frozen soil thermal conductivity for the typical soils in the QTEC are proposed. To further facilitate engineering applications, the prediction models of coarse and fine-grained soil categories have been proposed. Furthermore, the prediction models of two soil categories have also been used to predict the frozen soil thermal conductivity of small size soil samples in the QTEC.
2. Materials and Methods2.1. Soil Sample2.1.1. Collection of Soil Samples
The soil samples of thermal conductivity measurement were collected from the drilling specimens of Qinghai-Tibet expressway geological exploration project, which was implemented by the CCCC First Highway Consultants Co. Ltd. and Eco-Environment and Resources Northwest Institute of CAS from September 2017 to June 2018. As shown in Figure 1, the sampling spot mainly locates in the permafrost regions between Xidatan and Tanggula Mountain (the corresponding Qinghai-Tibet Highway mileage is K2870∼K3307). Every drilling spot is sampled for different depths, and the sampling depth varies from 0.5 m to 40 m. Total number of 638 unfrozen soil samples and 860 frozen soil samples are tested.
Schematic of the sampling spot along the QTEC [16].
Soil types of testing specimens are determined by the classification criteria of geotechnical engineering. Figure 2 shows the statistic results of different types of soil samples (only numbers more than 15 are given). It can be seen that silty clay is the most widely distributed frozen soil type in the QTEC, and the following soil types are sandy and gravel soil. Because some soil samples with large water (ice) content are impossible to implement unfrozen soil thermal conductivity testing, a total of 609 specimens tested for both frozen and unfrozen are selected as the research objects in this work, and the statistical distribution of various soil samples is shown in Table 1. It can be seen that the proportions of silty clay, silt, fine sand, gravel sand, boulder, breccia, and all weathered rock soils accounting for 86.2% are defined as the typical analytical soil types in the following predictive model research.
Statistic numbers of test samples for certain soil type.
Statistical distribution of various soil specimens.
Soil type
Numbers
Proportion (%)
Sampling depth (m)
≤1
1∼5
5∼10
10∼20
≥20
Fine-grained soil
Clay
24
3.94
1
6
3
9
5
Silty clay
202
33.17
10
63
40
69
20
Silt
52
8.54
6
23
10
8
5
Weathered slate
10
1.64
2
4
1
2
1
Weathered phyllite
29
4.76
0
7
10
10
2
Weathered sandstone
6
0.99
0
4
0
1
1
Weathered marl
5
0.82
0
1
2
1
1
Weathered mudstone
23
3.78
0
7
5
11
0
Coarse-grained soil
Silty sand
23
3.78
1
10
4
6
2
Fine sand
49
8.05
14
19
9
2
5
Medium sand
5
0.82
0
1
3
1
0
Coarse sand
3
0.49
0
3
0
0
0
Gravel sand
45
7.39
5
28
5
5
2
Boulder
49
8.05
3
31
5
5
5
Breccia
55
9.03
7
32
6
8
2
Gravel + pebble
10
1.64
0
4
1
3
2
Gravel soil
19
3.12
1
12
1
5
0
2.1.2. Water Content and Dry Density Distribution of Soil Samples
The dry density (ρd) and water content (w) of soil samples are statistically calculated, and the accumulative proportion distribution is shown in Figure 3. It shows that the water content and dry density distributions are significantly correlated with soil properties. The dry density is negatively correlated with soil particle size and basic orders of dry density of the samples are silt, silty clay, weathered rock, fine sand, gravel sand, breccia, and boulder soils. The average dry density of above soil samples are 1.64, 1.69, 1.79, 1.83, 1.87, 1.92, and 2.06 g/cm3, and the main distribution intervals (proportion within 10%∼90%) are 1.40∼1.81, 1.47∼1.91, 1.53∼2.04, 1.61∼2.03, 1.66∼2.15, 1.66∼2.20, and 1.81∼2.31 g/cm3. Likewise, the orders of water content are silty clay, silt, weathered rock, fine sand, breccia, gravel sand, and boulder. The corresponding average values are 19.33%, 19.05%, 16.04%, 13.27%, 10.88%, 10.66%, and 9.28%, and main distribution interval are 11.68%∼26.7%, 11.03%∼26.62%, 3.96%∼25.9%, 4.7%∼20.16%, 4.6%∼17.4%, 4.66%∼18.64%, and 4%∼15.2%, respectively, which basically covers all the water content ranges from dry to saturated.
(a) Dry density and (b) water content probability distribution of soil samples.
2.2. Experimental Methods and Apparatus
The TPS method is utilized as the thermal conductivity testing measure for all soil samples. The experimental apparatus is the Hot Disk 1500s Thermal Conductivity Analyzer (as shown in Figure 4), which has ±3% measurement accuracy. The principle of the TPS method for measuring the thermal properties of materials is based on the transient temperature response of a step-heated disc-shaped heat source in an infinite medium. No. 4922 Kapton film probe with radius of 29.4 mm is used to make specimen as large as possible, and the thermal conductivity (λ) is calculated by line equation of temperature variation values and dimensionless time constant which is measured by the film sensor:(1)ΔT¯τ=P0π3/2rλDτ,where ΔT¯τ is the temperature variation value, P0 is the total output of power and r is the radius of Kapton film sensor, and Dτ is the dimensionless time constant. During the test, experimental soil samples are formed into 3 cm height and 8 cm diameters column, and the film probe is sandwiched between two soil samples and fixed with a stainless steel sample holder. The detailed experimental procedures can be obtained in our previous work [46].
Thermal conductivity test system.
2.3. Nonlinear Regression Model
Previous studies have shown that dry density and water content have significant influence on soil’s thermal conductivity. Partial correlation analysis between frozen soil thermal conductivity of dry density, water content, and unfrozen soil thermal conductivity for two soil categories have been performed, and the results are shown in Table 2. The results exhibit that dry density and water content, as well as unfrozen soil thermal conductivity (λu), have significant positive correlations with frozen soil thermal conductivity (λf). In particular, it should be noted that the correlation between the frozen and unfrozen soil thermal conductivity for both coarse and fine-grained soils are very high (0.88 and 0.93, respectively). It can be inferred that, considering for the essential soil properties information reflection of unfrozen soil thermal conductivity, it would be possible to utilize the unfrozen soil thermal conductivity to predict frozen soil thermal conductivity.
Partial correlation analysis between frozen soil thermal conductivity, dry density, water content, and unfrozen soil thermal conductivity.
Soil categories
Related variables
Control variables
Unfrozen soil thermal conductivity
Frozen soil thermal conductivity
Fine-grained soil
Water content
Dry density
0.312
0.450
Unfrozen soil thermal conductivity
—
0.880
Dry density
Water content
0.456
0.385
Unfrozen soil thermal conductivity
—
0.891
Coarse-grained soil
Water content
Dry density
0.501
0.608
Unfrozen soil thermal conductivity
—
0.927
Dry density
Water content
0.530
0.502
Unfrozen soil thermal conductivity
—
0.930
Curve estimation of the fitting relationship between the frozen soil thermal conductivity and three factors has been taken, and it is found that the thermal conductivity of the frozen soil is linear with the unfrozen soil thermal conductivity, while its relation with dry density and water content are in the form of logarithmic functions. The ternary fitting formulas of frozen soil thermal conductivity for 7 typical soils in the QTEC are given as follows:(2)λf=a+bλu+clnw+dlnρd,where a, b, c, and d are the fitting coefficients of the equation. The detail fitting results of 7 typical soils are listed in Table 3. It can be seen that the determine coefficients (R2) of 7 typical soils range from 0.76∼0.93, while the analog effect of sandy soils are better than other soil types.
Ternary fitting result of frozen soil thermal conductivity for 7 typical soils in the QTEC.
Soil type
Fitting coefficients
R2
a
b
c
d
Silty clay
−0.703
1.291
0.238
0.007
0.759
Silt
−1.358
1.283
0.424
0.248
0.792
Weathered rock
−0.102
1.370
0.183
−0.854
0.895
Fine sand
−0.220
1.290
0.146
−0.429
0.908
Gravel sand
−0.617
1.272
0.254
−0.138
0.933
Boulder
−1.257
1.049
0.497
0.702
0.865
Breccia
−1.282
0.884
0.459
1.239
0.78
In order to analysis the predictive effect of the ternary fitting model, the estimated frozen soil thermal conductivity based on the Kersten model [28], Gangadhara Rao and Singh model [47], and ternary fitting model are plotted against the measured value of testing dataset. The empirical formulas proposed by Kersten and Gangadhara Rao and Singh used dry density and water content as fitting parameters, while the applicable soil types are silt, fine sand and gravel sand soils, etc. The comparison results of the estimation values of three models and the experimental values are shown in Figure 5. It can be seen that the proposed ternary fitting model performs best among all three models as most of the predictive values are within the ±10% relative error line. However, the predictive values of Kersten and Gangadhara Rao and Singh models are generally overestimated for all soil types and their prediction errors reach up to 40%. Therefore, it can be concluded that, compared with nonideal prediction accuracy of binary fitting models, the ternary fitting model using the unfrozen soil thermal conductivity as the fitting parameter has a higher thermal conductivity prediction accuracy.
A comparison of measured and estimated frozen soil thermal conductivity. (a) Slit. (b) Fine sand. (c) Gravel sand. (d) Fine-grained sand.
3. Development of SVR Models3.1. Theory and Performance Evaluation of SVR Model
SVR method is a type of supervised machine learning method and proposed by Vapnik et al. [48, 49], which is based on statistical learning theory. The SVR method uses kernel functions to map low-dimensional nonlinear problems to high-dimensional space to achieve linear separability and then to seek linear regression equations to fit sample data (as shown in Figure 6). The regression function can be expressed as follows:(3)fx=KTφx+m,where x is the input vector, K is the weight vector, m is the offset vector, and φx is eigenvector that maps input data to the high-dimensional space. The Gaussian radial basis function is chosen as the kernel function in the calculation.
Schematic diagram of the SVR method.
The SVR method is based on the principle of minimizing structural risk, transforming the linear regression problem into the following optimization problem, and then the values of K and m can be determined [48]:(4)Minimize12KTK+t∑i=1nξi+ξi∗,Subject toyi−KTφx−m≤ε+ξiKTφx+m−yi≤ε+ξi∗ξi,ξi∗≥0,where t is the penalty parameter which is greater than 0, ε is the insensitive loss coefficient, and ξi and ξi∗ are slack variables.
The predictive results of SVR models for silt and gravel sand soils are shown in Figure 7. It can be seen that the predictive results of the SVR model are in good agreement with the experimental results (R2 = 0.86 for silt and R2 = 0.94 for gravel sand). The SVR model exhibits remarkable prediction accuracy and most of predictive values are within the ±10% error bars, which validates their predictive availability and engineering suitability. Furthermore, it also should be noted that evaluation spots in Figure 7 is distributed uniformly throughout the thermal conductivity interval (λf ranges 0.6∼3.3 W/(m·K)), which proves the broader application range of the SVR model.
Predictive results of SVR models for (a) slit and (b) gravel sand.
3.2. Comparison of SVR Model with Ternary Fitting Model
The comparison of predictive results of the SVR model and ternary fitting model are shown in Figure 8. It can be seen that the fitting degree of the SVR model is much higher than that of the ternary fitting model. Especially for the silty clay and breccia soils, the predictive improvement effect of the SVR method is more obvious. For improvement at both high and low thermal conductivity intervals, the R2 increases to 0.84 and 0.85, respectively. It can be considered that, for the complex multiphase and porous structural characteristics of soil, the SVR model has more advantages than the traditional empirical formula model.
Comparison of predictive results of the SVR model and ternary fitting model. (a) Slit. (b) Fine sand. (c) Boulder. (d) Breccia.
The fitting results of the ternary fitting method and SVR method are statistically calculated. The distributions of R2, sample probability of mean absolute percent error (MAPE) less than 10%, and maximum relative error (δmax) are shown in Table 4. It can be seen that, compared with ternary fitting method, the minimum R2 of the SVR method is 0.84 (silty clay), and the R2 of 4 soil types is above 0.9. The 7 typical soils’ average δmax of the SVR model is 18.1% and ternary fitting model is 21.9%. Moreover, the proportion of MAPE less than 10% of SVR model are larger than that of ternary fitting model for all 7 typical soils, which further indicate that the SVR model has better performance than the ternary fitting model. The sample probability of the MAPE within 20% of the ternary fitting method is more than 98%, while the sample probability of the MAPE within 15% of the SVR method is more than 98%. Therefore, it can be considered that the ternary fitting method has a predictive accuracy of ±20% and the SVR method has a predictive accuracy of ±15%.
Statistics of prediction results of the ternary fitting model and SVR model.
Predictive model
Statistical values
Silty clay
Silt
Weathered rock
Fine sand
Gravel sand
Boulder
Breccia
Ternary fitting method
R2
0.76
0.79
0.9
0.91
0.93
0.87
0.78
MAPE ≤ 10%
63.30%
68%
76.36%
76%
76.74%
77.08%
65.38%
δmax
26.6%
22.88%
16.78%
19.62%
22.11%
19.73%
25.3%
SVR model
R2
0.84
0.86
0.91
0.92
0.94
0.94
0.85
MAPE ≤ 10%
74.74%
83.64%
72.55%
78.43%
79.07%
83.33%
69.81%
δmax
20.47%
19.89%
19.62%
17.49%
15.62%
13.23%
20.04%
4. Generalized Prediction Model for Coarse and Fine-Grained Soil
Based on the previous analysis, it can be found that some soil types, such as clay, silty sand, and gravel soil, cannot be analyzed by the nonlinear regression or SVR method for their small sample sizes. Furthermore, the overly detailed prediction models of subdivided soil types will also induce a certain degree of inconvenience in engineering application. Thus, generalized frozen soil thermal conductivity prediction models for coarse and fine-grained soil categories have been developed. The R2 and fitting coefficient values of the ternary fitting model and SVR model are listed in Table 5, and the comparison of predictive results of the above two models for coarse and fine-grained soil categories are shown in Figure 9. It shows that both the ternary fitting model and SVR model have acceptable fitting effect for two soil categories, and most evaluation spots of two models are within or near the ±10% error. The R2 of the ternary fitting model for the coarse and fine-grained soil categories is 0.80 and 0.89, respectively, while the SVR model is 0.85 and 0.91, which proves the feasibility of soil categories predictive models. Nevertheless, it also can be found that the SVR method still has a better performance despite the increase of variety and complexity of soil samples.
Fitting results of coarse and fine-grained soil categories.
Soil type
Ternary fitting
SVR
a
b
c
d
R2
R2
Fine-grained soil
−0.770
1.355
0.253
−0.102
0.80
0.85
Coarse-grained soil
−0.788
1.190
0.293
0.225
0.89
0.91
Comparison of predictive results of ternary fitting and SVR models for (a) fine-grained sand and (b) coarse-grained sand.
The MAPE distribution of two prediction models of coarse- and fine-grained soil categories is shown in Figure 10. It can be seen that the prediction accuracy of coarse-grained soil category is generally higher than that of fine-grained soil category. The MAPE of coarse-grained soil category distributes uniformly in each thermal conductivity range and its average value is lower than 10%, which means that the predictive results of two prediction models for coarse-grained soil category are more reliable in the entire frozen soil thermal conductivity distribution interval. Additionally, it should be mentioned that the large error intervals of fine-grained soil category are mainly concentrated in 0∼1.2 W/(m·K) and 2.7∼3.6 W/(m·K) ranges and the proportion of the two ranges are 6.67% and 2.6% (total of 9.27%), which is roughly excluded from its main thermal conductivity distribution range. Therefore, considering for applicability scope and convenience, the fine-grained soil category prediction models definitely have certain application values in engineering.
MAPE distribution of ternary fitting and SVR models for coarse- and fine-grained soil categories.
The prediction models of two soil categories are applied to small sample soils (clay, silty sand, and gravel soil), and the results are shown in Figure 11. It can be seen that both ternary fitting model and SVR model exhibit good predictive accuracy for all three small sample soils, which supports applicability of prediction models. Furthermore, it can be clearly noted that the SVR model has better predictive performance for clay soil, which testify that the machine learning method has broad application prospects for its capability of nonlinear and complex relationship capturing and imitating.
Prediction results of small sample soils.
5. Conclusions
In present work, a large-scale soil thermal conductivity test has been conducted by the TPS method. Correlation analysis of thermal conductivity of unfrozen and frozen soil has been adopted. Based on the measurement data of unfrozen soil thermal conductivity, the ternary fitting and SVR prediction models of frozen soil thermal conductivity for the typical soils in the QTEC are proposed considering for the essential soil properties information reflection of unfrozen soil thermal conductivity. Furthermore, to facilitate engineering applications, the prediction models of coarse and fine-grained soil categories have also been proposed and compared. The results show that
Compared with nonideal prediction accuracy of using water content and dry density as the fitting parameter, the ternary fitting model has a higher thermal conductivity prediction accuracy for typical soil types in QTEC, and more than 98% of soil specimens’ relative error are within 20%
With the capability of nonlinear and complex relationship capturing and imitating, the SVR model can further improve the frozen soil thermal conductivity prediction accuracy and more than 98% of the soil specimens’ relative error are within 15%
For coarse- and fine-grained soil categories, the above ternary fitting and SVR models still have reliable prediction accuracy and their R2 ranges from 0.8 to 0.91, which validates the applicability for small sample soils
Data Availability
All 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.
Acknowledgments
This research was supported by the National Science Foundation of China (Grant nos. 41502292 and 51574037), Applied Fundamental Research Project of China Communications Construction Co., Ltd. (nos. 2018-ZJKJ-PTJS03 and 2016-ZJKJ-02), and project funded by China Postdoctoral Science Foundation (no. 2014M560739).
DoréG.NiuF.BrooksH.Adaptation methods for transportation infrastructure built on degrading permafrost201627435236410.1002/ppp.19192-s2.0-84989314756DingY.ZhangS.ZhaoL.LiZ.KangS.Global warming weakening the inherent stability of glaciers and permafrost201964424525310.1016/j.scib.2018.12.0282-s2.0-85061656720RanY.LiX.ChengG.Climate warming over the past half century has led to thermal degradation of permafrost on the Qinghai-Tibet Plateau201812259560810.5194/tc-12-595-20182-s2.0-85042406833ChenS.LiX. F.WuT. H.Soil thermal regime alteration under experimental warming in permafrost regions of the central Tibetan Plateau202037211439710.1016/j.geoderma.2020.114397ZhangZ. Q.WuQ. B.JiangG. L.GaoS. R.ChenJ.LiuY. Z.Changes in the permafrost temperatures from 2003 to 2015 in the Qinghai-Tibet Plateau202016910290410.1016/j.coldregions.2019.102904WangT.WuT.WangP.LiR.XieC.ZouD.Spatial distribution and changes of permafrost on the Qinghai-Tibet Plateau revealed by statistical models during the period of 1980 to 20102019650166167010.1016/j.scitotenv.2018.08.3982-s2.0-85053065160NelsonF. E.AnisimovO. A.ShiklomanovN. I.Subsidence risk from thawing permafrost2001410683188989010.1038/350737462-s2.0-0035912205HjortJ.KarjalainenO.AaltoJ.Degrading permafrost puts Arctic infrastructure at risk by mid-century201891514710.1038/s41467-018-07557-42-s2.0-85058460618ShenY. J.WangY. Z.WeiX.JiaH. L.YanR. X.Investigation on meso-debonding process of the sandstone–concrete interface induced by freeze–thaw cycles using NMR technology202025211896210.1016/j.conbuildmat.2020.118962ShenY. J.HouX.YuanJ. Q. Z. H.Thermal deterioration of high-temperature granite after cooling shock: multiple-identification and damage mechanism2020711410.1007/s10064-020-01888-7YangM. X.NelsonF. E.ShiklomanovN. I.DongD. L.WanG. N.Permafrost degradation and its environmental effects on the Tibetan Plateau: a review of recent research20101031-2314410.1016/j.earscirev.2010.07.0022-s2.0-77958197920TuretskyM. R.AbbottB. W.JonesM. C.Permafrost collapse is accelerating carbon release20195697754323410.1038/d41586-019-01313-42-s2.0-85065029692SchuurE. A. G.McguireA. D.SchädelC.Climate change and the permafrost carbon feedback2015520754617117910.1038/nature143382-s2.0-84927514015YangM.WangX.PangG.WanG.LiuZ.The Tibetan Plateau cryosphere: observations and model simulations for current status and recent changes201919035336910.1016/j.earscirev.2018.12.0182-s2.0-85059664181WangQ. Z.FangJ. H.ZhaoX. Q.HuK.The influence of pavement type on the thermal stability of block-stone embankments in the warm permafrost region20202310033410.1016/j.trgeo.2020.100334ZhangM.LaiY.WuQ.A full-scale field experiment to evaluate the cooling performance of a novel composite embankment in permafrost regions2016951047105610.1016/j.ijheatmasstransfer.2015.12.0672-s2.0-84954450125LiuZ. Y.CuiF. Q.ChenJ. B.JinL.WangW.ZhangW.Study on the permafrost heat transfer mechanism and reasonable interval of separate embankment for the Qinghai-Tibet expressway202017010295210.1016/j.coldregions.2019.102952ShenY. J.YangH. W.XiJ. M.YangY.WangY. Z.WeiX.A novel shearing fracture morphology method to assess the influence of freeze–thaw actions on concrete–granite interface202016910290010.1016/j.coldregions.2019.102900ZhangY.ChenW.RiseboroughD. W.Transient projections of permafrost distribution in Canada during the 21st century under scenarios of climate change2008603-444345610.1016/j.gloplacha.2007.05.0032-s2.0-38349194801AlrtimiA.RouainiaM.HaighS.Thermal conductivity of a sandy soil201610655156010.1016/j.applthermaleng.2016.06.0122-s2.0-84991063468LiR.ZhaoL.WuT.Soil thermal conductivity and its influencing factors at the Tanggula permafrost region on the Qinghai-Tibet Plateau201926423524610.1016/j.agrformet.2018.10.0112-s2.0-85055639732DuY. Z.LiR.ZhaoL.Evaluation of 11 soil thermal conductivity schemes for the permafrost region of the central Qinghai-Tibet Plateau202019310460810.1016/j.catena.2020.104608Barry-MacaulayD.BouazzaA.SinghR. M.WangB.RanjithP. G.Thermal conductivity of soils and rocks from the Melbourne (Australia) region201316413113810.1016/j.enggeo.2013.06.0142-s2.0-84881238333CaiG.ZhangT.PuppalaA. J.LiuS.Thermal characterization and prediction model of typical soils in Nanjing area of China2015191233010.1016/j.enggeo.2015.03.0052-s2.0-84925373862ZhaoX.ZhouG.JiangX.Measurement of thermal conductivity for frozen soil at temperatures close to 0°C201914050451010.1016/j.measurement.2019.03.0692-s2.0-85064435294ZhangN.WangZ.XiaS. Q.HouX. Y.Review of soil thermal conductivity and predictive models201711717218310.1016/j.ijthermalsci.2017.03.0132-s2.0-85016399196DongY.McCartneyJ. S.LuN.Critical review of thermal conductivity models for unsaturated soils201533220722110.1007/s10706-015-9843-22-s2.0-84939993900KerstenM. S.1949Minneapolis, MN, USAUniversity of MinnesotaJohansenO.1975Trondheim, NorwayTrondheim UniversityLuS.RenT.GongY.HortonR.An improved model for predicting soil thermal conductivity from water content at room temperature200771181410.2136/sssaj2006.00412-s2.0-33846807222YanH.HeH.DyckM.A generalized model for estimating effective soil thermal conductivity based on the Kasubuchi algorithm201935322724210.1016/j.geoderma.2019.06.0312-s2.0-85068573812CôtéJ.KonradJ.-M.A generalized thermal conductivity model for soils and construction materials200542244345810.1139/t04-1062-s2.0-24044537597LuoS. Q.LuS. H.ZhangY.Soil thermal conductivity parameterization establishment and application in numerical model of central Tibetan Plateau200952433833910.3969/j.issn.0001-5733.2009.04.0082-s2.0-70049087712ZhangN.YuX. B.PradhanA.PuppalaA. J.Thermal conductivity of quartz sands by thermo-time domain reflectometry probe and model prediction201527121100401505910.1061/(asce)mt.1943-5533.00013322-s2.0-84947785853BallandV.ArpP. A.Modeling soil thermal conductivities over a wide range of conditions20054654955810.1139/s05-0072-s2.0-32144446105De VriesD. A.Thermal properties of soils1963Amsterdam, NetherlandsNorth-Holland Publishing Co.210235GoriF.A theoretical model for predicting the effective thermal conductivity of unsaturated frozen soilsProceedings of the 4th International Conference on Permafrost, FairbanksJuly 1983Fairbanks, AK, USA363368TongF.JingL.ZimmermanR. W.A fully coupled thermo-hydro-mechanical model for simulating multiphase flow, deformation and heat transfer in buffer material and rock masses201047220521710.1016/j.ijrmms.2009.11.0022-s2.0-76349100273FaroukiO. T.The thermal properties of soils in cold regions198151677510.1016/0165-232x(81)90041-02-s2.0-0019610030XuX. Z.WangJ. C.ZhangL. X.2010Beijing, ChinaScience PressHaighS. K.Thermal conductivity of sands201262761762510.1680/geot.11.p.0432-s2.0-84870993336ZhangT.WangC. J.LiuS. Y.ZhangN.ZhangT. W.Assessment of soil thermal conduction using artificial neural network models201916910290710.1016/j.coldregions.2019.102907ZhangN.ZouH.ZhangL.PuppalaA. J.LiuS.CaiG.A unified soil thermal conductivity model based on artificial neural network202015510641410.1016/j.ijthermalsci.2020.106414BangH. T.YoonS.JeonH.Application of machine learning methods to predict a thermal conductivity model for compacted bentonite202014210739510.1016/j.anucene.2020.107395BiJ.ZhangM.ChenW.LuJ.LaiY.A new model to determine the thermal conductivity of fine-grained soils201812340741710.1016/j.ijheatmasstransfer.2018.02.0352-s2.0-85042931174CuiF. Q.LiuZ. Y.ChenJ. B.DongY. H.JinL.PengH.Experimental test and prediction model of soil thermal conductivity in permafrost regions2020107247610.3390/app10072476Gangadhara RaoM.SinghD. N.A generalized relationship to estimate thermal resistivity of soils199936476777310.1139/t99-0372-s2.0-0033170208VapnikV.GolowichS. E.SmolaA. J.Support vector method for function, approximation, regression estimation and signal processing1997Cambridge, MA, USAMIT Press281287BreretonR. G.LloydG. R.Support vector machines for classification and regression2010135223026710.1039/b918972f2-s2.0-76749142536