The Earth's temperature field is one of the controlling factors of petrophysical properties, such as seismic velocity and rheological strength. Borehole geothermal measurements open a window into the thermal regime of deep Earth, allowing extrapolations to heat flow and geotherms, with the combination of thermal conductivity and heat production measurements (Furlong & Chapman, 2013). Furthermore, the temperatures within the shallow boreholes are often perturbed by groundwater movement (Guillou-Frottier et al., 2013;Jessop, 1993) and changes in ground surface temperature (GST) induced by paleoclimate variations and anthropogenic land use (Gosnold et al., 2011;Majorowicz et al., 2012). Paleoclimatic signatures could perturb equilibrium heat flow densities to more than 2 km (Beltrami et al., 2014). Therefore, detailed borehole geothermal studies have the potential to provide unique constraints on both deep-earth thermal processes and GST changes.Available global heat-flow and paleoclimate studies primarily rely on shallow (<1,000 m) boreholes (Cuesta-Valero et al., 2021;Davies, 2013;Huang et al., 2000). Long-term climate changes cause significant perturbations in the measured heat flow, leading to deviations in the global heat-flow data set. Climatic perturbations, which are considered noise to determine the heat flow, can be used for GST reconstruction (Cermak and Bodri, 2011).