Time-lapse and probabilistic inversion strategies for plane-wave electromagnetic data.

L'utilisation efficace des systèmes géothermaux, la séquestration du CO2 pour limiter le changement climatique et la prévention de l'intrusion d'eau salée dans les aquifères costaux ne sont que quelques exemples qui démontrent notre besoin en technologies nouvelles pour suivre l'évolution des processus souterrains à partir de la surface. Un défi majeur est d'assurer la caractérisation et l'optimisation des performances de ces technologies à différentes échelles spatiales et temporelles. Les méthodes électromagnétiques (EM) d'ondes planes sont sensibles à la conductivité électrique du sous-sol et, par conséquent, à la conductivité électrique des fluides saturant la roche, à la présence de fractures connectées, à la température et aux matériaux géologiques. Ces méthodes sont régies par des équations valides sur de larges gammes de fréquences, permettant détudier de manières analogues des processus allant de quelques mètres sous la surface jusqu'à plusieurs kilomètres de profondeur. Néanmoins, ces méthodes sont soumises à une perte de résolution avec la profondeur à cause des propriétés diffusives du champ électromagnétique. Pour cette raison, l'estimation des modèles du sous-sol par ces méthodes doit prendre en compte des informations a priori afin de contraindre les modèles autant que possible et de permettre la quantification des incertitudes de ces modèles de façon appropriée. Dans la présente thèse, je développe des approches permettant la caractérisation statique et dynamique du sous-sol à l'aide d'ondes EM planes. Dans une première partie, je présente une approche déterministe permettant de réaliser des inversions répétées dans le temps (time-lapse) de données d'ondes EM planes en deux dimensions. Cette stratégie est basée sur l'incorporation dans l'algorithme d'informations a priori en fonction des changements du modèle de conductivité électrique attendus. Ceci est réalisé en intégrant une régularisation stochastique et des contraintes flexibles par rapport à la gamme des changements attendus en utilisant les multiplicateurs de Lagrange. J'utilise des normes différentes de la norme l2 pour contraindre la structure du modèle et obtenir des transitions abruptes entre les régions du model qui subissent des changements dans le temps et celles qui n'en subissent pas. Aussi, j'incorpore une stratégie afin d'éliminer les erreurs systématiques de données time-lapse. Ce travail a mis en évidence l'amélioration de la caractérisation des changements temporels par rapport aux approches classiques qui réalisent des inversions indépendantes à chaque pas de temps et comparent les modèles. Dans la seconde partie de cette thèse, j'adopte un formalisme bayésien et je teste la possibilité de quantifier les incertitudes sur les paramètres du modèle dans l'inversion d'ondes EM planes. Pour ce faire, je présente une stratégie d'inversion probabiliste basée sur des pixels à deux dimensions pour des inversions de données d'ondes EM planes et de tomographies de résistivité électrique (ERT) séparées et jointes. Je compare les incertitudes des paramètres du modèle en considérant différents types d'information a priori sur la structure du modèle et différentes fonctions de vraisemblance pour décrire les erreurs sur les données. Les résultats indiquent que la régularisation du modèle est nécessaire lorsqu'on a à faire à un large nombre de paramètres car cela permet d'accélérer la convergence des chaînes et d'obtenir des modèles plus réalistes. Cependent, ces contraintes mènent à des incertitudes d'estimations plus faibles, ce qui implique des distributions a posteriori qui ne contiennent pas le vrai modèledans les régions ou` la méthode présente une sensibilité limitée. Cette situation peut être améliorée en combinant des méthodes d'ondes EM planes avec d'autres méthodes complémentaires telles que l'ERT. De plus, je montre que le poids de régularisation des paramètres et l'écart-type des erreurs sur les données peuvent être retrouvés par une inversion probabiliste. Finalement, j'évalue la possibilité de caractériser une distribution tridimensionnelle d'un panache de traceur salin injecté dans le sous-sol en réalisant une inversion probabiliste time-lapse tridimensionnelle d'ondes EM planes. Etant donné que les inversions probabilistes sont très coûteuses en temps de calcul lorsque l'espace des paramètres présente une grande dimension, je propose une stratégie de réduction du modèle ou` les coefficients de décomposition des moments de Legendre du panache de traceur injecté ainsi que sa position sont estimés. Pour ce faire, un modèle de résistivité de base est nécessaire. Il peut être obtenu avant l'expérience time-lapse. Un test synthétique montre que la méthodologie marche bien quand le modèle de résistivité de base est caractérisé correctement. Cette méthodologie est aussi appliquée à un test de trac¸age par injection d'une solution saline et d'acides réalisé dans un système géothermal en Australie, puis comparée à une inversion time-lapse tridimensionnelle réalisée selon une approche déterministe. L'inversion probabiliste permet de mieux contraindre le panache du traceur salin gr^ace à la grande quantité d'informations a priori incluse dans l'algorithme. Néanmoins, les changements de conductivités nécessaires pour expliquer les changements observés dans les données sont plus grands que ce qu'expliquent notre connaissance actuelle des phénomenès physiques. Ce problème peut être lié à la qualité limitée du modèle de résistivité de base utilisé, indiquant ainsi que des efforts plus grands devront être fournis dans le futur pour obtenir des modèles de base de bonne qualité avant de réaliser des expériences dynamiques. Les études décrites dans cette thèse montrent que les méthodes d'ondes EM planes sont très utiles pour caractériser et suivre les variations temporelles du sous-sol sur de larges échelles. Les présentes approches améliorent l'évaluation des modèles obtenus, autant en termes d'incorporation d'informations a priori, qu'en termes de quantification d'incertitudes a posteriori. De plus, les stratégies développées peuvent être appliquées à d'autres méthodes géophysiques, et offrent une grande flexibilité pour l'incorporation d'informations additionnelles lorsqu'elles sont disponibles. -- The efficient use of geothermal systems, the sequestration of CO2 to mitigate climate change, and the prevention of seawater intrusion in coastal aquifers are only some examples that demonstrate the need for novel technologies to monitor subsurface processes from the surface. A main challenge is to assure optimal performance of such technologies at different temporal and spatial scales. Plane-wave electromagnetic (EM) methods are sensitive to subsurface electrical conductivity and consequently to fluid conductivity, fracture connectivity, temperature, and rock mineralogy. These methods have governing equations that are the same over a large range of frequencies, thus allowing to study in an analogous manner processes on scales ranging from few meters close to the surface down to several hundreds of kilometers depth. Unfortunately, they suffer from a significant resolution loss with depth due to the diffusive nature of the electromagnetic fields. Therefore, estimations of subsurface models that use these methods should incorporate a priori information to better constrain the models, and provide appropriate measures of model uncertainty. During my thesis, I have developed approaches to improve the static and dynamic characterization of the subsurface with plane-wave EM methods. In the first part of this thesis, I present a two-dimensional deterministic approach to perform time-lapse inversion of plane-wave EM data. The strategy is based on the incorporation of prior information into the inversion algorithm regarding the expected temporal changes in electrical conductivity. This is done by incorporating a flexible stochastic regularization and constraints regarding the expected ranges of the changes by using Lagrange multipliers. I use non-l2 norms to penalize the model update in order to obtain sharp transitions between regions that experience temporal changes and regions that do not. I also incorporate a time-lapse differencing strategy to remove systematic errors in the time-lapse inversion. This work presents improvements in the characterization of temporal changes with respect to the classical approach of performing separate inversions and computing differences between the models. In the second part of this thesis, I adopt a Bayesian framework and use Markov chain Monte Carlo (MCMC) simulations to quantify model parameter uncertainty in plane-wave EM inversion. For this purpose, I present a two-dimensional pixel-based probabilistic inversion strategy for separate and joint inversions of plane-wave EM and electrical resistivity tomography (ERT) data. I compare the uncertainties of the model parameters when considering different types of prior information on the model structure and different likelihood functions to describe the data errors. The results indicate that model regularization is necessary when dealing with a large number of model parameters because it helps to accelerate the convergence of the chains and leads to more realistic models. These constraints also lead to smaller uncertainty estimates, which imply posterior distributions that do not include the true underlying model in regions where the method has limited sensitivity. This situation can be improved by combining planewave EM methods with complimentary geophysical methods such as ERT. In addition, I show that an appropriate regularization weight and the standard deviation of the data errors can be retrieved by the MCMC inversion. Finally, I evaluate the possibility of characterizing the three-dimensional distribution of an injected water plume by performing three-dimensional time-lapse MCMC inversion of planewave EM data. Since MCMC inversion involves a significant computational burden in high parameter dimensions, I propose a model reduction strategy where the coefficients of a Legendre moment decomposition of the injected water plume and its location are estimated. For this purpose, a base resistivity model is needed which is obtained prior to the time-lapse experiment. A synthetic test shows that the methodology works well when the base resistivity model is correctly characterized. The methodology is also applied to an injection experiment performed in a geothermal system in Australia, and compared to a three-dimensional time-lapse inversion performed within a deterministic framework. The MCMC inversion better constrains the water plumes due to the larger amount of prior information that is included in the algorithm. The conductivity changes needed to explain the time-lapse data are much larger than what is physically possible based on present day understandings. This issue may be related to the base resistivity model used, therefore indicating that more efforts should be given to obtain high-quality base models prior to dynamic experiments. The studies described herein give clear evidence that plane-wave EM methods are useful to characterize and monitor the subsurface at a wide range of scales. The presented approaches contribute to an improved appraisal of the obtained models, both in terms of the incorporation of prior information in the algorithms and the posterior uncertainty quantification. In addition, the developed strategies can be applied to other geophysical methods, and offer great flexibility to incorporate additional information when available.

