Spatial and Temporal Melt Variability at Helheim Glacier, East Greenland, and Its Effect on Ice Dynamics

Understanding the behavior of large outlet glaciers draining the Greenland Ice Sheet is critical for assessing the impact of climate change on sea level rise. The flow of marine-terminating outlet glaciers is partly governed by calving-related processes taking place at the terminus but is also influenced by the drainage of surface runoff to the bed through moulins, cracks, and other pathways. To investigate the extent of the latter effect, we develop a distributed surface-energy-balance model for Helheim Glacier, East Greenland, to calculate surface melt and thereby estimate runoff. The model is driven by data from an automatic weather station operated on the glacier during the summers of 2007 and 2008, and calibrated with independent measurements of ablation. Modeled melt varies over the deployment period by as much as 68% relative to the mean, with melt rates approximately 77% higher on the lower reaches of the glacier trunk than on the upper glacier. We compare melt variations during the summer season to estimates of surface velocity derived from global positioning system surveys. Near the front of the glacier, there is a significant correlation (on >95% levels) between variations in runoff (estimated from surface melt) and variations in velocity, with a 1 day delay in velocity relative to melt. Although the velocity changes are small compared to accelerations previously observed following some calving events, our findings suggest that the flow speed of Helheim Glacier is sensitive to changes in runoff. The response is most significant in the heavily crevassed, fast-moving region near the calving front. The delay in the peak of the cross-correlation function implies a transit time of 12-36 h for surface runoff to reach the bed.

Correction to Variations of Ice Bed Coupling Beneath and Beyond Ice Streams: The Force Balance

A geometrical force balance that links stresses to ice bed coupling along a flow band of an ice sheet was developed in 1988 for longitudinal tension in ice streams and published 4 years later. It remains a work in progress. Now gravitational forces balanced by forces producing tensile, compressive, basal shear, and side shear stresses are all linked to ice bed coupling by the floating fraction phi of ice that produces the concave surface of ice streams. These lead inexorably to a simple formula showing how phi varies along these flow bands where surface and bed topography are known: phi = h(O)/h(I) with h(O) being ice thickness h(I) at x = 0 for x horizontal and positive upslope from grounded ice margins. This captures the basic fact in glaciology: the height of ice depends on how strongly ice couples to the bed. It shows how far a high convex ice sheet (phi = 0) has gone in collapsing into a low flat ice shelf (phi = 1). Here phi captures ice bed coupling under an ice stream and h(O) captures ice bed coupling beyond ice streams.

Late Pleistocene Interactions of East and West Antarctic Ice-Flow Regimes: Evidence from the McMurdo Ice Shelf

We present new interpretations of deglaciation in McMurdo Sound and the western Ross Sea, with observationally based reconstructions of interactions between East and West Antarctic ice at the last glacial maximum (LGM), 16 000, 12 000, 8000 and 4000 sp. At the LGM? East Antarctic ice from Mulock Glacier split, one branch turned westward south of Ross Island but the other branch rounded Ross Island before flowing southwest into McMurdo Sound. This flow regime, constrained by an ice saddle north of Ross Island, is consistent with the reconstruction of Stuiver and others (1981a). After the LGM, grounding-line retreat was most rapid in areas with greatest water depth, especially along the Victoria Land coast. By 12 000 sp, the ice-now regime in McMurdo Sound changed to through-flowing Mulock Glacier ice, with lesser contributions from Koettlitz, Blue and Ferrar Glaciers, because the former ice saddle north of Ross Island was replaced by a dome. The modern flew regime was established similar to 4000 BP. Ice derived from high elevations on the Polar Plateau but now stranded on the McMurdo Ice Shelf, and the pattern of the Transantarctic Mountains erratics support our reconstructions of Mulock Glacier ice rounding Minna Bluff but with all ice from Skelton Glacier ablating south of the bluff. They are inconsistent with Drewry's (1979) LGM reconstruction that includes Skelton Glacier ice in the McMurdo-Sound through-flow. Drewry's (1979) model closely approximates our results for 12 000-4000 BP. Ice-sheet modeling holds promise for determining whether deglaciation proceeded by grounding-line retreat of an ice sheet that was largely stagnant, because it never approached equilibrium flowline profiles after the Ross Ice Shelf, grounded, or of a dynamic ice sheet with flowline profiles kept low by active ice streams that extended northward from present-day outlet glaciers after the Ross Ice Shelf grounded.

