MIGCLIM: Predicting plant distribution and dispersal in a changing climate

Many studies have forecasted the possible impact of climate change on plant distribution using models based on ecological niche theory. In their basic implementation, niche-based models do not constrain predictions by dispersal limitations. Hence, most niche-based modelling studies published so far have assumed dispersal to be either unlimited or null. However, depending on the rate of climatic change, the landscape fragmentation and the dispersal capabilities of individual species, these assumptions are likely to prove inaccurate, leading to under- or overestimation of future species distributions and yielding large uncertainty between these two extremes. As a result, the concepts of "potentially suitable" and "potentially colonisable" habitat are expected to differ significantly. To quantify to what extent these two concepts can differ, we developed MIGCLIM, a model simulating plant dispersal under climate change and landscape fragmentation scenarios. MIGCLIM implements various parameters, such as dispersal distance, increase in reproductive potential over time, barriers to dispersal or long distance dispersal. Several simulations were run for two virtual species in a study area of the western Swiss Alps, by varying dispersal distance and other parameters. Each simulation covered the hundred-year period 2001-2100 and three different IPCC-based temperature warming scenarios were considered. Our results indicate that: (i) using realistic parameter values, the future potential distributions generated using MIGCLIM can differ significantly (up to more than 95% decrease in colonized surface) from those that ignore dispersal; (ii) this divergence increases both with increasing climate warming and over longer time periods; (iii) the uncertainty associated with the warming scenario can be nearly as large as the one related to dispersal parameters; (iv) accounting for dispersal, even roughly, can importantly reduce uncertainty in projections.

Landscape structure affects dispersal in the greater white-toothed shrew: inference between genetic and simulated ecological distances

Animal dispersal in a fragmented landscape depends on the complex interaction between landscape structure and animal behavior. To better understand how individuals disperse, it is important to explicitly represent the properties of organisms and the landscape in which they move. A common approach to modelling dispersal includes representing the landscape as a grid of equal sized cells and then simulating individual movement as a correlated random walk. This approach uses a priori scale of resolution, which limits the representation of all landscape features and how different dispersal abilities are modelled. We develop a vector-based landscape model coupled with an object-oriented model for animal dispersal. In this spatially explicit dispersal model, landscape features are defined based on their geographic and thematic properties and dispersal is modelled through consideration of an organism's behavior, movement rules and searching strategies (such as visual cues). We present the model's underlying concepts, its ability to adequately represent landscape features and provide simulation of dispersal according to different dispersal abilities. We demonstrate the potential of the model by simulating two virtual species in a real Swiss landscape. This illustrates the model's ability to simulate complex dispersal processes and provides information about dispersal such as colonization probability and spatial distribution of the organism's path

Analysis, modelling and visualisation of spatiotemporal data for urban studies

The proportion of population living in or around cites is more important than ever. Urban sprawl and car dependence have taken over the pedestrian-friendly compact city. Environmental problems like air pollution, land waste or noise, and health problems are the result of this still continuing process. The urban planners have to find solutions to these complex problems, and at the same time insure the economic performance of the city and its surroundings. At the same time, an increasing quantity of socio-economic and environmental data is acquired. In order to get a better understanding of the processes and phenomena taking place in the complex urban environment, these data should be analysed. Numerous methods for modelling and simulating such a system exist and are still under development and can be exploited by the urban geographers for improving our understanding of the urban metabolism. Modern and innovative visualisation techniques help in communicating the results of such models and simulations. This thesis covers several methods for analysis, modelling, simulation and visualisation of problems related to urban geography. The analysis of high dimensional socio-economic data using artificial neural network techniques, especially self-organising maps, is showed using two examples at different scales. The problem of spatiotemporal modelling and data representation is treated and some possible solutions are shown. The simulation of urban dynamics and more specifically the traffic due to commuting to work is illustrated using multi-agent micro-simulation techniques. A section on visualisation methods presents cartograms for transforming the geographic space into a feature space, and the distance circle map, a centre-based map representation particularly useful for urban agglomerations. Some issues on the importance of scale in urban analysis and clustering of urban phenomena are exposed. A new approach on how to define urban areas at different scales is developed, and the link with percolation theory established. Fractal statistics, especially the lacunarity measure, and scale laws are used for characterising urban clusters. In a last section, the population evolution is modelled using a model close to the well-established gravity model. The work covers quite a wide range of methods useful in urban geography. Methods should still be developed further and at the same time find their way into the daily work and decision process of urban planners. La part de personnes vivant dans une région urbaine est plus élevé que jamais et continue à croître. L'étalement urbain et la dépendance automobile ont supplanté la ville compacte adaptée aux piétons. La pollution de l'air, le gaspillage du sol, le bruit, et des problèmes de santé pour les habitants en sont la conséquence. Les urbanistes doivent trouver, ensemble avec toute la société, des solutions à ces problèmes complexes. En même temps, il faut assurer la performance économique de la ville et de sa région. Actuellement, une quantité grandissante de données socio-économiques et environnementales est récoltée. Pour mieux comprendre les processus et phénomènes du système complexe "ville", ces données doivent être traitées et analysées. Des nombreuses méthodes pour modéliser et simuler un tel système existent et sont continuellement en développement. Elles peuvent être exploitées par le géographe urbain pour améliorer sa connaissance du métabolisme urbain. Des techniques modernes et innovatrices de visualisation aident dans la communication des résultats de tels modèles et simulations. Cette thèse décrit plusieurs méthodes permettant d'analyser, de modéliser, de simuler et de visualiser des phénomènes urbains. L'analyse de données socio-économiques à très haute dimension à l'aide de réseaux de neurones artificiels, notamment des cartes auto-organisatrices, est montré à travers deux exemples aux échelles différentes. Le problème de modélisation spatio-temporelle et de représentation des données est discuté et quelques ébauches de solutions esquissées. La simulation de la dynamique urbaine, et plus spécifiquement du trafic automobile engendré par les pendulaires est illustrée à l'aide d'une simulation multi-agents. Une section sur les méthodes de visualisation montre des cartes en anamorphoses permettant de transformer l'espace géographique en espace fonctionnel. Un autre type de carte, les cartes circulaires, est présenté. Ce type de carte est particulièrement utile pour les agglomérations urbaines. Quelques questions liées à l'importance de l'échelle dans l'analyse urbaine sont également discutées. Une nouvelle approche pour définir des clusters urbains à des échelles différentes est développée, et le lien avec la théorie de la percolation est établi. Des statistiques fractales, notamment la lacunarité, sont utilisées pour caractériser ces clusters urbains. L'évolution de la population est modélisée à l'aide d'un modèle proche du modèle gravitaire bien connu. Le travail couvre une large panoplie de méthodes utiles en géographie urbaine. Toutefois, il est toujours nécessaire de développer plus loin ces méthodes et en même temps, elles doivent trouver leur chemin dans la vie quotidienne des urbanistes et planificateurs.

