Pore-scale modeling of two-phase flow instabilities in porous media

Les problèmes d'écoulements multiphasiques en média poreux sont d'un grand intérêt pour de nombreuses applications scientifiques et techniques ; comme la séquestration de C02, l'extraction de pétrole et la dépollution des aquifères. La complexité intrinsèque des systèmes multiphasiques et l'hétérogénéité des formations géologiques sur des échelles multiples représentent un challenge majeur pour comprendre et modéliser les déplacements immiscibles dans les milieux poreux. Les descriptions à l'échelle supérieure basées sur la généralisation de l'équation de Darcy sont largement utilisées, mais ces méthodes sont sujettes à limitations pour les écoulements présentant de l'hystérèse. Les avancées récentes en terme de performances computationnelles et le développement de méthodes précises pour caractériser l'espace interstitiel ainsi que la distribution des phases ont favorisé l'utilisation de modèles qui permettent une résolution fine à l'échelle du pore. Ces modèles offrent un aperçu des caractéristiques de l'écoulement qui ne peuvent pas être facilement observées en laboratoire et peuvent être utilisé pour expliquer la différence entre les processus physiques et les modèles à l'échelle macroscopique existants. L'objet premier de la thèse se porte sur la simulation numérique directe : les équations de Navier-Stokes sont résolues dans l'espace interstitiel et la méthode du volume de fluide (VOF) est employée pour suivre l'évolution de l'interface. Dans VOF, la distribution des phases est décrite par une fonction fluide pour l'ensemble du domaine et des conditions aux bords particulières permettent la prise en compte des propriétés de mouillage du milieu poreux. Dans la première partie de la thèse, nous simulons le drainage dans une cellule Hele-Shaw 2D avec des obstacles cylindriques. Nous montrons que l'approche proposée est applicable même pour des ratios de densité et de viscosité très importants et permet de modéliser la transition entre déplacement stable et digitation visqueuse. Nous intéressons ensuite à l'interprétation de la pression capillaire à l'échelle macroscopique. Nous montrons que les techniques basées sur la moyenne spatiale de la pression présentent plusieurs limitations et sont imprécises en présence d'effets visqueux et de piégeage. Au contraire, une définition basée sur l'énergie permet de séparer les contributions capillaires des effets visqueux. La seconde partie de la thèse est consacrée à l'investigation des effets d'inertie associés aux reconfigurations irréversibles du ménisque causé par l'interface des instabilités. Comme prototype pour ces phénomènes, nous étudions d'abord la dynamique d'un ménisque dans un pore angulaire. Nous montrons que, dans un réseau de pores cubiques, les sauts et reconfigurations sont si fréquents que les effets d'inertie mènent à différentes configurations des fluides. A cause de la non-linéarité du problème, la distribution des fluides influence le travail des forces de pression, qui, à son tour, provoque une chute de pression dans la loi de Darcy. Cela suggère que ces phénomènes devraient être pris en compte lorsque que l'on décrit l'écoulement multiphasique en média poreux à l'échelle macroscopique. La dernière partie de la thèse s'attache à démontrer la validité de notre approche par une comparaison avec des expériences en laboratoire : un drainage instable dans un milieu poreux quasi 2D (une cellule Hele-Shaw avec des obstacles cylindriques). Plusieurs simulations sont tournées sous différentes conditions aux bords et en utilisant différents modèles (modèle intégré 2D et modèle 3D) afin de comparer certaines quantités macroscopiques avec les observations au laboratoire correspondantes. Malgré le challenge de modéliser des déplacements instables, où, par définition, de petites perturbations peuvent grandir sans fin, notre approche numérique apporte de résultats satisfaisants pour tous les cas étudiés. - Problems involving multiphase flow in porous media are of great interest in many scientific and engineering applications including Carbon Capture and Storage, oil recovery and groundwater remediation. The intrinsic complexity of multiphase systems and the multi scale heterogeneity of geological formations represent the major challenges to understand and model immiscible displacement in porous media. Upscaled descriptions based on generalization of Darcy's law are widely used, but they are subject to several limitations for flow that exhibit hysteric and history- dependent behaviors. Recent advances in high performance computing and the development of accurate methods to characterize pore space and phase distribution have fostered the use of models that allow sub-pore resolution. These models provide an insight on flow characteristics that cannot be easily achieved by laboratory experiments and can be used to explain the gap between physical processes and existing macro-scale models. We focus on direct numerical simulations: we solve the Navier-Stokes equations for mass and momentum conservation in the pore space and employ the Volume Of Fluid (VOF) method to track the evolution of the interface. In the VOF the distribution of the phases is described by a fluid function (whole-domain formulation) and special boundary conditions account for the wetting properties of the porous medium. In the first part of this thesis we simulate drainage in a 2-D Hele-Shaw cell filled with cylindrical obstacles. We show that the proposed approach can handle very large density and viscosity ratios and it is able to model the transition from stable displacement to viscous fingering. We then focus on the interpretation of the macroscopic capillary pressure showing that pressure average techniques are subject to several limitations and they are not accurate in presence of viscous effects and trapping. On the contrary an energy-based definition allows separating viscous and capillary contributions. In the second part of the thesis we investigate inertia effects associated with abrupt and irreversible reconfigurations of the menisci caused by interface instabilities. As a prototype of these phenomena we first consider the dynamics of a meniscus in an angular pore. We show that in a network of cubic pores, jumps and reconfigurations are so frequent that inertia effects lead to different fluid configurations. Due to the non-linearity of the problem, the distribution of the fluids influences the work done by pressure forces, which is in turn related to the pressure drop in Darcy's law. This suggests that these phenomena should be taken into account when upscaling multiphase flow in porous media. The last part of the thesis is devoted to proving the accuracy of the numerical approach by validation with experiments of unstable primary drainage in a quasi-2D porous medium (i.e., Hele-Shaw cell filled with cylindrical obstacles). We perform simulations under different boundary conditions and using different models (2-D integrated and full 3-D) and we compare several macroscopic quantities with the corresponding experiment. Despite the intrinsic challenges of modeling unstable displacement, where by definition small perturbations can grow without bounds, the numerical method gives satisfactory results for all the cases studied.

