The role of climate-driven chemical weathering on soil production

Climate plays an important role in controlling rates of weathering and weathered regolith production. Regolith production functions, however, seldom take climate parameters into account. Based on a climate-dependent weathered regolith production model, at low denudation rates, relative regolith thicknesses are less sensitive to changes in precipitation rates, while at high denudation rates, small changes in climatic parameters can result in complete stripping of hillslopes. This pattern is compounded by the long residence times and system response times associated with low denudation rates, and vice versa. As others have shown, the transition between regolith-mantled and bedrock slopes is dependent on the ratio of denudation to production. Here, we further suggest that this is itself a function of precipitation rate and temperature. We suggest that climatic parameters can be easily incorporated into existing soil production models and that such additions improve the predictive power of soil production models. (C) 2013 Elsevier B.V. All rights reserved.

Drilling into debris-rich basal ice at the bottom of the NEEM (Greenland) borehole

After the NEEM (Greenland) deep ice-core drilling was declared terminated with respect to developing stratigraphic climate reconstructions, efforts were turned toward collecting basal ice-sheet debris and, if possible, drilling into the bedrock itself. In 2010, several meters of banded debris-rich ice were obtained under normal ice-drilling operations with the NEEM version of the Hans Tausen (HT) drill, but further penetration was obstructed by a rock in the path of the drill head at 2537.36 m. During short campaigns in 2011 and 2012, attempts were made to penetrate further using various reinforced ice cutters mounted on the HT drill head, tailored to cut through rock. These had some success in penetrating coarse material, but produced severely damaged cutters. Additionally a 51 mm diameter diamond cutting tipped rock drill was adapted to fit the NEEM drill. With this device, several additional meters of core containing subglacial sediments, rocks and rock fragments were collected. With these tools 1.39 m of additional material were obtained during the 2011 field season, and 7.1 m during 2012. Subglacial water refreezing into the newly formed borehole hindered further penetration, and the bedrock interface was not reached before final closure of the NEEM Camp.

Ce(III) and Ce(IV) (re)distribution and fractionation in a laterite profile from Madagascar: Insights from in situ XANES spectroscopy at the Ce LIII-edge

The distribution of trivalent and tetravalent cerium, Ce(III) and Ce(IV) respectively, in a lateritic profile from Madagascar, has been characterized by X-ray-absorption near-edge structure (XANES) spectroscopy at the Ce LIII-edge on the LUCIA beamline (SOLEIL synchrotron, France). XANES spectra were acquired on bulk-rock samples as well as on specific lateritic minerals or polymineral zones (in-situ measurements) of the tonalite bedrock and the three overlying weathered horizons (C-, B- and A-horizons). Geochemically, the bedrock, and the A- and C-horizons show similar rare earth element content (REE = 363–405 mg/kg). They also display the same positive Ce-anomaly (CeCN/Ce∗ = 1.12–1.45), which is therefore likely to be inherited from the bedrock. In the B-horizon, the higher REE content (REE = 2194 mg/kg) and the larger Ce-anomaly (CeCN/Ce∗ = 4.26) are consistent with an accumulation zone caused by the evaporation of groundwater during the dry season. There is a good agreement between the Ce(III)/Cetotal ratio (XCe(III)) deduced from the positive Ce-anomaly (bulk-rock geochemical data) and that derived from XANES spectroscopy on the same bulk-rock samples (BR-XCe(III)-XANES) in the bedrock, and the C- and B-horizons. In the A-horizon, XANES measurements on bulk rock and minerals revealed a higher BR-XCe(III)-XANES (up to 100%) compared to the XCe(III) deduced from geochemical data (XCe(III) = 79%). The preservation of a positive Ce-anomaly in the A-horizon suggests that the Ce mobilization and redistribution during weathering occurred with no significant Ce fractionation from other trivalent REE. Remarkably, the only investigated sample where cerianite is observed belongs to the B-horizon. Within this horizon, Ce oxidation state varies depending on the microstructural position (porosity, cracks, clay-rich groundmass). The highest Ce(IV) concentrations are measured in cerianite (and aluminophosphates) localized in pores at the vicinity of Mn-rich domains (XCe(III)-XANES = 30–51%). Therefore, Ce fractionation from other REE is attributed to a Ce oxidation and precipitation potentially assisted by oxyhydroxide scavenging. In the C-horizon, Ce(III) and Ce(IV) are mainly distributed in REE-minerals of the rhabdophane group found in pores and cracks. The similarity between the Ce(III) proportion of rhabdophane grains (XCe(III)-XANES = 74–89%) with that of the bedrock (BR-XCe(III)-XANES = 79%) suggests no significant fractionation of Ce(III) and Ce(IV) between solution and mineral during the successive stages of primary REE-mineral alteration, transport in solution and secondary precipitation in the incipient stages of weathering. Overall, our novel spectroscopic approach shows that Ce is not necessarily oxidized nor fractionated from other REE during weathering in lateritic conditions. This implies that like Ce(III), Ce(IV) can be mobilized in aqueous fluids during weathering, possibly thanks to complexation with organic molecules, and can precipitate together with Ce(III) in secondary REE-bearing minerals. The corollary is that (paleo)redox reconstructions in soils and/or sediments based on Ce-anomaly in weathered rocks or minerals must be interpreted with caution.

