Self-potential investigations of a gravel bar in a restored river corridor

Self-potentials (SP) are sensitive to water fluxes and concentration gradients in both saturated and unsaturated geological media, but quantitative interpretations of SP field data may often be hindered by the superposition of different source contributions and time-varying electrode potentials. Self-potential mapping and close to two months of SP monitoring on a gravel bar were performed to investigate the origins of SP signals at a restored river section of the Thur River in northeastern Switzerland. The SP mapping and subsequent inversion of the data indicate that the SP sources are mainly located in the upper few meters in regions of soil cover rather than bare gravel. Wavelet analyses of the time-series indicate a strong, but non-linear influence of water table and water content variations, as well as rainfall intensity on the recorded SP signals. Modeling of the SP response with respect to an increase in the water table elevation and precipitation indicate that the distribution of soil properties in the vadose zone has a very strong influence. We conclude that the observed SP responses on the gravel bar are more complicated than previously proposed semi-empiric relationships between SP signals and hydraulic head or the thickness of the vadose zone. We suggest that future SP monitoring in restored river corridors should either focus on quantifying vadose zone processes by installing vertical profiles of closely spaced SP electrodes or by installing the electrodes within the river to avoid signals arising from vadose zone processes and time-varying electrochemical conditions in the vicinity of the electrodes.

Body surface area and body weight predict total liver volume in Western adults.

Computed tomography (CT) is used increasingly to measure liver volume in patients undergoing evaluation for transplantation or resection. This study is designed to determine a formula predicting total liver volume (TLV) based on body surface area (BSA) or body weight in Western adults. TLV was measured in 292 patients from four Western centers. Liver volumes were calculated from helical computed tomographic scans obtained for conditions unrelated to the hepatobiliary system. BSA was calculated based on height and weight. Each center used a different established method of three-dimensional volume reconstruction. Using regression analysis, measurements were compared, and formulas correlating BSA or body weight to TLV were established. A linear regression formula to estimate TLV based on BSA was obtained: TLV = -794.41 + 1,267.28 x BSA (square meters; r(2) = 0.46; P <.0001). A formula based on patient weight also was derived: TLV = 191.80 + 18.51 x weight (kilograms; r(2) = 0.49; P <.0001). The newly derived TLV formula based on BSA was compared with previously reported formulas. The application of a formula obtained from healthy Japanese individuals underestimated TLV. Two formulas derived from autopsy data for Western populations were similar to the newly derived BSA formula, with a slight overestimation of TLV. In conclusion, hepatic three-dimensional volume reconstruction based on helical CT predicts TLV based on BSA or body weight. The new formulas derived from this correlation should contribute to the estimation of TLV before liver transplantation or major hepatic resection.

The prelude of the end-Permian mass extinction predates a postulated bolide impact

The mass extinction at the Permian-Triassic Boundary (PTB) is said to have been abrupt and probably caused by an extraterrestrial impact. However, evidence from the Global Stratotype Section and Point (GSSP) of the base of the Induan at Meishan, China, shows that the biotic crisis began prior to the level, in beds 25 and 26 at which the postulated impact event occurred. Evidence of such an earlier biotic crisis occurs in other sections in South China, and in central and western Tethyan regions. This event is characterized by the extinction of a range of faunas, including corals, deep-water radiolarians, most fusulinids and pseudotirolitidammonoids, and many Permian brachiopods. In all sections, this extinction level is usually a few decimeters to meters below that of the main mass extinction in the event beds (25 and 26) at Meishan, and their correlatives elsewhere. This earlier extinction event happened before the postulated bolide impact at the level of beds 25 and 26, and constrains interpretation of the mechanisms that brought about this greatest mass extinction.

Numerical simulation of ambient seismic wavefield modification caused by pore-fluid effects in an oil reservoir

