Effect of tibial tray inclination on the femoro-tibial contact pattern of a mobile total knee arthroplasty

Introduction: The posterior inclination of the tibial component is an important factor that can affect the success of total knee arthroplasty. It can reduce the posterior impingement and thus increase the range of flexion, but it may also induce instability in flexion, anterior impingement between the polyethylene of postero-stabilizing knee prosthesis, and anterior conflict with the cortical bone and the stem. Although the problem is identified, there is still a debate on the ideal inclination angle and the surgical technique to avoid an excessive posterior inclination. The aim of this study was to predict the effect of a posterior inclination of the tibial component on the contact pattern on the tibial insert, using a numerical musculoskeletal model of the knee joint. Methods: A 3D finite element model of the knee joint was developed to simulate an active and loaded squat movement after total knee arthroplasty. Flexion was actively controlled by the quadriceps muscle and muscle activations were estimated from EMG data and were synchronized by a feedback algorithm. Two inclinations of the tibial tray were considered: a posterior inclination of 0° or 10°. During the entire range of flexion, the following quantities were calculated: the tibiofemoral and patello-femoral contact force, and the contact pattern on polyethylene insert. The antero-posterior displacement of the contact pattern was also measured. Abaqus 6.7 was used for all analyses. Results: The tibio-femoral and patello-femoral contact forces increased during flexion and reached respectively 4 and 7 BW (bodyweight) at 90° of flexion. They were slightly affected by the inclination of the tibial tray. Without posterior inclination, the contact pattern on the tibial insert remained centered. The contact pressure was lower than 5 MPa below 60° of flexion, but exceeded 20 MPa at 90° of flexion. The posterior inclination displaced the contact point posteriorly by 2 to 4 mm. Conclusion: The inclination of the tibial tray displaced the contactpattern towards the posterior border of the tibial insert. However, even for 10° of inclination, the contact center remained far from the posterior border (12 mm). There was no instability predicted for this movement.

Anatomic shoulder arthroplasty: consequences of humeral head malposition

Purpose Third generation anatomic total shoulder prostheses offer a wide range of adaptability (size, thickness, retroversion and offset of the humeral head, cervico-diaphyseal angle) in order to reproduce anatomy and biomechanics of the shoulder as normal as possible. The large variability of the implants may also induce malposition. Our goal was to analyse the consequences of a humeral head malposition, which is one of the most frequent placement errors. Material and Methods A 3D finite element model of the glenohumeral joint, including the rotator cuff muscles and the deltoid, was used with the Aequalis anatomic prosthesis. Active abduction was simulated. Three humeral head placements were compared : anatomic positioning (A), 5 mm inferior positioning (B), 5 mm superior positioning (C). The effect of humeral head malposition was evaluated through the following quantities : the range of motion free of impingements, the glenohumeral contact pattern, and the stress within the polyethylene and the cement. Results Inferior positioning (B) of the humeral head produced a superior impingement before 90° of abduction, an inferior eccentric contact point on the glenoid, and 165% increase of cement stress. Superior positioning (C) of the humeral head produced a postero-superior eccentric contact point on the glenoid, 300% increase of glenohumeral contact pressure, 450% increase of polyethylene stress, and 207% increase of cement stress. Conclusion Malposition of the humeral head of anatomic prostheses induces biomechanical consequences that may preclude the glenoid survival. Particular attention must be paid to reproduce the humeral anatomy as normal as possible.

Comparison of polyethylene wear in anatomical and reversed shoulder prostheses.

Wear of polyethylene is associated with aseptic loosening of orthopaedic implants and has been observed in hip and knee prostheses and anatomical implants for the shoulder. The reversed shoulder prostheses have not been assessed as yet. We investigated the volumetric polyethylene wear of the reversed and anatomical Aequalis shoulder prostheses using a mathematical musculoskeletal model. Movement and joint stability were achieved by EMG-controlled activation of the muscles. A non-constant wear factor was considered. Simulated activities of daily living were estimated from in vivo recorded data. After one year of use, the volumetric wear was 8.4 mm(3) for the anatomical prosthesis, but 44.6 mm(3) for the reversed version. For the anatomical prosthesis the predictions for contact pressure and wear were consistent with biomechanical and clinical data. The abrasive wear of the polyethylene in reversed prostheses should not be underestimated, and further analysis, both experimental and clinical, is required.