Efficacy of levofloxacin in the treatment of experimental endocarditis caused by viridans group streptococci.

Levofloxacin was investigated against viridans group streptococci in vitro and in rats with experimental aortic endocarditis. The MIC(90)s of levofloxacin and ciprofloxacin for 20 independent isolates of such bacteria were 1 and 8 mg/L, respectively. Rats were infected with two types of organism: either fully susceptible to levofloxacin MIC < or = 0.5 mg/L) or borderline susceptible (MIC 1-2 mg/L). Fully levofloxacin-susceptible bacteria comprised one penicillin-susceptible (MIC 0.004 mg/L) Streptococcus gordonii, and one penicillin-tolerant as well as one intermediate penicillin-resistant (MIC 0.125 mg/L) isogenic strains. Borderline levofloxacin-susceptible bacteria comprised one penicillin-susceptible Streptococcus sanguis and one highly penicillin-resistant Streptococcus mitis (MIC 2 mg/L). Rats were treated for 5 days with drug dosages simulating the following treatments in humans: (i) levofloxacin 500 mg orally once a day (q24 h), (ii) levofloxacin 500 mg orally twice a day (q12 h), (iii) levofloxacin 1 g orally q24 h, (iv) ciprofloxacin 750 mg orally q12 h, and (v) ceftriaxone 2 g iv q24 h. Levofloxacin was equivalent or superior to ceftriaxone, and was successful in treating experimental endocarditis irrespective of penicillin resistance. Nevertheless, standard levofloxacin treatment equivalent to 500 mg q24 h in human was less effective than twice daily 500 mg or once daily 1 g doses against borderline-susceptible organisms. Ciprofloxacin, used as a negative control, was ineffective and selected for resistant isolates. This underlines the importance of MIC determinations when treating severe streptococcal infection with quinolones. In the case of borderline-susceptible pathogens, total daily doses of 1 g of levofloxacin should be considered.

Sex differences in urinary levels of several biological indicators of exposure: a simulation study using a compartmental based toxicokinetic model

Toxicokinetic modeling is a useful tool to describe or predict the behavior of a chemical agent in the human or animal organism. A general model based on four compartments was developed in a previous study in order to quantify the effect of human variability on a wide range of biological exposure indicators. The aim of this study was to adapt this existing general toxicokinetic model to three organic solvents, which were methyl ethyl ketone, 1-methoxy-2-propanol and 1,1,1,-trichloroethane, and to take into account sex differences. We assessed in a previous human volunteer study the impact of sex on different biomarkers of exposure corresponding to the three organic solvents mentioned above. Results from that study suggested that not only physiological differences between men and women but also differences due to sex hormones levels could influence the toxicokinetics of the solvents. In fact the use of hormonal contraceptive had an effect on the urinary levels of several biomarkers, suggesting that exogenous sex hormones could influence CYP2E1 enzyme activity. These experimental data were used to calibrate the toxicokinetic models developed in this study. Our results showed that it was possible to use an existing general toxicokinetic model for other compounds. In fact, most of the simulation results showed good agreement with the experimental data obtained for the studied solvents, with a percentage of model predictions that lies within the 95% confidence interval varying from 44.4 to 90%. Results pointed out that for same exposure conditions, men and women can show important differences in urinary levels of biological indicators of exposure. Moreover, when running the models by simulating industrial working conditions, these differences could even be more pronounced. In conclusion, a general and simple toxicokinetic model, adapted for three well known organic solvents, allowed us to show that metabolic parameters can have an important impact on the urinary levels of the corresponding biomarkers. These observations give evidence of an interindividual variablity, an aspect that should have its place in the approaches for setting limits of occupational exposure.

Multiscale descriptions of density-driven flow instabilities in porous media