Tandem repeat sequence variation as causative cis-eQTLs for protein-coding gene expression variation: the case of CSTB.

Association studies have revealed expression quantitative trait loci (eQTLs) for a large number of genes. However, the causative variants that regulate gene expression levels are generally unknown. We hypothesized that copy-number variation of sequence repeats contribute to the expression variation of some genes. Our laboratory has previously identified that the rare expansion of a repeat c.-174CGGGGCGGGGCG in the promoter region of the CSTB gene causes a silencing of the gene, resulting in progressive myoclonus epilepsy. Here, we genotyped the repeat length and quantified CSTB expression by quantitative real-time polymerase chain reaction in 173 lymphoblastoid cell lines (LCLs) and fibroblast samples from the GenCord collection. The majority of alleles contain either two or three copies of this repeat. Independent analysis revealed that the c.-174CGGGGCGGGGCG repeat length is strongly associated with CSTB expression (P = 3.14 × 10(-11)) in LCLs only. Examination of both genotyped and imputed single-nucleotide polymorphisms (SNPs) within 2 Mb of CSTB revealed that the dodecamer repeat represents the strongest cis-eQTL for CSTB in LCLs. We conclude that the common two or three copy variation is likely the causative cis-eQTL for CSTB expression variation. More broadly, we propose that polymorphic tandem repeats may represent the causative variation of a fraction of cis-eQTLs in the genome.