We have modeled numerically the seismic response of a poroelastic inclusion with properties applicable to an oil reservoir that interacts with an ambient wavefield. The model includes wave-induced fluid flow caused by pressure differences between mesoscopic-scale (i.e., in the order of centimeters to meters) heterogeneities. We used a viscoelastic approximation on the macroscopic scale to implement the attenuation and dispersion resulting from this mesoscopic-scale theory in numerical simulations of wave propagation on the kilometer scale. This upscaling method includes finite-element modeling of wave-induced fluid flow to determine effective seismic properties of the poroelastic media, such as attenuation of P- and S-waves. The fitted, equivalent, viscoelastic behavior is implemented in finite-difference wave propagation simulations. With this two-stage process, we model numerically the quasi-poroelastic wave-propagation on the kilometer scale and study the impact of fluid properties and fluid saturation on the modeled seismic amplitudes. In particular, we addressed the question of whether poroelastic effects within an oil reservoir may be a plausible explanation for low-frequency ambient wavefield modifications observed at oil fields in recent years. Our results indicate that ambient wavefield modification is expected to occur for oil reservoirs exhibiting high attenuation. Whether or not such modifications can be detected in surface recordings, however, will depend on acquisition design and noise mitigation processing as well as site-specific conditions, such as the geologic complexity of the subsurface, the nature of the ambient wavefield, and the amount of surface noise.

Rate and processes of river network rearrangement during incipient faulting: The case of the Cahabon river, Guatemala

Deeply incised river networks are generally regarded as robust features that are not easily modified by erosion or tectonics. Although the reorganization of deeply incised drainage systems has been documented, the corresponding importance with regard to the overall landscape evolution of mountain ranges and the factors that permit such reorganizations are poorly understood. To address this problem, we have explored the rapid drainage reorganization that affected the Cahabon River in Guatemala during the Quaternary. Sediment-provenance analysis, field mapping, and electrical resistivity tomography (ERT) imaging are used to reconstruct the geometry of the valley before the river was captured. Dating of the abandoned valley sediments by the Be-10-Al-26 burial method and geomagnetic polarity analysis allow us to determine the age of the capture events and then to quantify several processes, such as the rate of tectonic deformation of the paleovalley, the rate of propagation of post-capture drainage reversal, and the rate at which canyons that formed at the capture sites have propagated along the paleovalley. Transtensional faulting started 1 to 3 million years ago, produced ground tilting and ground faulting along the Cahabon River, and thus generated differential uplift rate of 0.3 +/- 0.1 up to 0.7 +/- 0.4 mm . y(-1) along the river's course. The river responded to faulting by incising the areas of relative uplift and depositing a few tens of meters of sediment above the areas of relative subsidence. Then, the river experienced two captures and one avulsion between 700 ky and 100 ky. The captures breached high-standing ridges that separate the Cahabon River from its captors. Captures occurred at specific points where ridges are made permeable by fault damage zones and/or soluble rocks. Groundwater flow from the Cahabon River down to its captors likely increased the erosive power of the captors thus promoting focused erosion of the ridges. Valley-fill formation and capture occurred in close temporal succession, suggesting a genetic link between the two. We suggest that the aquifers accumulated within the valley-fills, increased the head along the subterraneous system connecting the Cahabon River to its captors, and promoted their development. Upon capture, the breached valley experienced widespread drainage reversal toward the capture sites. We attribute the generalized reversal to combined effects of groundwater sapping in the valley-fill, axial drainage obstruction by lateral fans, and tectonic tilting. Drainage reversal increased the size of the captured areas by a factor of 4 to 6. At the capture sites, 500 m deep canyons have been incised into the bedrock and are propagating upstream at a rate of 3 to 11 mm . y(-1) deepening at a rate of 0.7 to 1 5 mm . y(-1). At this rate, 1 to 2 million years will be necessary for headward erosion to completely erase the topographic expression of the paleovalley. It is concluded that the rapid reorganization of this drainage system was made possible by the way the river adjusted to the new tectonic strain field, which involved transient sedimentation along the river's course. If the river had escaped its early reorganization and had been given the time necessary to reach a new dynamic equilibrium, then the transient conditions that promoted capture would have vanished and its vulnerability to capture would have been strongly reduced.

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.