Les instabilités engendrées par des gradients de densité interviennent dans une variété d'écoulements. Un exemple est celui de la séquestration géologique du dioxyde de carbone en milieux poreux. Ce gaz est injecté à haute pression dans des aquifères salines et profondes. La différence de densité entre la saumure saturée en CO2 dissous et la saumure environnante induit des courants favorables qui le transportent vers les couches géologiques profondes. Les gradients de densité peuvent aussi être la cause du transport indésirable de matières toxiques, ce qui peut éventuellement conduire à la pollution des sols et des eaux. La gamme d'échelles intervenant dans ce type de phénomènes est très large. Elle s'étend de l'échelle poreuse où les phénomènes de croissance des instabilités s'opèrent, jusqu'à l'échelle des aquifères à laquelle interviennent les phénomènes à temps long. Une reproduction fiable de la physique par la simulation numérique demeure donc un défi en raison du caractère multi-échelles aussi bien au niveau spatial et temporel de ces phénomènes. Il requiert donc le développement d'algorithmes performants et l'utilisation d'outils de calculs modernes. En conjugaison avec les méthodes de résolution itératives, les méthodes multi-échelles permettent de résoudre les grands systèmes d'équations algébriques de manière efficace. Ces méthodes ont été introduites comme méthodes d'upscaling et de downscaling pour la simulation d'écoulements en milieux poreux afin de traiter de fortes hétérogénéités du champ de perméabilité. Le principe repose sur l'utilisation parallèle de deux maillages, le premier est choisi en fonction de la résolution du champ de perméabilité (grille fine), alors que le second (grille grossière) est utilisé pour approximer le problème fin à moindre coût. La qualité de la solution multi-échelles peut être améliorée de manière itérative pour empêcher des erreurs trop importantes si le champ de perméabilité est complexe. Les méthodes adaptatives qui restreignent les procédures de mise à jour aux régions à forts gradients permettent de limiter les coûts de calculs additionnels. Dans le cas d'instabilités induites par des gradients de densité, l'échelle des phénomènes varie au cours du temps. En conséquence, des méthodes multi-échelles adaptatives sont requises pour tenir compte de cette dynamique. L'objectif de cette thèse est de développer des algorithmes multi-échelles adaptatifs et efficaces pour la simulation des instabilités induites par des gradients de densité. Pour cela, nous nous basons sur la méthode des volumes finis multi-échelles (MsFV) qui offre l'avantage de résoudre les phénomènes de transport tout en conservant la masse de manière exacte. Dans la première partie, nous pouvons démontrer que les approximations de la méthode MsFV engendrent des phénomènes de digitation non-physiques dont la suppression requiert des opérations de correction itératives. Les coûts de calculs additionnels de ces opérations peuvent toutefois être compensés par des méthodes adaptatives. Nous proposons aussi l'utilisation de la méthode MsFV comme méthode de downscaling: la grille grossière étant utilisée dans les zones où l'écoulement est relativement homogène alors que la grille plus fine est utilisée pour résoudre les forts gradients. Dans la seconde partie, la méthode multi-échelle est étendue à un nombre arbitraire de niveaux. Nous prouvons que la méthode généralisée est performante pour la résolution de grands systèmes d'équations algébriques. Dans la dernière partie, nous focalisons notre étude sur les échelles qui déterminent l'évolution des instabilités engendrées par des gradients de densité. L'identification de la structure locale ainsi que globale de l'écoulement permet de procéder à un upscaling des instabilités à temps long alors que les structures à petite échelle sont conservées lors du déclenchement de l'instabilité. Les résultats présentés dans ce travail permettent d'étendre les connaissances des méthodes MsFV et offrent des formulations multi-échelles efficaces pour la simulation des instabilités engendrées par des gradients de densité. - Density-driven instabilities in porous media are of interest for a wide range of applications, for instance, for geological sequestration of CO2, during which CO2 is injected at high pressure into deep saline aquifers. Due to the density difference between the C02-saturated brine and the surrounding brine, a downward migration of CO2 into deeper regions, where the risk of leakage is reduced, takes place. Similarly, undesired spontaneous mobilization of potentially hazardous substances that might endanger groundwater quality can be triggered by density differences. Over the last years, these effects have been investigated with the help of numerical groundwater models. Major challenges in simulating density-driven instabilities arise from the different scales of interest involved, i.e., the scale at which instabilities are triggered and the aquifer scale over which long-term processes take place. An accurate numerical reproduction is possible, only if the finest scale is captured. For large aquifers, this leads to problems with a large number of unknowns. Advanced numerical methods are required to efficiently solve these problems with today's available computational resources. Beside efficient iterative solvers, multiscale methods are available to solve large numerical systems. Originally, multiscale methods have been developed as upscaling-downscaling techniques to resolve strong permeability contrasts. In this case, two static grids are used: one is chosen with respect to the resolution of the permeability field (fine grid); the other (coarse grid) is used to approximate the fine-scale problem at low computational costs. The quality of the multiscale solution can be iteratively improved to avoid large errors in case of complex permeability structures. Adaptive formulations, which restrict the iterative update to domains with large gradients, enable limiting the additional computational costs of the iterations. In case of density-driven instabilities, additional spatial scales appear which change with time. Flexible adaptive methods are required to account for these emerging dynamic scales. The objective of this work is to develop an adaptive multiscale formulation for the efficient and accurate simulation of density-driven instabilities. We consider the Multiscale Finite-Volume (MsFV) method, which is well suited for simulations including the solution of transport problems as it guarantees a conservative velocity field. In the first part of this thesis, we investigate the applicability of the standard MsFV method to density- driven flow problems. We demonstrate that approximations in MsFV may trigger unphysical fingers and iterative corrections are necessary. Adaptive formulations (e.g., limiting a refined solution to domains with large concentration gradients where fingers form) can be used to balance the extra costs. We also propose to use the MsFV method as downscaling technique: the coarse discretization is used in areas without significant change in the flow field whereas the problem is refined in the zones of interest. This enables accounting for the dynamic change in scales of density-driven instabilities. In the second part of the thesis the MsFV algorithm, which originally employs one coarse level, is extended to an arbitrary number of coarse levels. We prove that this keeps the MsFV method efficient for problems with a large number of unknowns. In the last part of this thesis, we focus on the scales that control the evolution of density fingers. The identification of local and global flow patterns allows a coarse description at late times while conserving fine-scale details during onset stage. Results presented in this work advance the understanding of the Multiscale Finite-Volume method and offer efficient dynamic multiscale formulations to simulate density-driven instabilities. - Les nappes phréatiques caractérisées par des structures poreuses et des fractures très perméables représentent un intérêt particulier pour les hydrogéologues et ingénieurs environnementaux. Dans ces milieux, une large variété d'écoulements peut être observée. Les plus communs sont le transport de contaminants par les eaux souterraines, le transport réactif ou l'écoulement simultané de plusieurs phases non miscibles, comme le pétrole et l'eau. L'échelle qui caractérise ces écoulements est définie par l'interaction de l'hétérogénéité géologique et des processus physiques. Un fluide au repos dans l'espace interstitiel d'un milieu poreux peut être déstabilisé par des gradients de densité. Ils peuvent être induits par des changements locaux de température ou par dissolution d'un composé chimique. Les instabilités engendrées par des gradients de densité revêtent un intérêt particulier puisque qu'elles peuvent éventuellement compromettre la qualité des eaux. Un exemple frappant est la salinisation de l'eau douce dans les nappes phréatiques par pénétration d'eau salée plus dense dans les régions profondes. Dans le cas des écoulements gouvernés par les gradients de densité, les échelles caractéristiques de l'écoulement s'étendent de l'échelle poreuse où les phénomènes de croissance des instabilités s'opèrent, jusqu'à l'échelle des aquifères sur laquelle interviennent les phénomènes à temps long. Etant donné que les investigations in-situ sont pratiquement impossibles, les modèles numériques sont utilisés pour prédire et évaluer les risques liés aux instabilités engendrées par les gradients de densité. Une description correcte de ces phénomènes repose sur la description de toutes les échelles de l'écoulement dont la gamme peut s'étendre sur huit à dix ordres de grandeur dans le cas de grands aquifères. Il en résulte des problèmes numériques de grande taille qui sont très couteux à résoudre. Des schémas numériques sophistiqués sont donc nécessaires pour effectuer des simulations précises d'instabilités hydro-dynamiques à grande échelle. Dans ce travail, nous présentons différentes méthodes numériques qui permettent de simuler efficacement et avec précision les instabilités dues aux gradients de densité. Ces nouvelles méthodes sont basées sur les volumes finis multi-échelles. L'idée est de projeter le problème original à une échelle plus grande où il est moins coûteux à résoudre puis de relever la solution grossière vers l'échelle de départ. Cette technique est particulièrement adaptée pour résoudre des problèmes où une large gamme d'échelle intervient et évolue de manière spatio-temporelle. Ceci permet de réduire les coûts de calculs en limitant la description détaillée du problème aux régions qui contiennent un front de concentration mobile. Les aboutissements sont illustrés par la simulation de phénomènes tels que l'intrusion d'eau salée ou la séquestration de dioxyde de carbone.