Effect of a pathological scapular tilt after total shoulder arthroplasty.

Total shoulder arthroplasty (TSA) is an accepted and most successfully used treatment for different shoulder pathologies. Different risk factors for the failure of the prosthesis are known. A pathological scapular orientation, observed in elderly people or in patients suffering from neuromuscular diseases, could be a cause of failure, which has not been investigated yet. To test this hypothesis, a numerical musculoskeletal model of the glenohumeral joint was used to compare two TSA cases: a reference normal case and a case with a pathological anterior tilt of the scapula. An active abduction of 150° was simulated. Joint force, contact pattern, polyethylene and cement stress were evaluated for both cases. The pathological tilt slightly increased the joint force and the contact pressure, but also shifted the contact pattern. This eccentric contact increased the stress level within the polyethylene of the glenoid component and within the surrounding cement layer. This adverse effect occurred mainly during the first 60° of abduction. Therefore, a pathological orientation of the scapula may increase the risk of a failure of the cement layer around the glenoid component. These preliminary numerical results should be confirmed by a clinical study.

Effect of the Scapulo-Humeral Rhythm on Anatomical and Reverse Shoulder Prostheses

Introduction: Several studies have reported significant alteration of the scapula-humeral rythm after total shoulder arthroplasty. However, the biomechanical and clinical effects, particularly on implants lifespan, are still unknown. The goal of this study was to evaluate the biomechanical consequences of an altered scapula-humeral rhythm. Methods: A numerical musculoskeletal model of the shoulder was used. The model included the scapula, the humerus and 6 scapulohumeral muscles: middle, anterior, and posterior deltoid, supraspinatus, subscapularis and infraspinatus combined with teres minor. Arm motion and joint stability were achieved by muscles. The reverse and anatomic Aequalis prostheses (Tornier Inc) were inserted. Two scapula-humeral rhythms were considered for each prosthesis: a normal 2:1 rhythm, and an altered 1:2 rhythm. For the 4 configurations, a movement of abduction in the scapular plane was simulated. The gleno-humeral force and contact pattern, but also the stress in the polyethylene and cement were evaluated. Results: With the anatomical prosthesis, the gleno-humeral force increased of 23% for the altered rhythm, with a more eccentric (posterior and superior) contact. The contact pressure, polyethylene stress, and cement stress increased respectively by 20%, 48% and 64%. With the reverse prosthesis, the gleno-humeral force increased of 11% for an altered rhythm. There was nearly no effect on the contact pattern on the polyethylene component surface. Conclusion: The present study showed that alteration oft the scapula-humeral rythm induced biomechanical consequences which could preclude the long term survival of the glenoid implant of anatomic prostheses. However,an altered scapula-humeral rhythm, even severe, should not be a contra indication for the use of a reverse prosthesis. 

Polyethylene wear in anatomic and reversed shoulder prostheses

Purpose: Polyethylene wear is a recurrent problem in joint arthroplasty. Small debris particles are also associated to inflammation reaction of the surrounding bone, eventually leading to the failure of the bound between the implant and the host bone, and implant loosening. The goal of this study was thus to estimate the volume of polyethylene wear of a reversed prosthesis, and compare it to an anatomic prosthesis, during one year of activities of daily living. Material and Methods: A numerical musculoskeletal model of the glenohumeral joint was used for this comparative study. The reversed (RP) and anatomic (AP) Aequalis prostheses were positioned in the numerical model. Eight levels of abduction were considered. Their daily frequency was estimated from in-vivo recorded data on healthy volunteers during activities of daily living. One year of use was simulated to predict the linear and volumetric wear. The volumetric wear was the difference of volume between the original and worn component. Results: With the AP, the contact pattern on the glenoid surface moved rapidly from the inferior to the superior side during the first 30 degrees of abduction, and then went back to the inferior side. With the RP, the contact pattern on the humeral cup surface remained at the inferior side. Contact pressure was 20 times lower with the RP than with the AP. One year of use produced a maximum linear wear of 0.2 mm with the AP, and 0.13 mm with the RP. However, the volumetric wear was 8.4 mm3 with the AP, but reached 44.6 mm3 with the RP. Conclusion: Polyethylene particles are a matter of concern with AP. Infiltration of these particles within the bone-implant interface can induce a implant loosening. This problem should not be underestimated with RP. It might be associated to a higher level of humeral stem loosening reported with RP. The long term survival of RP might be improved by using a highly cross-linked polyethylene, which has a better abrasion resistance but lower plastic resistance.