Gp74 a membrane glycoprotein of the cis-Golgi network that cycles through the endoplasmic reticulum and intermediate compartment

A monoclonal antibody CC92 (IgM), raised against a fraction of rat liver enriched in Golgi membranes, recognizes a novel Endo H-resistant 74-kD membrane glycoprotein (gp74). The bulk of gp74 is confined to the cis-Golgi network (CGN). Outside the Golgi gp74 is found in tubulovesicular structures and ER foci. In cells incubated at 37 degrees C the majority of gp74 is segregated from the intermediate compartment (IC) marker p58. However, in cells treated with organelle perturbants such as low temperature, BFA, and [AIF4]- the patterns of the two proteins become indistinguishable. Both proteins are retained in the Golgi complex at 20 degrees C and in the IC at 15 degrees C. Incubation of cells with BFA results in relocation of gp74 to p58 positive IC elements. [AIF4]- induces the redistribution of gp74 from the Golgi to p58-positive vesicles and does not retard the translocation of gp74 to IC elements in cells treated with BFA. Disruption of microtubules by nocodazol results in the rapid disappearance of the Golgi elements stained by gp74 and redistribution of the protein into vesicle-like structures. The responses of gp74 to cell perturbants are in sharp contrast with those of cis/middle and trans-Golgi resident proteins whose location is not affected by low temperatures or [AIF4]-, are translocated to the ER upon addition of BFA, and stay in slow disintegrating Golgi elements in cells treated with nocodazol. The results suggest that gp74 is an itinerant protein that resides most of the time in the CGN and cycles through the ER/IC following the pathway used by p58.

