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.

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.

Holtite and Dumortierite from the Szklary Pegmatite, Lower Silesia, Poland

The Szklary holtite is represented by three compositional varieties: (I) Ta-bearing (up to 14.66 wt.% Ta(2)O(5)), which forms homogeneous crystals and cores within zoned crystals; (2) Ti-bearing (up to 3.82 wt.% TiO(2)), found as small domains within the core; and (3) Nb-bearing (up to 5.30 wt.% Nb(2)O(5),) forming the rims of zoned crystals. All three varieties show variable Sb+As content, reaching 19.18 wt.% Sb(2)O(3) (0.87 Sb a.p.f.u.) and 3.30 wt.% As(2)O(3) (0.22 As a.p.f.u.) in zoned Ta-bearing holtite, which constitutes the largest Sb+As content reported for the mineral. The zoning in holtite is a result of Ta-Nb fractionation in the parental pegmatite-forming melt together with contamination of the relatively thin Szklary dyke by Fe, Mg and Ti. Holtite and the As- and Sb-bearing dumortierite, which in places overgrows the youngest Nb-bearing zone, suggest the following crystallization sequence: Ta-bearing holtite -> Ti-bearing holtite -> Nb-bearing holtite -> As- and Sb-bearing, (Ta,Nb,Ti)-poor dumortierite -> As- and Sb-dominant, (Ta,Nb,Ti)-free dumortierite-like mineral (16.81 wt.% As(2)O(3) and 10.23 wt.% Sb(2)O(3)) with (As+Sb) > Si. The last phase is potentially a new mineral species, Al(6)rectangle B(Sb,As)(3)O(15). or Al(5)rectangle(2)B(Sb,As)(3)O(12)(OH)(3), belonging to the dumortierite group. The Szklary holtite shows no evidence of clustering of compositions around 'holtite I' and 'holtite II'. Instead, the substitutions of Si(4+) by Sb(3+)+As(3+) at the Si/Sb sites and of Ta(5+) by Nb(5+) or Ti(4+) at the Al(l) site suggest possible solid solutions between: (1) (Sb,As)-poor and (Sb,As)-rich holtite; (2) dumortierite and the unnamed (As+Sb)-dominant dumortierite-like mineral; and (3) Ti-bearing dumortierite and holtite, i.e. our data provide further evidence for miscibility between holtite and dumortierite, but leave open the question of defining the distinction between them. The Szklary holtite crystallized from the melt along with other primary Ta-Nb-(Ti) minerals such as columbite-(Mn), tantalite-(Mn), stibiotantalite and stibiocolumbite as the availability of Ta decreased. The origin of the parental melt can be related to anatexis in the adjacent Sowie Mountains complex, leading to widespread migmatization and metamorphic segregation in pelitic-psammitic sediments metamorphosed at similar to 390-380 Ma.

Magma Extraction from the Mantle Wedge at Convergent Margins Through Dikes: A Parametric Sensitivity Analysis