Modelling the influence of change in fire regime on the local distribution of a Mediterranean pyrophytic plant species (Cistus salviifolius) at its northern range limit

Aims: To assess the potential distribution of an obligate seeder and active pyrophyte, Cistus salviifolius, a vulnerable species in the Swiss Red List; to derive scenarios by changing the fire return interval; and to discuss the results from a conservation perspective. A more general aim is to assess the impact of fire as a natural factor influencing the vegetation of the southern slopes of the Alps. Locations: Alps, southern Switzerland. Methods: Presence-absence data to fit the model were obtained from the most recent field mapping of C. salviifolius. The quantitative environmental predictors used in this study include topographic, climatic and disturbance (fire) predictors. Models were fitted by logistic regression and evaluated by jackknife and bootstrap approaches. Changes in fire regime were simulated by increasing the time-return interval of fire (simulating longer periods without fire). Two scenarios were considered: no fire in the past 15 years; or in the past 35 years. Results: Rock cover, slope, topographic position, potential evapotranspiration and time elapsed since the last fire were selected in the final model. The Nagelkerke R-2 of the model for C. salviifolius was 0.57 and the Jackknife area under the curve evaluation was 0.89. The bootstrap evaluation revealed model robustness. By increasing the return interval of fire by either up to 15 years, or 35 years, the modelled C. salviifolius population declined by 30-40%, respectively. Main conclusions: Although fire plays a significant role, topography and rock cover appear to be the most important predictors, suggesting that the distribution of C. salviifolius in the southern Swiss Alps is closely related to the availability of supposedly competition-free sites, such as emerging bedrock, ridge locations or steep slopes. Fire is more likely to play a secondary role in allowing C. salviifolius to extend its occurrence temporarily, by increasing germination rates and reducing the competition from surrounding vegetation. To maintain a viable dormant seed bank for C. salviifolius, conservation managers should consider carrying out vegetation clearing and managing wild fire propagation to reduce competition and ensure sufficient recruitment for this species.

Flow structures at an idealized bifurcation : a numerical experiment

River bifurcations are key nodes within braided river systems controlling the flow and sediment partitioning and therefore the dynamics of the river braiding process. Recent research has shown that certain geometrical configurations induce instabilities that lead to downstream mid-channel bar formation and the formation of bifurcations. However, we currently have a poor understanding of the flow division process within bifurcations and the flow dynamics in the downstream bifurcates, both of which are needed to understand bifurcation stability. This paper presents results of a numerical sensitivity experiment undertaken using computational fluid dynamics (CFD) with the purpose of understanding the flow dynamics of a series of idealized bifurcations. A geometric sensitivity analysis is undertaken for a range of channel slopes (0.005 to 0.03), bifurcation angles (22 degrees to 42 degrees) and a restricted set of inflow conditions based upon simulating flow through meander bends with different curvature on the flow field dynamics through the bifurcation. The results demonstrate that the overall slope of the bifurcation affects the velocity of flow through the bifurcation and when slope asymmetry is introduced, the flow structures in the bifurcation are modified. In terms of bifurcation evolution the most important observation appears to be that once slope asymmetry is greater than 0.2 the flow within the steep bifurcate shows potential instability and the potential for alternate channel bar formation. Bifurcation angle also defines the flow structures within the bifurcation with an increase in bifurcation angle increasing the flow velocity down both bifurcates. However, redistributive effects of secondary circulation caused by upstream curvature can very easily counter the effects of local bifurcation characteristics. Copyright (C) 2011 John Wiley & Sons, Ltd.

Soil Water Potentials and Capsicum annuum L. under Salinity

ABSTRACT Investigations into water potentials in the soil-plant system are of great relevance in environments with abiotic stresses, such as salinity and drought. An experiment was developed using bell pepper in a Neossolo Flúvico (Fluvent) irrigated with water of six levels of electrical conductivity (0, 1, 3, 5, 7 and 9 dS m-1) by using exclusively NaCl and by simulating the actual condition (using a mixture of salts). The treatments were arranged in a randomized block design, in a 6 × 2 factorial arrangement, with four replicates. Soil matric (Ψm) and osmotic (Ψo) potentials were determined 70 days after transplanting (DAT). Soil total potential was considered as the sum of Ψm and Ψo. Leaf water (obtained with the Scholander Chamber) and osmotic potentials were determined before sunrise (predawn) and at noon at 42 and 70 DAT. There were no significant differences between the salt sources used in the irrigation water for soil and plant water potentials. The supply of salts to the soil through irrigation water was the main factor responsible for the decrease in Ψo in the soil and in bell pepper leaves. The total potential of bell pepper at predawn reached values of -1.30 and -1.33 MPa at 42 and 70 DAT, respectively, when water of 9 dS m-1 was used in the irrigation. The total potential at noon reached -2.19 MPa. The soil subjected to the most saline treatment reached a water potential of -1.20 MPa at 70 DAT. There was no predawn equilibrium between the total water potentials of the soil and the plant, indicating that soil potential cannot be considered similar to that of the plant. The determination of the osmotic potential in the soil solution should not be neglected in saline soils, since it has strong influence on the calculation of the total potential.

Hydrogen isotope fractionations between amphiboles, micas, and fluids in alkaline igneous intrusions

