Evaluation of a mixed approach combining stationary and wearable systems to monitor gait over long distance.

Thanks to decades of research, gait analysis has become an efficient tool. However, mainly due to the price of the motion capture systems, standard gait laboratories have the capability to measure only a few consecutive steps of ground walking. Recently, wearable systems were proposed to measure human motion without volume limitation. Although accurate, these systems are incompatible with most of existing calibration procedures and several years of research will be necessary for their validation. A new approach consisting of using a stationary system with a small capture volume for the calibration procedure and then to measure gait using a wearable system could be very advantageous. It could benefit from the knowledge related to stationary systems, allow long distance monitoring and provide new descriptive parameters. The aim of this study was to demonstrate the potential of this approach. Thus, a combined system was proposed to measure the 3D lower body joints angles and segmental angular velocities. It was then assessed in terms of reliability towards the calibration procedure, repeatability and concurrent validity. The dispersion of the joint angles across calibrations was comparable to those of stationary systems and good reliability was obtained for the angular velocities. The repeatability results confirmed that mean cycle kinematics of long distance walks could be used for subjects' comparison and pointed out an interest for the variability between cycles. Finally, kinematics differences were observed between participants with different ankle conditions. In conclusion, this study demonstrated the potential of a mixed approach for human movement analysis.

Inference about the number of contributors to a DNA mixture: Comparative analyses of a Bayesian network approach and the maximum allele count method.

In the forensic examination of DNA mixtures, the question of how to set the total number of contributors (N) presents a topic of ongoing interest. Part of the discussion gravitates around issues of bias, in particular when assessments of the number of contributors are not made prior to considering the genotypic configuration of potential donors. Further complication may stem from the observation that, in some cases, there may be numbers of contributors that are incompatible with the set of alleles seen in the profile of a mixed crime stain, given the genotype of a potential contributor. In such situations, procedures that take a single and fixed number contributors as their output can lead to inferential impasses. Assessing the number of contributors within a probabilistic framework can help avoiding such complication. Using elements of decision theory, this paper analyses two strategies for inference on the number of contributors. One procedure is deterministic and focuses on the minimum number of contributors required to 'explain' an observed set of alleles. The other procedure is probabilistic using Bayes' theorem and provides a probability distribution for a set of numbers of contributors, based on the set of observed alleles as well as their respective rates of occurrence. The discussion concentrates on mixed stains of varying quality (i.e., different numbers of loci for which genotyping information is available). A so-called qualitative interpretation is pursued since quantitative information such as peak area and height data are not taken into account. The competing procedures are compared using a standard scoring rule that penalizes the degree of divergence between a given agreed value for N, that is the number of contributors, and the actual value taken by N. Using only modest assumptions and a discussion with reference to a casework example, this paper reports on analyses using simulation techniques and graphical models (i.e., Bayesian networks) to point out that setting the number of contributors to a mixed crime stain in probabilistic terms is, for the conditions assumed in this study, preferable to a decision policy that uses categoric assumptions about N.

Measuring and Explaining the Increase of Travel Distance : A Multilevel Analysis Using Repeated Cross Sectional Travel Surveys

We want to shed some light on the development of person mobility by analysing the repeated cross-sectional data of the four National Travel Surveys (NTS) that were conducted in Germany since the mid seventies. The above mentioned driving forces operate on different levels of the system that generates the spatial behaviour we observe: Travel demand is derived from the needs and desires of individuals to participate in spatially separated activities. Individuals organise their lives in an interactive process within the context they live in, using given infrastructure. Essential determinants of their demand are the individual's socio-demographic characteristics, but also the opportunities and constraints defined by the household and the environment are relevant for the behaviour which ultimately can be realised. In order to fully capture the context which determines individual behaviour, the (nested) hierarchy of persons within households within spatial settings has to be considered. The data we will use for our analysis contains information on these three levels. With the analysis of this micro-data we attempt to improve our understanding of the afore subsumed macro developments. In addition we will investigate the prediction power of a few classic sociodemographic variables for the daily travel distance of individuals in the four NTS data sets, with a focus on the evolution of this predictive power. The additional task to correctly measure distances travelled by means of the NTS is threatened by the fact that although these surveys measure the same variables, different sampling designs and data collection procedures were used. So the aim of the analysis is also to detect variables whose control corrects for the known measurement error, as a prerequisite to apply appropriate models in order to better understand the development of individual travel behaviour in a multilevel context. This task is complicated by the fact that variables that inform on survey procedures and outcomes are only provided with the data set for 2002 (see Infas and DIW Berlin, 2003).

Truncated robust distance for clinical laboratory safety data monitoring and assessment.

Laboratory safety data are routinely collected in clinical studies for safety monitoring and assessment. We have developed a truncated robust multivariate outlier detection method for identifying subjects with clinically relevant abnormal laboratory measurements. The proposed method can be applied to historical clinical data to establish a multivariate decision boundary that can then be used for future clinical trial laboratory safety data monitoring and assessment. Simulations demonstrate that the proposed method has the ability to detect relevant outliers while automatically excluding irrelevant outliers. Two examples from actual clinical studies are used to illustrate the use of this method for identifying clinically relevant outliers.

Sex differences in the computation of traveling distance

