Examining the information content of time-lapse crosshole GPR data collected under different infiltration conditions to estimate unsaturated soil hydraulic properties

Time-lapse geophysical data acquired during transient hydrological experiments are being increasingly employed to estimate subsurface hydraulic properties at the field scale. In particular, crosshole ground-penetrating radar (GPR) data, collected while water infiltrates into the subsurface either by natural or artificial means, have been demonstrated in a number of studies to contain valuable information concerning the hydraulic properties of the unsaturated zone. Previous work in this domain has considered a variety of infiltration conditions and different amounts of time-lapse GPR data in the estimation procedure. However, the particular benefits and drawbacks of these different strategies as well as the impact of a variety of key and common assumptions remain unclear. Using a Bayesian Markov-chain-Monte-Carlo stochastic inversion methodology, we examine in this paper the information content of time-lapse zero-offset-profile (ZOP) GPR traveltime data, collected under three different infiltration conditions, for the estimation of van Genuchten-Mualem (VGM) parameters in a layered subsurface medium. Specifically, we systematically analyze synthetic and field GPR data acquired under natural loading and two rates of forced infiltration, and we consider the value of incorporating different amounts of time-lapse measurements into the estimation procedure. Our results confirm that, for all infiltration scenarios considered, the ZOP GPR traveltime data contain important information about subsurface hydraulic properties as a function of depth, with forced infiltration offering the greatest potential for VGM parameter refinement because of the higher stressing of the hydrological system. Considering greater amounts of time-lapse data in the inversion procedure is also found to help refine VGM parameter estimates. Quite importantly, however, inconsistencies observed in the field results point to the strong possibility that posterior uncertainties are being influenced by model structural errors, which in turn underlines the fundamental importance of a systematic analysis of such errors in future related studies.

The effect of recalled previous work environment on return to work after a rehabilitation program including vocational aspects for trauma patients.

INTRODUCTION: The aim of the present study was to assess the association between remembered previous work place environment and return to work (RTW) after hospitalisation in a rehabilitation hospital. METHODS: A cohort of 291 orthopedic trauma patients discharged from hospital between 15 December 2004 and 31 December 2005 was included in a study addressing quality of life and work-related questions. Remembered previous work environment was measured by Karasek's 31-item Job Content Questionnaire (JCQ), given to the patients during hospitalisation. Post-hospitalisation work status was assessed 3 months, 1, and 2 years after discharge, using a questionnaire sent to the ex-patients. Logistic regression models were used to test the role of four JCQ variables on RTW at each time point while controlling for relevant confounders. RESULTS: Subjects perceiving a higher physical demand were less likely to return to work 1 year after hospital discharge. Social support at work was positively associated with RTW at all time points. A high job strain appeared to be positively associated with RTW 1 year after rehabilitation, with limitations due to large confidence intervals. CONCLUSIONS: Perceptions of previous work environment may influence the probability of RTW. In a rehabilitation setting, efforts should be made to assess those perceptions and, if needed, interventions to modify them should be applied.

Health at work and risk factors related to ergonomics: some results from the Swiss sample of the 4th European working conditions survey