RÉSUMÉ DE LA THÈSE Les teneurs des amphiboles en éléments majeurs et en isotopes stables ont été analysées dans plusieurs complexes ignés alcalins et hyperalcalins, dans le but de déterminer l'importance des variations de composition des minéraux pour le fractionnement isotopique de l'hydrogène dans un système naturel minéral-magma-fluide. Cette étude se concentre principalement sur les syénites néphéliniques de complexes intrusifs alcalins bien connus mais à chimie variable, dont les amphiboles, ainsi que d'autres silicates hydratés tels que micas et eudialytes, lorsque cela était possible, ont été séparés. L'intérêt principal s'est porté sur le complexe alcalin d'Ilímaussaq de la Province du Gardar, au Sud du Groenland. Dans une optique de comparaison, nous avons collecté et analysé d'autres échantillons provenant du complexe de Tugtutôq (Sud Groenland), des complexes de Khibina et Lovozero (Péninsule de Kola, Russie), du Mont St-Hilaire et du Mont Royal (Canada) et de 6 autres du nord-ouest de la Namibie (Cape Cross, Okenyenya, Messum, Etaneno, Kalkfeld,et Okorusu). Les compositions isotopiques de l'hydrogène des amphiboles des ces différentes zones présentent de grandes variations (-227 à -700/00), ce qui est atypique pour des magmas d'origine mantellique. Les valeurs comprises entre -80 et -400/00 indiquent une provenance du manteau. Ces larges variations de compositions ainsi que l'extrême appauvrissement en isotope lourd de l'hydrogène (D), en comparaison avec d'autres roches ignées, semblent être propres.aux roches alcalines et hyperalcalines de ce type, ce qui indiquerait un processus commun. Les différents complexes alcalins choisis présentent un large intervalle de composition chimique des amphiboles. La caractérisation des amphiboles par microscopie électronique et par spectroscopie Mössbauer contribuent à observer le contrôle du Fe sur le fractionnement des isotopes de l'hydrogène. En effet, cela a mis en évidence un contrôle du Fe sur le fractionnement et même, dans le cas du complexe hyperalcalin d'Ilímaussaq, une relation entre le rapport Fei3+/FeT et les variations du rapport D/H. Les complexes étudiés diffèrent de par leur index agpaïtique (Na+K/Al) et également de par leur contenu en fer. Les plus hautes valeurs en Fe (27-35 wt%) et en éléments alcalins dans les amphiboles, ainsi que les teneurs de D/H les plus basses et leur grande variation, sont celles du complexe d'Ilímaussaq. Les amphiboles de la Péninsule de Kola et du Canada sont similaires, mais toutefois moins appauvries en D. En ce qui concerne les amphiboles des complexes du NO de la Namibie, elles présentent des compositions isotopiques de l'hydrogène magmatiques normales (-73 à -100 0/00), contiennent moins de Fe (15-17 wt%) et sont fortement enrichies en Ca et moins en Na. Dans ce cas, l'alcalinité est moins importante en comparaison des autres complexes étudiés. En dehors des teneurs en éléments alcalins des amphiboles, l'alcalinité des fluides s'avère également un facteur important, ce qui est cohérent avec certaines suggestions à partir de systèmes expérimentaux. Afin de mieux contraindre ce facteur, des expériences d'échanges hydrothermaux entre les amphiboles et les fluides de salinité différente ont été effectuées en simulant des conditions naturelles. L'approximation d'amphiboles naturelles de complexes ignés alcalins, couplée aux expériences d'échange, aide à préciser les facteurs contrôlant le fractionnement des isotopes de l'hydrogène dans les roches alcalines. Les valeurs extrêmement basses de 3D des amphiboles de ces complexes alcalins peuvent être dues à une combinaison de différents facteurs, telles qu'une haute alcalinité, une haute teneur en Fe et une faible profondeur d'intrusion. Les grandes variations ainsi que les faibles valeurs de SD des amphiboles étudiées peuvent résulter d'un processus magmatique interne et il est peu probable que de l'eau météorique soit impliquée et/ou que le dégazage magmatique ait joué un rôle. THESIS ABSTRACT Major element and stable isotope compositions of amphiboles were analyzed from a number of alkaline and peralkaline igneous complexes in order to determine the importance of compositional variations in minerals to hydrogen isotope fractionations in natural mineral-melt-fluid systems. The thesis mainly focuses on nepheline syenites of well-studied, but chemically variable alkaline intrusive rocks, from which amphiboles and, if possible, other hydrous silicates such as micas and eudialytes were separated. The system of primary interest was the alkaline Ilímaussaq Complex of the Gardar Province of South Greenland. For the purpose of comparison additional samples were collected and examined from the Tugtutôq Complex (South Greenland), the Khibina and Lovozero Complexes (Kola Peninsular, Russia), Mount St-Hilaire and Mount Royal (Canada) and six further complexes from NW Namibia (Cape Cross, Okenyenya, Messum, Etaneno, Kalkfeld, and Okorusu). The hydrogen isotope compositions of amphiboles from the localities studied differ greatly, which is atypical for amphiboles from mantle, range between - 227 and - 700/00 (latter compatible with a simple mantle origin). As this wide range in compositions and the extreme depletion in the heavy hydrogen isotope (D) content relative to other igneous rocks appear to be unique to alkaline to peralkaline rocks of this type, a common process is indicated. The different alkaline complexes chosen cover a wide range of amphibole chemical compositions. Detailed chemical characterization of amphiboles by electron microprobe and Mössbauer spectroscopy analyses helped to constrain the control of Fe on the H-isotope fractionations. Complete characterization of the chemical compositions of the amphiboles support Fe-control on fractionations and at least for the peralkaline Ilímaussaq complex a relationship between Fe3+/FeT ratios and variations in D/H. The studied complexes differ in their agpaitic index (Na+K/Al) and also in their Fe-content. The most iron (27-35 wt. %) and alkaline element rich amphiboles, with the lowermost D/H ratio, as well with very wide range, are the ones from Ilímaussaq complex. Similar, but less D depleted amphiboles are from the Kola Peninsula and the Canadian localities. The complexes described from NW Namibia have amphiboles with normal magmatic hydrogen isotope composition (-730/00 to -1000/00), and have less Fe-content (15-17 wt. %), and are more Ca-and less Na-rich. In this case alkalinity is not that important in comparison to the other studied complexes. Beside the alkaline element contents in the amphiboles, the alkalinity of the fluids has been found to be an important factor, in conjunction with earlier suggestions from experimental systems. To further constrain this factor, hydrothermal exchange experiments between amphiboles and fluids of different salinity simulating natural conditions were performed. The approach of examining natural amphiboles from alkaline igneous complexes in parallel to performing exchange experiments - helped to further constrain the factors controlling the H-isotope fractionations in alkaline rocks. The observed changes between the hydrogen and oxygen isotope compositions of amphiboles and fluids before and after the experiments suggest that another phase was produced during the experiments, which influenced the final hydrogen isotope composition of the system. This presumably hydrous phase has also influenced the Fe3 +/Fe2+ ratio of the amphiboles, which became more oxidized. The extremely low SD values of amphiboles in these alkaline complexes may be due to a combination of different factors such as high alkalinity, high Fe-content, and shallow intrusion depths. This wide range and the low SD values of the amphiboles studied might be a result of internal, magmatic processes and it is unlikely that meteoric water was involved and/or magmatic degassing played an important role. RÉSUMÉ DE LA THÈSE (pour le grand public) Fractionnement isotopique de l'hydrogène entre amphiboles, micas et fluides dans des intrusions alcalines Zsófia Wáczek Directeur de thèse, Prof. Torsten W. Vennemann Institut de Minéralogie et Géochimie, Université de Lausanne Les roches alcalines et celles qui leurs sont associées sont des sources importantes de nombreux minéraux et minerais, tels l'apatite, le niobium, le diamant et autres pierres précieuses. Cette étude se concentre sur des complexes alcalins localisés dans le sud du Groenland, au Canada, dans la péninsule de Kola en Russie et au nord-ouest de la Namibie. Ces complexes sont composés de roches ayant cristallisé à partir de magmas et de fluides très enrichis en alcalins. Cet enrichissement permet la précipitation de minéraux inhabituels riches en potassium et/ou sodium, telles les amphiboles sodiques, également enrichies en fer. Les amphiboles étudiées ont des compositions calciques, sodi-calciques et sodiques, qui reflètent leurs différents environnements de formation. Des études précédentes ont révélé une large gamme de rapports isotopiques de l'hydrogène dans les amphiboles de roches hyperalcalines, dont certains extrêmement bas. Cette variation importante est très intrigante, sachant que des valeurs entre -40 et -800/00 correspondent à des silicates ignés hydratés et non altérés, alors que des valeurs descendant jusqu'a -1500/00 nécessiteraient une altération par de l'eau météorique et/ou une contamination par les roches environnantes ou des sédiments riches en matière organique. Dans lé cas précis du complexe d'Ilímaussaq (sud du Groenland), aucune de ces explications n'a pu être démontrée et des valeurs encore plus faibles ont été trouvées. Le complexe d'Ilímaussaq présente des valeurs de rapport isotopique de l'hydrogène entre -227 et -500/00 dans les amphiboles. Une origine mantellique permet d'expliquer les valeurs élevées, mais d'autres processus doivent entrer en jeu pour engendrer les valeurs les plus négatives. C'est à l'identification de ces processus que nous nous sommes attachés dans ce travail. Les grandes variations observées dans les teneurs en fer et dans le rapport Fe3+/FeT des roches et des minéraux de ces complexes sont corrélées avec d'autres paramètres chimiques, tels que la composition isotopique de l'hydrogène dans les amphiboles. Nous avons dès lors abordé les questions suivantes: quelle est la relation entre la teneur en fer des amphiboles et leur composition isotopique? Que nous apprennent les changements de la teneur en fer et les changements dans le rapport Fe3+/FeT sur les processus pétrologiques dans ces roches? Pour répondre à ces questions, nous avons analysé les compositions isotopiques de l'oxygène et de l'hydrogène dans les amphiboles et d'autres silicates hydratés. La composition chimique et le rapport Fe3+/FeT des amphiboles ont également été déterminés. Des expériences hydrothermales simulant des conditions naturelles ont été entreprises afin de mieux comprendre les processus de fractionnement isotopiques dans ces systèmes très alcalins. Nos conclusions sont les suivantes: (1) Les valeurs extrêmement faibles ainsi que les larges variations des rapports isotopiques de l'hydrogène des amphiboles de ces complexes alcalins sont dues à une combinaison de facteurs tels que la forte alcalinité, la haute teneur en fer et la profondeur très faible de l'intrusion. (2) Ces valeurs sont probablement le résultat de processus magmatiques internes. (3) Il est peu probable que les eaux météoriques et/ou le dégazage magmatique aient joué un rôle lors de la formation de ces amphiboles. (4) Certaines corrélations, en accord avec les études précédentes, ont pu être trouvées au niveau des concentrations en fer. (5) Dans le cas du complexe d'Ilímaussaq exclusivement, une relation a été trouvée entre le rapport Fe3+/FeT et la composition isotopique de l'hydrogène des amphiboles.

