Printer Friendly

On the prediction of the thermal conductivity of saturated clayey soils: Effect of the specific surface area.


The models for predicting thermal conductivity of soils are usually sensitive to the chosen value of the soil solid conductivity [[lambda].sub.s]. Existing approaches of estimating [[lambda].sub.s], in some cases, may lead to significant biases in predictions of overall thermal conductivities. It is postulated in the paper that the value of [[lambda].sub.s] can be estimated using the information on the microstructural, intrinsic soil property, namely the specific surface area (SSA). The new model is validated against thirty-four laboratory measurements of thermal conductivity performed on silt and clay soils. For validation purposes five widely known models, predicting overall soil thermal conductivity, are used. The analyses are performed with the use of four different values of solid thermal conductivity. The best agreement between predicted and measured conductivities is obtained when the new model is incorporated into the Johansen's method, i.e. root mean square error (RMSE) and bias are 0.172 W[m.sup.-1][K.sup.-1] and -0.040 W[m.sup.-1][K.sup.-1], respectively.


Heat transfer in solids is commonly described as the conductive flow. In some cases other mechanisms like convection or radiation also may contribute to the overall transfer. This is particularly common in fluids or gases, while in the case of soils, the conductive mechanism overwhelmingly controls (Johansen, 1975). The main parameters affecting the transfer of energy through a soil are the thermal conductivity and the specific heat (Farouki, 1981). The thermal conductivity of a soil [lambda] (being an interest of this work) can be defined as the amount of heat which is passing through a unit cross sectional area of a medium in a unit time due to a unit temperature gradient applied in the direction of this heat flow. This particular soil parameter can be determined either by in situ/laboratory measurements or by making use of existing empirical/semi-theoretical models.

In general, the methods for measuring thermal conductivity of soils can be classified into two categories, i.e. steady-state and transient heat transfer methods (Mohsenin, 1980). The steady state methods require a long time to complete the measurement. This is triggered by the fact that a soil portion under study should be in a steady state when the investigation is made. Transient methods usually require much less time and, in addition, are more versatile and easily performed (De Vries and Peck, 1958). Recently, the line source method, using a steel needle probe as a heating source, is most widely used (see Nicolas et al., 1993; Hanson et al., 2000; Rozahski and Sobotka, 2013; Lydzba et al., 2014a).

Thermal conductivity of soils (among other factors) strongly depends on mineral composition, texture, dry density [[rho].sub.d], moisture content [theta], porosity n, etc. Hence, much effort has been made to develop models incorporating these soil parameters which can be evaluated from laboratory investigations (Mickley, 1951; Gemant, 1952; Woodside and Messmer, 1961; De Vries, 1963; Johansen, 1975; Campbell, 1985; Cote and Konrad, 2005; Lu et al., 2007; Lu et al., 2014). The majority of these empirical or semi-theoretical relations need the soil solid thermal conductivity ([[lambda].sub.s]) as an input parameter. On the other hand, a proper evaluation of [[lambda].sub.s] requires cumbersome and time consuming laboratory measurements which are in contradiction with a quick prediction of the thermal conductivity. In practice, the conductivity of solids can be directly assumed to be a constant value (Farouki, 1981), it can be predicted from the quartz content of the total solids content (Johansen, 1975) or it can be estimated on the basis of the percent of clay in the soil solids (Gemant, 1952). However, as it is shown further in the paper, all existing models are sensitive to the chosen value of [[lambda].sub.s]. Therefore, a proper estimation of solid conductivity [[lambda].sub.s] is required to obtain a good prediction of the overall soil conductivity. The objective of this study is to formulate a simple model, that can be used for determination of [[lambda].sub.s], as a function of the specific surface area (SSA) which is a measure directly related to the clay fraction content, clayey mineral type as well as to the organic matter content. Furthermore, SSA is an intrinsic microstructure property which depends on both size and shape of soil particle, and various overall soil properties, used in engineering practice, are strongly affected by the magnitude of the specific surface area (e.g. Dolinar et al., 2007; Yukselen-Aksoy and Kaya, 2010; Koszela-Marek, 2014). The new model is validated against laboratory measurements of thermal conductivity of thirty-four saturated soils. For validation purposes five widely known models, predicting overall soil thermal conductivity, are used.


Thirty-four soils from three different mining regions of south-western Poland (Lower Silesia region) were used in the study. The soils were extracted from the depth ranging from approximately 77 m to 320 m below the ground level. The degree of saturation [S.sub.r] , of all investigated soils, was estimated to be in the range from 0.95 to 1. It is, therefore, assumed that the soils were fully saturated, i.e. [theta]=[[theta].sub.sat*] The soils were characterized as silt-clay or clay. The sample notations as well as basic characteristics - the contents of clay (Cl), silt (Si) and sand (Sa); specific surface area (SSA); dry bulk density ([[rho].sub.d]); soil particle density ([[rho].sub.s]); moisture content ([theta]); porosity (n) - are listed in Table 1. The water content, bulk and particle densities and soil particle size distribution were determined according to the guidelines of CEN ISO/TS 17892 standard (2004). In particular, the particle size distribution was evaluated using the aerometric and sedimentation methods. The specific surface area was measured by making use of the methylene blue adsorption method (e.g. Koszela-Marek, 2014). These measurements have been performed by extended team and the results have been archived in following reports: SPR 1 (2013), SPR 2 (2013) and SPR 1 (2014). For the purposes of this paper only a small part of all results (saturated clayey soils with negligible content of organic matter) was chosen.

