Méthodes sismiques actives adaptées à l'étude de milieux hétérogènes : application aux glissements de terrain

Résumé Les glissements de terrain représentent un des principaux risques naturels dans les régions montagneuses. En Suisse, chaque année les glissements de terrains causent des dégâts qui affectent les infrastructures et ont des coûts financiers importants. Une bonne compréhension des mécanismes des glissements peut permettre d'atténuer leur impact. Celle-ci passe notamment par la connaissance de la structure interne du glissement, la détermination de son volume et de son ou ses plans de glissement. Dans un glissement de terrain, la désorganisation et la présence de fractures dans le matériel déplacé engendre un changement des paramètres physiques et en particulier une diminution des vitesses de propagation des ondes sismiques ainsi que de la densité du matériel. Les méthodes sismiques sont de ce fait bien adaptées à l'étude des glissements de terrain. Parmi les méthodes sismiques, l'analyse de la dispersion des ondes de surface est une méthode simple à mettre en oeuvre. Elle présente l'avantage d'estimer les variations des vitesses de cisaillement avec la profondeur sans avoir spécifiquement recours à l'utilisation d'une source d'onde S et de géophones horizontaux. Sa mise en oeuvre en trois étapes implique la mesure de la dispersion des ondes de surface sur des réseaux étendus, la détermination des courbes de dispersion pour finir par l'inversion de ces courbes. Les modèles de vitesse obtenus à partir de cette procédure ne sont valides que lorsque les milieux explorés ne présentent pas de variations latérales. En pratique cette hypothèse est rarement vérifiée, notamment pour un glissement de terrain dans lequel les couches remaniées sont susceptibles de présenter de fortes hétérogénéités latérales. Pour évaluer la possibilité de déterminer des courbes de dispersion à partir de réseaux de faible extension des mesures testes ont été effectuées sur un site (Arnex, VD) équipé d'un forage. Un profil sismique de 190 m de long a été implanté dans une vallée creusée dans du calcaire et remplie par des dépôts glacio-lacustres d'une trentaine de mètres d'épaisseur. Les données acquises le long de ce profil ont confirmé que la présence de variations latérales sous le réseau de géophones affecte l'allure des courbes de dispersion jusqu'à parfois empêcher leur détermination. Pour utiliser l'analyse de la dispersion des ondes de surface sur des sites présentant des variations latérales, notre approche consiste à déterminer les courbes de dispersions pour une série de réseaux de faible extension, à inverser chacune des courbes et à interpoler les différents modèles de vitesse obtenus. Le choix de la position ainsi que de l'extension des différents réseaux de géophones est important. Il tient compte de la localisation des hétérogénéités détectées à partir de l'analyse de sismique réfraction, mais également d'anomalies d'amplitudes observées sur des cartes qui représentent dans le domaine position de tir - position du récepteur, l'amplitude mesurée pour différentes fréquences. La procédure proposée par Lin et Lin (2007) s'est avérée être une méthode efficace permettant de déterminer des courbes de dispersion à partir de réseaux de faible extension. Elle consiste à construire à partir d'un réseau de géophones et de plusieurs positions de tir un enregistrement temps-déports qui tient compte d'une large gamme de distances source-récepteur. Au moment d'assembler les différentes données une correction de phase est appliquée pour tenir compte des hétérogénéités situées entre les différents points de tir. Pour évaluer cette correction nous suggérons de calculer pour deux tir successif la densité spectrale croisée des traces de même offset: Sur le site d'Arnex, 22 courbes de dispersions ont été déterminées pour de réseaux de géophones de 10 m d'extension. Nous avons également profité du forage pour acquérir un profil de sismique verticale en ondes S. Le modèle de vitesse S déduit de l'interprétation du profil de sismique verticale est utilisé comme information à priori lors l'inversion des différentes courbes de dispersion. Finalement, le modèle en deux dimension qui a été établi grâce à l'analyse de la dispersion des ondes de surface met en évidence une structure tabulaire à trois couches dont les limites coïncident bien avec les limites lithologiques observées dans le forage. Dans celui-ci des argiles limoneuses associées à une vitesse de propagation des ondes S de l'ordre de 175 m/s surmontent vers 9 m de profondeur des dépôts de moraine argilo-sableuse caractérisés par des vitesses de propagation des ondes S de l'ordre de 300 m/s jusqu'à 14 m de profondeur et supérieur ou égal à 400 m/s entre 14 et 20 m de profondeur. Le glissement de la Grande Combe (Ballaigues, VD) se produit à l'intérieur du remplissage quaternaire d'une combe creusée dans des calcaires Portlandien. Comme dans le cas du site d'Arnex les dépôts quaternaires correspondent à des dépôts glacio-lacustres. Dans la partie supérieure la surface de glissement a été localisée à une vingtaine de mètres de profondeur au niveau de l'interface qui sépare des dépôts de moraine jurassienne et des dépôts glacio-lacustres. Au pied du glissement 14 courbes de dispersions ont été déterminées sur des réseaux de 10 m d'extension le long d'un profil de 144 m. Les courbes obtenues sont discontinues et définies pour un domaine de fréquence de 7 à 35 Hz. Grâce à l'utilisation de distances source-récepteur entre 8 et 72 m, 2 à 4 modes de propagation ont été identifiés pour chacune des courbes. Lors de l'inversion des courbes de dispersion la prise en compte des différents modes de propagation a permis d'étendre la profondeur d'investigation jusqu'à une vingtaine de mètres de profondeur. Le modèle en deux dimensions permet de distinguer 4 couches (Vs1 < 175 m/s, 175 m/s < Vs2 < 225 m/s, 225 m/s < Vs3 < 400 m/s et Vs4 >.400 m/s) qui présentent des variations d'épaisseur. Des profils de sismiques réflexion en ondes S acquis avec une source construite dans le cadre de ce travail, complètent et corroborent le modèle établi à partir de l'analyse de la dispersion des ondes de surface. Un réflecteur localisé entre 5 et 10 m de profondeur et associé à une vitesse de sommation de 180 m/s souligne notamment la géométrie de l'interface qui sépare la deuxième de la troisième couche du modèle établi à partir de l'analyse de la dispersion des ondes de surface. Abstract Landslides are one of the main natural hazards in mountainous regions. In Switzerland, landslides cause damages every year that impact infrastructures and have important financial costs. In depth understanding of sliding mechanisms may help limiting their impact. In particular, this can be achieved through a better knowledge of the internal structure of the landslide, the determination of its volume and its sliding surface or surfaces In a landslide, the disorganization and the presence of fractures in the displaced material generate a change of the physical parameters and in particular a decrease of the seismic velocities and of the material density. Therefoe, seismic methods are well adapted to the study of landslides. Among seismic methods, surface-wave dispersion analysis is a easy to implement. Through it, shearwave velocity variations with depth can be estimated without having to resort to an S-wave source and to horizontal geophones. Its 3-step implementation implies measurement of surface-wave dispersion with long arrays, determination of the dispersion curves and finally inversion of these curves. Velocity models obtained through this approach are only valid when the investigated medium does not include lateral variations. In practice, this assumption is seldom correct, in particular for landslides in which reshaped layers likely include strong lateral heterogeneities. To assess the possibility of determining dispersion curves from short array lengths we carried out tests measurements on a site (Arnex, VD) that includes a borehole. A 190 m long seismic profile was acquired in a valley carved into limestone and filled with 30 m of glacio-lacustrine sediments. The data acquired along this profile confirmed that the presence of lateral variations under the geophone array influences the dispersion-curve shape so much that it sometimes preventes the dispersion curves determination. Our approach to use the analysis of surface-wave dispersion on sites that include lateral variations consists in obtaining dispersion curves for a series of short length arrays; inverting each so obtained curve and interpolating the different obtained velocity model. The choice of the location as well as the geophone array length is important. It takes into account the location of the heterogeneities that are revealed by the seismic refraction interpretation of the data but also, the location of signal amplitude anomalies observed on maps that represent, for a given frequency, the measured amplitude in the shot position - receiver position domain. The procedure proposed by Lin and Lin (2007) turned out to be an efficient one to determine dispersion curves using short extension arrays. It consists in building a time-offset from an array of geophones with a wide offset range by gathering seismograms acquired with different source-to-receiver offsets. When assembling the different data, a phase correction is applied in order to reduce static phase error induced by lateral variation. To evaluate this correction, we suggest to calculate, for two successive shots, the cross power spectral density of common offset traces. On the Arnex site, 22 curves were determined with 10m in length geophone-arrays. We also took advantage of the borehole to acquire a S-wave vertical seismic profile. The S-wave velocity depth model derived from the vertical seismic profile interpretation is used as prior information in the inversion of the dispersion-curves. Finally a 2D velocity model was established from the analysis of the different dispersion curves. It reveals a 3-layer structure in good agreement with the observed lithologies in the borehole. In it a clay layer with a shear-wave of 175 m/s shear-wave velocity overlies a clayey-sandy till layer at 9 m depth that is characterized down to 14 m by a 300 m/s S-wave velocity; these deposits have a S-wave velocity of 400 m/s between depths of 14 to 20 m. The La Grand Combe landslide (Ballaigues, VD) occurs inside the Quaternary filling of a valley carved into Portlandien limestone. As at the Arnex site, the Quaternary deposits correspond to glaciolacustrine sediments. In the upper part of the landslide, the sliding surface is located at a depth of about 20 m that coincides with the discontinuity between Jurassian till and glacio-lacustrine deposits. At the toe of the landslide, we defined 14 dispersion curves along a 144 m long profile using 10 m long geophone arrays. The obtained curves are discontinuous and defined within a frequency range of 7 to 35 Hz. The use of a wide range of offsets (from 8 to 72 m) enabled us to determine 2 to 4 mode of propagation for each dispersion curve. Taking these higher modes into consideration for dispersion curve inversion allowed us to reach an investigation depth of about 20 m. A four layer 2D model was derived (Vs1< 175 m/s, 175 m/s <Vs2< 225 m/s, 225 m/s < Vs3 < 400 m/s, Vs4> 400 m/s) with variable layer thicknesses. S-wave seismic reflection profiles acquired with a source built as part of this work complete and the velocity model revealed by surface-wave analysis. In particular, reflector at a depth of 5 to 10 m associated with a 180 m/s stacking velocity image the geometry of the discontinuity between the second and third layer of the model derived from the surface-wave dispersion analysis.