Path integration is known to provide information to keep track of spatial location. Surprisingly, few investigations concerning sex differences in computation of the traveling distance have been done. This work was aimed at analyzing the reproduction of both passive and active linear displacements in women and men. To this end, the displacement of blindfolded subjects was done in a wheelchair, then on foot, three times in each condition for a fixed distance. Copies of passive and active traveling distance, distance estimations and pointing responses towards the starting point were analyzed. In passive condition and comparatively to men, women error was larger. Whereas traveling distance was generally underestimated in women, it was overestimated in men. In active condition, no sex differences were observed. When blindfolded subjects have to estimate the traveling distance, the female error was larger than the male one. But, when subjects were asked to indicate the visual cue corresponding to the traveling distance, the male error was larger than the female one. Finally, pointing to the starting point (0°) after a whole-body rotation showed a larger deviation from 0° in men than in women. These results suggest that sex of the subjects influence brain computation of path integration information.

Adaptively weighted maximum likelihood estimation of discrete distributions

SummaryDiscrete data arise in various research fields, typically when the observations are count data.I propose a robust and efficient parametric procedure for estimation of discrete distributions. The estimation is done in two phases. First, a very robust, but possibly inefficient, estimate of the model parameters is computed and used to indentify outliers. Then the outliers are either removed from the sample or given low weights, and a weighted maximum likelihood estimate (WML) is computed.The weights are determined via an adaptive process such that if the data follow the model, then asymptotically no observation is downweighted.I prove that the final estimator inherits the breakdown point of the initial one, and that its influence function at the model is the same as the influence function of the maximum likelihood estimator, which strongly suggests that it is asymptotically fully efficient.The initial estimator is a minimum disparity estimator (MDE). MDEs can be shown to have full asymptotic efficiency, and some MDEs have very high breakdown points and very low bias under contamination. Several initial estimators are considered, and the performances of the WMLs based on each of them are studied.It results that in a great variety of situations the WML substantially improves the initial estimator, both in terms of finite sample mean square error and in terms of bias under contamination. Besides, the performances of the WML are rather stable under a change of the MDE even if the MDEs have very different behaviors.Two examples of application of the WML to real data are considered. In both of them, the necessity for a robust estimator is clear: the maximum likelihood estimator is badly corrupted by the presence of a few outliers.This procedure is particularly natural in the discrete distribution setting, but could be extended to the continuous case, for which a possible procedure is sketched.RésuméLes données discrètes sont présentes dans différents domaines de recherche, en particulier lorsque les observations sont des comptages.Je propose une méthode paramétrique robuste et efficace pour l'estimation de distributions discrètes. L'estimation est faite en deux phases. Tout d'abord, un estimateur très robuste des paramètres du modèle est calculé, et utilisé pour la détection des données aberrantes (outliers). Cet estimateur n'est pas nécessairement efficace. Ensuite, soit les outliers sont retirés de l'échantillon, soit des faibles poids leur sont attribués, et un estimateur du maximum de vraisemblance pondéré (WML) est calculé.Les poids sont déterminés via un processus adaptif, tel qu'asymptotiquement, si les données suivent le modèle, aucune observation n'est dépondérée.Je prouve que le point de rupture de l'estimateur final est au moins aussi élevé que celui de l'estimateur initial, et que sa fonction d'influence au modèle est la même que celle du maximum de vraisemblance, ce qui suggère que cet estimateur est pleinement efficace asymptotiquement.L'estimateur initial est un estimateur de disparité minimale (MDE). Les MDE sont asymptotiquement pleinement efficaces, et certains d'entre eux ont un point de rupture très élevé et un très faible biais sous contamination. J'étudie les performances du WML basé sur différents MDEs.Le résultat est que dans une grande variété de situations le WML améliore largement les performances de l'estimateur initial, autant en terme du carré moyen de l'erreur que du biais sous contamination. De plus, les performances du WML restent assez stables lorsqu'on change l'estimateur initial, même si les différents MDEs ont des comportements très différents.Je considère deux exemples d'application du WML à des données réelles, où la nécessité d'un estimateur robuste est manifeste : l'estimateur du maximum de vraisemblance est fortement corrompu par la présence de quelques outliers.La méthode proposée est particulièrement naturelle dans le cadre des distributions discrètes, mais pourrait être étendue au cas continu.

Optimal robust estimates using the Hellinger distance

Optimal robust M-estimates of a multidimensional parameter are described using Hampel's infinitesimal approach. The optimal estimates are derived by minimizing a measure of efficiency under the model, subject to a bounded measure of infinitesimal robustness. To this purpose we define measures of efficiency and infinitesimal sensitivity based on the Hellinger distance.We show that these two measures coincide with similar ones defined by Yohai using the Kullback-Leibler divergence, and therefore the corresponding optimal estimates coincide too.We also give an example where we fit a negative binomial distribution to a real dataset of "days of stay in hospital" using the optimal robust estimates.

Phenotypic novelty in experimental hybrids is predicted by the genetic distance between species of cichlid fish.