The thermal conductivities were measured on intact soil cores by the needle probe test (at ambient temperature, approximately 20 [degrees]C) with the use of KD2 Pro meter with TR-1 sensor (KD2 Pro Operator's Manual, 2011). Before the measurement was conducted, first, the soil cores of length approximately 50 cm and of radius equal to 7.5 cm had been divided into three specimens (say I, II and III) of equal length, roughly 16 cm. Then, each specimen was measured three times; the TR-1 sensor (of length 10 cm) was pressed by hand and embedded parallel to the longer side of the specimen, one by one with required time shift, in three (say 1, 2, 3) different locations (see Fig. 1). Final result, overall thermal conductivity for a given soil core, was evaluated as a mean value from all performed measurements (three different locations for each specimen). The results are summarized in the last column of Table 1. It should be noted, that in the case of highly compacted soils, the sensor insertion was followed by the pilot hole bored by the drill of the diameter nearly equal to the diameter of the sensor. In this case, thermal grease was applied to reduce contact resistance errors with pre-drilled holes.

During each test a steel needle probe is heated for a time period [t.sub.h]=150 seconds (a half of the total measurement time). The temperature changes, for heating and cooling phases, are recorded and these data are fit to the following equations (in the majority of our measurements the range of temperatures was 20 [degrees]C [+ or -] 1.5 [degrees]C):

T = [m.sub.0]+[m.sub.1]t + [m.sub.3] ln t, (1)

for heating, and:

T = [m.sub.1]+[m.sub.2]t + [m.sub.3] ln-[t/[t-[t.sub.h]]], (2)

for cooling phase (KD2 Pro Operator's Manual 2011). In the equations above T is the temperature, t is the time, [m.sub.0] and [m.sub.1] is the ambient temperature for heating and cooling phase, respectively, [m.sub.2] is the rate of background temperature drift and [m.sub.3] is the slope of a line that relates temperature increase to ln t. Since (1) and (2) are log time approximations (KD2 Pro Operator's Manual 2011, Lydzba et al. 2014a) the early time data are ignored, i.e. only the final 2/3 of the data collected (during heating and cooling) are used for fitting. Finally, the thermal conductivities are evaluated using the relation:


where Q denotes the applied heat input rate. It should be emphasized that after each measurement the fitting quality was assessed by the relative error (err), which is a measure of how well the model (Eqs. (1) and (2)) fits the recorded data. According to the user's manual (KD2 Pro Operator's Manual, 2011) a good data set gives err values below 0.01. If the err value is larger than 0.01 the data has to be discarded and another reading is taken.



As mentioned in the Introduction, there exists a large number of soil thermal conductivity prediction methods. The majority of these models have some limits--other words, their usefulness is limited to a specific type of soil being under particular conditions, e.g. saturation level, texture and frozen or unfrozen state. A brief summary of existing models (used within this study) recognized by Farouli (1981) as the best models to apply to saturated, unfrozen fine soils will be presented below. Thermal conductivities are expressed in metric units, i.e. W[m.sup.-1][K.sup.-1], while the porosity n is dimensionless.


Mickley considered the unit cube of soil as the volume composed of subdomains (columns) being a solid, water and air (Mickley, 1951). For the given direction of heat flow the contributions of each subdomain are added and the thermal conductivity (with the assumption of full saturation, [S.sub.r]=1) is given by the following equation:


where [[lambda].sub.w] is the thermal conductivity of water and the length of the column l [-] (dimensionless) can be calculated, using the information on soil porosity n, i.e.:

[3l.sup.2]-[2l.sup.3]=n. (5)


Gemant considered the moist soil with the contact between soil particles as the point one (Gemant, 1952). The water was assumed to collect around these contact points forming a thermal bridge. Gemant's approach yields the following correlation:




h=9.988 * [10.sup.-3] [[rho].sub.d][theta]-[h.sub.0] [-], (8)

z = [(1-a/a).sup.2 3] [(h/2).sup.1/3] [-], (9)

[b.sup.2]=[[a/(1-a)].sup.2/3] [(h/2).sup.2/3] [-], (10)

In the equations above, [[rho].sub.d] is soil dry density [g/[cm.sub.3]], [theta] is moisture content [%], h represents the apex water (water collected around the contact points), [h.sub.0] [-] is a water adsorbed as a film around the soil particles. Adsorbed water [h.sub.0] is a function of temperature and should be recorded from the appropriate diagram. For the value of temperature 20[degrees]C, it is found that [h.sub.0]=0.01. The values of the function f{[b.sup.2]/a} have to be estimated using another diagram provided by Gemant (1952). For the purpose of this work, as well as for better efficiency of this method, it is found - by fitting the data from (Gemant 1952) - that this function can be approximated using the expression below:


