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.

Buckling Rate and Overhang Development at a Calving Face

Using the finite-element we have modeled the stress field near the calving face of an idealized tidewater glacier under a variety of assumptions about submarine calving-face height, subaerial calving-face height, and ice rheology These simulations all suggest that a speed maximum should be present at the calving face near the waterline. In experiments without crevassing, the decrease in horizontal velocity above this maximum culminates in a zone of longitudinal compression at the surface somewhat Up-glacier from the face. This zone of compression appears to be a consequence of the non-linear rheology of ice. It disappears when a linear rheology is assumed. Explorations of the near-surface stress field indicate that when pervasive crevassing of the surface ice is accounted for in the simulations (by rheological softening), the zone of compressive strain rates does not develop. Variations in the pattern of horizontal velocity with glacier thickness support the contention that calving rates should increase with water depth at the calving face. In addition, the height of the subaerial calving face may have an importance that is not visible ill Current field data owing to the lack of variation in height of such faces in nature. Glaciers with lower calving faces may not have sufficient tensile stress to calve actively, while tensile stresses in simulated higher faces are sufficiently high that such faces will be unlikely to build in nature.

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.

The Jakobshanvs Effect

The Jakobshavns Effect may have been a significant factor in hastening the collapse of palaeo ice sheets with the advent of climatic warming after 18,000 years ago and may precipitate partial collapse of the present‐day Greenland and Antarctic Ice Sheets following CO2‐induced climatic warming in the decades ahead. The Jakobshavns Effect is observed today on Jakobshavns Glacier, which is located at 69°10′N on the west coast of Greenland. The Jakobshavns Effect is a group of positive feedback mechanisms which allow Jakobshavns Glacier to literally pull ice out of the Greenland Ice Sheet at a rate exceeding 7 km/a across a floating terminus 800 m thick and 6 km wide. The pulling power results from an imbalance of horizontal hydrostatic forces in ice and water columns at the grounding line of the floating terminus. Positive feedback mechanisms that sustain the rapid ice discharge rate are ubiquitous surface crevassing, high summer rates of surface melting, extending creep flow, progressive basal uncoupling, progressive lateral uncoupling, and rapid iceberg calving.

Dynamics of the last four glacial terminations recorded in Lake Van, Turkey

A well-dated suite of Lake Van climate-proxy data covering the last 360 ka documents environmental changes over 4 glacial/interglacial cycles in Eastern Anatolia, Turkey. The picture of cold and dry glacials and warm and wet interglacials emerging from pollen, organic carbon, authigenic carbonate content, elemental profiling by XRF and lithological analyses is inconsistent with classical interpretation of ox- ygen isotopic composition of carbonates pointing to a more complex pattern in Lake Van region. Detailed analysis of glacial terminations allows for the constraining of a depositional model explaining different patterns observed in all the proxies. We hypothesize that variations in relative contribution of rainfall, snowmelt and glacier meltwater recharging the basin have a very important role for all sedimentary processes in Lake Van. Lake level of glacial Lake Van, predominantly fed by snowmelt, was low, the water column was oxic, and carbonates precipitating in the epilimnion recorded the light isotopic signature of inflow. During terminations, increasing rainfall and significant supply of mountain glaciers' meltwater contributed to lake level rise. Increased rainfall enhanced density gradients in the water column, and hindered mixing leading to development of bottom-water anoxia. Carbonates precipitating during terminations show large fluctuations in their isotopic composition. Full interglacial conditions in Lake Van are characterized by high or slowly falling lake level. Rainfall and snowmelt feed the lake but due to re-established mixing, the isotopic composition of authigenic carbonates is heavier and closer to that of evaporation-influenced lake water than that of runoff representing snowmelt and atmospheric precipitation.

Ice-Core Based Assessment of Historical Anthropogenic Heavy Metal (Cd, Cu, Sb, Zn) Emissions in the Soviet Union