Earthquake impacts in old-growth Nothofagus forests in New Zealand

Six stands located on different land forms in mixed old-growth Nothofagus forests in the Matiri Valley (northwest of South Island. New Zealand) were sampled to examine the effects of two recent large earthquakes on tree establishment and tree-ring growth, and how these varied across land forms. 50 trees were cor ed in each stand to determine age structure and the cores were cross-dated to precisely date unusual periods of radial growth. The 1968 earthquake (M = 7.1, epicentre 35 km from the study area) had no discernible impact on the sampled stands. The impact of the 1929 earthquake (M = 7.7, epicentre 20 kin from the study area) varied between stands, depending on whether or not they had been damaged by soil or rock movement. In all stands, the age structures showed a pulse of N. fusca establishment following the 1929 earthquake, with this species dominating establishment in large gaps created by landslides. Smaller gaps, created by branch or tree death, were closed by both N. fusca and N. menziesii. The long period of releases (1929-1945) indicates that direct earthquake damage was not the only cause of tree death, and that many trees died subsequently most likely of pathogen attack or a drought in the early 1930s. The impacts of the 1929 earthquake are compared to a storm in 1905 and a drought in 1974-1978 which also affected forests in the region. Our results confirm that earthquakes are an important factor driving forest dynamics in this tectonically active region, and that the diversity of earthquake impacts is a major source of heterogeneity in forest structure and regeneration.