Inflammation palpébrale pseudotumorale. A propos d'une observation anatomo-clinique [Pseudotumoral eyelid inflammation. Apropos of an anatomo-clinical case].

A clinicopathological case of a 76-year-old male patient with a chronic inflammatory change of the inferior left eyelid is reported. The inflammation appeared as a reddish area of the inner part of the eyelid, without sharp limits, but with loss of lashes. Numerous local treatments did not to cure this condition. As some true eyelid tumors may mimic an inflammation during growth and, for example, sebaceous carcinoma may clinically present as chronic unilateral blepharitis, a surgical excisional biopsy was performed on this left eyelid. Its histopathological study showed a granulomatous inflammation, which was typical of a simple chalazion. This case clearly illustrates that the chalazion may not always appear as a limited nodular inflammation of the eyelid, but may have a more diffuse clinical presentation.

Epitaxial growth of biferroic YMnO3(0001) on platinum electrodes

Epitaxial films of the biferroic YMnO3 (YMO) oxide have been grown on platinum-coated SrTiO3(1 1 1) and Al2O3(0 0 0 1) substrates. The platinum electrodes, (1 1 1) oriented, are templates for the epitaxy of the hexagonal phase of YMO with a (0 0 0 1) out-of-plane orientation, which is of interest as this is the polarization direction of YMO. X-ray diffractometry indicates the presence of two crystal domains, 60° rotated in-plane, in the Pt(1 1 1) layers which subsequently are transferred on the upperlaying YMO. Cross-section analysis by high-resolution transmission electron microscopy (HRTEM) of YMnO3/Pt/SrTiO3(1 1 1) shows high-quality epitaxy and sharp interfaces across the structure in the observed region. We present a detailed study of the epitaxial growth of the hexagonal YMO on the electrodes.

Extraction and concentration of biogenic calcium oxalate from plant leaves

The objective of this study was to extract and concentrate calcium oxalate (CaOx) crystals from plant leaves that form the above mentioned crystals. The chemical and physical studies of CaOx from plant to be performed depend on an adequate amount of the crystals. The plant used in this study was croton (Codiaeum variegatum). The leaves were ground in a heavy duty blender and sieved through a 0.20 mm sieve. The suspension obtained was suspended in distilled water. The crystals were concentrated at the bottom of a test tube. The supernatant must be washed until it is free of plant pigments and other organic substances. Biogenic CaOx crystals have well-defined and sharp peaks, indicating very high crystallinity. Moreover, the CaOx crystals were not damaged during the extraction procedure, as can be seen on the scanning electron microscope images. The porposed method can be considered efficient to extract and concentrate biogenic calcium oxalate.

Identification of heart rate-associated loci and their effects on cardiac conduction and rhythm disorders.

