An Introduction to mountain glaciers as climate indicators with spatial and temporal diversity

This article gives an introduction to the spatial and temporal diversity of mountain glaciers as climate indicators. Alongside some information about the present extent of mountain glaciation and available databases, some specific problems with the interpretation of mountain glacier changes are highlighted.

Early last glacial maximum in the southern Central Andes reveals northward shift of the westerlies at similar to 39 ka

The latitudinal position of the southern westerlies has been suggested to be a key parameter for the climate on Earth. According to the general notion, the southern westerlies were shifted equatorward during the global Last Glacial Maximum (LGM: ~24–18 ka), resulting in reduced deep ocean ventilation, accumulation of old dissolved carbon, and low atmospheric CO2 concentrations. In order to test this notion, we applied surface exposure dating on moraines in the southern Central Andes, where glacial mass balances are particularly sensitive to changes in precipitation, i.e. to the latitudinal position of the westerlies. Our results provide robust evidence that the maximum glaciation occurred already at ~39 ka, significantly predating the global LGM. This questions the role of the westerlies for atmospheric CO2, and it highlights our limited understanding of the forcings of atmospheric circulation.

Using Micro-Gravity Techniques to Map Alluvium Thickness and Pleistocene Location of the West Branch of the Susquehanna River Near Muncy, Pennsylvania

Laurentide glaciation during the early Pleistocene (~970 ka) dammed the southeast-flowing West Branch of the Susquehanna River (WBSR), scouring bedrock and creating 100-km-long glacial Lake Lesley near the Great Bend at Muncy, Pennsylvania (Ramage et al., 1998). Local drill logs and well data indicate that subsequent paleo-outwash floods and modern fluvial processes have deposited as much as 30 meters of alluvium in this area, but little is known about the valley fill architecture and the bedrock-alluvium interface. By gaining a greater understanding of the bedrock-alluvium interface the project will not only supplement existing depth to bedrock information, but also provide information pertinent to the evolution of the Muncy Valley landscape. This project determined if variations in the thickness of the valley fill were detectable using micro-gravity techniques to map the bedrock-alluvium interface. The gravity method was deemed appropriate due to scale of the study area (~30 km2), ease of operation by a single person, and the available geophysical equipment. A LaCoste and Romberg Gravitron unit was used to collect gravitational field readings at 49 locations over 5 transects across the Muncy Creek and Susquehanna River valleys (approximately 30 km2), with at least two gravity base stations per transect. Precise latitude, longitude and ground surface elevation at each location were measured using an OPUS corrected Trimble RTK-GPS unit. Base stations were chosen based on ease of access due to the necessity of repeat measurements. Gravity measurement locations were selected and marked to provide easy access and repeat measurements. The gravimeter was returned to a base station within every two hours and a looping procedure was used to determine drift and maximize confidence in the gravity measurements. A two-minute calibration reading at each station was used to minimize any tares in the data. The Gravitron digitally recorded finite impulse response filtered gravity measurements every 20 seconds at each station. A measurement period of 15 minutes was used for each base station occupation and a minimum of 5 minutes at all other locations. Longer or multiple measurements were utilized at some sites if drift or other externalities (i.e. train or truck traffic) were effecting readings. Average, median, standard deviation and 95% confidence interval were calculated for each station. Tidal, drift, latitude, free-air, Bouguer and terrain corrections were then applied. The results show that the gravitational field decreases as alluvium thickness increases across the axes of the Susquehanna River and Muncy Creek valleys. However, the location of the gravity low does not correspond with the present-day location of the West Branch of the Susquehanna River (WBSR), suggesting that the WBSR may have been constrained along Bald Eagle Mountain by a glacial lobe originating from the Muncy Creek Valley to the northeast. Using a 3-D inversion model, the topography of the bedrock-alluvium interface was determined over the extent of the study area using a density contrast of -0.8 g/cm3. Our results are consistent with the bedrock geometry of the area, and provide a low-cost, non-invasive and efficient method for exploring the subsurface and for supplementing existing well data.

Exposure dating of Late Glacial and pre-LGM moraines in the Cordon de Doña Rosa, Northern/Central Chile (~ 31°S)

