Effects of sample size on the performance of species distribution models

A wide range of modelling algorithms is used by ecologists, conservation practitioners, and others to predict species ranges from point locality data. Unfortunately, the amount of data available is limited for many taxa and regions, making it essential to quantify the sensitivity of these algorithms to sample size. This is the first study to address this need by rigorously evaluating a broad suite of algorithms with independent presence-absence data from multiple species and regions. We evaluated predictions from 12 algorithms for 46 species (from six different regions of the world) at three sample sizes (100, 30, and 10 records). We used data from natural history collections to run the models, and evaluated the quality of model predictions with area under the receiver operating characteristic curve (AUC). With decreasing sample size, model accuracy decreased and variability increased across species and between models. Novel modelling methods that incorporate both interactions between predictor variables and complex response shapes (i.e. GBM, MARS-INT, BRUTO) performed better than most methods at large sample sizes but not at the smallest sample sizes. Other algorithms were much less sensitive to sample size, including an algorithm based on maximum entropy (MAXENT) that had among the best predictive power across all sample sizes. Relative to other algorithms, a distance metric algorithm (DOMAIN) and a genetic algorithm (OM-GARP) had intermediate performance at the largest sample size and among the best performance at the lowest sample size. No algorithm predicted consistently well with small sample size (n < 30) and this should encourage highly conservative use of predictions based on small sample size and restrict their use to exploratory modelling.

Spontaneous viral clearance, viral load, and genotype distribution of hepatitis C virus (HCV) in HIV-infected patients with anti-HCV antibodies in Europe.

BACKGROUND: Variables influencing serum hepatitis C virus (HCV) RNA levels and genotype distribution in individuals with human immunodeficiency virus (HIV) infection are not well known, nor are factors determining spontaneous clearance after exposure to HCV in this population. METHODS: All HCV antibody (Ab)-positive patients with HIV infection in the EuroSIDA cohort who had stored samples were tested for serum HCV RNA, and HCV genotyping was done for subjects with viremia. Logistic regression was used to identify variables associated with spontaneous HCV clearance and HCV genotype 1. RESULTS: Of 1940 HCV Ab-positive patients, 1496 (77%) were serum HCV RNA positive. Injection drug users (IDUs) were less likely to have spontaneously cleared HCV than were homosexual men (20% vs. 39%; adjusted odds ratio [aOR], 0.36 [95% confidence interval {CI}, 0.24-0.53]), whereas patients positive for hepatitis B surface antigen (HBsAg) were more likely to have spontaneously cleared HCV than were those negative for HBsAg (43% vs. 21%; aOR, 2.91 [95% CI, 1.94-4.38]). Of patients with HCV viremia, 786 (53%) carried HCV genotype 1, and 53 (4%), 440 (29%), and 217 (15%) carried HCV genotype 2, 3, and 4, respectively. A greater HCV RNA level was associated with a greater chance of being infected with HCV genotype 1 (aOR, 1.60 per 1 log higher [95% CI, 1.36-1.88]). CONCLUSIONS: More than three-quarters of the HIV- and HCV Ab-positive patients in EuroSIDA showed active HCV replication. Viremia was more frequent in IDUs and, conversely, was less common in HBsAg-positive patients. Of the patients with HCV viremia analyzed, 53% were found to carry HCV genotype 1, and this genotype was associated with greater serum HCV RNA levels.

A multi-temporal approach to model endangered species distribution in Europe: the case of the Eurasian otter in Italy

We analysed the relationship between changes in land cover patterns and the Eurasian otter occurrence over the course of about 20 years (1985-2006) using multi-temporal Species Distribution Models (SDMs). The study area includes five river catchments covering most of the otter's Italian range. Land cover and topographic data were used as proxies of the ecological requirements of the otter within a 300-m buffer around river courses. We used species presence, pseudo-absence data, and environmental predictors to build past (1985) and current (2006) SDMs by applying an ensemble procedure through the BIOMOD modelling package. The performance of each model was evaluated by measuring the area under the curve (AUC) of the receiver-operating characteristic (ROC). Multi-temporal analyses of species distribution and land cover maps were performed by comparing the maps produced for 1985 and 2006. The ensemble procedure provided a good overall modelling accuracy, revealing that elevation and slope affected the otter's distribution in the past; in contrast, land cover predictors, such as cultivations and forests, were more important in the present period. During the transition period, 20.5% of the area became suitable, with 76% of the new otter presence data being located in these newly available areas. The multi-temporal analysis suggested that the quality of otter habitat improved in the last 20 years owing to the expansion of forests and to the reduction of cultivated fields in riparian belts. The evidence presented here stresses the great potential of riverine habitat restoration and environmental management for the future expansion of the otter in Italy

