(Supplement 3) Radar sounding observations including bedrock and surface elevation and ice thickness in the ASB region, Antarctica

The first Cenozoic ice sheets initiated in Antarctica from the Gamburtsev Subglacial Mountains and other highlands as a result of rapid global cooling ~34 million years ago. In the subsequent 20 million years, at a time of declining atmospheric carbon dioxide concentrations and an evolving Antarctic circumpolar current, sedimentary sequence interpretation and numerical modelling suggest that cyclical periods of ice-sheet expansion to the continental margin, followed by retreat to the subglacial highlands, occurred up to thirty times. These fluctuations were paced by orbital changes and were a major influence on global sea levels. Ice-sheet models show that the nature of such oscillations is critically dependent on the pattern and extent of Antarctic topographic lowlands. Here we show that the basal topography of the Aurora Subglacial Basin of East Antarctica, at present overlain by 2-4.5 km of ice, is characterized by a series of well-defined topographic channels within a mountain block landscape. The identification of this fjord landscape, based on new data from ice-penetrating radar, provides an improved under¬standing of the topography of the Aurora Subglacial Basin and its surroundings, and reveals a complex surface sculpted by a succession of ice-sheet configurations substantially different from today's. At different stages during its fluctuations, the edge of the East Antarctic Ice Sheet lay pinned along the margins of the Aurora Subglacial Basin, the upland boundaries of which are currently above sea level and the deepest parts of which are more than 1 km below sea level. Although the timing of the channel incision remains uncertain, our results suggest that the fjord landscape was carved by at least two ice- flow regimes of different scales and directions, each of which would have over-deepened existing topographic depressions, reversing valley floor slopes.

(Supplement 1) Selected chemical composition of Icelandic river waters and bedrock age

In this study we map the spatial distribution of selected dissolved constituents in Icelandic river waters using GIS methods to study and interpret the connection between river chemistry, bedrock, hydrology, vegetation and aquatic ecology. Five parameters were selected: alkalinity, SiO2, Mo, F and the dissolved inorganic nitrogen and dissolved inorganic phosphorus mole ratio (DIN/DIP). The highest concentrations were found in rivers draining young rocks within the volcanic rift zone and especially those draining active central volcanoes. However, several catchments on the margins of the rift zone also had high values for these parameters, due to geothermal influence or wetlands within their catchment area. The DIN/DIP mole ratio was higher than 16 in rivers draining old rocks, but lowest in rivers within the volcanic rift zone. Thus primary production in the rivers is limited by fixed dissolved nitrogen within the rift zone, but dissolved phosphorus in the old Tertiary catchments. Nitrogen fixation within the rift zone can be enhanced by high dissolved molybdenum concentrations in the vicinity of volcanoes. The river catchments in this study were subdivided into several hydrological categories. Importantly, the variation in the hydrology of the catchments cannot alone explain the variation in dissolved constituents. The presence or absence of central volcanoes, young reactive rocks, geothermal systems and wetlands is important for the chemistry of the river waters. We used too many categories within several of the river catchments to be able to determine a statistically significant connection between the chem¬istry of the river waters and the hydrological categories. More data are needed from rivers draining one single hydrological category. The spatial dissolved constituent distribution clearly revealed the difference between the two extremes, the young rocks of the volcanic rift zone and the old Tertiary terrain.

Hydrogeological investigations in the upper Ouémé catchment in Benin, West Africa

