High-resolution numerical modelling of flow-vegetation interactions

In this paper, we present and apply a new three-dimensional model for the prediction of canopy-flow and turbulence dynamics in open-channel flow. The approach uses a dynamic immersed boundary technique that is coupled in a sequentially staggered manner to a large eddy simulation. Two different biomechanical models are developed depending on whether the vegetation is dominated by bending or tensile forces. For bending plants, a model structured on the Euler-Bernoulli beam equation has been developed, whilst for tensile plants, an N-pendula model has been developed. Validation against flume data shows good agreement and demonstrates that for a given stem density, the models are able to simulate the extraction of energy from the mean flow at the stem-scale which leads to the drag discontinuity and associated mixing layer.

Geodynamic reconstructions of the South America-Antarctica plate system

The South America-Antarctica plate system shows many oceanic accretionary systems and subduction zones that initiated and then stopped. To better apprehend the evolution of the system, geodynamic reconstructions (global) have been created from Jurassic (165 Ma) to present, following the techniques used at the University of Lausanne. However, additional synthetic magnetic anomalies were used to refine the geodynamics between 33 Ma and present. The reconstructions show the break up of Gondwana with oceanisation between South America (SAM) and Antarctica (ANT), together with the break off of `Andean' geodynamical units (GDUs). We propose that oceanisation occurs also east and south of the Scotian GDUs. Andean GDUs collide with other GDUs crossing the Pacific. The west coast of SAM and ANT undergo a subsequent collision with all those GDUs between 103 Ma and 84 Ma, and the Antarctic Peninsula also collides with Tierra del Fuego. The SAM-ANT plate boundary experienced a series of extension and shortening with large strike-slip component, culminating with intra-oceanic subduction leading to the presence of the `V-' and anomalies in the Weddell Sea. From 84 Ma, a transpressive collision takes place in the Scotia region, with active margin to the east. As subduction propagates northwards into an old and dense oceanic crust, slab roll-back initiates, giving rise to the western Scotia Sea and the Powell Basin opening. The Drake Passage opens. As the Scotian GDUs migrate eastwards, there is enough space for them to spread and allow a north-south divergence with a spreading axis acting simultaneously with the western Scotia ridge. Discovery Bank stops the migration of South Orkney and `collides with' the SAM-ANT spreading axis, while the northern Scotian GDUs are blocked against the Falkland Plateau and the North-East Georgia Rise. The western and central Scotia and the Powell Basin spreading axes must cease, and the ridge jumps to create the South Sandwich Islands Sea. The Tierra del Fuego-Patagonia region has always experienced mid-oceanic ridge subduction since 84 Ma. Slab window location is also presented (57-0 Ma), because of its important implication for heat flux and magmatism. (C) 2011 Elsevier Ltd. All rights reserved.

Relating the permeability of quartz sands to their grain size and spectral induced polarization characteristics

Recently, Revil & Florsch proposed a novel mechanistic model based on the polarization of the Stern layer relating the permeability of granular media to their spectral induced polarization (SIP) characteristics based on the formation of polarized cells around individual grains. To explore the practical validity of this model, we compare it to pertinent laboratory measurements on samples of quartz sands with a wide range of granulometric characteristics. In particular, we measure the hydraulic and SIP characteristics of all samples both in their loose, non-compacted and compacted states, which might allow for the detection of polarization processes that are independent of the grain size. We first verify the underlying grain size/permeability relationship upon which the model of Revil & Florsch is based and then proceed to compare the observed and predicted permeability values for our samples by substituting the grain size characteristics by corresponding SIP parameters, notably the so-called Cole-Cole time constant. In doing so, we also asses the quantitative impact of an observed shift in the Cole-Cole time constant related to textural variations in the samples and observe that changes related to the compaction of the samples are not relevant for the corresponding permeability predictions. We find that the proposed model does indeed provide an adequate prediction of the overall trend of the observed permeability values, but underestimates their actual values by approximately one order-of-magnitude. This discrepancy in turn points to the potential importance of phenomena, which are currently not accounted for in the model and which tend to reduce the characteristic size of the prevailing polarization cells compared to the considered model, such as, for example, membrane polarization, contacts of double-layers of neighbouring grains, and incorrect estimation of the size of the polarized cells because of the irregularity of natural sand grains.

