Identification of glacial meltwater runoff in a karstic environment and its implication for present and future water availability

Glaciers all over the world are expected to continue to retreat due to the global warming throughout the 21st century. Consequently, future seasonal water availability might become scarce once glacier areas have declined below a certain threshold affecting future water management strategies. Particular attention should be paid to glaciers located in a karstic environment, as parts of the meltwater can be drained by underlying karst systems, making it difficult to assess water availability. In this study tracer experiments, karst modeling and glacier melt modeling are combined in order to identify flow paths in a high alpine, glacierized, karstic environment (Glacier de la Plaine Morte, Switzerland) and to investigate current and predict future downstream water availability. Flow paths through the karst underground were determined with natural and fluorescent tracers. Subsequently, geologic information and the findings from tracer experiments were assembled in a karst model. Finally, glacier melt projections driven with a climate scenario were performed to discuss future water availability in the area surrounding the glacier. The results suggest that during late summer glacier meltwater is rapidly drained through well-developed channels at the glacier bottom to the north of the glacier, while during low flow season meltwater enters into the karst and is drained to the south. Climate change projections with the glacier melt model reveal that by the end of the century glacier melt will be significantly reduced in the summer, jeopardizing water availability in glacier-fed karst springs.

Palaeoceanographic changes in the North Atlantic during the Mid-Pleistocene Transition (MIS 31-19) as inferred from planktonic foraminiferal and calcium carbonate records

Marine sediments from the Integrated Ocean Drilling Project (IODP) Site U1314 (56.36°N, 27.88°W), in the subpolar North Atlantic, were studied for their planktonic foraminifera, calcium carbonate content, and Neogloboqudrina pachyderma sinistral (sin.) δ13C records in order to reconstruct surface and intermediate conditions in this region during the Mid-Pleistocene Transition (MPT). Variations in the palaeoceanography and regional dynamics of the Arctic Front were estimated by comparing CaCO3 content, planktonic foraminiferal species abundances, carbon isotopes and ice-rafted debris (IRD) data from Site U1314 with published data from other North Atlantic sites. Site U1314 exhibited high abundances of the polar planktonic foraminifera N. pachyderma sin. and low CaCO3 content until Marine Isotope Stage (MIS) 26, indicating a relatively southeastward position of the Arctic Front (AF) and penetration of colder and low-salinity surface arctic water-masses. Changing conditions after MIS 25, with oscillations in the position of the AF, caused an increase in the northward export of the warmer North Atlantic Current (NAC), indicated by greater abundances of non-polar planktonic foraminifera and higher CaCO3. The N. pachyderma sin. δ13C data indicate good ventilation of the upper part of the intermediate water layer in the eastern North Atlantic during both glacial and interglacial stages, except during Terminations 24/23, 22/21 and 20/1. In addition, for N. pachyderma (sin.) we distinguished two morphotypes: non-encrusted and heavily encrusted test. Results indicate that increases in the encrusted morphotype and lower planktonic foraminiferal diversity are related to the intensification of glacial conditions (lower sea-surface temperatures, sea-ice formation) during MIS 22 and 20.

Biogeochemical variability during the past 3.6 million years recorded by FTIR spectroscopy in the sediment record of Lake El'gygytgyn, Far East Russian Arctic

Abstract. A number of studies have shown that Fourier transform infrared spectroscopy (FTIRS) can be applied to quantitatively assess lacustrine sediment constituents. In this study, we developed calibration models based on FTIRS for the quantitative determination of biogenic silica (BSi; n = 420; gradient: 0.9–56.5 %), total organic carbon (TOC; n = 309; gradient: 0–2.9 %), and total inorganic carbon (TIC; n = 152; gradient: 0–0.4 %) in a 318 m-long sediment record with a basal age of 3.6 million years from Lake El’gygytgyn, Far East Russian Arctic. The developed partial least squares (PLS) regression models yield high cross-validated (CV) R2 CV = 0.86–0.91 and low root mean square error of crossvalidation (RMSECV) (3.1–7.0% of the gradient for the different properties). By applying these models to 6771 samples from the entire sediment record, we obtained detailed insight into bioproductivity variations in Lake El’gygytgyn throughout the middle to late Pliocene and Quaternary. High accumulation rates of BSi indicate a productivity maximum during the middle Pliocene (3.6–3.3 Ma), followed by gradually decreasing rates during the late Pliocene and Quaternary. The average BSi accumulation during the middle Pliocene was �3 times higher than maximum accumulation rates during the past 1.5 million years. The indicated progressive deterioration of environmental and climatic conditions in the Siberian Arctic starting at ca. 3.3 Ma is consistent with the first occurrence of glacial periods and the finally complete establishment of glacial–interglacial cycles during the Quaternary.

