Hydrological system analysis and modelling of the Kara River basin (West Africa) using a lumped metric conceptual model

This paper discusses the analysis and modelling of the hydrological system of the basin of the Kara River, a transboundary river in Togo and Benin, as a necessary step towards sustainable water resources management. The methodological approach integrates the use of discharge parameters, flow duration curves and the lumped conceptual model IHACRES. A Sobol sensitivity analysis is performed and the model is calibrated by applying the shuffled complex evolution algorithm. Results show that discharge generation in three nested catchments of the basin is affected by landscape physical characteristics. The IHACRES model adequately simulates the rainfall–runoff dynamics in the basin with a mean modified Nash-Sutcliffe efficiency measure of 0.6. Modelling results indicate that parameters controlling rainfall transformation to effective rainfall are more sensitive than those routing the streamflow. This study provides insights into understanding the catchment’s hydrological system. Nevertheless, further investigations are required to better understand detailed runoff generation processes.

EDITOR M.C. Acreman; ASSOCIATE EDITOR N Verhoest     

High resolution regional climate model simulations for Germany: Part II—projected climate changes

The projected climate change signals of a five-member high resolution ensemble, based on two global climate models (GCMs: ECHAM5 and CCCma3) and two regional climate models (RCMs: CLM and WRF) are analysed in this paper (Part II of a two part paper). In Part I the performance of the models for the control period are presented. The RCMs use a two nest procedure over Europe and Germany with a final spatial resolution of 7 km to downscale the GCM simulations for the present (1971–2000) and future A1B scenario (2021–2050) time periods. The ensemble was extended by earlier simulations with the RCM REMO (driven by ECHAM5, two realisations) at a slightly coarser resolution. The climate change signals are evaluated and tested for significance for mean values and the seasonal cycles of temperature and precipitation, as well as for the intensity distribution of precipitation and the numbers of dry days and dry periods. All GCMs project a significant warming over Europe on seasonal and annual scales and the projected warming of the GCMs is retained in both nests of the RCMs, however, with added small variations. The mean warming over Germany of all ensemble members for the fine nest is in the range of 0.8 and 1.3 K with an average of 1.1 K. For mean annual precipitation the climate change signal varies in the range of ?2 to 9 % over Germany within the ensemble. Changes in the number of wet days are projected in the range of ±4 % on the annual scale for the future time period. For the probability distribution of precipitation intensity, a decrease of lower intensities and an increase of moderate and higher intensities is projected by most ensemble members. For the mean values, the results indicate that the projected temperature change signal is caused mainly by the GCM and its initial condition (realisation), with little impact from the RCM. For precipitation, in addition, the RCM affects the climate change signal significantly.     

Approximation of van Genuchten Parameter Ranges from Hydraulic Conductivity Data

The use of retention function and relative conductivity function is essential for the calculation of flow in a variably saturated media using the Richards equation. A widely used mathematical model for this is the Mualem-van Genuchten model which requires the shape parameters and . These, however, are difficult to obtain. When data is scarce, and are often taken from literature and may deviate largely from actual values. The current study presents a novel mathematical model for the approximation of $$ and $$ and for the further estimation of realistic value ranges, which may be used as parameter space, for example, for the calibration of a numerical model. The model was developed for cases where data is scarce and only values of saturated hydraulic conductivity are available. It is based on a large data set from literature and it is demonstrated that the model estimates mean values from that data set with a good accuracy. In order to show the applicability of the model, a second data set has been compiled anew (provided as Supporting Information). The model is incorporated into the current version of the freeware computer program HYPAGS, which enables easy usage.     

PGE fractionation in seafloor hydrothermal systems: examples from mafic- and ultramafic-hosted hydrothermal fields at the slow-spreading Mid-Atlantic Ridge