Ambient vertical flow in long-screen wells: a case study in the Fontainebleau Sands Aquifer (France)

A tritium (H-3) profile was constructed in a long-screened well (LSW) of the Fontainebleau Sands Aquifer (France), and the data were combined with temperature logs to gain insight into the potential effects of the ambient vertical flow (AVF) of water through the well on the natural aquifer stratification. AVF is commonly taken into account in wells located in fracture aquifers or intercepting two different aquifers with distinct hydraulic heads. However, due to the vertical hydraulic gradient of the flow lines intercepted by wells, AVF of groundwater is a common process within any type of aquifer. The detection of 3H in the deeper parts of the studied well ( approximate depth 50m), where H-3-free groundwater is expected, indicates that shallow young water is being transported downwards through the well itself. The temperature logs show a nearly zero gradient with depth, far below the mean geothermal gradient in sedimentary basins. The results show that the age distribution of groundwater samples might be biased in relation to the age distribution in the surroundings of the well. The use of environmental tracers to investigate aquifer properties, particularly in LSWs, is then limited by the effects of the AVF of water that naturally occurs through the well.

Urinary sphincter based on electronics and artificial muscles

Objectives: The AMS 800TM is the current artificial urinary sphincter (AUS) for incontinence due to intrinsic sphincter deficiency. Despite good clinical results, technical failures inherent to the hydraulic mechanism or urethral ischemic injury contribute to revisions up to 60%. We are developing an electronic AUS, called ARTUS to overcome the rigors of AMS. The objective of this study was to evaluate the technical efficacy and tissue tolerance of the ARTUS system in an animal model.Methods: The ARTUS is composed by three parts: the contractile unit, a series of rings and an integrated microprocessor. The contractile unit is made of Nitinol fibers. The rings are placed around the urethra to control the flow of urine by squeezing the urethra. They work in a sequential alternative mode and are controlled by a microprocessor. In the first phase a three-rings device was used while in the second phase a two-rings ARTUS was used. The device was implanted in 14 sheep divided in two groups of six and eight animals for study purpose. The first group aimed at bladder leak point pressure (BLPP) measurement and validation of the animal model; the second group aimed at verifying mid-term tissue tolerance by explants at twelve weeks. General animal tolerance was also evaluated.Results: The ARTUS system implantation was uneventful. When the system was activated, the BLPP was measured at 1.038±0.044 bar (mean±SD). Urethral tissue analysis did not show significant morphological changes. No infection and no sign of discomfort were noted in animals at 12 weeks.Conclusions: The ARTUS proved to be effective in continence achievement in this study. Histological results support our idea that a sequential alternative mode can avoid urethral atrophy and ischemia. Further technical developments are needed to verify long-term outcome and permit human use.

Canvis ambientals i relacions tròfiques en rius mediterranis

Els rius i rieres mediterranis són ecosistemes que es caracteritzen per fortes oscil•lacions de cabal i temperatura al llarg de l’any. Aquestes oscil•lacions provoquen canvis ambientals en l'hàbitat i en els recursos que afecten directament o indirecta la biota que habita aquests ecosistemes, la qual, per tant, ha de presentar adaptacions a aquestes oscil•lacions ambientals. L'escenari actual de canvi climàtic preveu una intensificació dels fenòmens de sequera i augment de temperatura. Entendre com la biota dels rius respon a aquestes fluctuacions és de gran importància per poder anticipar les respostes d'aquests sistemes als imminents canvis ambientals així com per gestionar adequadament els recursos hídrics en un futur. Els objectius principals d'aquesta tesi eren: caracteritzar estructural i funcionalment dues rieres intermitents mediterrànies al llarg dels diferents períodes característics del cicle anual i veure els efectes d'un augment de la sequera; veure com aquests efectes podien afectar l'ecosistema ripari circumdant i establir com diferències en la qualitat de la matèria orgànica derivades del canvi climàtic pot afectar el fitness i desenvolupament dels invertebrats. Aquests objectius s'han pogut complir només parcialment, ja que adversitats climàtiques van impedir finalitzar amb èxit la manipulació del cabal al camp i la resolució d'algunes dades no ha estat prou bona com per aplicar els models corresponents. Aquests contratemps s'han solucionat amb la incorporació de dos nous experiments (un encara s'ha de realitzar), fet que ha fet enlentir la finalització de la tesi.