Intraoperative angiography reloaded: a new hybrid operating theater for combined endovascular and surgical treatment of cerebral arteriovenous malformations: a pilot study on 25 patients.

BACKGROUND: Multimodality treatment suites for patients with cerebral arteriovenous malformations (AVM) have recently become available. This study was designed to evaluate feasibility, safety and impact on treatment of a new intraoperative flat-panel (FP) based integrated surgical and imaging suite for combined endovascular and surgical treatment of cerebral AVM. METHODS: Twenty-five patients with AVMs to treat with combined endovascular and surgical interventions were prospectively enrolled in this consecutive case series. The hybrid suite allows combined endovascular and surgical approaches with intraoperative scanner-like imaging (XperCT®) and intraoperative 3D rotational angiography (3D-RA). The impact of intraoperative multimodal imaging on feasibility, workflow of combined interventions, surgery, and unexpected imaging findings were analyzed. RESULTS: Twenty-five patients (mean age 38 ± 18.6 year) with a median Spetzler-Martin grade 2 AVM (range 1-4) underwent combined endovascular and surgical procedures. Sixteen patients presented with a ruptured AVM and nine with an unruptured AVM. In 16 % (n = 4) of cases, intraoperative imaging visualized AVM remnants ≤3 mm and allowed for completion of the resections in the same sessions. Complete resection was confirmed in all n = 16 patients who had follow-up angiography one year after surgery so far. All diagnostic and therapeutical steps, including angiographic control, were performed without having to move the patients CONCLUSION: The hybrid neurointerventional suite was shown to be a safe and useful setup which allowed for unconstrained combined microsurgical and neuroradiological workflow. It reduces the need for extraoperative angiographic controls and subsequent potential surgical revisions a second time, as small AVM remnants can be detected with high security.