Comparison of an EMG-based and a stress-based method to predict shoulder muscle forces.

The estimation of muscle forces in musculoskeletal shoulder models is still controversial. Two different methods are widely used to solve the indeterminacy of the system: electromyography (EMG)-based methods and stress-based methods. The goal of this work was to evaluate the influence of these two methods on the prediction of muscle forces, glenohumeral load and joint stability after total shoulder arthroplasty. An EMG-based and a stress-based method were implemented into the same musculoskeletal shoulder model. The model replicated the glenohumeral joint after total shoulder arthroplasty. It contained the scapula, the humerus, the joint prosthesis, the rotator cuff muscles supraspinatus, subscapularis and infraspinatus and the middle, anterior and posterior deltoid muscles. A movement of abduction was simulated in the plane of the scapula. The EMG-based method replicated muscular activity of experimentally measured EMG. The stress-based method minimised a cost function based on muscle stresses. We compared muscle forces, joint reaction force, articular contact pressure and translation of the humeral head. The stress-based method predicted a lower force of the rotator cuff muscles. This was partly counter-balanced by a higher force of the middle part of the deltoid muscle. As a consequence, the stress-based method predicted a lower joint load (16% reduced) and a higher superior-inferior translation of the humeral head (increased by 1.2 mm). The EMG-based method has the advantage of replicating the observed cocontraction of stabilising muscles of the rotator cuff. This method is, however, limited to available EMG measurements. The stress-based method has thus an advantage of flexibility, but may overestimate glenohumeral subluxation.

Low-pressure, water-assisted anatexis of basic dykes in a contact metamorphic aureole, Fuerteventura (Canary Islands): oxygen isotope evidence for a meteoric fluid origin

Migmatites produced by low-pressure anatexis of basic dykes are found in a contact metamorphic aureole around a pyroxenite-gabbro intrusion (PX2), on Fuerteventura. Dykes outside and inside the aureole record interaction with meteoric water, with low or negative delta O-18 whole-rock values (+0.2 to -3.4 parts per thousand), decreasing towards the contact. Recrystallised plagioclase, diopside, biotite and oxides, from within the aureole, show a similar evolution with lowest delta O-18 values (-2.8, -4.2, - 4.4 and -7.6 parts per thousand, respectively) in the migmatite zone, close to the intrusion. Relict clinopyroxene phenocrysts preserved in all dykes, retain typically magmatic delta O-18 values up to the anatectic zone, where the values are lower and more heterogeneous. Low delta O-18 values, decreasing towards the intrusion, can be ascribed 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.

Assessing the Functionality and Biocompatibility of a Wireless Contact Lens Sensor for 24 hour-intraocular Pressure Monitoring in Volunteers: Preliminary Results

