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.

Dinoflagellate cyst counts and age model from sediment samples of ODP Hole 175-1085A

Sediment samples from ODP Site 1085 were investigated in order to obtain more information on the initiation and development of the Benguela upwelling system during the middle and upper Miocene. In particular, our intent was to establish the causes of the upwelling as well as the response of the upwelling regime to the development of the Antarctic Circumpolar Current. Based on changes in the calcareous dinoflagellate cyst association, we found an initial increase of the dinoflagellate cyst productivity, probably related to the initiation of upwelling about 11.8 Ma ago. Two distinct increases in cyst productivity in conjunction with temperature decreases of the upper water masses reflect upwelling pulses off Namibia and occur at the end of the Miocene cooling events Mi5 (about 11.5 Ma) and Mi6 (about 10.5 Ma). Both cooling events are associated with an ice volume increase in Antarctica and are thought to have led to an increase in southeasterly winds, possibly causing these two upwelling pulses. We demonstrate a decrease in dinoflagellate cyst productivity and enhanced terrigenous input via the Orange River after the Mi5 event. At about 11.1 Ma, the dinoflagellate cyst productivity increases again. The polar cyst species Caracomia arctica occurs here for the first time. This implies an influence of subantarctic mode water and therefore a change in the quality of the upwelling water which allowed the Benguela upwelling to develop into modern conditions. From about 10.4 Ma, C. arctica forms a permanent part of the association, pointing to an establishment of the upwelling regime.

Planktonic foraminifera of the Mediterranean Sea

A high-resolution biochronology is presented for the Late Quaternary of the central Mediterranean. In the Late Pleistocene-Holocene successions three assemblage zones are distinguished on the basis of frequency patterns of planktic foraminifera. The age of these zones is determined by Accelerator Mass Spectrometry (AMS)14C dating. The zonal boundaries are dated at 12,700 yr B.P. (the end of Termination Ia) and 9600 yr B.P. (the start of Termination Ib), respectively. The AMS dates show that major changes in the planktic and benthic realms occurred synchronously over wide areas, although records of individual species may show important regional differences. In the studied areas, resedimentation processes revealed by anomalous successions of14C dates, play a far more important role than indicated by the sedimentological and micropaleontological data. Possibly these processes contribute to the very high accumulation rates in the glacial Zone III. Although the AMS technique has increased the accuracy of14C-measurements, admixture of older carbonate may still lead to substantial age differences between areas with different sedimentary regimes.

Carbon and oxygen isotopic compositions of Cibicidoides spp. and fine fraction separates (<10 µm) from ODP Holes 113-689B, 119-738B, and 120-748B

A prominent middle Eocene warming event is identified in Southern Ocean deep-sea cores, indicating that long-term cooling through the middle and late Eocene was not monotonic. At sites on Maud Rise and the Kerguelen Plateau, a distinct negative shift in d18O values (~1.0 per mil) is observed ca. 41.5 Ma. This excursion is interpreted as primarily a temperature signal, with a transient warming of 4°C over 600 k.y. affecting both surface and middle-bathyal deep waters in the Indian-Atlantic region of the Southern Ocean. This isotopic event is designated as the middle Eocene climatic optimum, and is interpreted to represent a significant climatic reversal in the midst of middle to late Eocene deep-sea cooling. The lack of a significant negative carbon isotope excursion, as observed during the Paleocene-Eocene thermal maximum, and the gradual rate of high-latitude warming suggest that this event was not triggered by methane hydrate dissociation. Rather, a transient rise in pCO2 levels is suspected, possibly as a result of metamorphic decarbonation in the Himalayan orogen or increased ridge/arc volcanism during the late middle Eocene.

Oxygen and carbon isotopic composition of shells from Radix sp. in lakes across the Tibetan Plateau