The development of strategies and policies aiming at the reduction of environmental exposure to air pollution requires the assessment of historical emissions. Although anthropogenic emissions from the extended territory of the Soviet Union (SU) considerably influenced concentrations of heavy metals in the Northern Hemisphere, Pb is the only metal with long-term historical emission estimates for this region available, whereas for selected other metals only single values exist. Here we present the first study assessing long-term Cd, Cu, Sb, and Zn emissions in the SU during the period 1935–1991 based on ice-core concentration records from Belukha glacier in the Siberian Altai and emission data from 12 regions in the SU for the year 1980. We show that Zn primarily emitted from the Zn production in Ust-Kamenogorsk (East Kazakhstan) dominated the SU heavy metal emission. Cd, Sb, Zn (Cu) emissions increased between 1935 and the 1970s (1980s) due to expanded non-ferrous metal production. Emissions of the four metals in the beginning of the 1990s were as low as in the 1950s, which we attribute to the economic downturn in industry, changes in technology for an increasing metal recovery from ores, the replacement of coal and oil by gas, and air pollution control.

Temperature and precipitation signal in two Alpine ice cores over the period 1961-2001

Water stable isotope ratios and net snow accumulation in ice cores are commonly interpreted as temperature or precipitation proxies. However, only in a few cases has a direct calibration with instrumental data been attempted. In this study we took advantage of the dense network of observations in the European Alpine region to rigorously test the relationship of the annual and seasonal resolved proxy data from two highly resolved ice cores with local temperature and precipitation. We focused on the time period 1961–2001 with the highest amount and quality of meteorological data and the minimal uncertainty in ice core dating (±1 year). The two ice cores were retrieved from the Fiescherhorn glacier (northern Alps, 3900 m a.s.l.), and Grenzgletscher (southern Alps, 4200 m a.s.l.). A parallel core from the Fiescherhorn glacier allowed assessing the reproducibility of the ice core proxy data. Due to the orographic barrier, the two flanks of the Alpine chain are affected by distinct patterns of precipitation. The different location of the two glaciers therefore offers a unique opportunity to test whether such a specific setting is reflected in the proxy data. On a seasonal scale a high fraction of δ18O variability was explained by the seasonal cycle of temperature (~60% for the ice cores, ~70% for the nearby stations of the Global Network of Isotopes in Precipitation – GNIP). When the seasonality is removed, the correlations decrease for all sites, indicating that factors other than temperature such as changing moisture sources and/or precipitation regimes affect the isotopic signal on this timescale. Post-depositional phenomena may additionally modify the ice core data. On an annual scale, the δ18O/temperature relationship was significant at the Fiescherhorn, whereas for Grenzgletscher this was the case only when weighting the temperature with precipitation. In both cases the fraction of interannual temperature variability explained was ~20%, comparable to the values obtained from the GNIP stations data. Consistently with previous studies, we found an altitude effect for the δ18O of −0.17‰/100 m for an extended elevation range combining data of the two ice core sites and four GNIP stations. Significant correlations between net accumulation and precipitation were observed for Grenzgletscher during the entire period of investigation, whereas for Fiescherhorn this was the case only for the less recent period (1961–1977). Local phenomena, probably related to wind, seem to partly disturb the Fiescherhorn accumulation record. Spatial correlation analysis shows the two glaciers to be influenced by different precipitation regimes, with the Grenzgletscher reflecting the characteristic precipitation regime south of the Alps and the Fiescherhorn accumulation showing a pattern more closely linked to northern Alpine stations.

Terrestrial and submarine evidence for the extent and timing of the Last Glacial Maximum and the onset of deglaciation on the maritime-Antarctic and sub-Antarctic islands

This paper is the maritime and sub–Antarctic contribution to the Scientific Committee for Antarctic Research (SCAR) Past Antarctic Ice Sheet Dynamics (PAIS) community Antarctic Ice Sheet reconstruction. The overarching aim for all sectors of Antarctica was to reconstruct the Last Glacial Maximum (LGM) ice sheet extent and thickness, and map the subsequent deglaciation in a series of 5000 year time slices. However, our review of the literature found surprisingly few high quality chronological constraints on changing glacier extents on these timescales in the maritime and sub–Antarctic sector. Therefore, in this paper we focus on an assessment of the terrestrial and offshore evidence for the LGM ice extent, establishing minimum ages for the onset of deglaciation, and separating evidence of deglaciation from LGM limits from those associated with later Holocene glacier fluctuations. Evidence included geomorphological descriptions of glacial landscapes, radiocarbon dated basal peat and lake sediment deposits, cosmogenic isotope ages of glacial features and molecular biological data. We propose a classification of the glacial history of the maritime and sub–Antarctic islands based on this assembled evidence. These include: (Type I) islands which accumulated little or no LGM ice; (Type II) islands with a limited LGM ice extent but evidence of extensive earlier continental shelf glaciations; (Type III) seamounts and volcanoes unlikely to have accumulated significant LGM ice cover; (Type IV) islands on shallow shelves with both terrestrial and submarine evidence of LGM (and/or earlier) ice expansion; (Type V) Islands north of the Antarctic Polar Front with terrestrial evidence of LGM ice expansion; and (Type VI) islands with no data. Finally, we review the climatological and geomorphological settings that separate the glaciological history of the islands within this classification scheme.