Joergensenium rotatile n. gen., n. sp (Entactinaria, Radiolaria): its distribution in west Norwegian fjords

A new radiolarian genus and species, Joergensenium rotatile, is described. This species is restricted to recent sediments and plankton samples from the North Atlantic. Its recent distribution in the Norwegian Sea and West Norwegian fjords shows a strong affinity to the neritic province and reaches almost 2% in Hryangerfjord. This species is only known from late Glacial and Holocene sediments in the Nordic seas. This genus shows, however, a patchy stratigraphic distribution with its first occurrence in the south-west Pacific within Palcocene, in the Middle to Late Miocene from the Norwegian Sea, and in the Labrador Sea at the base of biozone NN 21. Two fjords are compared for the general radiolarian distribution, but with special emphasis on the occurrence of J. rotatile in both sediment and plankton.

Comparison of plantar pressure distribution in adolescent runners at low vs. high running velocity

This study aimed to compare foot plantar pressure distribution while jogging and running in highly trained adolescent runners. Eleven participants performed two constant-velocity running trials either at jogging (11.2 ± 0.9 km/h) or running (17.8 ± 1.4 km/h) pace on a treadmill. Contact area (CA in cm(2)), maximum force (F(max) in N), peak pressure (PP in kPa), contact time (CT in ms), and relative load (force time integral in each individual region divided by the force time integral for the total plantar foot surface, in %) were measured in nine regions of the right foot using an in-shoe plantar pressure device. Under the whole foot, CA, F(max) and PP were lower in jogging than in running (-1.2% [p<0.05], -12.3% [p<0.001] and -15.1% [p<0.01] respectively) whereas CT was higher (+20.1%; p<0.001). Interestingly, we found an increase in relative load under the medial and central forefoot regions while jogging (+6.7% and +3.7%, respectively; [p<0.05]), while the relative load under the lesser toes (-8.4%; p<0.05) was reduced. In order to prevent overloading of the metatarsals in adolescent runners, excessive mileage at jogging pace should be avoided.

The DNA damage response to adeno-associated virus : centrosome overduplication and induction of cell death independently of the spindle checkpoint

