Assessing uncertainties in a conceptual water balance model using Bayesian methodology / Estimation bayésienne des incertitudes au sein d’une modélisation conceptuelle de bilan hydrologique

Abstract

Abstract The aim of this study was to estimate the uncertainties in the streamflow simulated by a rainfall–runoff model. Two sources of uncertainties in hydrological modelling were considered: the uncertainties in model parameters and those in model structure. The uncertainties were calculated by Bayesian statistics, and the Metropolis-Hastings algorithm was used to simulate the posterior parameter distribution. The parameter uncertainty calculated by the Metropolis-Hastings algorithm was compared to maximum likelihood estimates which assume that both the parameters and model residuals are normally distributed. The study was performed using the model WASMOD on 25 basins in central Sweden. Confidence intervals in the simulated discharge due to the parameter uncertainty and the total uncertainty were calculated. The results indicate that (a) the Metropolis-Hastings algorithm and the maximum likelihood method give almost identical estimates concerning the parameter uncertainty, and (b) the uncertainties in the simulated streamflow due to the parameter uncertainty are less important than uncertainties originating from other sources for this simple model with fewer parameters.     

Flood forecasting and infiltration modeling/Prévision de crue et modélisation de l’infiltration

Abstract

Abstract The role of accuracy in the representation of infiltration on the effectiveness of real-time flood forecasting models was investigated. A simple semi-distributed model of conceptual type with adaptive estimate of hydraulic characteristics included in the infiltration component was selected. Infiltration was described by a very accurate approach recently formulated for complex rainfall patterns, or alternatively through a simpler formulation known as an extension of the classical time compression approximation. The results indicated that, for situations involving a significant rainfall variability in space, the inaccuracy in the representation of infiltration cannot be corrected by the adaptive component of the rainfall–runoff model. A preliminary analysis of the role of an approximation of saturated hydraulic conductivity to be used in each homogeneous area of the semi-distributed model used both in non-adaptive version and in real-time is also presented.     

Identification of runoff generation processes using combined hydrometric,tracer and geophysical methods in a headwater catchment in South Africa / Identification des processus de génération de l'écoulement par combinaison de méthodes hydrométriques,de traçage et géophysiques dans un bassin versant sud-africain

Résumé

Des mesures hydrométriques classiques et l'imagerie 2-D de résistivité électrique (IRE) détaillée ont été combinées avec le prélèvement de traceurs pour identifier les processus hydrologiques en jeu dans un bassin versant semi-aride de la province Cap-Oriental, en Afrique du Sud. L'étude des événements de précipitation et d'écoulement a souligné la forte relation qui existe entre les précipitations et les caractéristiques de la génération des écoulements. L'observation du potentiel hydrique du sol ainsi que du niveau d'eau souterraine a mis en évidence le développement d'une nappe perchée dans le sol. Ces résultats sont confirmés par des séparations d'hydrogrammes à base de traceurs et illustrent le rôle important des écoulements souterrains superficiels. L'étude par IRE a permis de mieux connaître la structure du sous-sol. Finalement, l'étude par IRE, combinée avec des mesures obtenues par réflectométrie en domaine temporel (TDR), a permis d'extrapoler à des mesures sélectives de teneur en eau. Pour récapituler, l'utilisation et la combinaison de différentes méthodes de terrain ont mené au développement d'un modèle conceptuel du fonctionnement hydrologique de ce bassin versant. Le rôle dominant des mécanismes de subsurface a été démontré.     

The last deglaciation in the Picos de Europa National Park (Cantabrian Mountains,northern Spain)

A sedimentological and geochemical study of the Lago Enol sequence (Cantabrian Mountains, northern Spain), together with detailed geomorphological mapping, provides a first record of glacier evolution and climate change over the last 40 ka in the Picos de Europa National Park. The Enol glacier retreated from its maximum extent prior to 40 ka BP as demonstrated by the onset of proglacial lacustrine sedimentation in two glaciated depressions: the Comella hollow to the north (before 40 ka BP) and the Lago Enol (before 38 ka BP). These results support previous evidence that the maximum extent of southern European glaciers occurred earlier than in northern Europe. Alternation of homogeneous and laminated proglacial sediments during the glacier retreat illustrate a dynamic glacial evolution during the Marine Isotope Stage (MIS) 3 (40–26 ka BP). A slight warming is detected at 26 ka ago with the change from proglacial sediments (in a lake located in contact to the glacier) to glaciolacustrine sedimentation (in a non‐contact or distal lake). Finally, the onset of organic‐rich sediments took place at 18 ka ago. This last transition occurred in two phases, similarly to the North Atlantic Last Termination, suggesting a link between North Atlantic Deep Water formation oscillations and palaeohydrological variability in the Cantabrian Mountains. Copyright © 2009 John Wiley & Sons, Ltd.     

