Dose reconstruction in the context of tomotherapy

In vivo dosimetry is a way to verify the radiation dose delivered to the patient in measuring the dose generally during the first fraction of the treatment. It is the only dose delivery control based on a measurement performed during the treatment. In today's radiotherapy practice, the dose delivered to the patient is planned using 3D dose calculation algorithms and volumetric images representing the patient. Due to the high accuracy and precision necessary in radiation treatments, national and international organisations like ICRU and AAPM recommend the use of in vivo dosimetry. It is also mandatory in some countries like France. Various in vivo dosimetry methods have been developed during the past years. These methods are point-, line-, plane- or 3D dose controls. A 3D in vivo dosimetry provides the most information about the dose delivered to the patient, with respect to ID and 2D methods. However, to our knowledge, it is generally not routinely applied to patient treatments yet. The aim of this PhD thesis was to determine whether it is possible to reconstruct the 3D delivered dose using transmitted beam measurements in the context of narrow beams. An iterative dose reconstruction method has been described and implemented. The iterative algorithm includes a simple 3D dose calculation algorithm based on the convolution/superposition principle. The methodology was applied to narrow beams produced by a conventional 6 MV linac. The transmitted dose was measured using an array of ion chambers, as to simulate the linear nature of a tomotherapy detector. We showed that the iterative algorithm converges quickly and reconstructs the dose within a good agreement (at least 3% / 3 mm locally), which is inside the 5% recommended by the ICRU. Moreover it was demonstrated on phantom measurements that the proposed method allows us detecting some set-up errors and interfraction geometry modifications. We also have discussed the limitations of the 3D dose reconstruction for dose delivery error detection. Afterwards, stability tests of the tomotherapy MVCT built-in onboard detector was performed in order to evaluate if such a detector is suitable for 3D in-vivo dosimetry. The detector showed stability on short and long terms comparable to other imaging devices as the EPIDs, also used for in vivo dosimetry. Subsequently, a methodology for the dose reconstruction using the tomotherapy MVCT detector is proposed in the context of static irradiations. This manuscript is composed of two articles and a script providing further information related to this work. In the latter, the first chapter introduces the state-of-the-art of in vivo dosimetry and adaptive radiotherapy, and explains why we are interested in performing 3D dose reconstructions. In chapter 2 a dose calculation algorithm implemented for this work is reviewed with a detailed description of the physical parameters needed for calculating 3D absorbed dose distributions. The tomotherapy MVCT detector used for transit measurements and its characteristics are described in chapter 3. Chapter 4 contains a first article entitled '3D dose reconstruction for narrow beams using ion chamber array measurements', which describes the dose reconstruction method and presents tests of the methodology on phantoms irradiated with 6 MV narrow photon beams. Chapter 5 contains a second article 'Stability of the Helical TomoTherapy HiArt II detector for treatment beam irradiations. A dose reconstruction process specific to the use of the tomotherapy MVCT detector is presented in chapter 6. A discussion and perspectives of the PhD thesis are presented in chapter 7, followed by a conclusion in chapter 8. The tomotherapy treatment device is described in appendix 1 and an overview of 3D conformai- and intensity modulated radiotherapy is presented in appendix 2. - La dosimétrie in vivo est une technique utilisée pour vérifier la dose délivrée au patient en faisant une mesure, généralement pendant la première séance du traitement. Il s'agit de la seule technique de contrôle de la dose délivrée basée sur une mesure réalisée durant l'irradiation du patient. La dose au patient est calculée au moyen d'algorithmes 3D utilisant des images volumétriques du patient. En raison de la haute précision nécessaire lors des traitements de radiothérapie, des organismes nationaux et internationaux tels que l'ICRU et l'AAPM recommandent l'utilisation de la dosimétrie in vivo, qui est devenue obligatoire dans certains pays dont la France. Diverses méthodes de dosimétrie in vivo existent. Elles peuvent être classées en dosimétrie ponctuelle, planaire ou tridimensionnelle. La dosimétrie 3D est celle qui fournit le plus d'information sur la dose délivrée. Cependant, à notre connaissance, elle n'est généralement pas appliquée dans la routine clinique. Le but de cette recherche était de déterminer s'il est possible de reconstruire la dose 3D délivrée en se basant sur des mesures de la dose transmise, dans le contexte des faisceaux étroits. Une méthode itérative de reconstruction de la dose a été décrite et implémentée. L'algorithme itératif contient un algorithme simple basé sur le principe de convolution/superposition pour le calcul de la dose. La dose transmise a été mesurée à l'aide d'une série de chambres à ionisations alignées afin de simuler la nature linéaire du détecteur de la tomothérapie. Nous avons montré que l'algorithme itératif converge rapidement et qu'il permet de reconstruire la dose délivrée avec une bonne précision (au moins 3 % localement / 3 mm). De plus, nous avons démontré que cette méthode permet de détecter certaines erreurs de positionnement du patient, ainsi que des modifications géométriques qui peuvent subvenir entre les séances de traitement. Nous avons discuté les limites de cette méthode pour la détection de certaines erreurs d'irradiation. Par la suite, des tests de stabilité du détecteur MVCT intégré à la tomothérapie ont été effectués, dans le but de déterminer si ce dernier peut être utilisé pour la dosimétrie in vivo. Ce détecteur a démontré une stabilité à court et à long terme comparable à d'autres détecteurs tels que les EPIDs également utilisés pour l'imagerie et la dosimétrie in vivo. Pour finir, une adaptation de la méthode de reconstruction de la dose a été proposée afin de pouvoir l'implémenter sur une installation de tomothérapie. Ce manuscrit est composé de deux articles et d'un script contenant des informations supplémentaires sur ce travail. Dans ce dernier, le premier chapitre introduit l'état de l'art de la dosimétrie in vivo et de la radiothérapie adaptative, et explique pourquoi nous nous intéressons à la reconstruction 3D de la dose délivrée. Dans le chapitre 2, l'algorithme 3D de calcul de dose implémenté pour ce travail est décrit, ainsi que les paramètres physiques principaux nécessaires pour le calcul de dose. Les caractéristiques du détecteur MVCT de la tomothérapie utilisé pour les mesures de transit sont décrites dans le chapitre 3. Le chapitre 4 contient un premier article intitulé '3D dose reconstruction for narrow beams using ion chamber array measurements', qui décrit la méthode de reconstruction et présente des tests de la méthodologie sur des fantômes irradiés avec des faisceaux étroits. Le chapitre 5 contient un second article intitulé 'Stability of the Helical TomoTherapy HiArt II detector for treatment beam irradiations'. Un procédé de reconstruction de la dose spécifique pour l'utilisation du détecteur MVCT de la tomothérapie est présenté au chapitre 6. Une discussion et les perspectives de la thèse de doctorat sont présentées au chapitre 7, suivies par une conclusion au chapitre 8. Le concept de la tomothérapie est exposé dans l'annexe 1. Pour finir, la radiothérapie «informationnelle 3D et la radiothérapie par modulation d'intensité sont présentées dans l'annexe 2.

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.