Multiple-scale hydrological and geophysical data integration through non-linear Bayesian sequential simulation

Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale. However, extending the corresponding approaches to the scale of a field site represents a major, and as-of-yet largely unresolved, challenge. To address this problem, we have developed downscaling procedure based on a non-linear Bayesian sequential simulation approach. The main objective of this algorithm is to estimate the value of the sparsely sampled hydraulic conductivity at non-sampled locations based on its relation to the electrical conductivity logged at collocated wells and surface resistivity measurements, which are available throughout the studied site. The in situ relationship between the hydraulic and electrical conductivities is described through a non-parametric multivariatekernel density function. Then a stochastic integration of low-resolution, large-scale electrical resistivity tomography (ERT) data in combination with high-resolution, local-scale downhole measurements of the hydraulic and electrical conductivities is applied. The overall viability of this downscaling approach is tested and validated by comparing flow and transport simulation through the original and the upscaled hydraulic conductivity fields. Our results indicate that the proposed procedure allows obtaining remarkably faithful estimates of the regional-scale hydraulic conductivity structure and correspondingly reliable predictions of the transport characteristics over relatively long distances.

Maps as risk mitigation tools. Adaptation of the Swiss hazard assessment and mapping methodology to a Moroccan site: Beni Mellal

Dans le contexte climatique actuel, les régions méditerranéennes connaissent une intensification des phénomènes hydrométéorologiques extrêmes. Au Maroc, le risque lié aux inondations est devenu problématique, les communautés étant vulnérables aux événements extrêmes. En effet, le développement économique et urbain rapide et mal maîtrisé augmente l'exposition aux phénomènes extrêmes. La Direction du Développement et de la Coopération suisse (DDC) s'implique activement dans la réduction des risques naturels au Maroc. La cartographie des dangers et son intégration dans l'aménagement du territoire représentent une méthode efficace afin de réduire la vulnérabilité spatiale. Ainsi, la DDC a mandaté ce projet d'adaptation de la méthode suisse de cartographie des dangers à un cas d'étude marocain (la ville de Beni Mellal, région de Tadla-Azilal, Maroc). La méthode suisse a été adaptée aux contraintes spécifiques du terrain (environnement semi-aride, morphologie de piémont) et au contexte de transfert de connaissances (caractéristiques socio-économiques et pratiques). Une carte des phénomènes d'inondations a été produite. Elle contient les témoins morphologiques et les éléments anthropiques pertinents pour le développement et l'aggravation des inondations. La modélisation de la relation pluie-débit pour des événements de référence, et le routage des hydrogrammes de crue ainsi obtenus ont permis d'estimer quantitativement l'aléa inondation. Des données obtenues sur le terrain (estimations de débit, extension de crues connues) ont permis de vérifier les résultats des modèles. Des cartes d'intensité et de probabilité ont été obtenues. Enfin, une carte indicative du danger d'inondation a été produite sur la base de la matrice suisse du danger qui croise l'intensité et la probabilité d'occurrence d'un événement pour obtenir des degrés de danger assignables au territoire étudié. En vue de l'implémentation des cartes de danger dans les documents de l'aménagement du territoire, nous nous intéressons au fonctionnement actuel de la gestion institutionnelle du risque à Beni Mellal, en étudiant le degré d'intégration de la gestion et la manière dont les connaissances sur les risques influencent le processus de gestion. L'analyse montre que la gestion est marquée par une logique de gestion hiérarchique et la priorité des mesures de protection par rapport aux mesures passives d'aménagement du territoire. Les connaissances sur le risque restent sectorielles, souvent déconnectées. L'innovation dans le domaine de la gestion du risque résulte de collaborations horizontales entre les acteurs ou avec des sources de connaissances externes (par exemple les universités). Des recommandations méthodologiques et institutionnelles issues de cette étude ont été adressées aux gestionnaires en vue de l'implémentation des cartes de danger. Plus que des outils de réduction du risque, les cartes de danger aident à transmettre des connaissances vers le public et contribuent ainsi à établir une culture du risque. - Severe rainfall events are thought to be occurring more frequently in semi-arid areas. In Morocco, flood hazard has become an important topic, notably as rapid economic development and high urbanization rates have increased the exposure of people and assets in hazard-prone areas. The Swiss Agency for Development and Cooperation (SADC) is active in natural hazard mitigation in Morocco. As hazard mapping for urban planning is thought to be a sound tool for vulnerability reduction, the SADC has financed a project aimed at adapting the Swiss approach for hazard assessment and mapping to the case of Morocco. In a knowledge transfer context, the Swiss method was adapted to the semi-arid environment, the specific piedmont morphology and to socio-economic constraints particular to the study site. Following the Swiss guidelines, a hydro-geomorphological map was established, containing all geomorphic elements related to known past floods. Next, rainfall / runoff modeling for reference events and hydraulic routing of the obtained hydrographs were carried out in order to assess hazard quantitatively. Field-collected discharge estimations and flood extent for known floods were used to verify the model results. Flood hazard intensity and probability maps were obtained. Finally, an indicative danger map as defined within the Swiss hazard assessment terminology was calculated using the Swiss hazard matrix that convolves flood intensity with its recurrence probability in order to assign flood danger degrees to the concerned territory. Danger maps become effective, as risk mitigation tools, when implemented in urban planning. We focus on how local authorities are involved in the risk management process and how knowledge about risk impacts the management. An institutional vulnerability "map" was established based on individual interviews held with the main institutional actors in flood management. Results show that flood hazard management is defined by uneven actions and relationships, it is based on top-down decision-making patterns, and focus is maintained on active mitigation measures. The institutional actors embody sectorial, often disconnected risk knowledge pools, whose relationships are dictated by the institutional hierarchy. Results show that innovation in the risk management process emerges when actors collaborate despite the established hierarchy or when they open to outer knowledge pools (e.g. the academia). Several methodological and institutional recommendations were addressed to risk management stakeholders in view of potential map implementation to planning. Hazard assessment and mapping is essential to an integrated risk management approach: more than a mitigation tool, danger maps represent tools that allow communicating on hazards and establishing a risk culture.

