A Practicable Guideline for Predicting the Thermal Conductivity of Unconsolidated Soils

: For large infrastructure projects, such as high-voltage underground cables or for evaluating the very shallow geothermal potential (vSGP) of small-scale horizontal geothermal systems, large-scale geothermal collector systems (LSCs), and fifth generation low temperature district heating and cooling networks (5GDHC), the thermal conductivity ( λ ) of the subsurface is a decisive soil parameter in terms of dimensioning and design. In the planning phase, when direct measurements of the thermal conductivity are not yet available or possible, λ must therefore often be estimated. Various empirical literature models can be used for this purpose, based on the knowledge of bulk density, moisture content, and grain size distribution. In this study, selected models were validated using 59 series of thermal conductivity measurements performed on soil samples taken from different sites in Germany. By considering different soil texture and moisture categories, a practicable guideline in the form of a decision tree, employed by empirical models to calculate the thermal conductivity of unconsolidated soils, was developed. The Hu et al. (2001) model showed the smallest deviations from the measured values for clayey and silty soils, with an RMSE value of 0.20 W/(m · K). The Markert et al. (2017) model was determined to be the best-fitting model for sandy soils, with an RMSE value of 0.29 W/(m · K).


Introduction
The decarbonization of the energy supply in the heating sector is necessary to meet climate targets [1].Shallow geothermal systems can contribute significantly to achieving these goals.Particularly for very shallow closed-loop systems, as well as regarding heat extraction rates [2], the very shallow geothermal potential (vSGP) in the form of the thermal conductivity of soils (λ) is a decisive factor for efficient performance [3].
In addition to the generation of sustainable geothermal heating and cooling energy for individual residential and industrial buildings, there is also the possibility of using large-scale geothermal collector systems (LSCs) to cover significantly higher energy demands [4,5].In combination with fifth generation low-temperature district heating and cooling networks (5GDHC), these systems can meet the needs of entire residential areas [5][6][7][8].Due to the low temperature level and the non-insulated pipes [5] used in these systems, these networks also act as ground heat exchangers; thus, thermal conductivity is also a decisive factor for the energetic performance here.
Another key element of the energy transformation, particularly in regards to the SuedLink project in Germany [9], is the distribution of electricity from the site of production to the consumer over long distances via high-voltage direct current underground cables.The basic requirement for this transformation is that the function of the cables Soil Syst.2024, 8, 47 2 of 22 is guaranteed.For this reason, the thermal conductivity of the surrounding subsurface is of major relevance [10], as the heating of the cables due to the lack of heat dissipation determines their durability [11].For planning such large infrastructure projects, the thermal conductivity of the soil must therefore be considered when estimating the heat dissipation in the cable environment.
The thermal conductivity of unconsolidated soils is mainly influenced by the grain size distribution, bulk density, and moisture content of the soil physical parameters [2,[12][13][14].Direct measurements of λ are often performed in situ or in the laboratory using single-needle sensors [15][16][17][18][19].For example, Markert et al. [20] proposed an approach for measuring λ in a wide moisture content spectrum based on the evaporation method.
Thermal conductivity measurements in the planning phases of large-scale infrastructure projects such as LSCs, 5HDHC, or high-voltage underground cables are often not possible due to the high effort required or the early stage of the project.However, λ can be calculated using a variety of largely empirical models, if specific soil input parameters are known or assumed [21][22][23][24][25][26][27][28][29][30].This requires the estimation of the accuracy or error of the model used.For this reason, the aim of this study is to validate selected thermal conductivity models using a variety of laboratory measurement series.In addition, by consideration of subsets with different soil properties regarding grain size distribution or moisture content, a practicable guideline for the use of empirical models for calculating the thermal conductivity of soils can be developed, based on the knowledge of the grain size distribution, bulk density, and moisture content.
The soil properties differ from those of previous validations [31][32][33][34] due to different sampling sites, all in Germany, and the combination of two measurement methods using needle sensors.The selection of empirical thermal conductivity models for validation was based on the work of Wessolek et al. [34], with additional consideration of the model of Markert et al. [27].Therefore, the selection used in this study also differs from that of Dong et al. [31] and Tarnawski et al. [33].Moreover, while He et al. [32] only focused on models based on the model of Johansen [35,36], this study also validated models with different basic approaches.
Tarnawski et al. [33] only considered the measured thermal conductivities of soils in dry and saturated conditions, whereas in this study, a wide range of saturations was considered.Furthermore, the detailed consideration of subsets, in terms of moisture and grain size distribution, employing this amount of data and a wider range of model selection in terms of publication periods, was not performed in the previous validations and is therefore novel.

Thermal Conductivity Measurement Database Used for Validation
For the validation, a total of 59 measurement series were used.The data are taken from various research projects in Germany of the Shallow Geothermal Energy Working Group from the Friedrich-Alexander-University Erlangen-Nuremberg.A measurement series is defined as a certain number of measurements (measurement points) of the thermal conductivity at specific moisture contents on a soil sample with a defined bulk density, grain size distribution, and measurement method.
The thermal conductivity measurements were performed exclusively with a singleneedle sensor.The devices used and their measurement accuracy are listed in Table 1.The configurations of the used single-needle sensors comply with the specifications according to IEEE Std 442-1981 [37] and ASTM D5334 [38].
The measurements were carried out using two different methods.These are described in Table 2 and in Figure 1.The measurements were carried out using two different methods.These are described in Table 2 and in Figure 1.
Table 2. Measurement methods used for the thermal conductivity measurements.For each thermal conductivity measurement, the moisture content was determined.Here, for the number of measurement points per measurement series, the complete water content range is considered.

49
• Evaporation method based on the approach of Markert et al. [20], with continuous measurements, obtained at 15 min to 60 min intervals, or temporary measurements, e.g., performed at daily intervals; measurements performed in cylinders.
Figure 1.Schematic illustration of the (a) evaporation method and (b) the punctual method of thermal conductivity measurement used for validation, taken from the work of Rammler et al. [18], with ω 1/2 = medium water contents, ω dry = water content after oven drying, and ω sat = water content after capillary saturation.Also shown are three exemplary measurement series, which are distinguished by their different colors, for measurements (c) with the punctual method, (d) evaporation method with temporary measurements, and (e) evaporation method with continuous measurements (measuring points cannot be recognized as points due to the high data density).For the evaporation method, oven drying and the measurement of thermal conductivity at a moisture content θ = 0% were not carried out for each sample.In the case of high implausible deviations of λ within one measurement series, the deviating measurement points were removed manually.
All soil samples were installed in the cylinders with volumes of 397 cm 3 or 944 cm 3 with a specific bulk density using undisturbed sampling.Possible changes in bulk density due to shrinkage and swelling processes during the evaporation method [18,39] were neglected for validation.These changes in density are responsible for the fact that the maximum thermal conductivity in the evaporation method is sometimes not reached at saturation, but rather at an intermediate moisture content.In the case of slightly fluctuating bulk densities in measurement series measured using the punctual method, the average bulk density was assumed for validation.The individual bulk densities deviate from the mean value by a maximum of 0.04 g/cm 3 .
Due to the different methods or measurement intervals used, the number of measurement points per measurement series varies.However, for the validation, each measurement series was considered equally for authentication, independent of the number of measurement points.
As the measurement series originate from various research projects, the soil samples were taken from different locations in Germany; therefore, a wide range of soil texture classes was covered for the validation (Figure 2).Therefore, only the USDA soil texture classes sandy clay and silt were not considered for the validation.Very occasionally, the soil samples considered are mixed samples, topsoil samples, or anthropogenically influenced soils.
after capillary saturation.Also shown are three exemplary measurement series, which are distinguished by their different colors, for measurements (c) with the punctual method, (d) evaporation method with temporary measurements, and (e) evaporation method with continuous measurements (measuring points cannot be recognized as points due to the high data density).
For the evaporation method, oven drying and the measurement of thermal conductivity at a moisture content θ = 0% were not carried out for each sample.In the case of high implausible deviations of λ within one measurement series, the deviating measurement points were removed manually.
All soil samples were installed in the cylinders with volumes of 397 cm 3 or 944 cm 3 with a specific bulk density using undisturbed sampling.Possible changes in bulk density due to shrinkage and swelling processes during the evaporation method [18,39] were neglected for validation.These changes in density are responsible for the fact that the maximum thermal conductivity in the evaporation method is sometimes not reached at saturation, but rather at an intermediate moisture content.In the case of slightly fluctuating bulk densities in measurement series measured using the punctual method, the average bulk density was assumed for validation.The individual bulk densities deviate from the mean value by a maximum of 0.04 g/cm 3 .
Due to the different methods or measurement intervals used, the number of measurement points per measurement series varies.However, for the validation, each measurement series was considered equally for authentication, independent of the number of measurement points.
As the measurement series originate from various research projects, the soil samples were taken from different locations in Germany; therefore, a wide range of soil texture classes was covered for the validation (Figure 2).Therefore, only the USDA soil texture classes sandy clay and silt were not considered for the validation.Very occasionally, the soil samples considered are mixed samples, topsoil samples, or anthropogenically influenced soils.Grain size limits were determined according to the grain size limits of the DIN EN ISO 14688-1 [40] or USDA soil classification.The deviations at the silt-sand boundary at 0.063 mm, according to DIN EN ISO 14688-1, and at 0.05 mm, according to USDA soil classification, were neglected here.The red line illustrates the boundary of >50% and <50% sand content (fSand).Different measurement series can have the same grain size distribution.If present, particles coarser than 2 mm were sieved out before the measurements.For two soil samples with a sand content ≥ 93%, the remaining finegrained content was assumed to be half clay and half silt.
A statistical evaluation of the database regarding the decisive soil parameters is shown in the Table 3. Grain size limits were determined according to the grain size limits of the DIN EN ISO 14688-1 [40] or USDA soil classification.The deviations at the silt-sand boundary at 0.063 mm, according to DIN EN ISO 14688-1, and at 0.05 mm, according to USDA soil classification, were neglected here.The red line illustrates the boundary of >50% and <50% sand content (f Sand ).Different measurement series can have the same grain size distribution.If present, particles coarser than 2 mm were sieved out before the measurements.For two soil samples with a sand content ≥ 93%, the remaining fine-grained content was assumed to be half clay and half silt.
A statistical evaluation of the database regarding the decisive soil parameters is shown in the Table 3.
Table 3. Statistical analysis of the data base based on the work of Wessolek et al. [34], with λ = thermal conductivity, θ = moisture content, ρ b = bulk density, S r = saturation and grain size distributions regarding clay (f Clay ), silt (f Silt ), and sand content (f Sand ).For this overview, the complete water content range is considered here.

Approach for Validation of Selected Empirical Literature Models
For the validation, the focus was on two key indicators that demonstrate the match between the thermal conductivities calculated using the literature models and the measured values.Therefore, the RMSE (root mean square error) and the bias value were used as follows: The RMSE corresponds to the absolute error or deviation.Squaring the deviation prevents positive and negative deviations from neutralizing each other.The bias value corresponds to an average deviation.Here, positive and negative deviations neutralize each other; thus, a negative value indicates that a literature model underestimates the thermal conductivity, on average, compared to measured values.A positive value indicates an overestimation.
For the validation, the RMSE and the bias values were determined for each measurement series based on their measurement points (n).The values of both key indicators were then averaged for all measurement series.As a result, the measurement series were considered equivalent, regardless of their number of measurement points.
For the overall analysis, all measurement series were considered regardless of any soil properties.In addition, for the development of a practicable guideline, the measurement series were divided into two categories according to their grain size distribution, specifically their sand content (f sand ).The subdivision was based on the work of Kersten [25].A distinction was made between clay and silt soils, with a sand content f Sand ≤ 50%, and sandy soils, with f Sand > 50% (Table 4).Furthermore, the validation of the two soil texture categories was also carried out regarding the moisture content.For this differentiation, only measurement points in the measurement series with saturation degrees S r ≤ 0.5 were considered for dry conditions.For moist conditions, only measurement points with S r > 0.5 were used.The calculation of the saturation is described in Equation (4).
To ensure comparability, measuring points in the very dry range were not included in the validation, in accordance with the validity limits of the model suggested by Kersten [25].Therefore, for clay and silty soils, only the measured values with a gravimetric water content w ≥ 7% were used, and for sandy soils, only the measured values with w ≥ 1% were employed.As a result, in some cases the number of measuring points described in Table 2 was reduced.
As the aim was to develop a practical guideline for the use of literature models for calculating thermal conductivity, the evaluation focused mainly on the models showing the smallest deviations from the measured values.Moist Sand (S) 25

Validated Models for Calculating Thermal Conductivity
The basic principle for the selection of the calculation models was that not more than the basic soil input parameters of grain size distribution, bulk density, and moisture content must be known, or the parameters derived from and the equations will not contain any free parameters.The selection was based on the models used by Wessolek et al. [34].
From the basic input parameters, the porosity (Φ) and saturation (S r ) were calculated as follows, assuming a particle density (ρ p ) of 2.65 g/cm 3 : S r = saturation [-] θ = volumetric water content or moisture content [-] In some cases, the calculated maximum saturation of a measurement series can be greater than 1.0, e.g., due to a deviating particle density or swelling processes within the course of the evaporation method.In these cases, the maximum measured volumetric water content of the measurement series is set as the porosity.
Various models differentiate between fine and sandy soils.Unless otherwise stated, fine-grained soils were defined with a sand content of ≤50%.
The quartz content is calculated according to the method of Hu et al. [41]: For the validation, the Python programming language (Python Software Foundation, v3.11) was used.The script can be downloaded at https://github.com/wehnsdaefflae/fusl(accessed on 16 April 2024).
In the following section, the selected models are described in detail.

Kersten (1949) Model
With the model proposed by Kersten [25], for the calculation of the thermal conductivity, only the parameters of bulk density and gravimetric water content are required.Furthermore, the soils are divided into two categories, depending on their sand content (fine textured soils with f Sand ≤ 50% and granular soils with f Sand > 50%).Thus, the knowledge of the exact grain size distribution is not required.
In the work of Kersten [25], the equations are given in imperial units.Farouki [14] has rewritten these with a conversion factor for metric units.
For fine-grained soils, the thermal conductivity for w ≥ 7% is calculated as follows: Thermal conductivity is calculated as follows for granular soils of w ≤ 1%: Johansen [35,36] proposed a new model for calculating the λ of partial saturated soils as a function of the thermal conductivity in both the dry and saturated state.The thermal conductivity between these is interpolated using the degree of saturation and a normalized thermal conductivity (Kersten's number): For coarse soils, which are defined by f Clay ≤ 5%, the Kersten's number can be calculated for S r > 0.05, as follows: In this study, Ke = 0 was assumed when S r < 0.05.For fine soils, which are defined by f Clay > 5%, Ke can be calculated for the range where S r > 0.1, as follows: Ke = log S r + 1 (10) In the case of S r < 0.1, in this study, Ke = 0 was assumed.The calculation of the thermal conductivity of saturated and dry soils is described in Equations ( 11)-(13).Johansen [35] assumed a value of 0.57 W/(m•K) for the thermal conductivity of water.The thermal conductivity of quartz is assumed to be 7.7 W/(m•K).The remaining solid components have a thermal conductivity of 2.0 W/(m•K).For quartz contents lower than 20%, a thermal conductivity of 3.0 W/(m•K) is set for the remaining components.
The model of Brakelmann [10] is based on experimental results.The thermal conductivity is calculated as follows: For the calculation of λ s , the approach described in Wessolek et al. [34] is used in this study: λ s = 0.0812f sand + 0.054f silt + 0.02f clay (15) 2.3.4.Ewen and Thomas (1987) Model Based on the Johansen (1975) model, Ewen and Thomas [23] proposed a new approach for calculation of Kersten's number.The original Kersten's number, according to Johansen [35], cannot be applied for dry soils.In the approach of Ewen and Thomas [23], this was changed using Equation ( 16), and a new parameter with the value of −8.9 was obtained: ξ = fitted parameter of Ewen and Thomas (1987) model [−8.9]The further calculation steps correspond to those of the Johansen (1975) model, as shown in Equations ( 8) and ( 11)-(13).For Equation ( 12), a value of 2.0 W/(m•K) is used for the thermal conductivity of all minerals other than quartz.
Markle et al. [28] expanded the Ewen and Thomas (1987) model and used ξ as a fitting parameter.However, this model was not considered for validation, as it contains a free parameter.The thermal conductivity is calculated with Equation ( 8).Ke is determined using a new approach: Ke = 0.1811ln S r + 0.9878 (17) In the work of Hu et al.
[24], λ dry and λ sat were measured.In this study, Equations ( 11) and (13) were used for the calculation, with a thermal conductivity of water of 0.6 W/(m•K) and of solid particles of 3.35 W/(m•K), following the method of He et al. [32].

Côté and Konrad (2005) Model
According to the approach by Côté and Konrad [21], λ and the saturated thermal conductivity are calculated based on the Johansen (1975) model, using Equation ( 8) and (11).For the thermal conductivity of water, a value of 0.6 W/(m•K) is assumed here.
Instead of calculating the thermal conductivity of the solid particles, the values for different material are taken from Table 5.A new equation for the calculation of Kersten's number was introduced (Equation ( 18)).The soil texture dependent parameter κ can be taken from Table 6.The calculation of the thermal conductivity of dry soils is also based on two new parameters, whose values can be taken from Table 7:  The model by Yang et al. [30] is also based on the Johansen (1975) model, using Equations ( 8), (11), and (12), with a value of 0.5 W/(m•K) for λ w and 2.0 W/(m•K) for λ o .The equation for the calculation of the thermal conductivity of dry soils (Equation ( 20)) is a modification of Equation (13).The coefficient was originally set to 0.135 by Peters-Lidard et al. [42] and was adopted in this model.
ρ p = particle density [2700 kg/m 3 ] The Kersten's number is calculated as follows: The model of Lu et al. [26] expresses the nonlinear relationship between thermal conductivity, moisture content, soil texture, and bulk density, introducing the parameters α and β: λ dry = −0.56Φ+ 0.51 (23) α = 0.67f clay + 0.24 (24) Markert et al. [27] proposed a model based on the approach of Lu et al. [26], using Equation ( 22) and replacing the coefficients of Equations ( 23)-( 25) with parameters that depend on soil texture: The values for p 1-8 are listed in Table 8.For the sub-models "u + p", "u", and "p", a subdivision of three soil texture groups must be carried out.The "Sand" group is defined in Equation ( 29), and the "Silt" group in Equation (30).All soils which do not belong to the texture groups "Sand" or "Silt" correspond to the texture group "Loam".
f sand < 50% and f clay < 27% Table 8.Parameters depending on texture group and soil sample packing, with u = unpacked and p = packed, taken from the work of Markert et al. [27], and the designation of the sub-models used in this study.The diversity of the models, highlighting their reliance on different soil properties as input parameters, is shown in Table 9.

Sub-Model
The water content and bulk density are required for each selected model.The quartz content, porosity, and acordingly, the saturation parameters, can be derived from the basic input parameters (Equations (3)-( 5)).
With the approach of the Hu et al. ( 2001) model chosen in this study and the assumptions made, no knowledge of grain size distribution is required.However, for the Côté and Konrad (2005) model, this parameter must be known at least sufficiently to categorize the material as shown in Table 6.The Kersten (1949) model can be used without further derived parameters to calculate the thermal conductivity based on a rough estimate of the sand content (≤50% or >50%), so the exact grain size distribution is not required.For the other models, the grain size distribution-at least the exact sand content-must be known to perform the thermal conductivity calculation.

Validation Using the Entire Set of Measurement Data
The results of the validation considering the entire set of measurement data are shown in Table 10 and as scatter plots in Figures 3 and A1.
Table 10.Results of the validation in the form of the key indicator RMSE and bias considering the entire measurement dataset, with w ≥ 1% for sandy soils and w ≥ 7% for clay and silt soils, for all selected literature models.

RMSE (W/(m•K)) Bias
Kersten Overall, RMSE values between 0.28 and 0.47 W/(m•K) were determined.The bias values were almost exclusively in the negative range, meaning that the calculated thermal conductivity was underestimated, on average, by the literature models when compared to the measured values.
In the validation, the Hu et al. (2001) model, using only bulk density and water content-derived input parameters, showed the lowest RMSE value of 0.28 W/(m•K), with an average underestimation of the thermal conductivity with a bias value of −0.10.The RMSE value for the Markert et al. (2017) model, here the texture group unspecific submodel, was only marginally higher, with a value of 0.29 W/(m•K).With a bias value of −0.05, the thermal conductivity was underestimated, to a minor extent.The texture group specific sub-models "u + p" and "u" also showed comparably low deviations from the measured values, with RMSE values of 0.31 and 0.32 W/(m•K).Similarly, a low RMSE value of 0.31 W/(m•K) was determined for the Ewen and Thomas (1987) model, again underestimating thermal conductivity, with a bias value of −0.16.Overall, RMSE values between 0.28 and 0.47 W/(m•K) were determined.The bias values were almost exclusively in the negative range, meaning that the calculated thermal conductivity was underestimated, on average, by the literature models when compared to the measured values.
In the validation, the Hu et al. (2001) model, using only bulk density and water content-derived input parameters, showed the lowest RMSE value of 0.28 W/(m•K), with an average underestimation of the thermal conductivity with a bias value of −0.10.The RMSE value for the Markert et al. (2017) model, here the texture group unspecific submodel, was only marginally higher, with a value of 0.29 W/(m•K).With a bias value of −0.05, the thermal conductivity was underestimated, to a minor extent.The texture group specific sub-models "u + p" and "u" also showed comparably low deviations from the measured values, with RMSE values of 0.31 and 0.32 W/(m•K).Similarly, a low RMSE value The Kersten (1949) model, for which the knowledge of the exact grain size distribution is not required for the calculation of the thermal conductivity, exhibits a higher error, with an RMSE value of 0.37 W/(m•K) compared to the error rates of the models previously mentioned.However, the Johansen (1975)

Validation of the Models Considering Subsets in Terms of Soil Texture and Moisture
The results of the validations considering different subsets in terms of moisture content and soil texture are graphically shown in Figure 4 to express general patterns.−0.26, the Kersten (1949) model showed a significant underestimation of the thermal conductivity compared to the measured values.

Validation of the Models Considering Subsets in Terms of Soil Texture and Moisture
The results of the validations considering different subsets in terms of moisture content and soil texture are graphically shown in Figure 4 to express general patterns.In general, but not for all models, the RMSE values appeared to be higher for sandy soils in comparison to silt and clay soils.Sandy soils generally have higher thermal conductivity than fine-grained soils.The quartz content, and possibly, the grain size of the sand content, can be decisively influencing variables on the thermal conductivity, which are only assumed or not considered in detail in the models studied here.
With the bias values, almost exclusively negative values, i.e., an underestimation of the thermal conductivity, were achieved.In general, higher negative values were shown for dry soil conditions.Particularly high bias amounts resulted for dry clay and silt soils.The partially significant negative bias values for dry soils, especially for silty and clay soils could be due to shrinking processes that may occur during the evaporation method, along with the resulting increase in bulk density, which further positively influence the measured thermal conductivity [18].However, this correlation does not apply to all soil texture categories or models, including the Ewen and Thomas (1987)  In general, but not for all models, the RMSE values appeared to be higher for sandy soils in comparison to silt and clay soils.Sandy soils generally have higher thermal conductivity than fine-grained soils.The quartz content, and possibly, the grain size of the sand content, can be decisively influencing variables on the thermal conductivity, which are only assumed or not considered in detail in the models studied here.
With the bias values, almost exclusively negative values, i.e., an underestimation of the thermal conductivity, were achieved.In general, higher negative values were shown for dry soil conditions.Particularly high bias amounts resulted for dry clay and silt soils.The partially significant negative bias values for dry soils, especially for silty and clay soils could be due to shrinking processes that may occur during the evaporation method, along with the resulting increase in bulk density, which further positively influence the measured thermal conductivity [18].However, this correlation does not apply to all soil texture categories or models, including the Ewen and Thomas (1987)  The RMSE and bias values for the individual models and subsets are listed in detail in Tables 11 and 12.
For the Hu et al. ( 2001) model, a very low RMSE value of 0.20 W/(m•K) was obtained for clayey and silty soils, largely independent of the saturation, with only a slight underestimation of λ, on average.For this soil texture category, moisture had almost no influence on the quality of the model.However, the deviations of the calculated thermal conductivities compared to the measured values were significantly higher for sandy soils, with an RMSE value of 0.38 W/(m•K).The deviations were the highest for moist sandy soils, with an RMSE value of 0.41 W/(m•K) and a bias value of −0.26.
The texture group unspecific Markert et al. (2017) model showed low RMSE values of 0.22 and 0.30 W/(m•K), depending on the soil texture and moisture.The thermal conductivity for clay and silt soils was underestimated by a bias value of −0.12, on average.For sandy soils, the thermal conductivity was slightly overestimated compared to the measured values, with a bias value of 0.05.The "u + p" sub-model also showed similarly small deviations for sandy soils.organic soils, high errors using empirical models are to be expected.For example, thermal conductivity measurements of decomposed basalt showed significantly lower values in comparison to those for loamy valley deposits, despite similar grain size distribution and bulk density [18].This is mainly due to the low particle thermal conductivity of basalt [21].
In addition to the validation of models for the calculation of the thermal conductivity of soils based on the parameters of bulk density, moisture content, and grain size distribution, the aim of the study was to develop a practical guideline for the use of empirical literature models, additionally based on the validation results of the considered subsets.For the development of the guideline, the determined RMSE and bias values of the validation of this study were used, as well as the results for the selected moisture and soil texture categories.The guideline is shown in Figure 5.
In  2017) model is proposed here, if information on grain size distribution is available.As can be seen in Figure 5a and Table 11, soil moisture is not a decisive factor for the selection of the literature model.The Ewen and Thomas (1987) model showed the lowest RMSE value for dry soils in general, which was, however, neglected here for reasons of simplicity.The results can therefore be broken down into a simple decision tree (Figure 5b).When using the Markert et al. (2017) model to calculate the thermal conductivity of sandy soils, as proposed here, the slight overestimation of the thermal conductivity by the model should be considered.
If the exact grain size distribution is not known, the Hu et al. ( 2017) model, with the approach and assumptions applied in this study, can be used to generate particularly high accuracy for clayey and silty soils.In the applied approach, only porosity and saturation are required as input parameters for the calculation of the thermal conductivity, while other parameters are assumed and not calculated, as in comparable models.The Ewen and Thomas (1987) model showed differences between the validation results for dry and wet soil conditions.Low RMSE values of up to 0.19 and 0.30 W/(m•K) were obtained for dry soils, whereby the thermal conductivity for clay and silt soils is underestimated by a bias value of 0.10, on average.No underestimation or overestimation was observed for dry sandy soils.For moist soils, higher RMSE values of 0.31 and 0.35 W/(m•K) and a considerable underestimation, with bias values of −0.24 and −0.20, were achieved.
The Kersten (1949) model showed the highest deviation, with an RMSE value of 0.42 W/(m•K) for dry clay and silt soils and a clear average underestimation of the thermal conductivity compared to the measured values, with a bias value of −0.53.The smallest deviations were found for moist soils, with an RMSE value of 0.33 W/(m•K).

Discussion
For examples based on the Johansen (1975) model, the quartz content is required for the calculation of the thermal conductivity.In this study, the quartz content was assumed by deriving it from the sand content; thus, this parameter exhibits significant uncertainty.Changing the assumed quartz content, e.g., due to a different calculation equation, can lead to significantly different results [32].Therefore, actual measurements of the thermal conductivity are always preferable.Especially in the case of specific mineralogical compositions or organic soils, high errors using empirical models are to be expected.For example, thermal conductivity measurements of decomposed basalt showed significantly lower values in comparison to those for loamy valley deposits, despite similar grain size distribution and bulk density [18].This is mainly due to the low particle thermal conductivity of basalt [21].
In addition to the validation of models for the calculation of the thermal conductivity of soils based on the parameters of bulk density, moisture content, and grain size distribution, the aim of the study was to develop a practical guideline for the use of empirical literature models, additionally based on the validation results of the considered subsets.For the development of the guideline, the determined RMSE and bias values of the validation of this study were used, as well as the results for the selected moisture and soil texture categories.The guideline is shown in Figure 5.  2017) model is proposed here, if information on grain size distribution is available.As can be seen in Figure 5a and Table 11, soil moisture is not a decisive factor for the selection of the literature model.The Ewen and Thomas (1987) model showed the lowest RMSE value for dry soils in general, which was, however, neglected here for reasons of simplicity.The results can therefore be broken down into a simple decision tree (Figure 5b).When using the Markert et al. (2017) model to calculate the thermal conductivity of sandy soils, as proposed here, the slight overestimation of the thermal conductivity by the model should be considered.
If the exact grain size distribution is not known, the Hu et al. (2017) model, with the approach and assumptions applied in this study, can be used to generate particularly high accuracy for clayey and silty soils.In the applied approach, only porosity and saturation are required as input parameters for the calculation of the thermal conductivity, while other parameters are assumed and not calculated, as in comparable models., is considered to be less accurate here due to the higher RMSE values for the sandy soils compared to those of the approach by Markert et al. [27].
Furthermore, if the exact grain size distribution is not known, the Kersten (1949) model can be used to calculate the thermal conductivity, based on a rough estimate of the sand content (≤50% or >50%) and using a knowledge of the bulk density and water content, with no further assumptions.This simplicity enables the Kersten (1949) model to be used for preliminary thermal exploration, e.g., for derivation of thermal conductivity from geoelectric field measurements [18,43].The validation of this study showed an RMSE value of 0.37 W/(m•K).This significant underestimation would provide a large degree of certainty regarding the estimation of thermal conductivity for the dimensioning and design of thermal infrastructure projects.
For the sake of completeness, however, the results of other validation studies should also be briefly mentioned.In the study by He et al. [32], the models of Johansen [35,36], Côté and Konrad [21], Côté and Konrad [44], Lu et al. [45], Tarnawski et al. [46], and Yan et al. [47] were determined to be the best-fitting models without free parameters.In addition to the model of Markle et al. [28], which, however, requires a fitting parameter, Wessolek et al. [34] also found the Hu et al. (2001) model to be the model with the lowest RMSE value.However, the Markert et al. (2017) model was not included in both studies.