The distribution of platinum group elements (PGEs) in massive sulfides and hematite–magnetite±pyrite assemblages from the recently discovered basalt-hosted Turtle Pits hydrothermal field and in massive sulfides from the ultramafic-hosted Logatchev vent field both on the Mid-Atlantic Ridge was studied and compared to that from selected ancient volcanic-hosted massive sulfide (VHMS) deposits. Cu-rich samples from black smoker chimneys of both vent fields are enriched in Pd and Rh (Pd up to 227 ppb and Rh up to 149 ppb) when compared to hematite–magnetite-rich samples from Turtle Pits (Pd up to 10 ppb, Rh up to 1.9 ppb). A significant positive correlation was established between Cu and Rh in sulfide samples from Turtle Pits. PGE chondrite-normalized patterns (with a positive Rh anomaly and Pd and Au enrichment), Pd/Pt and Pd/Au ratios close to global MORB, and high values of Pd/Ir and Pt/Ir ratios indicate mafic source rock and seawater involvement in the hydrothermal system at Turtle Pits. Similarly shaped PGE chondrite-normalized patterns and high values of Pd/Pt and Pd/Ir ratios in Cu-rich sulfides at Logatchev likely reflect a similar mechanism of PGE enrichment but with involvement of ultramafic source rocks.     

Using subsurface metazoan fauna to indicate groundwater–surface water interactions in the Nakdong River floodplain, South Korea

Hydrological interactions between surface water and groundwater (GW) can be described using hydrochemical and biological methods. Surface water–groundwater interactions and their effects on groundwater invertebrate communities were studied in the Nakdong River floodplain in South Korea. Furthermore, the GW-Fauna-Index, a promising new index for assessing the strength of surface-water influence on groundwater, was tested. The influence of surface water on groundwater decreased with increasing depth and distance from the river. While hydrochemistry prevailingly reflected the origin of the waters in the study area (i.e. whether alluvial or from adjacent rock), faunal communities seemed to display an affinity to surface-water intrusion. Fauna reacted quickly to changes in hydrology, and temporal changes in faunal community structure were significantly linked to the hydrological situation in the floodplain. The metazoan faunal community and the GW-Fauna-Index allow a distinction between surface and subsurface waters with varying degrees of exchange. The results indicate that hydrological conditions are reflected by faunal assemblages on a high spatiotemporal resolution, and that surface-water intrusion can be estimated using the GW-Fauna-Index.     

Facies architecture of an isolated carbonate platform: tracing the cycles of the Latemàr (Middle Triassic, northern Italy)

The 720-m-thick succession of the Middle Triassic Latemàr Massif (Dolomites, Italy) was used to reconstruct the lagoonal facies architecture of a small atoll-like carbonate platform. Facies analysis of the lagoonal sediments yields a bathymetric interpretation of the lateral facies variations, which reflect a syndepositional palaeorelief. Based on tracing of lagoonal flooding surfaces, the metre-scale shallowing-upward cycles are interpreted to be of allocyclic origin. Short-term sea-level changes led to subaerial exposure of wide parts of the marginal zone, resulting in the development of a tepee belt of varying width. Occasional emergence of the entire lagoon produced lagoon-wide decimetre-thick red exposure horizons. The supratidal tepee belt in the backreef area represented the zone of maximum elevation, which circumscribed the sub- to peritidal lagoonal interior during most of the platform's development. This tepee rim, the subtidal reef and a sub- to peritidal transition zone in between stabilized the platform margin. The asymmetric width of facies belts within individual metre-scale cycles was caused by redistribution processes that reflect palaeowinds and storm paths from the present-day south and west. The overall succession shows stratigraphic changes on a scale of tens of metres from a basal subtidal unit, overlain by three tepee-rich intervals, separated by tepee-poor units composed of subtidal to peritidal facies. This stacking pattern reflects two third-order sequences during the late Anisian to early middle Ladinian.     

Current trends of space occupation by encrusting excavating sponges on Colombian coral reefs

Some sponges of the genus Cliona (Porifera, Hadromerida, Clionidae) simultaneously excavate and encrust calcareous substratum, competing aggressively for illuminated space with corals and other organisms. To interpret current trends of reef space occupation, the patterns of distribution and size of three Caribbean species were examined at San Andrés Island and Islas del Rosario in Colombia. While Cliona aprica was ubiquitous, C. caribbaea (= C. langae) preferred deep and protected reef zones, and C. tenuis shallow and wave‐exposed settings. In contrast to the effect on other excavating sponges, chronic exposure to raw sewage did not significantly increase the abundance of the studied sponges. Substratum occupation/availability ratios showed a positive tendency of the sponges toward certain coral skeletons, and a negative or neutral tendency toward calcareous rock, indicating that establishment may be easier on clean, recently dead coral than on older, heavily incrusted substratum. High relief generally limits sponge size to that of the illuminated portions of the substratum. A generally lower proportion of small individuals than of larger ones indicates currently low recruitment rates and low subsequent mortality. Successful events of higher recruitment seem to have occurred for C. tenuis. These are related to the massive acroporid coral die‐off in the early 1980s and to asexual dispersion during storms, resulting in a current 10% substratum cover. Reefs with high coral mortality were and/or are thus more susceptible to colonization and subsequent space occupation by these sponges, although relief may prevent space monopolization.     