Recent Increase in Black Carbon Concentrations from a Mt. Everest Ice Core Spanning 1860-2000 AD

A Mt. Everest ice core spanning 1860-2000 AD and analyzed at high resolution for black carbon (BC) using a Single Particle Soot Photometer (SP2) demonstrates strong seasonality, with peak concentrations during the winter-spring, and low concentrations during the summer monsoon season. BC concentrations from 1975-2000 relative to 1860-1975 have increased approximately threefold, indicating that BC from anthropogenic sources is being transported to high elevation regions of the Himalaya. The timing of the increase in BC is consistent with BC emission inventory data from South Asia and the Middle East, however since 1990 the ice core BC record does not indicate continually increasing BC concentrations. The Everest BC and dust records provide information about absorbing impurities that can contribute to glacier melt by reducing the albedo of snow and ice. There is no increasing trend in dust concentrations since 1860, and estimated surface radiative forcing due to BC in snow exceeds that of dust in snow. This suggests that a reduction in BC emissions may be an effective means to reduce the effect of absorbing impurities on snow albedo and melt, which affects Himalayan glaciers and the availability of water resources in major Asian rivers. Citation: Kaspari, S. D., M. Schwikowski, M. Gysel, M. G. Flanner, S. Kang, S. Hou, and P. A. Mayewski (2011), Recent increase in black carbon concentrations from a Mt. Everest ice core spanning 1860-2000 AD, Geophys. Res. Lett., 38, L04703, doi: 10.1029/2010GL046096.

Elastic Uplift in Southeast Greenland Due to Rapid Ice Mass Loss

The rapid unloading of ice from the southeastern sector of the Greenland ice sheet between 2001 and 2006 caused an elastic uplift of similar to 35 mm at a GPS site in Kulusuk. Most of the uplift results from ice dynamic-induced volume losses on two nearby outlet glaciers. Volume loss from Helheim Glacier, calculated from sequential digital elevation models, contributes about similar to 16 mm of the observed uplift, with an additional similar to 5 mm from volume loss of Kangerdlugssuaq Glacier. The remaining uplift signal is attributed to significant melt-induced ice volume loss from the ice sheet margin along the southeast coast between 62 degrees N and 66 degrees N.

Theoretical Calving Rates from Glaciers Along Ice Walls Grounded in Water of Variable Depths

Calving has been studied for glaciers ranging from slow polar glaciers that calve on dry land, such as on Deception Island (63.0-degrees-S, 60.6-degrees-W) in Antarctica, through temperate Alaskan tide-water glaciers, to fast outlet glaciers that float in fiords and calve in deep water, such as Jakobshavns Isbrae (69.2-degrees-N, 49.9-degrees-W) in Greenland. Calving from grounded ice walls and floating ice shelves is the main ablation mechanism for the Antarctic and Greenland ice sheets, as it was along marine and lacustrine margins of former Pleistocene ice sheets, and is for tide-water and polar glaciers. Yet, the theory of ice calving is underdeveloped because of inherent dangers in obtaining field data to test and constrain calving models. An attempt is made to develop a calving theory for ice walls grounded in water of variable depth, and to relate slab calving from ice walls to tabular calving from ice shelves. A calving law is derived in which calving rates from ice walls are controled by bending creep behind the ice wall, and depend on wall height h, forward bending angle-theta, crevasse distance c behind the ice wall and depth d of water in front of the ice wall. Reasonable agreement with calving rates reported by Brown and others (1982) for Alaskan tide-water glaciers is obtained when c depends on wall height, wall height above water and water depth. More data are needed to determine which of these dependencies is correct. A calving ratio c/h is introduced to understand the transition from slab calving to tabular calving as water deepens and the calving glacier becomes afloat.