(Re) deposition of loess in southern Limbourg,the Netherlands 3. Field evidence for conditions of deposition of the middle and upper silt loam complexes,and landscape evolution at Nagelbeek

This paper is about conditions and processes during the accumulation of the Upper Pleniglacial Middle and Upper Silt Loam complexes at Nagelbeek, Limbourg, a Weichselian stratigraphical type locality. The complexes correlate with Hesbayan and Brabantian loesses in Belgium, formed between 28,000 and 12,400 yrs BP. They are characterized by sediment features (cracks, mass deformations, and cut-and-fill features) that are ascribed to a cryogenic regimen. The main Middle Silt Loam subunit contains contraction cracks deformed and further modified by masswasting. Intra-Upper-Pleniglacial pedogenic intervals are inferred from three greyish bands, from a truncated weathered zone surmounting the Middle Silt Loam, and from an aeolian-cumulic palaeosol at the base of the Upper Silt Loam. Latter two zones and an intervening erosional unconformity were cryogenically reworked, less than 22,000 yrs BP, to form a complex stratigraphical marker-zone: the ‘Horizon à Langues de Nagelbeek’. Morphological variants of the marker-zone are described. The geomorphic evolution of the Nagelbeek upland, from Saalian to Holocene times, is summarized in unconformity maps and a time diagram.     

Lithosphere formation in the central Slave Craton (Canada): plume subcretion or lithosphere accretion?

Major-element compositions of minerals in peridotite xenoliths from the Lac de Gras kimberlites provide constraints on the mode of lithosphere formation beneath the central Slave Craton, Canada. Magnesia contents of reconstructed whole rocks correlate positively with NiO and negatively with CaO contents, consistent with variable partial melt extraction. Alumina and Cr₂O₃ contents are broadly positively correlated, suggestive of melt depletion in the absence of a Cr–Al phase. Garnet modes are high at a given Al₂O₃ content (a proxy for melt depletion), falling about a 7 GPa melt depletion model. These observations, combined with high olivine Mg# and major-element relationships of FeO-poor peridotites (<7.5 wt%) indicative of melt loss at pressures >3 GPa (residual FeO content being a sensitive indicator of melt extraction pressure), and similar high pressures of last equilibration (∼4.2 to 5.8 GPa), provide multiple lines of evidence that the mantle beneath the central Slave Craton has originated as a residue from high-pressure melting, possibly during plume subcretion. Apparent low melt depletion pressures for high-FeO peridotites (>7.5 wt%) could suggest formation in an oceanic setting, followed by subduction to their depth of entrainment. However, these rocks, which are characterised by low SiO₂ contents (<43 wt%), are more likely to be the result of post-melting FeO-addition, leading to spuriously low estimates of melt extraction pressures. They may have reacted with a silica-undersaturated melt that dissolved orthopyroxene, or experienced olivine injection by crystallising melts. A secular FeO-enrichment of parts of the deep mantle lithosphere is supported by lower average Mg# in xenolithic olivine (91.7) compared to olivine inclusions in diamond (92.6).     

The November 2002 eruption at Piton de la Fournaise volcano,La Réunion Island: ground deformation,seismicity, and pit crater collapse

An eruption on the eastern flank of Piton de la Fournaise volcano started on 16 November, 2002 after 10 months of quiescence. After a relatively constant level of activity during the first 13 days of the eruption, lava discharge, volcanic tremor and seismicity increased from 29 November to 3 December. Lava effusion suddenly ceased on 3 December while shallow earthquakes beneath the Dolomieu summit crater were still recorded at a rate of about one per minute. This unusual activity continued and increased in intensity over the next three weeks, ending with the formation of a pit crater within Dolomieu. Based on ground deformation, measured by rapid-static and continuous GPS and an extensometer, seismic data, and lava effusion patterns, the eruptive period is divided into five stages: 1) slow summit inflation and sporadic seismicity; 2) rapid summit inflation and a short seismic crisis; 3) rapid flank inflation, onset of summit deflation, sporadic seismicity, accompanied by stable effusion; 4) flank inflation, coupled with summit deflation, intense seismicity, and increased lava effusion; and finally 5) little deflation, intense shallow seismicity, and the end of lava effusion. We propose a model in which the pre-intrusive inflation of Stage 1 in the months preceding the eruption was caused by a magma body located near sea level. The magma reservoir was the source of an intrusion rising under the summit during Stage 2. In Stage 3, the magma ponded at a shallow level in the edifice while the lateral injection of a radial dike reached the surface on the eastern flank of the basaltic volcano, causing lava effusion. Pressure decrease in the magmatic plumbing system followed, resulting in upward migration of a collapse front, forming a subterranean column of debris by faulting and stoping. This caused intense shallow seismicity, increase in discharge of lava and volcanic tremor at the lateral vent in Stage 4 and, eventually the formation of a pit crater in Stage 5.     