Conclusions
By validating selected empirical models for calculating the thermal conductivity of unconsolidated soils using 59 series of measurements, it was possible to derive a practicable guideline in the form of a decision tree for the use of thermal conductivity models considering different soil texture categories.This is based on knowledge of the physical soil parameters of moisture content, bulk density, and if applicable, grain size distribution.
The measurement series each consist of a different number of measurement points.However, each measurement series received equal consideration for validation.The measurements were carried out on soil samples with a known grain size distribution, bulk density, and moisture content using single-needle sensors.As the measurement series originate from various research projects in Germany, a wide range of soil texture classes was included.
The correspondence between the calculated thermal conductivities and the measured values was evaluated using the RMSE and the bias value.As a result, it was determined that the thermal conductivity is mainly underestimated by the literature models, on average, with basically higher deviations for sandy soils.2001) model, with the approach and assumptions used in this study, can also be employed to calculate thermal conductivity when grain size distribution information is not available.
If the exact grain size distribution is not known, the Kersten (1949) model can also be used by roughly estimating the sand content.Although this model shows higher deviations, the significant underestimation of the thermal conductivity on average provides a large degree of certainty for the planning of thermal infrastructure projects.
As the soil samples were obtained exclusively from Germany, the applicability of the results or the guideline on a global scale is limited.Future research could therefore extend the validation carried out to international soil samples.In addition, the guideline could be extended by the integration of further models or by additional subset analyses, e.g., with regard to organic matter content or bulk density.
Due to the small number of input parameters and the high accuracy of the Hu et al. (2001) model, future attempts could also be made to perform thermal exploration using geoelectric field measurements by deriving thermal conductivities using this model.