3.1.3. W-M METHOD

Woodside and Messmer approach, which is based on the well-known three resistor model (originally formulated for the electrical conductivity), leads to the following correlation (Woodside and Messmer, 1961):



Within the De Vries model, saturated soil is treated as a two-phase material composed of ellipsoidal soil particles embedded in continuous medium of water. This prediction is given as (De Vries, 1963):


where F is the factor expressed by following relation:


In the equation above the g values, which were originally intended to be shape factors, are used rather as parameters to fit empirical data. Note that the g values sum to unity, i.e. [g.sub.a]+[g.sub.b]+[g.sub.c]=1.


Johansen expressed the thermal conductivity of an unsaturated soil as a function of both thermal conductivity at saturated state ([[lambda].sub.sat]) and the one at dry state ([[lambda].sub.dry]) (Johansen, 1975). The concept of normalized thermal conductivity [K.sub.e] (called Kersten's number) was introduced:

[K.sub.e]=([lambda]-[[lambda].sub.dry])/([[lambda].sub.sat]-[[lambda].sub.dry]) [-]. (15)

A semi-empirical relation, based on the dry density, was formulated to estimate [[lambda].sub.dry], while the thermal conductivity at saturated state [[lambda].sub.sat] was proposed to be estimated using the geometric mean, "average" of the parallel and series thermal conductivities, i.e.:


Kersten's number ([K.sub.e] ), for fine textured soils, was related to the degree of saturation [S.sub.r] by the following formula:

[K.sub.e] = log [S.sub.r]+1.0 [-]. (17)

Therefore, it is evident that in the case of fully saturated soils ([S.sub.r]=1), Kersten's number is also taking the value of 1, and using Eqs. (15) and (16), we obtain the formula for the thermal conductivity at saturation, i.e.[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]


Each correlation, briefly described above, needs the value of the thermal conductivity of soil solid [[lambda].sub.s] as an input parameter. In the literature, different approaches for evaluation of [[lambda].sub.s] are proposed. Farouki (1981), considering results of different researchers, pointed out that thermal conductivity of soil solid, in the case of silt-clay soils, can be directly assumed as [[lambda].sub.s]=2.0-3.0 W[m.sup.-1][K.sup.-1].

Johansen (1975) determined solid conductivity using the geometric mean equation from the quartz content q as a fraction of the total solids content and thermal conductivities of quartz [[lambda].sub.q] (7.7 W[m.sup.-1][K.sup.-1]) and other minerals [[lambda].sub.o] (2.0 W[m.sup.-1][K.sup.-1] for soils with q>0.2, and 3.0 W[m.sup.-1][K.sup.-1] for q[less than or equal to]0.2), i.e.: [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (18)

The estimation of solid conductivity using Eq. (18) requires, therefore, the measurement of quartz content. This test is not commonly performed. In practice, it is often assumed that quartz content q is equal to the content of sand (Peters-Lidard et al., 1998). However, it is obvious that in a case of fine soils (being the study of the paper), i.e. when the content of sand fraction is, in general, very low, the geometric mean does not deviate from constant value 3.0 W[m.sup.-1][K.sup.-1] (being the conductivity of other minerals) by more than a few percent.

Gemant (1952) suggested that thermal conductivity of soil solid can be estimated using the information on the percent of clay in the soil solids Cl:

[[lambda].sub.s]=5.84-0.033(Cl). (19)

The above equation can lead to underestimation or overestimation of solid conductivity. At one extreme case, i.e. when there is no clay, Farouki (1981) pointed out that Eq. (19) gives 5.84 W[m.sup.-1][K.sup.-1] which is too low for quartz sand by about 25 %. On the contrary, in the case of 100% of clay, Eq. (19) gives 2.54 W[m.sup.-1][K.sup.-1]--this value is greater than the minimum one suggested by Farouki (1981) for silt-clay soils, namely 2.0 W[m.sup.-1][K.sup.-1].


Soil mineral composition is an important factor affecting thermal conductivity of the soil solid due to the substantial diversity of thermal conductivity of particular minerals being its constituents (Farouki, 1981). For instance, the thermal conductivity of quartz minerals is approximately 2-3 times higher than the clay ones, e.g. 7.7 W[m.sup.-1][K.sup.-1] for quartz (Clauser and Huenges, 1995) or 2.9 W[m.sup.-1][K.sup.-1] for clay minerals (Bristow, 2002). Furthermore, the thermal conductivities of an organic matter (0.25 W[m.sup.-1][K.sup.-1]) and soil minerals differ by one order of magnitude. Hence, the recognition of soil mineral composition plays a vital role in proper evaluation of solid conductivity. In other words, if the mineral composition of soil is well recognized and, in addition, thermal conductivities of particular minerals are given, the conductivity of soil solid can be quite accurately estimated.

In practice, however, the investigation of soil mineral fabric is not commonly conducted. Therefore, there is a need for searching such measures that can directly exhibit the mineral composition of soil solid. It was, in some sense, the idea of Gemant (Eq. (19)) who assumed that as the content of clay fraction is increasing--and as a consequence, the contribution of clay minerals in the solid is also increasing--the conductivity of soil solid should decrease. As mentioned earlier, this approach has some disadvantages at extreme values of clay content.

Hence, for the purpose of solid thermal conductivity estimation, it is proposed to use another microstructural measure, namely the specific surface area (SSA), which is directly related to the clay fraction content, clayey mineral type as well as to the organic part content. The SSA of soil material can variously be defined (Hillel, 1998) but a reasonable definition is that SSA is the total surface area of particles per unit mass, expressed in terms of square meters per gram. Specific surface area is an intrinsic microstructure property which depends on both size and shape of soil particle. It is evident that flattened or elongated particles exhibit greater surface per unit mass or volume than spherical particles do. Although the clay minerals have small sizes, their shape is--in general--platy (Fig. 2) and therefore they mainly contribute to the overall value of SSA. For example, the SSA of sand may be no more than 1-2 [m.sub.2]/g, while that of clay may be at least 100 times higher (Hillel, 1998; Mitchell and Soga, 2005). In that sense, SSA appears to be a measure which reveals the mineral composition of soil solid.

According to the proposition of Gemant (Eq. (19)), the solid conductivity is decreasing as the content of clay fraction is increasing. As a consequence, the overall thermal conductivity of soil is also expected to be a decreasing function of Cl. Measured thermal conductivities, for all investigated soils, are plotted in Figure 3a against the clay content. The line of trend visualizes the decreasing thermal conductivity as Cl increases. The same behaviour is observed in the second graph (Fig. 3b) where measured thermal conductivity is plotted against the specific surface area. Thus, due to the fact that SSA (as well as volume fraction of clay particles) is a micro-structural property and the overall soil conductivity depends on its value, it is postulated that it is the solid conductivity [[lambda].sub.s] that is affected by the SSA (similarly as in the case of the proposition of Gemant (1952) where solid conductivity is a function of Cl). It has also been found in the paper of Lydzba et al. (2014b) where micromechanical approaches, namely Mori-Tanaka and Self-Consistent schemes are used to simulate the microstructure of clays.

Since high values of SSA exhibit large contribution of clay minerals in the solid, and bearing in mind the relative difference in conductivities of quartz and clay minerals, it is evident that the value of solid conductivity [[lambda].sub.s] is a decreasing function of specific surface area. The relation between [[lambda].sub.s] and the SSA is postulated to be characterized by the modified power function:


where [a.sup.1], [a.sup.2], [a.sup.3] are the parameters to be fit. Note that the choice of the above function is not arbitrary since Eq. (22) should hold the following conditions. As the SSA value tends to zero (no content of clay minerals) the conductivity should be close to that of quartz sand (the value of pure quartz mineral, 7 7 W[m.sup.-1][K.sup.-1], is prescribed). On the other hand for , the large value of SSA, e.g. 800 [m.sub.2]/g (such extreme value of SSA is reported for the minerals of the smectite group, e.g. montmorillonite (Mitchell and Soga, 2005; Koszela-Marek, 2014) the minimum value of solid conductivity, found in the literature (Farouki, 1981), namely 2.0 W[m.sup.-1][K.sup.-1] is assumed. Furthermore, as mentioned earlier, the typical range of SSA for clay soils is 100 , to 200 [m.sub.2]/g (the average value is 150 [m.sub.2]/g). Thus, taking this into account a third condition is defined: we assign the value of conductivity 3.0 W[m.sup.-1][K.sup.-1] (which is typical for silt-clay soils) to the value of SSA equal to 150 [m.sub.2]/g. Using the above conditions, a set of points, (SSA, [[lambda].sub.s]), i.e. {(0, 7.7);(150, 3.0);(800, 2.0)} creates the data which is used for fitting Eq. (22). The least square method is used and the following values of parameters [a.sub.1], [a.sub.2] and [a.sub.3] are found: [a.sub.1]= 5.7, [a.sub.2]=0.988, [a.sub.3]= 2.0. The data as well as the fitting function (Eq. (22)) is presented in Figure 4. It has to be emphasized that the proposed model, Eq. (22), is verified for all thirty-four soils, for which the range of SSA is: 55.8-171.7 [m.sub.2]/g (see Table 1). The range, in which the model is tested, is also indicated in Figure 4 by the continuous line and filled area under the curve.



The new model--Eq. (22) - was used to predict the solid conductivities of all investigated soils. The solid conductivity, estimated by the new model, is in the range from 2.72 to 4.91 W[m.sup.-1][K.sup.-1]. In Figures 5a and 5b estimations of the new model are plotted against the conductivities obtained with the use of geometric mean (Eq. (18)) and Gemant's method (Eq. (19)), respectively. Furthermore, for ease of comparison, 1:1 line is also provided.

It can be noticed that the results of solid conductivity obtained with the use of the new model are, in general, greater than the ones resulting from the geometric mean. Moreover, as it was expected, the geometric mean estimation, due to the low content of sand fraction, is close to the constant value 3.0 W[m.sup.-1][K.sup.-1], for all considered soils. The new model gives slightly lower values of solid conductivity, when compared to the Gemant's approach (only five of thirty-four results are above 1:1 line).


All models described in Section 3.1 were used to predict the overall soil conductivities. The results were obtained using four different values of solid conductivities, i.e. constant value according to suggestion of Farouki (1981), namely [[lambda].sub.s]=3.0 W[m.sup.-1][K.sup.-1], the one estimated as the geometric mean (Eq. (18)), the one determined with the use of Gemant's approach (Eq. (19)) and the value of conductivity resulting from the new model (Eq. (22)). The conductivity of water phase was assumed as [[lambda].sub.w]=0.585 W[m.sup.-1][K.sup.-1]. Predicted conductivities were compared with the laboratory results using the root mean square error (RMSE):


and the bias:


where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] and [[lambda].sub.predicted] represent measured and predicted values of thermal conductivity, respectively, and m is the number of measurements. The RMSE (Eq. (23)) and the bias (Eq. (24)) were calculated for all considered cases (four different values of [[lambda].sub.s]) and they are listed in Tables 2 and 3, respectively. The results are also graphically presented in Figure 6, where predicted thermal conductivities are plotted against the measured ones. For better visualization of results the linear trend as well as 1:1 line is also displayed. In addition, the values of the slope of the trend line are provided in Table 4 (note that as this value is equal to 1, the trend and 1:1 lines possesses the same slope). For clarity of presentation, in Figure 6 the results corresponding to the constant value of solid conductivity ([[lambda].sub.s]=3.0 W[m.sup.-1][K.sup.-1]) are not displayed--this is due to fact that these results almost coincide with the ones obtained using the geometric mean (Eq. (18)).

Performed analyses confirmed that the proper estimation of overall thermal conductivity is strongly dependent on the appropriate value of the solid conductivity that is used, as an input parameter, within the particular model, e.g. the W-M model gives a relatively good prediction (low value of RMSE=0.243 W[m.sup.-1][K.sup.-1] and bias=-0.005 W[m.sup.-1][K.sup.-1]) in the case of the constant value of [[lambda].sub.s]=3.0 W[m.sup.-1][K.sup.-1] while it becomes the worst prediction when the [[lambda].sub.s] determined according to Eq. (19) is used - the greatest values of RMSE=0.619 W[m.sup.-1][K.sup.-1] and bias= 0.580 W[m.sup.-1][K.sup.-1].

In general, the predicted thermal conductivities rather underestimate the measured ones when the calculations are performed with solid conductivity estimated as the geometric mean (Eq. (18)) and overestimate them when [[lambda].sub.s], evaluated by Gemant's expression (Eq. (19)), is used. The results obtained with the solid conductivity as a constant value or with [[lambda].sub.s] as the geometric mean almost coincide. Furthermore, when [[lambda].sub.s] is a constant value or it is estimated as the geometric mean, the slope of the trend line is significantly different than 1:1 line. This fact is the most evident in the case of W-M method--one can see (Fig. 6) that in spite the RMSE and bias are relatively small, some results strongly overestimate while others underestimate the measured conductivities (the dividing point is around [lambda]2.0 W[m.sup.-1][K.sup.-1])--the slope is 0.301 and 0.298 for [[lambda].sub.s]=3.0 W[m.sup.-1][K.sup.-1] and [[lambda].sub.s] determined as the geometric mean, respectively.

When the new model is not taken into consideration, the best agreement between predicted and measured thermal conductivities is obtained in the case of Gemant's method with the value of [[lambda].sub.s] estimated according to Eq. (19) --RMSE=0.218 W[m.sup.-1][K.sup.-1] and bias=0.006 W[m.sup.-1][K.sup.-1]. Nevertheless, this method is the most cumbersome one among all considered models. This is due to the fact that within this model the variation of adsorbed water ([h.sub.0]) with temperature as well as the variation of the function f{[b.sub.2]/a} with the soil dry density and moisture content (see Chapter 3.2) has to be estimated using diagrams provided by Gemant (1952). Furthermore, the value of the slope, 0.657, is still not satisfactory--some results are underestimated while others are overestimated.

Focusing now on the results obtained with the use of the new model (Eq. (22)), the following conclusions can be drawn. Considering RMSE and bias only, the new model gives better predictions for all methods, except four cases, namely: Gemant's method with [[lambda].sub.s] estimated by his own proposition (Eq. (18)), W-M model with [[lambda].sub.s]=3.0 W[m.sup.-1][K.sup.-1] or with [[lambda].sub.s] estimated as the geometric mean and De Vries model with solid conductivity evaluated as the geometric mean. On the other hand, it is most remarkable that when one considers the slope of the trend line it can be seen that the predictions obtained by the new model are significantly improved, i.e. for all models when [[lambda].sub.s]=3.0 W[m.sup.-1][K.sup.-1] the slope is in the range from 0.301 to 0.358, [[lambda].sub.s] evaluated according to Gemant proposition gives the slope from 0.640 to 0.707, in the case of geometric mean it is from 0.298 to 0.362, while using the new model the slope is in the range from 0.798 to 0.874.

Finally, the best prediction is obtained when the new model is used as an input parameter within the Johansen's method, i.e. RMSE=0.172 W[m.sup.-1][K.sup.-1], bias=-0.040 W[m.sup.-1][K.sup.-1] and the slope is 0.855.


It was shown that all correlations (existing in the literature), used to predict thermal conductivities of saturated clayey soils, are sensitive to the chosen value of the soil solid conductivity [[lambda].sub.s]. Laboratory measurement of the solid conductivity is a difficult task. As a consequence, different researchers have attempted to find the model for proper estimation of [[lambda].sub.s]. Nevertheless, the analyses performed in this work showed that using either one constant value of solid conductivity [[lambda].sub.s]=3.0 W[m.sup.-1][K.sup.-1] (as suggested by Farouki (1981)), estimating it as the geometric mean (Eq. (18)) or evaluating on the basis of clay content (Gemant's approach--Eq. (19)) may lead to significant biases in predictions of overall thermal conductivities.

It is obvious that soil solid conductivity is affected by its mineral composition due to the substantial diversity of thermal conductivity of particular minerals being its constituents and therefore the appropriate recognition of soil mineral composition plays a vital role in proper evaluation of solid conductivity. It was postulated, in the paper, that this can be performed, indirectly, by the investigation of soil specific surface area, i.e. as the value of SSA is increasing, the contribution of clay minerals (in soil solid) is also increasing, and as a consequence the thermal conductivity of solid should decrease. The relation between [[lambda].sub.s] and the SSA was postulated to be characterized by the modified power function--Eq. (22). The new model was validated against laboratory data, i.e. thirty-four soil samples (range of SSA: 55.8-171.7 [m.sub.2]/g) were used. The predictions of overall thermal conductivity were obtained using five widely known correlations, namely: Mickley model, Gemant model, W-M model, De Vries model and Johansen model. The best prediction was obtained when the new model was incorporated into the Johansen's expression and therefore this correlation is suggested to be used when the thermal conductivity of saturated clayey soils is estimated.


Bristow, K.L.: 2002, Thermal conductivity. Methods of Soil Analysis, Soil Science Society of America, Inc., Madison, 1209-1226.

Campbell, G.S.: 1985, Soil physics with BASIC: Transport models for soil-plant systems. Elsevier Sci. Publ. Co., New York.

CEN ISO/TS 17892: 2004, Geotechnical investigation and testing - Laboratory testing of soil. International Organization for Standardization, Geneva, Switzerland.

Clauser, C. and Huenges, E.: 1995, Thermal conductivity of rocks and minerals. Rock Physics & Phase Relations: A Handbook of Physical Constants (ed. T.J. Ahrens), American Geophysical Union, Washington, D. C. DOI: 10.1029/RF003p0105

Cote, J. and Konrad, J.M.: 2005, A generalized thermal conductivity model for soils and construction materials. Can. Geotech. J., 42, 443-458. DOI: 10.1139/cgj-2013-0388

De Vries, D.A. and Peck, A.J: 1958, On the cylindrical probe method of measuring thermal conductivity with special reference to soils. Austr. Journal of Physics, 11.

De Vries, D.A.: 1963, Thermal properties of soils. W.R. Van Wijk (ed.) Physics of plant environment, North-Holland Public Co., Amsterdam.

Dolinar, B., Misic, M. and Trauner, L.: 2007, Correlation between surface area and Atterberg limits of finegrained soils. Clays and Clay Miner., 55, 519-523. DOI: 10.1346/CCMN.2007.0550506

Farouki, O.T.: 1981, Thermal properties of soils. CRREL Monograph, 81-1, New Hampshire.

Gemant, A.: 1952, How to compute thermal soil conductivities. Heating, Piping and Air Conditioning, 24 (1), 122-123.

Hanson, J.L., Edil, T.B. and Yesiller, N.: 2000, Thermal properties of high water content materials. Geotechnics of High Water Content Materials, T.B. Edil and P.J. Fox (eds.), West Conshohocken, PA, 137-151.

Hillel, D.: 1998, Environmental soil physics: Fundamentals, applications, and environmental considerations. Academic press.

Johansen, O.: 1975, Thermal conductivity of soils. Ph.D. diss., Norwegian Univ. of Science and Technology, Trondheim.

KD2 Pro Thermal Properties Analyzer Operator's Manual: 2011, Version 10, Decagon Devices, Inc.

Kersten, M.S.: 1949, Laboratory research for the determination of the thermal properties of soils. ACFEL Tech. Tech. Rep., 23, Univ. of Minnesota, Minneapolis.

Koszela-Marek, E.: 2014, The effect of 200 MPa pressure on specific surface area of clay. Studia geotechnica et mechanica, 36 (4), 3-6. DOI: 10.2478/sgem-2014-0030

Lu, S., Ren, T., Gong, Y. and Horton, R: 2007, An improved model for predicting soil thermal conductivity from water content at room temperature. Soil Science Society of America Journal, 71 (1), 8-14. DOI: 10.2136/sssaj2006.0041

Lu, Y., Lu, S. Horton, R. and Ren, T.: 2014, An empirical model for estimating soil thermal conductivity from texture, water content, and bulk density. Soil Science Society of America Journal, 78 (6), 1859-1868. DOI: 10.2136/sssaj2014.05.0218

Lydzba, D., Rozanski, A., Rajczakowska, M. and Stefaniuk, D.: 2014a, Efficiency of the needle probe test for evaluation of thermal conductivity of composite materials: two-scale analysis. Studia geotechnica et mechanica, 36 (1), 55-62. DOI: 10.2478/sgem-2014-0007

Lydzba, D., Rajczakowska, M., Rozanski, A. and Stefaniuk, D.: 2014b, Influence of the moisture content and temperature on the thermal properties of soils: laboratory investigation and theoretical analysis. Procedia Engineering, 91, 298-303. DOI: 10.1016/j.proeng.2014.12.064

Mickley, A.S.: 1951, The thermal conductivity of moist soil. American Institute of Electrical Engineers Transactions, 70, 1789-1797.

Mitchell, J. and Soga, K.: 2005, Fundamentals of soil behavior. John Wiley and Sons, Inc., Hoboken, NJ.

Mohsenin, N.N.: 1980, Thermal properties of foods and agricultural materials. Gordon and Breach, New York.

Nicolas, J., Andr'e, P., Rivez, J.F. and Debbaut, V.: 1993, Thermal conductivity measurements in soil using an instrument based on the cylindrical probe method. Review of Scientific Instruments, 64, No. 3, 774-780. DOI: 10.1063/1.1144158

Peters-Lidard, C.D., Blackburn, E., Liang, X. and Wood, E.F.: 1998, The effect of soil thermal conductivity parameterization on surface energy fluxes and temperatures. Journal of the Atmospheric Sciences, 55(7), 1209-1224. DOI: 10.1175/1520-0469(1998)055<1209:TEOSTC>2.0.CO;2

Report SPR 1: 2013, Report from the measurements of mechanical and thermal properties of soils. Wroclaw University of Technology, I10/2013/S-001.

Report SPR 2: 2013, Report from the measurements of mechanical and thermal properties of soils. Wroclaw University of Technology, I10/2013/S-002.

Report SPR 1: 2014, Report from the measurements of mechanical and thermal properties of soils. Wroclaw University of Technology, I10/2014/S-001.

Rozanski, A. and Sobotka, M.: 2013, On the interpretation of the needle probe test results: thermal conductivity measurement of clayey soils. Studia Geotechnica et Mechanica, 35(1), 195-207, DOI: 10.2478/sgem-2013-0015

Woodside, W. and Messmer, J.M.: 1961, Thermal conductivity of porous media. Journal of Applied Physics, 32 (9), 1688-1706. DOI: 10.1063/1.1728420

Yukselen-Aksoy, Y. and Kaya, A.: 2010, Method dependency of relationships between specific surface area and soil physicochemical properties. Applied Clay Science, 50, 182-190.

Adrian ROZANSKI (*) and Damian STEFANIUK

Cite this article as: Rozanski A, Stefaniuk D: On the prediction of the thermal conductivity of saturated clayey soils: Effect of the specific surface area. Acta Geodyn. Geomater., 13, No. 4 (184), 339-349, 2016. DOI: 10.13168/AGG.2016.0016

Wroclaw University of Science and Technology, Faculty of Civil Engineering ul. Wybrzeze Wyspiahskiego 27, 50-370 Wroclaw, Poland

(*) Corresponding author's e-mail:


Article history:

Received 30 October 2015

Accepted 4 May 2016

Available online 6 June 2016




Thermal conductivity

Specific surface area Modelling

Table 1 Notations, basic characteristics and thermal conductivities
of soils (Report SPR 1 2013, Report SPR 2 2013 and Report SPR 1 2014).

Sample  Cl  Si  Sa      SSA      [[rho].sub.d]  [[rho].sub.s]
        %   %    %  [m.sub.2]/g  g/[cm.sub.3]   g/[cm.sub.3]

S_01    36  58   6     93.0          1.66           2.61
S_02    43  45  12    101.7          1.97           2.58
S_03    39  49  12    127.3          1.88           2.65
S_04    54  46   0    148.5          1.61           2.56
S_05    36  58   6    103.9          1.52           2.63
S_06    70  24   6     90.7          1.62           2.61
S_07    58  37   5    102.3          1.73           2.67
S_08    41  57   2     88.5          1.80           2.62
S_09    46  50   4    121.9          1.97           2.68
S_10     7  92   1     59.6          1.51           2.63
S_11    41  32  27     57.0          1.75           2.63
S_12    58  37   5    129.3          1.51           2.60
S_13    36  28  36     90.4          1.93           2.63
S_14    54  42   4    132.6          1.94           2.65
S_15    42  57   1     84.1          1.32           2.61
S_16    60  33   7    102.3          1.68           2.61
S_17    39  40  21     64.9          1.85           2.60
S_18    47  48   5    148.1          1.25           2.24
S_19    67  32   1    171.7          1.67           2.64
S_20    44  52   4    109.0          1.84           2.63
S_21    31  69   0    119.0          1.53           2.63
S_22    35  65   0     63.2          1.58           2.42
S_23    30  63   7     95.2          1.86           2.65
S_24    31  65   4     78.8          1.74           2.59
S_25    66  33   1    147.1          1.59           2.61
S_26    64  35   1    146.0          1.58           2.67
S_27    46  54   1     86.8          1.72           2.63
S_28    32  67   1    118.3          1.67           2.63
S_29    58  28  15    159.7          1.67           2.65
S_30    24  75   1     73.2          1.67           2.68
S_31    40  59   1    112.4          1.64           2.63
S_32    13  87   0     55.8          1.63           2.62
S_33    31  69   0     90.4          1.77           2.62
S_34    71  24   5     76.9          1.82           2.64

Sample  [theta]=[[theta].sub.sat]  n           measured
                   %               -      W[m.sup.-1][K.sup.-1]

S_01             21.9             0.364           1.897
S_02             11.9             0.235           2.505
S_03             15.6             0.292           1.904
S_04             22.9             0.370           1.642
S_05             27.9             0.423           1.494
S_06             23.4             0.379           1.780
S_07             20.3             0.351           1.943
S_08             17.4             0.313           2.230
S_09             13.4             0.264           2.103
S_10             28.2             0.426           1.960
S_11             19.2             0.335           2.325
S_12             27.6             0.418           1.721
S_13             13.9             0.268           2.611
S_14             13.8             0.267           2.098
S_15             37.5             0.495           1.425
S_16             21.1             0.355           2.147
S_17             15.5             0.287           2.553
S_18             35.2             0.441           1.439
S_19             21.9             0.366           1.426
S_20             16.3             0.300           2.004
S_21             27.2             0.417           1.916
S_22             22.0             0.347           1.871
S_23             16.1             0.299           2.042
S_24             18.8             0.328           1.779
S_25             24.6             0.391           1.511
S_26             25.7             0.406           1.404
S_27             20.1             0.345           1.817
S_28             21.9             0.365           2.083
S_29             22.3             0.372           1.505
S_30             22.5             0.376           1.895
S_31             23.1             0.378           1.603
S_32             23.0             0.376           2.190
S_33             18.2             0.323           1.926
S_34             17.0             0.310           2.136

Table 2 Root mean square error (RMSE) of the Mickley model,
Gemant model, W-M model, De Vries model and Johansen model
corresponding to four different values of solid conductivity

 Method                           RMSE [W[m.sup.-1][K.sup.-1]]
          [[lambda].sub.s] = 3.0   [[lambda].sub.s]--Eq. (18)

Mickley           0.296                      0.290
 Gemant           0.504                      0.483
  W-M             0.243                      0.267
De Vries          0.275                      0.274
Johansen          0.317                      0.308

          [[lambda].sub.s]--Eq. (19)  [[lambda].sub.s]--New model

Mickley             0.435                        0.217
 Gemant             0.218                        0.258
  W-M               0.619                        0.368
De Vries            0.401                        0.207
Johansen            0.332                        0.172

Table 3 Bias of the Mickley model, Gemant model, W-M model, De
Vries model and Johansen model corresponding to four different
values of solid conductivity [[lambda].sub.s].

 Method                           Bias [W[m.sup.-1][K.sup.-1]]
          [[lambda].sub.s] = 3.0   [[lambda].sub.s]--Eq. (18)

Mickley           0.176                      0.135
 Gemant           0.441                      0.407
  W-M             0.005                     -0.039
De Vries          0.141                      0.103
Johansen          0.215                      0.179

          [[lambda].sub.s]--Eq. (19) [[lambda].sub.s]--New model

Mickley             -0.378                     -0.128
 Gemant              0.006                      0.200
  W-M               -0.580                     -0.319
De Vries            -0.344                     -0.131
Johansen            -0.261                     -0.040

Table 4 The slope of the trend line of the Mickley model, Gemant
model, W-M model, De Vries model and Johansen model corresponding
to four different values of solid conductivity [[lambda].sub.s].

          [[lambda].sub.s] = 3.0          Slope [-]
          W[m.sup.-1][K.sup.-1]   [[lambda].sub.s]--Eq. (18)

Mickley           0.328               0.329
 Gemant           0.358               0.362
  W-M             0.301               0.298
De Vries          0.342               0.345
Johansen          0.357               0.362


         [[lambda].sub.s]--Eq. (19)  [[lambda].sub.s]--New model

Mickley            0.679                        0.874
 Gemant            0.657                        0.798
  W-M              0.640                        0.865
De Vries           0.676                        0.841
Johansen           0.707                        0.855
COPYRIGHT 2016 Akademie Ved Ceske Republiky, Ustav Struktury a Mechaniky Hornin
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2016 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Title Annotation:ORIGINAL PAPER
Author:Rozanski, Adrian; Stefaniuk, Damian
Publication:Acta Geodynamica et Geromaterialia
Article Type:Report
Date:Oct 1, 2016
Previous Article:Characterizing gas permeability and pore properties of Czech granitic rocks.
Next Article:Geodynamic studies in the Pieniny Klippen belt in 2004-2015.

Terms of use | Privacy policy | Copyright © 2022 Farlex, Inc. | Feedback | For webmasters |