Elevated resting heart rate is associated with greater risk of cardiovascular disease and mortality. In a 2-stage meta-analysis of genome-wide association studies in up to 181,171 individuals, we identified 14 new loci associated with heart rate and confirmed associations with all 7 previously established loci. Experimental downregulation of gene expression in Drosophila melanogaster and Danio rerio identified 20 genes at 11 loci that are relevant for heart rate regulation and highlight a role for genes involved in signal transmission, embryonic cardiac development and the pathophysiology of dilated cardiomyopathy, congenital heart failure and/or sudden cardiac death. In addition, genetic susceptibility to increased heart rate is associated with altered cardiac conduction and reduced risk of sick sinus syndrome, and both heart rate-increasing and heart rate-decreasing variants associate with risk of atrial fibrillation. Our findings provide fresh insights into the mechanisms regulating heart rate and identify new therapeutic targets.

Neurogenic subventricular zone stem/progenitor cells are notch1-dependent in their active but not quiescent state.

The adult mammalian forebrain contains neural stem/progenitor cells (NSCs) that generate neurons throughout life. As in other somatic stem cell systems, NSCs are proposed to be predominantly quiescent and proliferate only sporadically to produce more committed progeny. However, quiescence has recently been shown not to be an essential criterion for stem cells. It is not known whether NSCs show differences in molecular dependence based on their proliferation state. The subventricular zone (SVZ) of the adult mouse brain has a remarkable capacity for repair by activation of NSCs. The molecular interplay controlling adult NSCs during neurogenesis or regeneration is not clear but resolving these interactions is critical in order to understand brain homeostasis and repair. Using conditional genetics and fate mapping, we show that Notch signaling is essential for neurogenesis in the SVZ. By mosaic analysis, we uncovered a surprising difference in Notch dependence between active neurogenic and regenerative NSCs. While both active and regenerative NSCs depend upon canonical Notch signaling, Notch1-deletion results in a selective loss of active NSCs (aNSCs). In sharp contrast, quiescent NSCs (qNSCs) remain after Notch1 ablation until induced during regeneration or aging, whereupon they become Notch1-dependent and fail to fully reinstate neurogenesis. Our results suggest that Notch1 is a key component of the adult SVZ niche, promoting maintenance of aNSCs, and that this function is compensated in qNSCs. Therefore, we confirm the importance of Notch signaling for maintaining NSCs and neurogenesis in the adult SVZ and reveal that NSCs display a selective reliance on Notch1 that may be dictated by mitotic state.

Highway Relocation: U.S. 30/Lincoln Highway film production, 1960

The film depicts period traffic congestion, sharp and winding sections of road, steep hills making trucks slow to a crawl, and dangerous vehicle and pedestrian crossings, all important reasons why highway design and safety improvements, and highway relocation were needed. In fact, when the film was produced, U.S. 30 or the Lincoln Highway was the busiest primary road in Iowa; and the section between State Center and Boone was deemed “critical,” meaning it was considered dangerous by the ISHC’s Efficiency Standards.

Phase-field model for Hele-Shaw flows with arabitrary contrast II: numerical approach

We present a phase-field model for the dynamics of the interface between two inmiscible fluids with arbitrary viscosity contrast in a rectangular Hele-Shaw cell. With asymptotic matching techniques we check the model to yield the right Hele-Shaw equations in the sharp-interface limit, and compute the corrections to these equations to first order in the interface thickness. We also compute the effect of such corrections on the linear dispersion relation of the planar interface. We discuss in detail the conditions on the interface thickness to control the accuracy and convergence of the phase-field model to the limiting Hele-Shaw dynamics. In particular, the convergence appears to be slower for high viscosity contrasts.

Reaction-diffusion fronts under stochastic advection

We study front propagation in stirred media using a simplified modelization of the turbulent flow. Computer simulations reveal the existence of the two limiting propagation modes observed in recent experiments with liquid phase isothermal reactions. These two modes respectively correspond to a wrinkled although sharp propagating interface and to a broadened one. Specific laws relative to the enhancement of the front velocity in each regime are confirmed by our simulations.

Front dynamics in turbulent media

A study of a stable front propagating in a turbulent medium is presented. The front is generated through a reaction-diffusion equation, and the turbulent medium is statistically modeled using a Langevin equation. Numerical simulations indicate the presence of two different dynamical regimes. These regimes appear when the turbulent flow either wrinkles a still rather sharp propagating interfase or broadens it. Specific dependences of the propagating velocities on stirring intensities appropriate to each case are found and fitted when possible according to theoretically predicted laws. Different turbulent spectra are considered.