La prévision du climat : de l'échelle saisonnière à l'échelle décennale

The atmosphere and the ocean are subject to many dynamical instabilities, which limit the time during which their behaviour can be deterministically forecasted. At longer timescales, the atmosphere can be predicted at best using statistical methods, as a response to external forcing linked to sea- and land-surface anomalies. Climate being defined as the mean of atmospheric states, it appears that it can be predicted up to a few months in advance, which is the characteristic time of the so-called slow components of the climate system. Forecasting can sometimes be extended to longer time ranges, especially when the coupled ocean–atmosphere system exhibits internal variability modes, with characteristic times of a few years. Seasonal climate forecasting is most often based upon Monte-Carlo simulations, where the various realisations correspond to slightly different initial conditions. The present sate-of-the-art in Europe (ECMWF) and/or in the USA (IRI) allows to forecast such major phenomena, as El Niño, up to six months in advance. Finally, some parameters may exhibit predictability at still longer time-ranges (inter-annual to decadal), but only for certain regions. The example of electricity production is used to underline the potentially large economical benefit of seasonal climate forecasting. To cite this article: J.-C. André et al., C. R. Geoscience 334 (2002) 1115–1127.

Résumé

L'atmosphère et l'océan sont le siège d'instabilités dynamiques, qui limitent la durée pendant laquelle il est possible d'en prévoir l'évolution de façon déterministe. Au-delà, l'atmosphère n'est plus prévisible, au mieux, que de façon statistique, en fonction du forçage externe qu'exerce(nt) sur elle l'océan et/ou la surface des continents. Le climat (au sens d'une moyenne des états atmosphériques) se révèle ainsi prévisible jusqu'à des échéances temporelles de quelques mois, échelle de temps caractéristique des composantes dites « lentes » du système climatique. La prévision peut s'étendre à des échéances parfois plus longues, dans le cas où le système couplé océan–atmosphère posséderait des modes de variabilité temporelle de périodes caractéristiques de quelques années. La prévision climatique saisonnière est très souvent construite à partir de simulations de type Monte-Carlo, avec des ensembles de réalisations utilisant des conditions initiales légèrement différentes. Dans l'état actuel de ces prévisions, qu'elles soient réalisées en Europe (CEPMMT) ou aux États-Unis (IRI), il est possible de prévoir environ six mois à l'avance un certain nombre de phénomènes climatiques, en particulier ceux liés aux épisodes dits « El Niño », pour lesquels l'amplitude des variations est suffisamment importante. Il existe, par ailleurs, une prévisibilité à encore plus longue échéance (inter-annuelle à décennale), mais seulement pour certains paramètres et certaines régions. L'exemple de la production d'électricité montre l'importance économique potentielle très grande de la prévision climatique saisonnière. Pour citer cet article : J.-C. André et al., C. R. Geoscience 334 (2002) 1115–1127.     

Une dépression piézométrique naturelle en hausse au Sahel (Sud-Ouest du Niger)A rising piezometric depression in the Sahel (southwestern Niger)

In southwest Niger, the Continental Terminal water table displays a natural hollow shape about 10 m in depth over an area of 4000 km². A 10-year survey of this hollow aquifer has shown that current recharge is above $\text{20 }\text{mm}\text{yr}{}^{\mathrm{?1}}$. The water table has risen continuously since the 1950–1960s as a result of land clearance. This shows a disequilibrium in the aquifer balance. The long-term recharge rate is estimated by radioisotopes to be around $\text{1 }\text{mm}\text{yr}{}^{\mathrm{?1}}$. This figure fits with the only possible origin of the piezometric depression, i.e. evapotranspiration losses in its centre. To cite this article: G. Favreau et al., C. R. Geoscience 334 (2002) 395–401.     

Research on Autonomous Orbit Determination of Navigation Satellite Based on Crosslink Range and Orientation Parameters Constraining

Introduction Thepotentialvulnerabilityofsatellitenaviga tionsystemthatreliesongroundstationsisthat thesystemwouldbreakdownifgroundstations weredestroyed,whichcannotmeettherequire mentofnavigationwarfare[1].Withthedevelop mentofsuchspace basedsystemsasgrou…