Summary: Adeno-associated virus type 2 (AAV2) is a small virus containing single-stranded DNA of approximately 4.7kb in size. Both ends of the viral genome are flanked with inverted terminal repeat sequences (ITRs), which serve as primers for viral replication. Previous work in our laboratory has shown that AAV2 DNA with ultraviolet radiation-generated crosslinks (UV-AAV2) provokes a DNA damage response in the host cell by mimicking a stalled replication fork. Infection of cells with UV-AAV2 leads to a p53-and Chk1-mediated cell cycle arrest at the G2/M border of the cell cycle. However, tumour cells lacking the tumour suppressor protein p53 cannot sustain this arrest and enter a prolonged impaired mitosis, the outcome of which is cell death. The aim of my thesis was to investigate how UV-inactivated AAV2 kilts p53-deficient cancer cells. I found that the UV-AAV2-induced DNA damage signalling induces centriole overduplication in infected cells. The virus is able to uncouple the centriole duplication cycle from the cell cycle, leading to amplified centrosome numbers. Chk1 colocalises with centrosomes in the infected cells and the centrosome overduplication is dependent on the presence of Chk1, as well as on the activities of ATR and Cdk kinases and on the G2 arrest. The UV-AAV2-induced DNA damage signalling inhibits the degradation of cyclin B 1 and securin by the anaphase promoting complex, suggesting that the spindle checkpoint is activated in these mitotic cells. Interference with the spindle checkpoint components Mad2 and BubR1 revealed that the UV-AAV2-provoked mitotic catastrophe occurs independently of spindle checkpoint function, This work shows that, in the p53 deficient cells, UV-AAV2 triggers mitotic catastrophe associated with a dramatic Chk1-dependent overduplication of centrioles and the consequent formation of multiple spindle poles in mitosis. Résumé Le virus associé à l'adénovirus type 2 (AAV2) est un petit virus contenant un simple brin d'ADN d'environ 4.7kb. Des expériences antérieures dans notre laboratoire ont montré que les liens intramoléculaires sur l'ADN de AAV2 provoqués paz l'irradiation aux ultraviolets (UV) ressemblent à une fourche de réplication bloquée, ce qui provoque une réponse aux dommages à l'ADN dans la cellule hôte. L'infection des cellules avec UV-AAV2 résulte en un arrêt du cycle cellulaire à la transition G2/M entraîné par les protéines ATR et Chk1. Cependant, les cellules tumorales auxquelles il manque le suppresseur de tumeur p53 ne peuvent pas tenir cet arrêt et entrent dans une mitose anormale et prolongée qui se terminera par la mort cellulaire. Le but de ma thèse était d'étudier comment l'AAV2 inactivé par l'irradiation UV tue les cellules cancéreuses n'ayant pas p53. Je montre ici que le signal de dommages à l'ADN induit par UV-AAV2 génère une surduplication des centrioles dans les cellules infectées. Le virus est capable de dissocier le cycle de duplication du centriole du cycle cellulaire ce qui crée un nombre amplifié de centrosomes. Chk1 est co-localisé avec le centrosome dans les cellules infectées et la swduplication du centrosome est dépendante de la présence de Chk1, de l'activité des kinases ATR et Cdk et de l'arrêt en G2 de la cellule. Le signal d'ADN endommagé induit par UV-AAV2 réprime la dégradation des protéines cycline B1 et securine par le complexe promoteur de l'anaphase (APC), ce qui suggère que le point de contrôle du fuseau mitotique est activé dans ces cellules en mitose. L'étude d'interférence avec des éléments du point de contrôle du fuseau mitotique, Mad2 et BubR1, a révélé que la catastrophe mitotique provoquée paz UV-AAV2 survient indépendamment du point de contrôle du fuseau mitotique. Ce travail montre que dans les cellules déficientes en p53, UV-AAV2 induit une catastrophe mitotique associée à une surduplication des centrioles dépendant de Chk1 et ayant pour conséquence dramatique la formation de multiples fuseaux mitotiques dans la cellule en mitose.

Fibrosarcoma-like lipomatous neoplasm: a reappraisal of so-called spindle cell liposarcoma defining a unique lipomatous tumor unrelated to other liposarcomas

The term "spindle cell liposarcoma" has been applied to liposarcomas (LPSs) composed predominantly or exclusively of spindled cells. These tumors have been considered variants of well-differentiated LPS (WDL), myxoid LPS, and spindle cell lipoma, suggesting that this is a heterogenous group of lesions. Using strict morphologic criteria and molecular and immunohistochemical analyses, we have identified a homogenous group of spindle cell lipomatous tumors, histologically and genetically distinct from other forms of LPS, which we have called "fibrosarcoma-like lipomatous neoplasm." Cases classified as "spindle cell LPS" or "low-grade LPS with spindle cell features" were reviewed. Final selection criteria included: (1) an exclusive low-grade spindle cell component resembling fibrosarcoma; (2) a mixture of bland fibroblastic cells resembling the preadipocyte and early-adipocyte stage of embryonic fat; and (3) molecular-genetic analysis that excluded other forms of lipomatous tumors. Of the initial 25 cases identified, comparative genomic hybridization (CGH) was uninformative in 2 cases; 5 were reclassified as WDL on the basis of molecular data (MDM2 amplification) and 6 as spindle cell lipoma (CGH profiles with a few gains and losses including a constant loss of chromosome 13 and frequent losses of chromosomes 16 and 6). The 12 remaining cases showed flat CGH profiles; of these cases, 11 were negative for DDIT3 gene rearrangements, and 1 result was uninterpretable. Patients ranged in age from 15 to 82 years (mean 50 y); male patients were affected slightly more often (7:5). Tumors arose in the deep (6) and superficial (3) soft tissue of the groin (4), buttock (3), thigh (2), flank (1), shoulder (1), and paratesticular tissue (1) and ranged in size from 2 to 20 cm (mean 7.5 cm). Clinical follow-up in 11 patients (9 mo to 20 y; mean 68 mo) showed no recurrences or metastases. As defined above, "fibrosarcoma-like lipomatous neoplasm" is a unique lipomatous tumor that should be distinguished from WDL/(low-grade) dedifferentiated LPS and myxoid LPS on combined histologic/molecular features because of its better prognosis.