Meio-epifaunal wood colonization in the vicinity of methane seeps

In deep-sea environments, plant remains of several origins are found, including branches, twigs, leaves, and wood pieces, among others. As most of the deep-sea bottoms are oligotrophic and nutrient-limited, plant remains provide an oasis of localized organic enrichment and a substrate for colonization. Sunken wood was suggested to play an important evolutionary role in the diversification of chemosynthetic ecosystems, possibly representing stepping stones for the colonization between vent and seep ecosystems. In order to understand colonization processes of the Pacific Costa Rican meio-epifaunal assemblages associated with sunken wood, a field experiment was conducted on Mound 12 (8°55.778′N, 84°18.730′W) at ~1,000 m water depth. Woodblocks were placed in four different habitats (Mussel beds, tube worms, near mussel beds, rubble bottoms), and different local environmental conditions (seepage-active and seepage-inactive sites). Seven experimental Douglas fir woodblocks (each 1,047 cm² in surface area) were deployed from the R/V Atlantis using the manned submersible Alvin in February 2009 and recovered after 10.5 months in January 2010. Sample processing and analyses led to a data set of abundance (total 9,951 individuals) and spatial distribution of nine meio-epifaunal higher taxa/groups. Meio-epifaunal densities on individual woodblocks ranged from 3 to 26 ind.10 cm². Copepods accounted for the highest abundances (75.1%), followed by nauplii larvae (11.7%) and nematodes (9.8%). The maximum number of individuals (26.3 ind.10 cm⁻²) was found in blocks placed in seepage-inactive areas (near active mussel beds) in contrast to 2.9 ind.10 cm⁻² in active areas (within a mussel bed). A hierarchical cluster analysis grouped blocks according to seepage activity and not to habitat, but tests of similarity showed no significant differences in higher taxon composition and abundances, probably owing either to substrate homogeneity or low sample size. Copepods were the most abundant representatives, suggesting that this group is one of the most successful in colonizing in the early stage of succession, in this case while hardwood substrates are not yet decomposed or bored by bivalves.     

Seagrass ecophysiology meets ecological genomics: are we ready?

In March 2011, the Ecophysiology and Genetics Working Groups of the European Science Foundation COST Action ES 0906, entitled Seagrass Productivity: From Genes to Ecosystem Management, organized an exploratory workshop entitled “Linking Ecophysiology and Ecogenomics in Seagrass Systems”. The goal of the workshop was to discuss how to integrate comparative gene expression studies with physiological processes such as photosynthetic performance, carbon and nitrogen utilization and environmental adaptation. The main questions discussed for integrative research related to mechanisms of carbon utilization, light requirements, temperature effects and natural variation in pH and ocean acidification. It was concluded that the seagrass research community is still in the nascent stages of linking ecophysiology with genomic responses, as carbon and nitrogen metabolism of seagrasses have not been sufficiently well studied. Likewise, genomic approaches have only been able to assign meaningful interpretations to a handful of differentially expressed genes. Nevertheless, the way forward has been established.     

Treatment of time derivative and calculation of flow when solving groundwater flow problems by Galerkin finite element methods

Two techniques connected with the use of the finite element Galerkin method for solving the linear parabolic differential equation describing unsteady groundwater flow in an anisotropic non-homogeneous aquifer are introduced. The first is a mode superposition technique for dealing with the time derivative which involves computing the smallest eigenvalues and associated eigenvectors of the matrices arising from the Galerkin method. It is shown how such a technique allows us to interpret the response of the groundwater level to input in terms of parallel linear reservoirs. It is further argued that if properly implemented, the technique will have computational advantages over standard finite difference methods, e.g. in the case when the input function is constant over relatively large time subintervals. The second is a technique based on so-called generalized flow formulae for calculating flow values across external or internal boundaries, posterior to obtaining the groundwater level values. The implementation of the technique in the case of linear triangular elements on an irregular grid is discussed. It is finally argued from simplified cases that, apart from guaranteeing a match with prescribed input, the technique may often be expected to give more accurate flow values than those obtained directly from the groundwater gradients.