Diagnostic value of vestibulo-ocular reflex parameters in the detection and characterization of labyrinthine lesions.

OBJECTIVE: To evaluate the power of various parameters of the vestibulo-ocular reflex (VOR) in detecting unilateral peripheral vestibular dysfunction and in characterizing certain inner ear pathologies. STUDY DESIGN: Prospective study of consecutive ambulatory patients presenting with acute onset of peripheral vertigo and spontaneous nystagmus. SETTING: Tertiary referral center. PATIENTS: Seventy-four patients (40 females, 34 males) and 22 normal subjects (11 females, 11 males) were included in the study. Patients were classified in three main diagnoses: vestibular neuritis: 40; viral labyrinthitis: 22; Meniere's disease: 12. METHODS: The VOR function was evaluated by standard caloric and impulse rotary tests (velocity step). A mathematical model of vestibular function was used to characterize the VOR response to rotational stimulation. The diagnostic value of the different VOR parameters was assessed by uni- and multivariable logistic regression. RESULTS: In univariable analysis, caloric asymmetry emerged as the most powerful VOR parameter in identifying unilateral vestibular deficit, with a boundary limit set at 20%. In multivariable analysis, the combination of caloric asymmetry and rotational time constant asymmetry significantly improved the discriminatory power over caloric alone (p&lt;0.0001) and produced a detection score with a correct classification of 92.4%. In discriminating labyrinthine diseases, different combinations of the VOR parameters were obtained for each diagnosis (p&lt;0.003) supporting that the VOR characteristics differ between the three inner ear disorders. However, the clinical usefulness of these characteristics in separating the pathologies was limited. CONCLUSION: We propose a powerful logistic model combining the indices of caloric and time constant asymmetries to detect a peripheral vestibular loss, with an accuracy of 92.4%. Based on vestibular data only, the discrimination between the different inner ear diseases is statistically possible, which supports different pathophysiologic changes in labyrinthine pathologies.