Despite the important role of the Central Andes (15–30° S) for climate reconstruction, knowledge about the Quaternary glaciation is very limited due to the scarcity of organic material for radiocarbon dating. We applied 10Be surface exposure dating (SED) on 22 boulders from moraines in the Cordon de Doña Rosa, Northern/Central Chile (~31° S). The results show that several glacial advances in the southern Central Andes occurred during the Late Glacial between ~14.7±1.5 and 11.6±1.2 ka. A much more extensive glaciation is dated to ~32±3 ka, predating the temperature minimum of the global LGM (Last Glacial Maximum: ~20 ka). Reviewing these results in the paleoclimatic context, we conclude that the Late Glacial advances were most likely caused by an intensification of the tropical circulation and a corresponding increase in summer precipitation. High-latitude temperatures minima, e.g. the Younger Dryas (YD) and the Antarctic Cold Reversal (ACR) may have triggered individual advances, but current systematic exposure age uncertainties limit precise correlations. The absence of LGM moraines indicates that moisture advection was too limited to allow significant glacial advances at ~20 ka. The tropical circulation was less intensive despite the maximum in austral summer insolation. Winter precipitation was apparently also insufficient, although pollen and marine studies indicate a northward shift of the westerlies at that time. The dominant pre-LGM glacial advances in Northern/Central Chile at ~32 ka required lower temperatures and increased precipitation than today. We conclude that the westerlies were more intense and/or shifted equatorward, possibly due to increased snow and ice cover at higher southern latitudes coinciding with a minimum of insolation.

LGM and Late Glacial glacier advances in the Cordillera Real and Cochabamba (Bolivia) deduced from 10Be surface exposure dating

Surface exposure dating (SED) is an innovative tool already being widely applied for moraine dating and for Late Quaternary glacier and climate reconstruction. Here we present exposure ages of 28 boulders from the Cordillera Real and the Cordillera Cochabamba, Bolivia. Our results indicate that the local Last Glacial Maximum (LGM) in the Eastern Cordilleras occurred at ~22–25 ka and was thus synchronous to the global temperature minimum. We were also able to date several Late Glacial moraines to ~11–13 ka, which likely document lower temperatures and increased precipitation ("Coipasa" humid phase). Additionally, we recognize the existence of older Late Glacial moraines re-calculated to ~15 ka from published cosmogenic nuclide data. Those may coincide with the cold Heinrich 1 event in the North Atlantic region and the pronounced "Tauca" humid phase. We conclude that (i) exposure ages in the tropical Andes may have been overestimated so far due to methodological uncertainties, and (ii) although precipitation plays an important role for glacier mass balances in the tropical Andes, it becomes the dominant forcing for glaciation only in the drier and thus more precipitation-sensitive regions farther west and south.

The Geology of Lost Creek Canyon with Special Reference to Igneous Activity

This investigation was undertaken primarily as a problem in geologic mapping, coupled with a study of stratigraphy, glaciation, igneous phenomena, and structure. The area is admirably suited to a study of geology and geologic events. Because it is small in extent, the area was studied in some detail during the time which was devoted to field work. The record of igneous activity of past geological ages is re­markably well exposed, since Lost Creek Canyon was carved through the roof of a stock or batholith by the glaciers of the Pleistocene epoch.

The importance of glacier and forest change in hydrological climate-impact studies

