Antarctic Volcanic Flux Ratios from Law Dome Ice Cores

Explosive volcanic eruptions can inject large quantities of sulphur dioxide into the stratosphere. The aerosols that result from oxidation of the sulphur dioxide can produce significant cooling of the troposphere by reflecting or absorbing solar radiation. It is possible to obtain an estimate of the relative stratospheric sulphur aerosol concentration produced by different volcanoes by comparing sulphuric acid fluxes determined by analysis of polar ice cores. Here, we use a non-sea-salt sulphate time series derived from three well-dated Law Dome ice cores to investigate sulphuric acid flux ratios for major eruptions over the period AD 1301-1995. We use additional data from other cores to investigate systematic spatial variability in the ratios. Only for the Kuwae eruption (Law Dome ice date AD 1459.5) was the H2SO4 flux larger than that deposited by Tambora (Law Dome ice date AD 1816.7).

West Antarctic Ice Streams

Solar heat is the acknowledged driving force for climatic change. However, ice sheets are also capable of causing climatic change. This property of ice sheets derives from the facts that ice and rock are crystalline whereas the oceans and atmosphere are fluids and that ice sheets are massive enough to depress the earth's crust well below sea level. These features allow time constants for glacial flow and isostatic compensation to be much larger than those for ocean and atmospheric circulation and therefore somewhat independent of the solar variations that control this circulation. This review examines the nature of dynamic processes in ice streams that give ice sheets their degree of independent behavior and emphasizes the consequences of viscoplastic instability inherent in anisotropic polycrystalline solids such as glacial ice. Viscoplastic instability and subglacial topography are responsible for the formation of ice streams near ice sheet margins grounded below sea level. As a result the West Antarctic marine ice sheet is inherently unstable and can be rapidly carved away by calving bays which migrate up surging ice streams. Analyses of tidal flexure along floating ice stream margins, stress and velocity fields in ice streams, and ice stream boundary conditions are presented and used to interpret ERTS 1 photomosaics for West Antarctica in terms of characteristic ice sheet crevasse patterns that can be used to monitor ice stream surges and to study calving bay dynamics.

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.

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.

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.

On the Pulling Power of Ice Streams

Gravity wants to pull an ice sheet to the center of the Earth, but cannot because the Earth's crust is in the way, so ice is pushed out sideways instead. Or is it? The ice sheet "sees" nothing preventing it from spreading out except air, which is much less massive than ice. Therefore, does not ice rush forward to fill this relative vacuum; does not the relative vacuum suck ice into it, because Nature abhors a vacuum? If so, the ice sheet is not only pulled downward by gravity, it is also pulled outward by the relative vacuum. This pulling outward will be most rapid where the ice sheet encounters least resistance. The least resistance exists along the bed of ice streams, where ice-bed coupling is reduced by a basal water layer, especially if the ice stream becomes afloat and the floating part is relatively unconfined around its perimeter and unpinned to the sea floor. Ice streams are therefore fast currents of ice that develop near the margins of an ice sheet where these conditions exist. Because of these conditions, ice streams pull ice out of ice sheets and have pulling power equal to the longitudinal gravitational pulling force multiplied by the ice-stream velocity. These boundary conditions beneath and beyond ice streams can be quantified by a basal buoyancy factor that provides a life-cycle classification of ice streams into inception, growth, mature, declining and terminal stages, during which ice streams disintegrate the ice sheet. Surface profiles of ice streams are diagnostic of the stage in a life cycle and, hence, of the vitality of the ice sheet.

Stratigraphic Variation Within Polar Firn Caused by Differential Accumulation and Ice Flow: Interpretation of a 400 Mhz Short-Pulse Radar Profile from West Antarctica

We investigate causes of the stratigraphic variation revealed in a 177 km, 400 MHz short-pulse radar profile of firn from West Antarctica. The profile covers 56 m depth, and its direction was close to those of the ice flow and mean wind. The average, near-surface accumulation rates calculated from the time delays of one radar horizon consistently show minima on leeward slopes and maxima on windward slopes, confirming an earlier study based on stake observations. The stratigraphic variation includes up to 30 m depth variation in individual horizons over tens of km, fold limbs that become progressively steeper with depth, and fold-hinge loci that change direction or propagate down-ice with depth over distances far less than predicted by the ice speeds. We use an accumulation rate model to show how local rate anomalies and the effect of ice speed upon a periodic variation in accumulation rate cause these phenomena, and we reproduce two key features seen in the stratigraphic variations. We conclude that the model provides an explanation of changes in spatial stratigraphy and local measures of accumulation history given the constraints of surface topography, ice and wind velocities, and a general accumulation rate for an area.

Dating the Siple Dome (Antarctica) Ice Core By Manual and Computer Interpretation of Annual Layering

The Holocene portion of the Siple Dome (Antarctica) ice core was dated by interpreting the electrical, visual and chemical properties of the core. The data were interpreted manually and with a computer algorithm. The algorithm interpretation was adjusted to be consistent with atmospheric methane stratigraphic ties to the GISP2 (Greenland Ice Sheet Project 2) ice core, (BE)-B-10 stratigraphic ties to the dendrochronology C-14 record and the dated volcanic stratigraphy. The algorithm interpretation is more consistent and better quantified than the tedious and subjective manual interpretation.