A reassessment of models for hydrocarbon generation in the Khibiny nepheline syenite complex, Kola Peninsula, Russia

Although hydrocarbon-bearing fluids have been known from the alkaline igneous rocks of the Khibiny intrusion for many years, their origin remains enigmatic. A recently proposed model of post-magmatic hydrocarbon (HC) generation through Fischer-Tropsch (FT) type reactions suggests the hydration of Fe-bearing phases and release of H-2 which reacts with magmatically derived CO2 to form CH4 and higher HCs. However, new petrographic, microthermometric, laser Raman, bulk gas and isotope data are presented and discussed in the context of previously published work in order to reassess models of HC generation. The gas phase is dominated by CH4 with only minor proportions of higher hydrocarbons. No remnants of the proposed primary CO2-rich fluid are found in the complex. The majority of the fluid inclusions are of secondary nature and trapped in healed microfractures. This indicates a high fluid flux after magma crystallisation. Entrapment conditions for fluid inclusions are 450-550 degrees C at 2.8-4.5 kbar. These temperatures are too high for hydrocarbon gas generation through the FT reaction. Chemical analyses of rims of Fe-rich phases suggest that they are not the result of alteration but instead represent changes in magma composition during crystallisation. Furthermore, there is no clear relationship between the presence of Fe-rich minerals and the abundance of fluid inclusion planes (FIPs) as reported elsewhere. delta C-13 values for methane range from -22.4% to -5.4%, confirming a largely abiogenic origin for the gas. The presence of primary CH4-dominated fluid inclusions and melt inclusions, which contain a methane-rich gas phase, indicates a magmatic origin of the HCs. An increase in methane content, together with a decrease in delta C-13 isotope values towards the intrusion margin suggests that magmatically derived abiogenic hydrocarbons may have mixed with biogenic hydrocarbons derived from the surrounding country rocks. (C) 2006 Elsevier BV. All rights reserved.

Exposure to bioaerosols in poultry houses at different stages of fattening: use of real-time PCR for airborne bacterial quantification

Previous studies have demonstrated that poultry-house workers are exposed to very high levels of organic dust and consequently have an increased prevalence of adverse respiratory symptoms. However, the influence of the age of broilers, on bioaerosol concentrations has not been investigated. To evaluate the evolution of bioaerosol concentration during the fattening period, bioaerosol parameters (inhalable dust, endotoxin and bacteria) were measured in 12 poultry confinement buildings in Switzerland, at 3 different stages of the birds' growth; Samples of air taken from within the breathing zones of individual poultry-house employees as they caught the chickens ready to be transported for slaughter, were also analysed. Quantitative PCR (Q-PCR) was used to assess the quantity of total airborne bacteria and total airborne Staphylococcus species. Bioaerosol levels increased significantly during the fattening period of the chickens. During the task of catching mature birds, the mean inhalable dust concentration for a worker was 31 ± 4.7 mg/m3, and endotoxin concentration was 11'080 ± 3436 UE/m3 air, more than ten-fold higher than the Swiss occupational recommended value (1000 UE/m3). The mean exposure level of bird catchers to total bacteria and Staphylococcus species measured by Q-PCR is also very high, respectively reaching values of 72 (± 11) x107 cells/m3 air and 70 (± 16) x106/m3 air. It was concluded that in the absence of wearing protective breathing apparatus, chicken catchers in Switzerland risk exposure beyond recommended limits for all measured bioaerosol parameters. Moreover, the use of Q-PCR to estimate total and specific numbers of airborne bacteria is a promising tool for evaluating any modifications intended to improve the safety of current working practices.

Clinical effectiveness of direct class II restorations: a meta-analysis.