Localization and characterization of an active fault in an urbanized area in central Guatemala by means of geoelectrical imaging

The Polochic and Motagua faults define the active plate boundary between the North American and Caribbean plates in central Guatemala. A splay of the Polochic Fault traverses the rapidly growing city of San Miguel Uspantan that is periodically affected by destructive earthquakes. This fault splay was located using a 2D electrical resistivity tomography (ERT) survey that also characterized the fault damage zone and evaluated the thickness and nature of recent deposits upon which most of the city is built. ERT images show the fault as a similar to 50 m wide, near-vertical low-resistivity anomaly, bounded within a few meters by high resistivity anomalies. Forward modeling reproduces the key aspects of the observed electrical resistivity data with remarkable fidelity thus defining the overall location, geometry, and internal structure of the fault zone as well as the affected lithologies. Our results indicate that the city is constructed on a similar to 20 m thick surficial layer consisting of poorly consolidated, highly porous, water-logged pumice. This soft layer is likely to amplify seismic waves and to liquefy upon moderate to strong ground shaking. The electrical conductivity as well as the major element chemistry of the groundwater provides evidence to suggest that the local aquifer might, at least in part, be fed by water rising along the fault. Therefore, the potential threat posed by this fault splay may not be limited to its seismic activity per se, but could be compounded its potential propensity to enhance seismic site effects by injecting water into the soft surficial sediments. The results of this study provide the basis for a rigorous analysis of seismic hazard and sustainable development of San Miguel Uspantan and illustrate the potential of ERT surveying for paleoseismic studies.

Changes over time in risk factors for cardiovascular disease and use of lipid-lowering drugs in HIV-infected individuals and impact on myocardial infarction.

BACKGROUND: Because of the known relationship between exposure to combination antiretroviral therapy and cardiovascular disease (CVD), it has become increasingly important to intervene against risk of CVD in human immunodeficiency virus (HIV)-infected patients. We evaluated changes in risk factors for CVD and the use of lipid-lowering therapy in HIV-infected individuals and assessed the impact of any changes on the incidence of myocardial infarction. METHODS: The Data Collection on Adverse Events of Anti-HIV Drugs Study is a collaboration of 11 cohorts of HIV-infected patients that included follow-up for 33,389 HIV-infected patients from December 1999 through February 2006. RESULTS: The proportion of patients at high risk of CVD increased from 35.3% during 1999-2000 to 41.3% during 2005-2006. Of 28,985 patients, 2801 (9.7%) initiated lipid-lowering therapy; initiation of lipid-lowering therapy was more common for those with abnormal lipid values and those with traditional risk factors for CVD (male sex, older age, higher body mass index [calculated as the weight in kilograms divided by the square of the height in meters], family and personal history of CVD, and diabetes mellitus). After controlling for these, use of lipid-lowering drugs became relatively less common over time. The incidence of myocardial infarction (0.32 cases per 100 person-years [PY]; 95% confidence interval [CI], 0.29-0.35 cases per 100 PY) appeared to remain stable. However, after controlling for changes in risk factors for CVD, the rate decreased over time (relative rate in 2003 [compared with 1999-2000], 0.73 cases per 100 PY [95% CI, 0.50-1.05 cases per 100 PY]; in 2004, 0.64 cases per 100 PY [95% CI, 0.44-0.94 cases per 100 PY]; in 2005-2006, 0.36 cases per 100 PY [95% CI, 0.24-0.56 cases per 100 PY]). Further adjustment for lipid levels attenuated the relative rates towards unity (relative rate in 2003 [compared with 1999-2000], 1.06 cases per 100 PY [95% CI, 0.63-1.77 cases per 100 PY]; in 2004, 1.02 cases per 100 PY [95% CI, 0.61-1.71 cases per 100 PY]; in 2005-2006, 0.63 cases per 100 PY [95% CI, 0.36-1.09 cases per 100 PY]). CONCLUSIONS: Although the CVD risk profile among patients in the Data Collection on Adverse Events of Anti-HIV Drugs Study has decreased since 1999, rates have remained relatively stable, possibly as a result of a more aggressive approach towards managing the risk of CVD.

