the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Probabilistic estimation of depth-resolved profiles of soil thermal diffusivity from temperature time series
Carlotta Brunetti
John Lamb
Stijn Wielandt
Sebastian Uhlemann
Ian Shirley
Patrick McClure
Baptiste Dafflon
Improving the quantification of soil thermal and physical properties is key to achieving a better understanding and prediction of soil hydro-biogeochemical processes and their responses to changes in atmospheric forcing. Obtaining such information at numerous locations and/or over time with conventional soil sampling is challenging. The increasing availability of low-cost, vertically resolved temperature sensor arrays offers promise for improving the estimation of soil thermal properties from temperature time series, and the possible indirect estimation of physical properties. Still, the reliability and limitations of such an approach need to be assessed. In the present study, we develop a parameter estimation approach based on a combination of thermal modeling, sliding time windows, Bayesian inference, and Markov chain Monte Carlo simulation to estimate thermal diffusivity and its uncertainty over time, at numerous locations and at an unprecedented vertical spatial resolution (i.e., down to 5 to 10 cm vertical resolution) from soil temperature time series. We provide the necessary framework to assess under which environmental conditions (soil temperature gradient, fluctuations, and trend), temperature sensor characteristics (bias and level of noise), and deployment geometries (sensor number and position) soil thermal diffusivity can be reliably inferred. We validate the method with synthetic experiments and field studies. The synthetic experiments show that in the presence of median diurnal fluctuations ≥ 1.5 ∘C at 5 cm below the ground surface, temperature gradients > 2 ∘C m−1, and a sliding time window of at least 4 d the proposed method provides reliable depth-resolved thermal diffusivity estimates with percentage errors ≤ 10 % and posterior relative standard deviations ≤ 5 % up to 1 m depth. Reliable thermal diffusivity under such environmental conditions also requires temperature sensors to be spaced precisely (with accuracy to a few millimeters), with a level of noise ≤ 0.02 ∘C, and with a bias defined by a standard deviation ≤ 0.01 ∘C. Finally, the application of the developed approach to field data indicates significant repeatability in results and similarity with independent measurements, as well as promise in using a sliding time window to estimate temporal changes in soil thermal diffusivity, as needed to potentially capture changes in bulk density or water content.
- Article
(4401 KB) - Full-text XML
- BibTeX
- EndNote
Knowledge of soil thermal properties (i.e., soil thermal conductivity, thermal diffusivity, and specific heat capacity) is fundamental to solving problems in many fields, such as in engineering, agriculture, meteorology, and geology (Farouki, 1981). Thermal properties, which are controlled by the spatial arrangement and fraction of soil components, such as minerals, organic matter, water, ice, and air, modulate heat fluxes in soil and at the soil–surface boundary. As a consequence, an improved quantification of soil thermal properties is a cornerstone for advancing the indirect estimation of fraction of soil components (Al Nakshabandi and Kohnke, 1965; Ochsner et al., 2001; Abu-Hamdeh, 2003; Arkhangel'skaya, 2009; Tong et al., 2016; Arkhangel'skaya and Lukyashchenko, 2018; Xie et al., 2018; Jafarov et al., 2020) needed to parametrize hydro-biogeochemical models, as well as for understanding and predicting variability in the subsurface thermal regime and fluxes of water, carbon, and nutrients (Koven et al., 2013; Rasmussen et al., 2018; Nicolsky and Romanovsky, 2018; Oliva and Fritz, 2018).
Direct measurement of soil thermal and physical properties (Ochsner and Baker, 2008; Mengistu et al., 2017) obtained by sampling and analyzing soil from multiple depths and locations is highly valuable, but it is also time-consuming, invasive, difficult to repeat over time, and costly. In addition, the sampling and analysis are not error free, primarily because of the difficulty in minimizing error in the bulk density measurement. Similarly, this approach provides limited vertical and lateral resolution, which in some cases is insufficient to capture the highly heterogeneous distribution of soil physical and thermal properties. An assessment of alternative methods is needed that can potentially complement direct measurements and provide estimates of physical and thermal properties over space and time with unprecedented resolution. Soil temperature time series have been proven to be a valuable source of information for monitoring spatiotemporal changes in subsurface properties, and they have been used to infer (for instance) heat and water fluxes (e.g., Steele-Dunne et al., 2010; An et al., 2016; Irvine et al., 2017; Tabbagh et al., 2017), organic matter content (Tran et al., 2017), and thermal properties (e.g., Nicolsky et al., 2009; Rajeev and Kodikara, 2016).
The estimation of thermal properties from soil temperature time series can be achieved through analytical or inverse methods. With analytical methods, thermal diffusivity is computed from the analytical solution of the heat equation, under the assumption that soil surface temperature is a sinusoidal function. In the literature, thermal diffusivity has been estimated from the analytical solution of the heat conduction equation using, for instance, the harmonic (Fourier) method (Carson, 1963; Beardsmore et al., 2020) and the popular phase-shift and amplitude-ratio methods (Krzeminska et al., 2012; Hinkel, 1997; Jong van Lier and Durigon, 2013; Andújar Márquez et al., 2016) or from the analytical solution of the conduction–convection equation (Gao et al., 2017, 2008). The main advantage of analytical methods is that they provide daily or monthly variation in bulk thermal diffusivity without the need for computational resources. These estimates are important input in physics-based models that simulate land–atmosphere interactions (Gao et al., 2017). One limitation of these methods is that they cannot be applied when the temperature does not vary sinusoidally, such as on rainy or cloudy days. Moreover, since the amplitude of the soil temperature wave decreases exponentially with depth, the phase shift and amplitude methods might be unable to infer thermal diffusivity deeper than 20–30 cm (Farouki, 1981; Gao et al., 2017) unless very strong diurnal fluctuations are recorded at the surface of conductive soils. Additionally, most of the studies that rely on analytical methods are typically characterized by a low vertical spatial resolution, which often involves the estimation of thermal properties in a single layer in the top 20 cm (e.g., Gao et al., 2017), in two layers in the top 30 cm (e.g., Jong van Lier and Durigon, 2013), or in three layers in the top 40 cm (Krzeminska et al., 2012; Gao et al., 2008). A few studies have investigated thermal properties at greater depths, such as up to 1 m with a resolution of 20 cm (Beardsmore et al., 2020) or with a single layer up to 5 m (Andújar Márquez et al., 2016).
Some of the challenges inherent to the analytical methods can be overcome by the use of numerical approximations of the heat equation (which can take in input temperature time series of any form) and by embedding them in a deterministic or probabilistic inverse method. A deterministic inversion relies on optimization techniques to estimate the value of the unknown parameter(s) of interest (e.g., soil thermal properties) that minimize the difference between the simulated and measured data (e.g., soil temperature time series). Deterministic inversion of thermal properties has been achieved with the finite difference (Krzeminska et al., 2012) or the finite element scheme to approximate the heat conduction (Kim et al., 2019), the conduction–convection equation (Tabbagh et al., 2017), or the heat conduction equation with phase changes (Nicolsky et al., 2009). Comparison of analytical and numerical methods can be found in Horton et al. (1983) and Rajeev and Kodikara (2016). Deterministic inversions are computationally inexpensive, but the convergence to the global minimum is not guaranteed. Similar to the studies based on analytical methods, most of the investigations done with deterministic inverse approaches infer thermal properties at a low vertical resolution, such as in a single layer (e.g., Tabbagh et al., 2017) or two to three layers in the top 40 cm (e.g., Kim et al., 2019; Krzeminska et al., 2012), or, for instance, five layers up to 50 m (Nicolsky et al., 2009). Moreover, deterministic inversion, as with analytical methods, only provides a single solution of the system, with no means to explicitly model and take into account different sources of uncertainty.
Probabilistic inverse methods provide a set (distribution) of equally probable values of the unknown parameters by aiming at “capturing both the average response of the system and the variability due to uncertainties of any kind” (Renard et al., 2013). In the last decades, Bayesian inference has gained popularity among the probabilistic inverse approaches. Unlike classical statistical methods to estimate uncertainty (e.g., Beardsmore et al., 2020), Bayesian inference allows incorporation of prior knowledge about the parameter to be estimated as well as taking explicitly into account different sources of uncertainty, such as measurement errors. Moreover, a great variety of sampling algorithms are available to avoid nonconvergence issues. Though promising, the applications of the Bayesian framework to infer soil thermal properties are still limited. Huang et al. (2017) have derived the heat conductivity in the layer between the atmosphere and the soil surface, Choi et al. (2018) have inferred soil thermal conductivity from a thermal response test, and Tran et al. (2017) have inferred organic matter content from soil temperature, liquid water, and apparent resistivity data. Bayesian inference has been more widely applied in engineering (Kaipio and Fox, 2011) to estimate thermal properties of fins (Gnanasekaran and Balaji, 2013; Somasundharam and Reddy, 2017) and walls (De Simon et al., 2018; Rodler et al., 2019). While the above studies have shown promise in estimating thermal properties from time series of temperature, several challenges still remain, including the need for an approach that can assess the conditions under which such estimates are reliable. In addition, the recent increasing availability of dense, vertically resolved arrays of temperature sensors offers a new opportunity to increase the vertical and lateral resolution of the estimates and potentially quantify their temporal variability, although such a strategy has not been assessed yet.
The aim of our study is to evaluate the potential of estimating soil thermal diffusivity and its uncertainty at numerous locations and at an unprecedented vertical spatial resolution (i.e., 5–10 cm in the top 1 m), using vertically resolved time series of soil temperature. A second objective is to explore the possibility of estimating changes in thermal diffusivity over time by sequencing the soil temperature time series with a sliding time window. To this end, we developed a parameter estimation approach based on Bayesian inference that also allows us to assess the impact of different sources of uncertainty linked to various environmental conditions (e.g., soil temperature gradient, fluctuations, and trend), temperature sensor characteristics, and deployment geometries. In this study, the Bayesian inference approach is combined with a heat conduction thermal model to estimate thermal diffusivity under the assumption that conduction dominates advection and that water–ice phase change is absent. Though these assumptions limit the applicability of the approach to specific time windows where water fluxes are limited and temperatures are well above or below the freezing point, they avoid the need of constraining water fluxes or water content in or around the modeled domain, which would require measurements that are rarely available at numerous locations. While several studies have successfully accounted for advection and/or freeze–thaw processes when estimating soil thermal parameters, they relied on additional hydrological measurements (e.g., Jafarov et al., 2012; Tran et al., 2017) or involved fully saturated environments (e.g., Nicolsky et al., 2009). In our study, we aim at developing a method that can be applied at locations where only temperature data are available. While our approach reduces the number of parameters to be estimated or constrained (and hence the likelihood of non-uniqueness of solution), we emphasize that its applicability is limited to time windows where phase change is absent and heat advection is limited (i.e., period with little water infiltration).
Using the abovementioned approach, this study aims at answering the following research questions: (1) can soil thermal diffusivity be reliably estimated (i.e., with percentage errors ≤ 10 % and posterior relative standard deviations ≤ 5 %, which are representative of what can be presumably measured with a thermal analyzer on soil samples) at multiple depths and locations from solely vertically resolved soil temperature time series (i.e., without additional data), and if so, under which environmental conditions and sensor characteristics? (2) Can we capture temporal changes in thermal diffusivity by estimating diffusivity sequentially for a sliding time window of soil temperature time series? (3) Can thermal diffusivity, estimated by applying our method to several locations across a field site, be used to retrieve soil composition and other soil thermal–physical properties? To address these questions, we perform synthetic experiments in which we infer soil thermal diffusivity and assess its uncertainty under different environmental conditions (i.e., soil temperature gradients and fluctuations), length of sliding time window, level of measurement errors, and temperature sensor geometries (Sect. 3). In Sect. 4, we further evaluate the reliability and sensitivity of the proposed method with an in situ study that compares estimated thermal diffusivities for a silty–clayey soil in a warm summer Mediterranean climate (Berkeley, California, USA) with independent measurements obtained with a thermal-property analyzer. Finally, in Sect. 5, we apply the method at a field site in a discontinuous permafrost environment in a tundra climate (Nome, Alaska, USA) during a period of low vertical gradient in soil temperature and compare thermal diffusivity estimates at numerous locations across the site with soil sample measurements. The synthetic and in situ experiments enable evaluation and validation of the developed method and its applicability to numerous environmental conditions, including challenging situations such as the presence of low vertical gradients in soil temperature. Finally, it can be noted that the evaluation of the method in an Arctic environment is particularly important to potentially improve the parameterization of soil thermal parameters and organic matter content in ecosystem models simulating the feedback from Arctic ecosystem to climate warming.
2.1 Heat equation and thermal diffusivity
Heat is exchanged at the soil surface and within the soil through different processes, such as radiation, convection, conduction, and latent heat. However, conduction is the process that dominates the transport of heat in soil (Hillel, 1982). Assuming a purely heat conduction process in a heterogeneous medium of length L (m) under unsteady-state conditions over a time window of duration T (s), the heat equation in one dimension is (Hillel, 1982)
where temperature, u (∘C), varies with time, t (s), and depth, z (m), and thermal diffusivity, α (m s−2), changes with z. Thermal diffusivity characterizes unsteady heat conduction and defines how quickly a material transfers heat from a hot to a cold area due to a change in temperature (Farouki, 1981). High thermal diffusivity values indicate that a material is capable of a rapid transfer of heat. Thermal diffusivity is defined as the ratio of the thermal conductivity, κ (W m−1 K−1), to the density, ρ (kg m−3), and specific heat capacity, cp (J kg−1 K−1) at a constant pressure:
Thermal diffusivity of soil is influenced by (1) soil composition, (2) soil water content (through its impact on wet bulk density, specific heat capacity, and thermal conductivity), (3) soil bulk density (which depends mainly on the soil composition and its degree of compaction), and, to a lesser extent, (4) temperature changes (Al Nakshabandi and Kohnke, 1965; Ochsner et al., 2001; Abu-Hamdeh, 2003; Arkhangel'skaya, 2009; Tong et al., 2016; Arkhangel'skaya and Lukyashchenko, 2018; Xie et al., 2018; Zhu et al., 2019).
2.2 Forward modeling: finite differences and sliding time windows
The differential equation in Eq. (1) can be solved numerically by discretizing the space-time domain via finite difference or finite element methods such that an approximation of the exact solution at the grid nodes can be computed. We apply herein the most common and easy-to-implement explicit finite difference method based on the forward time and centered space scheme (Petter Langtangen and Linge, 2017; Praprotnik et al., 2004):
where i and j identify the grid node in space and time, respectively; Δt (s) is the length of the time step; Δz (m) is the distance in space between grid nodes; and represents a discrete approximation to u(z,t). Note that the thermal diffusivity is computed on a staggered grid. Both Δt and Δz are chosen so as to avoid numerical instabilities and to ensure the convergence of the explicit finite difference scheme:
which is known as the Courant–Friedrichs–Lewy condition (Courant et al., 1928). According to the literature (e.g., Farouki, 1981; Andújar Márquez et al., 2016), we consider soil thermal diffusivity values as large as m2 s−1 and Δz=0.05 m; therefore, we set Δt=300 s. While a uniform spacing of 0.05 m is used in the finite difference scheme, the actual spatial resolution at which soil thermal diffusivity is resolved by our method corresponds to the spacing between the temperature sensors. Similarly to Tabbagh et al. (2017) and Rodler et al. (2019), we sequence the soil temperature time series with a sliding time window of length T that moves every 24 h over the entire time period for which data are available. Thermal diffusivity is assumed to be constant in time within each time window. The optimal time window length T is investigated in detail later in this study, as it needs to be sufficiently small to limit the influence of hydrological processes (e.g., advection not represented in the heat conduction model), but long enough to contain sufficient information to reliably infer thermal diffusivity at multiple depths. Based on the measurement geometry assumed in this study, the soil temperature values at multiple depths at the beginning of each time window are used as initial conditions, whereas the soil temperature time series recorded over the time window by the shallowest and deepest sensors are employed as top and bottom boundary conditions (i.e., time-varying Dirichlet boundary conditions), respectively.
2.3 Inverse modeling: Bayesian inference with MCMC
Bayesian inference provides a probabilistic framework that enables us to derive from n data, , …, , the d-dimensional vector of the parameters of interest, θ, that are not directly measured or known while fully quantifying the associated uncertainty. This process is performed by the Bayes theorem:
which defines how the prior state of knowledge about the quantities of interest, p(θ), is updated by the information contained in the data through the likelihood function, . The denominator in Eq. (5) is a normalization factor called “evidence” that can be neglected if a single conceptual model is considered, as in this study. The Bayes theorem provides as output the posterior probability density function (pdf), , of the parameters of interest θ. Bayesian inference can be performed with likelihood functions of any form. However, we assume here, as it is often the case, uncorrelated and normally distributed measurement errors with constant standard deviation, , that define a Gaussian likelihood function as
The term ℱh(θ) is the forward model (e.g., Eq. 3, Sect. 2.2) used to simulate the observed data, . Larger likelihood values indicate that ℱh(θ) better predicts the data at hand. The posterior pdf in Eq. (5) is multidimensional and analytically intractable. Sampling schemes such as the popular Markov chain Monte Carlo (MCMC) (Gilks et al., 1995; Robert and Casella, 2013) algorithm are, therefore, applied to find an approximation of the posterior pdf. In particular, in this work, we make use of the DiffeRential Evolution Adaptive Metropolis, DREAM(ZS), algorithm (Laloy and Vrugt, 2012; Vrugt, 2016), which is a multi-chain MCMC sampling scheme based on the Metropolis acceptance ratio (Metropolis et al., 1953) but with an improved sampling efficiency. This algorithm allows us to avoid nonconvergence issues that might arise when inferring for more than one unknown parameter (i.e., thermal diffusivity at multiple depths). The convergence of the multiple Markov chains to the posterior pdf is assessed quantitatively with the Gelman–Rubin statistic (Gelman and Rubin, 1992). In this paper, we aim at inferring about 12 unknown parameters; therefore, following the guidelines provided by Vrugt (2016), we set the number of Markov chains to three and a total of 5×104 iterations were found to be enough to reach convergence and get acceptable acceptance rates (i.e., between 15 % and 40 %, according to Gelman et al., 1996). The unknown parameters to be inferred from the soil temperature time series are thermal diffusivity values at multiple depths that are drawn from the uniform prior distribution U[0.01, 3] mm2 s−1. The upper and lower limits of the thermal diffusivity prior range have been chosen based on typical values found in the literature (e.g., Farouki, 1981; Andújar Márquez et al., 2016). In the case of the in situ experiment (Sect. 4) and field case study (Sect. 5), we also infer the standard deviation of the measurement errors, , by drawing it from the uniform distribution U[0, 1] ∘C.
The uncertainty over the MCMC thermal diffusivity estimates is quantified with the percentage error and the posterior relative standard deviation. The percentage error (PE) is the distance of the MCMC posterior mean of thermal diffusivity, , from the synthetic true value, αtrue (in the case of synthetic experiments):
whereas the posterior relative standard deviation (RSD) measures the dispersion of the MCMC posterior distribution around the MCMC posterior mean of thermal diffusivity:
where σα is the standard deviation of the posterior MCMC distribution of thermal diffusivity.
The uncertainty regarding the inferred thermal diffusivity at multiple depths is influenced by a combination of many factors, including the temperature gradient at each depth, the length of the sliding time window used in the MCMC inversion, the amount of information contained in temperature trend and fluctuations, the level of noise and bias in the temperature time series, and the geometry of temperature sensors.
2.4 Data and field site description
2.4.1 Vertically resolved time series of temperature
Synthetic and field experiments performed in this study use an acquisition geometry in which temperature measurements are collected autonomously at numerous depths, with a vertical resolution of 0.05 to 0.1 m. This geometry corresponds to measurements obtained with vertically resolved temperature devices such as, for example, a distributed temperature profiling (DTP) system (Dafflon et al., 2022b). In this study, each DTP system has at least 13 digital temperature sensors (TMP117AIDRVR, http://www.ti.com/lit/ds/symlink/tmp117.pdf, last access: 23 June 2022) that are mounted on an interconnected printed circuit board and inserted in a ∼10 mm outer diameter plastic tube filled with epoxy. The DTP system configuration used in this study involves a spatial resolution of 0.05 m between the top five sensors and of 0.10 m between the bottom eight sensors, resulting in the acquisition of temperature time series at 0.05, 0.10, 0.15, 0.20, 0.25, 0.35, 0.45, 0.55, 0.65, 0.75, 0.85, 0.95, and 1.05 m below the ground surface. The sensors record the soil temperature every 15 min, with a resolution of 0.0078 ∘C and an accuracy (corresponding to 3 standard deviations) of ±0.1 or ±0.015 ∘C when relying on the manufacturer calibration or an additional in-house calibration, respectively (Dafflon et al., 2022b).
2.4.2 Synthetic experiment
The synthetic experiments (Sect. 3) are implemented to investigate the impact of different environmental conditions and sensor characteristics on the uncertainty of thermal diffusivity estimates. We generate synthetic soil temperature fields that mimic various types of plausible temperature gradients and fluctuations through forward modeling (Sect. 2.2). This is achieved by setting the initial, top, and bottom boundary conditions equal to in situ temperature time series, here specifically from two monitoring sites on the Seward Peninsula, Alaska. One site is used to evaluate summer conditions (Romanovsky et al., 2020) and the other (Brunetti et al., 2022) autumn conditions. Both sites are locations characterized by the absence of permafrost, the presence of tall shrubs, and decreasing organic matter content and porosity with depth. We then assume a soil column composed of three layers (i.e., top layer at 0.05–0.1 m, middle layer at 0.1–0.42 m, and bottom layer at 0.42–1.05 m) in each of which thermal diffusivity values are considered constant over time and assigned based on soil sample analysis performed at the same sites under similar environmental conditions, i.e., 0.16, 0.27, and 0.43 mm2 s−1 for the case of summer temperatures and 0.25, 0.75, and 0.6 mm2 s−1 for autumn. The summer case is characterized by median diurnal fluctuations (maximum minus minimum temperature over 1 d measured by the top sensor at 0.05 m below the ground surface) of 1.6 ∘C over the time period considered and temperature gradients all larger than 2 ∘C m−1 at each depth. The case for autumn has lower median diurnal fluctuations (i.e., 0.13 ∘C) and lower temperature gradients (mostly within −1 and 2 ∘C m−1). We perturb the synthetic temperature fields using a Gaussian noise with a standard deviation of ∘C.
2.4.3 Assessment of repeatability and temporal changes of in situ thermal diffusivity estimates
One of the two field experiments presented in this study is aimed at evaluating the repeatability of the estimated diffusivity and the method’s ability to detect changes over time (Sect. 4). To this end, measurements from four co-located (in a 25 cm radius area) DTP systems installed from 23 December 2020 to 17 February 2021, in a silt-dominated soil in Berkeley (CA, USA), were used to estimate thermal diffusivity. In this winter period, the median [min, max] temperature gradient among all the probes recorded between the two top sensors is 5.98 ∘C m−1 [−25.59, 37.07 ∘C m−1]; between the two bottom sensors it is 2.16 ∘C m−1 [1.16, 2.91 ∘C m−1]. The median [min, max] diurnal temperature variation measured by the top sensor at 0.05 m below the ground surface is 1.59 ∘C [0.52, 3.67 ∘C]. Accumulated precipitation every 15 min recorded from the Lawrence Berkeley National Laboratory meteorological station (lat: 37.8771, long: −122.2486) was obtained from https://mesowest.utah.edu (last access: 23 June 2022). Independent measurements of thermal diffusivity were collected on 24 December 2020, after a long dry period that lasted about 30 d, using a thermal-property analyzer (TEMPOS instrument with the SH-3 dual needle; METER Group), which measures conductivity, diffusivity, and volumetric heat capacity with an accuracy of ±10 %. Five measurements were recorded at each 0.10, 0.20, 0.30, 0.40, 0.50, and 0.55 m depth, by repetitively inserting the SH-3 dual needle along the walls of a hand-augured 8 cm diameter hole.
2.4.4 Field study in a cold, unfrozen soil with limited vertical temperature gradient
The second field experiment presented in this study involves the measurement and estimation of soil thermal diffusivity at numerous locations in a discontinuous permafrost environment, and the evaluation of the links between the estimated soil thermal diffusivity values and soil physical properties (Sect. 5). The study site is located along Teller Road about 40 km northwest of Nome, Alaska (64.72∘ N, 165.94∘ W). This site, referred to here as the “Teller site”, is characterized by discontinuous permafrost and a great variety of vegetation types (e.g., tall shrub, dwarf shrub, moss, graminoids; Léger et al., 2019; Uhlemann et al., 2021). Soil thermal diffusivity was inferred from temperature time series recorded in the dry period 7–27 October 2019 from 27 locations distributed over the 2.3 km2 area of the Teller site where during that period soil was entirely unfrozen along the temperature probes (i.e., temperature time series at all depths above 0.5 ∘C). At this time of the year, the median [min, max] temperature gradient among all the probes recorded between the two top sensors was 10.47 ∘C m−1 [−9.22, 35.78 ∘C m−1], and between the two bottom sensors it was as low as 0.63 ∘C m−1 [−0.94, 2.19 ∘C m−1]. The median [min, max] diurnal temperature variation measured by the top sensor at 0.05 m below the ground surface was 0.23 ∘C [0.03, 1.41 ∘C].
Moreover, 92 soil samples at 50 locations were retrieved during the first week of August 2019 and analyzed in the laboratory in order to measure thermal and physical properties (Dafflon et al., 2022a). Thermal conductivity, thermal diffusivity, volumetric heat capacity, and wet bulk density of the soil samples were measured the same day they were collected, using a weighing scale and the thermal-property analyzer (TEMPOS instrument with the SH-3 dual needle; METER Group). Later analyses in the lab involved drying the soil samples at 65 ∘C, recording dry weight, grinding soil to pass a 2 mm sieve, recording weight of the >2 and <2 mm portion, and measuring carbon concentration by combustion. These analyses were used to infer dry bulk density, water content, and carbon density. A total of 20 of these soil samples, retrieved at about 0.075, 0.20, and 0.80 m depth, were collocated with 13 out of the 27 locations monitored with the DTP systems.
Synthetic experiments were performed to investigate how the uncertainty of the inferred thermal diffusivity at multiple depths is affected by various factors, including (1) the amount of temperature gradient at each depth (Sect. 3.1 and 3.3), (2) the length of the sliding time window used in the MCMC inversion (Sect. 3.1), (3) the amount of information contained in temperature trends and fluctuations (Sect. 3.2), (4) the level of noise and bias in the temperature time series (Sect. 3.3), and (5) the geometry of the temperature sensors (Sect. 3.4).
3.1 Impact of time window length
In this section, we assess how the uncertainty of thermal diffusivity estimates is affected by the length of the sliding time window used in the MCMC inversion under two very different types of temperature fields. One temperature field has a strong gradient and diurnal fluctuations that generally occur in summer (Fig. 1a), whereas the other has much more limited gradient and fluctuations that are more typical for the autumn season (Fig. 1h). The diurnal signal in the top boundary condition in summer allows us to get PE ≤ 10 % and posterior RSD at maximum ±5 % on thermal diffusivity estimates in the top 0.25 m, even for very short time windows of 1 d (Fig. 1b, f, and g). Using a time window of 4 or more days allows us to infer thermal diffusivity with a PE ≤ 10 % and a RSD of ±2 % at all depths up to 1 m (Fig. 1c, f and g). A lower diurnal signal in the top boundary condition (between −1 and 2 ∘C m−1), as observed in autumn (Fig. 1h), implies much larger error. Indeed, in this case a time window of at least 21 d is required in order to infer thermal diffusivity at all depths with a PE ≤ 10 % and a RSD of ±5 % (Fig. 1k, m, and n). Moreover, thermal diffusivity estimates in the bottom part of the soil column, where temperature gradients and diurnal fluctuations are smaller than in the upper part, are characterized by higher uncertainty. In general, the optimal length of the time window depends on the information content (e.g., diurnal fluctuations and vertical gradient) in the temperature time series such that the less information content there is the longer the time window needs to be to reliably infer thermal properties.
This synthetic experiment further shows that soil heterogeneity is defined at a spatial resolution equal to the spacing of the temperature sensors. When an interface is located between two sensors, such as the one at 0.42 m in our synthetic soil column, the bulk thermal diffusivity inferred between 0.35 and 0.45 m corresponds to the weighted mean of the thermal diffusivity above and below the soil interface.
3.2 Impact of temperature trend and fluctuations
In this section, we assess the effect of temperature trend versus diurnal fluctuations over time (Fig. 2a–d) on the uncertainty of the estimated thermal diffusivity. For this purpose, we perform the MCMC inversion on various cases, which include (F0) a summer synthetic temperature field (Fig. 2a), (F1) a detrended summer synthetic temperature field (Fig. 2b), (F2) a summer synthetic temperature field with daily and smaller fluctuations smoothed out (Fig. 2c), and (F3) a summer synthetic temperature field without any fluctuations (Fig. 2d). All four cases preserve the same vertical temperature gradient, with a sliding time window of 7 d.
Results indicate that the information contained in the temperature trend has negligible value for estimating soil thermal diffusivity compared to the information provided by the fluctuations (F0 vs. F1, Fig. 2e and f). Indeed, the uncertainty of the thermal diffusivity estimates does not worsen when detrending the time series at each depth. Smoothing out the diurnal and smaller fluctuations can result in percentage errors higher than 10 % (F2 in Fig. 2e) and higher than 30 % when completely removing them (F3 in Fig. 2e). However, we should highlight that we can still get reliable thermal diffusivity estimates (PE ≤ 10 %, RSD ≤ 1 %) in F2 and F3 by increasing the time window length from 7 to 20 d (F2b and F3b in Fig. 2e and f).
3.3 Impact of measurement noise and bias
In this section, we investigate the impact of temperature sensor noise and bias on inferred thermal diffusivity estimates. We consider two distinct scenarios in which temperature data are perturbed with uncorrelated and normally distributed measurement errors, by using different combinations of the mean, μ, and standard deviation, . In the first scenario, soil temperature time series are assumed to be measured by temperature probes with unbiased sensors, μ=0 ∘C, that have different levels of noise equal to 0.01, 0.02, 0.03, 0.05, and 0.1 ∘C. The second scenario considers low levels of noise in the temperature time series, ∘C, and evaluates sensor bias, so that μ is randomly drawn from a zero-mean Gaussian distribution with standard deviation, σc, set equal to 0.005, 0.01, 0.02, and 0.04 ∘C. The MCMC inversion results (Fig. 3a–l) when using the summer temperature data of Fig. 1a show that thermal diffusivity can be inferred with PEs ≤ 10 % (Fig. 3e) and posterior RSDs ≤ 1 % (Fig. 3f) when the noise is no larger than 0.02 ∘C. The percent error may increase up to 20 % with ∘C. Bias in the temperature sensors has an overall larger impact on the PE of thermal diffusivity than a high level of noise in temperature sensors (Fig. 3g–j vs. Fig. 3a–d). Indeed, even with biases applied at each depth drawn randomly from a Gaussian distribution with σc as low as 0.01 ∘C, we can observe biased thermal diffusivity estimates (Fig. 3h) with PE > 10 % (Fig. 3k), and this effect worsens when increasing the bias in the temperature sensors (Fig. 3j).
3.4 Impact of sensor geometry
In this section, we evaluate how thermal diffusivity estimates are affected by the temperature sensor geometry, including the number and position of the sensors, and the effect of potential mispositioning of the sensors in soil. The sensor spatial configuration of the temperature probe used in the present work and described in Sect. 2.4.1 is considered the reference case and is indicated with P0 in Fig. 4.
When increasing the sensor spacing from 0.05 to 0.1 m in the top 0.25 m (P1), the PE of thermal diffusivity estimates can be as high as 20 % in the top 0.1 m, which is likely strengthened by the presence of an interface between the two top sensors (P0 vs. P1 in Fig. 4a). Moreover, a lower resolution in the top part of the probe increases the PE of the estimates in the bottom part (P1 in Fig. 4a). Increasing the sensor spacing from 0.1 to 0.2 m in the bottom part of the probe (i.e., below 0.4 m deep), where the soil is homogenous, does not significantly affect the PE of the thermal diffusivity estimates (P2 in Fig. 4a). A sensor spacing of 0.2 m along the entire probe yields PE larger than 10 % and ≤20 % (P3 in Fig. 4a). The sensor geometries in P0, P1, P2, and P3 do not significantly impact the RSDs of the estimates, which are ≤2 % (Fig. 4b). Decreasing the number of sensors to three or four leads to unreliable thermal diffusivity estimates (P4 and P5 in Fig. 4a), and the corresponding posterior RSDs are higher (Fig. 4b). Indeed, geometries with a significantly smaller number of sensors, such as probes P4 and P5, rely on much fewer data to compute the likelihood and constrain the MCMC inversion.
We further evaluate how an error in positioning a sensor at a specific depth or at a specific distance from other sensors impacts the estimated thermal diffusivity. Note that this risk is more strongly present when using discrete temperature sensors placed manually in the subsurface (Fig. 4c and d) instead of a distributed temperature profiling system. A downwards shift of 0.02 m at the top sensor leads to unreliable thermal diffusivity estimates in the uppermost part of the soil column and to PEs on thermal diffusivity larger than 10 % over the entire soil column (case S1 in Fig. 4c). Shifting the bottom sensor 0.02 m upwards has a smaller impact on the PE of thermal diffusivity estimates (case S2 in Fig. 4c). The main reasons for this are that the soil is homogeneous in the bottom part of the soil column below 0.42 m and that thermal diffusivity is more sensitive to changes in the top temperature time series compared to the bottom ones. Indeed, the top sensors record the temperature time series with the highest content of information (e.g., diurnal and seasonal fluctuations). Similarly, applying a shift at each sensor above 0.5 m depth, randomly drawn from the set −0.01, 0, and 0.01 m, degrades the PE of thermal diffusivity estimates more strongly than when the shifts are applied to the sensors below 0.5 m depth (case S3 vs. S4 in Fig. 4c). If the shifts are applied to all sensors, the PEs are higher than 30 % (case S5 in Fig. 4c). Reducing the number of sensors leads to higher uncertainty in thermal diffusivity estimates (i.e., increase in posterior RSDs), because fewer data are used to constrain the inversion. Instead, a malpositioning of one or more sensors does not affect the RSD of thermal diffusivity estimates (i.e., all ≤1 %, Fig. 4d), since the number of sensors and, hence, data used in the inversion do not change. Moreover, we find that increasing the sensor spacing from 5 to 10 cm in the top 25 cm of soil (probe P1 in Fig. 4a) has a smaller impact on the PE of thermal diffusivity estimates than an error of a few centimeters in positioning these sensors (cases S1 and S3 in Fig. 4c).
The estimation of soil thermal diffusivity from temperature time series was assessed by comparing MCMC inversion results from multiple depth profiles of temperature located close to each other, with independent measurements obtained using the thermal-property analyzer (Sect. 2.4.3). The potential of the MCMC method to infer thermal diffusivity over time was also further evaluated using the sliding time window approach over a 2-month period.
First, we compared the thermal diffusivity estimated with the developed method to those measured with the thermal-property analyzer. The replicated measurements from the thermal-property analyzer, while showing a high variability (large black error bars in Fig. 5a), suggest a fairly homogenous clayey–silty soil, which is in agreement with visual observations from the soil pit dug at the site and values found in the literature (e.g., Andújar Márquez et al., 2016). The thermal diffusivity inferred from the temperature time series shows very similar values, although one of the probes (probe 4) provides slightly higher values than the others at depth greater than 30 cm. Both the measurements from the thermal-property analyzer and the MCMC-inferred thermal diffusivity estimates show similar values and trends in the top 30 cm. Deeper, the measurements from the thermal analyzer suggest a slight decrease in thermal diffusivity that is not visible in the MCMC estimates. The thermal diffusivity values are still comparable, considering that their variability is within the range of typical clay–silt soils and that both methods tend to show the presence of spatial variability in soil thermal properties.
Further, we evaluated the variation in thermal diffusivity over time using the sliding time window approach. The MCMC inversion results from temperature data recorded by the four temperature probes are very similar and show that the estimated thermal diffusivity values in the top 30 cm of soil remain consistently lower than at deeper depths over the entire 2-month period. While the top soil does not show much temporal variability in the soil diffusivity over time, the deeper part shows larger changes, particularly at the time of precipitation events occurring between 22 January and 5 February 2021. These precipitations events are associated with a decrease in temperature (Fig. 5b) and an increase in thermal diffusivity, mainly in the bottom part of the soil column (Fig. 5c). Thermal diffusivities inferred from the sliding time window associated with the large precipitation events show posterior RSDs that are more than double those obtained from the MCMC inversion in the other time windows (Fig. 5d). Moreover, the MCMC inversion provides a worse fit to the temperature data in the time window corresponding to large precipitation events (Fig. 5e). Indeed, we find that the inferred standard deviation of the temperature measurement errors, (Sect. 2.4.3), is typically within 0.01 and 0.02 ∘C, but increases to 0.06 ∘C between 22 January and 5 February 2021 (Fig. 5e). Despite this increase in uncertainty in the bottom part of the soil column during the precipitation events, the observed changes in thermal diffusivity at these depths are consistent with the expected impact of an increase in water content in an initially relatively dry soil (e.g., Farouki, 1981; Arkhangel'skaya and Lukyashchenko, 2018). Unfortunately, we do not have soil moisture data and water level measurements to further investigate the controls on these changes in thermal diffusivity.
The MCMC method is used to estimate unfrozen soil thermal diffusivity at multiple depths and locations across a discontinuous permafrost environment along Teller Road (AK). Improving the estimation of thermal diffusivity in such an environment is critical to potentially improving the estimation of soil physical properties, including organic matter concentration and bulk density, that critically influence carbon cycle dynamics. A particular goal of this study is to infer soil thermal diffusivity at numerous locations and compare it with independent measurements of thermal and physical properties from soil samples, in order to evaluate the strength of these relationships for potential future estimation of soil physical properties from time series of temperature.
MCMC-based estimation of thermal diffusivity every 5 to 10 cm of depth is applied at 27 locations. The 1-month datasets in October (dry period) show very limited median diurnal fluctuations (i.e., less than 0.3 ∘C at 0.05 m below the ground surface), and temperature gradients are low, with values as small as 0.63 ∘C m−1 at 1 m deep (more details in Sect. 2.4). Based on the results from the synthetic experiments (Sect. 3), the use of longer time windows can counterbalance the negative impact of low temperature gradients on the reliability of thermal-property estimates. Hence, for the analysis at the Teller site, we increased the sliding time window from 7 to 10 d and inferred thermal diffusivity up to 0.85 m instead of 1 m depth. The level of fit achieved by the MCMC inversion reflects the accuracy of the temperature probes, since more than 94 % of inferred data errors among all probes and all time windows are ≤0.03 ∘C.
The distribution of estimated thermal diffusivity at each depth from the 27 temperature fields recorded at the Teller site (thin boxplots in Fig. 6) is consistent with the distribution of the thermal diffusivity measured from the 92 soil samples collected across the site and analyzed using the thermal analyzer (bold boxplots in Fig. 6). Estimated and measured thermal diffusivities suggest overall higher values and higher variability at depth larger than about 0.2 m. The thermal diffusivity values in the top 0.2 m are consistent with values observed for the peat layer in the literature (Farouki, 1981; Andújar Márquez et al., 2016). The transition from peat to a more mineral soil occurring at a depth between 0.10 and 0.20 m is also in agreement with the soil visual observations performed during the field campaign. Below 0.15–0.20 m, the soil is mainly composed of clay and silt (yellow shaded area in Fig. 6). The correlation coefficient between measured and inferred thermal diffusivity at 20 co-located locations is 0.79. This correlation is strong, considering that (i) we are comparing thermal diffusivities measured in August 2019 with those inferred from temperature fields recorded in October 2019, (ii) thermal diffusivity from soil samples can be influenced by changes in bulk density occurring during the sampling, and (iii) the scale of the two measurements is different, with soil samples being very local, sparse, and imperfectly co-located with the temperature measurements. Moreover, we note that the temporal variability (i.e., MCMC inference from each sliding time window over the time period of October 2019) in posterior mean thermal diffusivities from the deepest layers (not shown) is characterized by larger variability than those in the shallowest layers. This is in agreement with the results from the synthetic experiments for which smaller temperature gradients are found at deeper depths, causing larger uncertainties in the MCMC-inferred thermal diffusivities.
We further assessed the MCMC-inferred vertically resolved profiles of soil thermal diffusivity by evaluating the relationship between measured and MCMC-inferred thermal diffusivity with other soil thermal properties (i.e., thermal conductivity and volumetric heat capacity) and physical properties (i.e., bulk density, water content, and carbon density) retrieved from laboratory analyses (Fig. 7). For each pair of soil properties, we compute the distance correlation (dc) coefficient in order to capture linear and nonlinear dependencies. The strongest distance correlation coefficient (dc>0.9) is found between thermal diffusivity (measured from soil samples and inferred with the MCMC method) and thermal conductivity, wet bulk density, and dry bulk density measured from soil samples (Fig. 7a, c, and d). Also, thermal diffusivities inferred with the MCMC method and measured from soil samples are both similarly correlated (dc equals 0.78 and 0.76, respectively) to carbon density (Fig. 7f). The lowest correlation is observed between volumetric heat capacity and measured (dc=0.72) and estimated (dc=0.49) thermal diffusivity.
A large amount of the variability in each of the thermal and physical properties is linked to the vertical heterogeneity of the soil. All the shallow soil samples within the top 0.2 m are composed of at least 50 % water (Fig. 7e), are rich in organic carbon, and are characterized by lower thermal diffusivity values (Fig. 7f). Samples from depths between approximately 0.2 and 0.5 m show a large range of variability in all the properties. Samples from deeper than 0.5 m tend to show the highest values in wet and dry bulk density and in thermal diffusivity.
In this study, we developed a methodology for estimating soil thermal diffusivity and its uncertainty at an unprecedented vertical spatial resolution, at multiple locations, and over time, using depth-resolved time series of soil temperature, thermal modeling, and Bayesian inference. Through synthetic and field experiments, we assess how various sources of uncertainty impact thermal diffusivity estimates, and we evaluate the method’s performance and its potential.
6.1 Impact of environmental conditions, sensor characteristics, and deployment geometry on thermal diffusivity estimates
Results show that environmental conditions (i.e., temperature gradient, period and amplitude of temperature fluctuations, temperature trend) as well as measurement strategies affect the final amount of information contained in the temperature time series and, therefore, the quality of thermal diffusivity estimates. The outcome from our study is in agreement with the work of Jafarov et al. (2014), Rodler et al. (2019), and Jafarov et al. (2020), who concluded that the higher the temperature gradient, the better the results. In particular, we found that temperature gradients smaller than ±2 ∘C m−1 (Fig. 1) can be particularly problematic and lead to unreliable thermal-property estimates (i.e., percentage errors larger than 30 %), if an insufficiently long time window is used. Low information content in the temperature time series can be partly compensated for by increasing the size of the time window used to estimate soil thermal diffusivity. For example, by lengthening the time window (i.e., from 7 to 20 d in our synthetic experiment), the amount of data used to compute the likelihood in the MCMC inversion increases, and we can recover reliable thermal diffusivity values up to 1 m deep under temperature gradients smaller than 2 ∘C m−1 (Fig. 1h) or when diurnal fluctuations are absent (Fig. 2d). However, we need to keep in mind that in real case studies, the time window should be short enough to limit the influence of physical processes not represented in the heat diffusion model (e.g., advection, phase change). This is why in the field study at the Teller site we increased the time window from 7 to 10 d and investigated thermal diffusivity up to 0.85 m instead of 1 m depth.
Besides environmental conditions, additional sources of uncertainty for thermal diffusivity estimates are those related to acquisition strategies, including sensor quality, calibration, and deployment geometry. In Sect. 3.3, we investigated the impact of temperature fields perturbed with different levels of noise and bias. Results have shown that the use of sensors with even a small amount of bias is more detrimental to thermal diffusivity estimates than having sensors with some random noise (level of noise ≥ 0.05 ∘C), because the former might ultimately lead to a misleading interpretation of the soil layering (Fig. 3a–d vs. Fig. 3g–j). Indeed, the availability of high-accuracy sensors, with a bias defined by a standard deviation of 0.01 ∘C or less, is fundamental to ensure percentage errors on thermal diffusivity ≤ 10 %. This accuracy requirement underlines the challenge of this method and the importance of sensor-calibration approaches to increase sensor accuracy as much as possible.
Furthermore, results show that the deployment geometry is critical to potentially capturing heterogeneity present in the soil, and particularly the layer boundaries. The heterogeneity underlines the importance of inferring soil diffusivity with high spatial resolution, as done in this study. Note that the impact of sensor geometry on thermal diffusivity estimates is expected to be even larger in more heterogenous soils. Finally, while several environmental factors and measurement strategies strongly influence the estimate accuracy, the potential error in mispositioning a sensor in soil, which is less likely using a temperature probe than individual sensors, can easily lead to percentage errors in thermal diffusivity estimates larger than 30 % (Fig. 4).
In summary, we were able to identify five important aspects to consider at the design stage of temperature sensors: (1) having more sensors closer to the soil surface where temperature signal has the highest content of information (e.g., diurnal and seasonal fluctuations); (2) using temperature probes (i.e., sensors anchored on the same support) instead of discrete temperature sensors placed manually in the subsurface, which imply higher uncertainty in the distance between sensors; (3) ensuring a vertical spatial resolution that captures the soil layering as much as possible and thus decreases the sensor spacing in vertically heterogeneous environments; (4) collecting high-accuracy measurements; and (5) taking into consideration the environmental conditions at the sites of interest. Indeed, we concluded that, over a 1-month period with median diurnal fluctuations ≥ 1.5 ∘C at 5 cm below the ground surface, temperature gradients > 2 ∘C m−1 at each depth, and a sliding time window of at least 4 d, the proposed method provides reliable depth-resolved thermal diffusivity estimates, with a percentage error ≤ 10 % and posterior relative standard deviations ≤ 5 % up to 1 m depth. Reliable thermal diffusivity under such environmental conditions also requires temperature sensors with a level of noise ≤ 0.02 ∘C, a bias defined by a standard deviation of 0.01 ∘C or less, and a positioning accuracy of a few millimeters or less (i.e., temperature probe). These results were obtained from challenging synthetic experiments in which the boundary conditions used for the thermal modeling are soil temperature time series recorded in Alaska where temperature gradients and diurnal temperature variations are typically lower than in other regions. Therefore, the proposed method has the potential to provide reliable depth-resolved thermal diffusivity estimates at depths deeper than 1 m and over time windows shorter than 4 d.
6.2 Spatial variability of thermal diffusivity and link with soil physical properties
The reliability of the developed approach has been demonstrated for in situ estimation of vertically resolved thermal diffusivity profiles at field scale, including a site in Berkeley (CA, USA) and along Teller Road, northwest of Nome (AK, USA). At both sites, the estimated values compare well with independent measurements obtained using a thermal-property analyzer. Still, they do not match perfectly, because they are measured at different spatial scales. Moreover, soil temperatures recorded at the Teller site in autumn are characterized by very low temperature gradients and diurnal fluctuations, which make thermal diffusivity estimates extremely challenging. Overall, the obtained thermal diffusivity values were consistent with those reported in the literature for the various soil types encountered at these sites.
The field results underline the value of thermal diffusivity in quantifying soil physical properties, as well as challenges associated with interpreting spatial and temporal variations in thermal diffusivity. For example, the increase in water content in a particular dry material (<15 %) can drive an increase in thermal diffusivity, while a similar increase in a wetter material (>15 %) will barely impact the thermal diffusivity or drive a decrease in it (Farouki, 1981). This complexity explains why we observe a negative correlation between thermal diffusivity and soil moisture at the site along Teller Road, AK (Fig. 7e), while the thermal diffusivity of deep soil at the Berkeley site increases after a rain event (Fig. 5). In the latter case, the soil was drier in the deep subsurface as a consequence of a long dry period that lasted about 1 month.
The importance of estimating thermal diffusivity lies in the fact that this soil property shows strong correlation (i.e., between 0.78 and 0.92) with thermal conductivity, wet and dry bulk density, water content, and carbon density (Fig. 7). These results confirm previous laboratory studies (Farouki, 1981; Arkhangel'skaya and Lukyashchenko, 2018; Mengistu et al., 2017; Ochsner and Baker, 2008) investigating these relationships and offer new datasets to improve the integrated estimation of these properties. The recent work from Zhu et al. (2019) highlights the value of temperature measurements in quantifying thermal diffusivity and further evaluating the dominant effect of organic carbon on permafrost dynamics.
6.3 Temporal variability of thermal diffusivity
Finally, results from this study show promise in using temperature time series to estimate temporal changes in soil diffusivity. This is the first time, to our knowledge, that this concept has been introduced in environmental science at such a vertical spatial resolution and in a Bayesian framework. Tabbagh et al. (2017) have used a similar sliding time window approach to evaluate the temporal variability of thermal diffusivity, but for a single layer and with a deterministic inverse method. The assumptions of our method, based on which thermal diffusivity is constant in each time window, and on purely heat conduction processes taking place, can be violated if soil wetness changes considerably on timescales smaller than the length of the time window. However, we can still detect these time periods by looking at the MCMC outputs and at the inferred data errors, as shown in this study: the precipitation event at the Berkeley site caused at least a 2-fold increase in the MCMC-inferred thermal diffusivity estimates (Fig. 5c) and corresponding posterior relative standard deviations (Fig. 5d), as well as a larger inferred data error (Fig. 5e). These results indicate that the sliding time window approach can be used to detect when the assumption of a heat conduction model is not valid due to the presence of advective heat transfer and/or large change in water content over the used time windows. The detection of such processes is important to identify when the thermal diffusivity values estimated using a heat conduction model are reliable or not, as well as to uncover the presence of significant soil hydrological fluxes. While incorporating advection and/or freeze–thaw processes into models is feasible (e.g., Zhang et al., 2016; Zhao et al., 2016; Tabbagh et al., 2017; Jafarov et al., 2012; Romanovsky and Osterkamp, 2000; Nicolsky et al., 2009), it generally requires additional hydrological constraints. Thus the developed approach is particularly valuable for application at locations where hydrological measurements are generally absent or uncertain.
In this work, we have presented a parameter estimation approach based on the combination of thermal modeling (i.e., heat diffusion equation), sliding time windows, Bayesian inference, and MCMC simulation. The method enables us to estimate soil thermal diffusivity and its uncertainty from solely depth-resolved temperature time series at an unprecedented vertical spatial resolution (i.e., 5 to 10 cm resolution up to 1 m deep), at multiple locations, and over time. Results obtained through synthetic experiments indicate that reliable depth-resolved thermal diffusivity estimates with percentage errors ≤ 10 % and posterior relative standard deviations ≤ 5 % up to 1 m depth can be obtained in the presence of median diurnal fluctuations ≥ 1.5 ∘C at 5 cm below the ground surface, temperature gradients > 2 ∘C m−1, and a sliding time window of at least 4 d. Moreover, results indicate that reliable thermal diffusivity estimates in the above environmental conditions also require sensors to be spaced precisely (with accuracy to a few millimeters), with a level of noise ≤ 0.02 ∘C, and with a bias defined by a standard deviation ≤ 0.01 ∘C. These findings represent a critical step to guide the development of cost-effective and reliable methodologies to estimate soil thermal properties from time series of temperature at numerous locations. Importantly, the application of the method to field data shows results consistent with the synthetic experiments, as well as with independent measurements performed on soil samples. Finally, the field studies show promise in the application of a sliding time window to estimate temporal changes in soil thermal diffusivity which can be potentially used to capture the corresponding changes in water content.
The strength of our methodology lies in the fact that thermal diffusivity can be quantified from depth-resolved temperature time series, at numerous locations. This is of great value, considering that thermal diffusivity is strongly linked with thermal conductivity, wet and dry bulk density, water content, and carbon density, as shown by this and previous works. Hence, our approach opens the way for future research on the development of petrophysical relationships that can be integrated within the Bayesian inversion framework and used to derive from soil thermal diffusivity, inferred from temperature time series, the fraction of soil components, and in particular carbon density in organic-rich environments. The carbon density is a key property in the quantification of soil respiration rate. Similarly, an improved quantification of soil thermal parameters is needed to better parametrize ecosystem models, where soil thermal parameters modulate the effect of weather forcing on subsurface thermal and biogeochemical fluxes.
The forward model written in MATLAB that computes the 1D diffusion equation in a heterogeneous medium (Eq. 3 in this paper) was deposited in the NGEE Arctic (Next-Generation Ecosystem Experiments) portal at https://doi.org/10.5440/1433255 (Brunetti et al., 2022). The complete MATLAB code used for the inverse modeling based on MCMC is a proprietary software and more details can be found in Vrugt (2016). A similar implementation has been written in Python and is freely accessible at Shockley et al. (2021).
The data presented in this study are available from the NGEE Arctic (Next-Generation Ecosystem Experiments) data portal at https://doi.org/10.5440/1433255 (Brunetti et al., 2022) and https://doi.org/10.5440/1860505 (Dafflon et al., 2022a).
CB performed the formal analysis, methodology development, validation, visualization, writing and data collection. JL performed data collection. SW performed sensor development. SU performed data collection. IS performed data collection. PM performed sensor development. BD performed conceptualization, funding acquisition, review and editing, and data collection.
The contact author has declared that neither they nor their co-authors have any competing interests.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
We thank Vladimir Romanovsky for providing soil temperatures used for the synthetic examples. We thank Jasper Vrugt for making DREAM(ZS) available by providing an initial version of the code.
This material is based upon work supported by Next-Generation Ecosystem Experiments (NGEE Arctic), funded by the US Department of Energy Office of Science Office of Biological and Environmental Research (award no. DE-AC02-05CH11231).
This paper was edited by Michael Krautblatter and reviewed by Elchin Jafarov and one anonymous referee.
Abu-Hamdeh, N. H.: Thermal properties of soils as affected by density and water content, Biosyst. Eng., 86, 97–102, https://doi.org/10.1016/S1537-5110(03)00112-0, 2003. a, b
Al Nakshabandi, G. and Kohnke, H.: Thermal conductivity and diffusivity of soils as related to moisture tension and other physical properties, Agricult,. Meteorol., 2, 271–279, https://doi.org/10.1016/0002-1571(65)90013-0, 1965. a, b
An, K., Wang, W., Zhao, Y., Huang, W., Chen, L., Zhang, Z., Wang, Q., and Li, W.: Estimation from soil temperature of soil thermal diffusivity and heat flux in sub-surface layers, Bound.-Lay. Meteorol., 158, 473–488, https://doi.org/10.1007/s10546-015-0096-7, 2016. a
Andújar Márquez, J. M., Martínez Bohórquez, M. Á., and Gómez Melgar, S.: Ground thermal diffusivity calculation by direct soil temperature measurement. Application to very low enthalpy geothermal energy systems, Sensors, 16, 306, https://doi.org/10.3390/s16030306, 2016. a, b, c, d, e, f, g
Arkhangel'skaya, T.: Parameterization and mathematical modeling of the dependence of soil thermal diffusivity on the water content, Euras. Soil Sci., 42, 162–172, https://doi.org/10.1134/S1064229309020070, 2009. a, b
Arkhangel'skaya, T. and Lukyashchenko, K.: Estimating soil thermal diffusivity at different water contents from easily available data on soil texture, bulk density, and organic carbon content, Biosyst. Eng., 168, 83–95, https://doi.org/10.1016/j.biosystemseng.2017.06.011, 2018. a, b, c, d
Beardsmore, G., Egan, S., and Sandiford, M.: A Fourier Spectral Method to Measure the Thermal Diffusivity of Soil, Geotech. Test. J., 43, 3, https://doi.org/10.1520/GTJ20180300, 2020. a, b, c
Brunetti, C., Lamb, J., Wielandt, S., Uhlemann, S., Shirley, I., McClure, P., and Dafflon, B.: Probabilistic estimation of depth-resolved profiles of soil thermal diffusivity from temperature time series: Supporting Data, Next Generation Ecosystem Experiments Arctic Data Collection, Oak Ridge National Laboratory, US Department of Energy, Oak Ridge, Tennessee, USA [code and data set], https://doi.org/10.5440/1433255, 2022. a, b, c
Carson, J. E.: Analysis of soil and air temperatures by Fourier techniques, J. Geophys. Res., 68, 2217–2232, https://doi.org/10.1029/JZ068i008p02217, 1963. a
Choi, W., Kikumoto, H., Choudhary, R., and Ooka, R.: Bayesian inference for thermal response test parameter estimation and uncertainty assessment, Appl. Energ., 209, 306–321, https://doi.org/10.1016/j.apenergy.2017.10.034, 2018. a
Courant, R., Friedrichs, K., and Lewy, H.: Über die partiellen Differenzengleichungen der mathematischen Physik, Math. Annal., 100, 32–74, 1928. a
Dafflon, B., Brunetti, C., Lamb, J., Uhlemann, S., Shirley, I., McClure, P., and Wang, C.: Soil physical and thermal properties from soil samples at multiple locations at Teller Road Site, Seward Peninsula, Alaska, 2019, Next Generation Ecosystem Experiments Arctic Data Collection, Oak Ridge National Laboratory, US Department of Energy, Oak Ridge, Tennessee, USA [data set], https://doi.org/10.5440/1860505, 2022a. a, b
Dafflon, B., Wielandt, S., Lamb, J., McClure, P., Shirley, I., Uhlemann, S., Wang, C., Fiolleau, S., Brunetti, C., Akins, F. H., Fitzpatrick, J., Pullman, S., Busey, R., Ulrich, C., Peterson, J., and Hubbard, S. S.: A distributed temperature profiling system for vertically and laterally dense acquisition of soil and snow temperature, The Cryosphere, 16, 719–736, https://doi.org/10.5194/tc-16-719-2022, 2022b. a, b
De Simon, L., Iglesias, M., Jones, B., and Wood, C.: Quantifying uncertainty in thermophysical properties of walls by means of Bayesian inversion, Energ. Build., 177, 220–245, https://doi.org/10.1016/j.enbuild.2018.06.045, 2018. a
Farouki, O. T.: Thermal properties of soils, vol. 81-1, Cold Regions Research and Engineering Lab, Hanover, NH, USA, https://apps.dtic.mil/sti/pdfs/ADA111734.pdf (last access: 23 June 2022), 1981. a, b, c, d, e, f, g, h, i
Gao, Z., Lenschow, D. H., Horton, R., Zhou, M., Wang, L., and Wen, J.: Comparison of two soil temperature algorithms for a bare ground site on the Loess Plateau in China, J. Geophys. Res.-Atmos., 113, D18105, https://doi.org/10.1029/2008JD010285, 2008. a, b
Gao, Z., Tong, B., Horton, R., Mamtimin, A., Li, Y., and Wang, L.: Determination of desert soil apparent thermal diffusivity using a conduction-convection algorithm, J. Geophys. Res.-Atmos., 122, 9569–9578, https://doi.org/10.1002/2017JD027290, 2017. a, b, c, d
Gelman, A. and Rubin, D. B.: Inference from iterative simulation using multiple sequences, Stat. Sci., 7, 457–472, https://doi.org/10.1214/ss/1177011136, 1992. a
Gelman, A., Roberts, G. O., and Gilks, W. R.: Efficient Metropolis jumping rules, Bayesian Statistics, http://people.ee.duke.edu/~lcarin/baystat5.pdf (last access: 23 June 2022), 1996. a
Gilks, W. R., Richardson, S., and Spiegelhalter, D. J.: Introducing Markov chain Monte, Markov chain Monte Carlo in practice, Chapman and Hall/CRC, p. 1, https://doi.org/10.1201/b14835, 1995. a
Gnanasekaran, N. and Balaji, C.: Markov Chain Monte Carlo (MCMC) approach for the determination of thermal diffusivity using transient fin heat transfer experiments, Int. J. Therm. Sci, 63, 46–54, https://doi.org/10.1016/j.ijthermalsci.2012.07.005, 2013. a
Hillel, D.: Introduction to soil physics, Academic Press, New York, https://doi.org/10.1016/C2009-0-03052-9, 1982. a, b
Hinkel, K. M.: Estimating seasonal values of thermal diffusivity in thawed and frozen soils using temperature time series, Cold Reg. Sci. Technol., 26, 1–15, https://doi.org/10.1016/S0165-232X(97)00004-9, 1997. a
Horton, R., Wierenga, P., and Nielsen, D.: Evaluation of methods for determining the apparent thermal diffusivity of soil near the surface, Soil Sci. Soc. Am. J., 47, 25–32, https://doi.org/10.2136/sssaj1983.03615995004700010005x, 1983. a
Huang, Y., Jiang, J., Ma, S., Ricciuto, D., Hanson, P. J., and Luo, Y.: Soil thermal dynamics, snow cover, and frozen depth under five temperature treatments in an ombrotrophic bog: Constrained forecast with data assimilation, J. Geophys. Res.-Biogeo., 122, 2046–2063, https://doi.org/10.1002/2016JG003725, 2017. a
Irvine, D. J., Briggs, M. A., Lautz, L. K., Gordon, R. P., McKenzie, J. M., and Cartwright, I.: Using diurnal temperature signals to infer vertical groundwater-surface water exchange, Groundwater, 55, 10–26, https://doi.org/10.1111/gwat.12459, 2017. a
Jafarov, E., Nicolsky, D., Romanovsky, V., Walsh, J., Panda, S., and Serreze, M.: The effect of snow: How to better model ground surface temperatures, Cold Reg. Sci. Technol., 102, 63–77, https://doi.org/10.1016/j.coldregions.2014.02.007, 2014. a
Jafarov, E. E., Marchenko, S. S., and Romanovsky, V.: Numerical modeling of permafrost dynamics in Alaska using a high spatial resolution dataset, The Cryosphere, 6, 613–624, https://doi.org/10.5194/tc-6-613-2012, 2012. a, b
Jafarov, E. E., Harp, D. R., Coon, E. T., Dafflon, B., Tran, A. P., Atchley, A. L., Lin, Y., and Wilson, C. J.: Estimation of subsurface porosities and thermal conductivities of polygonal tundra by coupled inversion of electrical resistivity, temperature, and moisture content data, The Cryosphere, 14, 77–91, https://doi.org/10.5194/tc-14-77-2020, 2020. a, b
Jong van Lier, Q. D. and Durigon, A.: Soil thermal diffusivity estimated from data of soil temperature and single soil component properties, Revista Brasileira de Ciência do Solo, 37, 106–112, https://doi.org/10.1590/S0100-06832013000100011, 2013. a, b
Kaipio, J. P. and Fox, C.: The Bayesian framework for inverse problems in heat transfer, Heat Trans. Eng., 32, 718–753, https://doi.org/10.1080/01457632.2011.525137, 2011. a
Kim, H., Lee, J.-Y., and Kim, K.: Estimation of thermal diffusivity of soils in Antarctica using temperature time series data, Episod. J. Int. Geosci., 42, 245–252, https://doi.org/10.18814/epiiugs/2019/019020, 2019. a, b
Koven, C. D., Riley, W. J., and Stern, A.: Analysis of permafrost thermal dynamics and response to climate change in the CMIP5 Earth System Models, J. Climate, 26, 1877–1900, https://doi.org/10.1175/JCLI-D-12-00228.1, 2013. a
Krzeminska, D. M., Steele-Dunne, S. C., Bogaard, T. A., Rutten, M. M., Sailhac, P., and Geraud, Y.: High-resolution temperature observations to monitor soil thermal properties as a proxy for soil moisture condition in clay-shale landslide, Hydrol. Process., 26, 2143–2156, https://doi.org/10.1002/hyp.7980, 2012. a, b, c, d
Laloy, E. and Vrugt, J. A.: High-dimensional posterior exploration of hydrologic models using multiple-try DREAM(ZS) and high-performance computing, Water Resour. Res., 48, W01526, https://doi.org/10.1029/2011WR010608, 2012. a
Léger, E., Dafflon, B., Robert, Y., Ulrich, C., Peterson, J. E., Biraud, S. C., Romanovsky, V. E., and Hubbard, S. S.: A distributed temperature profiling method for assessing spatial variability in ground temperatures in a discontinuous permafrost region of Alaska, The Cryosphere, 13, 2853–2867, https://doi.org/10.5194/tc-13-2853-2019, 2019. a
Mengistu, A. G., van Rensburg, L. D., and Mavimbela, S. S.: The effect of soil water and temperature on thermal properties of two soils developed from aeolian sands in South Africa, Catena, 158, 184–193, https://doi.org/10.1016/j.catena.2017.07.001, 2017. a, b
Metropolis, N., Rosenbluth, A. W., Rosenbluth, M. N., Teller, A. H., and Teller, E.: Equation of state calculations by fast computing machines, J. Chem. Phys., 21, 1087–1092, https://doi.org/10.1063/1.1699114, 1953. a
Nicolsky, D. and Romanovsky, V. E.: Modeling long-term permafrost degradation, J. Geophys. Res.-Earth, 123, 1756–1771, https://doi.org/10.1029/2018JF004655, 2018. a
Nicolsky, D., Romanovsky, V., and Panteleev, G.: Estimation of soil thermal properties using in-situ temperature measurements in the active layer and permafrost, Cold Reg. Sci. Technol., 55, 120–129, https://doi.org/10.1016/j.coldregions.2008.03.003, 2009. a, b, c, d, e
Ochsner, T. E. and Baker, J. M.: In situ monitoring of soil thermal properties and heat flux during freezing and thawing, Soil Sci. Soc. Am. J., 72, 1025–1032, https://doi.org/10.2136/sssaj2007.0283, 2008. a, b
Ochsner, T. E., Horton, R., and Ren, T.: A new perspective on soil thermal properties, Soil Sci. Soc. Am. J., 65, 1641–1647, https://doi.org/10.2136/sssaj2001.1641, 2001. a, b
Oliva, M. and Fritz, M.: Permafrost degradation on a warmer Earth: Challenges and perspectives, Curr. Opin. Environ. Sci. Health, 5, 14–18, https://doi.org/10.1016/j.coesh.2018.03.007, 2018. a
Petter Langtangen, H. and Linge, S.: Finite difference computing with PDEs: a modern software approach, Springer Nature, https://doi.org/10.1007/978-3-319-55456-3, 2017. a
Praprotnik, M., Sterk, M., and Trobec, R.: Inhomogeneous heat-conduction problems solved by a new explicit finite difference scheme, Int. J. Pure Appl. Math., 13, 275–292, 2004. a
Rajeev, P. and Kodikara, J.: Estimating apparent thermal diffusivity of soil using field temperature time series, Geomech. Geoeng., 11, 28–46, https://doi.org/10.1080/17486025.2015.1006266, 2016. a, b
Rasmussen, L. H., Zhang, W., Hollesen, J., Cable, S., Christiansen, H. H., Jansson, P.-E., and Elberling, B.: Modelling present and future permafrost thermal regimes in Northeast Greenland, Cold Reg. Sci. Technol., 146, 199–213, https://doi.org/10.1016/j.coldregions.2017.10.011, 2018. a
Renard, P., Alcolea, A., and Gingsbourger, D.: Stochastic versus deterministic approaches, in: Environmental Modelling: Finding Simplicity in Complexity, 2nd Edn., edited by: Wainwright, J. and Mulligan, M., Wiley Online Library, 133–149, https://doi.org/10.1002/9781118351475.ch8, 2013. a
Robert, C. and Casella, G.: Monte Carlo statistical methods, Springer Science & Business Media, https://doi.org/10.1007/978-1-4757-4145-2, 2013. a
Rodler, A., Guernouti, S., and Musy, M.: Bayesian inference method for in situ thermal conductivity and heat capacity identification: Comparison to ISO standard, Construct. Build. Mater., 196, 574–593, https://doi.org/10.1016/j.conbuildmat.2018.11.110, 2019. a, b, c
Romanovsky, V., Cable, W., and Dolgikh, K.: Soil Temperature and Moisture, Teller Road Mile Marker 27, Seward Peninsula, Alaska, beginning 2016, Next Generation Ecosystem Experiments Arctic Data Collection, Oak Ridge National Laboratory, US Department of Energy, Oak Ridge, Tennessee, USA [data set], https://doi.org/10.5440/1581437, 2020. a
Romanovsky, V. E. and Osterkamp, T.: Effects of unfrozen water on heat and mass transport processes in the active layer and permafrost, Permafrost Periglac. Process., 11, 219–239, https://doi.org/10.1002/1099-1530(200007/09)11:3<219::AID-PPP352>3.0.CO;2-7, 2000. a
Shockley, E. M., Vrugt, J. A., and Lopez, C. F.: PyDREAM: high-dimensional parameter inference for biological models in python, Bioinformatics, 34, 695–697, https://doi.org/10.1093/bioinformatics/btx626, 2018. a
Somasundharam, S. and Reddy, K.: Inverse estimation of thermal properties using Bayesian inference and three different sampling techniques, Invers. Problem. Sci. Eng., 25, 73–88, https://doi.org/10.1080/17415977.2016.1138946, 2017. a
Steele-Dunne, S., Rutten, M., Krzeminska, D., Hausner, M., Tyler, S., Selker, J., Bogaard, T., and Van De Giesen, N.: Feasibility of soil moisture estimation using passive distributed temperature sensing, Water Resour. Res., 46, W03534, https://doi.org/10.1029/2009WR008272, 2010. a
Tabbagh, A., Cheviron, B., Henine, H., Guérin, R., and Bechkit, M.-A.: Numerical determination of vertical water flux based on soil temperature profiles, Adv. Water Resour., 105, 217–226, https://doi.org/10.1016/j.advwatres.2017.05.003, 2017. a, b, c, d, e, f
Tong, B., Gao, Z., Horton, R., Li, Y., and Wang, L.: An empirical model for estimating soil thermal conductivity from soil water content and porosity, J. Hydrometeorol., 17, 601–613, https://doi.org/10.1175/JHM-D-15-0119.1, 2016. a, b
Tran, A. P., Dafflon, B., and Hubbard, S. S.: Coupled land surface–subsurface hydrogeophysical inverse modeling to estimate soil organic carbon content and explore associated hydrological and thermal dynamics in the Arctic tundra, The Cryosphere, 11, 2089–2109, https://doi.org/10.5194/tc-11-2089-2017, 2017. a, b, c
Uhlemann, S., Dafflon, B., Peterson, J., Ulrich, C., Shirley, I., Michail, S., and Hubbard, S.: Geophysical monitoring shows that spatial heterogeneity in thermohydrological dynamics reshapes a transitional permafrost system, Geophys. Res. Lett., 48, e2020GL091149, https://doi.org/10.1029/2020gl091149, 2021. a
Vrugt, J. A.: Markov chain Monte Carlo simulation using the DREAM software package: Theory, concepts, and MATLAB implementation, Environ. Model. Softw., 75, 273–316, https://doi.org/10.1016/j.envsoft.2015.08.013, 2016. a, b, c
Xie, X., Lu, Y., Ren, T., and Horton, R.: An empirical model for estimating soil thermal diffusivity from texture, bulk density, and degree of saturation, J. Hydrometeorol., 19, 445–457, https://doi.org/10.1175/JHM-D-17-0131.1, 2018. a, b
Zhang, M., Wen, Z., Xue, K., Chen, L., and Li, D.: A coupled model for liquid water, water vapor and heat transport of saturated–unsaturated soil in cold regions: model formulation and verification, Environ. Earth Sci., 75, 701, https://doi.org/10.1007/s12665-016-5499-3, 2016. a
Zhao, Y., Si, B., He, H., Xu, J., Peth, S., and Horn, R.: Modeling of coupled water and heat transfer in freezing and thawing soils, Inner Mongolia, Water, 8, 424, https://doi.org/10.3390/w8100424, 2016. a
Zhu, D., Ciais, P., Krinner, G., Maignan, F., Puig, A. J., and Hugelius, G.: Controls of soil organic matter on soil thermal dynamics in the northern high latitudes, Nat. Commun., 10, 1–9, https://doi.org/10.1038/s41467-019-11103-1, 2019. a, b
- Abstract
- Introduction
- Theory and method
- Method assessment through synthetic experiments
- Assessment of repeatability and temporal changes of in situ thermal diffusivity estimates
- Estimation of soil thermal diffusivity in a cold, unfrozen soil with limited vertical temperature gradient
- Discussion
- Conclusions
- Code availability
- Data availability
- Author contributions
- Competing interests
- Disclaimer
- Acknowledgements
- Financial support
- Review statement
- References
- Abstract
- Introduction
- Theory and method
- Method assessment through synthetic experiments
- Assessment of repeatability and temporal changes of in situ thermal diffusivity estimates
- Estimation of soil thermal diffusivity in a cold, unfrozen soil with limited vertical temperature gradient
- Discussion
- Conclusions
- Code availability
- Data availability
- Author contributions
- Competing interests
- Disclaimer
- Acknowledgements
- Financial support
- Review statement
- References