GPR Reflection Profiles of Clark and Commonwealth Glaciers, Dry Valleys, Antarctica

Englacial horizons deeper than 100 m are absent within 100 MHz ground-penetrating radar (GPR) surface profiles we recorded on Clark and Commonwealth Glaciers in the Antarctic Dry Valleys region. Both glaciers show continuous bottom horizons to 280 m, with bottom signal-to-noise ratios near 30 dB. Density horizons should fade below 50 m depth because impermeable ice occurred by 36 m. Folding within Commonwealth Glacier could preclude radar strata beneath about 80 m depth, but there is no significant folding within Clark Glacier. Strong sulfate concentrations and contrasts exist in our shallow ice core. However, it appears that high background concentration levels, and possible decreased concentration contrasts with depth placed their corresponding reflection coefficients at the limit of, or below, our system sensitivity by about 77 m depth. Further verification of this conclusion awaits processing of our deep-core chemistry profiles.

Local Rates of Ice-Sheet Thickness Change in Greenland

The rate of ice-sheet thickness change is calculated for 10 sites in Greenland by comparing measured values of ice vertical velocity and snow-accumulation rate. Vertical velocities are derived from repeat surveys of markers using precision global positioning system techniques, and accumulation rates are determined from stratigraphic analysis of firn cores. The results apply to time-scales covered by the firn-core records, which in most cases are a few decades. A spectrum of thickness-change rates is obtained, ranging from substantial thinning to slow thickening. The sites where ice-sheet thinning is indicated are located near the ice-sheet margin or in outlet glacier catchments. Interior and high-elevation sites are predominantly in balance or thickening slowly. Uncertainties in the rates of thickness change are dominated by errors in the determination of accumulation rates. The results of this work are broadly comparable with regional estimates of mass balance obtained from the analysis of catchment input vs discharge.

Multi-Decadal Record of Ice Dynamics on Daugaard Jensen Gletscher, East Greenland, from Satellite Imagery and Terrestrial Measurements

The history of ice velocity and calving front position of Daugaard Jensen Gletscher, a large outlet glacier in East Greenland, is reconstructed from field measurements, aerial photography and satellite imagery for the period 1950-2001. The calving terminus of the glacier has remained in approximately the same position over the past similar to 50 years. There is no evidence of a change in ice motion between 1968 and 2001, based on a comparison of velocities derived from terrestrial surveying and feature tracking using sequential satellite images. Estimates of flux near the entrance to the fjord vs snow accumulation in the interior catchment show that Daugaard Jensen Gletscher has a small negative mass balance. This result is consistent with other mass-balance estimates for the inland region of the glacier.

Twentieth Century Increase of Atmospheric Ammonia Recorded in Mount Everest Ice Core

An NH4+ record covering the period A.D. 1845-1997 was reconstructed using an 80.4 m ice core from East Rongbuk Glacier at an elevation of 6450 m on the northern slope of Mount Everest. Variations in NH4+ are characterized by a dramatic increase since the 1950s. The highest NH4+ concentrations occur in the 1980s. They are about twofold more than those in the first half of twentieth century. Empirical orthogonal function (EOF) analysis on the eight major ion (Na+,K+,Mg2+,NH4+,Ca2+,NO3-,SO42- and Cl-) series from this core indicates that NH4+ is loaded mainly on EOF3 (60% of NH4+ variance), suggesting that NH4+ has a unique signature. Instrumental sea level pressure (SLP) and regional temperatures are used to explore the relationship between NH4+ variations and both atmospheric circulation and natural source strength over Asia. Higher NH4+ concentrations are associated with an enhanced winter Mongolian High and a deepened summer Mongolian Low. A positive relationship also exists between NH4+ concentrations and regional temperature changes of the GIS Box 36 (Indian subcontinent), indicating that an increase in temperature may contribute to the strengthening of natural ammonia emissions (e. g., from plants and soils). A close positive correlation between NH4+ and acidic species (SO42- plus NO3-) concentrations suggests that a portion of the increase in NH4+ concentrations could be contributed by enhanced atmospheric acidification. Anthropogenic ammonia emissions from enhanced agricultural activities and energy consumption over Asia in concert with population increase since the 1950s appear also to be a significant factor in the dramatic increase of NH4+ concentrations during the last few decades.