Characterization of gravitational rock slope deformations at different spatial scales based on field, remote sensing and numerical approaches

Rock slope instabilities such as rock slides, rock avalanche or deep-seated gravitational slope deformations are widespread in Alpine valleys. These phenomena represent at the same time a main factor that control the mountain belts erosion and also a significant natural hazard that creates important losses to the mountain communities. However, the potential geometrical and dynamic connections linking outcrop and slope-scale instabilities are often unknown. A more detailed definition of the potential links will be essential to improve the comprehension of the destabilization processes and to dispose of a more complete hazard characterization of the rock instabilities at different spatial scales. In order to propose an integrated approach in the study of the rock slope instabilities, three main themes were analysed in this PhD thesis: (1) the inventory and the spatial distribution of rock slope deformations at regional scale and their influence on the landscape evolution, (2) the influence of brittle and ductile tectonic structures on rock slope instabilities development and (3) the characterization of hazard posed by potential rock slope instabilities through the development of conceptual instability models. To prose and integrated approach for the analyses of these topics, several techniques were adopted. In particular, high resolution digital elevation models revealed to be fundamental tools that were employed during the different stages of the rock slope instability assessment. A special attention was spent in the application of digital elevation model for detailed geometrical modelling of past and potential instabilities and for the rock slope monitoring at different spatial scales. Detailed field analyses and numerical models were performed to complete and verify the remote sensing approach. In the first part of this thesis, large slope instabilities in Rhone valley (Switzerland) were mapped in order to dispose of a first overview of tectonic and climatic factors influencing their distribution and their characteristics. Our analyses demonstrate the key influence of neotectonic activity and the glacial conditioning on the spatial distribution of the rock slope deformations. Besides, the volumes of rock instabilities identified along the main Rhone valley, were then used to propose the first estimate of the postglacial denudation and filling of the Rhone valley associated to large gravitational movements. In the second part of the thesis, detailed structural analyses of the Frank slide and the Sierre rock avalanche were performed to characterize the influence of brittle and ductile tectonic structures on the geometry and on the failure mechanism of large instabilities. Our observations indicated that the geometric characteristics and the variation of the rock mass quality associated to ductile tectonic structures, that are often ignored landslide study, represent important factors that can drastically influence the extension and the failure mechanism of rock slope instabilities. In the last part of the thesis, the failure mechanisms and the hazard associated to five potential instabilities were analysed in detail. These case studies clearly highlighted the importance to incorporate different analyses and monitoring techniques to dispose of reliable and hazard scenarios. This information associated to the development of a conceptual instability model represents the primary data for an integrated risk management of rock slope instabilities. - Les mouvements de versant tels que les chutes de blocs, les éboulements ou encore les phénomènes plus lents comme les déformations gravitaires profondes de versant représentent des manifestations courantes en régions montagneuses. Les mouvements de versant sont à la fois un des facteurs principaux contrôlant la destruction progressive des chaines orogéniques mais aussi un danger naturel concret qui peut provoquer des dommages importants. Pourtant, les phénomènes gravitaires sont rarement analysés dans leur globalité et les rapports géométriques et mécaniques qui lient les instabilités à l'échelle du versant aux instabilités locales restent encore mal définis. Une meilleure caractérisation de ces liens pourrait pourtant représenter un apport substantiel dans la compréhension des processus de déstabilisation des versants et améliorer la caractérisation des dangers gravitaires à toutes les échelles spatiales. Dans le but de proposer un approche plus globale à la problématique des mouvements gravitaires, ce travail de thèse propose trois axes de recherche principaux: (1) l'inventaire et l'analyse de la distribution spatiale des grandes instabilités rocheuses à l'échelle régionale, (2) l'analyse des structures tectoniques cassantes et ductiles en relation avec les mécanismes de rupture des grandes instabilités rocheuses et (3) la caractérisation des aléas rocheux par une approche multidisciplinaire visant à développer un modèle conceptuel de l'instabilité et une meilleure appréciation du danger . Pour analyser les différentes problématiques traitées dans cette thèse, différentes techniques ont été utilisées. En particulier, le modèle numérique de terrain s'est révélé être un outil indispensable pour la majorité des analyses effectuées, en partant de l'identification de l'instabilité jusqu'au suivi des mouvements. Les analyses de terrain et des modélisations numériques ont ensuite permis de compléter les informations issues du modèle numérique de terrain. Dans la première partie de cette thèse, les mouvements gravitaires rocheux dans la vallée du Rhône (Suisse) ont été cartographiés pour étudier leur répartition en fonction des variables géologiques et morphologiques régionales. En particulier, les analyses ont mis en évidence l'influence de l'activité néotectonique et des phases glaciaires sur la distribution des zones à forte densité d'instabilités rocheuses. Les volumes des instabilités rocheuses identifiées le long de la vallée principale ont été ensuite utilisés pour estimer le taux de dénudations postglaciaire et le remplissage de la vallée du Rhône lié aux grands mouvements gravitaires. Dans la deuxième partie, l'étude de l'agencement structural des avalanches rocheuses de Sierre (Suisse) et de Frank (Canada) a permis de mieux caractériser l'influence passive des structures tectoniques sur la géométrie des instabilités. En particulier, les structures issues d'une tectonique ductile, souvent ignorées dans l'étude des instabilités gravitaires, ont été identifiées comme des structures très importantes qui contrôlent les mécanismes de rupture des instabilités à différentes échelles. Dans la dernière partie de la thèse, cinq instabilités rocheuses différentes ont été étudiées par une approche multidisciplinaire visant à mieux caractériser l'aléa et à développer un modèle conceptuel trois dimensionnel de ces instabilités. A l'aide de ces analyses on a pu mettre en évidence la nécessité d'incorporer différentes techniques d'analyses et de surveillance pour une gestion plus objective du risque associée aux grandes instabilités rocheuses.