Purpose: More than five hundred million direct dental restorations are placed each year worldwide. In about 55% of the cases, resin composites or compomers are used, and in 45% amalgam. The longevity of posterior resin restorations is well documented. However, data on resin composites that are placed without enamel/dentin conditioning and resin composites placed with self-etching adhesive systems are missing. Material and Methods: The database SCOPUS was searched for clinical trials on posterior resin composites without restricting the search to the year of publication. The inclusion criteria were: (1) prospective clinical trial with at least 2 years of observation; (2) minimum number of restorations at last recall = 20; (3) report on dropout rate; (4) report of operative technique and materials used; (5) utilization of Ryge or modified Ryge evaluation criteria. For amalgam, only those studies were included that directly compared composite resin restorations with amalgam. For the statistical analysis, a linear mixed model was used with random effects to account for the heterogeneity between the studies. P-values under 0.05 were considered significant. Results: Of the 373 clinical trials, 59 studies met the inclusion criteria. In 70% of the studies, Class II and Class I restorations had been placed. The overall success rate of composite resin restorations was about 90% after 10 years, which was not different from that of amalgam. Restorations with compomers had a significantly lower longevity. The main reason for replacement were bulk fractures and caries adjacent to restorations. Both of these incidents were infrequent in most studies and accounted only for about 6% of all replaced restorations after 10 years. Restorations with macrofilled composites and compomer suffered significantly more loss of anatomical form than restorations with other types of material. Restorations that were placed without enamel acid etching and a dentin bonding agent showed significantly more marginal staining and detectable margins compared to those restorations placed using the enamel-etch or etch-and-rinse technique; restorations with self-etching systems were between the other groups. Restorations with compomer suffered significantly more chippings (repairable fracture) than restorations with other materials, which did not statistically differ among each other. Restorations that were placed with a rubber-dam showed significantly fewer material fractures that needed replacement, and this also had a significant effect on the overall longevity. Conclusion: Restorations with hybrid and microfilled composites that were placed with the enamel-etching technique and rubber-dam showed the best overall performance; the longevity of these restorations was similar to amalgam restorations. Compomer restorations, restorations placed with macrofilled composites, and resin restorations with no-etching or self-etching adhesives demonstrated significant shortcomings and shorter longevity.

Subduction of continental crust in the Western Alps

As a result of recent deep reflection and refraction seismology the crustal structure of the Western Alps is now quite well-defined. However, this raises the question of what is present below the Moho, such as a crustal eclogitic root. This study attempts to estimate the volume of this eclogitic root on the basis of palinspastic reconstructions. Even with a minimum estimate of the crustal material involved in the subduction processes which took place during the Alpine orogeny, a significant eclogitized crustal root must be present down to depths of around 100 km below the Po plain. A maximum estimate suggests that a large part of this root could now be recycled in the asthenosphere.

Spatiotemporal adaptive multiscale multiphysics simulations of two-phase flow

We present a spatiotemporal adaptive multiscale algorithm, which is based on the Multiscale Finite Volume method. The algorithm offers a very efficient framework to deal with multiphysics problems and to couple regions with different spatial resolution. We employ the method to simulate two-phase flow through porous media. At the fine scale, we consider a pore-scale description of the flow based on the Volume Of Fluid method. In order to construct a global problem that describes the coarse-scale behavior, the equations are averaged numerically with respect to auxiliary control volumes, and a Darcy-like coarse-scale model is obtained. The space adaptivity is based on the idea that a fine-scale description is only required in the front region, whereas the resolution can be coarsened elsewhere. Temporal adaptivity relies on the fact that the fine-scale and the coarse-scale problems can be solved with different temporal resolution (longer time steps can be used at the coarse scale). By simulating drainage under unstable flow conditions, we show that the method is able to capture the coarse-scale behavior outside the front region and to reproduce complex fluid patterns in the front region.

Oxidation of methane at the CH4/H2O-(CO2) transition zone in the external part of the Central Alps, Switzerland: Evidence from stable isotope investigations

With the aim of understanding the mechanisms that control the metamorphic transition from the CH4- to the H2O-(CO2)-dominated fluid zone in the Helvetic domain of the Central Alps of Switzerland, fluid inclusions in quartz, illite ``crystallinity'' index, vitrinite reflectance, and the stable isotope compositions of vein and whole rock minerals and fluids trapped in quartz were investigated along four cross-sections. Increasing temperature during prograde metamorphism led to the formation of dry gas by hydrocarbon cracking in the CH4-zone. Fluid immiscibility in the H2O-CH4-(CO2)-NaCl system resulted in cogenetic, CH4- and H2O-dominated fluid inclusions. In the CH4-zone, fluids were trapped at temperatures <= 270 +/- 5 degrees C. The end of the CH4-zone is markedby a sudden increase of CO2 content in the gas phase of fluid inclusions. At temperatures > 270 +/- 5 degrees C, in the H2O-zone, the total amount of volatiles within the fluid decreased below 1 mol% with no immiscibility. This resulted m total homogenization temperatures of H2O-(CO2-CH4)-NaCl inclusions below 180 degrees C. Hydrogen isotope compositions of methane in fluid inclusion have delta D values of less than -100 parts per thousand in the CH4-zone, typical for an origin through cracking of higher hydrocarbons, but where the methane has not equilibrated with the pore water. delta D values of fluid inclusion water are around -40 parts per thousand., in isotopic equilibrium with phyllosilicates of the whole rocks. Within the CH4 to H2O(CO2) transition zone, delta D(H2O) values in fluid inclusions decrease to -130 parts per thousand interpreted to reflect the contribution of deuterium depleted water from methane oxidation. In the H2O-zone, delta D(H2O) values increase again towards an average of -30 parts per thousand which is again consistent with isotopic equilibrium with host-rock phyllosilicates. delta C-13 values of methane in fluid inclusions from the CH4-zone are around -27 parts per thousand in isotopic equilibrium with calcite in veins and whole rocks. The delta C-13(CH4) values decrease to less than -35 parts per thousand at the transition to the H2O-zone and are no longer in equilibrium with the carbonates in the whole rocks. delta C-13 values of CO, are variable but too low to be in equilibrium with the wall rock fluids, compatible with a contribution of CO2 from closed system oxidation of methane. Differences in isotopic composition between host-rock and Alpine fissure carbonate are generally small, suggesting that the amount of CO2 produced by oxidation of methane was small compared to the C-budget in the rocks and local pore fluids were buffered by the wall rocks during precipitation of calcite within the fissures. (c) 2006 Elsevier B.V. All rights reserved.

