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.

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.

The red island and the seven dwarfs: body size reduction in Cheirogaleidae

AimSmall body size in Madagascar's dwarf and mouse lemurs (Cheirogaleidae) is generally viewed as primitive. We investigated the evolution of body size in this family and in its sister-taxon, the Lepilemuridae, from phylogenetic, ontogenetic and adaptive perspectives. LocationMadagascar. MethodsWe used a phylogenetic method to reconstruct the evolution of body size in lemurs, and allometric regression models of gestation periods and static and growth allometries in Cheirogaleidae and Lepilemuridae to test the hypothesis that dwarfing occurred as a result of truncated ontogeny (progenesis). We also examined adaptive hypotheses relating body size to environmental variability, life history, seasonality of reproduction, hypothermy (use of torpor), and a diet rich in plant exudates. ResultsOur results indicated that cheirogaleids experienced at least four independent events of body size reduction from an ancestor as large as living Lepilemuridae, by means of progenesis. Our interpretation is supported by the paedomorphic appearance and parallel ontogenetic trajectories of the dwarf taxa, as well as their very short gestation periods and increased fecundity. Lepilemur species that occupy more predictable environments are significantly larger than those occupying unpredictable habitats. Main conclusionsCheirogaleidae appear to be paedomorphic dwarfs, a consequence of progenesis, probably as an adaptation to high environmental unpredictability. Although the capacity to use hypothermy is related to small body size, this advantage is unlikely to have driven dwarfing in cheirogaleids. We propose that gummmivory/exudativory co-evolved with body size reduction in this clade, probably from a folivorous ancestor. Their small size is derived, and their suitability as models for the ancestral primate' is therefore dubious.

Reduction of symmetry in grandite solid solution

The results of a crystal structure refinement of an anisotropic grandite garnet specimen with composition Gro36-4 And63-6 are given. The structure obtained has orthorrombic symmetry (space group Fddd) and is compared with similar results obtained by other authors. In all cases the reduction of symmetry is due to the ordering of Fe3+ and Al in octahedral sites. Non cubic structures of grandites are discussed in connection with optical, morphological an grou-th features of these minerals.

Cation distribution and magnetization of BaFe12-2xCoxSnxO19 (x=0.9,1.28) single crystals

The distribution of Sn4+ cations within the five crystallographic sites of the magnetoplumbite (M) ‐like compound BaFe12−2xCoxSnxO19 has been analyzed using single‐crystal x‐ray‐diffraction data. The species Fe3+ and Co2+ cannot be distinguished using x rays because of their very similar atomic numbers; however, the calculation of the apparent valencies for the different sites allows an insight into the Co2+ cation segregation. The use of previous data from neutron powder diffraction allows a precise picture of the cation distribution, which indicates a pronounced site selectivity for both Sn4+ and Co2+ cations. The Sn4+ cations prefer the 4f2 sites and to a much lower extent the 12k sites, while they do not enter the octahedral 2a sites at all. Co2+ cations are distributed among tetrahedral and octahedral sites displaying a clear preference for the tetrahedral 4f1 sites. Magnetic measurements indicate that the compound still exhibits uniaxial anisotropy with the easy direction parallel to the c axis. Nevertheless, the magnetic structure shows a considerable degree of noncolinearity. A strong reduction of the magnetic anisotropy regarding that of the undoped compound is also detected.

Pharmacokinetic and pharmacodynamic analysis of efavirenz dose reduction using an in vitro-in vivo extrapolation model.

The pharmacokinetics (PK) of efavirenz (EFV) is characterized by marked interpatient variability that correlates with its pharmacodynamics (PD). In vitro-in vivo extrapolation (IVIVE) is a "bottom-up" approach that combines drug data with system information to predict PK and PD. The aim of this study was to simulate EFV PK and PD after dose reductions. At the standard dose, the simulated probability was 80% for viral suppression and 28% for central nervous system (CNS) toxicity. After a dose reduction to 400 mg, the probabilities of viral suppression were reduced to 69, 75, and 82%, and those of CNS toxicity were 21, 24, and 29% for the 516 GG, 516 GT, and 516 TT genotypes, respectively. With reduction of the dose to 200 mg, the probabilities of viral suppression decreased to 54, 62, and 72% and those of CNS toxicity decreased to 13, 18, and 20% for the 516 GG, 516 GT, and 516 TT genotypes, respectively. These findings indicate how dose reductions might be applied in patients with favorable genetic characteristics.