Ambient seismic noise monitoring of a clay landslide: Toward failure prediction

Given that clay-rich landslides may become mobilized, leading to rapid mass movements (earthflows and debris flows), they pose critical problems in risk management worldwide. The most widely proposed mechanism leading to such flow-like movements is the increase in water pore pressure in the sliding mass, generating partial or complete liquefaction. This solid-to-liquid transition results in a dramatic reduction of mechanical rigidity in the liquefied zones, which could be detected by monitoring shear wave velocity variations. With this purpose in mind, the ambient seismic noise correlation technique has been applied to measure the variation in the seismic surface wave velocity in the Pont Bourquin landslide (Swiss Alps). This small but active composite earthslide-earthflow was equipped with continuously recording seismic sensors during spring and summer 2010. An earthslide of a few thousand cubic meters was triggered in mid-August 2010, after a rainy period. This article shows that the seismic velocity of the sliding material, measured from daily noise correlograms, decreased continuously and rapidly for several days prior to the catastrophic event. From a spectral analysis of the velocity decrease, it was possible to determine the location of the change at the base of the sliding layer. These results demonstrate that ambient seismic noise can be used to detect rigidity variations before failure and could potentially be used to predict landslides.

Lateglacial evolution of the Brenno glacier (Southern Swiss Alps) : preliminary results

The Lateglacial evolution of the Ticino glacier and tributaries is poorly known because of the lack of research by Quaternary geomorphologists during the last decades. In spite of the interest for the cryosphere reactions during the Lateglacial climate warming, only few scientific studies were carried out about the history of the northern valleys of the Ticino Alps during the deglaciation (e.g. Seiffert 1953, Renner 1982, Hantke 1983). Within the framework of geomorphological investigations on the Lateglacial and Holocene glacier/permafrost evolution in the Ticino Alps, the history of the Brenno glacier (Blenio Valley, Eastern Ticino Alps) during the end of the Pleistocene has been studied. The deglaciation sequence of the Blenio Valley is still not complete (Scapozza et al. 2009). Only the first glacial stadial of the Brenno glacier and the last Lateglacial stadials of the Greina region (northern Blenio valley, see Fontana et al. 2008) and of the upper Malvaglia Valley (eastern Blenio Valley, see Scapozza et al. 2008) have been unequivocally defined. For every stadial, the surface of the palaeoglacier and the depression of the Equilibrium Line Altitude (ELA) have been reconstructed on the base of geomorphological mapping. The first individual glacial stadial of the Brenno glacier corresponds to the Biasca stadial of the Ticino glacier defined by Hantke (1983). The ELA depression of 1100-1200 meters and its morphological and glaciological characteristics allow us to correlate this stadial with the Weissbad stadial defined by Keller (1988). In the Greina region, three stadials corresponding to the end of the Lateglacial have been identified, with an ELA depression of 110, 210 and 310-350 meters (Fontana et al. 2008). In the upper Malvaglia Valley, three stadials corresponding to the end of the Oldest Dryas and the Younger Dryas have been identified for the Orino glacier, with an ELA depression of 290, 400-420 and 470-560 meters (Scapozza et al. 2008). If we consider the other (fragmentary) glacial deposits of the Blenio Valley, it is possible to define a regression sequence of the Brenno glacier with 8 stadials, from the Biasca stadial to the end of the Younger Dryas. An attempt of correlation with the model "Gothard" developed by Renner (1982) and Hantke (1983) and with the model "Eastern Swiss Alps" developed by Maisch (1982) is proposed in Table 1. The following chronological conclusions are, therefore, proposed: (1) the Biasca stadial is probably the first stadial after the transition Pleniglacial - Lateglacial; (2) the stadials BRE 7 to BRE 3 are positioned between the beginning of the Lateglacial and the Bølling-Allerød interstadial; (3) the stadials BRE 2 and BRE 1 are assumed to be related to the Younger Dryas event.