Efficiency of frost-cracking processes through space and time: An example from the eastern Italian Alps

It is widely accepted that climate has a strong impact and exerts important feedbacks on erosional processes and sediment transport mechanisms. However, the extent at which climate influences erosion is still a matter of debate. In this paper we test whether frost-cracking processes and related temperature variations can influence the sediment production and surface erosion in a small catchment situated in the eastern Italian Alps. To this extent, we first present a geomorphic map of the region that we complement with published 10Be-based denudation rates. We then apply a preexisting heat-flow model in order to analyze the variations of the frost-cracking intensity (FCI) in the study area, which could have controlled the sediment production in the basin. Finally, we compare the model results with the pattern of denudation rates and Quaternary deposits in the geomorphic map. The model results, combined with field observations, mapping, and quantitative geomorphic analyses, reveal that frost-cracking processes have had a primary role in the production of sediment where the intensity of sediment supply has been dictated and limited by the combined effect of temperature variations and conditions of bedrock preservation. These results highlight the importance of a yet poorly understood process for the production of sediment in mountain areas.

The first major incision of the Swiss Deckenschotter landscape

The Swiss Deckenschotter (“cover gravels”) is the oldest Quaternary units in the northern Swiss Alpine Foreland. They are a succession of glaciofluvial gravel layers intercalated with glacial and/or overbank deposits. This lithostratigraphic sequence is called Deckenschotter because it “covers” Molasse or Mesozoic bedrock and forms mesa-type hill-tops. Deckenschotter occurs both within and beyond the extent of the Last Glacial Maximum glaciers. The Swiss Deckenschotter consist of two sub-units: Höhere (Higher) and Tiefere (Lower) Deckenschotter. Although the Höhere Deckenschotter sub-unit (HDS) is topographically higher than the Tiefere Deckenschotter, it is older. The only available age for the Swiss Deckenschotter is 2.5–1.8 Ma based on mammal remains found in HDS at the Irchel site. In this study, we present an exposure age for the topographically lowest HDS, calculated from a cosmogenic 10Be depth-profile. Our results show that the first phase of the Deckenschotter glaciations in the Swiss Alps terminated at least 1,020+80−120 ka ago, which is indicated by a significant fluvial incision. This line of evidence seems to be close to synchronous with the beginning of the Mid-Pleistocene Revolution, when the frequency of the glacial-interglacial cyclicity changed from 41 to 100 ka and the amplitude from low to high, between marine isotope stages 23 and 22.

Chronology of Lateglacial ice flow reorganization and deglaciation in the Gotthard Pass area, Central Swiss Alps, based on cosmogenic 10Be and in situ 14C

We reconstruct the timing of ice flow reconfiguration and deglaciation of the Central Alpine Gotthard Pass, Switzerland, using cosmogenic 10Be and in situ14C surface exposure dating. Combined with mapping of glacial erosional markers, exposure ages of bedrock surfaces reveal progressive glacier downwasting from the maximum LGM ice volume and a gradual reorganization of the paleoflow pattern with a southward migration of the ice divide. Exposure ages of ∼16–14 ka (snow corrected) give evidence for continuous early Lateglacial ice cover and indicate that the first deglaciation was contemporaneous with the decay of the large Gschnitz glacier system. In agreement with published ages from other Alpine passes, these data support the concept of large transection glaciers that persisted in the high Alps after the breakdown of the LGM ice masses in the foreland and possibly decayed as late as the onset of the Bølling warming. A younger group of ages around ∼12–13 ka records the timing of deglaciation following local glacier readvance during the Egesen stadial. Glacial erosional features and the distribution of exposure ages consistently imply that Egesen glaciers were of comparatively small volume and were following a topographically controlled paleoflow pattern. Dating of a boulder close to the pass elevation gives a minimum age of 11.1 ± 0.4 ka for final deglaciation by the end of the Younger Dryas. In situ14C data are overall in good agreement with the 10Be ages and confirm continuous exposure throughout the Holocene. However, in situ14C demonstrates that partial surface shielding, e.g. by snow, has to be incorporated in the exposure age calculations and the model of deglaciation.