Climatic variability during the last interglacial inferred from geochemical proxies in the Lake El'gygytgyn sediment record

The Last Interglacial Period (LIP) is often regarded as a good analogue for potential climatic conditions under predicted global warming scenarios. Despite this, there is still debate over the nature, duration and frequency of climatic changes during this period. One particularly contentious issue has been the apparent evidence of climatic instability identified in many marine cores but seemingly lacking from many terrestrial archives, especially within the Arctic, a key region for global climate change research. In this paper, geochemical records from Lake El'gygytgyn, north-eastern Russia, are used to infer past climatic changes during the LIP from within the high Arctic. With a sampling resolution of ~20–~90 years, these records offer the potential for detailed, high-resolution palaeoclimate reconstruction. This study shows that the LIP commenced in central Chukotka ~129 thousand years ago (ka), with the warmest climatic conditions occurring between ~128 and 127 ka before being interrupted by a short-lived cold reversal. Mild climatic conditions then persisted until ~122 ka when a marked reduction in the sedimentation rate suggests a decrease in precipitation. A further climatic deterioration at ~118 ka marks the return to glacial conditions. This study highlights the value of incorporating several geochemical proxies when inferring past climatic conditions, thus providing the potential to identify signals related to environmental change within the catchment. We also demonstrate the importance of considering how changes in sedimentation rate influence proxy records, in order to develop robust palaeoenvironmental reconstructions.

Lake Floor Morphology and Sediment Architecture of Lake Torneträsk, Northern Sweden

ABSTRACT. Here we present datasets from a hydroacoustic survey in July 2011 at Lake Torneträsk, northern Sweden. Our hydroacoustic data exhibit lake floor morphologies formed by glacial erosion and accumulation processes, insights into lacustrine sediment accumulation since the beginning of deglaciation, and information on seismic activity along the Pärvie Fault. Features of glacial scouring with a high-energy relief, steep slopes, and relative reliefs of more than 50 m are observed in the large W-basin. The remainder of the lacustrine subsurface appears to host a broad variety of well preserved formations from glacial accumulation related to the last retreat of the Fennoscandian ice sheet. Deposition of glaciolacustrine and lacustrine sediments is focused in areas situated in proximity to major inlets. Sediment accumulation in distal areas of the lake seldom exceeds 2 m or is not observable. We assume that lack of sediment deposition in the lake is a result of different factors, including low rates of erosion in the catchment, a previously high lake level leading to deposition of sediments in higher elevated paleodeltas, tributaries carrying low suspension loads as a result of sedimentation in upstream lakes, and an overall low productivity in the lake. A clear off-shore trace of the Pärvie Fault could not be detected from our hydroacoustic data. However, an absence of sediment disturbance in close proximity to the presumed fault trace implies minimal seismic activity since deposition of the glaciolacustrine and lacustrine sediments.

Environmental change within the Balkan region during the past ca. 50 ka recorded in the sediments from lakes Prespa and Ohrid

Abstract. Lakes Prespa and Ohrid, in the Balkan region, are considered to be amongst the oldest lakes in Europe. Both lakes are hydraulically connected via karst aquifers. From Lake Ohrid, several sediment cores up to 15m long have been studied over the last few years. Here, we document the first long sediment record from nearby Lake Prespa to clarify the influence of Lake Prespa on Lake Ohrid and the environmental history of the region. Radiocarbon dating and dated tephra layers provide robust age control and indicate that the 10.5m long sediment record from Lake Prespa reaches back to 48 ka. Glacial sedimentation is characterized by low organic matter content and absence of carbonates in the sediments, which indicate oligotrophic conditions in both lakes. Holocene sedimentation is characterized by particularly high carbonate content in Lake Ohrid and by particularly high organic matter content in Lake Prespa, which indicates a shift towards more mesotrophic conditions in the latter. Long-term environmental change and short-term events, such as related to the Heinrich events during the Pleistocene or the 8.2 ka cooling event during the Holocene, are well recorded in both lakes, but are only evident in certain proxies. The comparison of the sediment cores from both lakes indicates that environmental change affects particularly the trophic state of Lake Prespa due to its lower volume and water depth.