Figure 1 .
Figure 1.Schematic illustration of the (a) evaporation method and (b) the punctual method of thermal conductivity measurement used for validation, taken from the work of Rammler et al. [18], with ω1/2 = medium water contents, ωdry = water content after oven drying, and ωsat = water content

Figure 2 .
Figure 2. Grain size distributions and soil texture classes according to USDA soil classifications.

Figure 2 .
Figure 2. Grain size distributions and soil texture classes according to USDA soil classifications.Grain size limits were determined according to the grain size limits of the DIN EN ISO 14688-1[40] or USDA soil classification.The deviations at the silt-sand boundary at 0.063 mm, according to DIN EN ISO 14688-1, and at 0.05 mm, according to USDA soil classification, were neglected here.The red line illustrates the boundary of >50% and <50% sand content (f Sand ).Different measurement series can have the same grain size distribution.If present, particles coarser than 2 mm were sieved out before the measurements.For two soil samples with a sand content ≥ 93%, the remaining fine-grained content was assumed to be half clay and half silt.

2. 3 . 5 .
Hu et al. (2001) Model With the approach of Hu et al. [24], the thermal conductivity of unconsolidated porous media is calculated based on the Johansen (1975) model and the Leverett-Lewis equation for the relationship of capillary pressure saturation.

