Dynamique et modélisation des crues dans le bassin du Mono à Nangbéto (Togo/Bénin)

Résumé

Les inondations des dernières années en Afrique de l’Ouest rendent nécessaire la prévision hydrologique pour réduire la vulnérabilité des populations. Nous avons étudié la dynamique des crues dans le bassin du Mono en amont du barrage de Nangbéto (Togo et Bénin) par des méthodes statistiques et déterministes. L’augmentation des précipitations au cours de la décennie 1990 et la mise en eau du barrage de Nangbéto en 1988 ont modifié le régime du Mono, dont les débits d’étiages ont augmenté. Les débits entrants et sortants du barrage sont équivalents, même lors des crues. Les maximums annuels des débits entrants peuvent être modélisés par une distribution des valeurs extrêmes (GEV) de type Fréchet, tandis que ceux des débits sortants, suivent une loi GEV de type Weibull (bornée). Le mécanisme dominant de génération des crues dans le bassin est la montée progressive des eaux associée à l’augmentation de la saturation des sols. L’analyse des débits de crue entre 1988 et 2010 n’a pas montré d’augmentation de la magnitude ou de la fréquence des crues. Le modèle GR4J reproduit bien les débits de crues en période calage, mais ces critères de performance se dégradent lorsque l’on s’en éloigne. Dans ces conditions, l’utilisation de ce modèle pour évaluer les risques futurs de crues serait hasardeuse et il est nécessaire d’envisager d’autres types d’approches.     

The position of the lowest levels in the boundary layer of atmospheric circulation models

Abstract

The error associated with the position of the lowest wind level in atmospheric boundary‐layer modelling is studied in connection with vertical resolutions typical of parametrization schemes for atmospheric circulation models. The test case is the neutrally stratified steady‐state boundary layer. Finite‐difference and finite‐element schemes of two types are used: in one case the lowest wind level is centred with respect to the surface and the lowest internal level where the shear stress is calculated, in the other case the lowest wind level is set very close to the surface (5 m). It is found that schemes of the latter type underestimate the friction force at the lowest level and consequently overestimate the wind and the surface stress. This error is largest at low resolution since it is due to the uncentring of the lowest wind level with respect to the stress levels. The error in schemes of the former type is different, and is associated with the determination of the surface stress from a wind at a height that may exceed the extent of the surface layer. For all neutral cases, this error can be made small by adding a corrective term to the traditional logarithmic formulation. This paper shows that considerably more error is generated by uncentred differencing than by deepening the surface layer.     

Properties of unstable waves in the lower St Lawrence Estuary

Abstract

The lower St Lawrence Estuary is an interesting case amongst estuaries in that it is wide enough to accommodate the development of mesoscale unstable waves and eddies. These features are generated by the runoff‐driven jet along this body's south shore. We present data yielding estimates of the length, time and velocity scales of these unstable disturbances. To relate these quantities to the dynamics we employ a 2‐layer quasigeos‐trophic instability model featuring realistic lateral shear. All model runs show short time and length scales, e‐folding periods of less than 10 days and wavelengths less than 50 km.     

Experimental study of the K2ZrSi3O9 (wadeite)–K2TiSi3O9 and K2(Zr,Ti)Si3O9–phlogopite systems at 2–3 GPa

Experiments ranging from 2 to 3 GPa and 800 to 1300 °C and at 0.15 GPa and 770 °C were performed to investigate the stability and mutual solubility of the K₂ZrSi₃O₉ (wadeite) and K₂TiSi₃O₉ cyclosilicates under upper mantle conditions. The K₂ZrSi₃O₉–K₂TiSi₃O₉ join exhibits complete miscibility in the P–T interval investigated. With increasing degree of melting the solid solution becomes progressively enriched in Zr, indicating that K₂ZrSi₃O₉ is the more refractory end member. At 2 GPa, in the more complex K₂ZrSi₃O₉–K₂TiSi₃O₉–K₂Mg₆Al₂Si₆O₂₀(OH)₄ system, the presence of phlogopite clearly limits the extent of solid solution of the cyclosilicate to more Zr-rich compositions [Zr/(Zr + Ti) > 0.85], comparable to wadeite found in nature, with TiO₂ partitioning strongly into the coexisting mica and/or liquid. However, at 1200 °C, with increasing pressure from 2 to 3 GPa, the partitioning behaviour of TiO₂ changes in favour of the cyclosilicate, with Zr/(Zr + Ti) of the K₂(Zr,Ti)Si₃O₉ phase decreasing from ∼0.9 to ∼0.6. The variation in the Ti content of the coexisting phlogopite is related to its degree of melting to forsterite and liquid, following the major substitution ^VITi+^VI□=2^VIMg. Received: 26 January 1999 / Accepted: 10 January 2000     