Purpose: In vitro studies in porcine eyes have demonstrated a good correlation between induced intraocular pressure variations and corneal curvature changes, using a contact lens with an embedded microfabricated strain gauge. Continuous 24 hour-intraocular pressure (IOP) monitoring to detect large diurnal fluctuation is currently an unmet clinical need. The aims of this study is to evaluate precision of signal transmission and biocompatibility of 24 hour contact lens sensor wear (SENSIMED Triggerfish®) in humans. Methods: After full eye examination in 10 healthy volunteers, a 8.7 mm radius contact lens sensor and an orbital bandage containing a loop antenna were applied and connected to a portable recorder. Best corrected visual acuity and position, lubrication status and mobility of the sensor were assessed after 5 and 30 minutes, 4, 7 and 24 hours. Subjective comfort was scored and activities documented in a logbook. After sensor removal full eye examination was repeated, and the registration signal studied. Results: The comfort score was high and did not fluctuate significantly, except at the 7 hour-visit. The mobility of the contact lens was minimal but its lubrication remained good. Best corrected visual acuity was significantly reduced during the sensor wear and immediately after its removal. Three patients developed mild corneal staining. In all but one participant we obtained a registration IOP curve with visible ocular pulse amplitude. Conclusions: This 24 hour-trial confirmed the functionality and biocompatibility of SENSIMED Triggerfish® wireless contact lens sensor for IOP-fluctuation monitoring in volunteers. Further studies with a range of different contact lens sensor radii are indicated.

Shallow-level migmatization of gabbros in a metamorphic contact aureole, Fuerteventura Basal Complex, Canary Islands

Migmatization of gabbroic rocks at 2-3 kbar has occurred in the metamorphic contact aureole of a mafic pluton in the Fuerteventura Basal Complex (Canary Island;). Migmatites are characterized by a dense network: of closely spaced millimetre-wide leucocratic veins with perfectly preserved igneous textures. They are all relatively enriched in Al, Na I: Sr Ba, Nb, Y and the rare earth elements compared with the unaffected country rock beyond the aureole. Migmatization under such low-pressure conditions war possible because of the unusual tectonic and magmatic contact in which ii occurred. Multiple basic intrusions associated with extrusive volcanic activity created high heat flow in a small area. Alkaline and metasomatized rocks present in the country rock of the intruding pluton were leached by high-temperature fluids during contact metamorphism. These enriched fluids then favoured partial melting of the host gabbroic rocks, and contaminated both the leucosomes and melanosomes. A transpressive tectonic setting at the time of intrusion created shearing along the contact between the intrusion and its host rock. This shearing enhanced circulation of the fluids and allowed segregation of the nea-formed melts from their restite by opening tension veins into which the melts migrated. Depending on the relative timing of melt segregation and recrystallization leucosomes range in composition from a 40-60% mixture of clinopyroxene (+/- amphibole) and plagioclase to almost pure feldspathic veins. Comparable occurrences of gabbros migmatized at low pressure are expected only at a snail scale in localized areas of high heat flow in the presence of fluids, such as in. mid-ocean ridges or ocean-islands.

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

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

Devolatilization-induced pressure build-up: Implications for reaction front movement and breccia pipe formation

Generation of fluids during metamorphism can significantly influence the fluid overpressure, and thus the fluid flow in metamorphic terrains. There is currently a large focus on developing numerical reactive transport models, and with it follows the need for analytical solutions to ensure correct numerical implementation. In this study, we derive both analytical and numerical solutions to reaction-induced fluid overpressure, coupled to temperature and fluid flow out of the reacting front. All equations are derived from basic principles of conservation of mass, energy and momentum. We focus on contact metamorphism, where devolatilization reactions are particularly important owing to high thermal fluxes allowing large volumes of fluids to be rapidly generated. The analytical solutions reveal three key factors involved in the pressure build-up: (i) The efficiency of the devolatilizing reaction front (pressure build-up) relative to fluid flow (pressure relaxation), (ii) the reaction temperature relative to the available heat in the system and (iii) the feedback of overpressure on the reaction temperature as a function of the Clapeyron slope. Finally, we apply the model to two geological case scenarios. In the first case, we investigate the influence of fluid overpressure on the movement of the reaction front and show that it can slow down significantly and may even be terminated owing to increased effective reaction temperature. In the second case, the model is applied to constrain the conditions for fracturing and inferred breccia pipe formation in organic-rich shales owing to methane generation in the contact aureole.

Enregistrement des fluctuations de la pression intraoculaire sur 24 heures par une lentille de contact intégrant un capteur: tolérance et résultats chez 10 sujets sains