Figure 3 .
Figure 3. Results of the validation as exemplary scatter plots for the models of (a) Kersten [25], (b) Brakelmann [10], (c) Hu et al. [24] and (d) Markert et al. [27] (here: texture group independent submodel) considering the entire measurement data, with w ≥ 1% for sandy soils and w ≥ 7% for clay and silt soils.Measurement points of measurement series with a high data density (98 to 1882 measurement points per measurement series) are shown as blue dots.Measurement points from measurement series with low data density (3 to 10 measurement points per measurement series) are shown as black crosses.

Figure 3 .
Figure 3. Results of the validation as exemplary scatter plots for the models of (a) Kersten [25], (b) Brakelmann [10], (c) Hu et al. [24] and (d) Markert et al. [27] (here: texture group independent sub-model) considering the entire measurement data, with w ≥ 1% for sandy soils and w ≥ 7% for clay and silt soils.Measurement points of measurement series with a high data density (98 to 1882 measurement points per measurement series) are shown as blue dots.Measurement points from measurement series with low data density (3 to 10 measurement points per measurement series) are shown as black crosses.
model, the Côté and Konrad (2005) model, and the Yang et al. (2005) model showed significantly higher deviations.With a bias value of −0.26, the Kersten (1949) model showed a significant underestimation of the thermal conductivity compared to the measured values.