Heated Optical Fiber for Distributed Soil-Moisture Measurements: A Lysimeter Experiment

An Actively Heated Fiber Optics (AHFO) method to estimate soil moisture is tested and the analysis technique improved on. The measurements were performed in a lysimeter uniformly packed with loam soil with variable water content profiles. In the first meter of the soil profi le, 30 m of fiber optic cable were installed in a 12 loops coil. The metal sheath armoring the fiber cable was used as an electrical resistance heater to generate a heat pulse, and the soil response was monitored with a Distributed Temperature Sensing (DTS) system. We study the cooling following three continuous heat pulses of 120 s at 36 W m(-1) by means of long-time approximation of radial heat conduction. The soil volumetric water contents were then inferred from the estimated thermal conductivities through a specifically calibrated model relating thermal conductivity and volumetric water content. To use the pre-asymptotic data we employed a time correction that allowed the volumetric water content to be estimated with a precision of 0.01-0.035 (m(3) m(-3)). A comparison of the AHFO measurements with soil-moisture measurements obtained with calibrated capacitance-based probes gave good agreement for wetter soils [discrepancy between the two methods was less than 0.04 (m(3) m(-3))]. In the shallow drier soils, the AHFO method underestimated the volumetric water content due to the longertime required for the temperature increment to become asymptotic in less thermally conductive media [discrepancy between the two methods was larger than 0.1 (m(3) m(-3))]. The present work suggests that future applications of the AHFO method should include longer heat pulses, that longer heating and cooling events are analyzed, and, temperature increments ideally be measured with higher frequency.

Étude des sédiments quaternaires, de la molasse et sa tectonique, dans le Grand Lac (Léman) à partir de données sismiques 2D et 3D

