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.

PDMP blocks brefeldin a-induced retrograde membrane transport from Golgi to ER: evidence for involvement of calcium homeostasis and dissociation from sphingolipid metabolism

In this study, we show that an inhibitor of sphingolipid biosynthesis, d,l-threo-1-phenyl-2- decanoylamino-3-morpholino-1-propanol (PDMP), inhibits brefeldin A (BFA)-induced retrograde membrane transport from Golgi to endoplasmic reticulum (ER). If BFA treatment was combined with or preceded by PDMP administration to cells, disappearance of discrete Golgi structures did not occur. However, when BFA was allowed to exert its effect before PDMP addition, PDMP could not ¿rescue¿ the Golgi compartment. Evidence is presented showing that this action of PDMP is indirect, which means that the direct target is not sphingolipid metabolism at the Golgi apparatus. A fluorescent analogue of PDMP, 6-(N-[7-nitro-2,1,3-benzoxadiazol-4-yl]amino)hexanoyl-PDMP (C6-NBD-PDMP), did not localize in the Golgi apparatus. Moreover, the effect of PDMP on membrane flow did not correlate with impaired C6-NBD-sphingomyelin biosynthesis and was not mimicked by exogenous C6-ceramide addition or counteracted by exogenous C6-glucosylceramide addition. On the other hand, the PDMP effect was mimicked by the multidrug resistance protein inhibitor MK571. The effect of PDMP on membrane transport correlated with modulation of calcium homeostasis, which occurred in a similar concentration range. PDMP released calcium from at least two independent calcium stores and blocked calcium influx induced by either extracellular ATP or thapsigargin. Thus, the biological effects of PDMP revealed a relation between three important physiological processes of multidrug resistance, calcium homeostasis, and membrane flow in the ER/ Golgi system.

A link between eumelanism and calcium physiology in the barn owl.

In many animals, melanin-based coloration is strongly heritable and is largely insensitive to the environment and body condition. According to the handicap principle, such a trait may not reveal individual quality because the production of different melanin-based colorations often entails similar costs. However, a recent study showed that the production of eumelanin pigments requires relatively large amounts of calcium, potentially implying that melanin-based coloration is associated with physiological processes requiring calcium. If this is the case, eumelanism may be traded-off against other metabolic processes that require the same elements. We used a correlative approach to examine, for the first time, this proposition in the barn owl, a species in which individuals vary in the amount, size, and blackness of eumelanic spots. For this purpose, we measured calcium concentration in the left humerus of 85 dead owls. Results showed that the humeri of heavily spotted individuals had a higher concentration of calcium. This suggests either that plumage spottiness signals the ability to absorb calcium from the diet for both eumelanin production and storage in bones, or that lightly spotted individuals use more calcium for metabolic processes at the expense of calcium storage in bones. Our study supports the idea that eumelanin-based coloration is associated with a number of physiological processes requiring calcium.

Numerical analysis of wave-induced fluid flow effects on seismic data: Application to monitoring of CO2 storage at the Sleipner Field

In this work we analyze how patchy distributions of CO2 and brine within sand reservoirs may lead to significant attenuation and velocity dispersion effects, which in turn may have a profound impact on surface seismic data. The ultimate goal of this paper is to contribute to the understanding of these processes within the framework of the seismic monitoring of CO2 sequestration, a key strategy to mitigate global warming. We first carry out a Monte Carlo analysis to study the statistical behavior of attenuation and velocity dispersion of compressional waves traveling through rocks with properties similar to those at the Utsira Sand, Sleipner field, containing quasi-fractal patchy distributions of CO2 and brine. These results show that the mean patch size and CO2 saturation play key roles in the observed wave-induced fluid flow effects. The latter can be remarkably important when CO2 concentrations are low and mean patch sizes are relatively large. To analyze these effects on the corresponding surface seismic data, we perform numerical simulations of wave propagation considering reservoir models and CO2 accumulation patterns similar to the CO2 injection site in the Sleipner field. These numerical experiments suggest that wave-induced fluid flow effects may produce changes in the reservoir's seismic response, modifying significantly the main seismic attributes usually employed in the characterization of these environments. Consequently, the determination of the nature of the fluid distributions as well as the proper modeling of the seismic data constitute important aspects that should not be ignored in the seismic monitoring of CO2 sequestration problems.