BACKGROUND: Transgressive segregation describes the occurrence of novel phenotypes in hybrids with extreme trait values not observed in either parental species. A previously experimentally untested prediction is that the amount of transgression increases with the genetic distance between hybridizing species. This follows from QTL studies suggesting that transgression is most commonly due to complementary gene action or epistasis, which become more frequent at larger genetic distances. This is because the number of QTLs fixed for alleles with opposing signs in different species should increase with time since speciation provided that speciation is not driven by disruptive selection. We measured the amount of transgression occurring in hybrids of cichlid fish bred from species pairs with gradually increasing genetic distances and varying phenotypic similarity. Transgression in multi-trait shape phenotypes was quantified using landmark-based geometric morphometric methods. RESULTS: We found that genetic distance explained 52% and 78% of the variation in transgression frequency in F1 and F2 hybrids, respectively. Confirming theoretical predictions, transgression when measured in F2 hybrids, increased linearly with genetic distance between hybridizing species. Phenotypic similarity of species on the other hand was not related to the amount of transgression. CONCLUSION: The commonness and ease with which novel phenotypes are produced in cichlid hybrids between unrelated species has important implications for the interaction of hybridization with adaptation and speciation. Hybridization may generate new genotypes with adaptive potential that did not reside as standing genetic variation in either parental population, potentially enhancing a population's responsiveness to selection. Our results make it conceivable that hybridization contributed to the rapid rates of phenotypic evolution in the large and rapid adaptive radiations of haplochromine cichlids.

Proximal tibia volumetric bone mineral density is correlated to the magnitude of local acceleration in male long distance runners.

Introduction: The beneficial effect of physical exercise on bone mineral density (BMD) is at least partly explained by the forces exerted directly on the bones. Male runners present generally higher BMD than sedentary individuals. We postulated that the proximal tibia BMD is related to the running distance as well as to the magnitude of the shocks (while running) in male runners. Methods: A prospective study (three yearly measurements) included 81 healthy male subjects: 16 sedentary lean subjects and three groups of runners (5-30 km/week, n=19; 30-50 km/week, n=29; 50-100 km/week, n=17). Several measurements were performed at the proximal tibia level: volumetric BMD (vBMD), cortical index (CI) i.e. an index of cortical bone thickness and peak accelerations (an index of shocks during heel strike) while running (measured by a 3-D accelerometer). A general linear model assessed the prediction of vBMD or CI by a) simple effects (running distance, peak accelerations, time) and b) interactions (for instance if vBMD prediction by peak acceleration depends on running distance). Results: CI and vBMD a) increase with running distance to reach a plateau over 30 km/wk, b) are positively associated with peak accelerations over 30 km/week. Discussion: Running may be associated with high peak accelerations in order to have beneficial effects on BMD. More important strains are needed to be associated with the same increase in BMD during running sessions of short duration than those of long duration. Conclusion: CI and vBMD are associated with the magnitude of the shocks during heel strike in runners. Key words: Bone mineral density, strains, physical exercise, running distance.

Long-distance wound signalling in arabidopsis

