Ice Elevation Near the West Antarctic Ice Sheet Divide During the Last Glaciation

Interior ice elevations of the West Antarctic Ice Sheet (WAIS) during the last glaciation, which can serve as benchmarks for ice-sheet models, are largely unconstrained. Here we report past ice elevation data from the Ohio Range, located near the WAIS divide and the onset region of the Mercer Ice Stream. Cosmogenic exposure ages of glacial erratics that record a WAIS highstand similar to 125 m above the present surface date to similar to 11.5 ka. The deglacial chronology prohibits an interior WAIS contribution to meltwater pulse 1A. Our observational data of ice elevation changes compare well with predictions of a thermomechanical ice-sheet model that incorporates very low basal shear stress downstream of the present day grounding line. We conclude that ice streams in the Ross Sea Embayment had thin, low-slope profiles during the last glaciation and interior WAIS ice elevations during this period were several hundred meters lower than previous reconstructions.

Reduction in Northward Incursions of the South Asian Monsoon Since ~1400 AD Inferred from a Mt. Everest Ice Core

A highly resolved Mt. Everest ice core reveals a decrease in marine and increase in continental air masses related to relatively high summer surface pressure over Mongolia, and reduction in northward incursions of the summer South Asian monsoon since similar to 1400 AD. Previously published proxy records from lower sites south of the Himalayas indicate strengthening of the monsoon since this time. These regional differences are consistent with a south north seesaw in convective activity in the Asian monsoon region, and reflect a southward shift in the mean summer position of the monsoon trough since similar to 1400 AD. The change in monsoonal circulation at 1400 AD is synchronous with a reduction in solar irradiance and the onset of the LIA. This demonstrates a hemispheric scale circulation reorganization at this time, and the potential for future large shifts in monsoonal circulation.

The Polar Expression of ENSO and Sea-Ice Variability as Recorded in a South Pole Ice Core

An annually dated ice core recovered from South Pole (2850 in a.s.l.) in 1995, that covers the period 1487-1992, was analyzed for the marine biogenic sulfur species methanesulfonate (MS). Empirical orthogonal function analysis is used to calibrate the high-resolution MS series with associated environmental series for the period of overlap (1973-92). Utilizing this calibration we present a similar to500 year long proxy record of the polar expression of the El Nino-Southern Oscillation (ENSO) and southeastern Pacific sea-ice extent variations. These records reveal short-term periods of increased (1800-50, 1900-40) and decreased sea-ice extent (1550-1610., 1660-1710, 1760-1800). In general, increased (decreased) sea-ice extent is associated with a higher (lower) frequency of El Nino events.

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.

Bending Shear: The Rate-Controlling Mechanism for Calving Ice Walls

Bending shear was observed to produce nearly vertical shear bands in a calving ice wall standing on dry land on Deception Island (Iat. 63.0 oS., long. 60.6 W.), and slabs calved straight downward when shear rupture occurred along these shear bands (Hughes, 1989). A formula for the calving rate was developed from the Deception Island data, and we have attempted to justify generalizing this formula to include ice walls standing along beaches or in water. These are environments in which a wave-washed groove develops along the base of the ice wall or along a water line above the base. The rate of wave erosion provides an alternative mechanism for controlling the calving rate in these environments. We have determined that the rate at which bending creep produces nearly vertical shear bands, along which shear r upture occurs, controls the calving rate in all environments. Shear rupture occurs at a calving shear stress of about I bar. Our results justify using the calving formula to compute the calving rate of ice walls in computer models of ice-sheet dynamics. This is especially important in simulating retreat of Northern Hemisphere ice sheets during the last deglaciation, when marine and lacustrine environments were common along retreating ice margins. These margins would have been ice walls standing along beaches or in water, because floating ice shelves are not expected in the ablation zone of retreating ice sheets.

Solar Forcing Recorded by Aerosol Concentrations in Coastal Antarctic Glacier Ice, McMurdo Dry Valleys

Ice-core chemistry data from Victoria Lower Glacier, Antarctica, suggest, at least for the last 50 years, a direct influence of solar activity variations on the McMurdo Dry Valleys (MDV) climate system via controls on air-mass input from two competing environments: the East Antarctic ice sheet and the Ross Sea. During periods of increased solar activity, when total solar irradiance is relatively high, the MDV climate system appears to be dominated by air masses originating from the Ross Sea, leading to higher aerosol deposition. During reduced solar activity, the Antarctic interior seems to be the dominant air-mass source, leading to lower aerosol concentration in the ice-core record. We propose that the sensitivity of the MDV to variations in solar irradiance is caused by strong albedo differences between the ice-free MDV and the ice sheet.

The 1452 or 1453 A.D. Kuwae Eruption Signal Derived from Multiple Ice Core Records: Greatest Volcanic Sulfate Event of the Past 700 Years