Figure 4 .
Figure 4. Bar charts of (a) RMSE and (b) bias values for the considered subsets in terms of soil texture and moisture for all selected models [10,21,23-27,30,35,36].As a representative of the Markert et al. (2017) models, only the "unspecific" texture group sub-model was considered.

Figure 4 .
Figure 4. Bar charts of (a) RMSE and (b) bias values for the considered subsets in terms of soil texture and moisture for all selected models [10,21,23-27,30,35,36].As a representative of the Markert et al. (2017) models, only the "unspecific" texture group sub-model was considered.
this study, the Hu et al. (2017) model and the Markert et al. (2017) model (here: texture group unspecific sub-model) were identified as the generally best-fitting models.The guideline is based solely on the results of this study.The RMSE values for the Hu et al. (2001) model and for the Markert et al. (2017) model are almost identical for the soil texture-and moisture-independent approach.However, since the Hu et al. (2001) model shows significantly higher RMSE values for sandy soils, the Markert et al. (

Figure 5 .
Figure 5. (a) Basis of evaluation for the guideline derived from the results of the validation for the use of empirical literature models to calculate the thermal conductivity of unconsolidated soils (b) and the guideline in form of a decision tree.The basis for this is the knowledge (or assumption) of the bulk density and moisture content of the soil.The RMSE and bias values (RMSE in W/(m•K)/bias) are indicated in each case to assess the quality of the model.For sandy soils, the texture group independent sub-model, or the "u + p" sub-model of Markert et al. [27], can be used.In the figure, RMSE and bias values are shown for the texture group unspecific sub-model.Negative bias values correspond to an average underestimation of the thermal conductivity by the model.The validation was carried out for the range, with a gravimetric water content ≥ 1% for sandy soils and ≥7% for clay and silty soils.The calculation of the thermal conductivity with unknown grain size distribution, using the model of Hu et al. [24], is considered to be less accurate here due to the higher RMSE values for the sandy soils compared to those of the approach by Markert et al.[27].
The Hu et al. (2001) model, which is based on the Johansen (1975) model, and the Markert et al. (2017) model were identified as the best-fitting models.The Hu et al. (2001) model showed particularly low RMSE values for clay and silty soils, while for sandy soils, the Markert et al. (2017) model showed the lowest deviations.The Hu et al. (

Figure A1 .FigureFigure A1 .
Figure A1.Results of the validation as exemplary scatter plots for the models of (a) Johansen[35,36],

Figure A1 .
Figure A1.Results of the validation as exemplary scatter plots for the models of (a) Johansen[35,36], (b) Ewen and Thomas[23], (c) Côté and Konrad[21], (d) Yang, et al.[30], (e) Lu, et al.[26] and (f-h) Markert, et al.[27] (here: sub-models "u + p", "u" and "p") considering the entire measurement data, with w ≥ 1% for sandy soils and w ≥ 7% for clay and silt soils.Measurement points of measurement series with a high data density (98 to 1882 measurement points per measurement series) are shown as blue dots.Measurement points from measurement series with low data density (3 to 10 measurement points per measurement series) are shown as black crosses.

Table 1 .
Thermal properties analyzer and the specifications used for thermal conductivity measurements.

Table 2 .
Measurement methods used for the thermal conductivity measurements.For each thermal conductivity measurement, the moisture content was determined.Here, for the number of measurement points per measurement series, the complete water content range is considered.
[20]poration method based on the approach of Markert et al.[20], with continuous measurements, obtained at 15 min to 60 min intervals, or temporary measurements, e.g., performed at daily intervals; measurements performed in cylinders.

Table 1 .
Thermal properties analyzer and the specifications used for thermal conductivity measurements.

Table 4 .
Soil properties of the subsets considered in the validation with regard to soil texture and moisture, as well as their designation used for the specific soil texture categories.

Table 5 .
[21]age values for the thermal conductivity of solid particles (λ s ) of selected materials, taken from the work of Côté and Konrad[21].

Table 9 .
Basic input parameters and further derived parameters required (unless not measured) to calculate the thermal conductivity using each empirical model.