Petrology and chemistry of the Chalten Plutonic Complex and implications on the magmatic and tectonic evolution of the Southernmost Andes (Patagonia) during the Miocene

Résumé Scientific:Pétrologie et Géochimie du Complexe Plutonique de Chaltén et les conséquences pour l'évolution magmatique et tectonique du Andes du Sud (Patagonia) pendant le MiocèneLe sujet de cette thèse est le Complexe Plutonique de Chaltén (CHPC), situé à la frontière entre le Chili et l'Argentine, en Patagonie (49°15'S). Ce complexe s'est mis en place au début du Miocène, dans un contexte de changements tectoniques importants. La géométrie et la vitesse de migration des plaques en Patagonie a été modifiée suite l'ouverture de la plaque Farallon il y a 25Ma (Pardo-Casas and Molnar 1987) et la subduction de la ride active du Chili sous la plaque sud-américaine il y a 14Ma (Cande and Leslie 1986). Les effets de cette reconfiguration tectonique sur la morphologie et le magmatisme de la plaque supérieure sont encore sujets à discussion. Dans ce contexte, un groupe d'intrusions miocènes - telle que le CHPC - est particulièrement intriguant, car en position transitionnelle entre le batholithe patagonien et l'arc volcanique cénozoïque et récent à l'ouest, et les laves de plateau de Patagonie à l'est (Fig. 1). A cause de leur position tectonique transitoire, ces plutons isolés hors du batholithe représentent un endroit clé pour comprendre les interactions entre la tectonique à large échelle et le magmatisme en Patagonie. Ici, je présente de nouvelles données de terrain, petrologiques, géochimiques et géochronologiques dans le but de caractériser la nature du CHPC, qui était largement inconnu avant cette étude, dans le but de tester l'hypothèse de migration de l'arc et erosion par subduction.Les résultats de l'investigation géochimique (chapitre 2) montrent que le CHPC n'est qu'un exemple parmi les plutons isolés d'arrière arc ave une composition calco-alcaline caractéristique, c-à-d une signature d'arc. La plupart de ces plutons isolés ont une composition alcaline. Le CHPC, contrairement, a une signature calco-alcaline avec Κ intermédiaire, tel que le batholithe patagonien et la plupart des roches volcaniques quaternaires liées à l'arc le long des Andes.De nouvelles données géochronologiques U-Pb de haute précision sur des zircons, acquis par TIMS, sur le CHPC donnent des âges entre 17.0 et 16.4Ma. Les âges absolus sont en accord avec la séquence intrusive déduite des relations de terrain (chapitre 1). Ces données sont les premières contraintes d'âge U-Pb sur le CHPC. Elles montrent clairement que l'histoire magmatique du CHPC n'a pas de lien direct avec la subduction de la ride à cette latitude (Cande and Leslie 1986), car le complexe est au moins 6Ma plus ancien.Une comparaison en profondeur avec les autres intrusions d'âge Miocène en Patagonie révèlent - pour la première fois - une évolution temporelle intéressante. Il y a une tendance E-W distincte au magmatisme calco-alcalin entre 20-16Ma avec une diminution de l'âge vers l'est - le CHPC est l'expression la plus orientale de cette tendance. Je suggère que la relation espace-temps reflète une migration vers l'est (vers le continent) de l'arc magmatique. Je propose que le facteur principal contrôlant cette migration est la subduction rapide suite à la reconfiguration de la vitesse des plaques tectoniques après l'ouverture la plaque Farallon (à ~26Ma) qui résulterait en une déformation importante ainsi qu'à des taux élevés d'érosion dans la fosse de subduction.Les rapports d'isotopes radiogéniques (Pb, Sr, Nd) élevés, une signature 6018 basse et un rapport Th/La élevé sont des paramètres distinctifs pour les roches mafiques du CHPC. Le modèle isotopique présenté (chapitre 2) suggère que cette signature reflète une contamination de la source, dans le coin de manteau, plutôt qu'une contamination crustale. La signature des éléments en trace du CHPC indiquent que le coin de manteau a été contaminé par des composés terrigènes, le plus vraisemblablement par des sédiments paléozoïques.Les travaux de terrain, la pétrographie et la géothermobarométrie ont été utilisés dans le but de comprendre l'histoire interne du CHPC (chapitre 3). Ces données suggèrent deux niveaux distincts de cristallisation : l'un dans la croûte moyenne (6 à 4.5kbar) et l'autre à un niveau peu profond (3.5 à 2kbar). La modélisation isotopique AFC de la contamination crustale indique des taux variables d'assimilation, qui ne sont pas corrélés avec le degré de différenciation. Cela suggère que différents volumes de magma se sont différenciés en profondeur, de façon indépendante. Cela implique que le CHPC se serait formés en plusieurs puises de magmas provenant d'au moins trois sources différentes. Les textures des granodiorites et des granites indiquent des teneurs élevées en cristaux avant la mise en place et, par conséquent, des températures d'emplacement faibles. Les observations de terrain montrent que les roches mafiques sont déformées, alors que ce n'est pas le cas pour les granodiorites et granites (plus jeunes). La déformation des roches mafiques est encore sujet de recherche, afin de savoir si elle est liée à la déformation régionale en régime compressif ou à l'emplacement lui-même. Cependant, la mise en place de grand volume de magma felsique riche en cristaux suggère un régime d'extension.Scientific Abstract:Petrology and chemistry of the Chaltén Plutonic Complex and implications on the magmatic and tectonic evolution of the Southernmost Andes (Patagonia) during the MioceneThe subject of this thesis is the Chaltén Plutonic Complex (CHPC) located at the frontier between Chile and Argentina in Patagonia (at 49° 15 'Southern latitude). This complex intruded during early Miocene in a context of major tectonics changes. The plate geometry of Patagonia has been modified by changes in the plate motions after the break up of the Farallôn plate at 25Ma (Pardo-Casas and Molnar 1987) and by the subduction of the Chile spreading Ridge beneath South-America at 14 Ma (Cande and Leslie 1986). The effects of this tectonic setting on the morphology and the magmatism of the overriding plate are a matter of on-going discussion. Particularly intriguing in this context is a group of isolated Miocene intrusions - like the CHPC - which are located in a transitional position between the Patagonian Batholith and the Cenozoic and Recent volcanic arc in the West, and the Patagonian plateau lavas in the East (Fig. 1). Due to their transient tectonic position these isolated plutons outside the batholith represent a key to understanding the interaction between global-scale tectonics and magmatism in Patagonia. Here, I present new field, penological, geochemical and geochronological data to characterize the nature of the CHPC, which was largely unknown before this study, in order to test the hypothesis of time- transgressive magmatism.The results of the geochemical investigation (Chapter 2) show that the CHPC is only one among these isolated back-arc plutons with a characteristic calc-alkaline composition, i.e. arc signature. Most of these isolated intrusives have an alkaline character. The CHPC, in contrast, has a medium Κ calc-alkaline signature, like the Patagonian batholith and most of the Quaternary arc-related volcanic rocks along the Andes.New high precision TIMS U-Pb zircon dating of the CHPC yield ages between 17.0 to 16.4 Ma. The absolute ages support the sequence of intrusion relations established in the field (Chapter 1). These data are the first U-Pb age constraints on the CHPC, and clearly show that the magmatic history of CHPC has no direct link to the subduction of the ridge, since this complex is at least 6 Ma older than the time of collision of the Chile ridge at this latitude (Cande and Leslie 1986).An in-depth comparison with other intrusion of Miocene age in Patagonia reveals - for the first time - an interesting temporal pattern. There is a distinct E-W trend of calc-alkaline magmatism between 20-16 Ma with the younging of ages in the East - the CHPC is the easternmost expression of this trend. I suggest that this time-space relation reflects an eastward (landward) migration of the magmatic arc. I propose that main factor controlling this migration is the fast rates of subduction after the major reconfigurations of the plate tectonic motions after the break up of the Farallôn Plate (at -26 ) resulting in strong deformation and high rates of subduction erosion.High radiogenic isotope ratios (Pb, Sr, Nd) ratios, low 5018 signature and high Th/La ratios in mafic rocks are distinctive features of the CHPC. The presented isotopic models (Chapter 2) suggest that this signature reflects source contamination of the mantle wedge rather than crustal contamination. The trace element signature of the CHPC indicates that the mantle wedge was contaminated with a terrigenous component, most likely from Paleozoic sediments.Fieldwork, petrography and geothermobarometry were used to further unravel the internal history of the CHPC (Chapter 3). These data suggest two main levels of crystallization: one a mid crustal levels (6 to 4.5 kbar) and other a shallow level (3.5 to 2 kbar). Isotopic AFC modeling of crustal contamination indicate variable rates of assimilation, which are not correlated with the degree of differentiation. This suggests that different batches of magma differentiate independently at depths. This implies that the CHPC would have formed by several pulses of magmas from at least 3 different sources. Textures of granodiorites and granites indicate a high content of crystals previous to the emplacement and consequently low emplacement temperatures. Field observations show that the mafic rocks are deformed, whereas the (younger) granodiorites and granites are not. It is still subject of investigation whether the deformation of the mafic rocks is related to regional deformation during a compressional regime or to the emplacement it self. However, the emplacement of huge amount of crystal rich felsic magmas suggests an extensional regime.Résumé Grand PublicPétrologie et Géochimie du Complexe Plutonique de Chaltén et les conséquences pour l'évolution magmatique et tectonique du Andes du Sud (Patagonia) pendant le MiocèneLe Complexe Plutonique de Chaltén (CHPC) est un massif montagneux situé à 49°S à la frontière entre le Chili et l'Argentine, en Patagonie (région la plus au sud de l'Amérique du Sud). Il est composé de montagnes qui peuvent atteindre plus de 3000 mètres d'altitude, telles que le Cerro Fitz Roy (3400m) et le Cerro Torre (3100m). Ces montagnes sont composées de roches plutoniques, c.-à-d. des magmas qui se sont refroidis et ont cristallisés sous la surface terrestre.La composition chimique de ces roches montre que les magmas, qui ont formé ce complexe plutonique, font partie d'un volcanisme d'arc. Celui-ci se forme lorsqu'une plaque océanique plonge sous une plaque continentale. Les géologues appellent ce processus « subduction ». Dans un tel scénario, le manteau terrestre, qui se fait prendre entre ces deux plaques, fond pour former ainsi du magma. Ce magma remonte à travers la plaque continentale vers la surface. Si celui-ci atteint la surface, il forme les roches volcaniques, comme par exemple des laves. S'il n'atteint pas la surface, le magma se refroidit pour former finalement les roches plutoniques.Le long de la marge ouest d'Amérique du Sud, la plaque Nazca - qui se situe au sud-est de la plaque océanique pacifique - passe en dessous de la plaque d'Amérique du Sud. La bordure ouest du sud de la plaque sud-américaine a également été affectée par d'autres processus tectoniques, tels que des changements dramatiques dans les déplacements de plaques (il y a 25Ma) et la collision de la ride du Chili (depuis 15 Ma jusqu'à aujourd'hui). Ces caractéristiques tectoniques et magmatiques font de cette région un haut lieu pour les géologues. La plaque Nazca, s'est formée suite à l'ouverture d'une plaque océanique plus ancienne, il y a 25Ma. Cette ouverture est liée aux vitesses de subduction les plus rapides jamais connues. La ride du Chili est l'endroit où le sol de l'Océan Pacifique s'ouvre, formant deux plaques océaniques : les plaques Nazca et Antarctique. La ride du Chili subducte sous la plaque sud-américaine depuis 15Ma, en association avec la formation de grands volumes de magma ainsi que des changements morphologiques importants. La question de savoir lequel de ces changements tectoniques globaux affecte la géologie et la géographie de Patagonie a été, et est encore, discutée pendant de nombreuses années. De nombreux chercheurs suggèrent que la plupart des caractéristiques morphologiques et magmatiques en Patagonie sont liés à la subduction de la ride du Chili, mais cette suggestion est encore débattue comme le montre notre étude.Le batholithe de Patagonie du sud (SPB) est un énorme massif composé de roches plutoniques et il s'étend tout au long de la côte ouest de Patagonie (au sud de 47°S). Ces roches correspondent certainement aux racines d'un ancien arc volcanique, qui a été soulevé et érodé. Le CHPC, ainsi que d'autres petites intrusions dans la région, se situe dans une position exotique, à 100km à l'est du SPB. Certains chercheurs suggèrent que ces intrusions pourraient être liées à la subduction de la ride du Chili.Afin de débattre de cette problématique, nous avons utilisé différentes méthodes géochronologiques pour déterminer l'âge du CHPC et le comparer (a) à l'âge des roches intrusives similaires du SPB et (b) à l'âge de la collision de la ride du Chili. Dans ce travail, nous prouvons que le CHPC s'est formé au moins 7Ma avant la collision avec la ride du Chili. Sur la base des âges du CHPC et de la composition chimique de ses roches et minéraux, nous proposons que le CHPC fait partie d'un arc volcanique ancien. La migration de l'arc volcanique plus profondément dans le continent résulte de la grande vitesse de subduction entre 25 et lOMa. Des caractéristiques évidentes pour un tel processus - telles qu'une déformation importante et une vitesse d'érosion élevée - peuvent être rencontrées tout au long de la bordure ouest de l'Amérique du sud.

Right coronary artery flow velocity and volume assessment with spiral K-space sampled breathhold velocity-encoded MRI at 3 tesla: accuracy and reproducibility.

PURPOSE: To evaluate accuracy and reproducibility of flow velocity and volume measurements in a phantom and in human coronary arteries using breathhold velocity-encoded (VE) MRI with spiral k-space sampling at 3 Tesla. MATERIALS AND METHODS: Flow velocity assessment was performed using VE MRI with spiral k-space sampling. Accuracy of VE MRI was tested in vitro at five constant flow rates. Reproducibility was investigated in 19 healthy subjects (mean age 25.4 +/- 1.2 years, 11 men) by repeated acquisition in the right coronary artery (RCA). RESULTS: MRI-measured flow rates correlated strongly with volumetric collection (Pearson correlation r = 0.99; P < 0.01). Due to limited sample resolution, VE MRI overestimated the flow rate by 47% on average when nonconstricted region-of-interest segmentation was used. Using constricted region-of-interest segmentation with lumen size equal to ground-truth luminal size, less than 13% error in flow rate was found. In vivo RCA flow velocity assessment was successful in 82% of the applied studies. High interscan, intra- and inter-observer agreement was found for almost all indices describing coronary flow velocity. Reproducibility for repeated acquisitions varied by less than 16% for peak velocity values and by less than 24% for flow volumes. CONCLUSION: 3T breathhold VE MRI with spiral k-space sampling enables accurate and reproducible assessment of RCA flow velocity.

Geometry-based demosaicking

Demosaicking is a particular case of interpolation problems where, from a scalar image in which each pixel has either the red, the green or the blue component, we want to interpolate the full-color image. State-of-the-art demosaicking algorithms perform interpolation along edges, but these edges are estimated locally. We propose a level-set-based geometric method to estimate image edges, inspired by the image in-painting literature. This method has a time complexity of O(S) , where S is the number of pixels in the image, and compares favorably with the state-of-the-art algorithms both visually and in most relevant image quality measures.

Asymmetry of Functional Connectivity in Schizophrenia at Different Spatial Scales

Introduction: The interhemispheric asymmetries that originate from connectivity-related structuring of the cerebral cortex are compromised in schizophrenia (SZ). Recently, we have revealed the whole-head topography of EEG synchronization in SZ (Jalili et al. 2007; Knyazeva et al. 2008). Here we extended the analysis to assess the abnormality in the asymmetry of synchronization, which is further motivated by the evidence that the interhemispheric asymmetries suspected to be abnormal in SZ originate from the connectivity-related structuring of the cortex. Methods: Thirteen right-handed SZ patients and thirteen matched controls, participated in this study and the multichannel (128) EEGs were recorded for 3-5 minutes at rest. Then, Laplacian EEG (LEEG) were calculated using a 2-D spline. The LEEGs were analysis through calculating the power spectral density using Welch's average periodogram method. Furthermore, using a state-space based multivariate synchronization measure, S-estimator, we analyzed the correlate of the functional cortico-cortical connectivity in SZ patients compared to the controls. The values of S-estimator were obtained at three different special scales: first-order neighbors for each sensor location, second-order neighbors, and the whole hemisphere. The synchronization measures based on LEEG of alpha and beta bands were applied and tuned to various spatial scales including local, intraregional, and long-distance levels. To assess the between-group differences, we used a permutation version of Hotelling's T2 test. For correlation analysis, Spearman Rank Correlation was calculated. Results: Compared to the controls, who had rightward asymmetry at a local level (LEEG power), rightward anterior and leftward posterior asymmetries at an intraregional level (first- and second-order S-estimator), and rightward global asymmetry (hemispheric S-estimator), SZ patients showed generally attenuated asymmetry, the effect being strongest for intraregional synchronization. This deviation in asymmetry across the anterior-to-posterior axis is consistent with the cerebral form of the so-called Yakovlevian or anticlockwise cerebral torque. Moreover, the negative occipital and positive frontal asymmetry values suggest higher regional synchronization among the left occipital and the right frontal locations relative to their symmetrical counterparts. Correlation analysis linked the posterior intraregional and hemispheric abnormalities to the negative SZ symptoms, whereas the asymmetry of LEEG power appeared to be weakly coupled to clinical ratings. The posterior intraregional abnormalities of asymmetry were shown to increase with the duration of the disease. The tentative links between these findings and gross anatomical asymmetries, including the cerebral torque and gyrification pattern in normal subjects and SZ patients, are discussed. Conclusions: Overall, our findings reveal the abnormalities in the synchronization asymmetry in SZ patients and heavy involvement of the right hemisphere in these abnormalities. These results indicate that anomalous asymmetry of cortico-cortical connections in schizophrenia is amenable to electrophysiological analysis.

On fast-decodable space-time block codes

We focus on full-rate, fast-decodable space–time block codes (STBCs) for 2 x 2 and 4 x 2 multiple-input multiple-output (MIMO) transmission. We first derive conditions and design criteria for reduced-complexity maximum-likelihood (ML) decodable 2 x 2 STBCs, and we apply them to two families of codes that were recently discovered. Next, we derive a novel reduced-complexity 4 x 2 STBC, and show that it outperforms all previously known codes with certain constellations.

On high-rate full-diversity 2x2 space-time codes with low-complexity optimum detection

The 2×2 MIMO profiles included in Mobile WiMAX specifications are Alamouti’s space-time code (STC) fortransmit diversity and spatial multiplexing (SM). The former hasfull diversity and the latter has full rate, but neither of them hasboth of these desired features. An alternative 2×2 STC, which is both full rate and full diversity, is the Golden code. It is the best known 2×2 STC, but it has a high decoding complexity. Recently, the attention was turned to the decoder complexity, this issue wasincluded in the STC design criteria, and different STCs wereproposed. In this paper, we first present a full-rate full-diversity2×2 STC design leading to substantially lower complexity ofthe optimum detector compared to the Golden code with only a slight performance loss. We provide the general optimized form of this STC and show that this scheme achieves the diversitymultiplexing frontier for square QAM signal constellations. Then, we present a variant of the proposed STC, which provides a further decrease in the detection complexity with a rate reduction of 25% and show that this provides an interesting trade-off between the Alamouti scheme and SM.

A comparison of full-rate full-diversity 2x2 space-time codes for WiMAX systems

Multiple-input multiple-output (MIMO) techniques have become an essential part of broadband wireless communications systems. For example, the recently developed IEEE 802.16e specifications for broadband wireless access include three MIMOprofiles employing 2×2 space-time codes (STCs), and two of these MIMO schemes are mandatory on the downlink of Mobile WiMAX systems. One of these has full rate, and the other has full diversity, but neither of them has both of the desired features. The third profile, namely, Matrix C, which is not mandatory, is both a full rate and a full diversity code, but it has a high decoder complexity. Recently, the attention was turned to the decodercomplexity issue and including this in the design criteria, several full-rate STCs were proposed as alternatives to Matrix C. In this paper, we review these different alternatives and compare them to Matrix C in terms of performances and the correspondingreceiver complexities.

Silver space-time trellis-coded modulation

Silver Code (SilC) was originally discovered in [1–4] for 2×2 multiple-input multiple-output (MIMO) transmission. It has non-vanishing minimum determinant 1/7, slightly lower than Golden code, but is fast-decodable, i.e., it allows reduced-complexity maximum likelihood decoding [5–7]. In this paper, we present a multidimensional trellis-coded modulation scheme for MIMO systems [11] based on set partitioning of the Silver Code, named Silver Space-Time Trellis Coded Modulation (SST-TCM). This lattice set partitioning is designed specifically to increase the minimum determinant. The branches of the outer trellis code are labeled with these partitions. Viterbi algorithm is applied for trellis decoding, while the branch metrics are computed by using a sphere-decoding algorithm. It is shown that the proposed SST-TCM performs very closely to the Golden Space-Time Trellis Coded Modulation (GST-TCM) scheme, yetwith a much reduced decoding complexity thanks to its fast-decoding property.

A new hybrid summarizer based on vector Space model, statistical physics and linguistics

In this article we present a hybrid approach for automatic summarization of Spanish medical texts. There are a lot of systems for automatic summarization using statistics or linguistics, but only a few of them combining both techniques. Our idea is that to reach a good summary we need to use linguistic aspects of texts, but as well we should benefit of the advantages of statistical techniques. We have integrated the Cortex (Vector Space Model) and Enertex (statistical physics) systems coupled with the Yate term extractor, and the Disicosum system (linguistics). We have compared these systems and afterwards we have integrated them in a hybrid approach. Finally, we have applied this hybrid system over a corpora of medical articles and we have evaluated their performances obtaining good results.

The contribution of patch topology and demographic parameters to PVA predictions: the case of the European tree frog

Population viability analyses (PVA) are increasingly used in metapopulation conservation plans. Two major types of models are commonly used to assess vulnerability and to rank management options: population-based stochastic simulation models (PSM such as RAMAS or VORTEX) and stochastic patch occupancy models (SPOM). While the first set of models relies on explicit intrapatch dynamics and interpatch dispersal to predict population levels in space and time, the latter is based on spatially explicit metapopulation theory where the probability of patch occupation is predicted given the patch area and isolation (patch topology). We applied both approaches to a European tree frog (Hyla arborea) metapopulation in western Switzerland in order to evaluate the concordances of both models and their applications to conservation. Although some quantitative discrepancies appeared in terms of network occupancy and equilibrium population size, the two approaches were largely concordant regarding the ranking of patch values and sensitivities to parameters, which is encouraging given the differences in the underlying paradigms and input data.

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.

Detection of Optimal Models in Parameter Space with Support Vector Machines

The paper proposes an approach aimed at detecting optimal model parameter combinations to achieve the most representative description of uncertainty in the model performance. A classification problem is posed to find the regions of good fitting models according to the values of a cost function. Support Vector Machine (SVM) classification in the parameter space is applied to decide if a forward model simulation is to be computed for a particular generated model. SVM is particularly designed to tackle classification problems in high-dimensional space in a non-parametric and non-linear way. SVM decision boundaries determine the regions that are subject to the largest uncertainty in the cost function classification, and, therefore, provide guidelines for further iterative exploration of the model space. The proposed approach is illustrated by a synthetic example of fluid flow through porous media, which features highly variable response due to the parameter values' combination.

Distance education in neonatal nursing scenarios: a systematic review

OBJECTIVE Identify resources that support learning mediated by technology in the field of neonatal nursing. METHOD Systematic review with searches conducted in MEDLINE, LILACS and SciELO. Titles and abstracts were independently evaluated by two experts. RESULTS Of the 2,051 references, 203 full-text articles were analyzed, resulting in the inclusion of nine studies on semiotics and semiology, cardiopulmonary resuscitation, general aspects of neonatal care, diagnostic reasoning and assessment of pain. Only two articles addressed the development of educational strategies and seven papers described the assessment of these strategies by experts and/or users. CONCLUSION Distance education is an important resource for education, and its improvement and updating, and it particularly adds advantages for neonatal nursing by approximating teaching and real-life situations and by minimizing the exposure of newborns for teaching purposes. The lack of educational initiatives mediated by technology suggests the need for the development, evaluation and dissemination of educational resources focused on nursing care of newborns and their families.