Retrieving the paleoclimatic signal from the deeper part of the EPICA Dome C ice core

An important share of paleoclimatic information is buried within the lowermost layers of deep ice cores. Because improving our records further back in time is one of the main challenges in the near future, it is essential to judge how deep these records remain unaltered, since the proximity of the bedrock is likely to interfere both with the recorded temporal sequence and the ice properties. In this paper, we present a multiparametric study (δD-δ18Oice, δ18Oatm, total air content, CO2, CH4, N2O, dust, high-resolution chemistry, ice texture) of the bottom 60 m of the EPICA (European Project for Ice Coring in Antarctica) Dome C ice core from central Antarctica. These bottom layers were subdivided into two distinct facies: the lower 12 m showing visible solid inclusions (basal dispersed ice facies) and the upper 48 m, which we will refer to as the "basal clean ice facies". Some of the data are consistent with a pristine paleoclimatic signal, others show clear anomalies. It is demonstrated that neither large-scale bottom refreezing of subglacial water, nor mixing (be it internal or with a local basal end term from a previous/initial ice sheet configuration) can explain the observed bottom-ice properties. We focus on the high-resolution chemical profiles and on the available remote sensing data on the subglacial topography of the site to propose a mechanism by which relative stretching of the bottom-ice sheet layers is made possible, due to the progressively confining effect of subglacial valley sides. This stress field change, combined with bottom-ice temperature close to the pressure melting point, induces accelerated migration recrystallization, which results in spatial chemical sorting of the impurities, depending on their state (dissolved vs. solid) and if they are involved or not in salt formation. This chemical sorting effect is responsible for the progressive build-up of the visible solid aggregates that therefore mainly originate "from within", and not from incorporation processes of debris from the ice sheet's substrate. We further discuss how the proposed mechanism is compatible with the other ice properties described. We conclude that the paleoclimatic signal is only marginally affected in terms of global ice properties at the bottom of EPICA Dome C, but that the timescale was considerably distorted by mechanical stretching of MIS20 due to the increasing influence of the subglacial topography, a process that might have started well above the bottom ice. A clear paleoclimatic signal can therefore not be inferred from the deeper part of the EPICA Dome C ice core. Our work suggests that the existence of a flat monotonic ice–bedrock interface, extending for several times the ice thickness, would be a crucial factor in choosing a future "oldest ice" drilling location in Antarctica.

The 1987 debris flows in Switzerland: documentation and analysis

A great number of debris flows occurred during the flood catastrophes of the summer of 1987 in the Swiss Alps. Aerial photography, field investigations and eyewitness accounts documented and analysed the events. As an example of the reconstructed major events, the large debris flow in the Varuna valley involved an estimated peak discharge between 400 and 800 m3/s and an event magnitude of 200,000 m3. Several single pulses were observed; the duration of each of them appeared to be not more than a few minutes. Apart from incision into weak bedrock, the maximum erosion depth seemed to depend on the channel gradient. Based on approximately 600 events, typical starting zones and rainfall conditions are discussed with regard to the triggering conditions. Existing and new empirical formulae are proposed to estimate the most important flow parameters. These values are compared to debris flow data from Canada and Japan.

Reconstruction of in-situ porosity and porewater compositions of low-permeability crystalline rocks: Magnitude of artefacts induced by drilling and sample recovery

Geological site characterisation programmes typically rely on drill cores for direct information on subsurface rocks. However, porosity, transport properties and porewater composition measured on drill cores can deviate from in-situ values due to two main artefacts caused by drilling and sample recovery: (1) mechanical disruption that increases porosity and (2) contamination of the porewater by drilling fluid. We investigated the effect and magnitude of these perturbations on large drill core samples (12–20 cm long, 5 cmdiameter) of high-grade, granitic gneisses obtained from 350 to 600 m depth in a borehole on Olkiluoto Island (SW Finland). The drilling fluid was traced with sodium–iodide. By combining out-diffusion experiments, gravimetry, UV-microscopy and iodide mass balance calculations, we successfully quantified the magnitudes of the artefacts: 2–6% increase in porosity relative to the bulk connected porosity and 0.9 to 8.9 vol.% contamination by drilling fluid. The spatial distribution of the drilling-induced perturbations was revealed by numerical simulations of 2D diffusion matched to the experimental data. This showed that the rims of the samples have a mechanically disrupted zone 0.04 to 0.22 cm wide, characterised by faster transport properties compared to the undisturbed centre (1.8 to 7.7 times higher pore diffusion coefficient). Chemical contamination was shown to affect an even wider zone in all samples, ranging from 0.15 to 0.60 cm, inwhich iodide enrichmentwas up to 180 mg/kgwater, compared to 0.5 mg/kgwater in the uncontaminated centre. For all samples in the present case study, it turned out that the magnitude of the artefacts caused by drilling and sample recovery is so small that no correction is required for their effects. Therefore, the standard laboratory measurements of porosity, transport properties and porewater composition can be taken as valid in-situ estimates. However, it is clear that the magnitudes strongly depend on site- and drilling-specific factors and therefore our results cannot be transferred simply to other locations. We recommend the approach presented in this study as a route to obtain reliable values in future drilling campaigns aimed at characterising in-situ bedrock properties.