Are the light-harvesting I complexes from Rhodospirillum rubrum arranged around the reaction centre in a square geometry?

The basic photosynthetic unit containing the reaction centre and the light-harvesting I complex (RC-LHI) of the purple non-sulphur bacterium Rhodospirillum rubrum was purified and reconstituted into two-dimensional (2D) membrane crystals. Transmission electron microscopy using conventional techniques and cryoelectron microscopy of the purified single particles and of 2D crystals yielded a projection of the RC-LHI complex at a resolution of at least 1.6 nm. In this projection the LHI ring appears to have a square symmetry and packs in a square crystal lattice. The square geometry of the LHI ring was observed also in images of single isolated particles of the RC-LHI complex. However, although the LHI units are packed identically within the crystal lattice, a new rotational analysis developed here showed that the reaction centres take up one of four possible orientations within the ring. This fourfold disorder supports our interpretation of a square ring symmetry and suggests that a hitherto undetected component may be present within the photosynthetic unit.

The Lausanne experience with the Wichita fusion nail for arthrodesis of the knee

Background: Arthrodesis of the knee by intramedullary fixation hasbeen reported to have a higher rate of success than external fixationor compression plating. Antegrade nailing however can lead to complicationsdue to the different diameters of the medullary canals, fracturesduring insertion, poor rotational stability, breakage of the IM-nailand insufficient compression at the fusion site.Method: This retrospective study reports all knee fusions performedby the same orthopaedic surgeon with the Wichita (Stryker) fusion nail(WFN) from 2004 to 2010. The Wichita nail is a short nail with a deviceat the knee which allows for coupling of differently sized and interlockedfemoral and tibial components and at the same time for compression.Results: We report of 18 patients with a mean follow up of 28 months(range 3-71 months). Infected TKA was the most common indicationfor arthrodesis in 9 cases. The remaining reasons included asepticfailed TKA in 3 cases, 2 patients after fracture, 1 patient with neurologicalinstability after knee dislocation, 1 patient after tumoral resectionand 1 non union after failed arthrodesis with long antegrade nail.Finally 1 patient with bilateral congenital knee dislocation operated onboth sides. As expected, patients receiving the WFN had undergonea large number of previous knee surgeries with a mean of 3.8 (range0-8) procedures per patient. The complication rate was 27% (5 of 18).Two patients had persistent pain requiring revision surgery to increasestability with plating. One case of periprosthetic fracture needed openreduction and internal fixation. 2 patients with superficial hematomawere treated one with open drainage and the other with physiotherapy.Infection was erradicated in all septic cases, we found no new infectionand the fusion rate was 100%.Conclusion: The results in these often difficult cases are satisfyingand we think that this technique is a valid alternative to the otherknown techniques of knee fusion in patients with a poor bone stockand fragile soft tissues.