We combined 33 ice core records, 13 from the Northern Hemisphere and 20 from the Southern Hemisphere, to determine the timing and magnitude of the great Kuwae eruption in the mid-15th century. We extracted volcanic deposition signals by applying a high-pass loess filter to the time series and examining peaks that exceed twice the 31 year running median absolute deviation. By accounting for the dating uncertainties associated with each record, these ice core records together reveal a large volcanogenic acid deposition event during 1453 - 1457 A. D. The results suggest only one major stratospheric injection from the Kuwae eruption and confirm previous findings that the Kuwae eruption took place in late 1452 or early 1453, which may serve as a reference to evaluate and improve the dating of ice core records. The average total sulfate deposition from the Kuwae eruption was 93 kg SO4/km(2) in Antarctica and 25 kg SO4/km(2) in Greenland. The deposition in Greenland was probably underestimated since it was the average value of only two northern Greenland sites with very low accumulation rates. After taking the spatial variation into consideration, the average Kuwae deposition in Greenland was estimated to be 45 kg SO4/km(2). By applying the same technique to the other major eruptions of the past 700 years our result suggests that the Kuwae eruption was the largest stratospheric sulfate event of that period, probably surpassing the total sulfate deposition of the Tambora eruption of 1815, which produced 59 kg SO4/km(2) in Antarctica and 50 kg SO4/km(2) in Greenland.

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.

Antarctic Temperatures Over the Past Two Centuries from Ice Cores

We present a reconstruction of Antarctic mean surface temperatures over the past two centuries based on water stable isotope records from high-resolution, precisely dated ice cores. Both instrumental and reconstructed temperatures indicate large interannual to decadal scale variability, with the dominant pattern being anti-phase anomalies between the main Antarctic continent and the Antarctic Peninsula region. Comparative analysis of the instrumental Southern Hemisphere (SH) mean temperature record and the reconstruction suggests that at longer timescales, temperatures over the Antarctic continent vary in phase with the SH mean. Our reconstruction suggests that Antarctic temperatures have increased by about 0.2 degrees C since the late nineteenth century. The variability and the long-term trends are strongly modulated by the SH Annular Mode in the atmospheric circulation.

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.

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.

Ice Core Evidence for a Second Volcanic Eruption Around 1809 in the Northern Hemisphere

A volcanic signal observed in ice cores from both polar regions six years prior to Tambora is attributed to an unknown tropical eruption in 1809. Recovery of dacitic tephra from the 1809 horizon in a Yukon ice core ( Eclipse) that is chemically distinct from andesitic 1809 tephra found in Antarctic ice cores indicates a second eruption in the Northern Hemisphere at this time. Together with the similar magnitude and timing of the 1809 volcanic signal in the Arctic and Antarctic, this could suggest a large tropical eruption produced the sulfate and Antarctic tephra and a minor Northern Hemisphere eruption produced the Eclipse tephra. Nonetheless, the possibility that there were coincidental eruptions of similar magnitude in both hemispheres, rather than a single tropical eruption, should not be discounted. Correctly attributing the source of the 1809 volcanic signal has important implications for modeling the magnitude and latitudinal distribution of volcanic radiative forcing.

Controls on the Basal Water Pressure in Subglacial Channels Near the Margin of the Greenland Ice Sheet

Assuming a channelized drainage system in steady state, we investigate the influence of enhanced surface melting on the water pressure in subglacial channels, compared to that of changes in conduit geometry, ice rheology and catchment variations. The analysis is carried out for a specific part of the western Greenland ice-sheet margin between 66 degrees N and 66 degrees 30' N using new high-resolution digital elevation models of the subglacial topography and the ice-sheet surface, based on an airborne ice-penetrating radar survey in 2003 and satellite repeat-track interferometric synthetic aperture radar analysis of European Remote-sensing Satellite 1 and 2 (ERS-1/-2) imagery, respectively. The water pressure is calculated up-glacier along a likely subglacial channel at distances of 1, 5 and 9 km from the outlet at the ice margin, using a modified version of Rothlisberger's equation. Our results show that for the margin of the western Greenland ice sheet, the water pressure in subglacial channels is not sensitive to realistic variations in catchment size and mean surface water input compared to small changes in conduit geometry and ice rheology.

Polar MM5 Simulations of the Winter Climate of the Laurentide Ice Sheet at the LGM

Optimized regional climate simulations are conducted using the Polar MM5, a version of the fifth-generation Pennsylvania State University-NCAR Mesoscale Model (MM5), with a 60-km horizontal resolution domain over North America during the Last Glacial Maximum (LGM, 21 000 calendar years ago), when much of the continent was covered by the Laurentide Ice Sheet (LIS). The objective is to describe the LGM annual cycle at high spatial resolution with an emphasis on the winter atmospheric circulation. Output from a tailored NCAR Community Climate Model version 3 (CCM3) simulation of the LGM climate is used to provide the initial and lateral boundary conditions for Polar MM5. LGM boundary conditions include continental ice sheets, appropriate orbital forcing, reduced CO2 concentration, paleovegetation, modified sea surface temperatures, and lowered sea level. Polar MM5 produces a substantially different atmospheric response to the LGM boundary conditions than CCM3 and other recent GCM simulations. In particular, from November to April the upper-level flow is split around a blocking anticyclone over the LIS, with a northern branch over the Canadian Arctic and a southern branch impacting southern North America. The split flow pattern is most pronounced in January and transitions into a single, consolidated jet stream that migrates northward over the LIS during summer. Sensitivity experiments indicate that the winter split flow in Polar MM5 is primarily due to mechanical forcing by LIS, although model physics and resolution also contribute to the simulated flow configuration. Polar MM5 LGM results are generally consistent with proxy climate estimates in the western United States, Alaska, and the Canadian Arctic and may help resolve some long-standing discrepancies between proxy data and previous simulations of the LGM climate.