Estimation of In Situ Heat Capacity and Thermal Diffusivity from Undisturbed Ground Temperature Profile Measured in Ground Heat Exchangers

: Undisturbed ground temperature (UGT), thermal conductivity (TC) and heat capacity (HC) are essential parameters to design geothermal heat pumps and underground thermal energy storage systems, particularly borehole heat exchangers (BHE). However, field methods to assess the thermal state and properties of the subsurface are costly and time consuming. Moreover, HC is often not evaluated in situ but arbitrarily selected from literature considering the geological materials intercepted by boreholes. This work proposes an original empirical approach to reproduce a UGT profile and estimate in situ thermal diffusivity (TD) and HC in the scope of conventional thermal response tests (TRTs). Empirical equations were developed to reproduce a UGT profile measured along a BHE. Experimental coefficients are found with a non-linear least square solver optimization and used to calculate the damping depth, TD and HC. The suggested heat tracing method was verified and validated against other field methods demonstrating to be fast and reliable. The novelty of this new empirical approach relies on the use of a single temperature profile providing a simple way to better assess subsurface thermal properties.


Introduction
Thermal response tests (TRTs) made in borehole heat exchangers (BHEs) and analyzed with the infinite line source equation are commonly used to infer in situ thermal conductivity (TC) [1][2][3].However, the heterogeneity of the ground and the presence of groundwater flow are often neglected to simplify the analysis.Furthermore, estimation of in situ thermal diffusivity (TD) and heat capacity (HC) are also needed when TRTs are combined with the infinite line source equation to reproduce observed temperatures [4].
TRTs analyzed with numerical models or the moving line source equation, for example, are advanced options to characterize the thermal and hydrogeological properties of the subsurface [2,3].However, developing such numerical or analytical models is a complex task and only provides information on TC and not HC.Core pieces or chip samples can be analyzed in the laboratory in dry or water saturated conditions to better assess thermal storage properties and reduce uncertainty [4][5][6].However, samples are commonly disturbed and analysis may not be representative of the in situ conditions.Another common simple practice is to arbitrarily select a literature-based in situ HC that matches the geological description of materials found in a borehole [7][8][9][10].However, a study carried out by Giordano et al. [11] revealed that typical uncertainty associated to the in situ TD is about ±40% when a conventional TRT is analyzed with an approximated HC for common geological materials (1.50 to 3.20 MJ m −3 K −1 ).All these examples highlight uncertainty that prevails about HC when performing a conventional TRT [4,5].Therefore, alternative field methods to infer heat storage properties are still needed.Oscillatory TRT (OTRT) proposed and discussed by, for instance, Oberdorfer [12] and applied by Giordano et al. [11], have the potential to evaluate the in situ TD and HC.This methodology is an important step taken towards the in situ evaluation of this thermal property.However, analyzing the oscillating temperature response is complex and the method needs improvement to significantly reduce uncertainty (currently on the order of ±15%).Thus, improving the assessment of ground thermal properties with simple and efficient methods to infer the in situ TD and HC is necessary.Such improvements can help to better design geothermal systems relying on BHEs.Ultimately, this can have a significant positive impact on CO2 emission reduction and provide advantages to develop green energy alternatives [13].
Another available approach for the assessment of TD and HC is the analysis of ground temperature profiles to characterize the thermo-hydraulic properties of the geological materials [3,4,11,[14][15][16]. Ground temperature profiles can be acquired with a submersible probe that is lowered into the BHE during or after a TRT ⦋15,16⦌.These ground temperature profiles can be used to improve the evaluation of in situ TC and estimate horizontal groundwater fluxes [1,3,14,16].Moreover, ground temperature profiles measured in the BHE before the TRT have also commonly been used to evaluate the undisturbed ground temperature [10,[15][16][17].Numerical modeling can be used to reproduce an equilibrium temperature profile [18].The resulting profile can be used to extend the conventional in situ TC assessment of a TRT and to infer the terrestrial heat flow [18,19].Ground temperature profiles measured in observation wells are also used to evaluate other thermal properties such as in situ TD [20][21][22][23][24].A calculated ground temperature profile using such analytical or numerical approaches can be further matched to the observed ground temperature profile and used to evaluate vertical groundwater fluxes, hydraulic conductivities and thermal properties of the aquifer [22,25,26].The method is based on the evaluation of the annual amplitude temperature decay and the annual damping depth during long-term observation of the ground thermal disturbance diffusion resulting from the annual thermal flux at the ground surface [24,27,28].Often, one-dimensional semianalytical to analytical solutions or numerical simulations are used to infer the ground TD with various methods, such as the amplitude ratio, the phase lag and the harmonic method [10,[21][22][23][24][27][28][29][30][31][32].Despite the potential of these approaches to provide in situ evaluation of TD, they appear hardly applicable to the design of ground-coupled heat pump systems.The main reason is the time required to continuously measure the ground temperature profiles ranging from several days to a year to record a periodic cycle of heat diffusion in the subsurface [10].As a matter of fact, field measurements of prefeasibility studies for the design of ground-coupled heat pumps must be conducted within a few days, for instance, 2-3 days when considering the heating period of conventional TRTs or up to 5 days if the recovery phase of the TRT is included in the analysis [1][2][3][4]10,11,16,33].Recently, Márquez et al. [10] proposed a methodology for the indirect evaluation of the in situ TD.This approach assumes transient heat conduction in a semi-infinite medium and is based on the evaluation of the minimum and the maximum mean annual ground temperature measured in a shallow borehole and the depth where that mean annual ambient temperature is observed in the ground.In other words, the method proposed by Márquez et al. [10] is based on the assessment of the depth at which the annual ground temperature remains constant throughout the year (i.e., where there are no seasonal variations).That observed value is called the undisturbed ground temperature (TUGT) [34,35].The results obtained by Márquez et al. [10] were consistent with TD of mean geological materials identified from a reference borehole.However, this approach requires continuous ground temperature measurements during seven days to a year to identify and confirm the depth where TUGT is located.This field approach is time consuming for a ground-coupled heat pump project.Moreover, the result is sensitive to measurement errors since it is based on a single evaluation of TUGT [21,29,36].An entire temperature profile may provide more information from the BHE (in terms of geology and thermophysical properties of the geological materials) and thus potentially minimize uncertainty [10,15,16].
All the aforementioned difficulties in assessing accurately and quickly the in situ HC can eventually contribute to increasing errors in the design of geothermal heat pump systems and ultimately impact their installation cost [4,11].Therefore, evaluating the in situ HC can be useful for simulating the operation of BHEs used for both ground-coupled heat pump and underground thermal energy storage systems [11].In fact, the Giordano et al. [11] study indicated that the total drilling length of BHEs calculated when designing a ground-coupled heat pump system can be affected by ±6-7%, which influences the total system cost by 3-4%.This highlights that an accurate evaluation of the in situ HC can help to better design ground-coupled heat pump projects by accurately targeting their installation cost and, therefore, positively impacting the geothermal heating and cooling market [4,11,37,38].
Bearing in mind the importance of a quick and accurate assessment of in situ HC, this research study had the objective of developing an alternative heat tracing approach to evaluate in situ HC considering the main guidelines of the TRT [1,2,10,11].This study was carried out within the scope of a TRT performed in a pilot BHE [10,16].The resulting method relied on the measurement of a single equilibrium temperature profile that is not disturbed by the heat injection of a TRT or drilling of the BHE and can be recorded before the TRT.Analysis of this equilibrium temperature profile using heat tracing principles allowed the evaluation of in situ TD.Afterwards, HC was calculated based on the in situ TC evaluated with a TRT.This heat tracing method appears rather novel because the in situ HC is evaluated with a single ground temperature profile, which can be the same as evaluating TUGT before performing a conventional TRT.Moreover, using a single observed temperature profile determined at equilibrium and a TC assessment obtained from a conventional TRT analyzed with the slope method provide advantages in terms of simplicity.The methodology proposed does not need additional borehole or several temperature measurements.Finally, it does not depend on prior knowledge of the Earth heat flow since it is only based on an empirical approach to reproduce the observed undisturbed ground temperature profile measured in the BHE.

Theoretical Background
General concepts used in this heat tracing method are described below to provide the basis of the new field approach developed and applied in this study.The observed undisturbed ground temperature value (TUGT; °C) and the curve-fit between observed and calculated undisturbed ground temperature profiles are used as a criterion to constrain the analysis in order to infer the damping depth as well as the in situ TD and TC.
A practical approach to accurately evaluate TUGT is based on the graphical selection of the depth interval to be used for evaluation of the mean ground temperature and discarding near-surface data that visually appears to be affected by surface thermal disturbances [6,9,10,16,34,38].The depth at which TUGT is found or at which thermal disturbances from the surface are not perceived is called the depth of zero temperature amplitude (Figure 1) and it defines the boundary between thermocline and the thermostatic zones, where the geothermal gradient can be observed [34].Usually, estimation of the depth of zero temperature amplitude requires monitoring ground temperature profiles at different depths on a yearly basis [9,34], which we want to avoid here to fulfill the TRT practice [3,10,11,16].For this new method, the acquisition and analysis of ground temperature observations need to be short enough for the test to be reasonably implemented during prefeasibility studies of ground-coupled heat pumps.
Therefore, a new empirical equation is proposed to compute the damping depth (Equation ( 5)) from an observed ground temperature profile.It relies on empirical coefficients found by least square method with a solver applied to reproduce the ground temperature profile with calculations.This new empirical equation is described below and was obtained through analysis of the temperature profiles previously acquired during the field implementation of Pambou et al. [18].Furthermore, it was tested against the sites used by Kusuda et al. [20] and Xing [28].In this work, the application of the method is shown for TD and HC assessment in the scope of a TRT.
Thus, in this project, a ground temperature profile (Tobs (z)) is measured in the BHE before carrying out a TRT.Then, a heat transfer equation can be applied to the observed ground temperature profile, considering TD as an unknown parameter.With this assumption in mind, an empirical equation (Equation ( 1)) was developed to calculate a normalized ground temperature profile (Tg (z)) of the observed ground temperature (Tobs (z)), in which it is possible to evaluate a damping depth (Zdd) to infer TD by common analytical equation and the calculated undisturbed ground temperature (C1) at the depth of zero temperature amplitude.
A trial-and-error approach could be used to retrieve and approximate C1 and Zdd, but the analysis would be time consuming when considering the scope of a TRT.Therefore, a heat transfer equation (Equation ( 2)) related to the surface temperature variations and a solver optimization (Equation ( 3)) are suggested as the fastest method to accurately find C1 and Zdd values.These results are subsequently validated using Equation ( 4) and the curve-fit of the observed temperature against normalized depth-temperature profiles.
Detailed explanations of each equation used in this proposed heat tracing method are presented in the following subsection.

Calculated Undisturbed Ground Temperature Profile
A new empirical heat conduction equation (Equation ( 1)) was developed with the goal of calculating a normalized ground temperature profile at any depth z considering the entire length of the BHE.This proposed equation could be used to approximate the one-dimensional solution of the governing heat conduction equation considering a surface temperature diffusion in the half-infinite medium [9,21,23,[40][41][42].The proposed equation was described with an exponential form to normalize each observed ground temperature measurement that made an equilibrium (Tobs (z)) and it is defined as: where Tg (°C) is the normalized ground temperature at depth over the entire length of the BHE, C1 (°C) is the average ground temperature at the depth of zero temperature amplitude, Tobs (°C) and Tbot (°C) are the observed temperature at any depth and at the bottom of the borehole, respectively, Zdd (m) is the damping depth and Lobs (m) is the length of the borehole surveyed.Each ground temperature value calculated using Equation ( 1) is assumed to be a normalized value of the measured ground temperature at each depth over the entire length of the BHE.
Applying the heat balance concept, it is assumed that Zdd and C1 are integrated parameters of the thermocline zone (Figures 1 and 2) and the interpolated undisturbed ground temperature at the depth of zero temperature amplitude when considering onedimension heat conduction, respectively [19][20][21][22][23][24]27,28,[32][33][34].Therefore, an equation can be used as an upper boundary condition to describe the surface temperature variations transferred to the subsurface by oscillatory heat diffusion [28,33,[43][44][45][46].Thus, a new equation was defined to calculate a ground temperature, as sinusoidal function of the heat diffusion: where C1 (°C), C2 (°C), C3 (m) and C4 (m) are experimental coefficients that can be found by using a non-linear solver optimization (Equation ( 3)) related to the objective function (OF) and defined as: where OF (°C) is the sum of the squared residual computed from the difference between observed and calculated temperature at the same depth and 1 to N is the depth interval distribution covering the total length of the BHE.
The optimization function OF is validated when the bias error (BE) between C1 and the TUGT, inferred from the ground temperature profile measured in the BHE, is less than 5%, such that: where

Calculated Damping Depth
Using the experimental coefficients from Equation ( 2), Zdd is computed by a newly proposed field correlation defined as:

Calculated Subsurface Volumetric Heat Capacity
In situ HC is calculated directly from the analytical thermal diffusivity equation, with respect to the thermal conductivity inferred from a TRT done in the BHE ⦋4,11,47⦌: where HCcalc (J m −3 K −1 ) is the in situ heat capacity and λeff (W m −1 K −1 ) is the effective thermal conductivity inferred from the TRT.

Quality of Parameter Estimation
Statistical parameter analysis was used to evaluate accuracy and efficiency of the correlation between observed and calculated temperature.Relative error (RE) and root mean square error (RMSE) were calculated as: where yobs and ycalc are the observed and calculated parameters, respectively.N (−) is the total number of observations from the equilibrium temperature profile.RE (%) is an indicator of an overestimated (positive difference) versus an underestimated value (negative difference), while the RMSE (−) indicates the deviation magnitude in the range value.

Stepwise Procedure for Parameter Assessment
The following stepwise procedure (Figure 3) is suggested to summarize the parameter estimation analysis: 1. Accurately measure an equilibrium temperature profile (Tobs (z)) in a BHE and apply proper corrections for the rise of the water level in the U-pipe when using a wired probe as suggested by Pambou et al. [16]; 2. Perform a standard TRT and evaluate in situ TC; 3. Prepare the solver optimization to reproduce the normalized temperature profile (Tg (z)) through Equations ( 1) to (5); 4. Match the observed and calculated temperature profile using Equation ( 2) and refine the results by minimizing BE (Equation ( 4)), which describes the difference between TUGT and C1; 5. Evaluate the quality of parameter estimation using statistical analysis (Equations ( 8) and ( 9)) and proceed to the next step when the results are within the best value range and thus considered acceptable by statical analysis (Equations ( 8) and ( 9)); 6. Calculate the damping depth Zdd, in situ TD and HC (Equations ( 5)-( 7)).Microsoft Excel spreadsheet and its solver were used in this research to implement the equations, find the empirical coefficients by optimization and use the newly proposed method.The validation of the methodology is given in the next subsection.

Validation of the Proposed Method
The methodology proposed in this study was verified and validated by evaluating the in situ HC at the INRS geothermal experimental site in Quebec City (Figure 4).This geothermal experimental site was chosen for its scientific and technical interest due to the availability of BHEs and observation wells.In fact, several studies have been conducted at this site to improve the design of the BHEs and the methods for characterization of thermal properties [4,9,11,20,21,[40][41][42].For example, field measurements were performed with different equipment to accurately assess TUGT [4,10,16] and several types of TRTs were carried out in different seasons (fall, winter, summer) to infer TC and borehole thermal resistance [4,10,11,16,33] (see Figure 5, for an example of field setup).
Recently, other research of interest was made to evaluate heat flux density [19], subsurface thermostratigraphic log and groundwater flow [16], as well as in situ HC evaluation by analytical equations and numerical modeling [11].

Borehole Heat Exchanger and Site Description
The INRS geothermal experimental site has two BHEs (1U-pipe and 2U-pipe) and five observation wells (Obs) that were installed from 2015 [10] to 2020 [11].The single Upipe BHE used in this study has a diameter of 114 mm (4.5 in) and is grouted with a mix of bentonite and silica sand down to the entire depth of 154 m.The subsurface described at the site of INRS consists of shale bedrock under an overburden of 10 to 14 m in thickness (see Figure 6, for the stratigraphy at the location of the 1U-BHE).The shale bedrock is fractured and groundwater fluxes were inferred in the fractured zones [10,16].At the site, elevation of the water table varies from 14 to 16 m above sea level with northeast flow direction towards the Saint-Charles River.

Field Validation
Assessment of the mean TUGT is required to validate experimental coefficients C1 to C4 using the criterion described in the Equation ( 4).The first step involved estimating the mean TUGT from temperature profiles measured at equilibrium state.The ground temperature profiles were acquired in 2015 and 2016 before the TRTs, respectively [10,16;19].Two techniques were used for the temperature measurements.One was based on a submersible wired probe using a vertical spatial resolution of 1 m.The second was based on a fiber optic distributed temperature sensor with a spatial resolution of 0.25 m [16;19].The temperature measurements made with a wired probe were corrected for the rise of the water in the U-pipe of the BHE [16].In the second step, the coefficients C3 and C4 were used to calculate the damping depth Zdd and effective in situ TD and HC.The bulk TC previously estimated from a conventional TRT performed at this field site [10] was used for the evaluation of HC.
Additionally, the in situ HC estimated with the approach developed in this work was consecutively compared to the results obtained for the in situ HC determined by Giordano et al. [11].The approach developed by these authors to evaluate in situ HC is briefly described in the following lines to facilitate understanding and comparison.
Giordano et al. [11], at first, evaluated HC following the dual needle probe concept suggested in Raymond [4], using the 1-U BHE and the observation well (obs4; Figure 4).This allowed to independently assess the in situ HC and validate the OTRT method.The TRT was performed with a heating cable and temperature sensors in both the BHE and the observation well located 1.2 m apart.Temperature sensors were placed in the observation well at vertical distances varying from 2.5 to 5 m.The analysis was performed with the infinite line source equation and results from this test can be assumed as the most reliable field assessment of the in situ HC.As a second step, Giordano et al. [11] performed a sinusoidal heat injection for the OTRT based on Oberdorfer protocol [12].The oscillatory thermal response was analyzed with equations proposed by Eskilson [44].In situ TC was inferred using the infinite line source equation applied to the linear temperature component as in a conventional TRT [2,4,10,16].Then, in situ TD was calculated using the amplitude attenuation and the phase lag of the oscillatory component [20][21][22][23][24][27][28][29][30].HC was then evaluated similarly to what was done in this article using Equation ( 7).

Results
The validity and applicability of the new heat tracing method presented in this work to assess the in situ HC are presented and discussed below.A comparison between calculated and observed temperature was carried out, as well as between inferred thermal properties using results from various field methods applied at the same experimental site and previous geological characterization [10,11,16,18,19,45,46].

Observed Undisturbed Ground Temperature
Equilibrium ground temperature profiles (Tobs (z)) measured in warm and cold seasons were used to accurately analyze and evaluate the in situ TUTG (Figure 6).The evaluated TUTG was estimated to vary between 7.90 and 8.01 °C using temperature measurements from the depth interval 15 to 154 m and considering the temperature profiles measured at different times of the year with a vertical spatial resolution of 1 m (submersible sensor [10,16]) and 0.25 m (fiber optic [19]).
These profiles highlight the influence of the seasonal air temperature variations and the heat diffusion within the subsurface at the INRS experimental site.Two zones can be defined from the temperature profiles with analogy to Figures 1 and 2. These two zones are: the upper part which is influenced by the surface conditions (thermocline zone) and the lower part which is not affected by seasonal variations (thermostatic zone) but the geothermal gradient, on the order of 12.00 °C km −1 (Figure 7; [19,34,46]).The temperature profiles acquired show an inverted gradient in the upper thermostatic zone (Figure 7).This inverted gradient can be due to recent climate warming [19,34,45,47].Finally, it can be observed that in the thermostatic zone, both the temperature profile from the observation well and from the BHE have the same behavior.These results are in a good agreement with the assumption of temperature diffusion by heat conduction in a homogeneous and isotropic media [37,40,48,49].

Assessment of Empirical Coefficients
The TUTG was evaluated using a single equilibrium temperature profile (Tobs (z)) measured in 2015 before a conventional TRT was done (Figures 4 and 7).The calculated TUTG, which is assumed equal to C1, gave the value 7.96 °C.The absolute difference between measured mean ground temperature (7.90-8.01°C) and the calculated TUTG (or C1) ranges between 0.69% and 0.57%.Such results indicate a low bias error (Equation ( 4)).Furthermore, these results are used to validate the empirical coefficients C1 to C4 found by the solver (Equations ( 3) and (4); Table 1).
Table 1.Experimental coefficients estimated from solver optimization (Equations ( 2)-( 4)) applied on the temperature profile measured (Tobs (z)) in 2015 at the INRS experimental site before a TRT.

Empirical Coefficient Value
The results obtained suggest that calculated C1 is similar to TUTG evaluated from the measured temperature profile and is in the range of the validation criterion BE < 5%.Moreover, the value C2 equals to zero suggests that the calculated TUTG (or C1) is assessed close to the depth of zero temperature amplitude (Figures 1 and 2).Thus, these experimental coefficients C3 and C4 can be used for the next step to ultimately evaluate the in situ HC (Figure 3).The curve-fit between normalized ground temperature profile (Tg (z)) based on Equation (1) against the measured temperature profile (Tobs (z)) at equilibrium state during the TRT [16,39] is plotted in Figure 8.

HC Estimation Using the New Empirical Method
The conventional TRT done on this BHE revealed a bulk TC of about 1.75 W m −1 K −1 [10].Thus, the resulting subsurface HC (Equation ( 7)) is 2.40 MJ m −3 K −1 (Table 2) when considering the above TD.Temperature measurements taken at a depth of about 22 m in the observation well obs4 (Figure 4) with the heating cable TRT were used to evaluate both in situ TD and HC [11].The results, which are thought most accurate among Giordano field method [11], are presented for each submersible temperature sensor inserted at different depths in the BHE and ground layer encountered in the observation well obs4 (Table 3).Field results from dual needle concept were compared with those from this new heat tracing method using Equations ( 6) and (7) to infer the in situ HC (Tables 2 and 3).The average absolute discrepancy for the upper layer made of mixed unconsolidated sedimentary deposits and weathered shale (0 to 12.50 m) was 13.01%, while that for the lower layer made of shale (17 m to 22 m) was 4.48% (Table 3).The average absolute discrepancy considering both layers was 10.16% (Table 3; Figure 5).This difference between field results obtained at the same site using different methods was considered small enough, and therefore, this new heat tracing method was confirmed reliable (Table 3).

Comparison of Calculated Subsurface Heat Capacity with OTRT Method
Results from this new heating tracing method were also compared with results from the OTRT method of Giordano et al. [11] (Table 4).Absolute discrepancy when considering the HC evaluated with the oscillatory resistance method and phase shift method was found to be greater than 26% compared with this new empirical approach (Table 2 versus Table 4).Table 4. Subsurface HC inferred with the new empirical method compared to the results obtained by [11] using the OTRT method at the INRS experimental site.The results obtained from the recovery period revealed an underestimated value on the order of 10%; while the difference with the oscillatory resistance analysis and from the phase shift analysis are out of 15%, respectively (Table 2 versus Table 4).The differences found in these values suggest that the OTRT may be affected by uncertainties, and it may be useful to use corrected factors to adjust the range of the values when evaluating in situ HC using the OTRT (Table 4).The sources for such variability can be caused by, for instance, the HC of the grout filling the BHE [11].

Analysis
Despite some sources of uncertainty that need to be addressed, such as the influence of the backfilling material, the project made by Giordano et al. [11] proposed alternative avenues to improve in situ assessment of subsurface HC.Analysis of results reveals that the concept used by Giordano et al. [11] is straightforward but the implementation in the field of a sinusoidal heat injection can be complex and need specific analytical expertise, which is not always available.Finally, these results suggest that the new heat tracing method here developed and discussed below is, at present, more reliable than the OTRT and simpler than conducting a TRT with the dual probe concept requiring an observation well (Tables 2-4).

Discussion and Conclusions
A novel one-dimensional heat tracing field method was developed and applied to accurately evaluate the undisturbed ground temperature (TUGT = C1) and the damping depth (Zdd) of the surface temperature changes in the subsurface at the vicinity of the BHE.These values were used to infer the in situ TD (αcalc) and, subsequently, the in situ HC.The results were compared with other field methods undertaken at the same site.This comparison suggests that using this new empirical method in combination with the conventional TRT to assess TC can reduce uncertainty when characterizing heat storage properties.In other words, TUGT and effective subsurface TD and HC can now be inferred with the new approach in the scope of a TRT over the geological materials intercepted by a BHE.
This empirical heat tracing method assumes transient heat conduction mechanism of the surface temperature variations diffusing through an isotropic and homogeneous semiinfinite medium (i.e., the subsurface in which an equilibrium temperature profile is measured).This newly proposed method uses least squares and a nonlinear solver optimization to fit the observed temperature profile and to find the experimental coefficients (C1 to C4).These parameters are incorporated in the new semi-analytical ground temperature sinusoidal function which is assumed to be an upper boundary condition of the heat conduction equation [46,50,51].This method thus relies on an accurate equilibrium ground temperature profile measured in a BHE before a TRT [10,17].The calculated experimental coefficients are used to evaluate a damping depth and an effective TD.The bulk HC of the geological materials is afterwards estimated using a bulk in situ TC inferred from a conventional TRT.Evaluation of TC was made in the article with the infinite line source equation but could be made with a numerical approach as suggested by other authors [1,4,18,30,32].This can help to better evaluate TC but will not significantly change the assessment of TD and HC, mostly relying on the determination of the empirical coefficients from an undisturbed temperature profile, in this newly proposed heat tracing approach.
The in situ HC evaluated with the newly proposed method at the INRS experimental site in Quebec City was successfully compared to that inferred by other field methods (Tables 3 and 4) as a criterion of validation [11,51].Furthermore, these results were in the range of thermal properties for geological materials of the Quebec City area [46].Hence, the obtained results validate the model assumptions and the parameter estimation procedure.Consequently, the stepwise implementation of this new method (Figure 3) can be conducted in the scope of ground-coupled heat pump system design [1,3,4,10,11,19].Moreover, this method proposes a new damping depth equation which does not rely on the temperature amplitude as the previous methods adopted to evaluate the in situ TD [20][21][22][23][24][27][28][29][30].In addition, the newly proposed heat tracing method does not rely either on prior knowledge of the subsurface heat flux [18,19,51] or does not require time series of the annual temperature monitoring [9,14,[20][21][22][23][24][25][26][27][28][29][30].This is an advantage compared to previous damping depth methods used for the same purpose [34,35].Field measurements only require a single temperature profile that can be rapidly collected before a conventional TRT.This highlights the advantages and novelty of the proposed methodology when compared to the sampling steps required by other approaches using time series of ground temperature measurements [20][21][22][23][24][27][28][29][30]45,52,53].Another field benefit is that the new method does not need an additional borehole when compared to the dual probe method experimented in the field at the scale of a BHE by Giordano et al. [11].A practical advantage is also related to the mathematical formulation that can be easily implemented and optimized with a built-in solver found in a spreadsheet program.Considering one hour for the field setup and measurements with a wireline temperature probe for a BHE of 154 m depth and a single analysis that should not require more than half an hour for data processing, the method can be qualified as fast when compared against commonly used field and laboratory methods [4,11,18].
Moreover, results obtained using this newly proposed heat tracing method depend on experimental coefficients that, in turn, rely on an accurate ground temperature profile measurement that was successfully reproduced at the INRS experimental site.Current practice is to infer the mean ground temperature from the measurements of a ground temperature profile that is acquired before a TRT and by lowering a submersible temperature datalogger in the BHE [10,16,18].In some cases, measurements can be done with a 1 m spatial resolution over the length of the BHE that can reach 160 to 200 m.Measurements using a spatial resolution of 5 to 10 m are typically not good enough for this proposed method.Moreover, care should be taken with the field procedure by selecting an accurate temperature probe and by correcting the temperature profile for the water level rise when measured in a U-pipe [16].For example, the submersible temperature sensor and pressure probe used in this study had a ±2×10 −3 °C accuracy, <5×10 −5 °C resolution, ±5×10 −4 dbar accuracy and <1×10 −5 dbar resolution.
It is worth highlighting that these results are obtained with the assumption of conductive heat transfer in a homogeneous and isotropic subsurface.However, a validation process through field-testing complemented by numerical modeling is recommended in case of a strongly heterogeneous and anisotropic subsurface or the presence of significant groundwater flow [3,14,22,48,53,54].
As scientific contributions, this study puts forward a new damping depth equation and field estimation of TD and HC relying on a single measured equilibrium ground temperature profile.It confirms that ground temperature profiles measured in BHEs are an inexpensive source of data that can be analyzed to obtain more information on the subsurface thermophysical properties [4,9,11,[14][15][16][17][18][19][20][21][22][23].These contributions provide advantages for the design of ground-coupled heat pump systems by considering this heat tracing method as a complementary in situ tool for improving conventional TRTs [1,2,4].

Figure 1 .
Figure 1.Theoretical ground temperature distribution showing the depth of zero temperature amplitude.The blue triangle indicates the geothermal gradient.

Figure 2 .
Figure 2. Ground temperature distribution showing the damping depth and the depth of zero temperature amplitude, with (a) annual ground temperature profile span and (b) components of the daily oscillatory ground temperature influenced by surface temperature variations diffusing in the shallow subsurface.

Figure 3 .
Figure 3. Procedure for assessing the in situ HC from a temperature profile undisturbed by a TRT and the in situ TC inferred from the TRT.

Figure 4 .
Figure 4. Geothermal experimental site at INRS (Quebec City).Numbers beside boreholes (obs and U) and along yellow lines indicate elevation of the water table in meters above the sea level and local potentiometric level, respectively See field set up in 11Obs are four open observation wells, while 1-Uand 2-U are the borehole heat exchangers.

Figure 5 .
Figure 5. (a) Field equipment for ground temperature profile measurement at INRS geothermal experimental site (Quebec), with examples of (b) how to conduct measurements in the single U-pipe BHE and (c) results of observed ground temperature profiling; (2) and (3) are uncorrected for water level rise in the 1-U pipe, while (1) is corrected, respectively (see field protocol in Pambou et al.[16]).

Figure 6 .
Figure 6.Thermostratigraphic log of the 1UBHE at the INRS experimental site, with layered TC from Pambou et al. [19] and geological cross-correlation based on borehole logs from Raymond et al. [10].

Figure 7 .
Figure 7. Observed undisturbed ground temperatures at the INRS experimental site.

Table 2 .
In situ TD and HC estimated (Equations (2)-(6)) by applying the new heat tracing method on the temperature profile measured in 2015 at the INRS experimental site before a TRT.

Table 3 .
[11]urface HC inferred with the new empirical method compared to the results obtained by[11]using the dual needle concept at the INRS experimental site.Mean values of each method and mean relative difference in bold.