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.

Aortic vessel wall magnetic resonance imaging at 3.0 Tesla: a reproducibility study of respiratory navigator gated free-breathing 3D black blood magnetic resonance imaging.

The purpose of this study was to evaluate a free-breathing three-dimensional (3D) dual inversion-recovery (DIR) segmented k-space gradient-echo (turbo field echo [TFE]) imaging sequence at 3T for the quantification of aortic vessel wall dimensions. The effect of respiratory motion suppression on image quality was tested. Furthermore, the reproducibility of the aortic vessel wall measurements was investigated. Seven healthy subjects underwent 3D DIR TFE imaging of the aortic vessel wall with and without respiratory navigator. Subsequently, this sequence with respiratory navigator was performed twice in 10 healthy subjects to test its reproducibility. The signal-to-noise (SNR), contrast-to-noise ratio (CNR), vessel wall sharpness, and vessel wall volume (VWV) were assessed. Data were compared using the paired t-test, and the reproducibility of VWV measurements was evaluated using intraclass correlation coefficients (ICCs). SNR, CNR, and vessel wall sharpness were superior in scans performed with respiratory navigator compared to scans performed without. The ICCs concerning intraobserver, interobserver, and interscan reproducibility were excellent (0.99, 0.94, and 0.95, respectively). In conclusion, respiratory motion suppression substantially improves image quality of 3D DIR TFE imaging of the aortic vessel wall at 3T. Furthermore, this optimized technique with respiratory motion suppression enables assessment of aortic vessel wall dimensions with high reproducibility.

Mass conservative three-dimensional water tracer distribution from Markov chain Monte Carlo inversion of time-lapse ground-penetrating radar data

Time-lapse geophysical measurements are widely used to monitor the movement of water and solutes through the subsurface. Yet commonly used deterministic least squares inversions typically suffer from relatively poor mass recovery, spread overestimation, and limited ability to appropriately estimate nonlinear model uncertainty. We describe herein a novel inversion methodology designed to reconstruct the three-dimensional distribution of a tracer anomaly from geophysical data and provide consistent uncertainty estimates using Markov chain Monte Carlo simulation. Posterior sampling is made tractable by using a lower-dimensional model space related both to the Legendre moments of the plume and to predefined morphological constraints. Benchmark results using cross-hole ground-penetrating radar travel times measurements during two synthetic water tracer application experiments involving increasingly complex plume geometries show that the proposed method not only conserves mass but also provides better estimates of plume morphology and posterior model uncertainty than deterministic inversion results.

Mapping cell-matrix stresses during stretch reveals inelastic reorganization of the cytoskeleton.

The mechanical properties of the living cell are intimately related to cell signaling biology through cytoskeletal tension. The tension borne by the cytoskeleton (CSK) is in part generated internally by the actomyosin machinery and externally by stretch. Here we studied how cytoskeletal tension is modified during stretch and the tensional changes undergone by the sites of cell-matrix interaction. To this end we developed a novel technique to map cell-matrix stresses during application of stretch. We found that cell-matrix stresses increased with imposition of stretch but dropped below baseline levels on stretch release. Inhibition of the actomyosin machinery resulted in a larger relative increase in CSK tension with stretch and in a smaller drop in tension after stretch release. Cell-matrix stress maps showed that the loci of cell adhesion initially bearing greater stress also exhibited larger drops in traction forces after stretch removal. Our results suggest that stretch partially disrupts the actin-myosin apparatus and the cytoskeletal structures that support the largest CSK tension. These findings indicate that cells use the mechanical energy injected by stretch to rapidly reorganize their structure and redistribute tension.

Molecular imaging of matrix metalloproteinases in atherosclerotic plaques.

Ischaemic stroke and myocardial infarction often result from the sudden rupture of an atherosclerotic plaque. The subsequent arterial thrombosis occluding the vessel lumen has been widely indicated as the crucial acute event causing peripheral tissue ischaemia. A complex cross-talk between systemic and intraplaque inflammatory mediators has been shown to regulate maturation, remodeling and final rupture of an atherosclerotic plaque. Matrix metalloproteinases (MMPs) are proteolytic enzymes (released by several cell subsets within atherosclerotic plaques), which favour atherogenesis and increase plaque vulnerability. Thus, the assessment of intraplaque levels and activity of MMP might be of pivotal relevance in the evaluation of the risk of rupture. New imaging approaches, focused on the visualisation of inflammation in the vessel wall and plaque, may emerge as tools for individualised risk assessment and prevention of events. In this review, we summarize experimental findings of the currently available invasive and noninvasive imaging techniques, used to detect the presence and activity of MMPs in atherosclerotic plaques.