Post-Plutonic Magmatism: from the Source to the Emplacement of an Upper Crustal Dyke Swarm, Adamello Massif, Italy

Magmas of the arc-tholeiitic and calc-alkaline differentiation suites contribute substantially to the formation of continental crust in subduction zones. Different geochemical-petrological models have been put forward to achieve evolved magmas forming large volumes of tonalitic to granitic plutons, building an important part of the continental crust. Primary magmas produced in the mantle wedge overlying the subducted slab migrate through the mantle and the crust. During the transfer, magma can accumulate in intermediate reservoirs at different levels where crystallization leads to differentiation and the heat transfer from the magma, together with gained heat from solidification, lead to partial melting of the crust. Partial melts can be assimilated and mix with more primitive magma. Moreover, already formed crystal cumulates or crystal mushes can be recycled and reactivated to transfer to higher crustal levels. Magma transport in the crust involves fow through fractures within a brittle elastic rock. The solidified magma filled crack, a dyke, can crosscut previously formed geological structures and thus serves as a relative or absolute time marker. The study area is situated in the Adamello massif. The Adamello massif is a composite of plutons that were emplaced between 42 and 29 million years. A later dyke swarm intruded into the southern part of the Adamello Batholith. A fractionation model covering dyke compositions from picrobasalts to dacites results in the cummulative crystallization of 17% olivine, 2% Cr-rich spinel, 18% clinopyroxene, 41% amphibole, 4% plagioclase and 0.1% magnetite to achieve an andesitic composition out of a hydrous primitive picrobasalt. These rocks show a similar geochemical evolution as experimental data simulating fractional crystallization and associated magma differentiation at lower crustal depth (7-10 kbar). The peraluminous, corundum normative composition is one characteristic of more evolved dacitic magmas, which has been explained in a long lasting debate with two di_erent models. Melting of mafic crust or politic material provides one model, whereas an alternative is fractionation from primary mantle derived melts. Amphibole occurring in basaltic-andesitic and andesitic dyke rocks as fractionating cumulate phase extracted from lower crustal depth (6-7.5 kbar) is driving the magmas to peraluminous, corundum normative compositions, which are represented by tonalites forming most of the Adamello Batholith. Most primitive picrobasaltic dykes have a slightly steepened chondrite normalized rare earth elements (REE) pattern and the increased enrichment of light-REE (LREE) for andesites and dacites can be explained by the fractional crystallization model originating from a picrobasalt, taking the changing fractionating phase assemblage and temperature into account. The injection of hot basaltic magma (~1050°C) in a closely spaced dyke swarm increases the surface of the contact to the mainly tonalitic wallrock. Such a setting induces partial melting of the wall rock and selective assimilation. Partial melting of the tonalite host is further expressed through intrusion breccias from basaltic dykes. Heat conduction models with instantaneous magma injection for such a dyke swarm geometry can explain features of partial melting observed in the field. Geochemical data of minerals and bulk rock further underline the selective or bulk assimilation of the tonalite host rock at upper crustal levels (~2-3 kbar), in particular with regard to light ion lithophile elements (LILE) such as Sr, Ba and Rb. Primitive picrobasalts carry an immiscible felsic assimilant as enclaves that bring along refractory rutile and zircon with textures typically found in oceanic plagiogranites or high pressure/low-temperature metamorphic rocks in general. U-Pb data implies a lower Cretaceous age for zircon not yet described as assimilant in Eocene to Oligocene magmatic rocks of the Central Southern Alps. The distribution of post-plutonic dykes in large batholiths such as the Adamello is one of the key features for understanding the regional stress field during the post-batholith emplacement cooling history. The emplacement of the regional dyke swarm covering the southern part of the Adamello massif was associated with consistent left lateral strike-slip movement along magma dilatation planes, leading to en echelon segmentation of dykes. Through the dilation by magma of pre-existing weaknesses and cracks in an otherwise uniform host rock, the dyke propagation and according orientation in the horizontal plane adjusted continuously perpendicular to least compressive remote stress σ3, resulting in an inferred rotation of the remote principal stress field. Les magmas issus des zones de subduction contribuent substantiellement à la formation de la croûte continentale. Les plutons tonalitiques et granitiques représentent, en effet, une partie importante de la croûte continentale. Des magmas primaires produits dans le 'mantle wedge ', partie du manteau se trouvant au-dessus de la plaque plongeante dans des zones de subduction, migrent à travers le manteau puis la croûte. Pendant ce transfert, le magma peut s'accumuler dans des réservoirs intermédiaires à différentes profondeurs. Le stockage de magma dans ces réservoirs engendre, d'une part, la différentiation des magmas par cristallisation fractionnée et, d'autre part, une fusion partielle la croûte continentale préexistante associée au transfert de la chaleur des magmas vers l'encaissant. Ces liquides magmatiques issus de la croûte peuvent, ensuite, se mélanger avec des magmas primaires. Le transport du magma dans la croûte implique notamment un flux de magma à travers différentes fractures recoupant les roches encaissantes élastiques. Au cours de ce processus de migration, des cumulats de cristaux ou des agrégats de cristaux encore non-solidifiés, peuvent être recyclés et réactivés pour être transportés à des niveaux supérieures de la croûte. Le terrain d'étude est situé dans le massif d'Adamello. Celui-ci est composé de plusieurs plutons mis en place entre 42 et 29 millions d'années. Dans une phase tardive de l'activité magmatique liée à ce batholite, une série de filons de composition variable allant de picrobasalte à des compositions dacitiques s'est mise en place la partie sud du massif. Deux modèles sont proposés dans la littérature, pour expliquer la formation des magmas dacitiques caractérisés par des compositions peralumineux (i.e. à corindon normatif). Le premier modèle propose que ces magmas soient issus de la fusion de matériel mafique et pélitique présent dans la partie inférieur de la croûte, alors que le deuxième modèle suggère une évolution par cristallisation fractionnée à partir de liquides primaires issus du manteau. Un modèle de cristallisation fractionnée a pu être développé pour expliquer l'évolution des filons de l'Adamello. Ce modèle explique la formation des filons dacitiques par la cristallisation fractionnée de 17% olivine, 2% spinelle riche en Cr, 18% clinopyroxène, 41% amphibole, 4% plagioclase et 0.1% magnetite à partir de liquide de compositions picrobasaltiques. Ce modèle prend en considération les contraintes pétrologiques déduites de l'observation des différents filons ainsi que du champ de stabilité des différentes phases en fonction de la température. Ces roches montrent une évolution géochimique similaire aux données expérimentales simulant la cristallisation fractionnée de magmas évoluant à des niveaux inférieurs de la croûte (7-10 kbar). Le modèle montre, en particulier, le rôle prépondérant de l'amphibole, une phase qui contrôle en particulier le caractère peralumineux des magmas différentiés ainsi que leurs compositions en éléments en traces. Des phénomènes de fusion partielle de l'encaissant tonalitique lors de la mise en place de _lons mafiques sont observée sur le terrain. L'injection du magma basaltique chaud (~1050°C) sous forme de filons rapprochés augmente la surface du contact avec l'encaissante tonalitique. Une telle situation produit la fusion partielle des roches encaissantes nécessaire à l'incorporation d'enclaves mafiques observés au sein des tonalites. Pour comprendre les conditions nécessaires pour la fusion partielle des roches encaissantes, des modèles de conduction thermique pour une injection simultanée d'une série de filons ont été développées. Des données géochimiques sur les minéraux et sur les roches totales soulignent qu'au niveau supérieur de la croûte, l'assimilation sélective ou totale de l'encaissante tonalitique modifie la composition du liquide primaire pour les éléments lithophiles tel que le Sr, Ba et Rb. Un autre aspect important concernant la pétrologie des filons de l'Adamello est la présence d'enclaves felsiques dans les filons les plus primitifs. Ces enclaves montrent, en particulier, des textures proches de celles rencontrées dans des plagiogranites océaniques ou dans des roches métamorphiques de haute pression/basse température. Ces enclaves contiennent du zircon et du rutile. La datations de ces zircons à l'aide du géochronomètre U-Pb indique un âge Crétacé inférieur. Cet âge est important, car aucune roche de cet âge n'a été considérée comme un assimilant potentiel pour des roches magmatiques d'âge Eocène à Oligocène dans les Alpes Sud Centrales. La réparation spatiale des filons post-plutoniques dans des grands batholites tel que l'Adamello, est une caractéristique clé pour la compréhension des champs de contraintes lors du refroidissement du batholite. L'orientation des filons va, en particulier, indiqué la contrainte minimal au sein des roches encaissante. La mise en place de la série de filon recoupant la partie Sud du massif de l'Adamello est associée à un décrochement senestre, un décrochement que l'on peut lié aux contraintes tectoniques régionales auxquelles s'ajoutent l'effet de la dilatation produite par la mise en place du batholite lui-même. Ce décrochement senestre produit une segmentation en échelon des filons.