Changes in land cover alter the water balance components of a catchment, due to strong interactions between soils, vegetation and the atmosphere. Therefore, hydrological climate impact studies should also integrate scenarios of associated land cover change. To reflect two severe climate-induced changes in land cover, we applied scenarios of glacier retreat and forest cover increase that were derived from the temperature signals of the climate scenarios used in this study. The climate scenarios were derived from ten regional climate models from the ENSEMBLES project. Their respective temperature and precipitation changes between the scenario period (2074–2095) and the control period (1984–2005) were used to run a hydrological model. The relative importance of each of the three types of scenarios (climate, glacier, forest) was assessed through an analysis of variance (ANOVA). Altogether, 15 mountainous catchments in Switzerland were analysed, exhibiting different degrees of glaciation during the control period (0–51%) and different degrees of forest cover increase under scenarios of change (12–55% of the catchment area). The results show that even an extreme change in forest cover is negligible with respect to changes in runoff, but it is crucial as soon as changes in evaporation or soil moisture are concerned. For the latter two variables, the relative impact of forest change is proportional to the magnitude of its change. For changes that concern 35% of the catchment area or more, the effect of forest change on summer evapotranspiration is equally or even more important than the climate signal. For catchments with a glaciation of 10% or more in the control period, the glacier retreat significantly determines summer and annual runoff. The most important source of uncertainty in this study, though, is the climate scenario and it is highly recommended to apply an ensemble of climate scenarios in the impact studies. The results presented here are valid for the climatic region they were tested for, i.e., a humid, mid-latitude mountainous environment. They might be different for regions where the evaporation is a major component of the water balance, for example. Nevertheless, a hydrological climate-impact study that assesses the additional impacts of forest and glacier change is new so far and provides insight into the question whether or not it is necessary to account for land cover changes as part of climate change impacts on hydrological systems.

Mo isotopic composition of the mid-Neoproterozoic ocean: an iron formation perspective

The Neoproterozoic was a major turning point in Earth's surficial history, recording several widespread glaciations, the first appearance of complex metazoan life, and a major increase in atmospheric oxygen. Marine redox proxies have resulted in many different estimates of both the timing and magnitude of the increase in free oxygen, although the consensus has been that it occurred following the Marinoan glaciation, the second globally recorded “snowball Earth” event. A critically understudied rock type of the Neoproterozoic is iron formation associated with the Sturtian (first) glaciation. Samples from the <716 Ma Rapitan iron formation were analysed for their Re concentrations and Mo isotopic composition to refine the redox history of its depositional basin. Rhenium concentrations and Re/Mo ratios are consistently low throughout the bottom and middle of the iron formation, reflecting ferruginous to oxic basinal conditions, but samples from the uppermost jasper layers of the iron formation show significantly higher Re concentrations and Re/Mo ratios, indicating that iron formation deposition was terminated by a shift towards a sulfidic water column. Similarly, the δ98Mo values are close to 0.0‰ throughout most of the iron formation, but rise to ~+0.7‰ near the top of the section. The δ98Mo from samples of ferruginous to oxic basinal conditions are the product of adsorption to hematite, indicating that the Neoproterozoic open ocean may have had a δ98Mo of ~1.8‰. Together with the now well-established lack of a positive Eu anomaly in Neoproterozoic iron formations, these results suggest that the ocean was predominantly oxygenated at 700 Ma.

Water versus ice: The competing roles of modern climate and Pleistocene glacial erosion in the Central Alps of Switzerland

Recent studies have identified relationships between landscape form, erosion and climate in regions of landscape rejuvenation, associated with increased denudation. Most of these landscapes are located in non-glaciated mountain ranges and are characterized by transient geomorphic features. The landscapes of the Swiss Alps are likewise in a transient geomorphic state as seen by multiple knickzones. In this mountain belt, the transient state has been related to erosional effects during the Late Glacial Maximum (LGM). Here, we focus on the catchment scale and categorize hillslopes based on erosional mechanisms, landscape form and landcover. We then explore relationships of these variables to precipitation and extent of LGM glaciers to disentangle modern versus palaeo controls on the modern shape of the Alpine landscape. We find that in grasslands, the downslope flux of material mainly involves unconsolidated material through hillslope creep, testifying a transport-limited erosional regime. Alternatively, strength-limited hillslopes, where erosion is driven by bedrock failure, are covered by forests and/or expose bedrock, and they display oversteepened hillslopes and channels. There, hillslope gradients and relief are more closely correlated with LGM ice occurrence than with precipitation or the erodibility of the underlying bedrock. We relate the spatial occurrence of the transport- and strength-limited process domains to the erosive effects of LGM glaciers. In particular, strength-limited, rock dominated basins are situated above the equilibrium line altitude (ELA) of the LGM, reflecting the ability of glaciers to scour the landscape beyond threshold slope conditions. In contrast, transport-limited, soil-mantled landscapes are common below the ELA. Hillslopes covered by forests occupy the elevations around the ELA and are constrained by the tree line. We conclude that the current erosional forces at work in the Central Alps are still responding to LGM glaciation, and that the modern climate has not yet impacted on the modern landscape.

