Sudden Increase in Tidal Response Linked to Calving and Acceleration at a Large Greenland Outlet Glacier

Large calving events at Greenland's largest outlet glaciers are associated with glacial earthquakes and near instantaneous increases in glacier flow speed. At some glaciers and ice streams, flow is also modulated in a regular way by ocean tidal forcing at the terminus. At Helheim Glacier, analysis of geodetic data shows decimeter-level periodic position variations in response to tidal forcing. However, we also observe transient increases of more than 100% in the glacier's responsiveness to such tidal forcing following glacial-earthquake calving events. The timing and amplitude of the changes correlate strongly with the step-like increases in glacier speed and longitudinal strain rate associated with glacial earthquakes. The enhanced response to the ocean tides may be explained by a temporary disruption of the subglacial drainage system and a concomitant reduction of the friction at the ice-bedrock interface, and suggests a new means by which geodetic data may be used to infer glacier properties. Citation: de Juan, J., et al. (2010), Sudden increase in tidal response linked to calving and acceleration at a large Greenland outlet glacier, Geophys. Res. Lett., 37, L12501, doi: 10.1029/2010GL043289.

The Most Extensive Holocene Advance in the Stauning Alper, East Greenland, Occurred in the Little Ice Age

We present glacial geologic and chronologic data concerning the Holocene ice extent in the Stauning Alper of East Greenland. The retreat of ice from the late-glacial position back into the mountains was accomplished by at least 11 000 cal years B.P. The only recorded advance after this time occurred during the past few centuries (the Little Ice Age). Therefore, we postulate that the Little Ice Age event represents the maximum Holocene ice extent in this part of East Greenland.

Derived Bedrock Elevations, Strain Rates and Stresses from Measured Surface Elevations and Velocities - Jakobshavns-Isbrae, Greenland

Jakobshavns Isbrae (69 degrees 10'N, 49 degrees 5'W) drains about 6.5% of the Greenland ice sheet and is the fastest ice stream known. The Jakobshavns Isbrae basin of about 10 000 km(2) was mapped photogrammetrically from four sets of aerial photography, two taken in July 1985 and two in July 1986. Positions and elevations of several hundred natural features on the ice surface were determined for each epoch by photogrammetric block-aerial triangulation, and surface velocity vectors were computed from the positions. The two flights in 1985 yielded the best results and provided most common points (716) for velocity determinations and are therefore used in the modeling studies. The data from these irregularly spaced points were used to calculate ice elevations and velocity vectors at uniformly spaced grid paints 3 km apart by interpolation. The field of surface strain rates was then calculated from these gridded data and used to compute the field of surface deviatoric stresses, using the flow law of ice, for rectilinear coordinates, X, Y pointing eastward and northward. and curvilinear coordinates, L, T pointing longitudinally and transversely to the changing ice-flow direction. Ice-surface elevations and slopes were then used to calculate ice thicknesses and the fraction of the ice velocity due to basal sliding. Our calculated ice thicknesses are in fair agreement with an ice-thickness map based on seismic sounding and supplied to us by K. Echelmeyer. Ice thicknesses were subtracted from measured ice-surface elevations to map bed topography. Our calculation shows that basal sliding is significant only in the 10-15 km before Jakobshavns Isbrae becomes afloat in Jakobshavns IsfJord.

Evolution of Melt Pond Volume on the Surface of the Greenland Ice Sheet

The presence of surface meltwater on ice caps and ice sheets is an important glaciological and climatological characteristic. We describe an algorithm for estimating the depth and hence volume of surface melt ponds using multispectral ASTER satellite imagery. The method relies on reasonable assumptions about the albedo of the bottom surface of the ponds and the optical attenuation characteristics of the ponded meltwater. We apply the technique to sequences of satellite imagery acquired over the western margin of the Greenland Ice Sheet to derive changes in melt pond extent and volume during the period 2001 - 2004. Results show large intra- and interannual changes in ponded water volumes, and large volumes of liquid water stored in extensive slush zones.

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.

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.

Step-Wise Changes in Glacier Flow Speed Coincide With Calving and Glacial Earthquakes at Helheim Glacier, Greenland