Preliminary assessment of rockslide and rockfall hazards using a DEM (Oppstadhornet, Norway)

The increasing availability and precision of digital elevation model (DEM) helps in the assessment of landslide prone areas where only few data are available. This approach is performed in 6 main steps which include: DEM creation; identification of geomorphologic features; determination of the main sets of discontinuities; mapping of the most likely dangerous structures; preliminary rock-fall assessment; estimation of the large instabilities volumes. The method is applied to two the cases studies in the Oppstadhornet mountain (730m alt): (1) a 10 millions m3 slow-moving rockslide and (2) a potential high-energy rock falling prone area. The orientations of the foliation and of the major discontinuities have been determined directly from the DEM. These results are in very good agreement with field measurements. Spatial arrangements of discontinuities and foliation with the topography revealed hazardous structures. Maps of potential occurrence of these hazardous structures show highly probable sliding areas at the foot of the main landslide and potential rock falls in the eastern part of the mountain.

Bi-planar 2D-to-3D registration in Fourier domain for stereoscopic X-ray motion tracking

In this paper we present a new method to track bonemovements in stereoscopic X-ray image series of the kneejoint. The method is based on two different X-ray imagesets: a rotational series of acquisitions of the stillsubject knee that will allow the tomographicreconstruction of the three-dimensional volume (model),and a stereoscopic image series of orthogonal projectionsas the subject performs movements. Tracking the movementsof bones throughout the stereoscopic image series meansto determine, for each frame, the best pose of everymoving element (bone) previously identified in the 3Dreconstructed model. The quality of a pose is reflectedin the similarity between its simulated projections andthe actual radiographs. We use direct Fourierreconstruction to approximate the three-dimensionalvolume of the knee joint. Then, to avoid the expensivecomputation of digitally rendered radiographs (DRR) forpose recovery, we reformulate the tracking problem in theFourier domain. Under the hypothesis of parallel X-raybeams, we use the central-slice-projection theorem toreplace the heavy 2D-to-3D registration of projections inthe signal domain by efficient slice-to-volumeregistration in the Fourier domain. Focusing onrotational movements, the translation-relevant phaseinformation can be discarded and we only consider scalarFourier amplitudes. The core of our motion trackingalgorithm can be implemented as a classical frame-wiseslice-to-volume registration task. Preliminary results onboth synthetic and real images confirm the validity ofour approach.

Vestibulo-ocular reflex characteristics in patients with unilateral Ménière's disease.