Low viscosity contrast fingering in a rotating Hele-Shaw cell

We study the fingering instability of a circular interface between two immiscible liquids in a radial Hele-Shaw cell. The cell rotates around its vertical symmetry axis, and the instability is driven by the density difference between the two fluids. This kind of driving allows studying the interfacial dynamics in the particularly interesting case of an interface separating two liquids of comparable viscosity. An accurate experimental study of the number of fingers emerging from the instability reveals a slight but systematic dependence of the linear dispersion relation on the gap spacing. We show that this result is related to a modification of the interface boundary condition which incorporates stresses originated from normal velocity gradients. The early nonlinear regime shows nearly no competition between the outgrowing fingers, characteristic of low viscosity contrast flows. We perform experiments in a wide range of experimental parameters, under conditions of mass conservation (no injection), and characterize the resulting patterns by data collapses of two characteristic lengths: the radius of gyration of the pattern and the interface stretching. Deep in the nonlinear regime, the fingers which grow radially outwards stretch and become gradually thinner, to a point that the fingers pinch and emit drops. We show that the amount of liquid emitted in the first generation of drops is a constant independent of the experimental parameters. Further on there is a sharp reduction of the amount of liquid centrifugated, punctuated by periods of no observable centrifugation.

Asymmetric POSTS associated with unilateral EEG abnormalities.

OBJECTIVE: Positive occipital sharp transient of the sleep (POSTS) are considered a normal variant of non-REM sleep EEG. We describe a small series of patients with asymmetric POSTS and ipsilateral abnormal EEG findings. METHODS: Over a period of 30 weeks, we prospectively observed five consecutive subjects with strictly unilateral POSTS associated with ispilateral electrographic abnormalities. They represent 0.4% of all EEG performed over the same time lapse (5/1130), including inpatients, outpatients and long-term monitoring. RESULTS: Four women and one boy suffering from epileptic seizures (aged 7-76 years old) had unilateral POSTS, occurring only on the right side, during light sleep. They also presented ipsilateral epileptiform abnormalities. CONCLUSION: The fact that POSTS were asymmetric and found only on the same side as the abnormalities raises the question whether these transients should still be considered physiological or could be interpreted at times as markers of underlying electrical abnormalities, pointing to an increased cortical excitability on the more active side. Although larger samples are needed to confirm our preliminary results, this case study questions the interpretation of POSTS as a uniformly normal variant.

Origin of the neutron skin thickness of 208Pb in nuclear mean-field models

We study whether the neutron skin thickness Δrnp of 208Pb originates from the bulk or from the surface of the nucleon density distributions, according to the mean-field models of nuclear structure, and find that it depends on the stiffness of the nuclear symmetry energy. The bulk contribution to Δrnp arises from an extended sharp radius of neutrons, whereas the surface contribution arises from different widths of the neutron and proton surfaces. Nuclear models where the symmetry energy is stiff, as typical of relativistic models, predict a bulk contribution in Δrnp of 208Pb about twice as large as the surface contribution. In contrast, models with a soft symmetry energy like common nonrelativistic models predict that Δrnp of 208Pb is divided similarly into bulk and surface parts. Indeed, if the symmetry energy is supersoft, the surface contribution becomes dominant. We note that the linear correlation of Δrnp of 208Pb with the density derivative of the nuclear symmetry energy arises from the bulk part of Δrnp. We also note that most models predict a mixed-type (between halo and skin) neutron distribution for 208Pb. Although the halo-type limit is actually found in the models with a supersoft symmetry energy, the skin-type limit is not supported by any mean-field model. Finally, we compute parity-violating electron scattering in the conditions of the 208Pb parity radius experiment (PREX) and obtain a pocket formula for the parity-violating asymmetry in terms of the parameters that characterize the shape of the 208Pb nucleon densities.