Robust Estimators of the Generalized Log-Gamma Distribution

We propose robust estimators of the generalized log-gamma distribution and, more generally, of location-shape-scale families of distributions. A (weighted) Q tau estimator minimizes a tau scale of the differences between empirical and theoretical quantiles. It is n(1/2) consistent; unfortunately, it is not asymptotically normal and, therefore, inconvenient for inference. However, it is a convenient starting point for a one-step weighted likelihood estimator, where the weights are based on a disparity measure between the model density and a kernel density estimate. The one-step weighted likelihood estimator is asymptotically normal and fully efficient under the model. It is also highly robust under outlier contamination. Supplementary materials are available online.

A frequency distribution approach to hotspot identification

We present a new global method for the identification of hotspots in conservation and ecology. The method is based on the identification of spatial structure properties through cumulative relative frequency distributions curves, and is tested with two case studies, the identification of fish density hotspots and terrestrial vertebrate species diversity hotspots. Results from the frequency distribution method are compared with those from standard techniques among local, partially local and global methods. Our approach offers the main advantage to be independent from the selection of any threshold, neighborhood, or other parameter that affect most of the currently available methods for hotspot analysis. The two case studies show how such elements of arbitrariness of the traditional methods influence both size and location of the identified hotspots, and how this new global method can be used for a more objective selection of hotspots.

Measuring the relative effect of factors affecting species distribution model predictions

1. Species distribution models are increasingly used to address conservation questions, so their predictive capacity requires careful evaluation. Previous studies have shown how individual factors used in model construction can affect prediction. Although some factors probably have negligible effects compared to others, their relative effects are largely unknown. 2. We introduce a general "virtual ecologist" framework to study the relative importance of factors involved in the construction of species distribution models. 3. We illustrate the framework by examining the relative importance of five key factors-a missing covariate, spatial autocorrelation due to a dispersal process in presences/absences, sample size, sampling design and modeling technique-in a real study framework based on plants in a mountain landscape at regional scale, and show that, for the parameter values considered here, most of the variation in prediction accuracy is due to sample size and modeling technique. Contrary to repeatedly reported concerns, spatial autocorrelation has only comparatively small effects. 4. This study shows the importance of using a nested statistical framework to evaluate the relative effects of factors that may affect species distribution models.

Different concentrations of Mg++ ions affect nuclear matrix protein distribution during thermal stabilization of isolated nuclei.

The nuclear matrix, a proteinaceous network believed to be a scaffolding structure determining higher-order organization of chromatin, is usually prepared from intact nuclei by a series of extraction steps. In most cell types investigated the nuclear matrix does not spontaneously resist these treatments but must be stabilized before the application of extracting agents. Incubation of isolated nuclei at 37C or 42C in buffers containing Mg++ has been widely employed as stabilizing agent. We have previously demonstrated that heat treatment induces changes in the distribution of three nuclear scaffold proteins in nuclei prepared in the absence of Mg++ ions. We studied whether different concentrations of Mg++ (2.0-5 mM) affect the spatial distribution of nuclear matrix proteins in nuclei isolated from K562 erythroleukemia cells and stabilized by heat at either 37C or 42C. Five proteins were studied, two of which were RNA metabolism-related proteins (a 105-kD component of splicing complexes and an RNP component), one a 126-kD constituent of a class of nuclear bodies, and two were components of the inner matrix network. The localization of proteins was determined by immunofluorescent staining and confocal scanning laser microscope. Mg++ induced significant changes of antigen distribution even at the lowest concentration employed, and these modifications were enhanced in parallel with increase in the concentration of the divalent cation. The different sensitivity to heat stabilization and Mg++ of these nuclear proteins might reflect a different degree of association with the nuclear scaffold and can be closely related to their functional or structural role.

Spatial distribution of motor units recruited during electrical stimulation of the quadriceps muscle versus the femoral nerve.