Calcium homeostasis in the central nervous system: adaptation to neurodegeneration.

Here we review the results of our recent studies on neurodegeneration together with data on cerebral calcium precipitation in animal models and humans. A model that integrates the diversity of mechanisms involved in neurodegeneration is presented and discussed based on the functional relevance of calcium precipitation.

Accurate calculations of intermolecular interaction energies using explicitly correlated coupled cluster wave functions and a dispersion-weighted MP2 method.

Explicitly correlated coupled-cluster calculations of intermolecular interaction energies for the S22 benchmark set of Jurecka, Sponer, Cerny, and Hobza (Chem. Phys. Phys. Chem. 2006, 8, 1985) are presented. Results obtained with the recently proposed CCSD(T)-F12a method and augmented double-zeta basis sets are found to be in very close agreement with basis set extrapolated conventional CCSD(T) results. Furthermore, we propose a dispersion-weighted MP2 (DW-MP2) approximation that combines the good accuracy of MP2 for complexes with predominately electrostatic bonding and SCS-MP2 for dispersion-dominated ones. The MP2-F12 and SCS-MP2-F12 correlation energies are weighted by a switching function that depends on the relative HF and correlation contributions to the interaction energy. For the S22 set, this yields a mean absolute deviation of 0.2 kcal/mol from the CCSD(T)-F12a results. The method, which allows obtaining accurate results at low cost, is also tested for a number of dimers that are not in the training set.

The cervical spine in calcium pyrophosphate dihydrate deposition disease. A prevalent case-control study.

OBJECTIVE: To test the hypothesis that calcium pyrophosphate dihydrate (CPPD) deposition disease is a risk factor for neck pain. METHODS: A prevalent case-control study was conducted to assess cervical calcifications and neck pain between patients with and without known peripheral CPPD deposition disease. CPPD cases were included if diagnosed with CPPD deposition disease of peripheral joints, and excluded if their chief complaint was neck pain. Controls were randomly selected among consecutive patients, hospitalized for conditions unrelated to CPPD deposition disease or neck pain, and matched to CPPD cases by age and sex. Cervical calcifications were assessed by lateral cervical radiographs and computed tomography scans of the upper cervical spine; neck pain and cervical function were appraised by a validated questionnaire. RESULTS: Cervical calcifications were found in 24 out of 35 patients (69%) in the CPPD group compared to 4 out of 35 patients (11%) in the control group (p < 0.001). Patients with CPPD deposition disease reported significantly more neck pain and discomfort than controls (p < 0.001), and were 5 times more likely to report any neck pain (odds ratio 5.5; 95% confidence interval: 1.9, 21.9). Among male patients, more extensive cervical calcified deposits correlated with more severe neck pain (rs = 0.58, p = 0.03). CONCLUSION: These results suggest that CPPD deposition disease frequently involves the cervical spine and may be associated with the development of neck pain.

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.

Calcium reabsorption in the distal tubule: regulation by sodium, pH, and flow.

We developed a mathematical model of Ca transport along the late distal convoluted tubule (DCT2) and the connecting tubule (CNT) to investigate the mechanisms that regulate Ca reabsorption in the DCT2-CNT. The model accounts for apical Ca influx across transient receptor potential vanilloid 5 (TRPV5) channels and basolateral Ca efflux via plasma membrane Ca-ATPase pumps and type 1 Na/Ca exchangers (NCX1). Model simulations reproduce experimentally observed variations in Ca uptake as a function of extracellular pH, Na, and Mg concentration. Our results indicate that amiloride enhances Ca reabsorption in the DCT2-CNT predominantly by increasing the driving force across NCX1, thereby stimulating Ca efflux. They also suggest that because aldosterone upregulates both apical and basolateral Na transport pathways, it has a lesser impact on Ca reabsorption than amiloride. Conversely, the model predicts that full NCX1 inhibition and parathyroidectomy each augment the Ca load delivered to the collecting duct severalfold. In addition, our results suggest that regulation of TRPV5 activity by luminal pH has a small impact, per se, on transepithelial Ca fluxes; the reduction in Ca reabsorption induced by metabolic acidosis likely stems from decreases in TRPV5 expression. In contrast, elevations in luminal Ca are predicted to significantly decrease TRPV5 activity via the Ca-sensing receptor. Nevertheless, following the administration of furosemide, the calcium-sensing receptor-mediated increase in Ca reabsorption in the DCT2-CNT is calculated to be insufficient to prevent hypercalciuria. Altogether, our model predicts complex interactions between calcium and sodium reabsorption in the DCT2-CNT.