AbstractPurpose: to evaluate the tolerability, comfort and precision of the signal transmission of an ocular Sensor used for 24-hour intraocular pressure fluctuation monitoring in humans.Patients and methods: In this uncontrolled open trial involving 10 healthy volunteers an 8.7 mm radius prototype ocular telemetry Sensor (SENSIMED Triggerfish®, Lausanne, Switzerland) and an orbital bandage containing a loop antenna were applied and connected to a portable recorder after full eye examination. Best corrected visual acuity and position, surface wetting ability and mobility of the Sensor were assessed after 5 and 30 minutes, 4, 12 and 24 hours. Subjective wearing comfort was scored and activities documented in a logbook. After Sensor removal a full eye examination was repeated and the recorded signal analyzed.Results: The comfort score was high and did not fluctuate significantly over time. The mobility of the Sensor was limited across follow-up visits and its surface wetting ability remained good. Best corrected visual acuity was significantly reduced during Sensor wear and immediately after its removal (from 1.07 before, to 0.85 after, P-value 0.008). Three subjects developed a mild, transient corneal abrasion. In all but one participant we obtained usable data of a telemetric signal recording with sufficient sensitivity to depict ocular pulsation.Conclusions: This 24-hour- trial has encouraging results on the tolerability and functionality of the ocular telemetric Sensor for intraocular pressure fluctuation monitoring. Further studies with different Sensor radii conducted on a larger study population are needed to improve comfort, precision and interpretation of the telemetric signal.

Nucleation and growth of metamorphic minerals in contact aureoles