Introduction The European Foundation for the improvement of living and working conditions conducts a survey every 5 years since 1990. The foundation also offers the possibility to non-EU countries to be included in the survey: in 2005, Switzerland took part for the first time in the fourth edition of this survey. The Institute for Work and Health (IST) has been associated to the Swiss project conducted under the leadership of the SECO and the Fachhochschule Nordwestschweiz. The survey covers different aspects of work like job characteristics and employment conditions, health and safety, work organization, learning and development opportunities, and the balance between working and non-working life (Parent-Thirion, Fernandez Macias, Hurley, & Vermeylen, 2007). More particularly, one question assesses the worker's self-perception of the effects of work on health. We identified (for the Swiss sample) several factors affecting the risk to report health problems caused by work. The Swiss sample includes 1040 respondents. Selection of participants was based on a random multi-stage sampling and was carried out by M.I.S Trend S.A. (Lausanne). Participation rate was 59%. The database was weighted by household size, gender, age, region of domicile, occupational group, and economic sector. Specially trained interviewers carried out the interviews at the respondents home. The survey was carriedout between the 19th of September 2005 and the 30th of November 2005. As detailed in (Graf et al., 2007), 31% of the Swiss respondents identify work as the cause of health problems they experience. Most frequently reported health problems include back pain (18%), stress (17%), muscle pain (13%), and overall fatigue (11%). Ergonomic aspects associated with higher risk of reporting health problems caused by work include frequent awkward postures (odds ratio [OR] 4.7, 95% confidence interval [CI] 3.1 to 5.4), tasks involving lifting heavy loads (OR 2.7, 95% CI 2.0 to 3.6) or lifting people (OR 2.2, 95% CI 1.4 to 3.5), standing or walking (OR 1.4, 95% CI 1.1 to 1.9), as well as repetitive movements (OR 1.7, 95% CI 1.3 to 2.3). These results highlight the need to continue and intensify the prevention of work related health problems in occupations characterized by risk factors related to ergonomics.

Low dispersing species : population genetic processes in space and time : the genus "Trochulus" (Gastropoda: Hygromiidae) as a case study

RÉSUMÉ Une espèce est rarement composée d'une population unique. Parce que les individus ont des capacités de dispersion limitées et que les paysages sont des mosaïques d'habitats, la plupart des espèces sont plutôt composées de sous-populations connectées par la migration. Cette variation spatiale influence directement la distribution de la variabilité génétique dans et entre les populations. Durant ce travail, nous avons abordé certains des processus populationnels qui ont joué un rôle supposé dans l'apparition de nouvelles espèces au sein du genre Trochulus. Plus précisément, nous avons tenté d'évaluer les impacts respectifs de l'isolement passé (facteurs historiques) et présent (facteurs locaux). Nous avons d'abord pu montrer que les faibles capacités de dispersion des escargots terrestres ont directement influencé leur histoire évolutive à toutes les échelles spatiales et temporelles. En réduisant l'effet homogénéisant de la migration, une faible dispersion maintient dans les populations les traces génétiques d'évènements passés. A l'échelle de la distribution globale de Trochulus villosus, ces traces ont permis de reconstruire une histoire faite d'isolements et d'expansions de populations. En combinant des données génétiques avec une modélisation de la niche climatique passée, il a été possible de proposer un scénario significativement meilleur que toutes les hypothèses alternatives que nous avons testées. A l'échelle locale par contre, l'héritage historique est difficile à distinguer de la dynamique actuelle. Ce fut le cas des lignées mitochondriales du complexe sericeus-hispidus : les deux principales lignées étaient phylogénétiquement éloignées, avaient eu des démographies passées différentes et corrélaient avec des différences morphologiques. D'un autre côté, le flux de gène nucléaire était fort, contredisant l'idée de deux espèces cryptiques isolées reproductivement. Pour pouvoir conclure à la présence ou non de deux espèces, il nous a manqué des informations locales sur la dynamique des populations et les conditions écologiques que l'on trouve dans la région d'étude. Enfin, nous avons pu souligner que la connectivité entre populations d'escargots est soumise à la qualité des habitats et à leur organisation spatiale. Les escargots sont dépendants d'un habitat et s'y adaptent, comme l'indiquent la présence de «poils » uniquement sur la coquille d'espèces vivant dans des habitats humides ou la corrélation entre morphologie et habitat au sein du complexe sericeus-hispidus. Logiquement donc, les escargots migrent préférentiellement au travers d'habitats favorables comme l'a montré la réduction de flux de gènes au travers des prairies chez T. villosus (une espèce forestière). De ces données, nous pouvons supposer que les populations d'escargots en particulier, et des espèces à faible dispersion en général, ont de fortes chances d'être affectées par les changements climatiques, avec de probables implications pour leurs histoires évolutives. SUMMARY : Species rarely consists in a single population. Because individuals have limited dispersal abilities, because landscapes are habitat patchworks, most species are made of several subpopulations connected by migration. This spatial variation has consequences on the distribution of genetic diversity within and between populations, creating a structure among the populations. During the present work, we investigated some of the population processes assumed to have played an important role on the speciation within the genus Trochulus. More specifically, we questioned the respective impacts of past (historical factors) or present (local factors) population isolations. We first could show that the poor dispersal abilities of land snails have had profound impacts on their evolutionary histories at all spatial and temporal scales. Low dispersal maintains a strong signature of past events in the populations by minimising the homogenising effects of geneflow. At the scale of Trochulus villosus global distribution, they allowed to retrieve the detailed history of this species population isolations and expansions. Combining a large genetic dataset with paleo-climatic niche modelling ended up with a historical scenario significantly better than all traditional alternatives we tested. At local scale on the contrary, past events become difficult to tease apart from ongoing processes. This was the case for the divergent mitochondria) lineages within the sericeus-hispidus complex: the two principal lineages appeared to be phylogenetically distant, to have experienced different demographic histories and to correlate with morphological differences. On the other hand, nuclear (present day) geneflow was high, contradicting the idea of two reproductively isolated cryptic species. Information on the local population dynamics and environmental conditions are lacking to be able to decide whether past isolation has indeed resulted here in new species. Finally, we emphasised the importance of the habitat types present in a landscape as well as their spatial organisation for the population connectivity of land snails. These species are tightly dependent on a habitat and adapt to it as shown by thé occurrence of hair-like structures only in species living in humid environments or by the correlation between shell morphology and habitat in the sericeus-hispidus complex. As a result, land snails preferentially migrate through favourable habitats: Trochulus villosus, a forest species, had its geneflow significantly reduced across meadows. From these data, we can hypothesise that the populations of land snails in particular and of low dispersing species in general are likely to be strongly affected by the ongoing climate changes, with potential major consequences on their evolutionary histories.