The Tibetan Plateau (TP), including its surrounding mountain ranges, represents the largest store of ice outside the polar regions. It hosts numerous lakes as well as the head waters of major Asian rivers, on which billions of people depend, and it is particularly sensitive to climate change. The moisture transport to the TP is controlled by the Indian and Pacific monsoon and the Westerlies. Understanding the evolution of the interaction of these circulation systems requires studies on climate archives in different spatial and temporal contexts. The objective of this study is to learn more about the interannual variability of precipitation patterns across the TP and how different hydrologic systems react to different climatic factors. Aragonite shells of the aquatic gastropod Radix, which is widely distributed in the region, may represent suitable archives for inferring hydrologic and climatic signals in particularly high resolution. Therefore, sclerochronological studies of d18O and d13C ratios in Radix shells from seven lakes were conducted, each representing a different hydrologic and climatic setting, on a transect from the Pamirs across the TP. The shell patterns exhibit an increasing influence of precipitation and a decreasing influence of evaporation on the isotope compositions from west to east. d18O values of shells from lakes on the eastern and central TP (Donggi Cona, Yamdrok Yumco, Tarab Co) mirror monsoon signals, indicated by more negative values and higher variabilities compared to the more western lakes (Karakul, Bangong/Nyak, Manasarovar). In Yadang Co, located on the central southern TP, the monsoon rains did not reach the lake in the sampling year, although it is located in a region which is usually affected by monsoon circulation. The d18O values are used to differentiate the annual hydrological cycle into ice cover period, melt water period, precipitation period and evaporation period. d13C compositions in the shells particularly depend on specific habitats, which vary in biological productivity and in carbon sources. d18O and d13C patterns show a positive covariance in shells originating from large closed basins. The results show that Radix shells mirror general climatic differences between the seven lake regions. These differences reflect both regional and local climate signals in sub-seasonal resolution, without noticeable dependence on the particular lake system.

Geochemistry and oxygen and iron isotope ratios of mafic rocks

Recycling of oceanic crust into the deep mantle via subduction is a widely accepted mechanism for creating compositional heterogeneity in the upper mantle and for explaining the distinct geochemistry of mantle plumes. The oxygen isotope ratios (d18O) of some ocean island basalts (OIB) span values both above and below that of unmetasomatised upper mantle (5.5 ± 0.4 per mil) and provide support for this hypothesis, as it is widely assumed that most variations in d18O are produced by near-surface low-temperature processes. Here we show a significant linear relationship between d18O and stable iron isotope ratios (d57Fe) in a suite of pristine eclogite xenoliths. The d18O values of both bulk samples and garnets range from values within error of normal mantle to significantly lighter values. The observed range and correlation between d18O and d57Fe is unlikely to be inherited from oceanic crust, as d57Fe values determined for samples of hydrothermally altered oceanic crust do not differ significantly from the mantle value and show no correlation with d18O. It is proposed that the correlated d57Fe and d18O variations in this particular eclogite suite are predominantly related to isotopic fractionation by disequilibrium partial melting although modification by melt percolation processes cannot be ruled out. Fractionation of Fe and O isotopes by removal of partial melt enriched in isotopically heavy Fe and O is supported by negative correlations between bulk sample d57Fe and Cr content and bulk sample and garnet d18O and Sc contents, as Cr and Sc are elements that become enriched in garnet- and pyroxene-bearing melt residues. Melt extraction could take place either during subduction, where the eclogites represent the residues of melted oceanic lithosphere, or could take place during long-term residence within the lithospheric mantle, in which case the protoliths of the eclogites could be of either crustal or mantle origin. This modification of both d57Fe and d18O by melting processes and specifically the production of low-d18O signatures in mafic rocks implies that some of the isotopically light d18O values observed in OIB and eclogite xenoliths may not necessarily reflect near-surface processes or components.

Chronostratigraphic framework at sites GeoB14403 and GeoB14414 from the Gela Basin, Sicily Channel