Growth mechanisms and timing of emplacement of a vertically layered mafic intrusion in the feeder zone of an ocean island volcano (Fuerteventura, Canary islands)

Résumé pour le grand public L'île de Fuerteventura (Canaries) offre l'occasion rare d'observer les racines d'un volcan océanique édifié il y a 25 à 30 millions d'années et complètement érodé. On y voit de nombreux petits plutons de forme et composition variées, témoignant d'autant d'épisodes de l'activité magmatique. L'un de ces plutons, appelé PX1, présente une structure inhabituelle formée d'une alternance de bandes verticales d'épaisseur métrique à hectométrique de roches sombres de composition pyroxénilique ou gabbroïque. Les pyroxénites résultent clairement de l'accumulation de cristaux de pyroxènes et non de la simple solidification d'un magma? Se pose dès lors la question de la nature du processus qui a conduit à l'accumulation verticale de niveaux concentrés en pyroxènes. En effet, les litages pyroxénitiques classiques sont subhorizontaux, car ils résultent de l'accumulation gravitaire des cristaux séparés du magma dont ils cristalli¬sent par sédimentation. Cette étude vise à identifier et comprendre les mécanismes qui ont engendré ce Iitage minéralogique vertical et l'im¬portant volume de ces faciès cumulatifs. Nous nous sommes également intéressés aux conditions de pression et de température régnant au moment de la mise en place du pluton, ainsi qu'à sa durée de vie et à sa vitesse de refroidis¬sement. Enfin une approche géochimique nous a permis de préciser la nature de la source mantellique des magmas liés à cette activité magmatique. PX1 est en réalité un complexe filonien formé à des conditions de pression et de température de 1-2 kbar et 1050- 1100°C; sa construction a nécessité au moins 150 km3 de magma. L'alternance d'horizons gabbroïques et pyroxéniti¬ques représente des injections successives de magma sous la forme de filons verticaux, mis en place dans un contexte régional en extension. L'étude des orientations des minéraux dans ces faciès révèle que les horizons gabbroïques enregistrent l'extension régionale, alors que les pyroxénites sont générées par une compaction au sein du pluton. Ceci suggère que le régime des contraintes, qui était extensif lors de l'initiation de la mise en place de PX1, est pério¬diquement devenu compressif au sein même du pluton. Cette compression serait liée à des cycles de mise en place où la vitesse de croissance du pluton dépassait celle de l'extension régionale. La différenciation observée au sein de chaque horizon, depuis des pyroxénites riches en olivine jusqu'à des pyroxé¬nites à plagioclase interstitiel et des gabbros, ainsi que la composition géochimique des minéraux qui les constituent suggèrent que chaque filon vertical s'est mis en place à partir d'un magma de composition identique, puis a évolué indépendamment des autres en fonction du régime thermique et du régime des contraintes local. Lorsque le magma en train de cristalliser s'est trouvé en compression, le liquide résiduel a été séparé des cristaux déjà formés et extrait du système, laissant derrière lui une accumulation de cristaux dont la nature et les proportions dépendaient du stade de cristallisation atteint par le magma au moment de l'extraction. Ainsi, les niveaux de pyroxénites à olivine (premier minéral à cristalliser) ont été formés lorsque le magma correspondant était encore peu cristallisé; à l'inverse, les py¬roxénites riches en plagioclase (minéral plus tardif dans la séquence de cristallisation) et certains gabbros à caractère cumulatif résultent d'une compression tardive dans le processus de cristallisation du filon concerné. Les liquides résiduels extraits des niveaux pyroxénitiques sont rarement observés dans PX1, certaines poches et filonets de com¬position anorthositique pourraient en être les témoins. L'essentiel de ces liquides a probablement gagné des niveaux supérieurs du pluton, voire la surface du volcan. L'origine du régime compressif périodique affectant les filons en voie de cristallisation est attribuée aux injections suivantes de magma au sein du pluton, qui se sont succédées à un rythme plus rapide que la vitesse de consolidation des filons. Des datations U/Pb de haute précision sur des cristaux de zircon et de baddeleyite ainsi que40Ar/39Ar sur des cris¬taux d'amphibole révèlent une initiation de la mise en place de PX1 il y a 22.1 ± 0,7 Ma; celle-ci a duré quelque 0,48 ± 0,22 à 0,52 ± 0,29 Ma. Ce laps de temps est compatible avec celui nécessaire à la cristallisation des filons individuels, qui va de moins d'une année lors de l'initiation du magmatisme à 5 ans lors du maximum d'activité de PX1. La présence de cristaux résorbés enregistrant une cristallisation complexe suggère l'existence d'une chambre mag¬matique convective sous-jacente à PX1 et périodiquement rechargée. Les compositions isotopiques des roches étu¬diées révèlent une source mantellique profonde de type point chaud avec une contribution du manteau lithosphéri- que métasomatisé présent sous les îles Canaries. Résumé L'intrusion mafique Miocène PX1 fait partie du soubassement superficiel (0.15-0.2 GPa, 1100 °Q d'un volcan d'île océanique. La particularité de ce pluton est l'existence d'alternances d'unités de gabbros et de pyroxénites qui met¬tent en évidence un litage magmatique vertical (NNE-SSW). Les horizons gabbroiques et pyroxénitiques sont constitués d'unités de différenciation métriques qui suggèrent tine mise en place par injections périodiques de filons verticaux de magma formant un complexe filonien. Chaque filon vertical a subi une différenciation parallèle à un front de solidification sub-vertical parallèle aux bords du filon. Les pyroxénites résultent du fractionnement et de l'accumulation d'olivine ± clinopyroxene ± plagioclase à partir d'un magma basaltique faiblement alcalin et sont interprétées comme étant des imités de différenciation tronquées dont le liquide interstitiel a été extrait par compaction. L'orientation préférentielle des clinopyroxènes dans ces pyroxe- nites (obtenues par analyse EBSD et micro-tomographique) révèle une composante de cisaillement simple dans la genèse de ces roches, ce qui confirme cette interprétation. La compaction des pyroxénites est probablement causée par a mise en place de filons de magma suivants. Le liquide interstitiel expulsé est probablement par ces derniers. Les clinopyroxènes des gabbros, montrent une composante de cisaillement pure suggérant qu'ils sont affectés par une déformation syn-magmatique parallèle aux zones de cisaillement NNE-SSW observées autour de PX1 et liées au contexte tectonique Miocène d'extension régionale. Ceci suggère que les gabbros sont liés à des taux de mise en place faibles à la fin de cycles d'activité magmatique et sont peu ou pas affectés par la compaction. L'initiation et la géométrie de PX1 sont donc contrôlées par le contexte tectonique régional d'extension alors que les taux et les volumes de magma dépendent de facteurs liés à la source. Des taux d'injection élevés résultent probable¬ment en une croissance du pluton supérieure à la place crée par cette extension. Dans ce cas de figure, la propagation des nouveaux dykes et l'inaptitude du magma à circuler à travers les anciens dykes cristallisés pourrait causer une augmentation de la pression non-lithostatique sur ces derniers, exprimée par un cisaillement simple et l'expulsion du liquide interstitiel qu'ils contiennent (documenté par les zones de collecte anorthositiques). Les compositions en éléments majeurs et traces des gabbros et pyroxenites de PX1 sont globalement homogènes et dépendent de la nature cumulative des échantillons. Cependant, de petites variations des concentrations en éléments traces ainsi que les teneurs en éléments traces des bordures de clinopyroxenes suggèrent que ces derniers ont subi un processus de rééquilibrage et de cristallisation in situ. L'homogénéité des compositions chimiques des échantillons, ainsi que la présence de grains de clinopyroxene résorbés suggère que le complexe filonien PX1 s'est mis en place au dessus d'une chambre magmatique périodiquement rechargée dans laquelle la convection est efficace. Chaque filon est donc issu d'un même magma, mais a subi une différenciation par cristallisation in situ (jusqu'à 70% de fraction¬nement) indépendamment des autres. Dans ces filons cristallisés, les minéraux cumulatifs subissent un rééquilibrage partiel avec les liquide interstitiel avant que ce dernier ne soit expulsé lors de la compaction (mettant ainsi un terme à la différenciation). Ce modèle de mise en place signifie qu'un minimum de 150Km3 de magma est nécessaire à la genèse de PX1, une partie de ce volume ayant été émis par le 'Central Volcanic Complex' de Fuerteventura. Les rapports isotopiques radiogéniques mesurés révèlent la contribution de trois pôles mantelliques dans la genèse du magma formant PX1. Le mélange de ces pôles HIMU, DMM et EM1 refléterai l'interaction du point chaud Cana¬rien avec un manteau lithosphérique hétérogène métasomatisé. Les petites variations de ces rapports et des teneurs en éléments traces au sein des faciès pourrait refléter des taux de fusion partielle variable de la source, résultant en un échantillonnage variable du manteau lithosphérique métasomatisé lors de son interaction avec le point chaud. Des datations U/Pb de haute précision (TIMS) sur des cristaux de zircon et de baddeleyite extraits de gabbros de PX1 révèlent que l'initiation de la cristallisation du magma a eu lieu il y a 22.10±0.07 Ma et que l'activité magmatique a duré un minimum de 0.48 à 0.52 Ma. Des âges 40Ar/39Ar obtenus sur amphibole sont de 21.9 ± 0.6 à 21.8 ± 0.3 Ma, identiques aux âges U/Pb. La combinaison de ces méthodes de datations, suggère que le temps maximum nécessaire à PX1 pour se refroidir en dessous de la température de fermeture de l'amphibole est de 0.8Ma. Ceci signifie que la durée de vie de PX1 est de 520 000 à 800 000 ans. La coexistence de cristaux de baddeleyite et de zircon dans un gabbro est attribuée à son interaction avec un fluide riche en C02 relâché par les carbonatites encaissantes lors du métamorphisme de contact généré par la mise en place de PX1 environ 160 000 ans après le début de sa mise en place. Les durées de vie obtenue sont en accord avec le modèle de mise en place suggérant une durée de cristallisation poux chaque filon allant de 1 an à 5 ans. Abstract The Miocene PX1 gabbro-pyroxenite intrusion (Fuerteventura, Canary Islands), is interpreted as the shallow-level feeder-zone (0.15-0.2 GPa and 1100-1120°C), to an ocean island volcano. The particularity of PX1 is that it displays a NNE-SSW trending vertical magmatic banding expressed by alternating gabbro and pyroxeriite sequences. The gabbro and pyroxenite sequences consist of metre-thick differentiation units, which suggest emplacement by pe¬riodic injection of magma pulses as vertical dykes that amalgamated, similarly to a sub-volcanic sheeted dyke com¬plex. Individual dykes underwent internal differentiation following a solidification front (favoured by a significant lateral/horizontal thermal gradient) parallel to the dyke edges. Pyroxenitic layers result from the fractionation and accumulation of clinopyroxene ± olivine ± plagioclase crystals from a mildly alkaline basaltic liquid and are interpre¬ted as truncated differentiation sequences, from which residual melts were extracted by compaction. Clinopyroxene mineral orientation in pyroxenites (evidenced by EBSD and micro X-ray tomography analysis) display a marked pure shear component, supporting this interpretation. Compaction and squeezing of the crystal mush is ascribed to the incoming and inflating magma pulses. The resulting expelled interstitial liquid was likely collected and erupted along with the magma flowing through the newly injected dykes. Gabbro sequences represent crystallised coalesced magma batches, emplaced at lower rates at the end of eruptive cycles, and underwent minor melt extraction as evi¬denced by clinopyroxene orientations that record a simple shear component suggesting syn-magmatic deformation parallel to observed NNF.-SSW trending shear-zones induced by the regional tensional Miocene stress-field. The initiation and geometry of PX1 is controlled by the regional extensional tectonic regime whereas rates and vo¬lumes of magma depend on source-related factors. High injection rates are likely to induce intrusion growth rates larger than could be accommodated by the regional extension. In this case, dyke tip geometry and the inability of magma to circulate through previously emplaced and crystallised dykes could result in an increase of non-lithostatic pressure on previously emplaced mushy dyke walls; generating strong pure-shear compaction and interstitial melt expulsion within the feeder-zone as recorded by the cumulitic pyroxenite bands and anorthositic collection zones. The whole-rock major and trace-element chemistry of PX1 gabbros and pyroxenites is globally homogeneous and controlled by the cumulate nature of the samples (i.e. on the modal proportions of olivine, pyroxene, plagioclase and oxides). However, small variations of whole-rock trace-element contents as well as trace-element contents of clinopyroxene rims suggest that in-situ re-equilibration and crystallisation has occurred. Additionally, the global homogeneity and presence of complex zoning of rare resorbed clinopyroxene crystals suggest that the PX1 feeder- zone overlies a periodically replenished and efficiently mixed magma chamber. Each individual dyke of magma thus originated from a compositionally constant mildly alkaline magma and differentiated independently from the others reaching up to 70% fractionation. Following dyke arrest these are affected by interaction with the trapped interstitial liquid prior to its compaction-linked expulsion (thus stopping the differentiation process). This emplacement model implies that minimum amount of approximately 150 km3 of magma is needed to generate PX1, part of it having been erupted through the overlying Central Volcanic Complex of Fuerteventura. The radiogenic isotope ratios of PX1 samples reveal the contribution on three end-members during magma genesis. This mixing of the H1MU, EMI and DMM end-members could reflect the interaction of the deep-seated Canarian mantle plume with a heterogeneous metasomatic and sepentininsed lithospheric mantle. Additionally, the observed trace-element and isotopic variations within the same fades groups could reflect varying degrees of partial melting of the source region, thus tapping more or less large areas of the metasomatised lithospheric mantle during interac¬tion with the plume. High precision ID-TIMS U/Pb zircon and baddeleyite ages from the PX1 gabbro samples, indicate initiation of magma crystallisation at 22.10 ± 0.07 Ma. The magmatic activity lasted a minimum of 0.48 to 0.52 Ma. 40Ar/39Ar amphibole ages are of 21.9 ± 0.6 to 21.8 ± 0.3, identical within errors to the U/Pb ages. The combination of the 40Ar/39Ar and U/Pb datasets imply that the maximum amount of time PX1 took to cool below amphibole Tc is 0.8 Ma, suggesting PX1 lifetime of 520 000 to 800 000 years. On top of this, the coexistence of baddeleyite and zircon in a single sample is ascribed to the interaction of PX1 with C02-rich carbonatite-derived fluids released from the host-rock carbonatites during contact metamorphism 160 000 years after PX1 initiation. These ages are in agreement with the emplacement model, implying a crystallisation time of less than 1 to 5 years for individual dykes.

Streptabody, a high avidity molecule made by tetramerization of in vivo biotinylated, phage display-selected scFv fragments on streptavidin.

Phage display is a powerful method of isolating of antibody fragments from highly diverse naive human antibody repertoires. However, the affinity of the selected antibodies is usually low and current methods of affinity maturation are complex and time-consuming. In this paper, we describe an easy way to increase the functional affinity (avidity) of single chain variable fragments (scFvs) by tetramerization on streptavidin, following their site-specific biotinylation by the enzyme BirA. Expression vectors have been constructed that enable addition of the 15 amino acid biotin acceptor domain (BAD) on selected scFvs. Different domains were cloned at the C-terminus of scFv in the following order: a semi-rigid hinge region (of 16 residues), the BAD, and a histidine tail. Two such recombinant scFvs directed against the carcinoembryonic antigen (CEA) were previously selected from human non-immune and murine immune phage display libraries. The scFvs were first synthesized in Escherichia coli carrying the plasmid encoding the BirA enzyme, and then purified from the cytoplasmic extracts by Ni-NTA affinity chromatography. Purified biotinylated scFvs were tetramerized on the streptavidin molecule to create a streptabody (StAb). The avidity of various forms of anti-CEA StAbs, tested on purified CEA by competitive assays and surface plasmon resonance showed an increase of more than one log, as compared with the scFv monomer counterparts. Furthermore, the percentage of direct binding of 125I-labeled StAb or monomeric scFv on CEA-Sepharose beads and on CEA-expressing cells showed a dramatic increase for the tetramerized scFv (>80%), as compared with the monomeric scFv (<20%). Interestingly, the percentage binding of 125I-labeled anti-CEA StAbs to CEA-expressing colon carcinoma cells was definitely higher (>80%) than that obtained with a reference high affinity murine anti-CEA mAb (30%). Another advantage of using scFvs in a StAb format was demonstrated by Western blot analysis, where tetramerized anti-CEA scFv could detect a small quantity of CEA at a concentration 100-fold lower than the monomeric scFv.

Recurrent Macular Edema in Central Retinal Vein Occlusion Treated with Intravitreal Ranibizumab using a Modified Treat and Extend Regimen.

Background: The aim of this study was to evaluate the stability over time of the individually defined interval of intravitreal ranibizumab injection (IVR) for the treatment of recurrent macular edema (ME) in central retinal vein occlusion (CRVO). Patients and Methods: A case series of treatment naïve patients followed in the Jules Gonin Eye Hospital for macular edema due to central retinal vein occlusion is presented. Patients were treated monthly with IVR until complete absence of fluid on qualitative SD-OCT with a minimum of 5 monthly IVR. Thereafter, they were followed according to a modified treat and extend regimen (mTER). Results: Twelve eyes (12 patients) with ME due to CRVO were included. The mean follow-up period was 31.3 months. Analysis showed that best corrected visual acuity (BCVA), central macular thickness and qualitative spectral domain optical coherence tomography (SD-OCT) showed comparable results under monthly interval, after titration of an individualized interval and when performed in a series. 78 % of treating intervals were within ± 2 weeks of the first individually adjusted interval. The mean first defined interval was 4.3 weeks and the mean interval over time was 5.5 weeks (p = 0.003). There was a trend towards longer interval over time. Conclusion: The adjusted interval of retreatment of patients with ME due to CRVO showed a high stability with a trend toward longer duration over time. An mTER regimen seems to be valuable to follow patients with ME with good stabilization of VA.

A comparative study of ballpoint ink ageing parameters using GC/MS.

For more than a decade scientists tried to develop methods capable of dating ink by monitoring the loss of phenoxyethanol (PE) over time. While many methods were proposed in the literature, few were really used to solve practical cases and they still raise much concern within the scientific community. In fact, due to the complexity of ink drying processes it is particularly difficult to find a reliable ageing parameter to reproducibly follow ink ageing. Moreover, systematic experiments are required in order to evaluate how different factors actually influence the results over time. Therefore, this work aimed at evaluating the capacity of four different ageing parameters to reliably follow ink ageing over time: (1) the quantity of solvent PE in an ink line, (2) the relative peak area (RPA) normalising the PE results using stable volatile compounds present in the ink formulation, (3) the solvent loss ratio (R%) calculated from PE results obtained by the analyses of naturally and artificially aged samples, (4) a modified solvent loss ratio version (R%*) calculated from RPA results. After the determination of the limits of reliable measurements of the analytical method, the repeatability of the different ageing parameters was evaluated over time, as well as the influence of ink composition, writing pressure and storage conditions on the results. Surprisingly, our results showed that R% was not the most reliable parameter, as it showed the highest standard deviation. Discussion of the results in an ink dating perspective suggests that other proposed parameters, such as RPA values, may be more adequate to follow ink ageing over time.

A comparative study of ballpoint ink ageing parameters using GC/MS

For more than a decade scientists tried to develop methods capable of dating ink by monitoring the loss of phenoxyethanol (PE) over time. While many methods were proposed in the literature, few were really used to solve practical cases and they still raise much concern within the scientific community. In fact, due to the complexity of ink drying processes it is particularly difficult to find a reliable ageing parameter to reproducibly follow ink ageing. Moreover, systematic experiments are required in order to evaluate how different factors actually influence the results over time. Therefore, this work aimed at evaluating the capacity of four different ageing parameters to reliably follow ink ageing over time: (1) the quantity of solvent PE in an ink line, (2) the relative peak area (RPA) normalising the PE results using stable volatile compounds present in the ink formulation, (3) the solvent loss ratio (R%) calculated from PE results obtained by the analyses of naturally and artificially aged samples, (4) a modified solvent loss ratio version (R%*) calculated from RPA results. After the determination of the limits of reliable measurements of the analytical method, the repeatability of the different ageing parameters was evaluated over time, as well as the influence of ink composition, writing pressure and storage conditions on the results. Surprisingly, our results showed that R% was not the most reliable parameter, as it showed the highest standard deviation. Discussion of the results in an ink dating perspective suggests that other proposed parameters, such as RPA values, may be more adequate to follow ink ageing over time.

Obstructive Sleep Apnea Severity and Overnight Body Fluid Shift before and after Hemodialysis.

BACKGROUND AND OBJECTIVES: Obstructive sleep apnea is associated with significantly increased cardiovascular morbidity and mortality. Fluid overload may promote obstructive sleep apnea in patients with ESRD through an overnight fluid shift from the legs to the neck soft tissues. Body fluid shift and severity of obstructive sleep apnea before and after hemodialysis were compared in patients with ESRD. DESIGN, SETTING, PARTICIPANTS, & MEASUREMENTS: Seventeen patients with hemodialysis and moderate to severe obstructive sleep apnea were included. Polysomnographies were performed the night before and after hemodialysis to assess obstructive sleep apnea, and bioimpedance was used to measure fluid overload and leg fluid volume. RESULTS: The mean overnight rostral fluid shift was 1.27±0.41 L prehemodialysis; it correlated positively with fluid overload volume (r=0.39; P=0.02) and was significantly lower posthemodialysis (0.78±0.38 L; P<0.001). There was no significant difference in the mean obstructive apnea-hypopnea index before and after hemodialysis (46.8±22.0 versus 42.1±18.6 per hour; P=0.21), but obstructive apnea-hypopnea index was significantly lower posthemodialysis (-10.1±10.8 per hour) in the group of 12 patients, with a concomitant reduction of fluid overload compared with participants without change in fluid overload (obstructive apnea-hypopnea index +8.2±16.1 per hour; P<0.01). A lower fluid overload after hemodialysis was significantly correlated (r=0.49; P=0.04) with a lower obstructive apnea-hypopnea index. Fluid overload-assessed by bioimpedance-was the best predictor of the change in obstructive apnea-hypopnea index observed after hemodialysis (standardized r=-0.68; P=0.01) in multivariate regression analysis. CONCLUSIONS: Fluid overload influences overnight rostral fluid shift and obstructive sleep apnea severity in patients with ESRD undergoing intermittent hemodialysis. Although no benefit of hemodialysis on obstructive sleep apnea severity was observed in the whole group, the change in obstructive apnea-hypopnea index was significantly correlated with the change in fluid overload after hemodialysis. Moreover, the subgroup with lower fluid overload posthemodialysis showed a significantly lower obstructive sleep apnea severity, which provides a strong incentive to further study whether optimizing fluid status in patients with obstructive sleep apnea and ESRD will improve the obstructive apnea-hypopnea index.

Neonates with Bartter syndrome have enormous fluid and sodium requirements.

AIM: Managing neonatal Bartter syndrome by achieving adequate weight gain is challenging. We assessed the correlation between weight gain in neonatal Bartter syndrome and the introduction of fluid and sodium supplementations and indomethacin during the first 4 weeks of life. METHODS: Daily fluid and electrolytes requirements were analysed using linear regression and Spearman correlation coefficients. The weight gain coefficient was calculated as daily weight gain after physiological neonatal weight loss. RESULTS: We studied seven infants. The highest weight gain coefficients occurred between weeks two and four in the five neonates who either received prompt amounts of fluid (maximum 810 mL/kg/day) and sodium (maximum 70 mmol/kg/day) or were treated with indomethacin. For the two patients with the highest weight gain coefficient, water and sodium supplementations were decreased in weeks two to four leading to a significant negative Spearman correlation between weight gain and fluid supplements (r = -0.55 and -0.68) and weight gain and sodium supplementations (r = -0.96 and -0.72). The two patients with the lowest weight gain coefficient had positive Spearman correlation coefficients between weight gain and fluid and sodium supplementations. CONCLUSION: Infants with neonatal Bartter syndrome required rapid and enormous fluid and sodium supplementations or the early introduction of indomethacin treatment to achieve adequate weight gain during the early postnatal period.

External urethral sphincter electromyography in asymptomatic women and the influence of the menstrual cycle.

OBJECTIVE: To investigate by electromyography (EMG), the presence of complex repetitive discharges (CRDs) and decelerating bursts (DBs) in the striated external urethral sphincter during the menstrual cycle in female volunteers with no urinary symptoms and complete bladder emptying. SUBJECTS AND METHODS: Healthy female volunteers aged 20-40 years, with regular menstrual cycles and no urinary symptoms were recruited. Volunteers completed a menstruation chart, urinary symptom questionnaires, pregnancy test, urine dipstick, urinary free flow and post-void ultrasound bladder scan. Exclusion criteria included current pregnancy, use of hormonal medication or contraception, body mass index of >35 kg/m(2) , incomplete voiding and a history of pelvic surgery. Eligible participants underwent an external urethral sphincter EMG, using a needle electrode in the early follicular phase and the mid-luteal phase of their menstrual cycles. Serum oestradiol and progesterone were measured at each EMG test. RESULTS: In all, 119 women enquired about the research and following screening, 18 were eligible to enter the study phase. Complete results were obtained in 15 women. In all, 30 EMG tests were undertaken in the 15 asymptomatic women. Sphincter EMG was positive for CRDs and DBs at one or both phases of the menstrual cycle in eight (53%) of the women. Three had CRDs and DBs in both early follicular and mid-luteal phases. Five had normal EMG activity in the early follicular phase and CRDs and DBs in the mid-luteal phase. No woman had abnormal EMG activity in the early follicular phase and normal activity in the luteal phase. There was no relationship between EMG activity and age, parity or serum levels of oestradiol and progesterone. CONCLUSIONS: CRDs and DB activity in the external striated urethral sphincter is present in a high proportion of asymptomatic young women. This abnormal EMG activity has been shown for the first time to change during the menstrual cycle in individual women. CRDs and DBs are more commonly found in the luteal phase of the menstrual cycle. The importance of CRDs and DBs in the aetiology of urinary retention in young women remains uncertain. The distribution and or quantity of abnormal EMG activity in the external urethral sphincter may be important. In a woman with urinary retention the finding of CRDs and DBs by needle EMG does not automatically establish Fowler's syndrome as the explanation for the bladder dysfunction. Urethral pressure profilometry may be helpful in establishing a diagnosis. Opiate use and psychological stress should be considered in young women with urinary retention.