Late Quaternary palaeoenvironmental reconstruction from Lakes Ohrid and Prespa (Macedonia/Albania border) using stable isotopes

Abstract. Here we present stable isotope data from three sediment records from lakes that lie along the Macedonian- Albanian border (Lake Prespa: 1 core, and Lake Ohrid: 2 cores). The records only overlap for the last 40 kyr, although the longest record contains the MIS 5/6 transition (Lake Ohrid). The sedimentary characteristics of both lakes differ significantly between the glacial and interglacial phases. At the end of MIS 6 Lake Ohrid’s water level was low (high �18Ocalcite) and, although productivity was increasing (high calcite content), the carbon supply was mainly from inorganic catchment rock sources (high �13Ccarb). During the last interglacial, calcite and TOC production and preservation increased, progressively lower �18Ocalcite suggest increase in humidity and lake levels until around 115 ka. During ca. 80 ka to 11 ka the lake records suggest cold conditions as indicated by negligible calcite precipitation and low organic matter content. In Lake Ohrid, �13Corg are complacent; in contrast, Lake Prespa shows consistently higher �13Corg suggesting a low oxidation of 13C-depleted organic matter in agreement with a general deterioration of climate conditions during the glacial. From 15 ka to the onset of the Holocene, calcite and TOC begin to increase, suggesting lake levels were probably low (high �18Ocalcite). In the Holocene (11 ka to present) enhanced productivity is manifested by high calcite and organic matter content. All three cores show an early Holocene characterised by low �18Ocalcite, apart from the very early Holocene phase in Prespa where the lowest �18Ocalcite occurs at ca. 7.5 ka, suggesting a phase of higher lake level only in (the more sensitive) Lake Prespa. From 6 ka, �18Ocalcite suggest progressive aridification, in agreement with many other records in the Mediterranean, although the uppermost sediments in one core records low �18Ocalcite which we interpret as a result of human activity. Overall, the isotope data present here confirm that these two big lakes have captured the large scale, low frequency palaeoclimate variation that is seen in Mediterranean lakes, although in detail there is much palaeoclimate information that could be gained, especially small scale, high frequency differences between this region and the Mediterranean.

Denudation rates of small transient catchments controlled by former glaciation: The Hörnli nunatak in the northeastern Swiss Alpine Foreland

Extensive glaciers repeatedly occupied the northern Alpine Foreland during the Pleistocene and left a strongly glacially overprinted low slope landscape. Only few islands appeared as nunataks standing above the surface of the large piedmont glacier lobes. These nunatak areas kept their original shape, manifested in steep catchments with mean slopes up to 33 . Even though not glaciated, these catchments where significantly affected by base-level changes occurring as a consequence of phases of glacier advances and retreats. Both domains, the glacially eroded and non-eroded, are therefore prone to different mechanisms and time-scales of fluvial and colluvial re-adjustment. In this study we investigate these effects by exploring the spatial distribution and magnitude of denudation in the Hörnli region of the eastern Swiss Alpine Foreland in the present Interglacial. The area represents both domains in a relatively small area with largely uniform tectonic, lithologic and climatic conditions. The differences in Holocene andscape evolution are investigated using topographic analyses and catchment-averaged denudation rates derived from 10Be concentrations in fluvial quartz sand. We find that in formerly non-glaciated, fluvially dominated catchments close hillslope-channel coupling prevails and that these catchments yield high average denudation rates of 350 mm/ka. Glacially overprinted catchments yielded catchment-wide denudation rates an order of magnitude lower. These low denudation rates are hypothesized to be the consequence of both (i) a dominance of slow hillslope processes and (ii) admixture of high concentration, pre-LGM glacial sediment. This suggests that a) a careful field investigation must accompany the denudation rate studies and b) that the concept of area-weighted cosmogenic nuclide denudation rates must be considered in light of the predominant catchment processes.

Lake sediments as natural seismographs: A compiled record of Late Quaternary earthquakes in Central Switzerland and its implication for Alpine deformation