A multi-proxy chronological framework along with sequence-stratigraphic interpretations unveils composite Milankovitch cyclicity in the sedimentary records of the Last GlacialeInterglacial cycle at NE Gela Basin on the Sicilian continental margin. Chronostratigraphic data (including foraminifera-based eco-biostratigraphy and d18O records, tephrochronological markers and 14C AMS radiometric datings) was derived from the shallow-shelf drill sites GeoB14403 (54.6 m recovery) and GeoB14414 (27.5 m), collected with both gravity and drilled MeBo cores in 193 m and 146 m water depth, respectively. The recovered intervals record Marine Isotope Stages and Substages (MIS) from MIS 5 to MIS 1, thus comprising major stratigraphic parts of the progradational deposits that form the last 100-ka depositional sequence. Calibration of shelf sedimentary units with borehole stratigraphies indicates the impact of higher-frequency (20-ka) sea level cycles punctuating this 100-ka cycle. This becomes most evident in the alternation of thick interstadial highstand (HST) wedges and thinner glacial forced-regression (FSST) units mirroring seaward shifts in coastal progradation. Albeit their relatively short-lived depositional phase, these subordinate HST units form the bulk of the 100-ka depositional sequence. Two mechanisms are proposed that likely account for enhanced sediment accumulation ratios (SAR) of up to 200 cm/ka during these intervals: (1) intensified activity of deep and intermediate Levantine Intermediate Water (LIW) associated to the drowning of Mediterranean shelves, and (2) amplified sediment flux along the flooded shelf in response to hyperpycnal plumes that generate through extreme precipitation events during overall arid conditions. Equally, the latter mechanism is thought to be at the origin of undulated features resolved in the acoustic records of MIS 5 Interstadials, which bear a striking resemblance to modern equivalents forming on late-Holocene prodeltas of other Mediterranean shallow-shelf settings.

Chemobiostratigraphy of the Cretaceous/Paleocene boundary at ODP Hole 120-750A

An integrated biostratigraphic and stable isotope investigation was conducted on a high-latitude sequence across the Cretaceous/Paleogene (K/P) boundary recovered in Hole 750A in the southern Indian Ocean. The sequence consists of nannofossil chalk and is discontinuous across the boundary; missing is an estimated 0.3-m.y. late Maestrichtian and early Danian interval. Nonetheless, because calcareous nannofossil Zones NP1 and NP2 are well-developed, micropaleontological studies of the sequence have yielded a detailed record of Danian high-latitude microplankton evolution. In addition, stable carbon isotope analyses of planktonic and benthic foraminifer and bulk samples provide a record of late Maestrichtian and early Danian surface- and deep-water carbon isotope variations. Together, the carbon isotope and carbonate accumulation records serve as an index of regional marine net productivity across the boundary. Earliest Danian nannoplankton assemblages consisted mainly of persistent genera that were generally rare or absent in the Upper Cretaceous at Hole 750A. However, by 0.5-0.6 m.y. after the boundary, newly evolving Danian taxa became dominant. The turnover in nannofossil assemblages was accompanied by significant changes in rates of net productivity as gauged by carbon isotope distributions and carbonate accumulation rates. During the period dominated by persistent taxa, net productivity was extremely low, as reflected by the absence of vertical delta13C gradients and reduced carbonate accumulation rates. Later in the Danian, as new species evolved and flourished, vertical delta13C gradients reappeared and carbonate accumulation rates increased, signaling partial recovery of net productivity in this region. The absolute timing and magnitude of late Maestrichtian and early Danian biotic and geochemical changes in the southern Indian Ocean were similar to those recorded in other pelagic K/P boundary sequences from low- and mid-latitude Atlantic and Pacific sites, indicating that these events were ubiquitous.

Chronology for the EDML ice core from Dronning Maud Land, Antarctica, over the last 150 000 years

A chronology called EDML1 has been developed for the EPICA ice core from Dronning Maud Land (EDML). EDML1 is closely interlinked with EDC3, the new chronology for the EPICA ice core from Dome-C (EDC) through a stratigraphic match between EDML and EDC that consists of 322 volcanic match points over the last 128 ka. The EDC3 chronology comprises a glaciological model at EDC, which is constrained and later selectively tuned using primary dating information from EDC as well as from EDML, the latter being transferred using the tight stratigraphic link between the two cores. Finally, EDML1 was built by exporting EDC3 to EDML. For ages younger than 41 ka BP the new synchronized time scale EDML1/EDC3 is based on dated volcanic events and on a match to the Greenlandic ice core chronology GICC05 via 10Be and methane. The internal consistency between EDML1 and EDC3 is estimated to be typically ~6 years and always less than 450 years over the last 128 ka (always less than 130 years over the last 60 ka), which reflects an unprecedented synchrony of time scales. EDML1 ends at 150 ka BP (2417 m depth) because the match between EDML and EDC becomes ambiguous further down. This hints at a complex ice flow history for the deepest 350 m of the EDML ice core.