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.

(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.

Table 2 & Appendix A. Properties of operations injecting CO2 into saline aquifers

The experience from CO2 injection at pilot projects (Frio, Ketzin, Nagaoka, US Regional Partnerships) and existing commercial operations (Sleipner, Snøhvit, In Salah, acid-gas injection) demonstrates that CO2 geological storage in saline aquifers is technologically feasible. Monitoring and verification technologies have been tested and demonstrated to detect and track the CO2 plume in different subsurface geological environments. By the end of 2008, approximately 20 Mt of CO2 had been successfully injected into saline aquifers by existing operations. Currently, the highest injection rate and total storage volume for a single storage operation are approximately 1 Mt CO2/year and 25 Mt, respectively. If carbon capture and storage (CCS) is to be an effective option for decreasing greenhouse gas emissions, commercial-scale storage operations will require orders of magnitude larger storage capacity than accessed by the existing sites. As a result, new demonstration projects will need to develop and test injection strategies that consider multiple injection wells and the optimisation of the usage of storage space. To accelerate large-scale CCS deployment, demonstration projects should be selected that can be readily employed for commercial use; i.e. projects that fully integrate the capture, transport and storage processes at an industrial emissions source.

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.

(Table 2) Radiolaria species in nodule samples from the Clarion-Clipperton province, Pacific Ocean

Radiolaria were studied in 19 manganese nodules raised from the bottom. The nodules occurred mainly on the surface of thin Quaternary sediments covering Tertiary deposits of various ages (Middle Eocene to Early Miocene). Radiolaria in nodule cores and in inner and surface layers were studied. We found 85 radiolaria species and groups of species. Usually 1-4 to 6-19 radiolaria species were detected in each of the samples. Species belonging to Middle Eocene, Late Miocene to Early Oligocene, and Oligocene to Early Miocene were found. Rare Neogene species were revealed only in fractured surface layers. Age of the nodules is mainly Oligocene. Seismic waves cause sediment vibration, loosening disintegration, and removal of suspension by bottom currents. The vibration effect causes ancient nodules to float up to the surface of Quaternary sediment. This hypothesis suggests the reason for characteristics of the Clarion-Clipperton zone: regional stratigraphic hiatus, accumulation of residual fields of nodules, and the ''floating up'' of nodules to the surface of the Quaternary sediments.

Tab. 1 Mole % CO2 and 13C values of CO2 at baseline and at the final monitoring event for the eight observation wells

During CO2 storage operations in mature oilfields or saline aquifers it is desirable to trace the movement of injected CO2 for verification and safety purposes. We demonstrate the successful use of carbon isotope abundance ratios for tracing the movement of CO2 injected at the Cardium CO2 Storage Monitoring project in Alberta between 2005 and 2007. Injected CO2 had a d13C value of -4.6±1.1 per mil that was more than 10 per mil higher than the carbon isotope ratios of casing gas CO2 prior to CO2 injection with average d13C values ranging from -15.9 to -23.5 per mil. After commencement of CO2 injection, d13C values of casing gas CO2 increased in all observation wells towards those of the injected CO2 consistent with a two-source end-member mixing model. At four wells located in a NE-SW trend with respect to the injection wells, breakthrough of injected CO2 was registered chemically (>50 mol % CO2) and isotopically 1-6 months after commencement of CO2 injection resulting in cumulative CO2 fluxes exceeding 100000 m**3 during the observation period. At four other wells, casing gas CO2 contents remained below 5 mol % resulting in low cumulative CO2 fluxes (<2000 m**3) throughout the entire observation period, but carbon isotope ratios indicated contributions between <30 and 80% of injected CO2. Therefore, we conclude that monitoring the movement of CO2 in the injection reservoir with geochemical and isotopic techniques is an effective approach to determine plume expansion and to identify potential preferential flow paths provided that the isotopic composition of injected CO2 is constant and distinct from that of baseline CO2.

Composition of soil samples from injection location

CO2 leakage from subsurface storage sites is one of the main concerns connected with the CCS technology. As CO2 leakages into near surface formations appear to be very unlikely within pilot CCS projects, the aim of this work is to emulate a leakage by injecting CO2 into a near surface aquifer. The two main questions pursued by the injection test are (1) to investigate the impact of CO2 on the hydrogeochemistry of the groundwater as a base for groundwater risk assessment and (2) to develop and apply monitoring methods and monitoring concepts for detecting CO2 leakages in shallow aquifers. The presented injection test is planned within the second half of 2010, as a joint project of the University of Kiel (Germany), the Helmholtz-Centre for Environmental Research (Leipzig, Germany) and the Engineering Company GICON (Dresden, Germany). The test site has been investigated in detail using geophysical methods as well as direct-push soundings, groundwater well installation and soil and groundwater analyses. The present paper presents briefly the geological and hydrogeological conditions at the test site as well as the planned injection test design and monitoring concept.

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.