Polychlorinated Biphenyls in Glaciers. 2. Model Results of Deposition and Incorporation Processes

In previous work, Alpine glaciers have been identified as a secondary source of persistent organic pollutants (POPs). However, detailed understanding of the processes organic chemicals undergo in a glacial system was missing. Here, we present results from a chemical fate model describing deposition and incorporation of polychlorinated biphenyls (PCBs) into an Alpine glacier (Fiescherhorn, Switzerland) and an Arctic glacier (Lomonosovfonna, Norway). To understand PCB fate and dynamics, we investigate the interaction of deposition, sorption to ice and particles in the atmosphere and within the glacier, revolatilization, diffusion and degradation, and discuss the effects of these processes on the fate of individual PCB congeners. The model is able to reproduce measured absolute concentrations in the two glaciers for most PCB congeners. While the model generally predicts concentration profiles peaking in the 1970s, in the measurements, this behavior can only be seen for higher-chlorinated PCB congeners on Fiescherhorn glacier. We suspect seasonal melt processes are disturbing the concentration profiles of the lower-chlorinated PCB congeners. While a lower-chlorinated PCB congener is mainly deposited by dry deposition and almost completely revolatilized after deposition, a higher-chlorinated PCB congener is predominantly transferred to the glacier surface by wet deposition and then is incorporated into the glacier ice. The incorporated amounts of PCBs are higher on the Alpine glacier than on the Arctic glacier due to the higher precipitation rate and aerosol particle concentration on the former. Future studies should include the effects of seasonal melt processes, calculate the quantities of PCBs incorporated into the entire glacier surface, and estimate the quantity of chemicals released from glaciers to determine the importance of glaciers as a secondary source of organic chemicals to remote aquatic ecosystems.

210Pb dating of the Miaoergou ice core from the eastern Tien Shan, China

In 2005, two ice cores with lengths of 58.7 and 57.6 m respectively to bedrock were recovered from the Miaoergou flat-topped glacier (43 degrees 03 ' 19 '' N, 94 degrees 19 ' 21 '' E; 4512 m a.s.l.), eastern Tien Shan. Pb-210 dating of one of the ice cores (57.6 m) was performed, and an age of AD 1851 +/- 6 at a depth of 35.2 m w.e. was determined. For the period AD 1851-2005, a mean annual net accumulation of 229 +/- 7 mm w.e. a(-1) was calculated. At the nearby oasis city of Hami (similar to 80 km from the Miaoergou flat-topped glacier) the annual precipitation rate is 38 mm w.e. a(-1), hence glacial meltwater is a major water supply for local residents. The surface activity concentration of Pb-210(ex) was found to be similar to 400 mBq kg(-1), which is higher than observed at other continental sites such as Belukha, Russia, and Tsambagarav, Mongolia, which have surface activity concentrations of 280 mBq kg(-1). The Pb-210 dating agrees well with the chronological sequence deduced from the annual-layer counting resulting from the seasonalities of delta O-18 and trace metals for the period AD 1953-2005, and beta-activity horizons resulting from atmospheric nuclear testing during the period AD 1962-63. We conclude that Pb-210 analysis is a suitable method for obtaining a continuous dating of the Miaoergou ice core for similar to 160 years, which can also be applied to other ice cores recovered from the mountains of western China.

Isotopic constraints on marine and terrestrial N2O emissions during the last deglaciation