Silver-assisted laser desorption ionization for high spatial resolution imaging mass spectrometry of olefins from thin tissue sections.

Silver has been demonstrated to be a powerful cationization agent in mass spectrometry (MS) for various olefinic species such as cholesterol and fatty acids. This work explores the utility of metallic silver sputtering on tissue sections for high resolution imaging mass spectrometry (IMS) of olefins by laser desorption ionization (LDI). For this purpose, sputtered silver coating thickness was optimized on an assorted selection of mouse and rat tissues including brain, kidney, liver, and testis. For mouse brain tissue section, the thickness was adjusted to 23 ± 2 nm of silver to prevent ion suppression effects associated with a higher cholesterol and lipid content. On all other tissues, a thickness of at 16 ± 2 nm provided the best desorption/ionization efficiency. Characterization of the species by MS/MS showed a wide variety of olefinic compounds allowing the IMS of different lipid classes including cholesterol, arachidonic acid, docosahexaenoic acid, and triacylglyceride 52:3. A range of spatial resolutions for IMS were investigated from 150 μm down to the high resolution cellular range at 5 μm. The applicability of direct on-tissue silver sputtering to LDI-IMS of cholesterol and other olefinic compounds presents a novel approach to improve the amount of information that can be obtained from tissue sections. This IMS strategy is thus of interest for providing new biological insights on the role of cholesterol and other olefins in physiological pathways or disease.

Using matrix attachment regions to improve recombinant protein production.

Chinese hamster ovary (CHO) cells are the system of choice for the production of complex molecules, such as monoclonal antibodies. Despite significant progress in improving the yield from these cells, the process to the selection, identification, and maintenance of high-producing cell lines remains cumbersome, time consuming, and often of uncertain outcome. Matrix attachment regions (MARs) are DNA sequences that help generate and maintain an open chromatin domain that is favourable to transcription and may also facilitate the integration of several copies of the transgene. By incorporating MARs into expression vectors, an increase in the proportion of high-producer cells as well as an increase in protein production are seen, thereby reducing the number of clones to be screened and time to production by as much as 9 months. In this chapter, we describe how MARs can be used to increase transgene expression and provide protocols for the transfection of CHO cells in suspension and detection of high-producing antibody cell clones.

Building a Social Accounting Matrix within the ESA95 Framework: Obtaining a Dataset for Applied General Equilibrium Modelling

This research provides a description of the process followed in order to assemble a "Social Accounting Matrix" for Spain corresponding to the year 2000 (SAMSP00). As argued in the paper, this process attempts to reconcile ESA95 conventions with requirements of applied general equilibrium modelling. Particularly, problems related to the level of aggregation of net taxation data, and to the valuation system used for expressing the monetary value of input-output transactions have deserved special attention. Since the adoption of ESA95 conventions, input-output transactions have been preferably valued at basic prices, which impose additional difficulties on modellers interested in computing applied general equilibrium models. This paper addresses these difficulties by developing a procedure that allows SAM-builders to change the valuation system of input-output transactions conveniently. In addition, this procedure produces new data related to net taxation information.

The mean-variance model from the inverse of the variance-covariance matrix

En este artículo, a partir de la inversa de la matriz de varianzas y covarianzas se obtiene el modelo Esperanza-Varianza de Markowitz siguiendo un camino más corto y matemáticamente riguroso. También se obtiene la ecuación de equilibrio del CAPM de Sharpe.

Matrix metalloproteinase 9 (MMP-9/gelatinase B) proteolytically cleaves ICAM-1 and participates in tumor cell resistance to natural killer cell-mediated cytotoxicity.