The scope of this PhD thesis was the hydrogeological conceptualisation of the Upper Ouémé river catchment in Benin. The study area exceeds 14,500 km**2 and is underlain by a crystalline basement. At this setting the typical sequence of aquifers - a regolith aquifer at the top and a fractured bedrock aquifer at the bottom - is encountered, which is found in basement areas all over Africa and elsewhere in the world. The chosen regional approach revealed important information about the hydrochemistry and hydrogeology of this catchment. Based on the regional conceptual model a numerical groundwater flow model was designed. The numerical model was used to estimate the impact of climate change on the regional groundwater resources. This study was realised within the framework of the German interdisciplinary research project IMPETUS (English translation: "Integrated approach to the efficient management of scarce water resources in West Africa"), which is jointly managed by the German universities of Bonn and Cologne. Since the year 2000 the Upper Ouémé catchment was the principal target for investigations into the relevant processes of the regional water cycle. A first study from 2000 to 2003 (Fass, 2004, http://nbn-resolving.de/urn:nbn:de:hbz:5n-03849) focused on the hydrogeology of a small local catchment (~30 km**2). In the course of this thesis five field campaigns were underdone from the year 2004 to 2006. In the beginning of 2004 a groundwater monitoring net was installed based on 12 automatic data loggers. Manual piezometric measurements and the sampling of groundwater and surface water were realised for each campaign throughout the whole study area. Water samples were analysed for major ions, for a choice of heavy metals and for their composition by deuterium, oxygen-18 and tritium. The numerical model was performed with FEFLOW. The hydraulic and hydrochemical characteristics were described for the regolith aquifer and the bedrock aquifer. The regolith aquifer plays the role of the groundwater stock with low conductivity while the fractures of the bedrock may conduct water relatively fast towards extraction points. Flow in fractures of the bedrock depends on the connectivity of the fracture network which might be of local to subregional importance. Stable isotopes in combination with hydrochemistry proved that recharge occurs on catchment scale and exclusively by precipitation. Influx of groundwater from distant areas along dominant structures like the Kandi fault or from the Atacora mountain chain is excluded. The analysis of tritium in groundwater from different depths revealed the interesting fact of the strongly rising groundwater ages. Bedrock groundwater may possibly be much older than 50 years. Equilibrium phases of the silicate weathering products kaolinite and montmorillonite showed that the deeper part of the regolith aquifer and the bedrock aquifer feature either stagnant or less mobile groundwater while the shallow aquifer level is influenced by seasonal groundwater table fluctuations. The hydrochemical data characterised this zone by the progressive change of the hydrochemical facies of recently infiltrated rainwater on its flow path into deeper parts of the aquifers. Surprisingly it was found out that seasonal influences on groundwater hydrochemistry are minor, mainly because they affect only the groundwater levels close to the surface. The transfer of the hydrogeological features of the Upper Ouémé catchment into a regional numerical model demanded a strong simplification. Groundwater tables are a reprint of the general surface morphology. Pumping or other types of groundwater extraction would have only very local impact on the available groundwater resources. It was possible to integrate IMPETUS scenario data into the groundwater model. As a result it was shown that the impact of climate change on the groundwater resources until the year 2025 under the given conditions will be negligible due to the little share of precipitation needed for recharge and the low water needs for domestic use. Reason for concern is the groundwater quality on water points in the vicinity of settlements because of contamination by human activities as shown for the village of Dogué. Nitrate concentrations achieved in many places already alerting levels. Health risks from fluoride or heavy metals were excluded for the Upper Ouémé area.

Alpine permafrost index map

The objective of this study is the production of an Alpine Permafrost Index Map (APIM) covering the entire European Alps. A unified statistical model that is based on Alpine-wide permafrost observations is used for debris and bedrock surfaces across the entire Alps. The explanatory variables of the model are mean annual air temperatures, potential incoming solar radiation and precipitation. Offset terms were applied to make model predictions for topographic and geomorphic conditions that differ from the terrain features used for model fitting. These offsets are based on literature review and involve some degree of subjective choice during model building. The assessment of the APIM is challenging because limited independent test data are available for comparison and these observations represent point information in a spatially highly variable topography. The APIM provides an index that describes the spatial distribution of permafrost and comes together with an interpretation key that helps to assess map uncertainties and to relate map contents to their actual expression in terrain. The map can be used as a first resource to estimate permafrost conditions at any given location in the European Alps in a variety of contexts such as research and spatial planning. Results show that Switzerland likely is the country with the largest permafrost area in the Alps, followed by Italy, Austria, France and Germany. Slovenia and Liechtenstein may have marginal permafrost areas. In all countries the permafrost area is expected to be larger than the glacier-covered area.

Trace element concentrations and Li isotope composition in DSDP and ODP sediments from island arcs

We have studied the sedimentary and basaltic inputs of lithium to subduction zones. Various sediments from DSDP and ODP drill cores in front of the Mariana, South Sandwich, Banda, East Sunda and Lesser Antilles island arcs have been analysed and show highly variable Li contents and d7Li values. The sediment piles in front of the Mariana and South Sandwich arcs largely consist of pelagic sediments (clays and oozes). The pelagic clays have high Li contents (up to 57.3 ppm) and Li isotope compositions ranging from +1.3? to +4.1?. The oozes have lower Li contents (7.3-16 ppm) with d7Li values of the diatom oozes from the South Sandwich lower (+2.8? to +3.2?) than those of the radiolarian oozes from the Mariana arc (+8.1? to +14.5?). Mariana sediment also contains a significant portion of volcanogenic material, which is characterised by a moderate Li content (14 ppm) and a relatively heavy isotope composition (+6.4?). Sediments from the Banda and Lesser Antilles contain considerable amounts of continental detritus, and have high Li contents (up to 74.3 ppm) and low d7Li values (around 0?), caused by weathering of continental bedrock. East Sunda sediments largely consist of calcareous oozes. These carbonate sediments display intermediate to high Li contents (2.4-41.9 ppm) and highly variable d7Li values (-1.6? to +12.8?). Basaltic oceanic crust samples from worldwide DSDP and ODP drill cores are characterised by enrichment of Li compared to fresh MORB (6.6-33.1 vs. 3.6-7.5 ppm, respectively), and show a large range in Li isotope compositions (+1.7? to +11.8?). The elemental and isotopic enrichment of Li in altered basalts is due to the uptake of isotopically heavy seawater Li during weathering. However, old oceanic crust samples from Sites 417/418 exhibit lighter Li isotope compositions compared to young basaltic crust samples from Sites 332B and 504B. This lighter Li isotope signature in old crust is unexpected and further research is needed to explore this issue.

Crustal uplift rates of 4 GPS stations near Jakobshavn Isbrae (2006-2009)

We analyze 2006-2009 data from four continuous Global Positioning System (GPS) receivers located between 5 and 150 km from the glacier Jakobshavn Isbrae, West Greenland. The GPS stations were established on bedrock to determine the vertical crustal motion due to the unloading of ice from Jakobshavn Isbrae. All stations experienced uplift, but the uplift rate at Kangia North, only 5 km from the glacier front, was about 10 mm/yr larger than the rate at Ilulissat, located only ~45 km further away. This suggests that most of the uplift is due to the unloading of the Earth's surface as Jakobshavn thins and loses mass. Our estimate of Jakobshavn's contribution to uplift rates at Kangia North and Ilulissat are 14.6 ± 1.7 mm/yr and 4.9 ± 1.1 mm/yr, respectively. The observed rates are consistent with a glacier thinning model based on repeat altimeter surveys from NASA's Airborne Topographic Mapper (ATM), which shows that Jakobshavn lost mass at an average rate of 22 ± 2 km**3/yr between 2006 and 2009. At Kangia North and Ilulissat, the predicted uplift rates computed using thinning estimates from the ATM laser altimetry are 12.1 ± 0.9 mm/yr and 3.2 ± 0.3 mm/yr, respectively. The observed rates are slightly larger than the predicted rates. The fact that the GPS uplift rates are much larger closer to Jakobshavn than further away, and are consistent with rates inferred using the ATM-based glacier thinning model, shows that GPS measurements of crustal motion are a potentially useful method for assessing ice-mass change models.

Geometric parameters of shoreface profiles of Laptev Sea erosion coast (Table 1)

Field investigations of the Laptev Sea shoreface morphology were carried out (1) off erosional shores composed of unconsolidated sediments, (2) off the modern delta shores of the Lena River, and (3) off rocky shores. It was found that profiles off erosional shores had a concave shape. This shape is not well described by commonly applied power functions, a feature, which is in disagreement with the generally accepted concept of the equilibrium shape of shoreface profiles. The position of the lower shoreface boundary is determined by the elevation of the coastal lowland inundated during the last transgression (at -5 to -10 m) and may easily be recognized by a sharp, an order of magnitude decrease in the mean inclination of the sea floor. The mean shoreface inclination depends on sediment grain-size and ranges from 0.0022 to 0.033. The concave shape of the shoreface did not change substantially during the last 20-30 years, which indicates that shoreline retreat did not slow down and hence suggests continued intensive coastal erosion in the 21st century. The underwater part of the Lena River delta extends up to 35 km offshore. Its upper part is formed by a shallow and up to 18-km wide bench, which reaches depths of 2-3 m along the outer edge. The evolution of the delta was irregular. Whereas some parts of the delta are advancing rapidly (58 m/year), other parts are eroding. Comparison of measured profiles with older bathymetric data gave an opportunity to evaluate the changes of the underwater delta over past decades. Bathymetric surveys of the seabed around the delta can thus contribute towards a quantification of the sediment budget of the river-sea system. In addition, some sections of the Laptev Sea coast are composed of bedrock that has a comparatively low resistance to wave erosion. These sections may supply a considerable amount of sediment, especially if the cliffs are high. This source must therefore also be taken into account when assessing the contribution of shore erosion to the Laptev Sea sediment budget.

X-ray fluorescence (XRF) Ca intensity of sediment cores MD04-2760 and MD04-2788

Accelerator mass spectrometry (AMS) radiocarbon dating of ostracod and gastropod shells from the southwestern Black Sea cores combined with tephrochronology provides the basis for studying reservoir age changes in the lateglacial Black Sea. The comparison of our data with records from the northwestern Black Sea shows that an apparent reservoir age of ~1450 14C yr found in the glacial is characteristic of a homogenized water column. This apparent reservoir age is most likely due to the hardwater effect. Though data indicate that a reservoir age of ~1450 14C yr may have persisted until the Bølling-Allerød warm period, a comparison with the GISP2 ice-core record suggests a gradual reduction of the reservoir age to ~1000 14C yr, which might have been caused by dilution effects of inflowing meltwater. During the Bølling-Allerød warm period, soil development and increased vegetation cover in the catchment area of the Black Sea could have hampered erosion of carbonate bedrock, and hence diminished contamination by "old" carbon brought to the Black Sea basin by rivers. A further reduction of the reservoir age most probably occurred contemporary to the precipitation of inorganic carbonates triggered by increased phytoplankton activity, and was confined to the upper water column. Intensified deep water formation subsequently enhanced the mixing/convection and renewal of intermediate water. During the Younger Dryas, the age of the upper water column was close to 0 yr, while the intermediate water was ~900 14C yr older. The first inflow of saline Mediterranean water, at ~8300 14C yr BP, shifted the surface water age towards the recent value of ~400 14C yr.

Rocks and minerals from the Strakhov Fracture Zone, 4°N Mid-Atlantic Ridge

Geological-geophysical data obtained during Cruises 7, 11, and 12 of R/V Akademic Nikolay Strakhov (1989-1991) within the international project EQUARIDGE in the Strakhov Fracture Zone region (4°N) are presented. The trough of the fracture is interpreted as an open extension joint, a graben produced by stretching along the axis of the Mid-Atlantic Ridge. Bedrock studies showed that typical mid-ocean tholeiitic basalts occur within the narrow (60 nm wide) axial rift zone, whereas igneous rocks not typical for the ocean were found on the eastern and western flank plateaus. This allows to suppose that a reworked relict continental-type basement of pre-Upper Jurassic age possibly exists beneath the flank plateaus, within the segment under discussion. The above data correspond to the hypothesis of E. Bonatti about nonspreading nature of the basement of Mid-Atlantic Ridge within the equatorial segment and the Strakhov Fracture Zone.

Palynology on a profile from the Krusné Hory Mountains

Site details: The raised bog Fláje-Kiefern (50°429N, 13°329 E; 760 m a.s.l.; size ca. 500x500 m) lies in the Krusné Hory Mountains (Erzgebirge), Czech Republic, about 10 km from Georgenfelder Moor in Germany. Hejny and Slavík (1988) described the phytogeographic region of the Krusne Hory Mountains as 'a region of mountain flora and vegetation, with thermophilous species largely missing. In the natural forests, conifers, especially spruce (Picea excelsa) prevail. The deforested areas have been converted into meadows and pastures'. The climate is cool with annual average temperatures of about 5°C and annual precipitation of about 900 mm. The bedrock is Precambrian crystallinicum.