Provenance analysis using Raman spectroscopy of carbonaceous material: A case study in the Southern Alps of New Zealand

Detrital provenance analyses in orogenic settings, in which sediments are collected at the outlet of a catchment, have become an important tool to estimate how erosion varies in space and time. Here we present how Raman Spectroscopy on Carbonaceous Material (RSCM) can be used for provenance analysis. RSCM provides an estimate of the peak temperature (RSCM-T) experienced during metamorphism. We show that we can infer modern erosion patterns in a catchment by combining new measurements on detrital sands with previously acquired bedrock data. We focus on the Whataroa catchment in the Southern Alps of New Zealand and exploit the metamorphic gradient that runs parallel to the main drainage direction. To account for potential sampling biases, we also quantify abrasion properties using flume experiments and measure the total organic carbon content in the bedrock that produced the collected sands. Finally, we integrate these parameters into a mass-conservative model. Our results first demonstrate that RSCM-T can be used for detrital studies. Second, we find that spatial variations in tracer concentration and erosion have a first-order control on the RSCM-T distributions, even though our flume experiments reveal that weak lithologies produce substantially more fine particles than do more durable lithologies. This result implies that sand specimens are good proxies for mapping spatial variations in erosion when the bedrock concentration of the target mineral is quantified. The modeling suggests that highest present-day erosion rates (in Whataroa catchment) are not situated at the range front but around 10 km into the mountain belt.

Bedload composition, transport and modification in rivers of Westland, New Zealand, with implications for the distribution of alluvial pounamu (jade)

River bedload surveyed at 50 sites in Westland is dominated by Alpine Schist or Torlesse Greywacke from the Alpine Fault hanging wall, with subordinate Pounamu Ultramafics or footwall-derived Western Province rocks. Tumbling experiments found ultramafics to have the lowest attrition rates, compared with greywacke sandstone and granite (which abrade to produce silt to medium-sand), or incompetent schist (which fragments). Arahura has greater total concentrations (103–105 t/km2) and proportions (5–40%) of ultramafic bedload compared with Hokitika and Taramakau catchments (101–104 t/km2, mostly <10%), matching relative areas of mapped Pounamu Ultramafic bedrock, but enriched relative to absolute areal proportions. Western Province rocks downthrown by the Alpine Fault are under-represented in the bedload. Enriched concentrations of ultramafic bedload decrease rapidly with distance downstream from source rock outcrops, changing near prominent ice-limit moraines. Bedload evolution with transport involves both downstream fining and dilution from tributaries, in a sediment supply regime more strongly influenced by tectonics and the imprint of past glaciation. Treasured New Zealand pounamu (jade) is associated with ultramafic rocks. Chances of discovery vary between catchments, are increased near glacial moraines, and are highest near source-rock outcrops in remote mountain headwaters.

Constraining age distributions of groundwater from public supply wells in diverse hydrogeological settings in Scania, Sweden