Shedding of intercellular adhesion molecule 1 (ICAM-1) is believed to play a role in tumor cell resistance to cell-mediated cytotoxicity. However, the mechanism whereby ICAM-1 is shed from the surface of tumor cells remains unclear. In this study, we have addressed the possibility that matrix metalloproteinases are implicated in ICAM-1 shedding. Our observations suggest a functional relationship between ICAM-1 and matrix metalloproteinase 9 (MMP-9) whereby ICAM-1 provides a cell surface docking mechanism for proMMP-9, which, upon activation, proteolytically cleaves the extracellular domain of ICAM-1 leading to its release from the cell surface. MMP-9-dependent shedding of ICAM-1 is found to augment tumor cell resistance to natural killer (NK) cell-mediated cytotoxicity. Taken together, our observations propose a mechanism for ICAM-1 shedding from the cell surface and provide support for MMP involvement in tumor cell evasion of immune surveillance.

Circulating matrix γ-carboxyglutamate protein (MGP) species are refractory to vitamin K treatment in a new case of Keutel syndrome.

Summary.  Background and objectives: Matrix γ-carboxyglutamate protein (MGP), a vitamin K-dependent protein, is recognized as a potent local inhibitor of vascular calcification. Studying patients with Keutel syndrome (KS), a rare autosomal recessive disorder resulting from MGP mutations, provides an opportunity to investigate the functions of MGP. The purpose of this study was (i) to investigate the phenotype and the underlying MGP mutation of a newly identified KS patient, and (ii) to investigate MGP species and the effect of vitamin K supplements in KS patients. Methods: The phenotype of a newly identified KS patient was characterized with specific attention to signs of vascular calcification. Genetic analysis of the MGP gene was performed. Circulating MGP species were quantified and the effect of vitamin K supplements on MGP carboxylation was studied. Finally, we performed immunohistochemical staining of tissues of the first KS patient originally described focusing on MGP species. Results: We describe a novel homozygous MGP mutation (c.61+1G>A) in a newly identified KS patient. No signs of arterial calcification were found, in contrast to findings in MGP knockout mice. This patient is the first in whom circulating MGP species have been characterized, showing a high level of phosphorylated MGP and a low level of carboxylated MGP. Contrary to expectations, vitamin K supplements did not improve the circulating carboxylated MGP levels. Phosphorylated MGP was also found to be present in the first KS patient originally described. Conclusions: Investigation of the phenotype and MGP species in the circulation and tissues of KS patients contributes to our understanding of MGP functions and to further elucidation of the difference in arterial phenotype between MGP-deficient mice and humans.

Security strategies and equilibria in multiobjetives matrix games

Multiobjective matrix games have been traditionally analyzed from two different points of view: equiibrium concepts and security strategies. This paper is based upon the idea that both players try to reach equilibrium points playing pairs of security strategies, as it happens in scalar matrix games. We show conditions guaranteeing the existence of equilibria in security strategies, named security equilibria

Semiclassical evaluation of average nuclear one- and two-body matrix elements

Thomas-Fermi theory is developed to evaluate nuclear matrix elements averaged on the energy shell, on the basis of independent particle Hamiltonians. One- and two-body matrix elements are compared with the quantal results, and it is demonstrated that the semiclassical matrix elements, as function of energy, well pass through the average of the scattered quantum values. For the one-body matrix elements it is shown how the Thomas-Fermi approach can be projected on good parity and also on good angular momentum. For the two-body case, the pairing matrix elements are considered explicitly.

Density matrix functional theory that includes pairing correlations

The extension of density functional theory (DFT) to include pairing correlations without formal violation of the particle-number conservation condition is described. This version of the theory can be considered as a foundation of the application of existing DFT plus pairing approaches to atoms, molecules, ultracooled and magnetically trapped atomic Fermi gases, and atomic nuclei where the number of particles is conserved exactly. The connection with Hartree-Fock-Bogoliubov (HFB) theory is discussed, and the method of quasilocal reduction of the nonlocal theory is also described. This quasilocal reduction allows equations of motion to be obtained which are much simpler for numerical solution than the equations corresponding to the nonlocal case. Our theory is applied to the study of some even Sn isotopes, and the results are compared with those obtained in the standard HFB theory and with the experimental ones.