Echo-time independent signal modulations for strongly coupled systems in triple echo localization schemes: an extension of S-PRESS editing.

The double spin-echo point resolved spectroscopy sequence (PRESS) is a widely used method and standard in clinical MR spectroscopy. Existence of important J-modulations at constant echo times, depending on the temporal delays between the rf-pulses, have been demonstrated recently for strongly coupled spin systems and were exploited for difference editing, removing singlets from the spectrum (strong-coupling PRESS, S-PRESS). A drawback of this method for in vivo applications is that large signal modulations needed for difference editing occur only at relatively long echo times. In this work we demonstrate that, by simply adding a third refocusing pulse (3S-PRESS), difference editing becomes possible at substantially shorter echo times while, as applied to citrate, more favorable lineshapes can be obtained. For the example of an AB system an analytical description of the MR signal, obtained with this triple refocusing sequence (3S-PRESS), is provided.

Time-Lapse AFM Imaging of DNA Conformational Changes Induced by Daunorubicin.

Cancer is a major health issue that absorbs the attention of a large part of the biomedical research. Intercalating agents bind to DNA molecules and can inhibit their synthesis and transcription; thus, they are increasingly used as drugs to fight cancer. In this work, we show how atomic force microscopy in liquid can characterize, through time-lapse imaging, the dynamical influence of intercalating agents on the supercoiling of DNA, improving our understanding of the drug's effect.

Structural joint inversion of time-lapse crosshole ERT and GPR traveltime data

Time-lapse geophysical monitoring and inversion are valuable tools in hydrogeology for monitoring changes in the subsurface due to natural and forced (tracer) dynamics. However, the resulting models may suffer from insufficient resolution, which leads to underestimated variability and poor mass recovery. Structural joint inversion using cross-gradient constraints can provide higher-resolution models compared with individual inversions and we present the first application to time-lapse data. The results from a synthetic and field vadose zone water tracer injection experiment show that joint 3-D time-lapse inversion of crosshole electrical resistance tomography (ERT) and ground penetrating radar (GPR) traveltime data significantly improve the imaged characteristics of the point injected plume, such as lateral spreading and center of mass, as well as the overall consistency between models. The joint inversion method appears to work well for cases when one hydrological state variable (in this case moisture content) controls the time-lapse response of both geophysical methods. Citation: Doetsch, J., N. Linde, and A. Binley (2010), Structural joint inversion of time-lapse crosshole ERT and GPR traveltime data, Geophys. Res. Lett., 37, L24404, doi: 10.1029/2010GL045482.