Geodetic observations show several large, sudden increases in flow speed at Helheim Glacier, one of Greenland's largest outlet glaciers, during summer, 2007. These step-like accelerations, detected along the length of the glacier, coincide with teleseismically detected glacial earthquakes and major iceberg calving events. No coseismic offset in the position of the glacier surface is observed; instead, modest tsunamis associated with the glacial earthquakes implicate glacier calving in the seismogenic process. Our results link changes in glacier velocity directly to calving-front behavior at Greenland's largest outlet glaciers, on timescales as short as minutes to hours, and clarify the mechanism by which glacial earthquakes occur. Citation: Nettles, M., et al. (2008), Step-wise changes in glacier flow speed coincide with calving and glacial earthquakes at Helheim Glacier, Greenland, Geophys. Res. Lett., 35, L24503, doi: 10.1029/2008GL036127.

Experimental Determination of Salinity, Temperature, Growth, and Metabolic Effects on Shell Isotope Chemistry of Mytilus edulis Collected from Maine and Greenland

To study the effects of temperature, salinity, and life processes (growth rates, size, metabolic effects, and physiological/ genetic effects) on newly precipitated bivalve carbonate, we quantified shell isotopic chemistry of adult and juvenile animals of the intertidal bivalve Mytilus edulis (Blue mussel) collected alive from western Greenland and the central Gulf of Maine and cultured them under controlled conditions. Data for juvenile and adult M. edulis bivalves cultured in this study, and previously by Wanamaker et al. (2006), yielded statistically identical paleotemperature relationships. On the basis of these experiments we have developed a species-specific paleotemperature equation for the bivalve M. edulis [T degrees C = 16.28 (+/- 0.10) -4.57 (+/- 0.15) {delta(18)O(c) VPBD - delta(18)O(w) VSMOW} + 0.06 (+/- 0.06) {delta(18)O(c) VPBD - delta(18)O(w) VSMOW}(2); r(2) = 0.99; N = 323; p < 0.0001]. Compared to the Kim and O'Neil (1997) inorganic calcite equation, M. edulis deposits its shell in isotope equilibrium (delta(18)O(calcite)) with ambient water. Carbon isotopes (delta(13)C(calcite)) from sampled shells were substantially more negative than predicted values, indicating an uptake of metabolic carbon into shell carbonate, and delta(13)C(calcite) disequilibrium increased with increasing salinity. Sampled shells of M. edulis showed no significant trends in delta(18)O(calcite) based on size, cultured growth rates, or geographic collection location, suggesting that vital effects do not affect delta(18)O(calcite) in M. edulis. The broad modern and paleogeographic distribution of this bivalve, its abundance during the Holocene, and the lack of an intraspecies physiologic isotope effect demonstrated here make it an ideal nearshore paleoceanographic proxy throughout much of the North Atlantic Ocean.

Ice Dynamics Preceding Catastrophic Disintegration of the Floating Part of Jakobshavn Isbrie, Greenland

The floating terminal of Jakobshavn Isbr ae, the fastest Greenland ice stream, has disintegrated since 2002, resulting in a doubling of ice velocity and rapidly lowering inland ice elevations. Conditions prior to disintegration were modeled using control theory in a plane-stress solution, and the Missoula model of ice-shelf flow. Both approaches pointed to a mechanism that inhibits ice flow and that is not captured by either approach. Jamming of flow, an inherent property of granular materials passing through a constriction (Jakobshavn Isfjord), is postulated as the mechanism. Rapid disintegration of heavily crevassed floating ice accompanies break-up of the ice jam.

Boron-Bearing Kornerupine from Fiskenaesset, West Greenland: A Reexamination of Specimens from the Type Locality

In 1884, Lorenzen proposed the formula MgAI2SiO6 for his new mineral kornerupine from Fiskenæsset and did not suspect it to contain boron. Lacroix and de Gramont (1919) reported boron in Fiskenæsset kornerupine, while Herd (1973) found none. New analyses (ion microprobe mass analyser and spectrophotometric) of kornerupine in three specimens from the type locality, including the specimens analysed by Lorenzen and Herd, indicate the presence of boron in all three, in amounts ranging from 0.50 to 1.44 wt.% B203, e.g. (Li0.04 Na0.01 Ca0.01) (Mg3.49 Mn0.01 Fe0.17 Ti0.01 Al5.64)Σ9.30 (Si3.67 Al1.02 B0.31)Σ5 O21 (OH0.99 F0.01) for Lorenzen's specimen. Textures and chemical compositions suggest that kornerupine crystallized in equilibrium in the following assemblages, all with anorthite (An 92-95) and phlogopite (XFe = atomic Fe/(Fe + Mg) = 0.028-0.035): (1) kornerupine (0.045)-gedrite (0.067); (2) kornerupine (0.038-0.050)-sapphirine (0.032-0.035); and (3) kornerupine (0.050)-hornblende. Fluorine contents of kornerupine range from 0.01 to 0.06%, of phlogopite, from 0.09 to 0.10%. In the first assemblage, sapphirine (0.040) and corundum are enclosed in radiating bundles of kornerupine; additionally sapphirine, corundum, and/or gedrite occur with chlorite and pinite (cordierite?) as breakdown products of kornerupine. Kornerupine may have formed by reactions such as: gedrite + sapphirine + corundum + B203 (in solution) + H20 = kornerupine + anorthite + Na-phlogopite under conditions of the granulite facies. Boron for kornerupine formation was most likely remobilized by hydrous fluids from metasedimentary rocks occurring along the upper contact of the Fiskenæsset gabbro-anorthosite complex with amphibolite.

