Temperature sensors installed in a grid of nine full-depth boreholes drilled in the southwestern ablation zone of the Greenland Ice Sheet recorded cooling in discrete sections of ice over time within the lowest third of the ice column in most boreholes. Rates of temperature change outpace cooling expected from vertical conduction alone. Additionally, observed temperature profiles deviate significantly from the site-average thermal profile that is shaped by all thermomechanical processes upstream. These deviations imply recent, localized changes to the basal thermal state in the boreholes. Although numerous heat sources exist to add energy and warm ice as it moves from the central divide towards the margin such as strain heat from internal deformation, latent heat from refreezing meltwater, and the conduction of geothermal heat across the ice–bedrock interface, identifying heat sinks proves more difficult. After eliminating possible mechanisms that could cause cooling, we find that the observed cooling is a manifestation of previous warming in near-basal ice. Thermal decay after latent heat is released from freezing water in basal crevasses is the most likely mechanism resulting in the transient evolution of temperature and the vertical thermal structure observed at our site. We argue basal crevasses are a viable englacial heat source in the basal ice of Greenland's ablation zone and may have a controlling influence on the temperature structure of the near-basal ice.
Rates of ice deformation and the onset of basal sliding are regulated by the thermal regime of ice sheets (Hooke, 1981; Weertman, 1957), which makes ice temperature an important physical variable to be measured in the field and modeled to determine current and future ice motion. As ice moves from the central divide towards the margins, heat is conducted into the ice column from the surface and basal boundaries and is introduced through sources such as energy released by internal deformation and refreezing of meltwater (Cuffey and Paterson, 2010).
Mechanisms of heat transfer and heat sources in Greenland's ablation zone have been investigated by comparing sparsely collected temperature profiles to output from thermomechanical models (e.g., Iken et al., 1993; Funk et al., 1994; Lüthi et al., 2002; Harrington et al., 2015; Lüthi et al., 2015; Poinar et al., 2017). When attempting to replicate the vertical temperature structure, models consistently produced ice temperatures that are lower than collected temperatures. Previous modeling and observational studies highlight the appreciable influence of refreezing meltwater as a source of latent heat in Greenland's ablation zone (e.g., Phillips et al., 2013; Harrington et al., 2015; Lüthi et al., 2015; Poinar et al., 2017; Seguinot et al., 2020; Thomsen et al., 1991). Outside of Greenland, observations of englacial warming from latent heat sources have been demonstrated as early as the 1970s (Jarvis and Clarke, 1974).
Here we investigate heat transfer processes using an extensive dataset collected from over 300 temperature sensors emplaced in a grid of nine full-depth boreholes in southwestern Greenland. Our data provided sub-daily continuous records of ice temperature over a 2–3-year period, longer than many similar temperature time series datasets. The high temporal and spatial resolution allowed us to examine how the vertical thermal structure evolves in both time and space during the study period. Our continuous record of englacial temperatures enabled us to identify dominant drivers of heat flow at our field location.
Nine boreholes were drilled in the southwestern portion of the Greenland Ice Sheet in a sector that discharges in the land-terminating outlet glacier, Issunguata Sermia. The field site was located approximately 33 km from the margin and had a surface elevation of 980 m above sea level with ice depth ranging from 640 to 680 m. The dimensions of the borehole grid were approximately 250 m N–S and 500 m E–W. Figure 1 shows the layout of the instrumented field site.
Map of the field location showing the layout of the nine boreholes
(blue dots) and the five GPS stations (black triangles). Borehole nomenclature
is determined by the year that they were drilled and the direction from the
center borehole and GPS station. The boreholes are located
Borehole experiments conducted immediately after drilling indicated the ice
rests on bedrock with perhaps a veneer of sand and gravel approximately a
decimeter thick (Harper et al., 2017). The subglacial drainage system was
indicative of one overlying bedrock and governed by hard-bedded physics
(Harper et al., 2017). Mean annual air temperature measured at the site was
A hot water drill was used to melt boreholes
Our data allowed us to examine both the vertical temperature field within
the ice column and the temporal changes in temperature. The
temperature of a hot-water-drilled borehole decays towards the ambient ice
temperature over time. After 1 year
We also neglected the first year of measured temperatures in the temporal
records of temperature change to eliminate the drilling disturbance. After
this time, the measurements of temperature in the environment isolated from
electrical noise resulted in a time series of temperature changes that is
mostly noise-free and at the absolute limit of the sensor's resolution. The
data were cleaned of the few readings recording more than
We determined the rate of temperature change recorded by a sensor using a least-squares regression of the data recorded over time. Often the first and last digital steps of this record had fewer entries recorded due to the restrictions of the sampling window. Linearly regressing our stepped data to determine the rate of change of temperature required careful treatment to prevent a biased estimate from truncated observation periods. Assuming that temperature change is close to linear with time, we expected that with a longer record, the sensor would record approximately the same number of readings at each resolution step. Short data lengths occurring at either the beginning or end of the records were likely truncated by the restricted time period of our study. To correct for our restricted sampling time, we padded the first and last steps in our data to match the length of a fully recorded resolution step. This was performed by copying the required length of data from a full resolution step and pre- or post-appending it to the truncated data, while equating the temperatures. This procedure is illustrated in Fig. S2, which also shows the difference in derived rates of temperature change over time (Supplement).
The vertical temperature field collected from the nine boreholes 1 year
after temperature sensors were installed contains a cold central core with
temperatures increasing towards the surface and the bed (Fig. 2). Two
boreholes were only instrumented with temperature sensors in the lowest half
of the ice column. Temperatures are lowest in approximately the middle of
the ice column, ranging from
Temperature profiles for each borehole plotted with respect to each sensor's height above the bed. Sensor locations are marked with points, and temperatures between sensors are linearly interpolated. Dashed lines indicate the range of pressure melting temperatures given the various Clausius–Clapeyron gradients presented in Cuffey and Paterson (2010). The black box indicates the extent of the inset showing temperatures in the lowest 20 m.
Rates of temperature change are generally a function of position in the ice column (Fig. 3). We captured a clear cooling signature in the lowest third of the ice column in most boreholes, with some boreholes warming or cooling in locations outside of this region. Most boreholes cool at a rate of a few hundredths of a degree Celsius per year in the lowest portion of the ice column. Only five full-depth profiles are shown, as temperature sensors in boreholes 15N and 14W failed for most of the study period. Boreholes 15CB and 14SB were only instrumented in the lowest portion of the ice column.
Rates of temperature changes for each borehole plotted against each sensor's height above the bed. Boreholes drilled in 2014 show trends between July 2015 and July 2017, while temperature trends in boreholes drilled in 2015 represent temperature changes observed between July 2016 and July 2017. Error bars indicate 2 standard errors from the mean rate, determined by the linear regression.
While the observed temporal temperature changes are small, we believe the trends are a real signal and not an artifact of sensor drift or random noise. Analog-to-digital systems can exhibit both noise and drift over time; however, the expected magnitude of drift or aging is approximately an order of magnitude less than our observed trends. In addition, the coherence and lack of obvious bias in the data gives us confidence. We observe coherent sections of ice that cool (the bottom of 14SA, 14N, and 15CA) while other portions warm (top of 15CA and 14SA). Additionally, borehole pairs 14SA, 14SB and 15CA, 15CB are located within 10 m of each other and show similar rates and extent of cooling over time, establishing repeatability and viability of our measurements and methods.
We focus on observed basal cooling, since cooling of ice near the bed is enigmatic. Multiple heat sources are available to warm ice as it flows from the interior towards the margins. Heat sinks, however, are not obvious, other than vertical thermal conduction from the relatively warm ice at the bed into the cold central core. We proceed by investigating the potential causes of cooling in the lower part of the ice column.
We compare the magnitude of observed cooling over time to the temporal
temperature changes that would be produced from purely vertical thermal
diffusion given the shapes of our collected temperature profiles. We use the
1-D thermal diffusion equation in the vertical dimension:
Comparison of our data and expected temperature changes from the analytical
diffusion equation (Eq. 1) indicates the curvatures in our temperature
profiles cannot produce the observed temporal changes of temperature over
time through diffusive vertical heat flow alone. To illustrate the
discrepancies between calculated temperature changes caused by vertical
thermal diffusion and the observed temporal temperature change, we calculate
the additional sources or sinks of heat required throughout the ice column
to sustain the observed temperature change through time. Figure 4 presents a
vertical distribution of the required
Additional heat fluxes required to maintain observed rates of temperature change if the collected static temperature profiles evolve through vertical thermal conduction over time. The dashed line indicates the switch between a required heat source or sink.
In addition, we also highlight regions of the temperature profiles that are not only warming or cooling anomalously, but that may be warmer or cooler than would be expected given the average vertical temperature field at our site. The shape of the site-averaged temperature profile results from ice integrating changes in boundary conditions and all upstream thermomechanical processes. Any deviation from this averaged profile within each collected temperature profile will highlight thermal regions that reflect changes either to the local boundary conditions or to heat sources or sinks. These changes must have occurred recently for sensors within the closely spaced boreholes to capture temperature fields that do not match the average vertical temperature field at the site.
Discrepancies between each collected profile and the site-average temperature profile are shown in Fig. 5. Generally, there are positive deviations within the lowest portion of the ice column where there is a clear signal of cooling in boreholes 14SA, 14SB, 14N, 15CA, and 15CB, showing temperatures in the lowest portion of these boreholes are higher than those in the site-average profile and cool over time. Encouragingly for this analysis, the temperatures in the upper half of the 14N and 15CA profiles are warmer than the site-average profile, and cooling occurs over time. While the pattern of temporal temperature change in boreholes 15S and 15E is not as clear, 15E is colder throughout the full-ice depth than the mean thermal profile, but this borehole is also where sensors record the most warming throughout the profile. Spatially, the southeast portion of our field site does not have a cooling signature in relatively warm ice near the bed, which is evident throughout the other boreholes that have time series of temperature measurements.
Rates of temperature change for each borehole (black dots) plotted with temperature deviations within each borehole from the site-average temperature profile. Red shaded regions indicate that the collected static profile is warmer than the site-average profile, and blue shaded regions show the static profile is colder than the site-average profile. Error bars on the calculated rates are the same as in Fig. 3.
Pervasive cooling in the lowest portion of boreholes where ice is warmer than the site-average profile in some, but not all, of our borehole data is enigmatic, since multiple heat sources exist near the bed, but identifying heat sinks is problematic. In addition, the non-uniform nature of cooling near the bed requires explanation. We discuss several possible heat sinks below.
Ice near the bed at our field site is at the calculated pressure melting
temperature (Fig. 2). Some melting can occur in this temperate ice as the
pressure field varies due to ice flow, and transient increases in pressure
decrease the melting point. Melting may contribute to cooling near the bed
as heat is needed for the phase change. However, ice only reaches the
pressure melting temperature in the lowest
The ice sheet thins slightly along flow at our field site; the bed elevation
increases while the surface elevation remains relatively constant. With a
lower ice overburden pressure, the pressure melting temperature increases
(Cuffey and Paterson, 2010). Any interstitial water present in the ice would
then be colder than the new melting temperature and would refreeze.
Lüthi et al. (2002) estimated the water content of the temperate layer
of ice near the bed of Jakobshavn Isbræ has a moisture content of
All full-depth temperature profiles display a characteristic cold central core (Fig. 2). Vertical compression caused by longitudinal and transverse extension or basal melting could shift the cold ice downwards, causing temperatures in the lowest 100 m to decrease. Our sensors would fail to detect this cooling because they are mechanically coupled to the same package of ice through time after the borehole freezes. The sensors move with the ice and continue to record temperature of the same section of ice through time.
Alternatively, vertically compressing ice would steepen thermal gradients
and facilitate higher rates of heat flow into the cold core. With an average
surface vertical strain rate of
Variations in the boundary temperatures will propagate into the interior of the ice column. Here we address several of the most likely boundary temperature changes to assess how they could affect the temporal evolution of the ice at our site.
The Little Ice Age in Greenland occurred approximately 500–200 years ago,
and surface temperatures were
The length scale for changes in surface temperatures during the Little Ice
Age 500 years ago is
Equation (2) also provides an estimate of the timescale of diffusion,
illustrating changes in ice temperature from the end of the last
interglacial period would have diffused away through the entire ice column after
The most significant thermal change to the basal boundary is the transition from frozen to temperate, which occurs upstream of our site. This change introduces a warmer basal boundary relatively recently in terms of the ice flow; however, this would create a region of colder ice near the bed that is being warmed by conduction from the warm bed, not cooled as we observe.
The interior ice sheet, which is frozen to the bed, moves largely by internal deformation and accumulates dissipated heat near the bed where strain rates are highest. Ice at our site moves primarily by basal sliding (Maier et al., 2019) and this strain heat source is negligible, although frictional heating still exists at the base. Downstream from the region where this transition from deformation to sliding occurs, the ice near the bed would lose the strain-derived thermal energy and would cool as the heat diffused upward into the cold core. Although the location of the upstream transition from a frozen to thawed bed conducive to basal sliding is not precisely known, it is on the order of 200 km up-glacier from our field site (MacGregor et al., 2016), which translates to several thousand years for ice flow. Any heat anomalies would have decayed using Eq. (2). Additionally, we have shown that the shape of the temperature profiles that would result from previously accumulated strain heat near the base will not produce the magnitude of observed cooling from vertical thermal conduction alone. An additional heat sink must be required.
Although the frozen-to-sliding transition is probably the largest change to the strain heating term, analysis by Meierbachtol et al. (2016) shows a strong change in the basal shear stress regime approximately 75 km upstream of our site, where the basal shear stress decreases down-flow by a factor of 2. Again Eq. (2) illustrates that the thermal impact is negligible at our site, especially since the bulk of the strain, and therefore strain heating, occurs close to the basal boundary, which forms a fixed boundary condition at a short length scale for heat flow. However, this analysis indicates that variations in the basal boundary stress state can cause minor heating or cooling if strain, and therefore dissipated heat, is removed or reduced. We cannot completely eliminate this process, since basal conditions are not well known.
With the lack of obvious heat sinks, we investigate the possibility that the observed cooling in the lower depths of the ice occurs as thermal decay after latent heat is released from freezing water. Temporally and spatially uniform heating by refreezing of disseminated water would not lead to our observations, since horizontal thermal diffusion would not occur in laterally extensive sections of isothermal ice. However, widely spatially separated transient heating events, such as dispersed basal crevasses or moulins, only warm ice locally, leaving significant regions of colder, unheated ice to facilitate horizontal diffusion. Heat flow along horizontal temperature gradients allows the warmed ice to cool back to ambient temperature. Since heat would flow horizontally, not vertically, the temperature changes would not be revealed by our vertical thermal diffusion calculations governed by the shape of the temperature profiles.
A similar pattern of cooling in ice warmer than the site-average thermal profile is found in the upper portion of 14N (Fig. 4). Hills et al. (2017) discussed the thermal anomaly exhibited in borehole 14N; the authors recognized the relatively high temperatures of ice in the upper portion of the borehole when compared to others in the grid (Fig. 2) and noticed that temperature sensors recorded decreasing temperatures over time. Hills et al. (2017) attributed this localized thermal anomaly to a refrozen crevasse that was observed on the ice surface nearby. The lowest 100 m of ice in boreholes 14SA, 14N, and 15CA exhibits nearly identical thermal characteristics to the uppermost portion of 14N: ice is warmer than the site-average profile and it cools through time.
Basal crevasses can have a thermal influence on discrete but large regions of ice. Because these fractures are distinct features, the fact that we do not see a clear cooling signature in relatively warm ice throughout the base of the ice column at our site is not surprising. Slight warming in cold ice in borehole 15E and no clear signal in borehole 15S provide evidence that the ice is being affected by a thermal process that is not distributed in nature.
Multiple fractures can open similar to crevasses at the surface, and their growth is likely vertically restricted above the bed due to the cold, stiff central core of the ice column which reduces fracture propagation. Therefore, these features only affect the lower portion of the ice column. Additionally, basal crevasses likely are not simple planar features; strong deformational gradients near the bed can alter their shape and cause them to intersect more recently formed crevasses. Although basal crevasses may be more complex than an en echelon arrangement of vertical fractures, for simple modeling we focus on the thermal disturbance by basal crevasses that allow the influx of a plane of water from the bed up into the cold ice. The refreezing of this injected water releases heat into the surrounding ice. Accurately modeling the thermal disturbance created by a basal crevasse is hampered by our lack of knowledge of these features within grounded glaciers and ice sheets, since they have only been directly observed under mountain glaciers (Harper et al., 2010). Extensive basal crevasse fields have been identified based on strong hyperbolic reflections in radar surveys near the base of floating ice shelves (e.g., Jezek and Bentley, 1983; McGrath et al., 2012; Luckman et al., 2012) and in ice streams that were once a part of relict ice shelves (Catania et al., 2005). However, the base environment of ice shelves is different than that of grounded glaciers and ice sheets, with conditions characterized by deviatoric tensile stresses and basal water pressures at ice overburden pressure that are conducive for basal crevasse formation (Weertman, 1980), while tidal flexing at the grounding line also contributes to the formation and spacing of basal crevasses within floating ice (Jezek and Bentley, 1983).
We do not attempt detailed modeling of the thermal disturbance of basal
crevasses at our field site, but we can approximate the thermal disturbance
created by a single or by multiple basal crevasses to compare with our data.
One-dimensional thermal modeling in ice surrounding a crevasse requires
knowledge of the width of the crevasse and the presence of other nearby
crevasses, since closely spaced water-filled crevasses greatly reduce the
refreezing rate (Jarvis and Clarke, 1974). In addition, even the refreezing
problem is complicated by the Stefan boundary condition as the thickness of
the crevasse is reduced by refreezing (Carslaw and Jaeger, 1959). We are
interested in the thermal behavior of ice in the region of a basal crevasse
at times much longer than the crevasse remained open, so we model the
crevasse as an instantaneous planar heat source, considering only the total
latent heat released during the refreezing. We examine the thermal
disturbance at a long timescale in the half-space adjacent to the crevasse
using an analytic solution presented by Carslaw and Jaeger (1959):
Here
We assume a crevasse thickness of 1 m, matching the upper limit of previously observed basal crevasse widths (Harper et al., 2010). The release of heat initiates a wave of thermal energy propagating away from the source (Fig. 6a). The diffusion of the warming front causes ice temperature to increase before cooling asymptotically back to the ambient ice temperature (Fig. 6b). After the crevasse refreezes, the remaining period of warming is short compared to the timescale of cooling; thus the likelihood of temperature sensors installed in a previously crevassed area recording warming is low. To compare this theoretical result with our data, we average the rates of cooling observed in the lowest 100 m of the boreholes and show the range of these averages in Fig. 6c. The theoretical curves depend strongly on the original size of the crevasse and on the distance away from the crevasse; however, we have chosen values that show that this process can produce rates of cooling of the same magnitude as observed in the borehole array.
It is tempting to estimate the time since water-filled basal crevasses were
formed using our extensive temperature data. However, such an attempt is
under-constrained since similar thermal signatures could be produced by many
crevasse geometries and arrangements in time and space. Instead we examine
the same crevassing scenario above as a baseline to investigate the amount
of time elapsed since refreezing to produce observed cooling rates given
various distances between boreholes and hypothetical crevasses. Equation (3)
is solved for a range of time and linear distances from the heat source, and
rates of temperature change over time are calculated. Equivalent rates of
cooling could be produced from the dissipation of thermal energy released by
a refrozen crevasse after approximately 20–30 years if the boreholes
directly intersect the crevasse. The time elapsed is less if the boreholes
do not intersect the refrozen crevasse. Although speculative, it is
interesting that using a timescale of 10–30 years, assuming an ice velocity
of 100 m yr
Basal water pressures at or above ice overburden pressure promote the growth
of basal crevasses tens of meters above the bed (van der Veen, 1998).
Longitudinal extension further facilitates crevassing at distances over 100 m above the bed (van der Veen, 1998). Subglacial water pressures in
Greenland's ablation zone have consistently been measured at or near ice
overburden pressure (e.g., Iken et al., 1993; Meierbachtol et al., 2013;
Andrews et al., 2014; Ryser et al., 2014; van de Wal et al., 2015; Doyle et
al., 2018) and have peaked at
It is unlikely that a single basal crevasse is responsible for the cooling that we observe at the field site. More plausibly, the ice at our field site was fractured by multiple crevasses up-glacier. We cannot precisely determine the number or spacing of crevasses; however, a 1-D modeling analysis shows that three 1 m thick refrozen crevasses would lead to dispersed cooling over a large area of the bed after 20–30 years (Supplement). The 1-D temperature field evolution illustrates the warming that occurs due to the release of latent heat immediately as the crevasse completely refreezes, followed by dispersed cooling surrounding the crevasses over time (Fig. 7). Given negative heat fluxes are required in the lowest portion of all boreholes (Fig. 4) to create observed cooling, thermal energy likely diffuses away from our field site in three dimensions, which would result in cooling that we have not captured in our vertical 1-D diffusion calculations. It is plausible that boreholes intersect these regions of ice that have been thermally altered by basal crevassing.
Evolution of a synthetic 1-D temperature field from basal crevasses spaced 100 m apart. Temperatures are plotted every 5 years, and older temperature fields are shaded darker shades of grey. The red line shows the temperature field 30 years after refreezing. Distinct crevasses can cause cooling over large regions of the bed.
Moulins are unlikely to create the observed thermal observations. While the
internal architecture of moulins is likely complex, as a simple model, we
conceptualize moulins as linear pipes that release thermal energy from the
surface to the bed because of how these features are manifested at the
surface. Cooling in relatively warm ice is not observed throughout the
entire ice column in each borehole, which we would expect to see from a
linear heat source extending over the full ice depth. Moulins have been
shown to migrate horizontally with depth (Holmund, 1988), which could
possibly avoid this objection; however, numerous moulins would be needed to
generate the cooling we observe in multiple boreholes. In the ablation zone
surrounding Issunguata Sermia, active moulins have a mean density of 0.1 per
square kilometer (Smith et al., 2015), and Catania and Neumann (2010) found
relict moulins to be spaced
The existence of basal crevasses has not been widely considered outside of alpine or tidally influenced settings; however, these features likely impact the bulk water storage, basal water flow, and ice rheology in Greenland's ablation zone. Harrington et al. (2015) examined full-depth ice temperature data collected in boreholes at a site down-glacier, but on the same flow line as this study, and found temperate layer growth could not be reproduced by thermal modeling incorporating strain heat and longitudinal compression. The authors suggest basal crevasses allow latent heat to be released well above the bed and create the thick temperate layer.
The lack of a thick basal temperate layer indicates the transient thermal influence of basal crevasses. Ice viscosity is highly temperature dependent; thus basal crevasses could create regions of enhanced deformation, although previous observations in Greenland have shown that deformation rates in temperate ice are often lower than theory predicts (e.g., Ryser et al., 2014; Doyle et al., 2018; Maier et al., 2019). In addition, water content in basal temperate ice also influences rates of glacial ice deformation. A previous study at the down-glacier site found a higher water content in the basal temperate ice layer revealed by ground-based radar than is used in modeling studies (Brown et al., 2017). While unable to detect macro-spaces for bulk water storage in the radargrams, the authors speculated that basal crevasses provide macro-spaces for water storage based on the thermal evidence from Harrington et al. (2015), which results in the high water content observed.
Although basal crevasses have only been directly observed in alpine settings, their existence under the Greenland Ice Sheet has strong implications for the basal hydrology. Basal crevasses have been observed creating englacial water flow above the bed (Harper et al., 2010) and even allowing complex three-dimensional crossing water flow paths (Huzurbazar and Humphrey, 2008). These features have the potential to modulate subglacial discharge and efficiently transfer water through the subglacial drainage system.
Temperature measurements that began 1 year after sensor installation in nine boreholes and span 2–3 years show that ice cools in the lowest portion of the ice column, with all boreholes displaying cooling within the lowest third of the ice column. The main thermal sink is conduction of heat from this region up into the cold central core of the ice. Vertical diffusion calculations show this cooling to be insufficient to explain our data. In addition, the lowest section of ice in some boreholes is warmer than the mean thermal profile at our field site, which suggests relatively recent and localized changes to the boundary conditions or near-basal heat sources. The observed cooling regions of the boreholes are often the same regions of ice that are warmer than the site-average profile.
After eliminating various potential heat sinks, we conclude the observed thermal anomalies are created by the cooling of basal ice that has been previously warmed by basal crevassing occurring upstream of our observation site. Latent heat released by refrozen basal crevasses locally warms ice near the bed, but the effects are distributed as thermal energy propagates outwards, resulting in dispersed cooling at rates of the same magnitude as our observations. Our results provide additional thermal evidence for the existence of basal crevasses in Greenland's ablation zone, which has strong relevance for the temperature-dependent rheology of the basal ice and intriguing implications for the complexity of the basal hydrologic system.
Temperature measurements used in this paper are available at
The supplement related to this article is available online at:
NFH, JTH, and TWM designed the field study and installed temperature sensors. IEM performed the analysis and wrote the manuscript with contributions from all co-authors.
The authors declare that they have no conflict of interest.
We thank Megan Thompson-Munson for her helpful discussions throughout this project. We are grateful to Martin Lüthi and Samuel Doyle for their reviews that significantly improved the quality of this paper and to Ginny Catania, who served as the scientific editor.
This research has been supported by the National Science Foundation, Office of Polar Programs (grant nos. 1203451 and 0909495).
This paper was edited by Ginny Catania and reviewed by Martin Lüthi and Samuel Doyle.