INTRODUCTION: In this study we investigated differences in the spatial recruitment of motor units (MUs) in the quadriceps when electrical stimulation is applied over the quadriceps belly versus the femoral nerve. METHODS: M-waves and mechanical twitches were evoked using over-the-quadriceps and femoral nerve stimulation of gradually increasing intensity from 22 young, healthy subjects. Spatial recruitment was investigated using recruitment curves of M-waves recorded from the vastus medialis (VM) and vastus lateralis (VL) and of twitches recorded from the quadriceps. RESULTS: At maximal stimulation intensity (Imax), no differences were found between nerve and over-the-quadriceps stimulation. At submaximal intensities, VL M-wave amplitude was higher for over-the-quadriceps stimulation at 40% Imax, and peak twitch force was greater for nerve stimulation at 60% and 80% Imax. CONCLUSIONS: For the VM, MU spatial recruitment during nerve and over-the-quadriceps stimulation of increasing intensity occurred in a similar manner, whereas significant differences were observed for the VL. Muscle Nerve, 2013.

Le pergélisol dans les terrains sédimentaires à forte déclivité: distribution, régime thermique et instabilités

Résumé: Dans le contexte d'un climat de plus en plus chaud, la localisation du pergélisol dans les terrains sédimentaires à forte déclivité et l'évaluation des mouvements de terrain qui y ont cours s'avèrent primordiales. S'insérant dans cette problématique, ce travail de thèse s'articule autour de deux axes de recherche différents. D'un point de vue statique, cette recherche propose une étude de la distribution et des caractéristiques du pergélisol dans les éboulis de la zone périglaciaire alpine. D'un point de vue dynamique, une analyse de l'influence des caractéristiques du pergélisol (teneur en glace, température du pergélisol, etc.) et des variations des températures de l'air et du sol sur les vitesses de fluage des corps sédimentaires gelés est effectuée. Afin de répondre à ce double objectif, l'approche "terrain" a été privilégiée. Pour déterminer la répartition et les caractéristiques du pergélisol, les méthodes traditionnelles de prospection du pergélisol ont été utilisées, à savoir la mesure de la température du sol à la base du manteau neigeux (BTS), la mesure de la température du sol en continu ainsi que la méthode géoélectrique. Les mouvements de terrain ont pour leur part été mesurés à l'aide d'un GPS différentiel. L'étude de la distribution du pergélisol a été effectuée dans une quinzaine d'éboulis situés dans les régions du Mont Gelé (Verbier-Nendaz) et d'Arolla principalement. Dans la plupart des cas, un pergélisol a pu être mis en évidence dans la partie inférieure des accumulations sédimentaires, alors que la partie médiane des éboulis n'est, le plus souvent, pas gelée. Si cette absence de pergélisol se prolonge parfois dans les portions sommitales des pentes, les mesures réalisées montrent que dans d'autres cas des sédiments gelés y sont à nouveau présents. Les résistivités électriques mesurées dans les portions gelées des éboulis étudiés sont dans la plupart des cas nettement inférieures à celles mesurées sur les glaciers rocheux. Des études préalables ont montré que des circulations d'air internes sont responsables de l'anomalie thermique négative et, lorsqu'il existe, du pergélisol que l'on trouve dans la partie inférieure d'éboulis situés plus de 1000 m plus bas que la limite inférieure régionale du pergélisol discontinu. L'étude de quatre sites de basse altitude (1400-1900 m), et notamment l'équipement du site de Dreveneuse (Préalpes Valaisannes) avec deux forages, des capteurs de température de surface et un anémomètre a permis de vérifier et de préciser le mécanisme de ventilation actif au sein des éboulis froids de basse altitude. Ce mécanisme fonctionne de la manière suivante: en hiver, l'air contenu dans l'éboulis, plus chaud et plus léger que l'air extérieur, monte à l'intérieur de l'accumulation sédimentaire et est expulsé dans ses parties sommitales. Cet effet de cheminée provoque une aspiration d'air froid à l'intérieur de la partie inférieure de l'éboulis, causant ainsi un sur-refroidissement marqué du terrain. En été, le mécanisme s'inverse, l'éboulis étant plus froid que l'air environnant. De l'air froid est alors expulsé au bas de la pente. Une ventilation ascendante hivernale a pu être mise en évidence dans certains des éboulis de haute altitude étudiés. Elle est probablement en grande partie responsable de la configuration particulière des zones gelées observées. Même si l'existence d'un effet de cheminée n'a pu être démontrée dans tous les cas, du fait notamment de la glace interstitielle qui entrave le cheminement de l'air, des indices laissant présager son possible fonctionnement existent dans la quasi totalité des éboulis étudiés. L'absence de pergélisol à des altitudes qui lui sont favorables pourrait en tous les cas s'expliquer par un réchauffement du terrain lié à des expulsions d'air relativement chaud. L'étude des mouvements de terrain a été effectuée sur une dizaine de sites, principalement sur des glaciers rocheux, mais également sur une moraine de poussée et - II - Résumé ? abstract quelques éboulis. Plusieurs glaciers rocheux présentent des formes de déstabilisation récente (niches d'arrachement, blocs basculés, apparition de la matrice fine à la surface, etc.), ce qui témoigne d'une récente accélération des vitesses de déplacement. Ce phénomène, qui semble général à l'échelle alpine, est probablement à mettre sur le compte du réchauffement du pergélisol depuis une vingtaine d'années. Les vitesses mesurées sur ces formations sont souvent plus élevées que les valeurs habituellement proposées dans la littérature. On note par ailleurs une forte variabilité inter-annuelle des vitesses, qui semblent dépendre de la variation de la température moyenne annuelle de surface. Abstract: In the context of a warmer climate, the localisation of permafrost in steep sedimentary terrain and the measurement of terrain movements that occur in these areas is of great importance. With respect to these problems, this PhD thesis follows two different research axes. From a static point of view, the research presents a study of the permafrost distribution and characteristics in the talus slopes of the alpine periglacial belt. From a dynamic point of view, an analysis of the influence of the permafrost characteristics (ice content, permafrost temperature, etc.) and air and soil temperature variations on the creep velocities of frozen sedimentary bodies is carried out. In order to attain this double objective, the "field" approach was favoured. To determine the distribution and the characteristics of permafrost, the traditional methods of permafrost prospecting were used, i.e. ground surface temperature measurements at the base of the snow cover (BTS), year-round ground temperature measurements and DC-resistivity prospecting. The terrain movements were measured using a differential GPS. The permafrost distribution study was carried out on 15 talus slopes located mainly in the Mont Gelé (Verbier-Nendaz) and Arolla areas (Swiss Alps). In most cases, permafrost was found in the lower part of the talus slope, whereas the medium part was free of ice. In some cases, the upper part of the talus is also free of permafrost, whereas in other cases permafrost is present. Electrical resistivities measured in the frozen parts of the studied talus are in most cases clearly lower than those measured on rock glaciers. Former studies have shown that internal air circulation is responsible for the negative thermal anomaly and, when it exists, the permafrost present in the lower part of talus slopes located more than 1000 m below the regional lower limit of discontinuous permafrost. The study of four low-altitude talus slopes (1400-1900 m), and notably the equipment of Dreveneuse field site (Valais Prealps) with two boreholes, surface temperature sensors and an anemometer permitted to verify and to detail the ventilation mechanism active in low altitude talus slopes. This mechanism works in the following way: in winter, the air contained in the block accumulation is warmer and lighter than the surrounding air and therefore moves upward in the talus and is expelled in its upper part. This chimney effect induces an aspiration of cold air in the interior of the lower part of talus, that causes a strong overcooling of the ground. In summer, the mechanism is reversed because the talus slope is colder than the surrounding air. Cold air is then expelled in the lower part of the slope. Evidence of ascending ventilation in wintertime could also be found in some of the studied high-altitude talus slopes. It is probably mainly responsible for the particular configuration of the observed frozen areas. Even if the existence of a chimney effect could not be demonstrated in all cases, notably because of interstitial ice that obstructs Résumé ? abstract - III - the air circulation, indices of its presence exist in nearly all the studied talus. The absence of permafrost at altitudes favourable to its presence could be explained, for example, by the terrain warming caused by expulsion of relatively warm air. Terrain movements were measured at about ten sites, mainly on rock glaciers, but also on a push moraine and some talus slopes. Field observations reveal that many rock glaciers display recent destabilization features (landslide scars, tilted blocks, presence of fine grained sediments at the surface, etc.) that indicate a probable recent acceleration of the creep velocities. This phenomenon, which seems to be widespread at the alpine scale, is probably linked to the permafrost warming during the last decades. The measured velocities are often higher than values usually proposed in the literature. In addition, strong inter-annual variations of the velocities were observed, which seems to depend on the mean annual ground temperature variations.