Dense collection areas and terrestrial alteration of meteorites in the Atacama Desert

In the last 15 years, more than 2700 meteorites have been recovered and officially classified from the Atacama Desert. Although the number of meteorites collected in the Atacama has risen, the physical and climatic properties of the dense collection areas (DCAs) have not been fully characterized. In this article, we compiled the published data of all classified meteorites found in the Atacama Desert to (i) describe the distribution by meteorite groups, (ii) compare the weathering degree of chondrites among different Atacama DCAs and other hot and cold deserts, and (iii) determine the preservation conditions of chondrites in the main Atacama DCAs in relation with the local climatic conditions. The 35 DCAs so far identified in the Atacama Desert are located in three main morphotectonic units: The Coastal Range (CR), Central Depression (CD), and Pre-Andean Range/Basement. A comparison with reported weathering data from other cold and hot deserts indicates that the mean terrestrial weathering of Atacama chondrites (W1–2), displays less alteration than other hot deserts (W2–3) and resembles the weathering distribution of the Antarctic meteorites (W1–2). The highest abundance of Atacama chondrites with low weathering (≤W2) is localized in the CD (78.8%, N = 1435), which is protected from the coastal fog influence and seasonal rainfalls and displays the oldest surfaces in the Atacama Desert. The morphogenetic classification based on present-day temperatures and precipitations of the main Atacama DCAs reveals similar regional/subregional climatic conditions in the most productive areas and a truly productive surface for meteorite recovery between 5% and 58% of the quadrangles formally defined for each Atacama DCA. Our morphogenetic classification lacks consideration of some meteorological parameters such as the coastal fog, so it cannot fully explain the differences in weathering patterns among CR chondrites. Future studies of chondrite preservation in the Atacama DCAs should consider other meteorological variables such as relative humidity, specific humidity, or dew point, in combination with exposure ages of meteorites and its surfaces.     

Adaptation of the vertical resolution in the mixed layer for HYCOM

This paper focuses on the dynamics of the mixed layer. When the mixed layer depth increases, the vertical discretisation eventually becomes too sparse at the bottom of this layer to accurately resolve its evolution and strong numerical errors can appear. This is linked to the fact that the vertical resolution is concentrated in the upper part of the ocean and does not adapt to the deepening of the mixed layer.Knowing that the HYbrid Coordinate Ocean Model (HYCOM) is able to modify the distribution of the vertical levels, we suggest in this paper a method to adapt the resolution to the mixed layer extension. This method is tested in 1-D configurations for two academic atmospheric forcing conditions (strong convection and wind-mixing) and a realistic forcing extending over one year, with seasonal restratification following strong winter convection. The new method improves the results in all cases, and in particular when the mixed layer reaches deep layers.     

Soil cracking effects on hydrological and erosive processes: a study case in Mediterranean cultivated vertisols

Shrink–swell soils, such as those in a Mediterranean climate regime, can cause changes in terms of hydrological and erosive responses because of the changing soil water storage conditions. Only a limited number of long‐term studies have focused on the impacts on both hydrological and erosive responses and their interactions in an agricultural environment. In this context, this study aims to document the dynamics of cracks, runoff and soil erosion within a small Mediterranean cultivated catchment and to quantify the influence of crack processes on the water and sediment supplied to a reservoir located at the catchment outlet using water and sediment measurements at a cultivated field outlet as baseline. Detailed monitoring of the presence of topsoil cracks was conducted within the Kamech catchment (ORE OMERE, Tunisia), and runoff and suspended sediment loads were continuously measured over a long period of time (2005–2012) at the outlets of a field (1.32 ha) and a catchment (263 ha). Analysis of the data showed that topsoil cracks were open approximately half of the year and that the rainfall regime and water table level conditions locally control the seasonal cracking dynamics. Topsoil cracks appeared to seriously affect the generation of runoff and sediment concentrations and, consequently, sediment yields, with similar dynamics observed at the field and catchment outlets. A similar time lag in the seasonality between water and sediment delivery was observed at these two scales: although the runoff rates were globally low during the presence of topsoil cracks, most sediment transport occurred during this period associated with very high sediment concentrations. This study underlines the importance of a good prediction of runoff during the presence of cracks for reservoir siltation considerations. In this context, the prediction of cracking effects on runoff and soil erosion is a key factor for the development of effective soil and water management strategies and downstream reservoir preservation. Copyright © 2016 John Wiley & Sons, Ltd.     