Atmospheric Soluble Dust Records from a Tibetan Ice Core: Possible Climate Proxies and Teleconnection With the Pacific Decadal Oscillation

In autumn 2005, a joint expedition between the University of Maine and the Institute of Tibetan Plateau Research recovered three ice cores from Guoqu Glacier (33 degrees 34'37.80 '' N, 91 degrees 10'35.3 '' E, 5720 m above sea level) on the northern side of Mt. Geladaindong, central Tibetan Plateau. Isotopes ( delta(18)O), major soluble ions (Na(+), K(+), Mg(2+), Ca(2+), Cl(-), NO(3)(-), SO(4)(2-)), and radionuclide (beta-activity) measurements from one of the cores revealed a 70-year record (1935-2005). Statistical analysis of major ion time series suggests that atmospheric soluble dust species dominate the chemical signature and that background dust levels conceal marine ion species deposition. The soluble dust time series have interspecies relations and common structure (empirical orthogonal function (EOF) 1), suggesting a similar soluble dust source or transport route. Annual and seasonal correlations between the EOF 1 time series and National Centers for Environmental Prediction/National Center for Atmospheric Research reanalysis climate variables (1948-2004) suggest that the Mt. Geladaindong ice core record provides a proxy for local and regional surface pressure. An approximately threefold decrease of soluble dust concentrations in the middle to late 1970s, accompanied by regional increases in pressure and temperature and decreases in wind velocity, coincides with the major 1976-1977 shift of the Pacific Decadal Oscillation (PDO) from a negative to a positive state. This is the first ice core evidence of a potential teleconnection between central Asian atmospheric soluble dust loading and the PDO. Analysis of temporally longer ice cores from Mt. Geladaindong may enhance understanding of the relationship between the PDO and central Asian atmospheric circulation and subsequent atmospheric soluble dust loading.

Association Between Atmospheric Circulation Patterns and Firn-Ice Core Records from the Inilchek Glacierized Area, Central Tien Shan, Asia

Glacioclimatological research in the central Tien Shan was performed in the summers of 1998 and 1999 on the South Inilchek Glacier at 5100 - 5460 m. A 14.36 m firn-ice core and snow samples were collected and used for stratigraphic, isotopic, and chemical analyses. The firn-ice core and snow records were related to snow pit measurements at an event scale and to meteorological data and synoptic indices of atmospheric circulation at annual and seasonal scales. Linear relationships between the seasonal air temperature and seasonal isotopic composition in accumulated precipitation were established. Changes in the delta(18)O air temperature relationship, in major ion concentration and in the ratios between chemical species, were used to identify different sources of moisture and investigate changes in atmospheric circulation patterns. Precipitation over the central Tien Shan is characterized by the lowest ionic content among the Tien Shan glaciers and indicates its mainly marine origin. In seasons of minimum precipitation, autumn and winter, water vapor was derived from the arid and semiarid regions in central Eurasia and contributed annual maximal solute content to snow accumulation in Tien Shan. The lowest content of major ions was observed in spring and summer layers, which represent maximum seasonal accumulation when moisture originates over the Atlantic Ocean and Mediterranean and Black Seas.

Pliocene Paleoenvironment and Antarctic Ice Sheet Behavior: Evidence from Wright Valley