The leaves of all plants use elaborate and inducible defence systems to protect themselves. A wide variety of such defences are known and they include defence chemicals such as alkaloids, phenolics and terpenes, physical structures ranging from fibre cells to silica deposits, and a wide variety of defence proteins many of which target digestive processes in herbivores. It has long been known that the defence responses of plants under attack by insects are not restricted to the site of attack. Instead, if a leaf is damaged, defence can be triggered in other parts of the plant body, for example in distal leaves or even in roots and flowers. This raises the question of what are the organ-to-organ signals that coordinate this process. Several hypotheses have been proposed. These include the long-distance transfer of chemical signals through the plant vasculature, hydraulic signals that may transit through the xylem, and electrical signals that would move through living tissues such as the phloem. Much evidence for each of these scenarios has been published. In this thesis we took advantage of the fact that many plant defence responses are regulated by a signal transduction pathway based on a molecule called jasmonic acid. We used this molecule, one of its derivatives (jasmonoyl-isoleucine), and some of the genes it regulates as markers. Using these we investigated the possible role of the electrical signals in the leaf- to-leaf activation of the jasmonate pathway. We found that feeding insects stimulate easily detected electrical activity in the leaves of Arabidopsis thaliana and we used non-invasive surface electrodes to record this activity. This approach showed that jasmonate pathway activity and the electrical activity provoked by mechanical wounding occurred within identical spatial boundaries. Measurements of the apparent speed of surface potentials agreed well with previous velocity estimates for the speed of leaf-to-leaf signals that activate the jasmonate pathway. Using this knowledge we were able to investigate the effects of current injection into Arabidopsis leaves. This resulted in the strong expression of many jasmonate-regulated genes. All these results showed that electrical activity and the activation of jasmonate signalling were highly correlated. In order to test for possible causal links between the two processes, we conducted a small-scale reverse genetic screen on a series of T-DNA insertion mutants in ion channel genes and in other genes encoding proteins such as proton pumps. This screen, which was based on surface potential measurements, revealed that mutations in genes related to ionotropic glutamate receptors in animals had impaired electrical activity after wounding. Combining mutation of two of these glutamate-receptor-like genes in a double mutant reduced the response of leaves to current injection. When a leaf of this double mutant was wounded it failed to transmit a long-distance signal to a distal leaf. This result distinguished the double mutant from the wild-type plant and provides the first genetic evidence that electrical signalling is necessary to coordinate defence responses between organs in plants. - Les feuilles des plantes disposent de systèmes de défense inductibles très élaborés. Un grand nombre de ces systèmes de défenses sont connus et sont basés sur des composés chimiques comme les alcaloïdes, les composés phénoliques ou les terpènes, des systèmes physiques allant de la production de cellules fibreuses aux cristaux de silice ainsi qu'un grand nombre de protéines de défense ciblant le processus digestif des herbivores. Il est connu dépuis longtemps que la réponse défensive de la plante face à l'attaque pas un insecte n'est pas seulement localisée au niveau de la zone d'attaque. A la place, si une feuille est attaquée, les systèmes de défense peuvent être activés ailleurs dans la plante, comme par exemple dans d'autres feuilles, les racines ou même les fleurs. Ces observations soulèvent la question de la nature des signaux d'organes à organes qui régulent ces systèmes. Plusieurs hypothèses ont été formulées; une ou plusieures molécules pourraient être véhiculées dans la plante grâce au système vasculaire, un signal hydraulique transmis au travers du xylème ou encore des signaux électriques transmis par les cellules comme dans le phloème par exemple. De nombreuses études ont été publiées sur ces différentes hypothèses. Dans ce travail de thèse, nous avons choisi d'utiliser à notre avantage le fait que de nombreuses réponses de défense de la plante sont régulées par une même voie de signalisation utilisant l'acide jasmonique. Nous avons utilisé comme marqueurs cette molécule, un de ses dérivés (le jasmonoyl-isoleucine) ainsi que certains des gènes que l'acide jasmonique régule. Nous avons alors testé l'implication de la transmission de signaux électriques dans l'activation de la voie du jasmonate de feuille à feuille. Nous avons découvert que les insectes qui se nourrissent de feuilles d'Arabidopsis thaliana activent un signal électrique que nous avons pu mesurer grâce à une technique non invasive d'électrodes de surface. Les enregistrements ont montré que la génération de signaux électriques et l'activation de la voie du jasmonate avaient lieu aux mêmes endroits. La mesure de la vitesse de déplacement des impulsions électriques correspond aux estimations faites concernant l'activation de la voie du jasmonate. Grâce à cela, nous avons pu tester l'effet d'injection de courant électrique dans les feuilles d'Arabidopsis. La conséquence a été une forte expression de nombreux gènes de la voie du jasmonate, suggérant une forte corrélation entre l'activité électrique et l'activation de la voie du jasmonate. Afin de tester le lien de cause entre ces deux phénomènes, nous avons entrepris un criblage génétique sur une série de mutants d'insertion à l'ADN-T dans des gènes de canaux ioniques et d'autres gènes d'intérêt comme les gènes des pompes à protons. Ce criblage, basé sur la mesure de potentiels de surface, a permis de montrer que plusieurs mutations de gènes liés aux récepteurs au glutamate ionotropique présentent une baisse drastique de leurs activités électriques après une blessure mécanique des feuilles par rapport au type sauvage. Par la combinaison de deux mutations de ces récepteurs au glutamate en un double mutant, on obtient une réponse à la stimulation électrique encore plus faible. Quand une feuille du double mutant est blessée, elle est incapable de transmettre un signal à longue distance vers une feuille éloignée. Ce résultat permet de distinguer le double mutant de la plante sauvage et amène la première preuve génétique que l'activité électrique est nécessaire pour coordonner les réponses de défense entre les organes chez les plantes.

Kinetics of crystal growth and equilibrium domains in eclogite of the Sesia Zone, Western Alps