Regional-scale integration of multi-scale hydrological and geophysical data using a two-step Bayesian sequential simulation approach

Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale for the purpose of improving predictions of groundwater flow and solute transport. However, extending corresponding approaches to the regional scale still represents one of the major challenges in the domain of hydrogeophysics. To address this problem, we have developed a regional-scale data integration methodology based on a two-step Bayesian sequential simulation approach. Our objective is to generate high-resolution stochastic realizations of the regional-scale hydraulic conductivity field in the common case where there exist spatially exhaustive but poorly resolved measurements of a related geophysical parameter, as well as highly resolved but spatially sparse collocated measurements of this geophysical parameter and the hydraulic conductivity. To integrate this multi-scale, multi-parameter database, we first link the low- and high-resolution geophysical data via a stochastic downscaling procedure. This is followed by relating the downscaled geophysical data to the high-resolution hydraulic conductivity distribution. After outlining the general methodology of the approach, we demonstrate its application to a realistic synthetic example where we consider as data high-resolution measurements of the hydraulic and electrical conductivities at a small number of borehole locations, as well as spatially exhaustive, low-resolution estimates of the electrical conductivity obtained from surface-based electrical resistivity tomography. The different stochastic realizations of the hydraulic conductivity field obtained using our procedure are validated by comparing their solute transport behaviour with that of the underlying ?true? hydraulic conductivity field. We find that, even in the presence of strong subsurface heterogeneity, our proposed procedure allows for the generation of faithful representations of the regional-scale hydraulic conductivity structure and reliable predictions of solute transport over long, regional-scale distances.