A high proportion of users of low-threshold facilities with needle exchange programmes in Switzerland are currently on methadone treatment: implications for new approaches in harm reduction and care.

BACKGROUND: Increasingly, patients receiving methadone treatment are found in low threshold facilities (LTF), which provide needle exchange programmes in Switzerland. This paper identifies the characteristics of LTF attendees receiving methadone treatment (MT) compared with other LTF attendees (non-MT). METHODS: A national cross-sectional survey was conducted in 2006 over five consecutive days in all LTF (n=25). Attendees were given an anonymous questionnaire, collecting information on socio-demographic indicators, drug consumption, injection, methadone treatment, and self-reported HIV and HCV status. Univariate analysis and logistic regression were performed to compare MT to non-MT. The response rate was 66% (n=1128). RESULTS: MT comprised 57.6% of the sample. In multivariate analysis, factors associated with being on MT were older age (OR: 1.38), being female (OR: 1.60), having one's own accommodation (OR: 1.56), receiving public assistance (OR: 2.29), lifetime injecting (OR: 2.26), HIV-positive status (OR: 2.00), and having consumed cocaine during the past month (OR: 1.37); MT were less likely to have consumed heroin in the past month (OR: 0.76, not significant) and visited LTF less often on a daily basis (OR: 0.59). The number of injections during the past week was not associated with MT. CONCLUSIONS: More LTF attendees were in the MT group, bringing to light an underappreciated LTF clientele with specific needs. The MT group consumption profile may reflect therapeutic failure or deficits in treatment quality and it is necessary to acknowledge this and to strengthen the awareness of LTF personnel about potential needs of MT attendees to meet their therapeutic goals.

Rotação de culturas no sistema plantio direto em Tibagi (PR): II - Emissões de CO2 e N2O

A atividade agrícola pode alterar a quantidade e qualidade da matéria orgânica do solo (MOS), resultando em emissões de dióxido de carbono (CO2) e óxido nitroso (N2O) do solo para a atmosfera. O sistema plantio direto (SPD) com a utilização de leguminosas em sistemas de rotação é uma estratégia que deve ser considerada tanto para o aumento da quantidade de MOS como para seu efeito na redução das emissões dos gases de efeito estufa. Com o objetivo de determinar os fluxos de gases do efeito estufa (CO2 e N2O) do solo, um experimento foi instalado em Tibagi (PR), em um Latossolo Vermelho distroférrico textura argilosa. Os tratamentos, dispostos em faixas não casualizadas com parcelas subdivididas, foram: sistema plantio direto por 12 anos com sucessões milho/trigo e soja/trigo (PD12 M/T e PD12 S/T, respectivamente) e por 22 anos (PD22 M/T e PD22 S/T, respectivamente). As emissões de CO2 do solo foram aproximadamente 20 % mais elevadas no PD22 em relação ao PD12. As emissões de CO2 apresentaram correlação significativa (R² = 0,85; p < 0,05) com a temperatura do solo, com emissões médias 40 % menores, registradas nos meses com temperaturas mais baixas. As emissões mais elevadas de N2O foram observadas após a colheita das culturas de verão, sobretudo na sucessão milho/trigo, em relação à sucessão soja/trigo. As emissões de N2O foram aproximadamente 25 % maiores após aplicação do fertilizante nitrogenado na cultura do trigo nas duas sucessões e apresentaram correlação significativa (R² = 0,88; p < 0,01) com o grau de saturação de água no solo (Sr %).

Improving Reconstituted HDL Composition for Efficient Post-Ischemic Reduction of Ischemia Reperfusion Injury.

