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     

Geochemistry of igneous rocks associated with ultramafic–mafic-hosted Cu (Co, Ni, Au) VMS deposits from the Main Uralian Fault (Southern Urals, Russia)

Ultramafic–mafic- and ultramafic-hosted Cu (Co, Ni, Au) volcanogenic massive sulfide (VMS) deposits from ophiolite complexes of the Main Uralian Fault, Southern Urals, are associated with island arc-type igneous rocks. Trace element analyses show that these rocks are geochemically analogous to Early Devonian boninitic and island arc tholeiitic rocks found at the base of the adjacent Magnitogorsk volcanic arc system, while they are distinguished both from earlier, pre-subduction volcanic rocks and from later volcanic products that were erupted in progressively more internal arc settings. The correlation between the sulfide host-rocks and the earliest volcanic units of the Magnitogorsk arc suggests a connection between VMS formation and infant subduction-driven intraoceanic magmatism.     

Organic matter in eolian dusts over the Atlantic Ocean

The elemental and mineralogical composition and the microfossil and detritus content of particulate fallout from the lower troposphere over the Atlantic Ocean have been extensively documented in earlier work, and it was possible to ascribe terrigenous source areas to such fallout. A brief review of the organic geochemistry of eolian dust is also presented here. The lipids of eolian dusts sampled from the air mass over the eastern Atlantic from about 35°N to 30°S were analyzed here.These lipids consisted mainly of normal alkanes, carboxylic acids and alcohols. The n-alkanes were found to range from n-C₂₃ to n-C₃₅, with high CPI values and maximizing at n-C₂₇ in the North Atlantic, at n-C₂₉ in the equatorial Atlantic and at n-C₃₁ in the South Atlantic. The n-fatty acids had mostly bimodal distributions, ranging from n-C₁₂ to n-C₃₀ (high CPI), with maxima at n-C₁₆ and in the northern samples at n-C₂₄ and in the southern samples at n-C₂₆. The n-alcohols ranged from n-C₁₂ to n-C₃₂, with high CPI values and maxima mainly at n-C₂₈. The compositions of these lipids indicated that their terrigenous sources were comprised mainly of higher plant vegetation and desiccated lacustrine mud flats on the African continent.     

Optimizing a layered and laterally constrained 2D inversion of resistivity data using Broyden's update and 1D derivatives

Dense profile-oriented resistivity data allows for 2D and 3D inversions. However, huge amounts of data make it practically impossible to do full 2D or 3D inversions on a routine basis. Therefore, a number of approximations have been suggested over the years to speed up computations. We suggest using a combination of Broyden's update on the Jacobian matrix with derivatives calculated using a 1D formulation on a parameterized 2D model of locally 1D layered models. The approximations bring down the effective number of 2D forward responses to a minimum, which again gives us the ability to invert very large sections. Broyden's update is not as useful with a parameterized problem as is the case with a smooth minimum structure problem that has been the usual application. 1D derivatives, however, seem to be very effective when initiating a full 2D solution with Broyden's update. We compare the different methods using two different kinds of data on two synthetic models and on two field examples. The most effective and reliable optimization combines 1D derivatives with a full 2D solution and Broyden's update. When using Broyden's update the Jacobian matrix needs to be reset every once in a while. We do this whenever the difference in data residual from the previous iteration is less than 5%. This combined inversion method reduces the computation time approximately a factor of 3 without losing model resolution.     

Heat capacity of wadeite-type K2Si4O9 and the pressure-induced stable decomposition of K-feldspar

The heat capacity of synthetic, stoichiometric wadeite-type K₂Si₄O₉ has been measured by DSC in the 195≤T(K)≤598 range. Near the upper temperature limit of our data, the heat capacity observed by DSC agrees with that reported by Geisinger et al. (1987) based on a vibrational model of their infrared and Raman spectroscopic data. However, with decreasing temperature, the Cp observed by DSC is progressively higher than that predicted from the vibrational model, suggesting that the standard entropy of K₂Si₄O₉ is likely to be larger than 198.9 ± 4.0 J/K · mol computed from the spectroscopic data. A fit to the DSC data gave: Cp(T) = 499.13 (±1.87) − 4.35014 · 10³(±3.489 · 10¹) · T ^−0.5, with T in K and average absolute percent deviation of 0.37%. The room-temperature compressibilities of kalsilite and leucite, hitherto unknown, have been measured as well. The data, fitted to the Murnaghan equation of state, gave K _o = 58.6 GPa, K _o ^′ = 0.1 for kalsilite and K _o = 45 GPa, K _o ^′ = 5.7 for α-leucite. Apart from the above mentioned data on the properties of the individual phases, we have also obtained reaction-reversals on four equilibria in the system K₂O-Al₂O₃-SiO₂. The Bayesian method has been used simultaneously to process the properties of 13 phases and 15 reactions between them to derive an internally consistent thermodynamic dataset for the K₂O-Al₂O₃-SiO₂ ternary. The enthalpy of formation of K₂Si₄O₉ wadeite is in perfect agreement with its revised calorimetric value, the standard entropy is 232.1 ± 10.4 J/K · mol, ∼15% higher than that implied by vibrational modeling. The phase diagram, generated from our internally consistent thermodynamic dataset, shows that for all probable P-T trajectories in the subduction regime, the stable pressure-induced decomposition of K-feldspar will produce coesite + kalsilite rather than coesite + kyanite + K₂Si₄O₉ (cf. Urakawa et al. 1994). Received: 11 June 1997 / Accepted: 2 December 1997     