Imaging and quantifying salt-tracer transport in a riparian groundwater system by means of 3D ERT monitoring

Determining groundwater flow paths of infiltrated river water is necessary for studying biochemical processes in the riparian zone, but their characterization is complicated by strong temporal and spatial heterogeneity. We investigated to what extent repeat 3D surface electrical resistance tomography (ERT) can be used to monitor transport of a salt-tracer plume under close to natural gradient conditions. The aim is to estimate groundwater flow velocities and pathways at a site located within a riparian groundwater system adjacent to the perialpine Thur River in northeastern Switzerland. Our ERT time-lapse images provide constraints on the plume's shape, flow direction, and velocity. These images allow the movement of the plume to be followed for 35 m. Although the hydraulic gradient is only 1.43 parts per thousand, the ERT time-lapse images demonstrate that the plume's center of mass and its front propagate with velocities of 2x10(-4) m/s and 5x10(-4) m/s, respectively. These velocities are compatible with groundwater resistivity monitoring data in two observation wells 5 m from the injection well. Five additional sensors in the 5-30 m distance range did not detect the plume. Comparison of the ERT time-lapse images with a groundwater transport model and time-lapse inversions of synthetic ERT data indicate that the movement of the plume can be described for the first 6 h after injection by a uniform transport model. Subsurface heterogeneity causes a change of the plume's direction and velocity at later times. Our results demonstrate the effectiveness of using time-lapse 3D surface ERT to monitor flow pathways in a challenging perialpine environment over larger scales than is practically possible with crosshole 3D ERT.

Construcció i explotació d'un magatzem de dades de planificació hidrològica

El projecte que constitueix la base de aquest TFC, consisteix en la construcció d'un magatzem de dades, com a eina de treball, per als tècnics de la Confederació Hidrogràfica Nord-Est (CHNE), que és l'entitat encarregada de la gestió de l'aigua i infraestructures hidràuliques del nord-est de la península ibèrica. Aquest magatzem ha de permetre automatitzar la recollida de la informació provinent de fulls Excel i explotar aquesta informació amb eines de Business Intelligence adients.

A direct approach to the discounted penalty function

This paper provides a new and accessible approach to establishing certain results concerning the discounted penalty function. The direct approach consists of two steps. In the first step, closed-form expressions are obtained in the special case in which the claim amount distribution is a combination of exponential distributions. A rational function is useful in this context. For the second step, one observes that the family of combinations of exponential distributions is dense. Hence, it suffices to reformulate the results of the first step to obtain general results. The surplus process has downward and upward jumps, modeled by two independent compound Poisson processes. If the distribution of the upward jumps is exponential, a series of new results can be obtained with ease. Subsequently, certain results of Gerber and Shiu [H. U. Gerber and E. S. W. Shiu, North American Actuarial Journal 2(1): 48–78 (1998)] can be reproduced. The two-step approach is also applied when an independent Wiener process is added to the surplus process. Certain results are related to Zhang et al. [Z. Zhang, H. Yang, and S. Li, Journal of Computational and Applied Mathematics 233: 1773–1 784 (2010)], which uses different methods.

Conditioning of stochastic 3-D fracture networks to hydrological and geophysical data