Twenty-five public supply wells throughout the hydrogeologically diverse region of Scania, southern Sweden are subjected to environmental tracer analysis (³H–³He,⁴He, CFCs, SF₆ and for one well only also ⁸⁵Kr and ³⁹Ar) to study well and aquifer vulnerability and evaluate possibilities of groundwater age distribution assessment. We find CFC and SF₆ concentrations well above solubility equilibrium with modern atmosphere, indicating local contamination, as well as indications of CFC degradation. The tracer-specific complications considerably constrain possibilities for sound quantitative regional ground- water age distribution assessment and demonstrate the importance of initial qualitative assessment of tracer-specific reliability, as well a need for additional, complementary tracers (e.g. ⁸⁵Kr,³⁹Ar and potentially also ¹⁴C). Lumped parameter modelling yields credible age distribution assessments for representative wells in four type aquifers. Pollution vulnerability of the aquifer types was based on the selected LPM models and qualitative age characterisation. Most vulnerable are unconfined dual porosity and fractured bedrock aquifers, due to a large component of very young groundwater. Unconfined sedimentary aquifers are vulnerable due to young groundwater and a small pre-modern component. Less vulnerable are semi-confined sedimentary or dual-porosity aquifers, due to older age of the modern component and a larger pre-modern component. Confined aquifers appear least vulnerable, due an entirely pre-modern groundwater age distribution (recharged before 1963). Tracer complications aside, environmental tracer analyses and lumped parameter modelling aid in vulnerability assessment and protection of regional groundwater resources.

Patterns of landscape form in the upper Rhône basin, Central Swiss Alps, predominantly show lithologic controls despite multiple glaciations and variations in rock uplift rates

The development of topography depends mainly on the interplay between uplift and erosion. These processes are controlled by various factors including climate, glaciers, lithology, seismic activity and short-term variables, such as anthropogenic impact. Many studies in orogens all over the world have shown how these controlling variables may affect the landscape's topography. In particular, it has been hypothesized that lithology exerts a dominant control on erosion rates and landscape morphology. However, clear demonstrations of this influence are rare and difficult to disentangle from the overprint of other signals such as climate or tectonics. In this study we focus on the upper Rhône Basin situated in the Central Swiss Alps in order to explore the relation between topography, possible controlling variables and lithology in particular. The Rhône Basin has been affected by spatially variable uplift, high orographically driven rainfalls and multiple glaciations. Furthermore, lithology and erodibility vary substantially within the basin. Thanks to high-resolution geological, climatic and topographic data, the Rhône Basin is a suitable laboratory to explore these complexities. Elevation, relief, slope and hypsometric data as well as river profile information from digital elevation models are used to characterize the landscape's topography of around 50 tributary basins. Additionally, uplift over different timescales, glacial inheritance, precipitation patterns and erodibility of the underlying bedrock are quantified for each basin. Results show that the chosen topographic and controlling variables vary remarkably between different tributary basins. We investigate the link between observed topographic differences and the possible controlling variables through statistical analyses. Variations of elevation, slope and relief seem to be linked to differences in long-term uplift rate, whereas elevation distributions (hypsometry) and river profile shapes may be related to glacial imprint. This confirms that the landscape of the Rhône Basin has been highly preconditioned by (past) uplift and glaciation. Linear discriminant analyses (LDAs), however, suggest a stronger link between observed topographic variations and differences in erodibility. We therefore conclude that despite evident glacial and tectonic conditioning, a lithologic control is still preserved and measurable in the landscape of the Rhône tributary basins.

Influence of catchment quality and altitude on the water and sediment composition of 68 small lakes in Central Europe

68 lakes (63 Swiss, 2 French and 3 Italian) located in an altitudinal range between 334 and 2339m spanning a wide range of land-use have been investigated. The aim of the study was to discuss influences of geographic location, vegetation and land-use in the catchment area on the water and sediment chemistry of small lakes. Detailed quantitative description of land-use, vegetation, and climate in the watershed of all lakes was established. Surface and bottom water samples collected from each lake were analyzed for major ions and nutrients. Correlations were interpreted using linear regression analysis. Chemical parameters of water and sediment reflect the characteristics of the catchment areas. All lakes were alkaline since they were situated on calcareous bedrock. Concentrations of nitrogen and phosphorus strongly increase with increasing agricultural land-use. Na and K, however, are positively correlated with the amount of urbanization within the catchment area. These elements as well as dissolved organic carbon (DOC), Mg, Ca, and alkalinity, increase when the catchment is urbanized or used for agriculture. Total nitrogen and organic carbon in the sediments decrease distinctly if large parts of the catchment consist of bare land. No correlations between sediment composition and maximum water depth or altitude of the lakes were found.¶Striking differences in the water compositions of lakes above and below approximately 700 m of altitude were observed. Concentrations of total nitrogen and nitrate, total phosphorus, DOC, Na, K, Mg, Ca, and alkalinity are distinctly higher in most lakes below 700 m than above, and the pH of the bottom waters of these lakes is generally lower. Estimates of total nitrogen concentrations, even in remote areas, indicate that precipitation is responsible for increased background concentrations. At lower altitudes nitrogen concentrations in lakes is explained by the nitrogen loaded rain from urban areas deposited on the catchment, and with high percentages of agricultural land-use in the watershed.

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.