Geological transect across the Northwestern Himalaya in eastern Ladakh and Lahul (A model for the continental collision of India and Asia)

The detailed geological mapping and structural study of a complete transect across the northwestern Himalaya allow to describe the tectonic evolution of the north Indian continental margin during the Tethys ocean opening and the Himalayan Orogeny. The Late Paleozoic Tethys rifting is associated with several tectonomagmatic events. In Upper Lahul and SE Zanskar, this extensional phase is recorded by Lower Carboniferous synsedimentary transtensional faults, a Lower Permian stratigraphic unconformity, a Lower Permian granitic intrusion and middle Permian basaltic extrusions (Panjal Traps). In eastern Ladakh, a Permian listric normal fault is also related to this phase. The scarcity of synsedimentary faults and the gradual increase of the Permian syn-rift sediment thickness towards the NE suggest a flexural type margin. The collision of India and Asia is characterized by a succession of contrasting orogenic phases. South of the Suture Zone, the initiation of the SW vergent Nyimaling-Tsarap Nappe corresponds to an early phase of continental underthrusting. To the S, in Lahul, an opposite underthrusting within the Indian plate is recorded by the NE vergent Tandi Syncline. This structure is associated with the newly defined Shikar Beh Nappe, now partly eroded, which is responsible for the high grade (amphibolite facies) regional metamorphism of South Lahul. The main thrusting of the Nyimaling-Tsarap Nappe followed the formation of the Shikar Beh Nappe. The Nyimaling-Tsarap Nappe developed by ductile shear of the upper part of the subducted Indian continental margin and is responsible for the progressive regional metamorphism of SE Zanskar, reaching amphibolite facies below the frontal part of the nappe, near Sarchu. In Upper Lahul, the frontal parts of the Nyimaling-Tsarap and Shikar Beh nappes are separated by a zone of low grade metamorphic rocks (pumpellyite-actinolite facies to lower greenschist facies). At high structural level, the Nyimaling-Tsarap Nappe is characterized by imbricate structures, which grade into a large ductile shear zone with depth. The related crustal shortening is about 87 km. The root zone and the frontal part of this nappe have been subsequently affected by two zones of dextral transpression and underthrusting: the Nyimaling Shear Zone and the Sarchu Shear Zone. These shear zones are interpreted as consequences of the counterclockwise rotation of the continental underthrusting direction of India relative to Asia, which occurred some 45 and 36 Ma ago, according to plate tectonic models. Later, a phase of NE vergent `'backfolding'' developed on these two zones of dextral transpression, creating isoclinal folds in SE Zanskar and more open folds in the Nyimaling Dome and in the Indus Molasse sediments. During a late stage of the Himalayan Orogeny, the frontal part of the Nyimaling-Tsarap Nappe underwent an extension of about 15 km. This phase is represented by two types of structures, responsible for the tectonic unroofing of the amphibolite facies rocks of the Sarchu area: the Sarchu high angle Normal Fault, cutting a first set of low angle normal faults, which have been created by reactivation of older thrust planes related to the Nyimaling-Tsarap Nappe.

Physical considerations on discrepancies in target volume delineation.

BACKGROUND AND PURPOSE: To compare the delineations and interpretations of target volumes by physicians in different radio-oncology centers. MATERIALS AND METHODS: Eleven Swiss radio-oncology centers delineated volumes according to ICRU 50 recommendations for one prostate and one head and neck case. In order to evaluate the consistency of the volume delineations, the following parameters were determined: 1) the target volumes (GTV, CTV and manually expanded PTV) and their extensions in the three main axes and 2) the correlation of the volume delineated by each pair of centers using the ratio of the intersection to the union (called proximity index). RESULTS: The delineated prostate volume was 105+/-55cm(3) for the CTV and 218+/-44cm(3) for the PTV. The delineated head and neck volume was 46+/-15cm(3) for the GTV, 327+/-154cm(3) for the CTV and 528+/-106cm(3) for the PTV. The mean proximity index for the prostate case was 0.50+/-0.13 for the CTV and 0.57+/-0.11 for the PTV. The proximity index for the head and neck case was 0.45+/-0.09 for the GTV, 0.42+/-0.13 for the CTV and 0.59+/-0.06 for the PTV. CONCLUSIONS: Large discrepancies between all the delineated target volumes were observed. There was an inverse relationship between the CTV volume and the margin between CTV and PTV, leading to less discrepancies in the PTV than is the CTV delineations. There was more spread in the sagittal and frontal planes due to CT pixel anisotropy, which suggests that radiation oncologists should delineate the target volumes not only in the transverse plane, but also in the sagittal and frontal planes to improve the delineation by allowing a consistency check.

3D seismic facies characterization and geological patterns recognition (Australian North West Shelf)

EXECUTIVE SUMMARY This PhD research, funded by the Swiss Sciences Foundation, is principally devoted to enhance the recognition, the visualisation and the characterization of geobodies through innovative 3D seismic approaches. A series of case studies from the Australian North West Shelf ensures the development of reproducible integrated 3D workflows and gives new insight into local and regional stratigraphic as well as structural issues. This project was initiated in year 2000 at the Geology and Palaeontology Institute of the University of Lausanne (Switzerland). Several collaborations ensured the improvement of technical approaches as well as the assessment of geological models. - Investigations into the Timor Sea structural style were carried out at the Tectonics Special Research Centre of the University of Western Australia and in collaboration with Woodside Energy in Perth. - Seismic analysis and attributes classification approach were initiated with Schlumberger Oilfield Australia in Perth; assessments and enhancements of the integrated seismic approaches benefited from collaborations with scientists from Schlumberger Stavanger Research (Norway). Adapting and refining from "linear" exploration techniques, a conceptual "helical" 3D seismic approach has been developed. In order to investigate specific geological issues this approach, integrating seismic attributes and visualisation tools, has been refined and adjusted leading to the development of two specific workflows: - A stratigraphic workflow focused on the recognition of geobodies and the characterization of depositional systems. Additionally, it can support the modelling of the subsidence and incidentally the constraint of the hydrocarbon maturity of a given area. - A structural workflow used to quickly and accurately define major and secondary fault systems. The integration of the 3D structural interpretation results ensures the analysis of the fault networks kinematics which can affect hydrocarbon trapping mechanisms. The application of these integrated workflows brings new insight into two complex settings on the Australian North West Shelf and ensures the definition of astonishing stratigraphic and structural outcomes. The stratigraphic workflow ensures the 3D characterization of the Late Palaeozoic glacial depositional system on the Mermaid Nose (Dampier Subbasin, Northern Carnarvon Basin) that presents similarities with the glacial facies along the Neotethys margin up to Oman (chapter 3.1). A subsidence model reveals the Phanerozoic geodynamic evolution of this area (chapter 3.2) and emphasizes two distinct mode of regional extension for the Palaeozoic (Neotethys opening) and Mesozoic (abyssal plains opening). The structural workflow is used for the definition of the structural evolution of the Laminaria High area (Bonaparte Basin). Following a regional structural characterization of the Timor Sea (chapter 4.1), a thorough analysis of the Mesozoic fault architecture reveals a local rotation of the stress field and the development of reverse structures (flower structures) in extensional setting, that form potential hydrocarbon traps (chapter 4.2). The definition of the complex Neogene structural architecture associated with the fault kinematic analysis and a plate flexure model (chapter 4.3) suggest that the Miocene to Pleistocene reactivation phases recorded at the Laminaria High most probably result from the oblique normal reactivation of the underlying Mesozoic fault planes. This episode is associated with the deformation of the subducting Australian plate. Based on these results three papers were published in international journals and two additional publications will be submitted. Additionally this research led to several communications in international conferences. Although the different workflows presented in this research have been primarily developed and used for the analysis of specific stratigraphic and structural geobodies on the Australian North West Shelf, similar integrated 3D seismic approaches will have applications to hydrocarbon exploration and production phases; for instance increasing the recognition of potential source rocks, secondary migration pathways, additional traps or reservoir breaching mechanisms. The new elements brought by this research further highlight that 3D seismic data contains a tremendous amount of hidden geological information waiting to be revealed and that will undoubtedly bring new insight into depositional systems, structural evolution and geohistory of the areas reputed being explored and constrained and other yet to be constrained. The further development of 3D texture attributes highlighting specific features of the seismic signal, the integration of quantitative analysis for stratigraphic and structural processes, the automation of the interpretation workflow as well as the formal definition of "seismo-morphologic" characteristics of a wide range of geobodies from various environments would represent challenging examples of continuation of this present research. The 21st century will most probably represent a transition period between fossil and other alternative energies. The next generation of seismic interpreters prospecting for hydrocarbon will undoubtedly face new challenges mostly due to the shortage of obvious and easy targets. They will probably have to keep on integrating techniques and geological processes in order to further capitalise the seismic data for new potentials definition. Imagination and creativity will most certainly be among the most important quality required from such geoscientists.

Partial melting of basic dykes in the contact aureoles of two shallow-level, gabbro-pyroxenite intrusions, Fuerteventura

THESIS ABSTRACT Low-pressure anatexis of basic dykes gave rise to unusual, zebra-like migmatites, in the contact metamorphic aureoles of two layered gabbro-pyroxenite intrusions, PXl and PX2, in the root zone of an ocean island, Fuerteventura Basal Complex (Canary Islands). This thesis focuses on the understanding of processes attributing to the partial melting and formation of these migmatites, characterised by a dense network of closely spaced, millimetre-wide leucocratic segregations with perfectly preserved igneous textures. The presence of fluids are required to decrease the solidus of basic igneous lithologies, to allow partial melting in such aloes-pressure (1-2 kb) environment. An oxygen isotope study was thus carried out on dykes inside and beyond the PX2 aureole, in order to decipher the nature and origin of such fluids. Low or negative δ18O values were obtained for whole rocks and mineral-separates, decreasing towards the contact, with the intrusion itself retaining fairly high values. This trend has been attributed to the advection of meteoric water during magma emplacement, with increasing fluid/rock ratios (higher dyke intensities towards the intrusion acting as fluid-pathways) and higher temperatures promoting increasing exchange during recrystallisation. A comparison of whole rock and mineral major- and trace- element data allowed the redistribution of elements .between different mineral phases and generations, during contact metamorphism and partial melting to be assessed. Certain trace-elements, e.g. Zr, Hf, Y, and REEs, were internally redistributed during contact metamorphic recrystallisation, causing- the enrichment of neocrystallised diopsides compared to relict phenocrysts. This has been assigned to the liberation of trace elements on the breakdown of primary minerals, kaersutite and sphene, on entering the thermal aureole. Major and trace element compositions of minerals in migmatite melanosomes and leucosomes are almost identical, pointing to a syn- or post- solidus reequilibration on cooling of the migmatite terrain. The mineralogical, textural and geochemical evolution of dykes in a contact metamorphic aureole, is recorded around an apophysis of the PX1 intrusion, where there is evidence of incipient partial melting. Hydrothermal mineral pseudomorphs in the outer parts of the aureole are progressively replaced by dry mineral assemblages, with increasingly recrystallised diopside and evidence of partial melting -the extent of which varies from one lithology to another. The appearance of more mafic lithologies towards the intrusion, with lower whole rock SiO2 and mobile element abundances, e.g. Rb, Cs, K, has been explained by the migration and accumulation of feldspathic material into leucosomes outside the samples. A micro-structural study of leucosomes and leucocratic pods, with the aid of high-resolution X-ray computed micro-tomography (HRXµCT), allowing the visualization and quantification of shapes and orientations, was carried out in order to better understand the processes of melt segregation in the PX1 aureole. Leucocratic pods, representing former amygdales, are considered as natural strain ellipsoids. Their short axes are oriented perpendicular to leucosome planes, which sub-parallel the intrusive contact. Leucosomes thus effectively represent foliation planes. This implies that the direction of maximum shortening, during migmatisation, was perpendicular to the orientation of leucosomes, contradicting earlier models that suggest leucosomes represent tension veins. RESUME DE LA THESE Un phénomène rare de fusion partielle de filons basiques à basse pression a été étudié dans les auréoles de contact de deux intrusions litées de gabbro-pyroxénite, PX1 et PX2, localisées dans le soubassement de l'île volcanique de Fuerteventura aux Canaries. Cette anatexie a engendré des migmatites finement zébrées d'aspect très inhabituel, dont les processus de formation ont été étudiés dans le présent travail. Ces roches sont caractérisées par un réseau dense de veinules leucocrates d'épaisseur millimétrique, dont les textures ignées sont parfaitement préservées. La fusion partielle de roches basiques à basse pression (1-2 kbar) requiert la présence d'eau afin d'abaisser le solidus du système à des températures géologiquement réalistes. Une étude comparative des isotopes de l'oxygène a ainsi été menée sur des filons respectivement affectés et non affectés par le métamorphisme de contact, afin de confirmer la présence de ces fluides, de déterminer l'importance de leur interaction avec les roches et leur origine. Des valeurs de δ180 basses ou négatives ont été mesurées sur roche totale et minéraux séparés, décroissantes en direction du contact, alors que l'intrusion elle-même a conservé des valeurs élevées. Ce gradient a été attribué à l'advection d'eau météorique durant la mise en place du magma, les températures les plus élevées favorisant d'autant plus la circulation des fluides et les échanges isotopiques durant la recristallisation des roches. Cette recristallisation engendré une redistribution chimique complète des éléments entre les différentes générations de minéraux résultant du métamorphisme de contact et de l'anatexie, mise en évidence par microanalyse. Certains éléments traces comme Zr, Hf, Y et les REE ont été concentrés dans le diopside néoformé consécutivement à la déstabilisation de minéraux primaires riches en ces éléments comme la kaersutite ou le sphène. Les compositions en éléments majeurs et traces des minéraux des mélanosomes et leucosomes des migmatites sont pratiquement identiques, indiquant une rééquilibration syn- à postsolidus lors du refroidissement de l'auréole de contact. La transformation progressive des filons basiques au niveau de leur minéralogie, textures et composition chimique a pu être observée en détail à l'approche du contact d'une apophyse de l'intrusion PX1. La paragenèse magmatique initiale n'est jamais préservée, les faciès les plus distants du contact étant constitués d'un assemblage pseudomorphique hydrothermal. Ce dernier est progressivement remplacé par des assemblages anhydres incluant du diopside néoformé, puis apparaissent les premiers signes de fusion partielle, dont l'importance varie fortement d'une lithologie à l'autre. L'apparition de faciès plus basiques en direction du contact, avec des teneurs réduites en SiO2 et en éléments incompatibles tels Rb, Cs, K, a été attribuée à l'échappement de leucosomes feldspathiques hors du système. Une étude microstructurale de la distribution spatiale du matériel leucocrate au sein des migmatites par microtomographie X de haute résolution (HRXµCT) a été menée pour mieux comprendre les processus de ségrégation des liquides dans l'auréole de PX1. De petites entités ovoïdes, représentant d'anciennes structures amygdalaires au sein des filons, ont été considérées comme des ellipsoïdes marqueurs de la déformation finie. Leur petit axe est orienté perpendiculairement aux plans définis par les leucosomes, eux-mêmes subparallèles au contact intrusif. Les leucosomes matérialisent donc des plans de clivage. Ainsi, la direction de raccourcissement maximum durant la fusion partielle était perpendiculaire à l'orientation des leucosomes, contrairement à ce qui a été dit dans de précédentes publications, qui suggéraient que les leucosomes représentaient des veines de tension. RESUME DE LA THESE (POUR LE GRAND PUBLIC) L'observation directe du soubassement d'une île volcanique est une occasion rare, accessible dans le «complexe de base » de l'île canarienne de Fuerteventura. Ce dernier a enregistré divers phénomènes magmatiques, métamorphiques et de fusion partielle induits par l'intrusion répétée de magmas alimentant des appareils volcaniques sus jacents, sous forme de petits plutons, essaims de filons et complexes annulaires de gabbros alcalins, pyroxénites, syénites et carbonatites. Dans ce contexte de flux de chaleur élevé, des filons basiques ont subi une fusion partielle au contact de deux intrusions de gabbro-pyroxénite, un phénomène extrêmement rare à une profondeur aussi réduite, estimée à quelque 3-6 km. Les produits de cette fusion partielle sont des liquides très riches en feldspath, concentrés en un réseau dense de veinules blanches (leucosomes) au sein du matériau résiduel sombre non fondu (mélanosome) pour former ce qu'on appelle des migmatites. Outre les aspects pétrologiques liés à la formation de ces migmatites, l'intérêt majeur du phénomène réside dans le fait qu'il puisse représenter la source des magmas évolués parfois observés sur les îles océaniques. A des pressions aussi faibles que dans le soubassement de Fuerteventura (1-2 Kbar), la présence de fluides abondants est nécessaire pour abaisser la température de début de fusion des roches (solidus) à des valeurs géologiquement réalistes. Des expériences ont montré que même en présence de plusieurs %-poids d'eau, une température de 1000°C était encore nécessaire pour obtenir une proportion de liquide équivalente à celle observée sur le terrain, soit 25%. Or les magmas alcalins des îles océaniques, bien qu'hydratés, n'en contiennent de loin pas autant, ce qui implique une source d'eau externe. Une étude isotopique de l'oxygène a été entreprise afin de tester cette hypothèse. Les valeurs obtenues en 5180 sont basses ou négatives et indiquent l'influence d'eau d'origine météorique. Cette eau de pluie se serait infiltrée le long des filons depuis la surface du volcan et les aurait complètement hydrothermalisés en profondeur (situation encore visible à l'extérieur de l'auréole de contact), leur permettant ainsi de stocker l'eau nécessaire à leur fusion partielle ultérieure. L'interaction entre eau de pluie et filons a été d'autant plus importante que ces derniers étaient proches du contact avec l'intrusion, ce qui suggère que la circulation de ces eaux et leur interaction avec les roches a été favorisée par la chaleur fournie par l'intrusion elle-même. Un autre aspect de ce travail s'est focalisé sur la redistribution des éléments traces au sein des minéraux des filons basiques durant le métamorphisme de contact et la fusion partielle. Ainsi, le pyroxène de seconde génération est-il sensiblement enrichi en traces telles Zr, Hf, Y et les terres rares, par rapport au pyroxène magmatique originel, en relation avec la déstabilisation de minéraux primaires riches en ces éléments tels le sphène et la kaersutite. Cependant, les compositions en éléments majeurs et traces des minéraux recristallisés des migmatites sont pratiquement les mêmes dans les leucosomes et les mélanosomes, suggérant une rééquilibration chimique complète durant le refroidissement de ces lithologies. Si certaines migmatites se sont comportées en système chimiquement fermé (hormis l'eau météorique), d'autres filons ont manifestement perdu une partie de leurs leucosomes, ainsi qu'en témoigne leur composition progressivement appauvrie en silice et autres éléments incompatibles mobiles, tels K et Rb à l'approche du contact de l'intrusion. Parallèlement à cette évolution chimique, les paragenèses hydrothermales distantes du contact sont progressivement remplacées par des paragenèses anhydres, puis par l'apparition des premiers leucosomes, tandis que les textures magmatiques initiales sont complètement effacées au profit d'une combinaison de textures magmatique dans les leucosomes et en mosaïque dans les mélanosomes. Enfin, la distribution spatiale des liquides de fusion partielle a été étudiée par microtomographie X de haute résolution, dans des filons contenant des entités ovoïdes leucocrates, sans doute d'anciennes amygdales à zéolites. Ces dernières ont été considérées comme des ellipsoïdes de la déformation finie. L'orientation de leur petit axe, perpendiculaire au plan défini par les veinules de leucosomes, indique que ces derniers représentent des plans de clivage perpendiculaires à la direction de raccourcissement maximum. Ainsi, la ségrégation des liquides de fusion partielle se serait faite dans les plans de compression et non dans des plans de dilatation, contrairement à ce que laisserait penser le sens commun.

A Miocene mud volcano and its plumbing system: A chaotic complex revisited (Monferrato, NW Italy)

Chaotic deposits are frequently reported in the geological literature and are commonly interpreted as olistostromes or tectonic melanges. A chaotic complex in the Cenozoic succession of Monferrato (NW Italy) consists of interbedded mud breccia and burrowed silty clays that are pierced by sheared mud breccias and embed carbonate-cemented blocks. These may be represented by microcrystalline limestones or strongly cemented matrix-supported breccias locally containing remains of chemosymbiotic organisms (lucinid bivalves). Moreover, cylindrical concretions, up to 15 cm in diameter and 1 m long, occur in the chaotic complex and crosscut bedding planes at high angles. The cement of all these lithified portions is mainly dolomite characterized by low delta(13)C values (from -10.3 to -23parts per thousand PDB) and delta(18)O values up to + 7parts per thousand PDB. The delta(13)C values testify to precipitation of carbonates induced by microbial oxidation of methane, whereas the markedly positive delta(18)C signature, ubiquitous in the cylindrical concretions, is the evidence for the presence and destabilization of gas hydrates. The studied section provides a well-exposed example of the geological record of the birth, life, and death of a mud volcano. Unsheared, soft mud breccias represent mud flows along the flanks of the volcano, whereas sheared mud breccias are the result of the injection of unconsolidated overpressured fine-grained sediments, both taking place during ``eruptive'' phases. They were followed by more quiet stages of hemipelagic sedimentation, burrowing, and CH4 seeping. The cylindrical concretions represent the first described ancient example of the chimneys observed in present-day mud-volcano settings. They are the remnants of a cold-seep plumbing network that crosscut the mud volcano edifice. The chimneys were the pathway for the expulsion toward the sea floor of gas- and sediment-charged fluids likely originated from destabilization of methane gas hydrates. The association of mud breccias and methane-derived carbonates may not be due to mass gravity flows but can be primary and, therefore, is a diagnostic criterion for recognizing chaotic deposits due to mud volcano activity in the geological record.

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.

BIO-INSPIRED COMPUTATIONAL TECHNIQUES APPLIED TO THE CLUSTERING AND VISUALIZATION OF SPATIO-TEMPORAL GEOSPATIAL DATA

The coverage and volume of geo-referenced datasets are extensive and incessantly¦growing. The systematic capture of geo-referenced information generates large volumes¦of spatio-temporal data to be analyzed. Clustering and visualization play a key¦role in the exploratory data analysis and the extraction of knowledge embedded in¦these data. However, new challenges in visualization and clustering are posed when¦dealing with the special characteristics of this data. For instance, its complex structures,¦large quantity of samples, variables involved in a temporal context, high dimensionality¦and large variability in cluster shapes.¦The central aim of my thesis is to propose new algorithms and methodologies for¦clustering and visualization, in order to assist the knowledge extraction from spatiotemporal¦geo-referenced data, thus improving making decision processes.¦I present two original algorithms, one for clustering: the Fuzzy Growing Hierarchical¦Self-Organizing Networks (FGHSON), and the second for exploratory visual data analysis:¦the Tree-structured Self-organizing Maps Component Planes. In addition, I present¦methodologies that combined with FGHSON and the Tree-structured SOM Component¦Planes allow the integration of space and time seamlessly and simultaneously in¦order to extract knowledge embedded in a temporal context.¦The originality of the FGHSON lies in its capability to reflect the underlying structure¦of a dataset in a hierarchical fuzzy way. A hierarchical fuzzy representation of¦clusters is crucial when data include complex structures with large variability of cluster¦shapes, variances, densities and number of clusters. The most important characteristics¦of the FGHSON include: (1) It does not require an a-priori setup of the number¦of clusters. (2) The algorithm executes several self-organizing processes in parallel.¦Hence, when dealing with large datasets the processes can be distributed reducing the¦computational cost. (3) Only three parameters are necessary to set up the algorithm.¦In the case of the Tree-structured SOM Component Planes, the novelty of this algorithm¦lies in its ability to create a structure that allows the visual exploratory data analysis¦of large high-dimensional datasets. This algorithm creates a hierarchical structure¦of Self-Organizing Map Component Planes, arranging similar variables' projections in¦the same branches of the tree. Hence, similarities on variables' behavior can be easily¦detected (e.g. local correlations, maximal and minimal values and outliers).¦Both FGHSON and the Tree-structured SOM Component Planes were applied in¦several agroecological problems proving to be very efficient in the exploratory analysis¦and clustering of spatio-temporal datasets.¦In this thesis I also tested three soft competitive learning algorithms. Two of them¦well-known non supervised soft competitive algorithms, namely the Self-Organizing¦Maps (SOMs) and the Growing Hierarchical Self-Organizing Maps (GHSOMs); and the¦third was our original contribution, the FGHSON. Although the algorithms presented¦here have been used in several areas, to my knowledge there is not any work applying¦and comparing the performance of those techniques when dealing with spatiotemporal¦geospatial data, as it is presented in this thesis.¦I propose original methodologies to explore spatio-temporal geo-referenced datasets¦through time. Our approach uses time windows to capture temporal similarities and¦variations by using the FGHSON clustering algorithm. The developed methodologies¦are used in two case studies. In the first, the objective was to find similar agroecozones¦through time and in the second one it was to find similar environmental patterns¦shifted in time.¦Several results presented in this thesis have led to new contributions to agroecological¦knowledge, for instance, in sugar cane, and blackberry production.¦Finally, in the framework of this thesis we developed several software tools: (1)¦a Matlab toolbox that implements the FGHSON algorithm, and (2) a program called¦BIS (Bio-inspired Identification of Similar agroecozones) an interactive graphical user¦interface tool which integrates the FGHSON algorithm with Google Earth in order to¦show zones with similar agroecological characteristics.

Composition and growth mode of MoSx sputtered films

MoS(x) lubricating thin films were deposited by nonreactive, reactive, and low energy ion-assisted radio-frequency (rf) magnetron sputtering from a MoS2 target. Depending on the total and reactive gas pressures, the film composition ranges between MoS0.7 and MoS2.8. A low working pressure was found to have effects similar to those of low-energy ion irradiation. Films deposited at high pressure have (002) planes preferentially perpendicular to the substrate, whereas films deposited at low pressure or under low-energy ion irradiation have (002) mainly parallel to it. Parallel films are sulfur deficient (MoS1.2-1.4). Their growth is explained in terms of an increased reactivity of the basal surfaces, itself a consequence of the creation of surface defects due to ion irradiation. The films exhibit a lubricating character for all compositions above MoS1.2. The longest lifetime in ball-on-disk wear test was found for MoS1.5.

Detection of colorectal carcinoma by emission-computerized tomography after injection of 123I-labeled Fab or F(ab')2 fragments from monoclonal anti-carcinoembryonic antigen antibodies.

This clinical study was based on experimental results obtained in nude mice grafted with human colon carcinoma, showing that injected 131I-labeled F(ab')2 and Fab fragments from high affinity anti-carcinoembryonic antigen (CEA) monoclonal antibodies (MAb) gave markedly higher ratios of tumor to normal tissue localization than intact MAb. 31 patients with known colorectal carcinoma, including 10 primary tumors, 13 local tumor recurrences, and 21 metastatic involvements, were injected with 123I-labeled F(ab')2 (n = 14) or Fab (n = 17) fragments from MAb anti-CEA. The patients were examined by emission-computerized tomography (ECT) at 6, 24, and sometimes 48 h after injection using a rotating dual head scintillation camera. All 23 primary tumors and local recurrences except one were clearly visualized on at least two sections of different tomographic planes. Interestingly, nine of these patients had almost normal circulating CEA levels, and three of the visualized tumors weighed only 3-5 g. Among 19 known metastatic tumor involvements, 14 were correctly localized by ECT. Two additional liver and several bone metastases were discovered by immunoscintigraphy. Altogether, 86% of the tumor sites were detected, 82% with F(ab')2 and 89% with Fab fragments. The contrast of the tumor images obtained with Fab fragments suggests that this improved method of immunoscintigraphy has the potential to detect early tumor recurrences and thus to increase the survival of patients. The results of this retrospective study, however, should be confirmed in a prospective study before this method can be recommended for the routine diagnosis of cancer.

Epiphysenfugenverletzungen der distalen Tibia. Sinnvolle Mehrinformation durch die MRT? [Epiphyseal injuries of the distal tibia. Does MRI provide useful additional information?]

Plain film radiography often underestimates the extent of injury in children with epiphyseal fracture. Especially Salter-Harris V fractures (crush fracture of the epiphyseal plate) are often primarily not detected. MRI of the ankle was performed in 10 children aged 9-17 (mean 14) years with suspected epiphyseal injury using 1.0-T Magnetom Expert. The fractures were classified according to the Salter-Harris-Rang-Odgen classification and compared with the results of plain radiography. In one case MRI could exclude epiphyseal injury; in four cases the MRI findings changed the therapeutic management. The visualisation of the fracture in three orthogonal planes and the possibility of detection of cartilage and ligamentous injury in MR imaging makes this method superior to conventional radiography and CT. With respect to radiation exposure MRI instead of CT should be used for the diagnosis of epiphyseal injuries in children.

The 2006 Eiger rockslide, European Alps

In July 2006, approximately 2 million m3 of massive limestone began to move on the east flank of the Eiger in central Switzerland. For more than two years after the initial failure, the rock mass moved at rates of up to 70 cm per day. A detailed analysis of the structures and velocities of the different moving blocks was conducted with the aid of terrestrial laser scanning. The moving rock mass included a rear block that subsided, pushing a frontal block forward. Movement directions were controlled by discontinuity sets that formed wedges bounded on one side by sub-vertical bedding planes. The instability was, until recently, buttressed by a glacier. Slope observations and results of continuum and discontinuum modeling indicate that the structure of the rock mass and topography were the main causes of the instability. Progressive weathering and mechanical fatigue of the rock mass appear to have led to the failure. A dynamic analytical model further indicates that the rockslide was primarily controlled by a reduction in the strength of discontinuities, the effects of ice deformation, and ? to a limited extent ? groundwater flow. This study shows that realistic and simple instability models can be constructed for rock-slope failures if high-resolution data are available.

Neutron structure of type-III antifreeze protein allows the reconstruction of AFP-ice interface.

Antifreeze proteins (AFPs) inhibit ice growth at sub-zero temperatures. The prototypical type-III AFPs have been extensively studied, notably by X-ray crystallography, solid-state and solution NMR, and mutagenesis, leading to the identification of a compound ice-binding surface (IBS) composed of two adjacent ice-binding sections, each which binds to particular lattice planes of ice crystals, poisoning their growth. This surface, including many hydrophobic and some hydrophilic residues, has been extensively used to model the interaction of AFP with ice. Experimentally observed water molecules facing the IBS have been used in an attempt to validate these models. However, these trials have been hindered by the limited capability of X-ray crystallography to reliably identify all water molecules of the hydration layer. Due to the strong diffraction signal from both the oxygen and deuterium atoms, neutron diffraction provides a more effective way to determine the water molecule positions (as D(2) O). Here we report the successful structure determination at 293 K of fully perdeuterated type-III AFP by joint X-ray and neutron diffraction providing a very detailed description of the protein and its solvent structure. X-ray data were collected to a resolution of 1.05 Å, and neutron Laue data to a resolution of 1.85 Å with a "radically small" crystal volume of 0.13 mm(3). The identification of a tetrahedral water cluster in nuclear scattering density maps has allowed the reconstruction of the IBS-bound ice crystal primary prismatic face. Analysis of the interactions between the IBS and the bound ice crystal primary prismatic face indicates the role of the hydrophobic residues, which are found to bind inside the holes of the ice surface, thus explaining the specificity of AFPs for ice versus water.

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

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

Measurement of multi-segment foot joint angles during gait using a wearable system.

Usually the measurement of multi-segment foot and ankle complex kinematics is done with stationary motion capture devices which are limited to use in a gait laboratory. This study aimed to propose and validate a wearable system to measure the foot and ankle complex joint angles during gait in daily conditions, and then to investigate its suitability for clinical evaluations. The foot and ankle complex consisted of four segments (shank, hindfoot, forefoot, and toes), with an inertial measurement unit (3D gyroscopes and 3D accelerometers) attached to each segment. The angles between the four segments were calculated in the sagittal, coronal, and transverse planes using a new algorithm combining strap-down integration and detection of low-acceleration instants. To validate the joint angles measured by the wearable system, three subjects walked on a treadmill for five minutes at three different speeds. A camera-based stationary system that used a cluster of markers on each segment was used as a reference. To test the suitability of the system for clinical evaluation, the joint angle ranges were compared between a group of 10 healthy subjects and a group of 12 patients with ankle osteoarthritis, during two 50-m walking trials where the wearable system was attached to each subject. On average, over all joints and walking speeds, the RMS differences and correlation coefficients between the angular curves obtained using the wearable system and the stationary system were 1 deg and 0.93, respectively. Moreover, this system was able to detect significant alteration of foot and ankle function between the group of patients with ankle osteoarthritis and the group of healthy subjects. In conclusion, this wearable system was accurate and suitable for clinical evaluation when used to measure the multi-segment foot and ankle complex kinematics during long-distance walks in daily life conditions.

Preliminary description and slope stability analyses of the 2008 Little Salmon lake and 2007 Mt. Steele landslides, Yukon

In August 2008, reactivation of the Little Salmon Lake landslide occurred. During this event, hundreds of conical mounds of variable size and composition formed in the deposition zone. The characteristics of these landforms are described and a potential mechanism for their formation is proposed. A preliminary slope stability analysis of the 2007 Mount Steele rock and ice avalanche was also undertaken. The orientation of very high persistence (>20 m long) structural planes (e.g., faults, joints and bedding) within bedrock in the source zone was obtained using an airborne-LiDAR digital elevation model and the software COLTOP-3D. Using these discontinuity orientation measurements, kinematic, surface wedge and simple three-dimensional distinct element slope stability analyses were performed.