Central Switzerland lies tectonically in an intraplate area and recurrence rates of strong earthquakes exceed the time span covered by historic chronicles. However, many lakes are present in the area that act as natural seismographs: their continuous, datable and high-resolution sediment succession allows extension of the earthquake catalogue to pre-historic times. This study reviews and compiles available data sets and results from more than 10 years of lacustrine palaeoseismological research in lakes of northern and Central Switzerland. The concept of using lacustrine mass-movement event stratigraphy to identify palaeo-earthquakes is showcased by presenting new data and results from Lake Zurich. The Late Glacial to Holocene mass-movement units in this lake document a complex history of varying tectonic and environmental impacts. Results include sedimentary evidence of three major and three minor, simultaneously triggered basin-wide lateral slope failure events interpreted as the fingerprints of palaeoseismic activity. A refined earthquake catalogue, which includes results from previous lake studies, reveals a non-uniform temporal distribution of earthquakes in northern and Central Switzerland. A higher frequency of earthquakes in the Late Glacial and Late Holocene period documents two different phases of neotectonic activity; they are interpreted to be related to isostatic post-glacial rebound and relatively recent (re-)activation of seismogenic zones, respectively. Magnitudes and epicentre reconstructions for the largest identified earthquakes provide evidence for two possible earthquake sources: (i) a source area in the region of the Alpine or Sub-Alpine Front due to release of accumulated north-west/south-east compressional stress related to an active basal thrust beneath the Aar massif; and (ii) a source area beneath the Alpine foreland due to reactivation of deep-seated strike-slip faults. Such activity has been repeatedly observed instrumentally, for example, during the most recent magnitude 4.2 and 3.5 earthquakes of February 2012, near Zug. The combined lacustrine record from northern and Central Switzerland indicates that at least one of these potential sources has been capable of producing magnitude 6.2 to 6.7 events in the past.

A TEX 86 lake record suggests simultaneous shifts in temperature in Central Europe and Greenland during the last deglaciation

High-resolution quantitative temperature records from continents covering glacial to interglacial transitions are scarce but important for understanding the climate system. We present the first decadal resolution record of continental temperatures in Central Europe during the last deglaciation (similar to 14,60010,600cal. yrB.P.) based on the organic geochemical palaeothermometer TEX86. The TEX86-inferred temperature record from Lake Lucerne (Vierwaldstattersee, Switzerland) reveals typical oscillations during the Late Glacial Interstadial, followed by an abrupt cooling of 2 degrees C at the onset of Younger Dryas and a rapid warming of 4 degrees C at the onset of the Holocene, within less than 350years. The remarkable resemblance with the Greenland and regional stable oxygen isotope records suggests that temperature changes in continental Europe were dominated by large-scale reorganizations in the northern hemispheric climate system.

Deciphering the driving forces of erosion rates on millennial to million-year timescales in glacially impacted landscapes: An example from the Western Alps

In many regions, tectonic uplift is the main driver of erosion over million-year (Myr) timescales, but climate changes can markedly affect the link between tectonics and erosion, causing transient variations in erosion rates. Here we study the driving forces of millennial to Myr-scale erosion rates in the French Western Alps, as estimated from in situ produced cosmogenic 10Be and a newly developed approach integrating detrital and bedrock apatite fission-track thermochronology. Millennial erosion rates from 10Be analyses vary between ~0.27 and ~1.33 m/kyr, similar to rates measured in adjacent areas of the Alps. Significant positive correlations of millennial erosion rates with geomorphic measures, in particular with the LGM ice thickness, reveal a strong transient morphological and erosional perturbation caused by repeated Quaternary glaciations. The perturbation appears independent of Myr-scale uplift and erosion gradients, with the effect that millennial erosion rates exceed Myr-scale erosion rates only in the internal Alps where the latter are low (<0.4 km/Myr). These areas, moreover, exhibit channels that clearly plot above a general linear positive relation between Myr-scale erosion rates and normalized steepness index. Glacial erosion acts irrespective of rock uplift and thus not only leads to an overall increase in erosion rates but also regulates landscape morphology and erosion rates in regions with considerable spatial gradients in Myr-scale tectonic uplift. Our study demonstrates that climate change, e.g., through occurrence of major glaciations, can markedly perturb landscape morphology and related millennial erosion rate patterns, even in regions where Myr-scale erosion rates are dominantly controlled by tectonics.

Paleoglaciations in Anatolia: A Schematic Review and First Results

