Search for: All records

Linear elastic fracture mechanics suggests that short‐lived flow accelerations, such as the one initiated by a flooding event beneath Byrd Glacier in 2006, can form abnormally large basal crevasses at the grounding line. Airborne radar measurements acquired in 2011 reveal hundreds of basal crevasses ranging in height from40–335 m. Particle tracking results show that formation of the largest basal crevasse occurred at the grounding line during the 2006 flooding event. Very large basal crevasses form distinctive surface depressions directly overhead, which are observed along the Byrd Glacier flowline to the terminus of the Ross Ice Shelf. By using these surface depressions as proxies for abnormally large basal crevasses, we create a timeline of past subglacial flooding events on Byrd Glacier. Understanding the frequency of flooding events and their effect on glacier dynamics will help inform models of ice sheet stability and subglacial hydrology.

The one‐dimensional steady state analytical solution of the energy conservation equation obtained by Robin (1955, https://doi.org/10.3189/002214355793702028) is frequently used in glaciology. This solution assumes a linear change in surface velocity from a minimum value equal to minus the mass balance at the surface to zero at the bed. Here we show that this assumption of a linear velocity profile leads to large errors in the calculated temperature profile and especially in basal temperature. By prescribing a nonlinear power function of elevation above the bed for the vertical velocity profile arising from use of the Shallow Ice Approximation, we derive a new analytical solution for temperature. We show that the solution produces temperature profiles identical to numerical temperature solutions with the Shallow Ice Approximation vertical velocity near ice divides. We quantify the importance of strain heating and demonstrate that integrating the strain heating and adding it to the geothermal heat flux at the bed is a reasonable approximation for the interior regions. Our analytical solution does not include horizontal advection components, so we compare our solution with numerical solutions of a two‐dimensional advection‐diffusion model and assess the applicability and errors of the analytical solution away from the ice divide. We show that several parameters and assumptions impact the spatial extent of applicability of the new solution including surface mass balance rate and surface temperature lapse rate. We delineate regions of Greenland and Antarctica within which the analytical solution at any depth is likely within 2 K of the actual temperatures with horizontal advection.