Displacements and transformations of nitrate-rich and nitrate-poor water masses in the tropical Pacific during the 1997 El Niño

A Lagrangian analysis was applied to the outputs of a coupled physical-biogeochemical model to describe the redistribution of nitrate-rich and nitrate-poor surface water masses in the tropical Pacific throughout the major 1997 El Niño. The same tool was used to analyze the causes of nitrate changes along trajectories and to investigate the consequences of the slow nitrate uptake in the high nutrient low chlorophyll (HNLC) region during the growth phase of the event. Three patterns were identified during the drift of water masses. The first mechanism is well known along the equator: oligotrophic waters from the western Pacific are advected eastward and retain their oligotrophic properties along their drift. The second concerns the persistent upwelling in the eastern basin. Water parcels have complex trajectories within this retention zone and remain mesotrophic. This study draws attention to the third process which is very specific to the HNLC region and to the El Niño period. During the 1997 El Niño, horizontal and vertical inputs of nitrate decreased so dramatically that nitrate uptake by phytoplankton became the only mechanism driving nitrate changes along pathways. The study shows that because of the slow nitrate uptake characteristic of the tropical Pacific HNLC system, nitrate in the pre-El Niño photic layer can support biological production for a period of several months. As a consequence, the slow nitrate uptake delays the gradual onset of oligotrophic conditions over nearly all the area usually occupied by upwelled waters. Owing to this process, mesotrophic conditions persist in the tropical Pacific during El Niño events.     

An automated salt‐tracing gauge for flow‐velocity measurement

This article introduces the SVG (salt‐velocity gauge), a novel automated technique for measuring flow velocity by means of salt tracing. SVG allows a high measuring rate (up to one every 2 seconds), short control section length (down to 10 cm), high accuracy (+[sol ]?1·5 cm s^?1), and unbiased calculation of the mean velocity in experimental conditions with turbulent, supercritical flow. A few cubic centimetres of saturated salt solution (NaCl) are injected into the flow at regular time intervals using a programmable solenoid valve. The tracer successively passes two conductivity probes placed a short distance downstream. The transformation of the signal between the two probes is modelled as a one‐dimensional diffusion wave equation. Model calibration gives an estimation of the mean velocity and the diffusion for each salt plume. Two implementations of the SVG technique are described. The first was an outdoors simulated rainfall experiment in Senegal (conductivity probes at 40 cm apart, 8 Hz measurement rate, salt injections at 10 second intervals). Mean velocity was estimated to range between 0·1 and 0·3 m s^?1. The second was a laboratory‐based flume experiment (conductivity probes at 10 cm apart, 32 Hz, salt injections at 2 second intervals). Another SVG with probes at 34 cm apart was used for comparison. An acoustic Doppler velocimeter (ADV) was also used to give an independent assessment of velocity. Using the 10 cm salt gauge, estimated mean velocity ranged from 0·6 to 0·9 m s^?1 with a standard deviation of 1·5 cm s^?1. Comparisons between ADV, 10 cm SVG and 34 cm SVG were consistent and demonstrated that the salt‐tracing results were unbiased and independent of distance between probes. Most peaks were modelled with r² > 90 per cent. The SVG technology offers an alternative to the dye‐tracing technique, which has been severely criticized in the literature because of the wide interval of recommended values for the correction factor α to be applied to the timings. This article demonstrates that a fixed value of α is inappropriate, since the correction factor varies with velocity, diffusion and the length of the control section. Copyright © 2005 John Wiley & Sons, Ltd.