Anatolia is situated in the Eastern Mediterranean region between 36 – 42N and 26 – 45E. The geological records of paleoglaciations in the high terrains of Anatolia are key archives to quantify paleoclimate change in the Eastern Mediterranean area. The climate of the Eastern Mediterranean region is influenced by three main atmospheric systems: the main middle to high latitude westerlies, the mid-latitude subtropical high-pressure systems, and the monsoon climate. Glacial geological studies in Turkey have started in the late 19th century. Glacial deposits are found mainly in the eastern, northeastern and southern part of the Anatolian Peninsula. Anatolia is the fundamental element to understand the interactions between paleoenvironment, climatic variations, and development of the human societies. As the Taurus and Black Sea Mountains are sensitively situated for the paleoclimatic reconstructions, a chronostratigraphic framework on the paleoglaciation should be elaborated. The timing of the Last Glacial Maximum (LGM) in Anatolia is still unknown. Our first results from Kavron Valley (Kaçkar Mountains, NE Turkey) are encouraging for the reconstruction of paleoglaciations in Turkey and related paleoclimatological interpretations although it is presently difficult to pinpoint the classical Last Glacial Maximum – Younger Dryas – Little Ice Age moraine sequences in the field.

River loads and modern denudation of the Alps — A review

This paper presents the first comprehensive analysis of sediment and dissolved load across an entire mountain range. We investigate patterns and rates of modern denudation of the European Alps based on a compilation of data about river loads and reservoir sedimentation from 202 drainage basins that are between ca. 1 to 10,000 km2 large. The study basins cover about 50% of the total area of the Alps. Modern glaciated basins have the highest sediment yields of up to 7000 t km− 2 a− 1, which are on average 5 to 10 times higher than in non-glaciated basins. Likewise sediment yield and glacial cover are positively correlated. Instead, relief is a relatively weak predictor of sediment yield. The strong glacial impact in the correlations is due to glacier recession since the 19th century as well as due to glacial conditioning during repeated Quaternary glaciations which have produced the strong transient state of the Alpine landscape. We suggest that this is the major cause for ca. 3 fold enhanced denudation of the western compared to the eastern Alps. Chemical denudation rates are highest in the external Alps dominated by carbonate sedimentary rocks, where they make up about one third of total denudation. The high rates cannot be explained without anhydrite dissolution. We estimated that only 45% of the sediments mobilized in headwaters are exported out off the Alps, most sediments being trapped in artificial reservoirs. The total amount of sediment annually trapped within the Alps equates to 43 Mt. When corrected for sediment storage, we obtain an area-weighted mean total denudation rate for the Alps of about 0.32 mm a− 1. The pre-dam rate might be as high as 0.42 mm a− 1. In total, ca. 35 plus 23 Mt of mass are exported each year out of the Alps as solids and solutes, respectively. These rates are not enough to out pace modern rock uplift. Nevertheless, pattern of sediment yield across the Alps coincides roughly with the intensity of glacial conditioning and modern rock uplift, supporting the hypothesis of an erosion-driven uplift of the Alps.

Identification of erosional mechanisms during past glaciations based on a bedrock surface model of the central European Alps

The bedrock topography beneath the Quaternary cover provides an important archive for the identification of erosional processes during past glaciations. Here, we combined stratigraphic investigations of more than 40,000 boreholes with published data to generate a bedrock topography model for the entire plateau north of the Swiss Alps including the valleys within the mountain belt. We compared the bedrock map with data about the pattern of the erosional resistance of Alpine rocks to identify the controls of the lithologic architecture on the location of overdeepenings. We additionally used the bedrock topography map as a basis to calculate the erosional potential of the Alpine glaciers, which was related to the thickness of the LGM ice. We used these calculations to interpret how glaciers, with support by subglacial meltwater under pressure, might have shaped the bedrock topography of the Alps. We found that the erosional resistance of the bedrock lithology mainly explains where overdeepenings in the Alpine valleys and the plateau occur. In particular, in the Alpine valleys, the locations of overdeepenings largely overlap with areas where the underlying bedrock has a low erosional resistance, or where it was shattered by faults. We also found that the assignment of two end-member scenarios of erosion, related to glacial abrasion/plucking in the Alpine valleys, and dissection by subglacial meltwater in the plateau, may be adequate to explain the pattern of overdeepenings in the Alpine realm. This most likely points to the topographic controls on glacial scouring. In the Alps, the flow of LGM and previous glaciers were constrained by valley flanks, while ice flow was mostly divergent on the plateau where valley borders are absent. We suggest that these differences in landscape conditioning might have contributed to the contrasts in the formation of overdeepenings in the Alpine valleys and the plateau.