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.

Table 1 Microorganisms, identified with DOGE analysis. 16S rDNA sequences, obtained from sequencing of DOGE bands, were assigned to the corresponding microorganisms

Microorganisms play an important role in the transformation of material within the earth's crust. The storage of CO2 could affect the composition of inorganic and organic components in the reservoir, consequently influencing microbial activities. To study the microbial induced processes together with geochemical, petrophysical and mineralogical changes, occurring during CO2 storage, long-term laboratory experiments under simulated reservoir P-T conditions were carried out. Clean inner core sections, obtained from the reservoir region at the CO2 storage site in Ketzin (Germany) from a depth of about 650 m, were incubated in high pressure vessels together with sterile synthetic formation brine under in situ P-T conditions of 5.5 MPa and 40°C. A 16S rDNA based fingerprinting method was used to identify the dominant species in DNA extracts of pristine sandstone samples. Members of the alpha- and beta-subdivisions of Proteobacteria and the Actinobacteria were identified. So far sequences belonging to facultative anaerobic, chemoheterotrophic bacteria (Burkholderia fungorum, Agrobacterium tumefaciens) gaining their energy from the oxidation of organic molecules and a genus also capable of chemolithoautotrophic growth (Hydrogenophaga) was identified. During CO2 incubation minor changes in the microbial community composition were observed. The majority of microbes were able to adapt to the changed conditions. During CO2 exposure increased concentrations of Ca**2+, K**+, Mg**2+ and SO4**2- were observed. Partially, concentration rises are (i) due to equilibration between rock pore water and synthetic brine, and (ii) between rock and brine, and are thus independent on CO2 exposure. However, observed concentrations of Ca**2+, K**+, Mg**2+ are even higher than in the original reservoir fluid and therefore indicate mineral dissolution due to CO2 exposure.

Icebergs in the Amundsen Sea (Antarctica) determined from ENVISAT ASAR WSM data using object-based image analysis approach, link to ESRI GDB files

ENVISAT ASAR WSM images with pixel size 150 × 150 m, acquired in different meteorological, oceanographic and sea ice conditions were used to determined icebergs in the Amundsen Sea (Antarctica). An object-based method for automatic iceberg detection from SAR data has been developed and applied. The object identification is based on spectral and spatial parameters on 5 scale levels, and was verified with manual classification in four polygon areas, chosen to represent varying environmental conditions. The algorithm works comparatively well in freezing temperatures and strong wind conditions, prevailing in the Amundsen Sea during the year. The detection rate was 96% which corresponds to 94% of the area (counting icebergs larger than 0.03 km**2), for all seasons. The presented algorithm tends to generate errors in the form of false alarms, mainly caused by the presence of ice floes, rather than misses. This affects the reliability since false alarms were manually corrected post analysis.

Ocean acidification not likely to affect the survival and fitness of two yemperate benthic foraminiferal species: results from culture experiments

Specimens of Bolivina argentea and Bulimina marginata, two widely distributed temperate benthic foraminiferal species, were cultured at constant temperature and controlled pCO2 (ambient, 1000 ppmv, and 2000 ppmv) for six weeks to assess the effect of elevated atmospheric CO2 concentrations on survival and fitness using Adenosine Triphosphate (ATP) analyses and on shell microfabric using high-resolution SEM and image analysis. To characterize the carbonate chemistry of the incubation seawater, total alkalinity and dissolved inorganic carbon were measured approximately every two weeks. Survival and fitness were not directly affected by elevated pCO2 and the concomitant decrease in seawater pH and calcite saturation states (Omega c), even when seawater was undersaturated with respect to calcite. These results differ from some previous observations that ocean acidification can cause a variety of effects on benthic foraminifera, including test dissolution, decreased growth, and mottling (loss of symbiont color in symbiont-bearing species), suggesting that the benthic foraminiferal response to ocean acidification may be species specific. If so, this implies that ocean acidification may lead to ecological winners and losers even within the same taxonomic group.

Combined metabolome and proteome analysis of the mantle tissue from Pacific oyster Crassostrea gigas exposed to elevated pCO2

Ocean acidification (OA) has been found to affect an array of normal physiological processes in mollusks, especially posing a significant threat to the fabrication process of mollusk shell. In the current study, the impact of exposure to elevated pCO2 condition was investigated in mantle tissue of Crassostrea gigas by an integrated metabolomic and proteomic approach. Analysis of metabolome and proteome revealed that elevated pCO2 could affect energy metabolism in oyster C. gigas, marked by differentially altered ATP, succinate, MDH, PEPCK and ALDH levels. Moreover, the up-regulated calponin-2, tropomyosins and myosin light chains indicated that elevated pCO2 probably caused disturbances in cytoskeleton structure in mantle tissue of oyster C. gigas. This work demonstrated that a combination of proteomics and metabolomics could provide important insights into the effects of OA at molecular levels.