Wave-vector dependence of spin and density multipole excitations in quantum dots

We have employed time-dependent local-spin density-functional theory to analyze the multipole spin and charge density excitations in GaAs-AlxGa1-xAs quantum dots. The on-plane transferred momentum degree of freedom has been taken into account, and the wave-vector dependence of the excitations is discussed. In agreement with previous experiments, we have found that the energies of these modes do not depend on the transferred wave vector, although their intensities do. Comparison with a recent resonant Raman scattering experiment [C. Schüller et al., Phys. Rev. Lett. 80, 2673 (1998)] is made. This allows us to identify the angular momentum of several of the observed modes as well as to reproduce their energies

Different PTH response to oral peptone load and oral calcium load in patients with normocalcemic primary hyperparathyroidism, primary hyperparathyroidism, and healthy subjects.

BACKGROUND: Normocalcemic primary hyperparathyroidism (PHPT-N) is a condition that may have similar long-term implications to primary hyperparathyroidism (PHPT); however, differential diagnosis and treatment for parathyroid disorders are not clearly defined. We investigated the effect of an oral peptone and an oral calcium load on calcium-regulating hormones in PHPT-N compared with PHPT and healthy controls to provide a new potential diagnostic tool. DESIGN: Case-control study. METHODS: We evaluated serum gastrin, PTH, ionized calcium, and phosphate responses to oral calcium (1 g) and peptone (10 g) load in 22 PHPT and 20 PHPT-N patients matched for PTH serum values. Moreover, 30 healthy subjects were enrolled as controls. In 12 patients for each group, we also performed the oral peptone test adding aluminum hydroxide (AH) to suppress phosphate absorption. RESULTS: In PHPT patients, PTH increased significantly 30 min after the oral peptone load, while no significant increase was found in PHPT-N and controls. After oral calcium load, PTH remained stable in PHPT while it decreased dramatically in PHPT-N patients, and ionized calcium increased significantly in each of the three groups. Peptones plus AH induced a blunted PTH increase in the three groups. CONCLUSIONS: Considering the marked difference in PTH response elicited by peptones in PHPT compared with PHPT-N, we suggest that the oral peptone test could be added to the diagnostic evaluation of PHPT patients. In case of absent response to peptones, patients should have their serum calcium levels assessed twice a year in accordance with recent guidelines.

Patterns of calcium-binding proteins in human inferior colliculus: identification of subdivisions and evidence for putative parallel systems.

The subdivisions of human inferior colliculus are currently based on Golgi and Nissl-stained preparations. We have investigated the distribution of calcium-binding protein immunoreactivity in the human inferior colliculus and found complementary or mutually exclusive localisations of parvalbumin versus calbindin D-28k and calretinin staining. The central nucleus of the inferior colliculus but not the surrounding regions contained parvalbumin-positive neuronal somata and fibres. Calbindin-positive neurons and fibres were concentrated in the dorsal aspect of the central nucleus and in structures surrounding it: the dorsal cortex, the lateral lemniscus, the ventrolateral nucleus, and the intercollicular region. In the dorsal cortex, labelling of calbindin and calretinin revealed four distinct layers.Thus, calcium-binding protein reactivity reveals in the human inferior colliculus distinct neuronal populations that are anatomically segregated. The different calcium-binding protein-defined subdivisions may belong to parallel auditory pathways that were previously demonstrated in non-human primates, and they may constitute a first indication of parallel processing in human subcortical auditory structures.

Complexometric titration with potenciometric indicator to determination of calcium and magnesium in soil extracts¹

This study proposes a method of direct and simultaneous determination of the amount of Ca2+ and Mg2+ present in soil extracts using a Calcium Ion-Selective Electrode and by Complexometric Titration (ISE-CT). The results were compared to those obtained by conventional analytical techniques of Complexometric Titration (CT) and Flame Atomic Absorption Spectrometry (FAAS). There were no significant differences in the determination of Ca2+ and Mg2+ in comparison with CT and FAAS, at a 95 % confidence level. Additionally, results of this method were more precise and accurate than of the Interlaboratorial Control (IC).