Colonic movements in healthy subjects as monitored by a Magnet Tracking System.

The Magnet Tracking System (MTS) is a minimally-invasive technique of continuous evaluation of gastrointestinal motility. In this study, MTS was used to analyse colonic propulsive dynamics and compare the transit of a magnetic pill with that of standard radio-opaque markers. MTS monitors the progress in real time of a magnetic pill through the gut. Ten men and 10 women with regular daily bowel movements swallowed this pill and 10 radio-opaque markers at 8 pm. Five hours of recordings were conducted during 2 following mornings. Origin, direction, amplitude and velocity of movements were analysed relative to space-time plots of the pill trajectory. Abdominal radiographs were taken to compare the progress of both pill and markers. The magnetic pill lay idle for 90% of its sojourn in the colon; its total retrograde displacement accounted for only 20% of its overall movement. Analysis of these movements showed a bimodal distribution of velocities: around 1.5 and 50 cm min(-1), the latter being responsible for 2/3 of distance traversed. There were more movements overall and more mass movements in males. Net hourly forward progress was greater in the left than right colon, and greater in males. The position of the magnetic pill correlated well with the advancement of markers. MTS showed patterns and propulsion dynamics of colonic segments with as yet unmet precision. Detailed analysis of slow and fast patterns of colonic progress makes it possible to specify the motility of colonic segments, and any variability in gender. Such analysis opens up promising avenues in studies of motility disorders.

Adaptive frequency tracking of the baseline ECG identifies the site of atrial fibrillation termination by catheter ablation

Multiple organization indices have been used to predict the outcome of stepwise catheter ablation in long-standing persistent atrial fibrillation (AF), however with limited success. Our study aims at developinginnovative organization indices from baseline ECG (i.e. during the procedure, before ablation) in orderto identify the site of AF termination by catheter ablation. Seventeen consecutive male patients (age60 ± 5 years, AF duration 7 ± 5 years) underwent a stepwise catheter ablation. Chest lead V6 was placedin the back (V6b). QRST cancelation was performed from chest leads V1 to V6b. Using an innovativeadaptive harmonic frequency tracking, two measures of AF organization were computed to quantify theharmonics components of ECG activity: (1) the adaptive phase difference variance (APD) between theAF harmonic components as a measure of AF regularity, and (2) and adaptive organization index (AOI)evaluating the cyclicity of the AF oscillations. Both adaptive indices were compared to indices computedusing a time-invariant approach: (1) ECG AF cycle length (AFCL), (2) the spectrum based organizationindex (OI), and (3) the time-invariant phase difference TIPD. Long-standing persistent AF was terminatedinto sinus rhythm or atrial tachycardia in 13/17 patients during stepwise ablation, 11 during left atriumablation (left terminated patients - LT), 2 during the right atrium ablation (right terminated patients -RT), and 4 were non terminated (NT) and required electrical cardioversion. Our findings showed that LTpatients were best separated from RT/NT before ablation by the duration of sustained AF and by AOI onchest lead V1 and APD from the dorsal lead V6b as compared to ECG AFCL, OI and TIPD, respectively. Ourresults suggest that adaptive measures of AF organization computed before ablation perform better thantime-invariant based indices for identifying patients whose AF will terminate during ablation within theleft atrium. These findings are indicative of a higher baseline organization in these patients that could beused to select candidates for the termination of AF by stepwise catheter ablation.© 2013 Elsevier Ltd. All rights reserved.

Effects of the transition time between muscle-tendon stretch and shortening on mechanical efficiency.