Where to find 1.5 million yr old ice for the IPICS "Oldest-Ice" ice core

The recovery of a 1.5 million yr long ice core from Antarctica represents a keystone of our understanding of Quaternary climate, the progression of glaciation over this time period and the role of greenhouse gas cycles in this progression. Here we tackle the question of where such ice may still be found in the Antarctic ice sheet. We can show that such old ice is most likely to exist in the plateau area of the East Antarctic ice sheet (EAIS) without stratigraphic disturbance and should be able to be recovered after careful pre-site selection studies. Based on a simple ice and heat flow model and glaciological observations, we conclude that positions in the vicinity of major domes and saddle position on the East Antarctic Plateau will most likely have such old ice in store and represent the best study areas for dedicated reconnaissance studies in the near future. In contrast to previous ice core drill site selections, however, we strongly suggest significantly reduced ice thickness to avoid bottom melting. For example for the geothermal heat flux and accumulation conditions at Dome C, an ice thickness lower than but close to about 2500 m would be required to find 1.5 Myr old ice (i.e., more than 700 m less than at the current EPICA Dome C drill site). Within this constraint, the resolution of an Oldest-Ice record and the distance of such old ice to the bedrock should be maximized to avoid ice flow disturbances, for example, by finding locations with minimum geothermal heat flux. As the geothermal heat flux is largely unknown for the EAIS, this parameter has to be carefully determined beforehand. In addition, detailed bedrock topography and ice flow history has to be reconstructed for candidates of an Oldest-Ice ice coring site. Finally, we argue strongly for rapid access drilling before any full, deep ice coring activity commences to bring datable samples to the surface and to allow an age check of the oldest ice.

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.

Bedload composition, transport and modification in rivers of Westland, New Zealand, with implications for the distribution of alluvial pounamu (jade)

River bedload surveyed at 50 sites in Westland is dominated by Alpine Schist or Torlesse Greywacke from the Alpine Fault hanging wall, with subordinate Pounamu Ultramafics or footwall-derived Western Province rocks. Tumbling experiments found ultramafics to have the lowest attrition rates, compared with greywacke sandstone and granite (which abrade to produce silt to medium-sand), or incompetent schist (which fragments). Arahura has greater total concentrations (103–105 t/km2) and proportions (5–40%) of ultramafic bedload compared with Hokitika and Taramakau catchments (101–104 t/km2, mostly <10%), matching relative areas of mapped Pounamu Ultramafic bedrock, but enriched relative to absolute areal proportions. Western Province rocks downthrown by the Alpine Fault are under-represented in the bedload. Enriched concentrations of ultramafic bedload decrease rapidly with distance downstream from source rock outcrops, changing near prominent ice-limit moraines. Bedload evolution with transport involves both downstream fining and dilution from tributaries, in a sediment supply regime more strongly influenced by tectonics and the imprint of past glaciation. Treasured New Zealand pounamu (jade) is associated with ultramafic rocks. Chances of discovery vary between catchments, are increased near glacial moraines, and are highest near source-rock outcrops in remote mountain headwaters.

Paleoglaciated landscapes in Simen and other high-mountain areas of Ethiopia

In the present-day Ethiopia, glaciated landscapes do not exist, but paleoglaciated landscapes have been documented on a few mountain tops, which have altitudes higher than about 4,350 m asl in northern Ethiopia (Simen Mountains) and about 4,100 m asl in southern Ethiopia (Arsi and Bale Mountains). Glaciers were associated with the Late Pleistocene cold stages and reached as far down as 3,760 m asl in northern and 3,200 m asl in southern Ethiopia. Bale Mountains had the most extensive Late Pleistocene glaciation, covering over 190 km2, followed by Arsi Mountains (about 85 km2). In Simen, the Late Pleistocene glaciers covered merely 13 km2. In addition, paleo-periglacial slope deposits are found on all above-mentioned paleoglaciated mountains and in further mountain systems which did not host glaciers. This allows the reconstruction of the Late Pleistocene paleoclimate as being about 8 °C colder than at present (2014), much more dry, and probably without monsoon, at least in northern Ethiopia. Most probably in the Early Holocene, the re-emergence of monsoonal rains led to a strong erosion phase, which was followed by an extended stable phase with soil formation, building up about 70-cm-deep A-horizons (Andosol) on the paleo-periglacial slope deposits. These soils have been heavily degraded due to human-induced soil erosion up to about 3800 m asl since agriculture started several decades to millennia ago.