Sharp spatially constrained inversion with applications to transient electromagnetic data

Time‐domain electromagnetic data are conveniently inverted by using smoothly varying 1D models with fixed vertical discretization. The vertical smoothness of the obtained models stems from the application of Occam‐type regularization constraints, which are meant to address the ill‐posedness of the problem. An important side effect of such regularization, however, is that horizontal layer boundaries can no longer be accurately reproduced as the model is required to be smooth. This issue can be overcome by inverting for fewer layers with variable thicknesses; nevertheless, to decide on a particular and constant number of layers for the parameterization of a large survey inversion can be equally problematic. Here, we present a focusing regularization technique to obtain the best of both methodologies. The new focusing approach allows for accurate reconstruction of resistivity distributions using a fixed vertical discretization while preserving the capability to reproduce horizontal boundaries. The formulation is flexible and can be coupled with traditional lateral/spatial smoothness constraints in order to resolve interfaces in stratified soils with no additional hypothesis about the number of layers. The method relies on minimizing the number of layers of non‐vanishing resistivity gradient, instead of minimizing the norm of the model variation itself. This approach ensures that the results are consistent with the measured data while favouring, at the same time, the retrieval of horizontal abrupt changes. In addition, the focusing regularization can also be applied in the horizontal direction in order to promote the reconstruction of lateral boundaries such as faults. We present the theoretical framework of our regularization methodology and illustrate its capabilities by means of both synthetic and field data sets. We further demonstrate how the concept has been integrated in our existing spatially constrained inversion formalism and show its application to large‐scale time‐domain electromagnetic data inversions.     

Evaluation of recent hydro-climatic changes in four tributaries of the Niger River Basin (West Africa)

West Africa experienced severe drought during the 1970s and 1980s, posing a threat to water resources. A wetter climate more recently suggests recovery from the drought. The Mann-Kendall trend and Theil-Sen’s slope estimator were applied to detect probable trends in weather elements in four sub-basins of the Niger River Basin between 1970 and 2010. The cross-entropy method was used to detect breakpoints in rainfall and runoff, Spearman’s rank test for correlation between the two, and cross-correlation analysis for possible lags. Results showed an overall increase in rainfall and runoff and a decrease in sunshine duration. Spearman’s coefficients suggest significant (5%) moderate to strong rainfall–runoff correlation for three sub-basins. A significant lower runoff was observed around 1979, with a rainfall break around 1992, indicating possible cessation of the drought. Temperatures increased significantly, at 0.02–0.05°C year^-1, with a negative wind speed trend for most stations. Half of the stations exhibited an increase in potential evapotranspiration.

EDITOR M.C. Acreman

ASSOCIATE EDITOR Not assigned     

Organic geochemistry and pore water chemistry of sediments from Mangrove Lake,Bermuda

Mangrove Lake, Bermuda, is a small coastal, brackish-water lake that has accumulated 14 m of banded, gelatinous, sapropelic sediments in less than 10⁴ yr. Stratigraphic evidence indicates that Mangrove Lake's sedimentary environment has undergone three major depositional changes (peat, freshwater gel, brackish-water gel) as a result of sea level changes. The deposits were examined geochemically in an effort to delineate sedimentological and diagenetic changes. Gas and pore water studies include measurements of sulfides, ammonia, methane, nitrogen gas, calcium, magnesium, chloride, alkalinity, and pH. Results indicate that sulfate reduction is complete, and some evidence is presented for bacterial denitrification and metal sulfide precipitation. The organic-rich sapropel is predominantly algal in origin, composed mostly of carbohydrates and insoluble macromolecular organic matter called humin with minor amounts of proteins, lipids, and humic acids. Carbohydrates and proteins undergo hydrolysis with depth in the marine sapropel but tend to be preserved in the freshwater sapropel. The humin, which has a predominantly aliphatic structure, increases linearly with depth and composes the greatest fraction of the organic matter. Humic acids are minor components and are more like polysaccharides than typical marine humic acids. Fatty acid distributions reveal that the lipids are of an algal and/or terrestrial plant source. Normal alkanes with a total concentration of 75 ppm exhibit two distribution maxima. One is centered about n-C₂₂ with no odd/even predominance, suggestive of a degraded algal source. The other is centered at n-C₃₁ with a distinct odd/even predominance indicative of a vascular plant origin. Stratigraphic changes in the sediment correlate to observed changes in the gas and pore water chemistry and the organic geochemistry.