The geometry and connectivity of fractures exert a strong influence on the flow and transport properties of fracture networks. We present a novel approach to stochastically generate three-dimensional discrete networks of connected fractures that are conditioned to hydrological and geophysical data. A hierarchical rejection sampling algorithm is used to draw realizations from the posterior probability density function at different conditioning levels. The method is applied to a well-studied granitic formation using data acquired within two boreholes located 6 m apart. The prior models include 27 fractures with their geometry (position and orientation) bounded by information derived from single-hole ground-penetrating radar (GPR) data acquired during saline tracer tests and optical televiewer logs. Eleven cross-hole hydraulic connections between fractures in neighboring boreholes and the order in which the tracer arrives at different fractures are used for conditioning. Furthermore, the networks are conditioned to the observed relative hydraulic importance of the different hydraulic connections by numerically simulating the flow response. Among the conditioning data considered, constraints on the relative flow contributions were the most effective in determining the variability among the network realizations. Nevertheless, we find that the posterior model space is strongly determined by the imposed prior bounds. Strong prior bounds were derived from GPR measurements and helped to make the approach computationally feasible. We analyze a set of 230 posterior realizations that reproduce all data given their uncertainties assuming the same uniform transmissivity in all fractures. The posterior models provide valuable statistics on length scales and density of connected fractures, as well as their connectivity. In an additional analysis, effective transmissivity estimates of the posterior realizations indicate a strong influence of the DFN structure, in that it induces large variations of equivalent transmissivities between realizations. The transmissivity estimates agree well with previous estimates at the site based on pumping, flowmeter and temperature data.

Use of high-resolution geophysical data to characterize heterogeneous aquifers: influence of data integration method on hydrological predictions

The integration of geophysical data into the subsurface characterization problem has been shown in many cases to significantly improve hydrological knowledge by providing information at spatial scales and locations that is unattainable using conventional hydrological measurement techniques. The investigation of exactly how much benefit can be brought by geophysical data in terms of its effect on hydrological predictions, however, has received considerably less attention in the literature. Here, we examine the potential hydrological benefits brought by a recently introduced simulated annealing (SA) conditional stochastic simulation method designed for the assimilation of diverse hydrogeophysical data sets. We consider the specific case of integrating crosshole ground-penetrating radar (GPR) and borehole porosity log data to characterize the porosity distribution in saturated heterogeneous aquifers. In many cases, porosity is linked to hydraulic conductivity and thus to flow and transport behavior. To perform our evaluation, we first generate a number of synthetic porosity fields exhibiting varying degrees of spatial continuity and structural complexity. Next, we simulate the collection of crosshole GPR data between several boreholes in these fields, and the collection of porosity log data at the borehole locations. The inverted GPR data, together with the porosity logs, are then used to reconstruct the porosity field using the SA-based method, along with a number of other more elementary approaches. Assuming that the grid-cell-scale relationship between porosity and hydraulic conductivity is unique and known, the porosity realizations are then used in groundwater flow and contaminant transport simulations to assess the benefits and limitations of the different approaches.

Candida parapsilosis complex water isolates from a haemodialysis unit: biofilm production and in vitro evaluation of the use of clinical antifungals

Candida parapsilosis, currently divided into three distinct species, proliferates in glucose-rich solutions and has been associated with infections resulting from the use of medical devices made of plastic, an environment common in dialysis centres. The aims of this study were (i) to screen for Candida orthopsilosis and Candida metapsilosis (100 environmental isolates previously identified as C. parapsilosis), (ii) to test the ability of these isolates to form biofilm and (iii) to investigate the in vitro susceptibility of Candida spp biofilms to the antifungal agents, fluconazole (FLC) and amphotericin B (AMB). Isolates were obtained from a hydraulic circuit collected from a haemodialysis unit. Based on molecular criteria, 47 strains were re-identified as C. orthopsilosis and 53 as C. parapsilosis. Analyses using a formazan salt reduction assay and total viable count, together with microscopy studies, revealed that 72 strains were able to form biofilm that was structurally similar, but with minor differences in morphology. A microtitre-based colorimetric assay used to test the susceptibility of fungal biofilms to AMB and FLC demonstrated that the C. parapsilosis complex displayed an increased resistance to these antifungal agents. The results from these analyses may provide a basis for implementing quality controls and monitoring to ensure the microbiological purity of dialysis water, including the presence of yeast.