RESUME L'Institut de Géophysique de l'Université de Lausanne a développé au cours de ces dernières années un système d'acquisition de sismique réflexion multitrace à haute résolution 2D et 3D. L'objectif de cette thèse était de poursuivre ce développement tout améliorant les connaissances de la géologie sous le lac Léman, en étudiant en particulier la configuration des grands accidents sous-lacustres dans la Molasse (Tertiaire) qui forme l'essentiel du substratum des formations quaternaires. En configuration 2D, notre système permet d'acquérir des profils sismiques avec une distance inter-CDP de 1,25 m. La couverture varie entre 6 et 18 selon le nombre de traces et la distance inter-tir. Le canon à air (15/15 eu. in.), offre une résolution verticale de 1,25 ni et une pénétration maximale de 300 m sous le fond de l'eau. Nous avons acquis au total plus de 400 km de sections 2D dans le Grand Lac et le Haut Lac entre octobre 2000 et juillet 2004. Une campagne de sismique 3D a fourni des données au large d'Evian sur une surface de 442,5 m sur 1450 m, soit 0,64 km2. La navigation ainsi que le positionnement des hydrophones et de la source ont été réalisés avec des GPS différentiels. Nous avons utilisé un traitement sismique conventionnel, sans appliquer d'AGC et en utilisant une migration post-stack. L'interprétation du substratum antéquaternaire est basée sur l'identification des sismofaciès, sur leurs relations avec les unités géologiques adjacentes au lac, ainsi que sur quelques données de forages. Nous obtenons ainsi une carte des unités géologiques dans le Grand Lac. Nous précisons la position du chevauchement subalpin entre la ville de Lausanne, sur la rive nord, et le bassin de Sciez, sur la rive sud. Dans la Molasse de Plateau, nous avons identifié les décrochements de Pontarlier et de St. Cergue ainsi que plusieurs failles non reconnues jusqu'ici. Nous avons cartographié les accidents qui affectent la Molasse subalpine ainsi que le plan de chevauchement du flysch sur la Molasse près de la rive sud du lac. Une nouvelle carte tectonique de la région lémanique a ainsi pu être dressée. L'analyse du substratum ne montre pas de failles suggérant une origine tectonique de la cuvette lémanique. Par contre, nous suggérons que la forme du creusement glaciaire, donc de la forme du lac Léman, a été influencée par la présence de failles dans le substratum antéquaternaire. L'analyse des sédiments quaternaires nous a permis de tracer des cartes des différentes interfaces ou unités qui les composent. La carte du toit du substratum antéquaternaire montre la présence de chenaux d'origine glaciaire dont la profondeur maximale atteint la cote -200 ni. Leur pente est dirigée vers le nord-est, à l'inverse du sens d'écoulement actuel des eaux. Nous expliquons cette observation par l'existence de circulations sous-glaciaires d'eau artésienne. Les sédiments glaciaires dont l'épaisseur maximale atteint 150 ni au centre du lac ont enregistré les différentes récurrences glaciaires. Dans la zone d'Evian, nous mettons en évidence la présence de lentilles de sédiments glaciolacustres perchées sur le flanc de la cuvette lémanique. Nous avons corrélé ces unités avec des données de forage et concluons qu'il s'agit du complexe inférieur de la pile sédimentaire d'Evian. Celui-ci, âgé de plus de 30 000 ans, serait un dépôt de Kame associé à un lac périglaciaire. La sismique réflexion 3D permet de préciser l'orientation de l'alimentation en matériel détritique de l'unité. La finesse des images obtenues nous permet également d'établir quels types d'érosion ont affecté certaines unités. Les sédiments lacustres, dont l'épaisseur maximale imagée atteint plus de 225 m et sans doute 400 ni sous le delta du Rhône, indiquent plusieurs mécanismes de dépôts. A la base, une mégaturbidite, épaisse d'une trentaine de mètres en moyenne, s'étend entre l'embouchure de la Dranse et le delta du Rhône. Au-dessus, la décantation des particules en suspension d'origine biologique et détritique fournit l'essentiel des sédiments. Dans la partie orientale du lac, les apports détritiques du Rhône forment un delta qui prograde vers l'ouest en s'imbriquant avec les sédiments déposés par décantation. La structure superficielle du delta a brutalement évolué, probablement à la suite de l'évènement catastrophique du Tauredunum (563 A.D.). Sa trace probable se marque par la présence d'une surface érosive que nous avons cartographiée. Le delta a ensuite changé de géométrie, avec notamment un déplacement des chenaux sous-lacustres. Sur l'ensemble de nos sections sismiques, nous n'observons aucune faille dans les sédiments quaternaires qui attesterait d'une tectonique postglaciaire du substratum. ABSTRACT During the last few years the institute of Geophysics of the University of Lausanne cleveloped a 2D and 3D high-resolution multichannel seismic reflection acquisition system. The objective of the present work was to carry on this development white improving our knowledge of the geology under Lake Geneva, in particular by studying the configuration of the large accidents affecting the Tertiary Molasse that makes up the basement of most Quaternary deposits. In its 2D configuration, our system makes it possible to acquire seismic profiles with a CDP interval of 1.25 m. The fold varies from 6 to 18 depending on the number of traces and the shooting interval. Our air gun (15/15 cu. in.) provides a vertical resolution of 1.25 m and a maximum penetration depth of approximately 300 m under water bottom. We acquired more than 400 km of 2D sections in the Grand Lac and the Haut Lac between October 2000 and July 2004. A 3D seismic survey off the city of Evian provided data on a surface of 442.5 m x 1450 m (0.64 km2). Ship's navigation as well as hydrophone- and source positioning were carried out with differential GPS. The seismic data were processed following a conventional sequence without .applying AGC and using post-stack migration. The interpretation of the pre-Quaternary substratum is based on sismofacies, on their relationships with terrestrial geological units and on some borehole data. We thus obtained a map of the geological units in the Grand Lac. We defined the location of the subalpine thrust from Lausanne, on the north shore, to the Sciez Basin, on the south shore. Within the Molasse de Plateau, we identified the already know Pontarlier and St Cergue transforms Fault as well as faults. We mapped faults that affect subalpine Molasse as well as the thrust fault plane between alpine flysch and Molasse near the lake's south shore. A new tectonic map of the Lake Geneva region could thus be drawn up. The substratum does not show faults indicating a tectonic origin for the Lake Geneva Basin. However, we suggest that the orientation of glacial erosion, and thus the shape of Lake Geneva, vas influenced by the presence of faults in the pre-Quaternary basement. The analysis of Quaternary sediments enabled us to draw up maps of various discontinuities or internal units. The top pre-Quaternary basement map shows channels of glacial origin, the deepest of them reaching an altitude of 200 m a.s.l. The channel's slopes are directed to the North-East, in opposite direction of the present water flow. We explain this observation by the presence of artesian subglacial water circulation. Glacial sediments, the maximum thickness of which reaches 150 m in the central part of the lake, record several glacial recurrences. In the Evian area, we found lenses of glacio-lacustrine sediments set high up on the flank of the Lake Geneva Bassin. We correlated these units with on-land borehole data and concluded that they represent the lower complex of the Evian sedimentary pile. The lower complex is aider than 30 000 years, and it could be a Kame deposit associated with a periglacial lake. Our 3D seismic reflexion survey enables us to specify the supply direction of detrital material in this unit. With detailed seismic images we established how some units were affected by different erosion types. The lacustrine sediments we imaged in Lake Geneva are thicker than 225 m and 400 m or more Linder the Rhone Delta. They indicate several depositional mechanisms. Their base is a major turbidite, thirty meters thick on average, that spreads between the Dranse mouth and the Rhone delta. Above this unit, settling of suspended biological and detrital particles provides most of the sediments. In the eastern part of the lake, detrital contribution from the Rhone builds a delta that progrades to the west and imbricates with the settling sediments. The shallow structure of the Rhone delta abruptly evolved, probably after the catastrophic Tauredunum event (563 A.D.). It probably coincides with an erosive surface that we mapped. As a result, the delta geometry changed, in particular associated with a displacement of water bottom channels. In all our seismic sections, we do not observe fault in the Quaternary sediments that would attest postglacial tectonic activity in the basement.

