In shallow geothermal systems, the main equivalent underground thermal properties are commonly calculated with a thermal response test (TRT). This is a borehole heat exchanger production test where the temperature of a heat transfer fluid is recorded over time at constant power heat injection/extraction. The equivalent thermal parameters (thermal conductivity, heat capacity) are simply deduced from temperature data regression analysis that theoretically is a logarithmic function in the time domain, or else a linear function in the log-time domain. By interpreting the recorded temperatures as a regionalized variable whose drift is the regression function, in both cases the formal problem is a linear estimation of the mean. If the autocorrelation function (variogram, covariance) of residuals is known, coefficient variance can be directly deduced. Coefficient estimates are independent of the drift form adopted, and the residuals are the same in the same points. The random function is different in the time domain, however, and in the log-time domain. In fact, residual The work was presented to GeoEnv 2014 Conference on Geostatistics for Environmental Applications, held in Paris, 9-11/07/2014.