Sea surface temperature and age analysis from sediment core KT94-15_PC-9, XP98-PC-4 and XP98-PC-2

The Japan Sea experienced bottom water anoxia at the last glacial maximum (LGM) since it is surrounded by four shallow straits, the sill depths of which are close to, or shallower than, the drop in sea level (~120 m) that occurred then. A distinctive negative d18O excursion of planktonic foraminifera also took place during the LGM. This excursion has been interpreted from foraminiferal data as recording a drop in the paleosalinity of surface waters on the assumption of a constant low sea surface temperatures between 34 and 11 ka. We present here a profile of alkenone-based sea surface temperatures (alkenone-SSTs) over the past 36 kyr. Our results suggest that SSTs during the LGM were much higher than those previously assumed. After considering the factors that might affect estimation of alkenone-SSTs and comparisons of core-top alkenone-SSTs values with values for modern seawater we conclude that the higher alkenone-SSTs during the LGM are reliable and reasonable. These warm SSTs were probably caused by radiative equilibrium associated with the development of stable water stratification in the Japan Sea during the LGM.

Analysis of minerals from Arctic Ocean sediments

The Arctic Ocean System is a key player regarding the climatic changes of Earth. Its highly sensitive ice Cover, the exchange of surface and deep water masses with the global ocean and the coupling with the atmosphere interact directly with global climatic changes. The output of cold, polar water and sea ice influences the production of deep water in the North Atlantic and controls the global ocean circulation ("the conveyor belt"). The Arctic Ocean is surrounded by the large Northern Hemisphere ice sheets which not only affect the sedimentation in the Arctic Ocean but also are supposed to induce the Course of glacials and interglacials. Terrigenous sediment delivered from the ice sheets by icebergs and meltwater as well as through sea ice are major components of Arctic Ocean sediments. Hence, the terrigenous content of Arctic Ocean sediments is an outstanding archive to investigate changes in the paleoenvironment. Glazigenic sediments of the Canadian Arctic Archipelago and surface samples of the Arctic Ocean and the Siberian shelf regions were investigated by means of x-ray diffraction of the bulk fraction. The source regions of distinct mineral compositions were to be deciphered. Regarding the complex circumpolar geology stable christalline shield rocks, active and ancient fold belts including magmatic and metamorphic rocks, sedimentary rocks and wide periglacial lowlands with permafrost provide a complete range of possible mineral combinations. Non- glaciated shelf regions mix the local input from a possible point source of a particular mineral combination with the whole shelf material and function as a sampler of the entire region draining to the shelf. To take this into account, a literature research was performed. Descriptions of outcropping lithologies and Arctic Ocean sediments were scanned for their mineral association. The analyses of glazigenic and shelf sediments yielded a close relationship between their mineral composition and the adjacent source region. The most striking difference between the circumpolar source regions is the extensive outcrop of carbonate rocks in the vicinity of the Canadian Arctic Archipelago and in N Greenland while siliciclastic sediments dominate the Siberian shelves. In the Siberian shelf region the eastern Kara Sea and the western Laptev Sea form a destinct region defined by high smectite, (clino-) pyroxene and plagioclase input. The source of this signal are the extensive outcrops of the Siberian trap basalt in the Putorana Plateau which is drained by the tributaries of the Yenissei and Khatanga. The eastern Laptev Sea and the East Siberian Sea can also be treated as one source region containing a feldspar, quartz, illite, mica, and chlorite asscciation combined with the trace minerals hornblende and epidote. Franz Josef Land provides a mineral composition rich in quartz and kaolinite. The diverse rock suite of the Svalbard archipelago distributes specific mineral compositions of highly metamorphic christalline rocks, dolomite-rich carbonate rocks and sedimentary rocks with a higher diagenetic potential manifested in stable newly built diagenetic minerals and high organic maturity. To reconstruct the last 30,000 years as an example of the transition between glacial and interglacial conditions a profile of sediment cores, recovered during the RV Polarstern" expedition ARK-VIIIl3 (ARCTIC '91), and additional sediment cores around Svalbard were investigated. Besides the mineralogy of different grain size fractions several additional sedimentological and organo-geochemical Parameterswere used. A detailed stratigraphic framework was achieved. By exploiting this data set changes in the mineral composition of the Eurasian Basin sediments can be related to climatic changes. Certain mineral compositions can even be associated with particular transport processes, e.g. the smectitel pyroxene association with sea ice transport from the eastern Kara Sea and the western Laptev Sea. Hence, it is possible to decipher the complex interplay between the influx of warm Atlantic waters into the Southwest of the Eurasian Basin, the waxing and waning of the Svalbard1Barents- Sea- and Kara-Sea-Ice-Sheets, the flooding of the Siberian shelf regions and the surface and deep water circulation. Until now the Arctic Ocean was assumed to be a rather stable System during the last 30,000 years which only switched from a completely ice covered situation during the glacial to seasonally Open waters during the interglacial. But this work using mineral assemblages of sediment cores in the vicinity of Svalbard revealed fast changes in the inflow of warm Atlantic water with the Westspitsbergen Current (< 1000 years), short periods of advances and retreats of the marine based Eurasian ice sheets (1000-3000 years), and short melting phases (400 years?). Deglaciation of the marine-based Eurasian and the land-based north American and Greenland ice sheets are not simultaneous. This thesis postulates that the Kara Sea Ice Sheet released an early meltwater signal prior to 15,000 14C years leading the Barents Sea Ice Sheet while the western land-based ice sheets are following later than 13,500 14C years. The northern Eurasian Basin records the shift between iceberg and sea-ice material derived from the Canadian Arctic Archipelago and N-Greenland and material transported by sea-ice and surface currents from the Siberian shelf region. The phasing of the deglaciation becomes very obvious using the dolomite and quartd phyllosilicate record. It is also supposed that the flooding of the Laptev Sea during the Holocene is manifested in a stepwise increase of sediment input at the Lomonosov Ridge between the Eurasian and Amerasian Basin. Depending on the strength of meltwater pulses from the adjacent ice sheets the Transpolar Drift can probably be relocated. These movements are traceable by the distribution of indicator minerals. Based on the outcome of this work the feasibility of bulk mineral determination can be qualified as excellent tool for paleoenvironmental reconstructions in the Arctic Ocean. The easy preparation and objective determination of bulk mineralogy provided by the QUAX software bears the potential to use this analyses as basic measuring method preceding more time consuming and highly specialised mineralogical investigations (e.g. clay mineralogy, heavy mineral determination).