THESIS ABSTRACT Nucleation and growth of metamorphic minerals are the consequence of changing P-T-X-conditions. The thesis presented here focuses on processes governing nucleation and growth of minerals in contact metamorphic environments using a combination of geochemical analytics (chemical-, isotope-, and trace element composition), statistical treatments of spatial data, and numerical models. It is shown, that a combination of textural modeling and stable isotope analysis allows a distinction between several possible reaction paths for olivine growth in a siliceous dolomite contact aureole. It is suggested that olivine forms directly from dolomite and quartz. The formation of olivine from this metastable reaction implies metamorphic crystallization far from equilibrium. As a major consequence, the spatial distribution of metamorphic mineral assemblages in a contact aureole cannot be interpreted as a proxy for the temporal evolution of a single rock specimen, because each rock undergoes a different reaction path, depending on temperature, heating rate, and fluid-infiltration rate. A detailed calcite-dolomite thermometry study was initiated on multiple scales ranging from aureole scale to the size of individual crystals. Quantitative forward models were developed to evaluate the effect of growth zoning, volume diffusion and the formation of submicroscopic exsolution lamellae (<1 µm) on the measured Mg-distribution in individual calcite crystals and compare the modeling results to field data. This study concludes that Mg-distributions in calcite grains of the Ubehebe Peak contact aureole are the consequence of rapid crystal growth in combination with diffusion and exsolution. The crystallization history of a rock is recorded in the chemical composition, the size and the distribution of its minerals. Near the Cima Uzza summit, located in the southern Adamello massif (Italy), contact metamorphic brucite bearing dolomite marbles are exposed as xenoliths surrounded by mafic intrusive rocks. Brucite is formed retrograde pseudomorphing spherical periclase crystals. Crystal size distributions (CSD's) of brucite pseudomorphs are presented for two profiles and combined with geochemistry data and petrological information. Textural analyses are combined with geochemistry data in a qualitative model that describes the formation periclase. As a major outcome, this expands the potential use of CSD's to systems of mineral formation driven by fluid-infiltration. RESUME DE LA THESE La nucléation et la croissance des minéraux métamorphiques sont la conséquence de changements des conditions de pression, température et composition chimique du système (PT-X). Cette thèse s'intéresse aux processus gouvernant la nucléation et la croissance des minéraux au cours d'un épisode de métamorphisme de contact, en utilisant la géochimie analytique (composition chimique, isotopique et en éléments traces), le traitement statistique des données spatiales et la modélisation numérique. Il est montré que la combinaison d'un modèle textural avec des analyses en isotopes stables permet de distinguer plusieurs chemins de réactions possibles conduisant à la croissance de l'olivine dans une auréole de contact riche en Silice et dolomite. Il est suggéré que l'olivine se forme directement à partir de la dolomie et du quartz. Cette réaction métastable de formation de l'olivine implique une cristallisation métamorphique loin de l'équilibre. La principale conséquence est que la distribution spatiale des assemblages de minéraux métamorphiques dans une auréole de contact ne peut pas être considérée comme un témoin de l'évolution temporelle d'un type de roche donné, puisque chaque type de roche suit différents chemins de réactions, en fonction de la température, la vitesse de réchauffement et le taux d'infiltration du fluide. Une étude thermométrique calcite-dolomite détaillée a été réalisée à diverses échelles, depuis l'échelle de l'auréole de contact jusqu'à l'échelle du cristal. Des modèles numériques quantitatifs ont été développés pour évaluer l'effet des zonations de croissance, de la diffusion volumique et de la formation de lamelles d'exsolution submicroscopiques (<1µm) sur la distribution du magnésium mesuré dans des cristaux de calcite individuels. Les résultats de ce modèle ont été comparés ä des échantillons naturels. Cette étude montre que la distribution du Mg dans les grains de calcite de l'auréole de contact de l'Ubehebe Peak (USA) résulte d'une croissance cristalline rapide, associée aux processus de diffusion et d'exsolution. L'histoire de cristallisation d'une roche est enregistrée dans la composition chimique, la taille et la distribution de ses minéraux. Près du sommet Cima Uzza situé au sud du massif d'Adamello (Italie), des marbres dolomitiques à brucite du métamorphisme de contact forment des xénolithes dans une intrusion mafique. La brucite constitue des pseudomorphes rétrogrades du périclase. Les distributions de taille des cristaux (CSD) des pseudomorphes de brucite sont présentées pour deux profiles et sont combinées aux données géochimiques et pétrologiques. Les analyses textorales sont combinées aux données géochimiques dans un modèle qualitatif qui décrit la formation du périclase. Ceci élargit l'utilisation potentielle de la C5D aux systèmes de formation de minéraux controlés par les infiltrations fluides. THESIS ABSTRACT (GENERAL PUBLIC) Rock textures are essentially the result of a complex interaction of nucleation, growth and deformation as a function of changing physical conditions such as pressure and temperature. Igneous and metamorphic textures are especially attractive to study the different mechanisms of texture formation since most of the parameters like pressure-temperature-paths are quite well known for a variety of geological settings. The fact that textures are supposed to record the crystallization history of a rock traditionally allowed them to be used for geothermobarometry or dating. During the last decades the focus of metamorphic petrology changed from a static point of view, i.e. the representation of a texture as one single point in the petrogenetic grid towards a more dynamic view, where multiple metamorphic processes govern the texture formation, including non-equilibrium processes. This thesis tries to advance our understanding on the processes governing nucleation and growth of minerals in contact metamorphic environments and their dynamic interplay by using a combination of geochemical analyses (chemical-, isotope-, and trace element composition), statistical treatments of spatial data and numerical models. In a first part the thesis describes the formation of metamorphic olivine porphyroblast in the Ubehebe Peak contact aureole (USA). It is shown that not the commonly assumed succession of equilibrium reactions along a T-t-path formed the textures present in the rocks today, but rather the presence of a meta-stable reaction is responsible for forming the olivine porphyroblast. Consequently, the spatial distribution of metamorphic minerals within a contact aureole can no longer be regarded as a proxy for the temporal evolution of a single rock sample. Metamorphic peak temperatures for samples of the Ubehebe Peak contact aureole were determined using calcite-dolomite. This geothermometer is based on the temperature-dependent exchange of Mg between calcite and dolomite. The purpose of the second part of this thesis was to explain the interfering systematic scatter of measured Mg-content on different scales and thus to clarify the interpretation of metamorphic temperatures recorded in carbonates. Numerical quantitative forward models are used to evaluate the effect of several processes on the distribution of magnesium in individual calcite crystals and the modeling results were then compared to measured field. Information about the crystallization history is not only recorded in the chemical composition of grains, like isotope composition or mineral zoning. Crystal size distributions (CSD's) provide essential information about the complex interaction of nucleation and growth of minerals. CSD's of brucite pseudomorphs formed retrograde after periclase of the southern Adamello massif (Italy) are presented. A combination of the textural 3D-information with geochemistry data is then used to evaluate reaction kinetics and to constrain the actual reaction mechanism for the formation of periclase. The reaction is shown to be the consequence of the infiltration of a limited amount of a fluid phase at high temperatures. The composition of this fluid phase is in large disequilibrium with the rest of the rock resulting in very fast reaction rates. RESUME DE LA THESE POUR LE GRAND PUBLIC: La texture d'une roche résulte de l'interaction complexe entre les processus de nucléation, croissance et déformation, en fonction des variations de conditions physiques telles que la pression et la température. Les textures ignées et métamorphiques présentent un intérêt particulier pour l'étude des différents mécanismes à l'origine de ces textures, puisque la plupart des paramètres comme les chemin pression-température sont relativement bien contraints dans la plupart des environnements géologiques. Le fait que les textures soient supposées enregistrer l'histoire de cristallisation des roches permet leur utilisation pour la datation et la géothermobarométrie. Durant les dernières décennies, la recherche en pétrologie métamorphique a évolué depuis une visualisation statique, c'est-à-dire qu'une texture donnée correspondait à un point unique de la grille pétrogénétique, jusqu'à une visualisation plus dynamique, où les multiples processus métamorphiques qui gouvernent 1a formation d'une texture incluent des processus hors équilibre. Cette thèse a pour but d'améliorer les connaissances actuelles sur les processus gouvernant la nucléation et la croissance des minéraux lors d'épisodes de métamorphisme de contact et l'interaction dynamique existant entre nucléation et croissance. Pour cela, les analyses géochimiques (compositions chimiques en éléments majeurs et traces et composition isotopique), le traitement statistique des données spatiales et la modélisation numérique ont été combinés. Dans la première partie, cette thèse décrit la formation de porphyroblastes d'olivine métamorphique dans l'auréole de contact de l'Ubehebe Peak (USA). Il est montré que la succession généralement admise des réactions d'équilibre le long d'un chemin T-t ne peut pas expliquer les textures présentes dans les roches aujourd'hui. Cette thèse montre qu'il s'agirait plutôt d'une réaction métastable qui soit responsable de la formation des porphyroblastes d'olivine. En conséquence, la distribution spatiale des minéraux métamorphiques dans l'auréole de contact ne peut plus être interprétée comme le témoin de l'évolution temporelle d'un échantillon unique de roche. Les pics de température des échantillons de l'auréole de contact de l'Ubehebe Peak ont été déterminés grâce au géothermomètre calcite-dolomite. Celui-ci est basé sur l'échange du magnésium entre la calcite et la dolomite, qui est fonction de la température. Le but de la deuxième partie de cette thèse est d'expliquer la dispersion systématique de la composition en magnésium à différentes échelles, et ainsi d'améliorer l'interprétation des températures du métamorphisme enregistrées dans les carbonates. Des modèles numériques quantitatifs ont permis d'évaluer le rôle de différents processus sur la distribution du magnésium dans des cristaux de calcite individuels. Les résultats des modèles ont été comparés aux échantillons naturels. La composition chimique des grains, comme la composition isotopique ou la zonation minérale, n'est pas le seul témoin de l'histoire de la cristallisation. La distribution de la taille des cristaux (CSD) fournit des informations essentielles sur les interactions entre nucléation et croissance des minéraux. La CSD des pseudomorphes de brucite retrograde formés après le périclase dans le sud du massif Adamello (Italie) est présentée dans la troisième partie. La combinaison entre les données textorales en trois dimensions et les données géochimiques a permis d'évaluer les cinétiques de réaction et de contraindre les mécanismes conduisant à la formation du périclase. Cette réaction est présentée comme étant la conséquence de l'infiltration d'une quantité limitée d'une phase fluide à haute température. La composition de cette phase fluide est en grand déséquilibre avec le reste de la roche, ce qui permet des cinétiques de réactions très rapides.

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.