Measuring ecological niche overlap from occurrence and spatial environmental data

1. Aim - Concerns over how global change will influence species distributions, in conjunction with increased emphasis on understanding niche dynamics in evolutionary and community contexts, highlight the growing need for robust methods to quantify niche differences between or within taxa. We propose a statistical framework to describe and compare environmental niches from occurrence and spatial environmental data.¦2. Location - Europe, North America, South America¦3. Methods - The framework applies kernel smoothers to densities of species occurrence in gridded environmental space to calculate metrics of niche overlap and test hypotheses regarding niche conservatism. We use this framework and simulated species with predefined distributions and amounts of niche overlap to evaluate several ordination and species distribution modeling techniques for quantifying niche overlap. We illustrate the approach with data on two well-studied invasive species.¦4. Results - We show that niche overlap can be accurately detected with the framework when variables driving the distributions are known. The method is robust to known and previously undocumented biases related to the dependence of species occurrences on the frequency of environmental conditions that occur across geographic space. The use of a kernel smoother makes the process of moving from geographical space to multivariate environmental space independent of both sampling effort and arbitrary choice of resolution in environmental space. However, the use of ordination and species distribution model techniques for selecting, combining and weighting variables on which niche overlap is calculated provide contrasting results.¦5. Main conclusions - The framework meets the increasing need for robust methods to quantify niche differences. It is appropriate to study niche differences between species, subspecies or intraspecific lineages that differ in their geographical distributions. Alternatively, it can be used to measure the degree to which the environmental niche of a species or intraspecific lineage has changed over time.

Scale-dependent adaptive evolution and morphological convergence to climatic niche in Californian eriogonoids (Polygonaceae)