BACKGROUND: New evidence shows that high density lipoproteins (HDL) have protective effects beyond their role in reverse cholesterol transport. Reconstituted HDL (rHDL) offer an attractive means of clinically exploiting these novel effects including cardioprotection against ischemia reperfusion injury (IRI). However, basic rHDL composition is limited to apolipoprotein AI (apoAI) and phospholipids; addition of bioactive compound may enhance its beneficial effects. OBJECTIVE: The aim of this study was to investigate the role of rHDL in post-ischemic model, and to analyze the potential impact of sphingosine-1-phosphate (S1P) in rHDL formulations. METHODS AND RESULTS: The impact of HDL on IRI was investigated using complementary in vivo, ex vivo and in vitro IRI models. Acute post-ischemic treatment with native HDL significantly reduced infarct size and cell death in the ex vivo, isolated heart (Langendorff) model and the in vivo model (-48%, p<0.01). Treatment with rHDL of basic formulation (apoAI + phospholipids) had a non-significant impact on cell death in vitro and on the infarct size ex vivo and in vivo. In contrast, rHDL containing S1P had a highly significant, protective influence ex vivo, and in vivo (-50%, p<0.01). This impact was comparable with the effects observed with native HDL. Pro-survival signaling proteins, Akt, STAT3 and ERK1/2 were similarly activated by HDL and rHDL containing S1P both in vitro (isolated cardiomyocytes) and in vivo. CONCLUSION: HDL afford protection against IRI in a clinically relevant model (post-ischemia). rHDL is significantly protective if supplemented with S1P. The protective impact of HDL appears to target directly the cardiomyocyte.

CO2, CH4 and N2O fluxes in an Ultisol treated with sewage sludge and cultivated with castor bean

Organic residue application into soil alter the emission of gases to atmosphere and CO2, CH4, N2O may contribute to increase the greenhouse effect. This experiment was carried out in a restoration area on a dystrophic Ultisol (PVAd) to quantify greenhouse gas (GHG) emissions from soil under castor bean cultivation, treated with sewage sludge (SS) or mineral fertilizer. The following treatments were tested: control without N; FertMin = mineral fertilizer; SS5 = 5 t ha-1 SS (37.5 kg ha-1 N); SS10 = 10 t ha-1 SS (75 kg ha-1 N); and SS20 = 20 t ha-1 SS (150 kg ha-1 N). The amount of sludge was based on the recommended rate of N for castor bean (75 kg ha-1), the N level of SS and the mineralization fraction of N from SS. Soil gas emission was measured for 21 days. Sewage sludge and mineral fertilizers altered the CO2, CH4 and N2O fluxes. Soil moisture had no effect on GHG emissions and the gas fluxes was statistically equivalent after the application of FertMin and of 5 t ha-1 SS. The application of the entire crop N requirement in the form of SS practically doubled the Global Warming Potential (GWP) and the C equivalent emissions in comparison with FertMin treatments.

Embodied pollution in Spanish household consumption: a disaggregate analysis

[spa] En este artículo aplicamos un modelo input-output ampliado medioambientalmente para analizar un aspecto específico de la hipótesis de la curva de Kuznets ambiental. El propósito del estudio es analizar si las estructuras de consumo de los hogares con una mejor ‘posición económica’ pueden tener un efecto positivo para reducir las presiones medioambientales. Para ello combinamos información de diferentes bases de datos para analizar el impacto de la contaminación atmosférica del consumo de diferentes hogares españoles en el año 2000. Consideramos nueve gases, i.e. los seis gases de efecto invernadero (CO2, CH4, N2O, SF6, HFCs, y PFCs) y otros tres gases (SO2, NOx, y NH3). Clasificamos los hogares en quintiles de gasto per capita y quintiles de gasto equivalente. Los resultados obtenidos muestran que hay una relación positiva y elevada entre el nivel de gasto y las emisiones directas e indirectas generadas por el consumo de los hogares; sin embargo, las intensidades de emisión tienden a disminuir con el nivel de gasto para los diferentes gases, con la excepción de SF6, HFCs, y PFCs.