Outcome evaluation of ankle osteoarthritis treatments using spatio-temporal gait parameters and plantar pressure during unconstrained long distance walking : P25

Introduction: Ankle arthrodesis (AD) and total ankle replacement (TAR) are typical treatments for ankle osteoarthritis (AO). Despite clinical interest, there is a lack of their outcome evaluation using objective criteria. Gait analysis and plantar pressure assessment are appropriate to detect pathologies in orthopaedics but they are mostly used in lab with few gait cycles. In this study, we propose an ambulatory device based on inertial and plantar pressure sensors to compare the gait during long-distance trials between healthy subjects (H) and patients with AO or treated by AD and TAR. Methods: Our study included four groups: 11 patients with AO, 9 treated by TAR, 7 treated by AD and 6 control subjects. An ambulatory system (Physilog®, CH) was used for gait analysis; plantar pressure measurements were done using a portable insole (Pedar®-X, DE). The subjects were asked to walk 50 meters in two trials. Mean value and coefficient of variation of spatio-temporal gait parameters were calculated for each trial. Pressure distribution was analyzed in ten subregions of foot. All parameters were compared among the four groups using multi-level model-based statistical analysis. Results: Significant difference (p <0.05) with control was noticed for AO patients in maximum force in medial hindfoot and forefoot and in central forefoot. These differences were no longer significant in TAR and AD groups. Cadence and speed of all pathologic groups showed significant difference with control. Both treatments showed a significant improvement in double support and stance. TAR decreased variability in speed, stride length and knee ROM. Conclusions: In spite of a small sample size, this study showed that ankle function after AO treatments can be evaluated objectively based on plantar pressure and spatio-temporal gait parameters measured during unconstrained walking outside the lab. The combination of these two ambulatory techniques provides a promising way to evaluate foot function in clinics.