We address under what conditions a magma generated by partial melting at 100 km depth in the mantle wedge above a subduction zone can reach the crust in dikes before stalling. We also address under what conditions primitive basaltic magma (Mg # >60) can be delivered from this depth to the crust. We employ linear elastic fracture mechanics with magma solidification theory and perform a parametric sensitivity analysis. All dikes are initiated at a depth of 100 km in the thermal core of the wedge, and the Moho is fixed at 35 km depth. We consider a range of melt solidus temperatures (800-1100 degrees C), viscosities (10-100 Pa s), and densities (2400-2700 kg m(-3)). We also consider a range of host rock fracture toughness values (50-300 MPa m(1/2)) and dike lengths (2-5 km) and two thermal structures for the mantle wedge (1260 and 1400 degrees C at 100 km depth and 760 and 900 degrees C at 35 km depth). For the given parameter space, many dikes can reach the Moho in less than a few hundred hours, well within the time constraints provided by U series isotope disequilibria studies. Increasing the temperature in the mantle wedge, or increasing the dike length, allows additional dikes to propagate to the Moho. We conclude that some dikes with vertical lengths near their critical lengths and relatively high solidus temperatures will stall in the mantle before reaching the Moho, and these may be returned by corner flow to depths where they can melt under hydrous conditions. Thus, a chemical signature in arc lavas suggesting partial melting of slab basalts may be partly influenced by these recycled dikes. Alternatively, dikes with lengths well above their critical lengths can easily deliver primitive magmas to the crust, particularly if the mantle wedge is relatively hot. Dike transport remains a viable primary mechanism of magma ascent in convergent tectonic settings, but the potential for less rapid mechanisms making an important contribution increases as the mantle temperature at the Moho approaches the solidus temperature of the magma.

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.

Summer Temperature Trend Over the Past Two Millennia Using Air Content in Himalayan Ice

Two Himalayan ice cores display a factor-two decreasing trend of air content over the past two millennia, in contrast to the relatively stable values in Greenland and Antarctica ice cores over the same period. Because the air content can be related with the relative frequency and intensity of melt phenomena, its variations along the Himalayan ice cores provide an indication of summer temperature trend. Our reconstruction point toward an unprecedented warming trend in the 20th century but does not depict the usual trends associated with "Medieval Warm Period" (MWP), or "Little Ice Age" (LIA).

Oxygen Isotopic and Soluble Ionic Composition of a Shallow Firn Core, Inilchek Glacier, Central Tien Shan

Oxygen isotopic and soluble ionic measurements made on snow-pit (2 in depth) and firn-core (12.4 m depth samples recovered from the accumulation zone 5100 m) of Inilchek glacier 43degrees N, 79degrees E) provide information on recent (1992-98) climatic and environmental conditions in the central Tien Shan region of central Asia. The combined 14.4 m snow-pit/firn-core profile lies within the firn zone, arid contains only one observed melt feature (10 m temperature = - 12 degreesC), Although some post-depositional attenuation of the sub-seasonal delta(18)O record is possible, annual cycles are apparent throughout the isotope profile. We therefore use the preserved delta(18)O record to establish a depth/age scale for the core. Mean delta(18)O values for the entire core and for summer periods are consistent with delta(18)O/temperature observations, and suggest the delta(18)O record provides a means to reconstruct past changes in summer surface temperature at the site. Major-ion (Na(+), K(+), Mg(2+), Ca(2+), NH(4)(+), Cl(-), NO(3)(-), SO(4)(2-)) data from the core demonstrate the dominant influence of dust deposition on the soluble chemistry at the site, arid indicate significant interannual variability in atmospheric-dust loading during the 1900s. Anthropogenic impacts oil NH(4)(+) concentrations are observed at the site, and suggest a summer increase in atmospheric NH(4)(+) that may be related to regional agricultural (nitrogen-rich fertilizer use activities.

Evidence for a Frozen Bed, Byrd Glacier, Antarctica

Ice thickness, computed within the fjord region of Byrd Glacier on the assumptions that Byrd Glacier is in mass-balance equilibrium and that ice velocity is entirely due to basal sliding, are on average 400 m less than measured ice thicknesses along a radio-echo profile. We consider four explanations for these differences: (1) active glacier ice is separated from a zone of stagnant ice near the base of the glacier by a shear zone at depth; (2) basal melting rates are some 8 m/yr; (3) internal shear occurs with no basal sliding in much of the region above the grounding zone; or (4) internal creep and basal sliding contribute to the flow velocity in varying proportions above the grounding zone. Large gradients of surface strain rate seem to invalidate the first explanation. Computed values of basal shear stress (140 to 200 kPa) provide insufficient frictional heat to melt the ice demanded by the second explanation. Both the third and fourth explanations were examined by making simplifying assumptions that prevented a truly quantitative evaluation of their merit. Nevertheless, there is no escaping the qualitative conclusion that internal shear contributes strongly to surface velocities measured on Byrd Glacier, as is postulated in both these explanations.

On the Disintegration of Ice Shelves: The Role of Fracture

Crevasses can be ignored in studying the dynamics of most glaciers because they are only about 20 m deep, a small fraction of ice thickness. In ice shelves, however, s urface crevasses 20 m deep often reach sealevel and bottom crevasses can move upward to sea-level (Clough, 1974; Weertman, 1980). The ice shelf is fractured completely through if surface and basal crevasses meet (Barrett, 1975; Hughes, 1979). This is especially likely if surface melt water fills surface crevasses (Weertman, 1973; Pfeffer, 1982; Fastook and Schmidt, 1982). Fracture may therefore play an important role i n the disintegration of ice shelves. Two fracture criteria which can be evaluated experimentally and applied to ice shelves, are presented. Fracture is then examined for the general strain field of an ice shelf and for local strain fields caused by shear rupture alongside ice streams entering the ice shelf, fatigue rupture along ice shelf grounding lines, and buckling up-stream from ice rises. The effect of these fracture patterns on the stability of Antarctic ice shelves and the West Antarctic ice sheet is then discussed.

Seasonal Geophysical Monitoring of Biogenic Gases in a Northern Peatland: Implications for Temporal and Spatial Variability in Free Phase Gas Production Rates

A set of high resolution surface ground penetrating radar (GPR) surveys, combined with elevation rod ( to monitor surface deformation) and gas flux measurements, were used to investigate in situ biogenic gas dynamics within a northern peatland (Caribou Bog, Maine). Gas production rates were directly estimated from the time series of GPR measurements. Spatial variability in gas production was also investigated by comparing two sites with different geological and ecological attributes, showing differences and/or similarities depending on season. One site characterized by thick highly humified peat deposits (5-6 m), wooded heath vegetation and open pools showed large ebullition events during the summer season, with estimated emissions (based on an assumed range of CH(4) concentration) between 100 and 172 g CH(4) m(-2) during a single event. The other site characterized by thinner less humified peat deposits (2-3 m) and shrub vegetation showed much smaller ebullition events during the same season (between 13 and 23 g CH(4) m(-2)). A consistent period of free-phase gas (FPG) accumulation during the fall and winter, enhanced by the frozen surficial peat acting as a confining layer, was followed by a decrease in FPG after the snow/ice melt that released estimated fluxes between 100 and 200 g CH(4) m(-2) from both sites. Estimated FPG production rates during periods of biogenic gas accumulation ranged between 0.22 and 2.00 g CH(4) m(3) d(-1) and reflected strong seasonal and spatial variability associated with differences in temperature, peat soil properties, and/or depositional attributes (e. g., stratigraphy). Periods of decreased atmospheric pressure coincided with short-period increases in biogenic gas flux, including a very rapid decrease in FPG content associated with an ebullition event that released an estimated 39 and 67 g CH(4) m(-2) in less than 3.5 hours. These results provide insights into the spatial and seasonal variability in production and emission of biogenic gases from northern peatlands.

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.

Sliding of Ice Past an Obstacle at Engabreen, Norway

At Engabreen, Norway, an instrumented panel containing a decimetric obstacle was mounted flush With the bed surface beneath 210 m of ice. Simultaneous measurements of normal and shear stresses, ice velocity and temperature were obtained as dirty basal ice flowed past the obstacle. Our measurements were broadly consistent with ice thickness, flow conditions and bedrock topography near the site of the experiment. Ice speed 0.45 m above the bed was about 130 mm d(-1), much less than the surface velocity of 800 mm d(-1) Average normal stress on the panel was 1.0-1.6 MPa, smaller than the expected ice overburden pressure. Normal stress was larger and temperature was lower on the stoss side than on the lee side, in accord with flow dynamics and equilibrium thermodynamics. Annual differences in normal stresses were correlated with changes in sliding speed and ice-flow direction. These temporal variations may have been caused by changes in ice rheology associated with changes in sediment concentration, water content or both. Temperature and normal stress were generally correlated, except when clasts presumably collided with the panel. Temperature gradients in the obstacle indicated that regelation was negligible, consistent with the obstacle size. Melt rate was about 10% of the sliding speed. Despite high sliding speed, no significant ice/bed separation was observed in the lee of the obstacle. Frictional forces between sediment particles in the ice and the panel, estimated from Hallet's (1981) model, indicated that friction accounted for about 5% of the measured bed-parallel force. This value is uncertain, as friction theories are largely untested. Some of these findings agree with sliding theories, others do not.

Prismatine and Ferrohogbomite-2N2S in Granulite-Facies Fe-Oxide Lenses in the Eastern Ghats Belt at Venugopalapuram, Vizianagaram District, Andhra Pradesh, India: Do Such Lenses Have a Tourmaline-Enriched Lateritic Precursor?

Fluorine-rich prismatine, (square,Fe,Mg)(Mg,Al,Fe)(5)Al-4(Si,B,Al)(5)O-21(OH,F), with F/(OH+F) = 0.36-0.40 and hercynite are major constituents of a Fe-Al-B-rich lens in ultrahigh-temperature granulite-facies quartz-sillimanite-hypersthene-cordierite gneisses of the Eastern Ghats belt, Andhra Pradesh, India. Hemo-ilmenite. sapphirine, magnetite, biotite and sillimanite are subordinate. Lithium, Be and B are concentrated in prismatine (140 ppm Li, 170 ppm Be, and 2.8-3.0 wt.% B2O3). Another Fe-rich lens is dominantly magnetite, which encloses fine-grained zincian ferrohogbomite-2N2S, (Fe2+ Mg,Zn,Al)(6) (Al,Fe3+,Ti)(16)O-30(OH)(2), containing minor Ga2O3 (0.30-0.92 wt.%). Fe-Al-B-rich lenses with prismatine (or kornerupine) constitute a distinctive type of B-enrichment in granulite-facies rocks and have been reported from four other localities worldwide. A scenario involving a tourmaline-enriched lateritic precursor affected by dehydration melting is our preferred explanation for the origin of the Fe-Al-B-rich lenses at the five localities. Whole-rock analyses and field relationships at another of these localities, Bok se Puts, Namaqualand, South Africa, are consistent with this scenario. Under granulite-facies conditions, tourmaline would have broken down to give korner-upine-prismatine ( other borosilicates) plus a sodic melt containing H2O and B. Removal of this melt depleted the rock in Na and B, but the formation of ferromagnesian borosilicate phases in the restite prevented total loss of B.

Tertiary Minette and Melanephelinite Dikes, Wasatch Plateau, Utah - Records of Mantle Heterogeneities and Changing Tectonics

A swarm of minette and melanephelinite dikes is exposed over 2500 km2 in and near the Wasatch Plateau, central Utah, along the western margin of the Colorado Plateaus in the transition zone with the Basin and Range province. To date, 110 vertical dikes in 25 dike sets have been recognized. Strikes shift from about N80-degrees-W for 24 Ma dikes, to about N60-degrees-W for 18 Ma, to due north for 8-7 m.y. These orientations are consistent with a shift from east-west Oligocene compression associated with subduction to east-west late Miocene crustal extension. Minettes are the most common rock type; mica-rich minette and mica-bearing melanephelinite occurs in 24 Ma dikes, whereas more ordinary minette is found in 8-7 Ma dikes. One melanephelinite dike is 18 Ma. These mafic alkaline rocks are transitional to one another in modal and major element composition but have distinctive trace element patterns and isotopic compositions; they appear to have crystallized from primitive magmas. Major, trace element, and Nd-Sr isotopic data indicate that melanephelinite, which has similarities to ocean island basalt, was derived from small degree melts of mantle with a chondritic Sm/Nd ratio probably located in the asthenosphere, but it is difficult to rule out a lithospheric source. In contrast, mica-bearing rocks (mica melanephelinite and both types of minette) are more potassic and have trace element patterns with strong Nb-Ta depletions and Sr-Nd isotopic compositions caused by involvement with a component from heterogeneously enriched lithospheric mantle with long-term enrichment of Rb or light rare earth elements (REE) (epsilon Nd as low as - 15 in minette). Light REE enrichment must have occurred anciently in the mid-Proterozoic when the lithosphere was formed and is not a result of Cenozoic subduction processes. After about 25 Ma, foundering of the subducting Farallon plate may have triggered upwelling of warm asthenospheric mantle to the base of the lithosphere. Melanephelinite magma may have separated from the asthenosphere and, while rising through the lithosphere, provided heat for lithospheric magma generation. Varying degrees of interaction between melanephelinite and small potassic melt fractions derived from the lithospheric mantle can explain the gradational character of the melanephelinite to minette suite.