Borehole thermometry is an important tool for reconstructing past climate conditions, assessing changes in land energy storage, and understanding subsurface thermal regimes such as permafrost and glacial dynamics. Optimizing the temperature sensor placement within boreholes allows us to maximize the informativeness of temperature measurements, particularly in polar regions where operational constraints necessitate cost-effective solutions. Traditional sensor placement methods such as linear or exponential spacing, often overlook site-specific subsurface ice heat distribution characteristics, potentially limiting the accuracy of the measured temperature profile. In this paper, we propose a greedy optimal sampling technique for strategically placing temperature sensors in ice boreholes. Utilizing heat transfer model simulations, this method selects sensor locations that minimize interpolation errors in reconstructed temperature profiles. We apply our approach to two distinct ice borehole sites: EPICA Dronning Maud Land site in East Antarctica and the Greenland Ice Core Project site, each with unique surface conditions. Our results demonstrate that the greedy optimal sensor placement significantly outperforms conventional linear and exponential spacing methods, reducing sampling errors by up to a factor of ten and thus achieving similar informativeness with fewer sensors. This strategy offers a cost–effective means to maximize the information obtained from ice borehole temperature measurements, thereby potentially enhancing the precision of climate reconstructions.
Horizon 2020716092Introduction
The Earth's climate system is experiencing an imbalance, with continuous accumulation of heat over the past decades, leading to warming across various components including the ocean, land, cryosphere, and atmosphere . High precision measurements of subsurface ice and soil temperatures are essential for assessing changes in land energy storage, refining ice models , and determining the thermal state of permafrost . Subsurface temperature profiles obtained through borehole thermometry provide valuable information about the surface temperature evolution, glacial thermal regimes, and geothermal heat flux . Borehole temperature logs used for climate reconstruction must balance the precision of temperature measurements with the depth resolution . In addition, temporal resolution plays an important role: while many borehole temperature measurements are performed at a single point in time by continuously or stepwise lowering a single sensor, continuous measurements over time have several advantages. These continuous measurements help average out seasonal effects or fast temporal fluctuations, utilize transient data to better constrain the inverse temperature reconstruction problem , and allow for measurements in boreholes that are still disturbed from the drilling process. For example, placed a temperature chain with 16 temperature sensors in 80–90 m deep boreholes. Continuous measurements were then performed for over a year and sent via satellite with the objective of reconstructing multi-decadal temperature trends in East Antarctica. This strategy allowed the team to discard the initial months when the borehole was still affected by the drilling process and to average out the seasonal cycle, all without requiring a second expedition to return to the borehole. Thus, high precision temperature measurement technology for continuously monitoring borehole temperatures is an important tool for climate reconstructions. The system's ability to withstand extremely cold temperatures, its ultra low power consumption, and its high precision and long-term stability are crucial . To minimize costs and allow this temperature sensing and data transmission system to operate over multiple years on limited energy supplies, it may be necessary to compromise on the number of sensing devices. Therefore, one needs to be strategic about where to place the limited number of sensors in the borehole to maximize the informativeness of the measurements with minimal effort.
The distribution of heat in boreholes is mostly determined by diffusion and advection of surface temperature anomalies as well as the geothermal heat flux . Temperature logs from terrestrial boreholes used to reconstruct surface temperature histories vary widely in the number of measurements, ranging from 10 to over a hundred for 200–600 m deep borehole profiles . For surface temperature reconstruction, and used an ad hoc spacing of temperature sensors in the borehole, with the distance between the sensors increasing down the borehole. suggested that by exponential positioning of sensors, each data point contributes an identical amount of information to a reconstructed surface temperature history. Beyond sensor placement, uncertainties in borehole-based reconstructions can also arise from factors such as insufficient borehole depth, which hinders the separation of climatic and geothermal signals , and heterogeneous thermal properties of soils and bedrock, which vary with depth and location . These challenges are typically less significant in ice boreholes, where the thermal properties of snow and ice are homogeneous and well known and the boundary conditions are generally better constrained. Although the topic holds significance and potential, sensor placement in boreholes has not received much attention in the field of ice borehole thermometry.
At the same time, the sensor placement problem has been extensively studied in other research areas such as weather forecasting and environmental monitoring. A study conducted to control the environment inside the greenhouse, used error based and entropy based methods for optimal sensor placement to track temperature variations . Active data selection and test point rejection based on error variance , Gaussian Process based model to reduce uncertainty , and convolutional Gaussian neural process are few of the machine learning related approaches for sensor placement. Finding effective sensor placements in the context of borehole thermometry can be accomplished by applying similar ideas of selecting sensor positions to minimize error.
The knowledge of the heat distribution in ice is necessary for determining the best location of sensors for ice borehole thermometry. Snow falls on top of the ice sheet and is then densified into firn, eventually becoming ice at a depth of approximately 80–100 m. Over time, heat gradually diffuses into the ice, with stronger fluctuations in temperature occurring closer to the surface and smoothing out deeper within the ice sheet. This phenomenon primarily motivates the ad hoc placement of borehole sensors discussed above, where the distance between measurement points increases further down the borehole .
In this paper, we propose an approach to strategically install a set of sensors to measure the temperature of an ice borehole to gather information on borehole temperature profile as effectively as possible. Using heat transfer model simulations, a greedy sampling technique is employed to optimize the sensor placement. Greedy algorithms solve problems by selecting the solution that is locally optimal at each phase . We start by formularizing a greedy optimal sampling algorithm for finding sensor placements in the borehole. We evaluate the approach of greedy optimal sensor placement for two distinct boreholes with distinct local surface conditions and examine our findings by contrasting it with the linear and exponential sensor placements. To demonstrate how to find an optimal and cost-effective arrangement of sensors, we analyze how the informativeness of borehole thermometry depends on the number of sensors used. We further discuss the uncertainty in the algorithm and data used in this study.
Method and data
We discuss the underlying subsurface ice heat transfer model, which is given as a one-dimensional diffusion advection equation. The sensor placement problem is solved by using a greedy optimal sampling algorithm that utilizes the heat transfer model applied to a set of possible past surface temperature evolutions to include the physics of the problem into the sensor optimization. We also describe the physical parameters of the two example sites considered in this study.
Heat transfer model
The heat transfer of the surface temperature θ(t) into the ice sheet is modeled by the one-dimensional heat diffusion advection equation
∂Tθ∂t=k∂2Tθ∂z2-w∂Tθ∂z,
where Tθ is the temperature as a function of time t and depth z (positive downwards). As a convention, the model is defined for time t∈[t0,tpr], hence heat transfer is considered from t0 years in the past to the present time (tpr). w is the vertical velocity and k is the thermal diffusivity. The thermal diffusivity is calculated as
k=Kρc,
where c is the specific heat capacity, K is the thermal conductivity, and ρ is the density of the firn/ice. Details of the parametrization of density, thermal properties and vertical velocity are given in Appendix . The (vertical) spatial domain is [0,H], where H is the thickness of the ice sheet and z=0 corresponds to the surface, see Fig. . Note that we consider [0,L] (L<H) as the subdomain of the borehole in ice, while still modeling the full vertical temperature profile on domain [0,H], which includes the borehole and the below remaining ice sheet.
To solve Eq. (), we use the forward Euler method in time, central differences for the diffusion term, and forward differences for the advection term. The spatial grid (in z direction) with uniform step size Δz is ZF⊂[0,H]. We employ Dirichlet boundary conditions for the simulations in this study. The boundary conditions are given as
3aTθ(t,0)=θ(t),3bTθ(t,H)=θb,
representing the top and bottom boundary conditions of the model. The top boundary condition models the input of the surface temperature into the ice sheet. The bottom boundary condition is set to the basal temperature of the ice sheet θb. In order to determine the initial condition (Tθ(t0,z)), the model is run for 50 000 years by setting the top boundary condition to a constant mean temperature Tθ(t,0)=θm, until the profile reaches equilibrium. The result of the temperature profile simulation for surface temperature θ(t) is Tθ(tpr,z) which can be denoted as TF. This represents the vector of approximated temperatures at all grid points ZF of the discretized heat transfer model.
Sensor placement problem
Our goal is to select the set of n best sensor locations ZS in the borehole of depth L (Fig. ), so that the sampling error in the reconstruction of the borehole temperature profile is minimized. The sensor locations ZS are selected from a set of candidate sensor locations ZL. Here, this set is selected as a uniform grid of size l, over the interval [0,L], where l is the number of possible candidate sensor locations.
Schematic diagram representing space vectors used in the algorithm. ZS is the set of selected sensor locations, ZL is the grid of candidate sensor locations, and ZF is the spatial grid of the model.
Sampling error calculation
Our aim is to develop a measure of error or uncertainty, for a given choice of sensor placements ZS. To this end, we evaluate simulated temperature profiles for different surface temperature evolutions at locations ZS to find the respective set of temperature values TS. Then, we estimate the (interpolation) error that we introduce if we aim to recover the full temperature profiles from this artificial data (ZS and TS for each profile). We use cubic spline interpolation for this study (Fig. in Appendix ). Sampling error is calculated as the root mean square over these interpolation errors obtained for each of the given temperature profiles of possible surface temperature histories, thereby incorporating a prior assumption on possible past surface temperature evolutions. This procedure is formulated in Algorithm .
Note that in our algorithm, the final sampling error (ϵs) is influenced by the error due to sensor placement (i.e., interpolation error) and by the error due to uncertainties in the sensor device, which is called device error (ϵd). Here, ϵd is drawn from a normal distribution N(0,σd2) and is considered independent between borehole simulations and across sensors. The artificial temperature measurements TS for sensors ZS are generated by adding an error term ϵd.
Sampling_error.
Input: Sensor locations ZS, standard deviation of sensor device error σd, a set of P temperature profiles {TFp}p=1P, candidate sensor locations ZL.
Output: Sampling error ϵs
begin:
p←1
whilep≤Pdo
Create an interpolating function T̃F,p for the pth simulated temperature profile, T̃F,p(z)←I(ZF,TFp).
Draw device error ϵd from normal distribution with given standard deviation N(0,σd2I).
Create artificial measurements TSp for sensors in ZS using the function T̃F,p and device error, TSp←T̃F,p(ZS)+ϵd.
Create interpolation function T̃S,p for sensor locations and artificial sensor measurements, T̃S,p(z)←I(ZS,TSp).
Calculate error grid as the difference between two functions evaluated at all candidate locations, ϵsp=[|T̃S,p(z)-T̃F,p(z)|]z∈ZL.
p←p+1
endwhile
Calculate final sampling error as the root mean square across error grids with respect to all P profiles, ϵs=∑p=1P(ϵsp)2P.
end
Algorithm thus requires the given choice of sensor locations, standard deviation of sensor device error, a set of borehole temperature profiles, and all possible (candidate) sensor locations for the calculation of sampling error. For each of the borehole temperature profiles, the algorithm creates an interpolation function. Artificial sensor measurements are created by evaluating given sensor locations using this interpolation function and adding device error to them. A separate interpolation function is then created for this artificial measurements. Both these functions are evaluated at all possible candidate locations, and the difference between them is recorded as the error grid for each borehole temperature profile. The final sampling error is calculated as the root mean square of the error grids of all borehole temperature profiles.
Sensor placement using greedy optimal sampling
Greedy optimal sampling is an adaptive method of adding sensors to lower the sampling error . The sampling error is computed after each sensor is added, and the subsequent sensor is positioned at the location on the grid of candidate sensor locations ZL, which has the maximum error.
This procedure is summarized in Algorithm . The algorithm begins by initializing ZS with a set of sensor positions ZS*, in which one sensor is fixed to the top and one to the bottom position of the borehole, while the others are randomly selected using a quasi-random sequence with upper bound as L and lower bound as 0. The sampling error ϵs which represents the accumulated error over all candidate sensor location ZL (refer Sect. ) is calculated with respect to the current set of sensor positions in ZS. The algorithm identifies the next sensor position zopt as the location in ZL which has the maximum error. The optimal sensor position zopt thus found from the set of candidate sensor locations is added to the existing set of selected sensor positions in ZS, and this process is repeated with updated ZS each time, until all the remaining sensor positions are found.
There is a possibility of bias in the sensor locations ZS found using Algorithm resulting from ZS* being fixed at the beginning. This is reduced by averaging the various sets of sensor locations ZS obtained with respect to different ZS* in the greedy optimal sampling algorithm (Algorithm ). We average these different sets of ZS giving averaged sensor placements ZS‾. It is assumed that this averaging procedure allows to reduce the bias due to ZS* being fixed (also see Sect. ). ZS‾ is referred to the greedy optimal sensor placement in this study, and unless otherwise specified, ZS‾ is obtained by averaging over 1000 sets of ZS.
Greedy optimal sampling.
Input: No. of available sensors n, initial set of sensor locations ZS*, standard deviation of sensor device error σd, a set of P temperature profiles {TFp}p=1P, candidate sensor locations ZL.
Output: Set of sensor locations ZS
begin:
Initialize the set of sensor locations, ZS←ZS*.
Set n* with the size of initial set of sensor locations ZS*.
i←n*
whilei<(n-n*)do
Find the sampling error ϵs for the current set of sensor locations, ϵs←sampling_error(ZS,σd,ZF,{TFp}p=1P,ZL).
Select the location from ZL that corresponds to maximum error in ϵs as the next sensor location in-order to minimize sampling error,
j^←argmaxjϵs(j),
zopt←ZL(j^).
Add zopt to the current set of sensor locations, ZS←ZS∪zopt.
Sort the set ZS.
i←i+1
endwhile
end
Other sensor placement schemes
Linear and exponential spacing of sensor locations are used for comparative analysis in this paper. In what we denominate as linear sensor placement, sensors positions zi∈ZS are assumed to be equally spaced along the sample space, hence
zi=(i-1)d1≤i<nLi=n
where, d=Ln-1. In exponential sensor placement, sensors positions zi∈ZS are assumed to be spaced evenly on a log scale along the sample space, hence
zi=(r)i-1-11≤i<nLi=n
where, r is the common ratio between adjacent sensor positions.
DataParameters of the study sites
We study two known ice coring sites chosen to represent distinct boundary conditions (Table ): EPICA Dronning Maud Land drilling site at Kohnen Station, East Antarctica (EDML) (75.00° S, 0.07° E) , and the site of the Greenland Ice Core Project (GRIP) (72.35° N, 38.30° W). GRIP has an accumulation rate of 230 mm yr−1 and a mean surface snow temperature of -31.7 °C (borehole temperature at 10 m) , while EDML has a lower accumulation rate of 64 mm yr−1 and a mean surface snow temperature of -44.1 °C (borehole temperature at 10 m, measured in 2024).
Accumulation rates: GRIP and EDML , Snow surface temperature: GRIP and own measurements 2024, Basal temperature and ice sheet thickness: GRIP and EDML , Surface pressure: , Snow surface density: GRIP and EDML .
The heat transfer model is discretized with a temporal grid of size 15 984 (Δt=2-4 years) and the spatial grid of size 695 for EDML and 757 for GRIP leading to a spatial resolution of Δz∼4 m for both sites.
Surface temperature time series
We generate an ensemble of surrogate local surface temperature time series designed to capture the full range of possible temperature histories. Each surrogate is a random time-series (η(t,β)) following a power law spectral density with β as its scaling exponent and the variance of the local observed surface temperature, as well as a global warming component based on anomalies in global mean temperatures (GMT) multiplied by a factor PA, representing the potential amount of polar amplification at the borehole site.
θ(t)=η(t,β)+PA×GMT(t)
In particular, this approach considers the temporal covariance structure of the stochastic component, which mimics natural climate variability (such as weather), and accounts for the uncertain amplitude of the deterministic warming trend. A simple parametrization to describe observed climate variability across a wide range of timescales is to assume a power–law relationship for the power spectral density (PSD) of the signal, S(f). This relationship is expressed as
S(f)∝f-β,
where β represents the scaling exponent. In state-of-the-art climate models, a scaling exponent of β=0.2 is typically found for Antarctica . However, data from the EDML ice core and nearby ice cores for the past 1000 years suggest a scaling exponent of approximately β=0.6 after adjusting for local non-climate variability . In contrast, marine and terrestrial temperature proxy data indicates a scaling exponent closer to β=1. In this study, we select values of β between 0 and 1 to encompass the range of plausible scaling behaviors for surface temperature.
We generate 1000 year long surrogate time series as the sum of a realization of a stochastic process mimicking natural variability and the global mean surface temperature time series to represent the global warming component. The random time series follow the prescribed power-law spectral density and have the variance of the observed 2 m air temperature time-series from the ice core locations from NOAA 20th century reanalysis spanning 1836–2015. The global warming component are the anomalies in the NOAA 20th century reanalysis global mean temperature (GMT) with respect to mean value in 1836–1950 and multiplied by factor PA=1,2,3 to represent the uncertain polar amplification at the borehole site. This component is set to zero before 1836. A scenario without any warming component is referred to as PA=0.
The simulations of borehole temperature profiles for each site were carried out using their respective site specific parameters (Table ). For each site, 90 borehole temperature profiles are simulated with respect to different aforementioned surrogate local surface temperature time series θ(t). Ten samples for each of the nine therein discussed, different types of θ(t), obtained by various combinations of warming factor (PA=1,2,3) and stochastic component (β=0,0.6,1), are used in the borehole simulations. They serve as prior (knowledge) that we apply, hence narrowing down the optimization process to realistic past surface temperature histories.
Sensor sampling parameters
The settings of input parameters for sampling error calculation (Algorithm ) and greedy optimal sampling (Algorithm ) are based on few basic values including length of the borehole, length of the error grid, total and number of sensors to be placed in borehole. In this study, these values (Table ) are taken to be the same for both sites.
Values used for generating input parameters of Algorithms and .
ParameterSymbolValueBorehole length (in m)L200Length of the error gridl105No. of sensorsn20No. of initial sensorsn*4
Note: n* is used only in Algorithm .
ResultsComparative analysis of linear, exponential and greedy optimal sensor placements for borehole thermometry
We first compare the performance of linear, exponential, and greedy optimal sensor placements for n=20 sensors for the EDML site, assuming no device error (Fig. ). Naturally, the sampling error is near zero at the depths where sensors are placed and maximal between sensors. For the linear placement, the maximum sampling error is larger than 100 mK, which occurs between the top two sensors (Fig. b, c). For the exponential spacing, the error exceeds 2 mK, which occurs between the deepest two sensors (Fig. e, f). The greedy optimal spacing lies between the linear and exponential cases, with larger spacing deeper in the borehole but less than the exponential increase. For the greedy optimal sensor placement, the maximum sampling error is in the top 50 m but is less than 0.2 mK (Fig. g–i), and thus an improvement of a factor over 500 compared to the linear spacing.
Strong fluctuations in heat distribution occur closer to the surface and decrease deeper into the ice sheet. Therefore, many studies considered exponential and ad hoc sensor placements with sensors being placed farther apart down the borehole to recover heat distribution effectively. The greedy optimal sensor placements (Fig. g) manifests more or less naturally this trend of sensors being placed farther apart down the borehole. The last sensor is an exception due to its fixed location at the very bottom of the borehole.
Performance comparison of linear, exponential, and greedy optimal sensor placements in case of EDML, assuming no device error. Panels (a)–(c) are based on linear spacing, (d)–(f) are based on exponential spacing and, (g)–(i) are based on greedy optimal sampling of a set of 20 sensors in a 200 m borehole. The horizontal axis of (a), (d) and (g) shows the index of the sensors, sensor with index-1 is placed at the top and index-20 is placed at the bottom of the borehole. Panels (b), (e) and (h) show the sampling error due to linear, exponential, and greedy optimal sensor placement respectively. Panels (c), (f) and (i) show sampling error in logarithmic scale. The dashed line marks 1 mK in (c), (f), and (i).
Comparison of the results of greedy optimal sensor placements for the cases of EDML and GRIP, assuming no device error. (a) Depicts the greedy optimal sensor placements for both sites, for a set of 20 sensors in the 200 m boreholes and (b) depicts their respective sampling error.
Interestingly, generally similar results are also obtained for the GRIP ice core site that is characterized by very different site conditions. Comparing the results for both sites in detail, Fig. shows that due to higher accumulation rate of GRIP compared to EDML, the signals are drifted slightly deeper, and therefore the GRIP sensors are placed slightly deeper in the borehole than the corresponding sensors for EDML. The difference between the corresponding sensor locations of the two sites is noticeable (>0.5m) for sensors placed between ∼ 10 and 130 m depths (Fig. a). The sampling error curve appears to be slightly shifted from one another, and the maximum sampling errors of EDML and GRIP are ∼ 0.2 mK (Fig. b).
Dependence of sampling error on the number of sensors used for borehole thermometry
To evaluate the dependence of sampling error on the number of sensors used in the borehole sensor placement, we generate linear, exponential, and greedy optimal sensor placements for n=5 to 20 sensors, assuming initially no device error. For the EDML site, as expected, the maximum sampling error decreases with the increasing number of sensors for all three types of sensor placement (Fig. a). Similar results are also observed in the case of GRIP (Fig. a in Appendix ).
The sampling of linear, exponential and greedy optimal sensor placements with and without device error were computed for 5 to 20 sensors for EDML. Maximal sampling errors without device error (σd=0) (a), with 5 mK device error (σd=5 mK) (b) and with 10 mK device error (σd=10 mK) (c) are given. The number of sensors used in sensor placement is denoted on the horizontal axis of (a), (b) and (c).
We further investigate the influence of a device error on the sampling error. We again generate linear, exponential, and greedy optimal sensor placements for n=5 to 20 sensors, assuming device error of standard deviation σd=5 and 10 mK (Fig. b, c). The trend of maximum sampling error decreasing monotonically with increasing number of sensors is observed for linear and exponential sensor placement, but not reduced to the device error (σd line) with 20 sensors. In the case of greedy optimal sensor placement, this trend is observed until the error is in the range of the device error and then the error continues to be in that range with further increase in number of sensors (Fig. b, c). Considering that the sensors have device error, it is clear that the device error of the sensors influences the sampling error. The greedy optimal sensor placements outperformed linear and exponential spacing of sensors as it reduces the maximum sampling error to the range of device error with fewer sensors.
Assuming the sensors have device error, it is important to understand if we can reduce the number of sensors and get good results. As mentioned before, both linear and exponential sensor placement was not able to reduce the sampling error to the range of device error with 20 sensors. The sampling error is not reduced below 100 mK with linear sensor placements considering device error of standard deviation σd=5 and 10 mK. For exponential sensor placement the sampling error is not reduced below 8 and 16 mK respectively, when device error of standard deviation σd=5 and 10 mK is considered. It is evident that only half the number of sensors is required by greedy optimal sensor placement to reduce the sampling error to the close vicinity of the given σd, when compared to respective linear and exponential placements (Fig. b, c). GRIP also showed similar results when analyzing the influences of device errors (Fig. b, c).
To assess the reliability of greedy optimal sensor placement with fewer sensors, we tested one-by-one sensor failures for sets ranging from 11 to 20 sensors. For each set, we recorded the highest sampling error resulting from the failure of any single sensor, and plotted this maximum error to evaluate performance (Fig. in Appendix ). The greedy optimal placement resulted in the lowest error for all cases and the sampling errors remained within 30 mK when 15 or more sensors (with a sensor failure) are used. In contrast, for the other placement strategies, even with 20 sensors the maximum error under a single sensor failure reached about 60 mK for exponential placements and exceeded 250 mK for linear placements (Fig. ).
Discussion
We found that the greedy optimal sensor placement for borehole temperature measurements provides more informative data than linear or exponential sensor placements, as it enables measurement of the full borehole temperature profile with higher precision and fewer sensors. Here, we review the factors that affect the optimal sensor placement, examine the associated uncertainties, choice of device error, and discuss the broader implications of our findings.
Numerical uncertainty of greedy optimal sampling algorithm
The greedy optimal sensor placements using Algorithm determines the final sensor placement ZS‾ by an averaging over a larger number of sensor placements that correspond to different sets of initial sensor positions, see Sect. . This averaging is done for robustness reasons. While each individual sensor placement ZS fulfills optimality properties, their geometric average does not have to be optimal. Therefore, it is of interest to analyze how strongly the averaged sensor placement deviates from the best possible, non-averaged sensor placement and whether averaging over a larger set influences the result.
Numerical uncertainty handling of greedy optimal sampling. The obtained maximum sampling error is plotted against the number of greedy optimal sensor placement sets (ZS), over which the final placement is averaged. The dashed lines represents the best-case sampling error for EDML and GRIP which is the minimum of the maximum sampling errors across the 1000 distinct sensor placement.
We analyze the maximum sampling error of the greedy optimal sensor placement averaged with 1, 20, 250, 500, and 1000 sets of sensor placements. The best case of sampling error is calculated as the minimum of the maximum sampling errors of each of the 1000 distinct sensor placement sets. The maximum sampling error of averaged greedy optimal sensor placement is observed to be decreasing (towards the best case) with increasing the number of sets used in averaging (Fig. ). Although the maximum sampling error remains significantly higher than the best-case error, the absolute differences, at less than 0.3 mK, are negligible compared to the variations between different sampling strategies or the device error.
Sensitivity on past climate and ice borehole site properties
To better understand the factors influencing optimal sensor spacing, we conducted a sensitivity study to examine the impact of extreme surface temperature histories and site conditions on our results. In addition to the “realistic” parameter sets used previously, we include one case where the surface temperature history exhibits no trend, only stochastic variability (PA=0, β=0.6), and another case dominated by a warming trend (PA=3, β=0.6), which is stronger than trends currently assumed in the literature or found in state-of-the-art climate model simulations .
Greedy optimal sensor placements' sensitivity to the surface temperature time series and the accumulation rate that determines the advection. The distribution of the 20 sensor locations computed for different extreme cases of input temperature time series with respect to EDML are shown in (a). Greedy optimal sensor placements for EDML with ws=64 mm yr−1 and without advection (w=0), and GRIP with its given ws as 230 mm yr−1 and very high advection (ws=1000 mm yr−1) are shown in (b). Impact of advection is analyzed using surface temperature time series with PA=2 and β=0.6 (see Sect. ).
In contrast to our earlier analysis, we present the full distributions of sensor placements for each of the different types of surface temperature time series. We modified the Algorithm such that we generate greedy optimal sensor placements (for computational efficiency reasons, we only averaged over 20 sets of ZS in this analysis) for each of the 1000 surface temperature time series of a particular type and thus the distribution of sensor placements with respect to that particular type of surface temperature time series is obtained.
Interestingly, the distributions of 1000 sensor locations generated for the six types of surface temperature time series mostly overlap at the respective sensor indices (Fig. a). Even the extreme parameter sets for surface temperature have minimal influence on the results. This insensitivity indicates that sensor placement is primarily determined by the properties of the advection-diffusion equation rather than by any specific surface temperature history.
The primary factor varying between different ice core locations is the snow accumulation rate, which ranges from as little as 30 mm yr−1 on the East Antarctic Plateau to more than 1 m yr−1 in coastal regions. This directly affects the advection term. Additionally, site conditions influence the diffusion term by affecting the density and temperature profiles.
To investigate these influences, we conduct another sensitivity study. In addition to the actual cooler Antarctic EDML and warmer Greenlandic GRIP conditions, we examine two extreme parameter sets: EDML without advection and GRIP with very high advection (ws=1000 mm yr−1). For surface temperature, we use the standard parameter set (β=0.6, PA=2.0), as we have shown above that the results are not sensitive to this choice.
We find that for all surface conditions studied, the overall shape of the sensor placement is similar (Fig. b). However, the accumulation rate has a noticeable effect: the optimal sensor locations, as determined by the greedy algorithm in Sect. , shift deeper into the ice sheet as the accumulation rate increases in the borehole simulations (Fig. b). In the extreme case of optimal sensor placement with no advection and extremely high advection, the maximum difference reaches up to 6 m. This behavior aligns with the observation that GRIP sensors are placed slightly deeper in the borehole compared to the corresponding EDML sensors (Fig. a).
The uncertainty in bottom boundary condition parameterization of ice sheet have minimal impact in borehole temperature at 200 m. The borehole temperature profile simulated using a Dirichlet boundary condition differs slightly from that obtained using a Neumann boundary, but the effect on optimal sensor placements are undetectable. Greedy optimal sensor placements do not depend on the absolute temperature, but on the variations of the temperature across borehole depth. Therefore, according to our methodology, the difference in borehole temperature profile simulations due to different basal boundary conditions of the ice sheet, does not affect the sampling error, as it is not contributing to the uncertainty budget we are assessing.
Device error range in ice borehole thermistors
The device error values assumed in this study (5 and 10 mK) fall within the typical range reported for borehole thermistors in glaciological applications. For example, precision values range from around 0.5 mK when using a single thermistor on a winch in high-end systems to approximately 30 mK in simpler thermistor chain setups while accuracy typically lies in the range of 3–30 mK. Our results show that greedy optimal sampling can reduce the sampling error to a level comparable to these typical device errors, particularly when only a limited number of sensors are used.
Implications
The proposed ice borehole sensor placement method, utilizing greedy optimal sampling, utilizes the understanding of how heat is distributed in the ice borehole and incorporates our knowledge on potential past climate states.
This approach has broad applicability and can directly be applied for shallow ice borehole thermometry spanning the Holocene period. Given the robustness of our results to variations in site-specific conditions and climate history (e.g., stochastic variability or pronounced trends), already using the proposed sensor spacing is likely to offer substantial improvements over ad hoc placements, even when applied to locations other than those explicitly analyzed in this study. Rerunning our algorithm for specific sites and specific prior knowledge on potential past climate histories is straightforward and has the advantage of also providing the uncertainty caused by using only a finite set of sensors.
For boreholes extending to the bottom of an ice sheet, for example for the estimation of geothermal heat flux , the spacing will change near the basal regions and the sensitivity to prior information on the bottom heat flux needs to be tested. For studies encompassing time periods extending beyond the Holocene, such as glacial-interglacial cycles and orbital-scale variations, a revised set of surface temperature histories should be used to account for these longer-term climatic oscillations and depending on the site conditions, the forward model has to be extended with an ice-flow model .
The greedy optimal approach to sensor placement can further be extended to other problems such as terrestrial boreholes or boreholes in permafrost areas or the measurement of temperature profiles in sea ice using adapted forward heat transfer models. Expected phase changes (e.g., in permafrost) or vertically heterogeneous thermal properties could be incorporated into the forward model, allowing sensor positions to be optimized nonetheless.
Conclusion
We propose a greedy strategy for placing sensors in ice boreholes to maximize the informativeness of temperature sensor measurements. Unlike ad hoc, linear and exponential sensor placements, where sensors are positioned as a function of depth, our method determines sensor locations by accounting for the subsurface heat distribution influenced by surface temperature diffusion and advection. Exponential sensor placements, which space sensors farther apart down the ice borehole, are found to outperform linear sensor placements. This is because stronger vertical gradients in heat distribution occur closer to the surface and gradually diminish deeper into the ice sheet. Compared to exponential sensor placement, our simulations demonstrate that the greedy optimal sensor placement increases the informativeness of ice borehole measurements, reducing error by a factor of approximately 10.
Additionally, we analyzed sensor placement while considering the inherent measurement uncertainty (device error). This shows that, under realistic settings, greedy optimal sampling requires only half the number of sensors to minimize sampling error compared to linear and exponential placements. Therefore, we argue that our proposed approach could enable more cost-effective and accurate measurement of ice borehole profiles, making it a promising contribution to the establishment of a continuous subsurface temperature monitoring system.
Heat transfer model parameters
The thermal properties and vertical velocity of the firn/ice is calculated with respect to site specific parameters, which are mean surface temperature (θm), basal temperature (θb), basal velocity (wb), accumulation rate (ws) and thickness of the ice sheet (H).
Thermal properties
The thermal properties include density (ρ), specific heat capacity (c), thermal conductivity (K), and thermal diffusivity (k). The density profile is simulated using the Herron–Langway model , which requires surface snow density (ρs) in addition to θm, ws and H. The unit of density is kg m−3 and it is one of the key inputs to the heat transfer model to calculate heat capacity c and thermal conductivity K. The equations for their calculations are listed below .
Specific heat capacity (<inline-formula><mml:math id="M255" display="inline"><mml:mi>c</mml:mi></mml:math></inline-formula>)
Specific heat capacity of ice (cice) in J kg−1 K−1 is given by the following equation from ,
cice=152.5+7.122T,
where T is the temperature in Kelvin. Specific heat capacity of firn (cfirn) is calculated from the percentage of ice and air in firn which is given by,
cfirn=ciceρρice+ca1-ρρice,
where ρice is the density of ice (917 kg m−3) and ca (1005 J kg−1 K−1) is the specific heat capacity of dry air.
The temperature-dependent thermal conductivity of ice (Kice) in W m−1 K−1 is given by
Kice=2.22(1-0.0067T),
where T is temperature in °C. Thermal conductivity of firn (Kfirn) is given by
Kfirn=Kice(ρρice)α′-β′(ρρice),
where α′ and β′ are site specific coefficients. We used α′=2.4634 and β′=0 for EDML which was used by for the site NUS07-2, and α′=2 and β′=0.5 for GRIP .
As in , the velocity profile is computed using the equation from , which is based on the ice velocity model by . Using a relative vertical coordinate ζ=z/H, where H is the ice sheet thickness, w(ζ) is given as
w(ζ)=ρρice[ws-(ws-wb)(m+2m+1ζ)⋅(1-ζm+1m+2)].
The constants ws and wb are the vertical velocity at the surface and at the base of the ice sheet, respectively. ws is assumed to be equal to the accumulation rate and wb is the basal melting rate. m is the shape parameter of the vertical velocity profile. We use m=11 for EDML, and m=10 for GRIP.
Sensor placements for GRIP
Performance comparison of linear, exponential, and greedy optimal sensor placements in case of GRIP, assuming no device error. Panels (a)–(c) are based on linear spacing, (d)–(f) are based on exponential spacing and, (g)–(i) are based on greedy optimal sampling of a set of 20 sensors in a 200 m borehole. The horizontal axis of (a), (d) and (g) shows the index of the sensors, sensor with index-1 is placed at the top and index-20 is placed at the bottom of the borehole. Panels (b), (e) and (h) show the sampling error due to linear, exponential, and greedy optimal sensor placement respectively. Panels (c), (f) and (i) show sampling error in logarithmic scale. The dashed line marks 1 mK in (c), (f), and (i).
The sampling of linear, exponential and greedy optimal sensor placements with and without device error were computed for 5 to 20 sensors for GRIP. Maximal sampling errors without device error (σd=0) (a), with 5 mK device error (σd=5 mK) (b) and with 10 mK device error (σd=10 mK) (c) are given. The number of sensors used in sensor placement is denoted on the horizontal axis of (a), (b) and (c).
Interpolation function
The maximum sampling error of exponential and greedy optimal sensor placements with respect to different types of interpolation functions used is shown here. EDML site parameters are used for this analysis.
Sensor failure test
This figure shows the maximal sampling error observed due to single sensor failure for all the three sensor placement strategies of 11 to 20 sensors, each with a device error of standard deviation σd=10mK. The horizontal axis shows the no. of sensors including the failed sensor.
Code and data availability
The code and data used in this study are available on Zenodo at https://doi.org/10.5281/zenodo.17849760 .
Author contributions
PZ, KS, and TL conceptualized and designed the study. NH and TL provided surface temperature surrogates and glaciological expertise. PZ supervised the mathematical aspects of the project, while TL supervised the climate-related components. KS conducted the analysis and drafted the manuscript, with contributions from all co-authors.
Competing interests
The contact author has declared that none of the authors has any competing interests.
Disclaimer
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.
Acknowledgements
We would like to acknowledge the support of the ”Interdisciplinary Center for Machine Learning and Data Analytics (IZMD)” at the University of Wuppertal and we would like to thank Andrew Dolman and Thomas Münch for fruitful discussions.
Financial support
KS is funded through the Helmholtz School for Marine Data Science (MarDATA) under grant no. HIDSS-0005. KS and NH received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement no. 716092)The article processing charges for this open-access publication were covered by the Alfred-Wegener-Institut. Helmholtz-Zentrum für Polar- und Meeresforschung.
Review statement
This paper was edited by Lev Eppelbaum and reviewed by two anonymous referees.
ReferencesAndersson, T. R., Bruinsma, W. P., Markou, S., Requeima, J., Coca-Castro, A., Vaughan, A., Ellis, A.-L., Lazzara, M. A., Jones, D., Hosking, S., and Turner, R. E.: Environmental sensor placement with convolutional Gaussian neural processes, Environmental Data Science, 2, e32, 10.1017/eds.2023.22, 2023.Arthern, R. J., Vaughan, D. G., Rankin, A. M., Mulvaney, R., and Thomas, E. R.: In situ measurements of Antarctic snow compaction compared with predictions of models, Journal of Geophysical Research: Earth Surface, 115, 10.1029/2009JF001306, 2010.Bartos, M. and Kerkez, B.: Observability-Based Sensor Placement Improves Contaminant Tracing in River Networks, Water Resources Research, 57, 10.1029/2020WR029551, 2021.Beltrami, H.: Earth's Long-Term Memory, Science, 297, 206–207, 10.1126/science.1074027, 2002.Beltrami, H., Matharoo, G. S., and Smerdon, J. E.: Impact of borehole depths on reconstructed estimates of ground surface temperature histories and energy storage, Journal of Geophysical Research: Earth Surface, 120, 763–778, 10.1002/2014JF003382, 2015.Biskaborn, B. K., Smith, S. L., Noetzli, J., Matthes, H., Vieira, G., Streletskiy, D. A., Schoeneich, P., Romanovsky, V. E., Lewkowicz, A. G., Abramov, A., Allard, M., Boike, J., Cable, W. L., Christiansen, H. H., Delaloye, R., Diekmann, B., Drozdov, D., Etzelmüller, B., Grosse, G., Guglielmin, M., Ingeman-Nielsen, T., Isaksen, K., Ishikawa, M., Johansson, M., Johannsson, H., Joo, A., Kaverin, D., Kholodov, A., Konstantinov, P., Kröger, T., Lambiel, C., Lanckman, J.-P., Luo, D., Malkova, G., Meiklejohn, I., Moskalenko, N., Oliva, M., Phillips, M., Ramos, M., Sannel, A. B. K., Sergeev, D., Seybold, C., Skryabin, P., Vasiliev, A., Wu, Q., Yoshikawa, K., Zheleznyak, M., and Lantuit, H.: Permafrost is warming at a global scale, Nat. Commun., 10, 264, 10.1038/s41467-018-08240-4, 2019.Black, P. E.: Greedy algorithm, Dictionary of Algorithms and Data Structures, 2, 62, US National Institute of Standards and Technology, https://www.nist.gov/dads/HTML/greedyalgo.html (last access: 18 November 2025), 2005.Bolzan, J. F. and Strobel, M.: Accumulation-rate variations around Summit, Greenland, Journal of Glaciology, 40, 56–66, 10.3189/S0022143000003798, 1994.Casado, M., Münch, T., and Laepple, T.: Climatic information archived in ice cores: impact of intermittency and diffusion on the recorded isotopic signal in Antarctica, Clim. Past, 16, 1581–1598, 10.5194/cp-16-1581-2020, 2020.Clow, G. D.: The extent of temporal smearing in surface-temperature histories derived from borehole temperature measurements, Global and Planetary Change, 6, 81–86, 10.1016/0921-8181(92)90027-8, 1992.Clow, G. D., Saltus, R. W., and Waddington, E. D.: A new high-precision borehole-temperature logging system used at GISP2, Greenland, and Taylor Dome, Antarctica, Journal of Glaciology, 42, 576–584, 10.3189/S0022143000003555, 1996.Colgan, W., Shields, C., Talalay, P., Fan, X., Lines, A. P., Elliott, J., Rajaram, H., Mankoff, K., Jensen, M., Backes, M., Liu, Y., Wei, X., Karlsson, N. B., Spanggård, H., and Pedersen, A. Ø.: Design and performance of the Hotrod melt-tip ice-drilling system, Geosci. Instrum. Method. Data Syst., 12, 121–140, 10.5194/gi-12-121-2023, 2023.Compo, G., Whitaker, J., Sardeshmukh, P., Allan, R., McColl, C., Yin, X., Giese, B., Vose, R., Matsui, N., Ashcroft, L., Auchmann, R., Benoy, M., Bessemoulin, P., Brandsma, T., Brohan, P., Brunet, M., Comeaux, J., Cram, T., Crouthamel, R., Groisman, P., Hersbach, H., Jones, P., Jonsson, T., Jourdain, S., Kelly, G., Knapp, K., Kruger, A., Kubota, H., Lentini, G., Lorrey, A., Lott, N., Lubker, S., Luterbacher, J., Marshall, G., Maugeri, M., Mock, C., Mok, H., Nordli, O., Przybylak, R., Rodwell, M., Ross, T., Schuster, D., Srnec, L., Valente, M., Vizi, Z., Wang, X., Westcott, N., Woollen, J., and Worley, S.: NOAA/CIRES Twentieth Century Global Reanalysis Version 2c, NSF National Center for Atmospheric Research [data set], 10.5065/D6N877TW, 2015. Cuffey, K. M. and Paterson, W. S. B.: The Physics of Glaciers, Academic Press, ISBN 978-0-08-091912-6, 2010.Cuffey, K. M., Clow, G. D., Steig, E. J., Buizert, C., Fudge, T. J., Koutnik, M., Waddington, E. D., Alley, R. B., and Severinghaus, J. P.: Deglacial temperature history of West Antarctica, Proceedings of the National Academy of Sciences, 113, 14249–14254, 10.1073/pnas.1609132113, 2016.Eppelbaum, L. V.: The Possibility of Estimating the Permafrost's Porosity In Situ in the Hydrocarbon Industry and Environment, Geosciences, 14, 72, 10.3390/geosciences14030072, 2024.García-Pereira, F., González-Rouco, J. F., Schmid, T., Melo-Aguilar, C., Vegas-Cañas, C., Steinert, N. J., Roldán-Gómez, P. J., Cuesta-Valero, F. J., García-García, A., Beltrami, H., and de Vrese, P.: Thermodynamic and hydrological drivers of the soil and bedrock thermal regimes in central Spain, SOIL, 10, 1–21, 10.5194/soil-10-1-2024, 2024.Goujon, C., Barnola, J.-M., and Ritz, C.: Modeling the densification of polar firn including heat diffusion: Application to close-off characteristics and gas isotopic fractionation for Antarctica and Greenland sites, Journal of Geophysical Research: Atmospheres, 108, 10.1029/2002JD003319, 2003.Groenke, B., Langer, M., Miesner, F., Westermann, S., Gallego, G., and Boike, J.: Robust Reconstruction of Historical Climate Change From Permafrost Boreholes, Journal of Geophysical Research: Earth Surface, 129, e2024JF007734, 10.1029/2024JF007734, 2024.Hammer, C. U. and Dahl-Jensen, D.: GRIP Accumulation Rates, PANGAEA [data set], 10.1594/PANGAEA.55084, 1999.Hébert, R., Herzschuh, U., and Laepple, T.: Millennial-scale climate variability over land overprinted by ocean temperature fluctuations, Nat. Geosci., 15, 899–905, 10.1038/s41561-022-01056-4, 2022.Herron, M. M. and Langway Jr., C. C.: Firn Densification: An Empirical Model, Journal of Glaciology, 25, 373–385, 10.3189/S0022143000015239, 1980.Hersbach, H., Bell, B., Berrisford, P., Hirahara, S., Horányi, A., Muñoz-Sabater, J., Nicolas, J., Peubey, C., Radu, R., Schepers, D., Simmons, A., Soci, C., Abdalla, S., Abellan, X., Balsamo, G., Bechtold, P., Biavati, G., Bidlot, J., Bonavita, M., De Chiara, G., Dahlgren, P., Dee, D., Diamantakis, M., Dragani, R., Flemming, J., Forbes, R., Fuentes, M., Geer, A., Haimberger, L., Healy, S., Hogan, R. J., Hólm, E., Janisková, M., Keeley, S., Laloyaux, P., Lopez, P., Lupu, C., Radnoti, G., de Rosnay, P., Rozum, I., Vamborg, F., Villaume, S., and Thépaut, J.-N.: The ERA5 global reanalysis, Quarterly Journal of the Royal Meteorological Society, 146, 1999–2049, 10.1002/qj.3803, 2020.Johnsen, S. J.: GRIP Temperature Profile, PANGAEA [data set], 10.1594/PANGAEA.89007, 2003.Jones, J. M., Gille, S. T., Goosse, H., Abram, N. J., Canziani, P. O., Charman, D. J., Clem, K. R., Crosta, X., de Lavergne, C., Eisenman, I., England, M. H., Fogt, R. L., Frankcombe, L. M., Marshall, G. J., Masson-Delmotte, V., Morrison, A. K., Orsi, A. J., Raphael, M. N., Renwick, J. A., Schneider, D. P., Simpkins, G. R., Steig, E. J., Stenni, B., Swingedouw, D., and Vance, T. R.: Assessing recent trends in high-latitude Southern Hemisphere surface climate, Nature Climate Change, 6, 917–926, 10.1038/nclimate3103, 2016.Krause, A., Singh, A., and Guestrin, C.: Near-Optimal Sensor Placements in Gaussian Processes: Theory, Efficient Algorithms and Empirical Studies, Journal of Machine Learning Research, 9, 235–284, 2008. Laepple, T. and Huybers, P.: Ocean surface temperature variability: Large model–data differences at decadal and longer periods, Proceedings of the National Academy of Sciences, 111, 16682–16687, 10.1073/pnas.1412077111, 2014.Lee, S.-y., Lee, I.-b., Yeo, U.-h., Kim, R.-w., and Kim, J.-g.: Optimal sensor placement for monitoring and controlling greenhouse internal environments, Biosystems Engineering, 188, 190–206, 10.1016/j.biosystemseng.2019.10.005, 2019.Ligtenberg, S. R. M., Helsen, M. M., and van den Broeke, M. R.: An improved semi-empirical model for the densification of Antarctic firn, The Cryosphere, 5, 809–819, 10.5194/tc-5-809-2011, 2011. Lliboutry, L.: A critical review of analytical approximate solutions for steady state velocities and temperature in cold ice-sheets,, Z. Gletscherkde. Glazialgeol, 15, 135–148, 1979.Løkkegaard, A., Mankoff, K. D., Zdanowicz, C., Clow, G. D., Lüthi, M. P., Doyle, S. H., Thomsen, H. H., Fisher, D., Harper, J., Aschwanden, A., Vinther, B. M., Dahl-Jensen, D., Zekollari, H., Meierbachtol, T., McDowell, I., Humphrey, N., Solgaard, A., Karlsson, N. B., Khan, S. A., Hills, B., Law, R., Hubbard, B., Christoffersen, P., Jacquemart, M., Seguinot, J., Fausto, R. S., and Colgan, W. T.: Greenland and Canadian Arctic ice temperature profiles database, The Cryosphere, 17, 3829–3845, 10.5194/tc-17-3829-2023, 2023.Montelli, A. and Kingslake, J.: Geothermal heat flux is the dominant source of uncertainty in englacial-temperature-based dating of ice rise formation, The Cryosphere, 17, 195–210, 10.5194/tc-17-195-2023, 2023.Münch, T. and Laepple, T.: What climate signal is contained in decadal- to centennial-scale isotope variations from Antarctic ice cores?, Clim. Past, 14, 2053–2070, 10.5194/cp-14-2053-2018, 2018.Muto, A.: Multi-decadal surface temperature trends in East Antarctica inferred from borehole firn temperature measurements and geophysical inverse methods, PhD Thesis, University of Colorado at Boulder, https://scholar.colorado.edu/concern/graduate_thesis_or_dissertations/td96k2585 (last access: 8 December 2025), 2010.Muto, A., Scambos, T. A., Steffen, K., Slater, A. G., and Clow, G. D.: Recent surface temperature trends in the interior of East Antarctica from borehole firn temperature measurements and geophysical inverse methods, Geophysical Research Letters, 38, 10.1029/2011GL048086, 2011.Oerter, H., Wilhelms, F., Jung-Rothenhäusler, F., Göktas, F., Miller, H., Graf, W., and Sommer, S.: Accumulation rates in Dronning Maud Land, Antarctica, as revealed by dielectric-profiling measurements of shallow firn cores, Annals of Glaciology, 30, 27–34, 10.3189/172756400781820705, 2000.Oerter, H., Graf, W., Meyer, H., and Wilhelms, F.: The EPICA ice core from Dronning Maud Land: first results from stable-isotope measurements, Annals of Glaciology, 39, 307–312, 10.3189/172756404781814032, 2004.Orsi, A. J., Cornuelle, B. D., and Severinghaus, J. P.: Little Ice Age cold interval in West Antarctica: Evidence from borehole temperature at the West Antarctic Ice Sheet (WAIS) Divide: WAIS DIVIDE TEMPERATURE, Geophys. Res. Lett., 39, 10.1029/2012GL051260, 2012.Paterson, W.: Physics of Glaciers, Pergamon/Elsevier Science, 3rd edn., ISBN 9780080379449, 10.1016/C2009-0-14802-X, 1994.Pollack, H. N., Huang, S., and Shen, P.-Y.: Climate Change Record in Subsurface Temperatures: A Global Perspective, Science, 282, 279–281, 10.1126/science.282.5387.279, 1998.Salamatin, A.: Paleoclimatic reconstructions based on borehole temperature measurements in ice sheets: possibilities and limitations, Physics of Ice Core Records, 243–282, http://hdl.handle.net/2115/32471 (last access: 18 November 2025), 2000.Seo, S., Wallat, M., Graepel, T., and Obermayer, K.: Gaussian Process Regression: Active Data Selection and Test Point Rejection, in: Mustererkennung 2000, 27–34, Springer, Berlin, Heidelberg, ISBN 978-3-642-59802-9, 10.1007/978-3-642-59802-9_4, 2000.Shaju, K.: Ice borehole thermometry: Sensor placement using greedy optimal sampling, Zenodo [code, data set], 10.5281/zenodo.17849760, 2025.Shen, P. Y., Pollack, H. N., Huang, S., and Wang, K.: Effects of subsurface heterogeneity on the inference of climate change from borehole temperature data: Model studies and field examples from Canada, Journal of Geophysical Research: Solid Earth, 100, 6383–6396, 10.1029/94JB03136, 1995.Singh, A., Krause, A. R., Guestrin, C., Kaiser, W. J., and Batalin, M.: Efficient Planning of Informative Paths for Multiple Robots, in: 20th International Joint Conference on Artifical Intelligence, 2204–2211, https://escholarship.org/uc/item/4r48w3bb (last access: 18 November 2025), 2006.Thorsteinsson, T., Kipfstuhl, J., and Miller, H.: Textures and fabrics in the GRIP ice core, Journal of Geophysical Research: Oceans, 102, 26583–26599, 10.1029/97JC00161, 1997.von Schuckmann, K., Minière, A., Gues, F., Cuesta-Valero, F. J., Kirchengast, G., Adusumilli, S., Straneo, F., Ablain, M., Allan, R. P., Barker, P. M., Beltrami, H., Blazquez, A., Boyer, T., Cheng, L., Church, J., Desbruyeres, D., Dolman, H., Domingues, C. M., García-García, A., Giglio, D., Gilson, J. E., Gorfer, M., Haimberger, L., Hakuba, M. Z., Hendricks, S., Hosoda, S., Johnson, G. C., Killick, R., King, B., Kolodziejczyk, N., Korosov, A., Krinner, G., Kuusela, M., Landerer, F. W., Langer, M., Lavergne, T., Lawrence, I., Li, Y., Lyman, J., Marti, F., Marzeion, B., Mayer, M., MacDougall, A. H., McDougall, T., Monselesan, D. P., Nitzbon, J., Otosaka, I., Peng, J., Purkey, S., Roemmich, D., Sato, K., Sato, K., Savita, A., Schweiger, A., Shepherd, A., Seneviratne, S. I., Simons, L., Slater, D. A., Slater, T., Steiner, A. K., Suga, T., Szekely, T., Thiery, W., Timmermans, M.-L., Vanderkelen, I., Wjiffels, S. E., Wu, T., and Zemp, M.: Heat stored in the Earth system 1960–2020: where does the energy go?, Earth Syst. Sci. Data, 15, 1675–1709, 10.5194/essd-15-1675-2023, 2023.Weikusat, I., Jansen, D., Binder, T., Eichler, J., Faria, S. H., Wilhelms, F., Kipfstuhl, S., Sheldon, S., Miller, H., Dahl-Jensen, D., and Kleiner, T.: Physical analysis of an Antarctic ice core—towards an integration of micro- and macrodynamics of polar ice*, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 375, 20150347, 10.1098/rsta.2015.0347, 2017.Wesche, C., Weller, R., König-Langlo, G., Fromm, T., Eckstaller, A., Nixdorf, U., and Kohlberg, E.: Neumayer III and Kohnen Station in Antarctica operated by the Alfred Wegener Institute, Journal of large-scale research facilities JLSRF, 2, 1–6, 10.17815/jlsrf-2-152, 2016.Zagorodnov, V., Nagornov, O., Scambos, T. A., Muto, A., Mosley-Thompson, E., Pettit, E. C., and Tyuflin, S.: Borehole temperatures reveal details of 20th century warming at Bruce Plateau, Antarctic Peninsula, The Cryosphere, 6, 675–686, 10.5194/tc-6-675-2012, 2012. Zuo, G., Dou, Y., Chang, X., Chen, Y., and Ma, C.: Design and Performance Analysis of a Multilayer Sea Ice Temperature Sensor Used in Polar Region, Sensors, 18, 10.3390/s18124467, 2018.