Aim Macroevolutionary patterns and processes change substantially depending on levels of taxonomic and ecological organization, and the resolution of environmental and spatial variability. In comparative methods, the resolution of environmental and spatial variability often defines the number of selective regimes used to test whether phenotypic characteristics are adaptively correlated with the environment. Here, we examine how investigator choice of the number of selective regimes, determined by varying the resolution of among-species variability in the species climatic niche (hereafter called ecological scale'), influences trait morphological diversification among Eriogonoideae species. We assess whether adaptive or neutral processes drive the evolution of several morphological traits in these species. Location South-western North America. Methods We applied a phylogenetic framework of three evolutionary models to four morphological traits and the climatic niches of Eriogonoideae (in the buckwheat family, Polygonaceae). We tested whether morphological traits evolve in relation to climate by adaptive or neutral process, and whether the resulting patterns of morphological variability are conserved or convergent across the clade. We inspected adaptive models of evolution under different levels of resolution of among-species variability of the climatic niche. Results We show that morphological traits and climate niches of Eriogonoideae species are not phylogenetically conserved. Further, adaptive evolution of phenotypic traits is specific to climatic niche occupancy across this clade. Finally, the likely evolutionary process and the level of detectable niche conservatism change depending on the resolution of environmental variability of the climatic niche. Main conclusions Our study demonstrates the need to consider both the resolution of environmental variability and alternative evolutionary models to understand the morphological diversification that accompanies divergent adaptive evolution of lineages to climatic conditions.

Time-lapse and probabilistic inversion strategies for plane-wave electromagnetic data.

L'utilisation efficace des systèmes géothermaux, la séquestration du CO2 pour limiter le changement climatique et la prévention de l'intrusion d'eau salée dans les aquifères costaux ne sont que quelques exemples qui démontrent notre besoin en technologies nouvelles pour suivre l'évolution des processus souterrains à partir de la surface. Un défi majeur est d'assurer la caractérisation et l'optimisation des performances de ces technologies à différentes échelles spatiales et temporelles. Les méthodes électromagnétiques (EM) d'ondes planes sont sensibles à la conductivité électrique du sous-sol et, par conséquent, à la conductivité électrique des fluides saturant la roche, à la présence de fractures connectées, à la température et aux matériaux géologiques. Ces méthodes sont régies par des équations valides sur de larges gammes de fréquences, permettant détudier de manières analogues des processus allant de quelques mètres sous la surface jusqu'à plusieurs kilomètres de profondeur. Néanmoins, ces méthodes sont soumises à une perte de résolution avec la profondeur à cause des propriétés diffusives du champ électromagnétique. Pour cette raison, l'estimation des modèles du sous-sol par ces méthodes doit prendre en compte des informations a priori afin de contraindre les modèles autant que possible et de permettre la quantification des incertitudes de ces modèles de façon appropriée. Dans la présente thèse, je développe des approches permettant la caractérisation statique et dynamique du sous-sol à l'aide d'ondes EM planes. Dans une première partie, je présente une approche déterministe permettant de réaliser des inversions répétées dans le temps (time-lapse) de données d'ondes EM planes en deux dimensions. Cette stratégie est basée sur l'incorporation dans l'algorithme d'informations a priori en fonction des changements du modèle de conductivité électrique attendus. Ceci est réalisé en intégrant une régularisation stochastique et des contraintes flexibles par rapport à la gamme des changements attendus en utilisant les multiplicateurs de Lagrange. J'utilise des normes différentes de la norme l2 pour contraindre la structure du modèle et obtenir des transitions abruptes entre les régions du model qui subissent des changements dans le temps et celles qui n'en subissent pas. Aussi, j'incorpore une stratégie afin d'éliminer les erreurs systématiques de données time-lapse. Ce travail a mis en évidence l'amélioration de la caractérisation des changements temporels par rapport aux approches classiques qui réalisent des inversions indépendantes à chaque pas de temps et comparent les modèles. Dans la seconde partie de cette thèse, j'adopte un formalisme bayésien et je teste la possibilité de quantifier les incertitudes sur les paramètres du modèle dans l'inversion d'ondes EM planes. Pour ce faire, je présente une stratégie d'inversion probabiliste basée sur des pixels à deux dimensions pour des inversions de données d'ondes EM planes et de tomographies de résistivité électrique (ERT) séparées et jointes. Je compare les incertitudes des paramètres du modèle en considérant différents types d'information a priori sur la structure du modèle et différentes fonctions de vraisemblance pour décrire les erreurs sur les données. Les résultats indiquent que la régularisation du modèle est nécessaire lorsqu'on a à faire à un large nombre de paramètres car cela permet d'accélérer la convergence des chaînes et d'obtenir des modèles plus réalistes. Cependent, ces contraintes mènent à des incertitudes d'estimations plus faibles, ce qui implique des distributions a posteriori qui ne contiennent pas le vrai modèledans les régions ou` la méthode présente une sensibilité limitée. Cette situation peut être améliorée en combinant des méthodes d'ondes EM planes avec d'autres méthodes complémentaires telles que l'ERT. De plus, je montre que le poids de régularisation des paramètres et l'écart-type des erreurs sur les données peuvent être retrouvés par une inversion probabiliste. Finalement, j'évalue la possibilité de caractériser une distribution tridimensionnelle d'un panache de traceur salin injecté dans le sous-sol en réalisant une inversion probabiliste time-lapse tridimensionnelle d'ondes EM planes. Etant donné que les inversions probabilistes sont très coûteuses en temps de calcul lorsque l'espace des paramètres présente une grande dimension, je propose une stratégie de réduction du modèle ou` les coefficients de décomposition des moments de Legendre du panache de traceur injecté ainsi que sa position sont estimés. Pour ce faire, un modèle de résistivité de base est nécessaire. Il peut être obtenu avant l'expérience time-lapse. Un test synthétique montre que la méthodologie marche bien quand le modèle de résistivité de base est caractérisé correctement. Cette méthodologie est aussi appliquée à un test de trac¸age par injection d'une solution saline et d'acides réalisé dans un système géothermal en Australie, puis comparée à une inversion time-lapse tridimensionnelle réalisée selon une approche déterministe. L'inversion probabiliste permet de mieux contraindre le panache du traceur salin gr^ace à la grande quantité d'informations a priori incluse dans l'algorithme. Néanmoins, les changements de conductivités nécessaires pour expliquer les changements observés dans les données sont plus grands que ce qu'expliquent notre connaissance actuelle des phénomenès physiques. Ce problème peut être lié à la qualité limitée du modèle de résistivité de base utilisé, indiquant ainsi que des efforts plus grands devront être fournis dans le futur pour obtenir des modèles de base de bonne qualité avant de réaliser des expériences dynamiques. Les études décrites dans cette thèse montrent que les méthodes d'ondes EM planes sont très utiles pour caractériser et suivre les variations temporelles du sous-sol sur de larges échelles. Les présentes approches améliorent l'évaluation des modèles obtenus, autant en termes d'incorporation d'informations a priori, qu'en termes de quantification d'incertitudes a posteriori. De plus, les stratégies développées peuvent être appliquées à d'autres méthodes géophysiques, et offrent une grande flexibilité pour l'incorporation d'informations additionnelles lorsqu'elles sont disponibles. -- The efficient use of geothermal systems, the sequestration of CO2 to mitigate climate change, and the prevention of seawater intrusion in coastal aquifers are only some examples that demonstrate the need for novel technologies to monitor subsurface processes from the surface. A main challenge is to assure optimal performance of such technologies at different temporal and spatial scales. Plane-wave electromagnetic (EM) methods are sensitive to subsurface electrical conductivity and consequently to fluid conductivity, fracture connectivity, temperature, and rock mineralogy. These methods have governing equations that are the same over a large range of frequencies, thus allowing to study in an analogous manner processes on scales ranging from few meters close to the surface down to several hundreds of kilometers depth. Unfortunately, they suffer from a significant resolution loss with depth due to the diffusive nature of the electromagnetic fields. Therefore, estimations of subsurface models that use these methods should incorporate a priori information to better constrain the models, and provide appropriate measures of model uncertainty. During my thesis, I have developed approaches to improve the static and dynamic characterization of the subsurface with plane-wave EM methods. In the first part of this thesis, I present a two-dimensional deterministic approach to perform time-lapse inversion of plane-wave EM data. The strategy is based on the incorporation of prior information into the inversion algorithm regarding the expected temporal changes in electrical conductivity. This is done by incorporating a flexible stochastic regularization and constraints regarding the expected ranges of the changes by using Lagrange multipliers. I use non-l2 norms to penalize the model update in order to obtain sharp transitions between regions that experience temporal changes and regions that do not. I also incorporate a time-lapse differencing strategy to remove systematic errors in the time-lapse inversion. This work presents improvements in the characterization of temporal changes with respect to the classical approach of performing separate inversions and computing differences between the models. In the second part of this thesis, I adopt a Bayesian framework and use Markov chain Monte Carlo (MCMC) simulations to quantify model parameter uncertainty in plane-wave EM inversion. For this purpose, I present a two-dimensional pixel-based probabilistic inversion strategy for separate and joint inversions of plane-wave EM and electrical resistivity tomography (ERT) data. I compare the uncertainties of the model parameters when considering different types of prior information on the model structure and different likelihood functions to describe the data errors. The results indicate that model regularization is necessary when dealing with a large number of model parameters because it helps to accelerate the convergence of the chains and leads to more realistic models. These constraints also lead to smaller uncertainty estimates, which imply posterior distributions that do not include the true underlying model in regions where the method has limited sensitivity. This situation can be improved by combining planewave EM methods with complimentary geophysical methods such as ERT. In addition, I show that an appropriate regularization weight and the standard deviation of the data errors can be retrieved by the MCMC inversion. Finally, I evaluate the possibility of characterizing the three-dimensional distribution of an injected water plume by performing three-dimensional time-lapse MCMC inversion of planewave EM data. Since MCMC inversion involves a significant computational burden in high parameter dimensions, I propose a model reduction strategy where the coefficients of a Legendre moment decomposition of the injected water plume and its location are estimated. For this purpose, a base resistivity model is needed which is obtained prior to the time-lapse experiment. A synthetic test shows that the methodology works well when the base resistivity model is correctly characterized. The methodology is also applied to an injection experiment performed in a geothermal system in Australia, and compared to a three-dimensional time-lapse inversion performed within a deterministic framework. The MCMC inversion better constrains the water plumes due to the larger amount of prior information that is included in the algorithm. The conductivity changes needed to explain the time-lapse data are much larger than what is physically possible based on present day understandings. This issue may be related to the base resistivity model used, therefore indicating that more efforts should be given to obtain high-quality base models prior to dynamic experiments. The studies described herein give clear evidence that plane-wave EM methods are useful to characterize and monitor the subsurface at a wide range of scales. The presented approaches contribute to an improved appraisal of the obtained models, both in terms of the incorporation of prior information in the algorithms and the posterior uncertainty quantification. In addition, the developed strategies can be applied to other geophysical methods, and offer great flexibility to incorporate additional information when available.

Safe and effective variability-a criterion for dose individualization.

BACKGROUND: : A primary goal of clinical pharmacology is to understand the factors that determine the dose-effect relationship and to use this knowledge to individualize drug dose. METHODS: : A principle-based criterion is proposed for deciding among alternative individualization methods. RESULTS: : Safe and effective variability defines the maximum acceptable population variability in drug concentration around the population average. CONCLUSIONS: : A decision on whether patient covariates alone are sufficient, or whether therapeutic drug monitoring in combination with target concentration intervention is needed, can be made by comparing the remaining population variability after a particular dosing method with the safe and effective variability.

Sleep function: current questions and new approaches.

The mammalian brain oscillates through three distinct global activity states: wakefulness, non-rapid eye movement (NREM) sleep and REM sleep. The regulation and function of these 'vigilance' or 'behavioural' states can be investigated over a broad range of temporal and spatial scales and at different levels of functional organization, i.e. from gene expression to memory, in single neurons, cortical columns or the whole brain and organism. We summarize some basic questions that have arisen from recent approaches in the quest for the functions of sleep. Whereas traditionally sleep was viewed to be regulated through top-down control mechanisms, recent approaches have emphasized that sleep is emerging locally and regulated in a use-dependent (homeostatic) manner. Traditional markers of sleep homeostasis, such as the electroencephalogram slow-wave activity, have been linked to changes in connectivity and plasticity in local neuronal networks. Thus waking experience-induced local network changes may be sensed by the sleep homeostatic process and used to mediate sleep-dependent events, benefiting network stabilization and memory consolidation. Although many questions remain unanswered, the available data suggest that sleep function will best be understood by an analysis which integrates sleep's many functional levels with its local homeostatic regulation.

Activation and inhibition of bimanual movements in school-aged children

The development of motor activation and inhibition was compared in 6-to-12 year-olds. Children had to initiate or stop the externally paced movements of one hand, while maintaining that of the other hand. The time needed to perform the switching task (RT) and the spatio-temporal variables show different agerelated evolutions depending on the coordination pattern (inor anti-phase) and the type of transition (activation, selective inhibition, non selective inhibition) required. In the anti-phase mode, activation perturbs the younger subjects' responses while temporal and spatial stabilities transiently decrease around 9 years when activating in the in-phase mode. Aged-related changes differed between inhibition and activation in the antiphase mode, suggesting either the involvement of distinct neural networks or the existence of a single network that is reorganized. In contrast, stopping or adding one hand in the in-phase mode shows similar aged-related improvement. We suggest that selectively stopping or activating one arm during symmetrical coordination rely on the two faces of a common processing in which activation could be the release of inhibition

Unravelling the interaction between tectonic and sedimentary processes during lithospheric thinning in the Alpine Tethys margins

The discovery of exhumed continental mantle and hyper-extended crust in present-day magma-poor rifted margins is at the origin of a paradigm shift within the research field of deep-water rifted margins. It opened new questions about the strain history of rifted margins and the nature and composition of sedimentary, crustal and mantle rocks in rifted margins. Thanks to the benefit of more than one century of work in the Alps and access to world-class outcrops preserving the primary relationships between sediments and crustal and mantle rocks from the fossil Alpine Tethys margins, it is possible to link the subsidence history and syn-rift sedimentary evolution with the strain distribution observed in the crust and mantle rocks exposed in the distal rifted margins. In this paper, we will focus on the transition from early to late rifting that is associated with considerable crustal thinning and a reorganization of the rift system. Crustal thinning is at the origin of a major change in the style of deformation from high-angle to low-angle normal faulting which controls basin-architecture, sedimentary sources and processes and the nature of basement rocks exhumed along the detachment faults in the distal margin. Stratigraphic and isotopic ages indicate that this major change occurred in late Sinemurian time, involving a shift of the syn-rift sedimentation toward the distal domain associated with a major reorganization of the crustal structure with exhumation of lower and middle crust. These changes may be triggered by mantle processes, as indicated by the infiltration of MOR-type magmas in the lithospheric mantle, and the uplift of the Brianconnais domain. Thinning and exhumation of the crust and lithosphere also resulted in the creation of new paleogeographic domains, the Proto Valais and Liguria-Piemonte domains. These basins show a complex, 3D temporal and spatial evolution that might have evolved, at least in the case of the Liguria-Piemonte basin, in the formation of an embryonic oceanic crust. The re-interpretation of the rift evolution and the architecture of the distal rifted margins in the Alps have important implications for the understanding of rifted margins worldwide, but also for the paleogeographic reconstruction of the Alpine domain and its subsequent Alpine compressional overprint.

What drives invasibility? A multi-model inference test and spatial modelling of alien plant species richness patterns in northern Portugal

Understanding and anticipating biological invasions can focus either on traits that favour species invasiveness or on features of the receiving communities, habitats or landscapes that promote their invasibility. Here, we address invasibility at the regional scale, testing whether some habitats and landscapes are more invasible than others by fitting models that relate alien plant species richness to various environmental predictors. We use a multi-model information-theoretic approach to assess invasibility by modelling spatial and ecological patterns of alien invasion in landscape mosaics and testing competing hypotheses of environmental factors that may control invasibility. Because invasibility may be mediated by particular characteristics of invasiveness, we classified alien species according to their C-S-R plant strategies. We illustrate this approach with a set of 86 alien species in Northern Portugal. We first focus on predictors influencing species richness and expressing invasibility and then evaluate whether distinct plant strategies respond to the same or different groups of environmental predictors. We confirmed climate as a primary determinant of alien invasions and as a primary environmental gradient determining landscape invasibility. The effects of secondary gradients were detected only when the area was sub-sampled according to predictions based on the primary gradient. Then, multiple predictor types influenced patterns of alien species richness, with some types (landscape composition, topography and fire regime) prevailing over others. Alien species richness responded most strongly to extreme land management regimes, suggesting that intermediate disturbance induces biotic resistance by favouring native species richness. Land-use intensification facilitated alien invasion, whereas conservation areas hosted few invaders, highlighting the importance of ecosystem stability in preventing invasions. Plants with different strategies exhibited different responses to environmental gradients, particularly when the variations of the primary gradient were narrowed by sub-sampling. Such differential responses of plant strategies suggest using distinct control and eradication approaches for different areas and alien plant groups.

Analysis and spatial modeling of the indoor radon Swiss data

The present research deals with an important public health threat, which is the pollution created by radon gas accumulation inside dwellings. The spatial modeling of indoor radon in Switzerland is particularly complex and challenging because of many influencing factors that should be taken into account. Indoor radon data analysis must be addressed from both a statistical and a spatial point of view. As a multivariate process, it was important at first to define the influence of each factor. In particular, it was important to define the influence of geology as being closely associated to indoor radon. This association was indeed observed for the Swiss data but not probed to be the sole determinant for the spatial modeling. The statistical analysis of data, both at univariate and multivariate level, was followed by an exploratory spatial analysis. Many tools proposed in the literature were tested and adapted, including fractality, declustering and moving windows methods. The use of Quan-tité Morisita Index (QMI) as a procedure to evaluate data clustering in function of the radon level was proposed. The existing methods of declustering were revised and applied in an attempt to approach the global histogram parameters. The exploratory phase comes along with the definition of multiple scales of interest for indoor radon mapping in Switzerland. The analysis was done with a top-to-down resolution approach, from regional to local lev¬els in order to find the appropriate scales for modeling. In this sense, data partition was optimized in order to cope with stationary conditions of geostatistical models. Common methods of spatial modeling such as Κ Nearest Neighbors (KNN), variography and General Regression Neural Networks (GRNN) were proposed as exploratory tools. In the following section, different spatial interpolation methods were applied for a par-ticular dataset. A bottom to top method complexity approach was adopted and the results were analyzed together in order to find common definitions of continuity and neighborhood parameters. Additionally, a data filter based on cross-validation was tested with the purpose of reducing noise at local scale (the CVMF). At the end of the chapter, a series of test for data consistency and methods robustness were performed. This lead to conclude about the importance of data splitting and the limitation of generalization methods for reproducing statistical distributions. The last section was dedicated to modeling methods with probabilistic interpretations. Data transformation and simulations thus allowed the use of multigaussian models and helped take the indoor radon pollution data uncertainty into consideration. The catego-rization transform was presented as a solution for extreme values modeling through clas-sification. Simulation scenarios were proposed, including an alternative proposal for the reproduction of the global histogram based on the sampling domain. The sequential Gaussian simulation (SGS) was presented as the method giving the most complete information, while classification performed in a more robust way. An error measure was defined in relation to the decision function for data classification hardening. Within the classification methods, probabilistic neural networks (PNN) show to be better adapted for modeling of high threshold categorization and for automation. Support vector machines (SVM) on the contrary performed well under balanced category conditions. In general, it was concluded that a particular prediction or estimation method is not better under all conditions of scale and neighborhood definitions. Simulations should be the basis, while other methods can provide complementary information to accomplish an efficient indoor radon decision making.

Reticulocyte dynamic and hemoglobin variability in hemodialysis patients treated with Darbepoetin alfa and C.E.R.A.: a randomized controlled trial.

BACKGROUND: In a simulation based on a pharmacokinetic model we demonstrated that increasing the erythropoiesis stimulating agents (ESAs) half-life or shortening their administration interval decreases hemoglobin variability. The benefit of reducing the administration interval was however lessened by the variability induced by more frequent dosage adjustments. The purpose of this study was to analyze the reticulocyte and hemoglobin kinetics and variability under different ESAs and administration intervals in a collective of chronic hemodialysis patients. METHODS: The study was designed as an open-label, randomized, four-period cross-over investigation, including 30 patients under chronic hemodialysis at the regional hospital of Locarno (Switzerland) in February 2010 and lasting 2 years. Four subcutaneous treatment strategies (C.E.R.A. every 4 weeks Q4W and every 2 weeks Q2W, Darbepoetin alfa Q4W and Q2W) were compared with each other. The mean square successive difference of hemoglobin, reticulocyte count and ESAs dose was used to quantify variability. We distinguished a short- and a long-term variability based respectively on the weekly and monthly successive difference. RESULTS: No difference was found in the mean values of biological parameters (hemoglobin, reticulocytes, and ferritin) between the 4 strategies. ESAs type did not affect hemoglobin and reticulocyte variability, but C.E.R.A induced a more sustained reticulocytes response over time and increased the risk of hemoglobin overshooting (OR 2.7, p = 0.01). Shortening the administration interval lessened the amplitude of reticulocyte count fluctuations but resulted in more frequent ESAs dose adjustments and in amplified reticulocyte and hemoglobin variability. Q2W administration interval was however more favorable in terms of ESAs dose, allowing a 38% C.E.R.A. dose reduction, and no increase of Darbepoetin alfa. CONCLUSIONS: The reticulocyte dynamic was a more sensitive marker of time instability of the hemoglobin response under ESAs therapy. The ESAs administration interval had a greater impact on hemoglobin variability than the ESAs type. The more protracted reticulocyte response induced by C.E.R.A. could explain both, the observed higher risk of overshoot and the significant increase in efficacy when shortening its administration interval.Trial registrationClinicalTrials.gov NCT01666301.

Enzyme polymorphism and clinical variability of diseases: a study of diabetes mellitus.

We investigated possible relations among four common neonatal manifestations of diabetic pregnancy (macrosomia, hypoglycemia, hypocalcemia, jaundice) and four enzyme polymorphisms (PGM1, ADA, AK1, ACP1 in a sample of infants born of diabetic mothers. The pattern of associations observed between the two sets of variables is consistent with known differences in enzymatic activity within phenotypes of each system, suggesting that low enzymatic activity may have unfavorable effects on fetal development and on adaptability of the neonate to the extrauterine environment, Some of the polymorphic enzymes studied influence fetal growth in normal pregnancy as well. Analysis of relations between genetic polymorphisms and the clinical pattern of common diseases may provide a better understanding of the genetic basis of the clinical variability of diseases within and between human populations.

How do the different components of episodic memory develop? Role of executive functions and short-term feature-binding abilities.

This study investigated the development of all 3 components of episodic memory (EM), as defined by Tulving, namely, core factual content, spatial context, and temporal context. To this end, a novel, ecologically valid test was administered to 109 participants aged 4-16 years. Results showed that each EM component develops at a different rate. Ability to memorize factual content emerges early, whereas context retrieval abilities continue to improve until adolescence, due to persistent encoding difficulties (isolated by comparing results on free recall and recognition tasks). Exploration of links with other cognitive functions revealed that short-term feature-binding abilities contribute to all EM components, and executive functions to temporal and spatial context, although ability to memorize temporal context is predicted mainly by age.

Inter- and intraobserver variability of measurements with Pentacam HR® in patients with mild to moderate keratoconus

Purpose: To evaluate inter- and intraobserver variability of indices crucial for detection of keratoconus progression derived from the Pentacam HR® (high-resolution) tomographer (OCULUS Optikgeräte GmbH, Wetzlar, Germany) in patients with mild to moderate keratoconus. Methods: Three repeated corneal topography measurements in the 25-picture mode by two independent observers were performed. The extent of variability across a large range of measurement parameters was analyzed including anterior and posterior corneal surface measurements, pachymetry values, corneal volume, anterior chamber volume and depth, and iridocorneal angle. The intraclass correlation coefficient (ICC) between and within each investigator was calculated to assess reproducibility and repeatability, respectively. Results: 31 eyes of 20 patients (mean age 31.6, SD ± 8.6) were included. Overall, the repeatability and reproducibility were excellent. The range of variability was reported by calculating the standard deviation of measurements. The detailed results are shown in Table 1. Conclusions: This study shows that the Pentacam HR® tomographer provides reliable measurements in patients with mild to moderate keratoconus. However, all parameters showed a certain range of variability. This should be taken into account when assessing keratoconus progression in order to distinguish true progression from variability in measurements. In addition, the excellent reproducibility suggests that the measurements can be reliably performed by different individuals from one visit to another.