The net mechanical efficiency of positive work (eta(pos)) has been shown to increase if it is immediately preceded by negative work. This phenomenon is explained by the storage of elastic energy during the negative phase and its release during the subsequent positive phase. If a transition time (T) takes place, the elastic energy is dissipated into heat. The aim of the present study was to investigate the relationship between eta(pos) and T, and to determine the minimal T required so that eta(pos) reached its minimal value. Seven healthy male subjects were tested during four series of lowering-raising of the body mass. In the first series (S (0)), the negative and positive phases were executed without any transition time. In the three other series, T was varied by a timer (0.12, 0.24 and 0.56 s for series S (1), S (2) and S (3), respectively). These exercises were performed on a force platform sensitive to vertical forces to measure the mechanical work and a gas analyser was used to determine the energy expenditure. The results indicated that eta(pos) was the highest (31.1%) for the series without any transition time (S (0)). The efficiencies observed with transition times (S (1), S (2) and S (3)) were 27.7, 26.0 and 23.8%, respectively, demonstrating that T plays an important role for mechanical efficiency. The investigation of the relationship between eta(pos) and T revealed that the minimal T required so that eta(pos) reached its minimal value is 0.59 s.

Work motivation and pro-social behavior in the delivery of public services : theoretical and empirical insights

This Ph.D. dissertation seeks to study the work motivation of employees in the delivery of public services. The questioning on work motivation in public services in not new but it becomes central for governments which are now facing unprecedented public debts. The objective of this research is twofold : First, we want to see if the work motivation of employees in public services is a continuum (intrinsic and extrinsic motivations cannot coexist) or a bi-dimensional construct (intrinsic and extrinsic motivations coexist simultaneously). The research in public administration literature has focused on the concept of public service motivation, and considered motivation to be uni-dimensional (Perry and Hondeghem 2008). However, no study has yet tackled both types of motivation, the intrinsic and extrinsic ones, in the same time. This dissertation proposes, in Part I, a theoretical assessment and an empirical test of a global work motivational structure, by using a self-constructed Swiss dataset with employees from three public services, the education sector, the security sector and the public administrative services sector. Our findings suggest that work motivation in public services in not uni-dimensional but bi-dimensional, the intrinsic and extrinsic motivations coexist simultaneously and can be positively correlated (Amabile et al. 1994). Our findings show that intrinsic motivation is as important as extrinsic motivation, thus, the assumption that employees in public services are less attracted by extrinsic rewards is not confirmed for this sample. Other important finding concerns the public service motivation concept, which, as theoretically predicted, represents the major motivational dimension of employees in the delivery of public services. Second, the theory of public service motivation makes the assumption that employees in public services engage in activities that go beyond their self-interest, but never uses this construct as a determinant for their pro-social behavior. In the same time, several studies (Gregg et al. 2011 and Georgellis et al. 2011) bring evidence about the pro-social behavior of employees in public services. However, they do not identify which type of motivation is at the origin of this behavior, they only make the assumption of an intrinsically motivated behavior. We analyze the pro-social behavior of employees in public services and use the public service motivation as determinant of their pro-social behavior. We add other determinants highlighted by the theory of pro-social behavior (Bénabou and Tirole 2006), by Le Grand (2003) and by fit theories (Besley and Ghatak 2005). We test these determinants on Part II and identify for each sector of activity the positive or the negative impact on pro-social behavior of Swiss employees. Contrary to expectations, we find, for this sample, that both intrinsic and extrinsic factors have a positive impact on pro-social behavior, no crowding-out effect is identified in this sample. We confirm the hypothesis of Le Grand (2003) about the positive impact of the opportunity cost on pro-social behavior. Our results suggest a mix of action-oriented altruism and out-put oriented altruism of employees in public services. These results are relevant when designing incentives schemes for employees in the delivery of public services.

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.

Thalamic connectivity in patients with essential tremor treated with MR imaging-guided focused ultrasound: in vivo fiber tracking by using diffusion-tensor MR imaging.