Kinetics of crystal growth and equilibrium domains in eclogite of the Sesia Zone, Western Alps Darbellay Bastien Institut de Minéralogie et Géochimie Résumé grand public Comme toute matière, la roche est sensible à son environnement et cherche à s'adapter pour acquérir un état stable (état d'équilibre). Les changements des conditions physiques (température et pression) vont ainsi impliquer des modifications dans la roche. Le métamorphisme est l'étude de ces changements. Les minéraux qui constituent la roche peuvent modifier, leur structure, leur chimie ou être remplacer par d'autres minéraux plus stables. Il est ainsi crucial de déterminer les processus responsables et limitant de la croissance minérale. Trois processus permettent la croissance ; (1) la dissolution des éléments du réactant, (2) le transport de ces éléments vers le site de croissance, (3) l'incorporation de ces éléments dans la nouvelle structure. Cette thèse se focalise sur les structures des minéraux de haute pression (forme, zonation chimique, structure interne) pour essayer de déterminer les facteurs importants à l'origine de leur état final. Les zones d'étude se situent dans la zone de Sésia. La première partie traite de la problématique liée à l'incorporation d'un élément dans une structure minérale. A l'image de la croissance humaine, les irrégularités minéralogiques permettent de mettre en lumière un dysfonctionnement de la croissance due à un excès ou à une carence d'un élément. Bien dosé, cet élément est cependant essentiel à la croissance. Les zoisites (épidotes) des métabasites de la région de Cima di Bonze montrent une zonation chimique en sablier. Dans cette zonation la teneur en fer excède la capacité maximum que peut contenir la structure orthorhombique de la zoisite. Des défauts de structure permettent l'accommodation de cet excès. La zoisite peut ainsi adapter sa structure pour permettre l'incorporation d'une relativement grande quantité de fer. Les études précédentes montraient, pour des conditions similaires, la formation de deux épidotes distinctes. La deuxième partie se penche sur la compétition entre le minéral qui fait sa croissance et les minéraux (réactants) qui l'entourent. Les métapélites de la région du Monte Mucrone contiennent des grenats atollaires. Des études détaillées de la texture et de la zonation chimique du grenat ainsi qu'une modélisation thermodynamique ont permis de mieux cerner les facteurs importants responsables de la forme atollaire. Cette structure est obtenue par un changement du comportement de la croissance du grenat le long d'un chemin P-T hercynien. Dans un premier stade, le grenat croît rapidement et consume peu le quartz de la matrice. La croissance se fait ainsi le long des jointures des grains de quartz. Dans un second temps, les changements de conditions PT donnent une croissance lente du grenat et une forte consommation du quartz. Le grenat peut ainsi développer sa forme dodécaédrale classique. La troisième partie s'intéresse aux distances de transport par diffusion d'un élément (ici l'argon) durant la haute pression. Pour ce faire, un profile d'âges 40Ar/39Ar sur biotite a été mesuré depuis un veine de haute pression riche en argon jusque dans son encaissant (granitoïd du Monte Mucrone). Le profile montre une répartition des âges suivant une courbe de diffusion. Le transport se fait sur une longueur de deux centimètre avec l'aide d'un fluide. Il est réduit à une échelle millimétrique quand la phase fluide disparaît. Cette étude montre ainsi les difficultés de transport des éléments durant la haute pression ne permettant pas un rééquilibrage de la roche à grande échelle. Kinetics of crystal growth and equilibrium domains in eclogite of the Sesia Zone, Western Alps Darbellay Bastien Institut de Minéralogie et Géochimie Résumé de thèse Les processus de croissance (diffusion des éléments et les réactions d'interface) et les conditions dans lesquelles les minéraux grandissent (température, pression, fluide, composition chimique de la roche), déterminent la texture ainsi que la zonation des minéraux. Cette thèse se focalise, par le biais de textures peu communes, sur trois différents processus impliqués dans la croissance minérale à haute pression (Zone de Sésia, Alpes de l'Ouest, Italie). L'incorporation d'un élément dans une structure minérale ne peut se faire que dans des sites en accord avec la taille et la charge ionique de l'élément. De plus, la balance de charge doit être maintenue dans le minéral. La régularité de la structure cristalline fixe ainsi une limite maximum de concentration d'un élément donné. Les zoisites provenant des métabasites de la région de Cima di Bonze montrent des zonations en sablier caractérisées par une concentration anormale en fer. La zonation se marque par une différente teinte de biréfringence et par un plus grand angle d'extinction que le reste de la zoisite. Une inter-croissance de clinozoisite à l'intérieur de la structure orthorhombique de la zoisite peut ainsi être suspectée. Les analyses XRD (diffraction des rayons x) ainsi que les analyses Raman ne confirment pas cette suspicion. Seules les analyses TEM (microscope à électrons transmis) montrent des défauts de structure pouvant être interprétés comme des modules de clinozoisite. Ils ne peuvent cependant pas être considérés comme une phase thermodynamique. Un nouveau trou d'immiscibilité entre deux zoisite (X ep= 0.1 and Xep = 0.15) a ainsi pu être établi. Dans les métapélites la région du Monte Mucrone, des grenats fortement zonés montrent une évolution texturale singulière. Ils présentent une forme initiale de `champignon' qui se développe pour former une structure atollaire finale. L'étude conjuguée de la structure 3D et des zonations, ainsi que l'établissement d'un model thermodynamique, indiquent que ces structures proviennent de deux épisodes de croissances : (1) La croissance du grenat durant un chemin prograde hercynien (de 525 °C et 6.2 kbar à 640 °C et 9 kbar) permet la formation des textures atollaires. Elles sont le résultat d'une croissance poecilitique initiale suivie d'une croissance idiomorphique du grenat. (2) La structure est rendue plus complexe par la cristallisation d'un grenat homogène tout autour ainsi qu'à l'intérieur du grenat hercynien durant la haute pression alpine (550 °C and 20 kbar). L'arrivée de l'eau durant la haute pression facilite le transport d'éléments et permet une cristallisation rapide du grenat. La diffusion peut être un facteur limitant de la croissance minéralogique. Elle a aussi une grande importance pour la géochronologie. Une veine de haute pression à l'intérieur du granitoïde du Monte Mucrone a été étudiée dans le but de déterminer la distance de diffusion de l'argon. Le profile d'âges 40Ar/39Ar sur biotites, établi de la veine vers le métagranitoïde, suit une courbe de diffusion. Les âges sont élevés proche de la veine (800 Ma) puis décroissent jusqu'à des âges homogènes (170-150 Ma) à deux centimètres de la veine. La présence de fluide, marqué par de hautes concentrations en chlore, permet une diffusion centimétrique. Cependant, la distance est réduite à une échelle millimétrique quand le fluide est absent. Les très faibles distances de diffusion préservent les âges pré-alpins et impliquent un événement géologique pour les âges de 170-150 Ma. Kinetics of crystal growth and equilibrium domains in eclogite of the Sesia Zone, Western Alps Darbellay Bastien Institut de Minéralogie et Géochimie Thesis abstract Rock textures and zonings are the consequence of growth processes (element diffusion and interface reaction) steered by the environment in which they grew (pressure, differential stress, temperature, fluid and rock composition). The thesis presented here focuses on three different topics, each of it dealing with aspects of mineral growth processes during subduction, in a high-pressure environment. All studies were conducted in the Sesia Zone of the Western European Alps, Italy. The first study addresses the crystallography and geochemistry of element incorporation in zoisite, one of the major hydrous minerals found in subduction zone rocks. Elements can be incorporated into a mineral structure only on crystallographic sites that offer enough space for the ion and the overall charge balance has to be maintained. Element concentrations are hence limited. Incorporation of some elements produces complex zoning, including hourglass like patterns, which are the focus of the first contribution. Zoisites from Cima di Bonze (Sesia Zone) show spectacular hourglass zoning defined by Fe-content variations. The hourglass zones have a distinct birefringence and a different extinction angle than the regular part of the zoisite. We show by detailed XRD (X-ray diffraction) and confocal Raman analyses that the high Fe-zones are nevertheless zoisite, and not clinozoisite as one might expect. High resolution TEM (transmission electron microscopy) analyses show planar defects on (100) that can be interpreted as small-scale clinozoisite modules. However, these clinozoisites cannot be interpreted as a distinctive thermodynamic phase and the entire mineral has to be considered as zoisite. The miscibility gap between two zoisites (Xep = 0.1 and Xep = 0.15) can be then definite at 550 ± 50°C and 14 to 20 Kbar. Strongly zoned garnets in quartz rich metapelite from the Monte Mucrone area (Sesia Zone) show evolution form 3D mushroom to atoll structure. The second contribution presents textural investigations, garnet zoning and thermodynamic modeling that demonstrate that atoll garnets are the result of two distinctive growth events. (1) Garnet atoll structure is already formed during a prograde Hercynian path from 525 °C and 6.2 kbar to 640 °C and 9 kbar. It results in an initial poikilitic growth followed by a final idiomorphic growth event. (2) Alpine HP garnet are homogenous (550 °C and 20 kbar) and grew around and also inside the Hercynian garnet. Lack of prograde Alpine garnet and fast growth of the HP garnet is explained by the absence of water during much of the prograde path. Water saturation was only observed towards the end, close towards the peak metamorphic conditions. Diffusion could be a limiting factor for crystal growth. It has also a great importance in geochronology. HP vein inside the metagranitoide of the Monte Mucrone (~300 Ma) was investigated to determine argon diffusion scales during high-pressure metamorphism. 40Ar/39Ar biotite ages profile from the vein toward the metagranodiorite show a diffusion curve: old ages (800 Ma) located close to the vein decrease until homogenous 170-150 Ma ages are obtained, two centimeter away from the vein. Centimeter-scale diffusion occurs with help of a fluid phase marked by high chlorine concentrations. Argon diffusion is reduced to a millimeter scale when free fluid is absent. Very short diffusion distance permits to preserve pre-Alpine ages. The 170-150 Ma ages are considered to be geologic meaningful, probably resulting from the extensional tectonics linked to opening of the Tethian ocean.

Oxygen uptake kinetics and middle distance swimming performance.

OBJECTIVE: The aim of this study was to determine whether V˙O(2) kinetics and specifically, the time constant of transitions from rest to heavy (τ(p)H) and severe (τ(p)S) exercise intensities, are related to middle distance swimming performance. DESIGN: Fourteen highly trained male swimmers (mean ± SD: 20.5 ± 3.0 yr; 75.4 ± 12.4 kg; 1.80 ± 0.07 m) performed an discontinuous incremental test, as well as square wave transitions for heavy and severe swimming intensities, to determine V˙O(2) kinetics parameters using two exponential functions. METHODS: All the tests involved front-crawl swimming with breath-by-breath analysis using the Aquatrainer swimming snorkel. Endurance performance was recorded as the time taken to complete a 400 m freestyle swim within an official competition (T400), one month from the date of the other tests. RESULTS: T400 (Mean ± SD) (251.4 ± 12.4 s) was significantly correlated with τ(p)H (15.8 ± 4.8s; r=0.62; p=0.02) and τ(p)S (15.8 ± 4.7s; r=0.61; p=0.02). The best single predictor of 400 m freestyle time, out of the variables that were assessed, was the velocity at V˙O(2max)vV˙O(2max), which accounted for 80% of the variation in performance between swimmers. However, τ(p)H and V˙O(2max) were also found to influence the prediction of T400 when they were included in a regression model that involved respiratory parameters only. CONCLUSIONS: Faster kinetics during the primary phase of the V˙O(2) response is associated with better performance during middle-distance swimming. However, vV˙O(2max) appears to be a better predictor of T400.

Metamorphism and kinetics in the Torres del Paine contact aureole

Contact aureoles provide an excellent geologic environment to study the mechanisms of metamorphic reactions in a natural system. The Torres del Paine (TP) intrusion is one of the most spectacular natural laboratories because of its excellent outcrop conditions. It formed in a period from 12.59 to 12.43 Ma and consists of three large granite and four smaller mafic batches. The oldest granite is on top, the youngest at the bottom of the granitic complex, and the granites overly the mafic laccolith. The TP intruded at a depth of 2-3 km into regional metamorphic anchizone to greenschist facies pelites, sandstones, and conglomerates of the Cerro Toro and Punta Barrosa formations. It formed a thin contact aureole of 150-400 m width. This thesis focuses on the reaction kinetics of the mineral cordierite in the contact aureole using quantitative textural analysis methods. First cordierite was formed from chlorite break¬down (zone I, ca. 480 °C, 750 bar). The second cordierite forming reaction was the muscovite break-down, which is accompanied by a modal decrease in biotite and the appearance of k- feldspar (zone II, 540-550 °C, 750 bar). Crystal sizes of the roundish, poikiloblastic cordierites were determined from microscope thin section images by manually marking each crystal. Images were then automatically processed with Matlab. The correction for the intersection probability of each crystal radius yields the crystal size distribution in the rock. Samples from zone I below the laccolith have the largest crystals (0.09 mm). Cordierites from zone II are smaller, with a maximum crystal radius of 0.057 mm. Rocks from zone II have a larger number of small cordierite crystals than rocks from zone I. A combination of these quantitative analysis with numerical modeling of nucleation and growth, is used to infer nucleation and growth parameters which are responsible for the observed mineral textures. For this, the temperature-time paths of the samples need to be known. The thermal history is complex because the main body of the intrusion was formed by several intrusive batches. The emplacement mechanism and duration of each batch can influence the thermal structure in the aureole. A possible subdivision of batches in smaller increments, so called pulses, will focus heat at the side of the intrusion. Focusing all pulses on one side increases the contact aureole size on that side, but decreases it on the other side. It forms a strongly asymmetric contact aureole. Detailed modeling shows that the relative thicknesses of the TP contact aureole above and below the intrusion (150 and 400 m) are best explained by a rapid emplacement of at least the oldest granite batch. Nevertheless, temperatures are significantly too low in all models, compared to observed mineral assemblages in the hornfelses. Hence, an other important thermal mechanisms needs to take place in the host rock. Clastic minerals in the immature sediments outside the contact aureole are hydrated due to small amounts of expelled fluids during contact metamorphism. This leads to a temperature increase of up to 50 °C. The origin of fluids can be traced by stable isotopes. Whole rock stable isotope data (6D and δ180) and chlorine concentrations in biotite document that the TP intrusion induced only very small amounts of fluid flow. Oxygen whole rock data show δ180 values between 9.0 and 10.0 %o within the first 5 m of the contact. Values increase to 13.0 - 15.0 %o further away from the intrusion. Whole rock 6D values display a more complex zoning. First, host rock values (-90 to -70 %o) smoothly decrease towards the contact by ca. 20 %o, up to a distance of ca. 150 m. This is followed by an increase of ca. 20 %o within the innermost 150 m of the aureole (-97.0 to -78 %o at the contact). The initial decrease in 6D values is interpreted to be due to Rayleigh fractionation accompanying the dehydration reactions forming cordierite, while the final increase reflects infiltration of water-rich fluids from the intrusion. An over-estimate on the quantity and the corresponding thermal effect yields a temperature increase of less than 30 °C. This suggests that fluid flow might have contributed only for a small amount to the thermal evolution of the system. A combination of the numerical growth model with the thermal model, including the hydration reaction enthalpies but neglecting fluid flow and incremental growth, can be used to numerically reproduce the observed cordierite textures in the contact aureole. This yields kinetic parameters which indicate fast cordierite crystallization before the thermal peak in the inner aureole, and continued reaction after the thermal peak in the outermost aureole. Only small temperature dependencies of the kinetic parameters seem to be needed to explain the obtained crystal size data. - Les auréoles de contact offrent un cadre géologique privilégié pour l'étude des mécanismes de réactions métamorphiques associés à la mise en place de magmas dans la croûte terrestre. Par ses conditions d'affleurements excellentes, l'intrusion de Torres del Paine représente un site exceptionnel pour améliorer nos connaissances de ces processus. La formation de cette intrusion composée de trois injections granitiques principales et de quatre injections mafiques de volume inférieur couvre une période allant de 12.50 à 12.43 Ma. Le plus vieux granite forme la partie sommitale de l'intrusion alors que l'injection la plus jeune s'observe à la base du complexe granitique; les granites recouvrent la partie mafique du laccolite. L'intrusion du Torres del Paine s'est mise en place a 2-3 km de profondeur dans un encaissant métamorphique. Cet encaissant est caractérisé par un métamorphisme régional de faciès anchizonal à schiste vert et est composé de pélites, de grès, et des conglomérats des formations du Cerro Toro et Punta Barrosa. La mise en place des différentes injections granitiques a généré une auréole de contact de 150-400 m d'épaisseur autour de l'intrusion. Cette thèse se concentre sur la cinétique de réaction associée à la formation de la cordiérite dans les auréoles de contact en utilisant des méthodes quantitatives d'analyses de texture. On observe plusieurs générations de cordiérite dans l'auréole de contact. La première cordiérite est formée par la décomposition de la chlorite (zone I, environ 480 °C, 750 bar), alors qu'une seconde génération de cordiérite est associée à la décomposition de la muscovite, laquelle est accompagnée par une diminution modale de la teneur en biotite et l'apparition de feldspath potassique (zone II, 540-550 °C, 750 bar). Les tailles des cristaux de cordiérites arrondies et blastic ont été déterminées en utilisant des images digitalisées des lames minces et en marquant individuellement chaque cristal. Les images sont ensuite traitées automatiquement à l'aide du programme Matlab. La correction de la probabilité d'intersection en fonction du rayon des cristaux permet de déterminer la distribution de la taille des cristaux dans la roche. Les échantillons de la zone I, en dessous du lacolite, sont caractérisés par de relativement grands cristaux (0.09 mm). Les cristaux de cordiérite de la zone II sont plus petits, avec un rayon maximal de 0.057 mm. Les roches de la zone II présentent un plus grand nombre de petits cristaux de cordiérite que les roches de la zone I. Une combinaison de ces analyses quantitatives avec un modèle numérique de nucléation et croissance a été utilisée pour déduire les paramètres de nucléation et croissance contrôlant les différentes textures minérales observées. Pour développer le modèle de nucléation et de croissance, il est nécessaire de connaître le chemin température - temps des échantillons. L'histoire thermique est complexe parce que l'intrusion est produite par plusieurs injections successives. En effet, le mécanisme d'emplace¬ment et la durée de chaque injection peuvent influencer la structure thermique dans l'auréole. Une subdivision des injections en plus petits incréments, appelés puises, permet de concentrer la chaleur dans les bords de l'intrusion. Une mise en place préférentielle de ces puises sur un côté de l'intrusion modifie l'apport thermique et influence la taille de l'auréole de contact produite, auréole qui devient asymétrique. Dans le cas de la première injection de granite, une modélisation détaillée montre que l'épaisseur relative de l'auréole de contact de Torres del Paine au-dessus et en dessous de l'intrusion (150 et 400 m) est mieux expliquée par un emplacement rapide du granite. Néanmoins, les températures calculées dans l'auréole de con¬tact sont trop basses pour que les modèles thermiques soient cohérants par rapport à la taille de cette auréole. Ainsi, un autre mecanisme exothermique est nécessaire pour permettre à la roche encais¬sante de produire les assemblages observés. L'observation des roches encaissantes entourant les granites montre que les minéraux clastiques dans les sédiments immatures au-dehors de l'auréole sont hydratés suite à la petite quantité de fluide expulsée durant le métamorphisme de contact et/ou la mise en place des granites. Les réactions d'hydratation peuvent permettre une augmentation de la température jusqu'à 50 °C. Afin de déterminer l'origine des fluides, une étude isotopique de roches de l'auréole de contact a été entreprise. Les isotopes stables d'oxygène et d'hydrogène sur la roche totale ainsi que la concentration en chlore dans la biotite indiquent que la mise en place des granites du Torres del Paine n'induit qu'une circulation de fluide limitée. Les données d'oxygène sur roche totale montrent des valeurs δ180 entre 9.0 et 10.0%o au sein des cinq premiers mètres du contact. Les valeurs augmentent jusqu'à 13.0 - 15.0 plus on s'éloigne de l'intrusion. Les valeurs 5D sur roche totale montrent une zonation plus complexe. Les valeurs de la roche encaissante (-90 à -70%o) diminuent progressivement d'environ 20%o depuis l'extérieur de l'auréole jusqu'à une distance d'environ 150 m du granite. Cette diminution est suivie par une augmentation d'environ 20%o au sein des 150 mètres les plus proches du contact (-97.0 à -78%o au contact). La diminution initiale des valeurs de 6D est interprétée comme la conséquence du fractionnement de Rayleigh qui accompagne les réactions de déshydratation formant la cordiérite, alors que l'augmentation finale reflète l'infiltration de fluide riche en eau venant de l'intrusion. A partir de ces résultats, le volume du fluide issu du granite ainsi que son effet thermique a pu être estimé. Ces résultats montrent que l'augmentation de température associée à ces fluides est limitée à un maximum de 30 °C. La contribution de ces fluides dans le bilan thermique est donc faible. Ces différents résultats nous ont permis de créer un modèle thermique associé à la for¬mation de l'auréole de contact qui intègre la mise en place rapide du granite et les réactions d'hydratation lors du métamorphisme. L'intégration de ce modèle thermique dans le modèle numérique de croissance minérale nous permet de calculer les textures des cordiérites. Cepen¬dant, ce modèle est dépendant de la vitesse de croissance et de nucléation de ces cordiérites. Nous avons obtenu ces paramètres en comparant les textures prédites par le modèle et les textures observées dans les roches de l'auréole de contact du Torres del Paine. Les paramètres cinétiques extraits du modèle optimisé indiquent une cristallisation rapide de la cordiérite avant le pic thermique dans la partie interne de l'auréole, et une réaction continue après le pic thermique dans la partie la plus externe de l'auréole. Seules de petites dépendances de température des paramètres de cinétique semblent être nécessaires pour expliquer les don¬nées obtenues sur la distribution des tailles de cristaux. Ces résultats apportent un éclairage nouveau sur la cinétique qui contrôle les réactions métamorphiques.