Investigations in Wright Valley, adjacent to the Transantarctic Mountains in East Antarctica, shed light on the question of whether high-latitude Pliocene climate was warm enough to cause widespread deglaciation of the East Antarctic craton with a concurrent Magellanic moorland-like environment. If Pliocene age diatoms, presently in glaciogenic deposits high in the Transantarctic Mountains, had come from seaways on the East Antarctic craton, an expanding Late Pliocene ice sheet must have first eroded them from marine sediments and then deposited the diatoms at their present high-altitude locations. This hypothetical expanding glacier would have had to have come through Wright Valley. Glacial drift sediments from the central Wright Valley were mapped, sampled, analyzed, and Ar-40/Ar-39 whole rock dated. Our evidence indicates that an East Antarctic outlet glacier has not expanded through Wright Valley, and hence cannot have overridden the Dry Valleys sector of the Transantarctic Mountains, any time in the past 3.8 myr. Rather, there was only moderate Pliocene expansion of local cola-based alpine glaciers and continuous cold-desert conditions in Wright Valley. Persistence of a cold-desert paleoenvironment implies that the sector of the East Antarctic Ice Sheet adjacent to Wright Valley has remained relatively stable without melting ablation zones since at least 3.8 Ma, in Early Pliocene time. A further implication is that Antarctic Ice Sheet behavior in the Pliocene was much like that in the Quaternary, when the ice sheet consisted of a stable, terrestrial core in East Antarctica and a dynamic, marine-based appendage in West Antarctica.

Thermal Conditions at the Bed of the Laurentide Ice Sheet in Maine During Deglaciation: Implications for Esker Formation

The University of Maine Ice Sheet Model was used to study basal conditions during retreat of the Laurentide ice sheet in Maine. Within 150 km of the margin, basal melt rates average similar to 5 mm a(-1) during retreat. They decline over the next 100km, so areas of frozen bed develop in northern Maine during retreat. By integrating the melt rate over the drainage area typically subtended by an esker, we obtained a discharge at the margin of similar to 1.2 m(3) s(-1). While such a discharge could have moved the material in the Katahdin esker, it was likely too low to build the esker in the time available. Additional water from the glacier surface was required. Temperature gradients in the basal ice increase rapidly with distance from the margin. By conducting upward into the ice all of the additional viscous heat produced by any perturbation that increases the depth of flow in a flat conduit in a distributed drainage system, these gradients inhibit the formation of sharply arched conduits in which an esker can form. This may explain why eskers commonly seem to form near the margin and are typically segmented, with later segments overlapping onto earlier ones.

Calculating Basal Thermal Zones Beneath the Antarctic Ice Sheet

A procedure is presented for using a simple flowline model to calculate the fraction of the bed that is thawed beneath present-day ice sheets, and therefore for mapping thawed, frozen, melting and freezing basal thermal zones. The procedure is based on the proposition, easily demonstrated, that variations in surface slope along ice flowlines are due primarily to variations in bed topography and ice-bed coupling, where ice-bed coupling for sheet flow is represented by the basal thawed fraction. This procedure is then applied to the central flowlines of flow bands on the Antarctic ice sheet where accumulation rates, surface elevations and bed topography are mapped with sufficient accuracy, and where sheet flow rather than stream flow prevails. In East Antarctica, the usual condition is a low thawed fraction in subglacial highlands, but a high thawed fraction in subglacial basins and where ice converges on ice streams. This is consistent with a greater depression of the basal melting temperature and a slower rate of conducting basal heat to the surface where ice is thick, and greater basal frictional heat production where ice flow is fast, as expected for steady-state flow. This correlation is reduced or even reversed where steady-state flow has been disrupted recently, notably where ice-stream surges produced the Dibble and Dalton Iceberg Tongues, both of which are now stagnating. In West Antarctica, for ice draining into the Pine Island Bay polynya of the Amundsen Sea, the basal thawed fraction is consistent with a prolonged and ongoing surge of Pine Island Glacier and with a recently initiated surge of Thwaites Glacier. For ice draining into the Ross Ice Shelf, long ice streams extend nearly to the West Antarctic ice divide. Over the rugged bed topography near the ice divide, no correlation consistent with steady-state sheet flow exists between ice thickness and the basal thawed fraction. The bed is wholly thawed beneath ice streams, even where stream flow is slow. This is consistent with ongoing gravitational collapse of ice entering the Ross Sea embayment and with unstable flow in the ice streams.