PURPOSE: To use diffusion-tensor (DT) magnetic resonance (MR) imaging in patients with essential tremor who were treated with transcranial MR imaging-guided focused ultrasound lesion inducement to identify the structural connectivity of the ventralis intermedius nucleus of the thalamus and determine how DT imaging changes correlated with tremor changes after lesion inducement. MATERIALS AND METHODS: With institutional review board approval, and with prospective informed consent, 15 patients with medication-refractory essential tremor were enrolled in a HIPAA-compliant pilot study and were treated with transcranial MR imaging-guided focused ultrasound surgery targeting the ventralis intermedius nucleus of the thalamus contralateral to their dominant hand. Fourteen patients were ultimately included. DT MR imaging studies at 3.0 T were performed preoperatively and 24 hours, 1 week, 1 month, and 3 months after the procedure. Fractional anisotropy (FA) maps were calculated from the DT imaging data sets for all time points in all patients. Voxels where FA consistently decreased over time were identified, and FA change in these voxels was correlated with clinical changes in tremor over the same period by using Pearson correlation. RESULTS: Ipsilateral brain structures that showed prespecified negative correlation values of FA over time of -0.5 or less included the pre- and postcentral subcortical white matter in the hand knob area; the region of the corticospinal tract in the centrum semiovale, in the posterior limb of the internal capsule, and in the cerebral peduncle; the thalamus; the region of the red nucleus; the location of the central tegmental tract; and the region of the inferior olive. The contralateral middle cerebellar peduncle and bilateral portions of the superior vermis also showed persistent decrease in FA over time. There was strong correlation between decrease in FA and clinical improvement in hand tremor 3 months after lesion inducement (P < .001). CONCLUSION: DT MR imaging after MR imaging-guided focused ultrasound thalamotomy depicts changes in specific brain structures. The magnitude of the DT imaging changes after thalamic lesion inducement correlates with the degree of clinical improvement in essential tremor.

Adaptive tracking of EEG oscillations.

Neuronal oscillations are an important aspect of EEG recordings. These oscillations are supposed to be involved in several cognitive mechanisms. For instance, oscillatory activity is considered a key component for the top-down control of perception. However, measuring this activity and its influence requires precise extraction of frequency components. This processing is not straightforward. Particularly, difficulties with extracting oscillations arise due to their time-varying characteristics. Moreover, when phase information is needed, it is of the utmost importance to extract narrow-band signals. This paper presents a novel method using adaptive filters for tracking and extracting these time-varying oscillations. This scheme is designed to maximize the oscillatory behavior at the output of the adaptive filter. It is then capable of tracking an oscillation and describing its temporal evolution even during low amplitude time segments. Moreover, this method can be extended in order to track several oscillations simultaneously and to use multiple signals. These two extensions are particularly relevant in the framework of EEG data processing, where oscillations are active at the same time in different frequency bands and signals are recorded with multiple sensors. The presented tracking scheme is first tested with synthetic signals in order to highlight its capabilities. Then it is applied to data recorded during a visual shape discrimination experiment for assessing its usefulness during EEG processing and in detecting functionally relevant changes. This method is an interesting additional processing step for providing alternative information compared to classical time-frequency analyses and for improving the detection and analysis of cross-frequency couplings.

Improvement of cardiovascular risk prediction: time to review current knowledge, debates, and fundamentals on how to assess test characteristics.

Cardiovascular risk assessment might be improved with the addition of emerging, new tests derived from atherosclerosis imaging, laboratory tests or functional tests. This article reviews relative risk, odds ratios, receiver-operating curves, posttest risk calculations based on likelihood ratios, the net reclassification improvement and integrated discrimination. This serves to determine whether a new test has an added clinical value on top of conventional risk testing and how this can be verified statistically. Two clinically meaningful examples serve to illustrate novel approaches. This work serves as a review and basic work for the development of new guidelines on cardiovascular risk prediction, taking into account emerging tests, to be proposed by members of the 'Taskforce on Vascular Risk Prediction' under the auspices of the Working Group 'Swiss Atherosclerosis' of the Swiss Society of Cardiology in the future.