Results of analysis and model simulation of the bed profiles

Bedforms such as dunes and ripples are ubiquitous in rivers and coastal seas, and commonly described as triangular shapes from which height and length are calculated to estimate hydrodynamic and sediment dynamic parameters. Natural bedforms, however, present a far more complicated morphology; the difference between natural bedform shape and the often assumed triangular shape is usually neglected, and how this may affect the flow is unknown. This study investigates the shapes of natural bedforms and how they influence flow and shear stress, based on four datasets extracted from earlier studies on two rivers (the Rio Paraná in Argentina, and the Lower Rhine in The Netherlands). The most commonly occurring morphological elements are a sinusoidal stoss side made of one segment and a lee side made of two segments, a gently sloping upper lee side and a relatively steep (6 to 21°) slip face. A non-hydrostatic numerical model, set up using Delft3D, served to simulate the flow over fixed bedforms with various morphologies derived from the identified morphological elements. Both shear stress and turbulence increase with increasing slip face angle and are only marginally affected by the dimensions and positions of the upper and lower lee side. The average slip face angle determined from the bed profiles is 14°, over which there is no permanent flow separation. Shear stress and turbulence above natural bedforms are higher than above a flat bed but much lower than over the often assumed 30° lee side angle.

Core scans and lithologic analysis of Subglacial Lake Whillans sediment cores

The hydrologic system beneath the Antarctic Ice Sheet is thought to influence both the dynamics and distribution of fast flowing ice streams, which discharge most of the ice lost by the ice sheet. Despite considerable interest in understanding this subglacial network and its affect on ice flow, in situ observations from the ice sheet bed are exceedingly rare. Here we describe the first sediment cores recovered from an active subglacial lake. The lake, known as Subglacial Lake Whillans, is part of a broader, dynamic hydrologic network beneath the Whillans Ice Stream in West Antarctica. Even though "floods" pass through the lake, the lake floor shows no evidence of erosion or deposition by flowing water. By inference, these floods must have insufficient energy to erode or transport significant volumes of sediment coarser than silt. Consequently, water flow beneath the region is probably incapable of incising continuous channels into the bed and instead follows preexisting subglacial topography and surface slope. Sediment on the lake floor consists of till deposited during intermittent grounding of the ice stream following flood events. The fabrics within the till are weaker than those thought to develop in thick deforming beds suggesting subglacial sediment fluxes across the ice plain are currently low and unlikely to have a large stabilizing effect on the ice stream's grounding zone.