Treeline Fluctuations Recorded for 12,500 Years by Soil Profiles, Pollen, and Plant Macrofossils in the Central Swiss Alps

Past treelines can rarely be recorded by pollen percentages alone, but pollen concentration, pollen influx, and plant macrofossils (including stomata of conifers) are more reliable indicators. In addition, ancient forest soils above today's treeline may trace the maximum upper expansion of the forest since the last glaciation. Charcoal in such soil profiles may be radiocarbon dated. Our example from the Central Swiss Alps at the Alpe d'Essertse consists of a plant-macrofossil diagram and pollen diagrams of the pond Gouille Rion at 2343 m a.s.l. and a sequence of soil profiles from 1780 m to 2600 m a.s.l. The area around the pond was forested with LariJc decidua and Pinus cembra between 9500 and 3600 BP. After 4700 BP the forest became more open and Juniperus nana and Alnus viridis expanded (together with Picea abies in the subalpine forest). Between 1700 and 900 BP the Juniperus nana and Alnus viridis scrubs declined while meadows and pastures took over, so that the pond Gouille Rion was definitively above timber­ line. The highest Holocene treeline was at 2400 to 2450 m a.s.l. (i.e. 50 to 100 m higher than the uppermost single specimen of Pinus cembra today) between 9000 and 4700 BP, but it is not yet dated in more detail. The highest charcoal of Pinus cembra at 2380 m a.s.l. has a radiocarbon date of 6010 ± 70 BP. Around 6900 BP a strong climatic deterioration caused an opening of timberline forest. First indicators of anthropogenic influence occurred at 4700 BP, when the forest limit started to move down. The lowering of timberline after 4700 BP was probably due to combined effects of human and climatic impact.

Patterns of landscape form in the upper Rhône basin, Central Swiss Alps, predominantly show lithologic controls despite multiple glaciations and variations in rock uplift rates

The development of topography depends mainly on the interplay between uplift and erosion. These processes are controlled by various factors including climate, glaciers, lithology, seismic activity and short-term variables, such as anthropogenic impact. Many studies in orogens all over the world have shown how these controlling variables may affect the landscape's topography. In particular, it has been hypothesized that lithology exerts a dominant control on erosion rates and landscape morphology. However, clear demonstrations of this influence are rare and difficult to disentangle from the overprint of other signals such as climate or tectonics. In this study we focus on the upper Rhône Basin situated in the Central Swiss Alps in order to explore the relation between topography, possible controlling variables and lithology in particular. The Rhône Basin has been affected by spatially variable uplift, high orographically driven rainfalls and multiple glaciations. Furthermore, lithology and erodibility vary substantially within the basin. Thanks to high-resolution geological, climatic and topographic data, the Rhône Basin is a suitable laboratory to explore these complexities. Elevation, relief, slope and hypsometric data as well as river profile information from digital elevation models are used to characterize the landscape's topography of around 50 tributary basins. Additionally, uplift over different timescales, glacial inheritance, precipitation patterns and erodibility of the underlying bedrock are quantified for each basin. Results show that the chosen topographic and controlling variables vary remarkably between different tributary basins. We investigate the link between observed topographic differences and the possible controlling variables through statistical analyses. Variations of elevation, slope and relief seem to be linked to differences in long-term uplift rate, whereas elevation distributions (hypsometry) and river profile shapes may be related to glacial imprint. This confirms that the landscape of the Rhône Basin has been highly preconditioned by (past) uplift and glaciation. Linear discriminant analyses (LDAs), however, suggest a stronger link between observed topographic variations and differences in erodibility. We therefore conclude that despite evident glacial and tectonic conditioning, a lithologic control is still preserved and measurable in the landscape of the Rhône tributary basins.