CO(2) Diffusion in Polar Ice: Observations from Naturally Formed CO(2) Spikes in the Siple Dome (Antarctica) Ice Core

One common assumption in interpreting ice-core CO(2) records is that diffusion in the ice does not affect the concentration profile. However, this assumption remains untested because the extremely small CO(2) diffusion coefficient in ice has not been accurately determined in the laboratory. In this study we take advantage of high levels of CO(2) associated with refrozen layers in an ice core from Siple Dome, Antarctica, to study CO(2) diffusion rates. We use noble gases (Xe/Ar and Kr/Ar), electrical conductivity and Ca(2+) ion concentrations to show that substantial CO(2) diffusion may occur in ice on timescales of thousands of years. We estimate the permeation coefficient for CO(2) in ice is similar to 4 x 10(-21) mol m(-1) s(-1) Pa(-1) at -23 degrees C in the top 287 m (corresponding to 2.74 kyr). Smoothing of the CO(2) record by diffusion at this depth/age is one or two orders of magnitude smaller than the smoothing in the firn. However, simulations for depths of similar to 930-950m (similar to 60-70 kyr) indicate that smoothing of the CO(2) record by diffusion in deep ice is comparable to smoothing in the firn. Other types of diffusion (e.g. via liquid in ice grain boundaries or veins) may also be important but their influence has not been quantified.

Seasonal Deuterium Excess in a Tien Shan Ice Core: Influence of Moisture Transport and Recycling in Central Asia

Stable water isotope (delta(18)O, deltaD) data from a high elevation (5100 masl) ice core recovered from the Tien Shan Mountains, Kyrgyzstan, display a seasonal cycle in deuterium excess (d = deltaD - 8* delta(18)O) related to changes in the regional hydrologic cycle during 1994 - 2000. While there is a strong correlation (r(2) = 0.98) between delta(18)O and dD in the ice core samples, the regression slope (6.9) and mean d value (23.0) are significantly different than the global meteoric water line values. The resulting time-series ice core d profile contains distinct winter maxima and summer minima, with a yearly d amplitude of similar to 15 - 20parts per thousand. Local-scale processes that may affect d values preserved in the ice core are not consistent with the observed seasonal variability. Data from Central Asian monitoring sites in the Global Network of Isotopes in Precipitation (GNIP) have similar seasonal d changes. We suggest that regional-scale hydrological conditions, including seasonal changes in moisture source, transport, and recycling in the Caspian/ Aral Sea region, are responsible for the observed spatial and temporal d variability.

A 12,000 Year Record of Explosive Volcanism in the Siple Dome Ice Core, West Antarctica

Air mass trajectories in the Southern Hemisphere provide a mechanism for transport to and deposition of volcanic products on the Antarctic ice sheet from local volcanoes and from tropical and subtropical volcanic centers. This study extends the detailed record of Antarctic, South American, and equatorial volcanism over the last 12,000 years using continuous glaciochemical series developed from the Siple Dome A (SDMA) ice core, West Antarctica. The largest volcanic sulfate spike ( 280 mu g/L) occurs at 5881 B. C. E. Other large signals with unknown sources are observed around 325 B. C. E. ( 270 mu g/L) and 2818 B. C. E. ( 191 mu g/L). Ages of several large equatorial or Southern Hemisphere volcanic eruptions are synchronous with many sulfate peaks detected in the SDMA volcanic ice chemistry record. The microprobe "fingerprinting'' of glass shards in the SDMA core points to the following Antarctic volcanic centers as sources of tephra found in the SDMA core: Balenny Island, Pleiades, Mount Berlin, Mount Takahe, and Mount Melbourne as well as Mount Hudson and possibly Mount Burney volcanoes of South America. Identified volcanic sources provide an insight into the poorly resolved transport history of volcanic products from source volcanoes to the West Antarctic ice sheet.

A 200 Year Sulfate Record from Sixteen Antarctic Ice Cores and Associations With Southern Ocean Sea-Ice Extent

Chemistry data from 16, 50-115 m deep, sub-annually dated ice cores are used to investigate spatial and temporal concentration variability of sea-salt (ss) SO42- and excess (xs) SO42- over West Antarctica and the South Pole for the last 200 years. Low-elevation ice-core sites in western West Antarctica contain higher concentrations Of SO42- as a result of cyclogenesis over the Ross Ice Shelf and proximity to the Ross Sea Polynya. Linear correlation analysis of 15 West Antarctic ice-core SO42- time series demonstrates that at several sites concentrations Of ssSO(4)(2-) are higher when sea-ice (SIE) extent is greater, and the inverse for XSS04. Concentrations Of XSS04 from the South Pole site (East Antarctica) are associated with SIE from the Weddell region, and West Antarctic XSSO42- concentrations are associated with SIE from the Bellingshausen-Amundsen-Ross region. The only notable rise of the last 200 years in xsSO(4)(2-), around 1940, is not related to SIE fluctuations and is most likely a result of increased xsSO(4)(2-) production in the mid-low latitudes and/or an increase in transport efficiency from the mid-low latitudes to central West Antarctica. These high-resolution records show that the source types and source areas Of ssSO(4)(2-) and xsSO(4)(2-) delivered to eastern and western West Antarctica and the South Pole differ from site to site but can best be resolved using records from spatial ice-core arrays such as the International Trans-Antarctic Scientific Expedition (ITASE).

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.