X-ray photoelectron spectroscopy of oxygen adsorbates on Al(111): Theory experiment

We present the result of polar angle resolved x¿ray photoemission spectroscopy on Al(111)/O and cluster calculations of the O(1s) binding energy (BE) for various model situations. In the experimental data two O(1s) peaks are observed, separated by 1.3 eV. The angular behavior (depth¿resolution) could indicate that the lower BE peak is associated with an O atom under the surface, and the higher BE peak with an O atom above the surface. Equally, it could indicate oxygen islands on the surface where the perimeter atoms have a higher O(1s) BE than the interior atoms. The cluster calculations show that the former interpretation cannot be correct, since an O ads below the surface has a higher calculated O(1s) BE than one above. Cluster calculations simulating oxygen islands are, however, consistent with the experimental data.

Analysis of joint laxity after total ankle arthroplasty: cadaver study.

BACKGROUND: Clinical results of total ankle arthroplasty with early designs were disappointing. Recently-developed ankle prostheses have good mid-term results; however, limited information is available regarding effects of total ankle arthroplasty on ankle laxity. METHODS: Eight cadaveric lower extremities were tested with a custom device which enabled measurement of multi-axial forces, moments, and displacement during applied axial, shear, and rotational loading. Tests consisted of anterior-posterior and medial-lateral translation and internal-external rotation of the talus relative to the tibia during axial loads on the tibia simulating body weight (700 N) and an unloaded condition (5 N). Tests were performed in neutral, dorsiflexion, and plantarflexion. Laxity was determined for the intact ankle, and following insertion of an unconstrained total ankle implant, comparing load-displacement curve. FINDINGS: Laxity after total ankle arthroplasty did not approximate the normal ankle in most conditions tested. Displacement was significantly greater for total ankle arthroplasty in both posterior and lateral translation, and internal rotation, with 5 N axial loading, and anterior-posterior, medial-lateral translation, and internal-external rotation for 700 N axial loading. For the 700 N axial load condition, in the neutral ankle position, total anterior-posterior translation averaged 0.4 mm (SD 0.2 mm), but 6.0 mm (SD 1.5 mm) after total ankle arthroplasty (P<0.01). This study demonstrated more laxity in the replaced ankle than normal ankle for both unloaded and 700 N axially loaded conditions. INTERPRETATION: These data indicate the increased responsibility of the ligaments for ankle laxity after total ankle arthroplasty and suggest the importance of meticulous ligament reconstruction with total ankle arthroplasty operations.