OBJECTIVE: The pathophysiologic concept of Ménière's disease assumes that endolymphatic hydrops is the cause of the symptoms via increased pressure in the endolymphatic space and/or ionic disorder due to ruptured membrane. The goal of this study was to assess whether the vestibulo-ocular reflex (VOR) properties in patients with Ménière's disease were consistent with the classical theory. MATERIAL AND METHODS: We studied 34 patients (19 women, 15 men) presenting with unilateral Ménière's disease divided into 2 groups according to the duration of the disease: 18 were in the early stage (< 12 mo), and 16 were in the late stage (> or = 12 mo). Nineteen patients were examined during an attack. Eight of them and the 15 other patients were tested during the interval between attacks. Their characteristics were compared with those of a group of 22 normal subjects. The VOR function was evaluated via standard caloric and impulse rotatory tests (velocity step). A mathematic model of vestibular function was used to characterize the VOR response to rotational stimulation. Dynamic VOR parameters (sensitivity coefficient, time constant, gain, and asymmetries between the 2 directions of rotation) were statistically compared between the 2 groups of patients during and between attacks and between the patients and controls. RESULTS: All dynamic VOR parameters showed no statistically significant difference both with normal controls and among the patients during and between attacks (p > 0.05) except for gain asymmetry (p < or = 0.008). During attacks, patients with early Ménière's disease displayed a higher gain in rotation toward the affected ear, the opposite being observed in patients with late disease. Caloric test revealed a moderate canal paresis on the affected side during the crisis and a slight asymmetry between attacks. CONCLUSION: During attacks in patients with early Ménière's disease, the VOR gain toward the affected side is higher than that toward the intact side, supporting the fact that that the sensitivity of the cupuloendolymphatic system on the affected ear is increased. This observation is in agreement with a mechanical pressure effect of hydrops on the vestibular organs.

Geochemical and H-O-Sr-Nd isotope evidence for magmatic processes and meteoric-water interactions in the basal complex of La Gomera, Canary Islands

The plutonic rocks of the Basal Complex of La Gomera, Canary Islands, Spain, were studied by means of major and trace element contents and by H-O-Sr-Nd isotope compositions in order to distinguish primary magmatic characteristics and late-stage alteration products. Deciphering the effects of alteration allowed us to determine primary, plume-related compositions that indicated D- and (18)O-depletion relative to normal upper mantle, supporting the conclusions of earlier studies on the plutonic rocks of Fuerteventura and La Palma. Late-stage alteration took place during the formation of the intrusive series induced by interaction with meteoric water. Inferred isotopic compositions of the meteoric water indicate that the water infiltrated into the rock edifice at a height of about 1500 m above sea level, suggesting the existence of a subaerial volcano which was active during the intrusive activity and that it has been either distroyed or remain buried by later volcanic and landslide events.

Analysis of joint laxity after total ankle arthroplasty: cadaver study.

BACKGROUND: Clinical results of total ankle arthroplasty with early designs were disappointing. Recently-developed ankle prostheses have good mid-term results; however, limited information is available regarding effects of total ankle arthroplasty on ankle laxity. METHODS: Eight cadaveric lower extremities were tested with a custom device which enabled measurement of multi-axial forces, moments, and displacement during applied axial, shear, and rotational loading. Tests consisted of anterior-posterior and medial-lateral translation and internal-external rotation of the talus relative to the tibia during axial loads on the tibia simulating body weight (700 N) and an unloaded condition (5 N). Tests were performed in neutral, dorsiflexion, and plantarflexion. Laxity was determined for the intact ankle, and following insertion of an unconstrained total ankle implant, comparing load-displacement curve. FINDINGS: Laxity after total ankle arthroplasty did not approximate the normal ankle in most conditions tested. Displacement was significantly greater for total ankle arthroplasty in both posterior and lateral translation, and internal rotation, with 5 N axial loading, and anterior-posterior, medial-lateral translation, and internal-external rotation for 700 N axial loading. For the 700 N axial load condition, in the neutral ankle position, total anterior-posterior translation averaged 0.4 mm (SD 0.2 mm), but 6.0 mm (SD 1.5 mm) after total ankle arthroplasty (P<0.01). This study demonstrated more laxity in the replaced ankle than normal ankle for both unloaded and 700 N axially loaded conditions. INTERPRETATION: These data indicate the increased responsibility of the ligaments for ankle laxity after total ankle arthroplasty and suggest the importance of meticulous ligament reconstruction with total ankle arthroplasty operations.