Tidal Flexure of Jakobshavns Glacier, West Greenland

Jakobshavns Glacier, a floating outlet glacier on the West Greenland coast, was surveyed during July 1976. The vertical displacements of targets along two profiles perpendicular to the fjord wall bounding the north margin of the glacier were analyzed to determine the effect of flexure caused by tidal oscillations within the fjord. An analysis based on the assumption that vertical displacements of the glacier reflected pure elastic bending yielded the conclusion that the effective thickness of the ice (i.e., the thickness which remained unaffected by surface and basal cracking and which behaved as a continuum) was ∼160 m 2.6 km upglacier from the calving front and ∼110 m 0.6 km from the calving front. An analysis based on the more realistic assumption that observed bending reflected elastic and viscoplastic deformation yielded the conclusion that the average effective thickness of the ice was 316 ± 74 m (∼40% of the estimated 800-m total thickness) 2.6 km from the calving front and 160 ± 48 m (∼21% of the estimated 750-m total) 0.6 km from the calving front. A constitutive relationship appropriate for hard glide during flexure was used.

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.

Rapid Volume Loss from Two East Greenland Outlet Glaciers Quantified Using Repeat Stereo Satellite Imagery

The coastal portions of Kangerdlugssuaq and Helheim glaciers in southeast Greenland lost at least 51 +/- 8 km(-3) yr(-1) of ice between 2001-2006 due to thinning and retreat, according to an analysis of sequential digital elevation models (DEMs) derived from stereo ASTER satellite imagery. The dominant contribution to this ice loss was dynamic thinning caused by the acceleration in flow of both glaciers. Peak rates of change, including thinning rates of similar to 90 m yr(-1), coincided with the rapid increases in flow speed. Extrapolation of the measured data to the ice divides yields an estimated combined catchment volume loss of similar to 122 +/- 30 km(-3) yr(-1), which accounts for half the total mass loss from the ice sheet reported in recent studies. These catchment-wide volume losses contributed similar to 0.31 +/- 0.07 mm yr(-1) to global sea level rise over the 5-year observation period with the coastal regions alone contributing at least 0.1 +/- 0.02 mm yr(-1).

Relating Crevassing to Non-Linear Strain in the Floating Part of Jakobshavn Isbrae, West Greenland

Jakobshavn Isbrae is a major ice stream that drains the west-central Greenland ice sheet and becomes afloat in Jakobshavn Isfiord (69degreesN, 49degreesW), where it has maintained the world's fastest-known sustained velocity and calving rate (7 km a(-1)) for at least four decades. The floating portion is approximately 12 km long and 6 km wide. Surface elevations and motion vectors were determined photogrammetrically for about 500 crevasses on the floating ice, and adjacent grounded ice, using aerial photographs obtained 2 weeks apart in July 1985. Surface strain rates were computed from a mesh of 399 quadrilateral elements having velocity measurements at each corner. It is shown that heavy crevassing of floating ice invalidates the assumptions of linear strain theory that (i) surface strain in the floating ice is homogeneous in both space and time, (ii) the squares and products of strain components are nil, and (iii) first- and second-order rotation components are small compared to strain components. Therefore, strain rates and rotation rates were also computed using non-linear strain theory. The percentage difference between computed linear and non-linear second invariants of strain rate per element were greatest (mostly in the range 40-70%) where crevassing is greatest. Isopleths of strain rate parallel and transverse to flow and elevation isopleths relate crevassing to known and inferred pinning points.