Performance of correlation functionals in ab initio chemisorption cluster-model calculations: Alkali metals on Si(111)

The performance of different correlation functionals has been tested for alkali metals, Li to Cs, interacting with cluster models simulating different active sites of the Si(111) surface. In all cases, the ab initio Hartree-Fock density has been obtained and used as a starting point. The electronic correlation energy is then introduced as an a posteriori correction to the Hartree-Fock energy using different correlation functionals. By making use of the ionic nature of the interaction and of different dissociation limits we have been able to prove that all functionals tested introduce the right correlation energy, although to a different extent. Hence, correlation functionals appear as an effective and easy way to introduce electronic correlation in the ab initio Hartree-Fock description of the chemisorption bond in complex systems where conventional configuration interaction techniques cannot be used. However, the calculated energies may differ by some tens of eV. Therefore, these methods can be employed to get a qualitative idea of how important correlation effects are, but they have some limitations if accurate binding energies are to be obtained.

Prothèses totales de genoux. Analyse clinique et numérique des micromouvements de l'implant tibial [Total knee prosthesis. Clinical and numerical study of micromovements of the tibial implant].

INTRODUCTION: The importance of the micromovements in the mechanism of aseptic loosening is clinically difficult to evaluate. To complete the analysis of a series of total knee arthroplasties (TKA), we used a tridimensional numerical model to study the micromovements of the tibial implant. MATERIAL AND METHODS: Fifty one patients (with 57 cemented Porous Coated Anatomic TKAs) were reviewed (mean follow-up 4.5 year). Radiolucency at the tibial bone-cement interface was sought on the AP radiographs and divided in 7 areas. The distribution of the radiolucency was then correlated with the axis of the lower limb as measured on the orthoradiograms. The tridimensional numerical model is based on the finite element method. It allowed the measurement of the cemented prosthetic tibial implant's displacements and the micromovements generated at bone-ciment interface. A total load (2000 Newton) was applied at first vertically and asymetrically on the tibial plateau, thereby simulating an axial deviation of the lower limbs. The vector's posterior inclination then permitted the addition of a tangential component to the axial load. This type of effort is generated by complex biomechanical phenomena such as knee flexion. RESULTS: 81 per cent of the 57 knees had a radiolucent line of at least 1 mm, at one or more of the tibial cement-epiphysis jonctional areas. The distribution of these lucent lines showed that they came out more frequently at the periphery of the implant. The lucent lines appeared most often under the unloaded margin of the tibial plateau, when axial deviation of lower limbs was present. Numerical simulations showed that asymetrical loading on the tibial plateau induced a subsidence of the loaded margin (0-100 microns) and lifting off at the opposite border (0-70 microns). The postero-anterior tangential component induced an anterior displacement of the tibial implant (160-220 microns), and horizontal micromovements with non homogenous distribution at the bone-ciment interface (28-54 microns). DISCUSSION: Comparison of clinical and numerical results showed a relation between the development of radiolucent lines and the unloading of the tibial implant's margin. The deleterious effect of lower limbs' axial deviation is thereby proven. The irregular distribution of lucent lines under the tibial plateau was similar of the micromovements' repartition at the bone-cement interface when tangential forces were present. A causative relation between the two phenomenaes could not however be established. Numerical simulation is a truly useful method of study; it permits to calculate micromovements which are relative, non homogenous and of very low amplitude. However, comparative clinical studies remain as essential to ensure the credibility of results.