Fluid migration through thrust faults in the Helvetic nappes (Western Swiss Alps)

The delta(18)O, delta(13)C and Sr-87/Sr-86 values of calcite and organic matter were measured for carbonate mylonites from numerous thrusts in the Helvetic Alps. Carbonate mylonites in most of the thrusts retain essentially unaltered protolith delta(18)O and delta(13)C values, consistent with there having been little to no advection of isotopically distinct fluid through these faults. Only carbonate mylonites from the basal thrusts of the largest nappes have delta(18)O and/or delta(13)C values that differ from those of their protoliths. The zone of isotopic alteration/exchange is confined to c. 10 to 20 meters of these fault contacts. We propose the fluids that migrated through these faults contained variable amounts of organically derived carbon and radiogenic strontium, and were probably derived from dewatering of the sedimentary rocks and prograde metamorphic reactions in the nappes' root zones. Apart from the basal thrusts of the largest nappes that behaved as narrow, laterally extensive conduits for fluids, there is little isotopic evidence that large quantities of fluids passed through most of the carbonate-hosted thrusts in the Helvetic Alps.

Incidence and risk factors for chronic elevation of alanine aminotransferase levels in HIV-infected persons without hepatitis b or c virus co-infection.

BACKGROUND: Chronic liver disease in human immunodeficiency virus (HIV)-infected patients is mostly caused by hepatitis virus co-infection. Other reasons for chronic alanine aminotransferase (ALT) elevation are more difficult to diagnose. METHODS: We studied the incidence of and risk factors for chronic elevation of ALT levels (greater than the upper limit of normal at 2 consecutive semi-annual visits) in participants of the Swiss HIV Cohort Study without hepatitis B virus (HBV) or hepatitis C virus (HCV) infection who were seen during the period 2002-2008. Poisson regression analysis was used. RESULTS: A total of 2365 participants were followed up for 9972 person-years (median age, 38 years; male sex, 66%; median CD4+ cell count, 426/microL; receipt of antiretroviral therapy [ART], 56%). A total of 385 participants (16%) developed chronic elevated ALT levels, with an incidence of 3.9 cases per 100 person-years (95% confidence interval [CI], 3.5-4.3 cases per 100 person-years). In multivariable analysis, chronic elevated ALT levels were associated with HIV RNA level >100,000 copies/mL (incidence rate ratio [IRR], 2.23; 95% CI, 1.45-3.43), increased body mass index (BMI, defined as weight in kilograms divided by the square of height in meters) (BMI of 25-29.9 was associated with an IRR of 1.56 [95% CI, 1.24-1.96]; a BMI 30 was associated with an IRR of 1.70 [95% CI, 1.16-2.51]), severe alcohol use (1.83 [1.19-2.80]), exposure to stavudine (IRR per year exposure, 1.12 [95% CI, 1.07-1.17]) and zidovudine (IRR per years of exposure, 1.04 [95% CI, 1.00-1.08]). Associations with cumulative exposure to combination ART, nucleoside reverse-transcriptase inhibitors, and unboosted protease inhibitors did not remain statistically significant after adjustment for exposure to stavudine. Black ethnicity was inversely correlated (IRR, 0.52 [95% CI, 0.33-0.82]). Treatment outcome and mortality did not differ between groups with and groups without elevated ALT levels. CONCLUSIONS: Among patients without hepatitis virus co-infection, the incidence of chronic elevated ALT levels was 3.9 cases per 100 person-years, which was associated with high HIV RNA levels, increased BMI, severe alcohol use, and prolonged stavudine and zidovudine exposure. Long-term follow-up is needed to assess whether chronic elevation of ALT levels will result in increased morbidity or mortality.