Nitrous oxide (N2O) is an important greenhouse gas and ozone-depleting substance that has anthropogenic as well as natural marine and terrestrial sources. The tropospheric N2O concentrations have varied substantially in the past in concert with changing climate on glacial–interglacial and millennial timescales. It is not well understood, however, how N2O emissions from marine and terrestrial sources change in response to varying environmental conditions. The distinct isotopic compositions of marine and terrestrial N2O sources can help disentangle the relative changes in marine and terrestrial N2O emissions during past climate variations. Here we present N2O concentration and isotopic data for the last deglaciation, from 16,000 to 10,000 years before present, retrieved from air bubbles trapped in polar ice at Taylor Glacier, Antarctica. With the help of our data and a box model of the N2O cycle, we find a 30 per cent increase in total N2O emissions from the late glacial to the interglacial, with terrestrial and marine emissions contributing equally to the overall increase and generally evolving in parallel over the last deglaciation, even though there is no a priori connection between the drivers of the two sources. However, we find that terrestrial emissions dominated on centennial timescales, consistent with a state-of-the-art dynamic global vegetation and land surface process model that suggests that during the last deglaciation emission changes were strongly influenced by temperature and precipitation patterns over land surfaces. The results improve our understanding of the drivers of natural N2O emissions and are consistent with the idea that natural N2O emissions will probably increase in response to anthropogenic warming.

Pb pollution from leaded gasoline in South America in the context of a 2000-year metallurgical history

Exploitation of the extensive polymetallic deposits of the Andean Altiplano in South America since precolonial times has caused substantial emissions of neurotoxic lead (Pb) into the atmosphere; however, its historical significance compared to recent Pb pollution from leaded gasoline is not yet resolved. We present a comprehensive Pb emission history for the last two millennia for South America, based on a continuous, high-resolution, ice core record from Illimani glacier. Illimani is the highest mountain of the eastern Bolivian Andes and is located at the northeastern margin of the Andean Altiplano. The ice core Pb deposition history revealed enhanced Pb enrichment factors (EFs) due to metallurgical processing for silver production during periods of the Tiwanaku/Wari culture (AD 450–950), the Inca empires (AD 1450–1532), colonial times (AD 1532–1900), and tin production at the beginning of the 20th century. After the 1960s, Pb EFs increased by a factor of 3 compared to the emission level from metal production, which we attribute to gasoline-related Pb emissions. Our results show that anthropogenic Pb pollution levels from road traffic in South America exceed those of any historical metallurgy in the last two millennia, even in regions with exceptional high local metallurgical activity.

The deglaciation history of the Simplon region (southern Swiss Alps) constrained by 10Be exposure dating of ice-molded bedrock surfaces

The deglaciation history of the Swiss Alps after the Last Glacial Maximum involved the decay of several ice domes and the subsequent disintegration of valley glaciers at high altitude. Here we use bedrock exposure dating to reconstruct the temporal and spatial pattern of ice retreat at the Simplon Pass (altitude: ∼2000 m) located 40 km southwest of the ‘Rhône ice dome’. Eleven 10Be exposure ages from glacially polished quartz veins and ice-molded bedrock surfaces cluster tightly between 13.5 ± 0.6 ka and 15.4 ± 0.6 ka (internal errors) indicating that the Simplon Pass depression became ice-free at 14.1 ± 0.4 ka (external error of mean age). This age constraint is interpreted to record the melting of the high valley glaciers in the Simplon Pass region during the warm Bølling–Allerød interstadial shortly after the Oldest Dryas stadial. Two bedrock samples collected a few hundred meters above the pass depression yield older 10Be ages of 17.8 ± 0.6 ka and 18.0 ± 0.6 ka. These ages likely reflect the initial downwasting of the Rhône ice dome and the termination of the ice transfluence from the ice dome across the Simplon Pass toward the southern foreland. There, the retreat of the piedmont glacier in Val d’Ossola was roughly synchronous with the decay of the Rhône ice dome in the interior of the mountain belt, as shown by 10Be ages of 17.7 ± 0.9 ka and 16.1 ± 0.6 ka for a whaleback at ∼500 m elevation near Montecrestese in northern Italy. In combination with well-dated paleoclimate records derived from lake sediments, our new age data suggest that during the deglaciation of the European Alps the decay of ice domes was approximately synchronous with the retreat of piedmont glaciers in the foreland and was followed by the melting of high-altitude valley glaciers after the transition from the Oldest Dryas to the Bølling–Allerød, when mean annual temperatures rose rapidly by ∼3 °C.