Tectonic subsidence of the Lomonosov Ridge

The Cenozoic sedimentary record revealed by the Integrated Ocean Drilling Program's Arctic Coring Expedition (ACEX) to the Lomonosov Ridge microcontinent in 2004 is characterized by an unconformity attributed to the period 44-18 Ma. According to conventional thermal kinematic models, the microcontinent should have subsided to >1 km depth owing to rifting and subsequent separation from the Barents-Kara Sea margin at 56 Ma. We propose an alternative model incorporating a simple pressure-temperature (P-T) relation for mantle density. Using this model, we can explain the missing stratigraphic section by post-breakup uplift and erosion. The pattern of linear magnetic anomalies and the spreading geometry imply that the generation of oceanic crust in the central Eurasia Basin could have been restricted and confined by non-volcanic thinning of the mantle lithosphere at an early stage (ca. 56-40 Ma). In response to a rise in temperature, the mantle mineral composition may have changed through breakdown of spinet peridotite and formation of less dense plagioclase peridotite. The consequence of lithosphere heating and related mineral phase transitions would be post-breakup uplift followed by rapid subsidence to the deep-water environment observed on the Lomonosov Ridge today.

Effect of long-term climbing training on cerebellar ataxia: a case series.

Background. Efficient therapy for both limb and gait ataxia is required. Climbing, a complex task for the whole motor system involving balance, body stabilization, and the simultaneous coordination of all 4 limbs, may have therapeutic potential. Objective. To investigate whether long-term climbing training improves motor function in patients with cerebellar ataxia. Methods. Four patients suffering from limb and gait ataxia underwent a 6-week climbing training. Its effect on ataxia was evaluated with validated clinical balance and manual dexterity tests and with a kinematic analysis of multijoint arm and leg pointing movements. Results. The patients increased their movement velocity and achieved a more symmetric movement speed profile in both arm and leg pointing movements. Furthermore, the 2 patients who suffered the most from gait ataxia improved their balance and 2 of the 4 patients improved manual dexterity. Conclusion. Climbing training has the potential to serve as a new rehabilitation method for patients with upper and lower limb ataxia.

Plate tectonics of the Altaids

Abstract: The Altaids consist in a huge accretionary-type belt extending from Siberia through Mon-golia, northern China, Kyrgyzstan and Kazakhstan. They were formed from the Vendian through the Jurassic by the accretion of numerous displaced and exotic terranes (e.g. island arc, ribbon microcontinent, seamount, basaltic plateau, back-arc basin). The number, nature and origin of the terranes differ according to the palaeotectonic models of the different authors. Thanks to a geo- dynamic study (i.e. definition of tectonic settings and elaboration of geodynamic scenarios) and plate tectonics modelling, this work aims to present an alternative model explaining the Palaeozoic palaeotectonic evolution of the Altaids. Based on a large set of compiled geological data related to palaeogeography and geodyna¬mic (e.g. sedimentology, stratigraphy, palaeobiogeography, palaeomagnetism, magmatism, me- tamorphism, tectonic...), a partly new classification of the terranes and sutures implicated in the formation of the Altaids is proposed. In the aim to elaborate plate tectonics reconstructions, it is necessary to fragment the present arrangement of continents into consistent geological units. To avoid confusion with existing terminology (e.g. tectonic units, tectono-stratigraphic units, micro- continents, terranes, blocks...), the new concept of "Geodynamic Units (GDU)" was introduced. A terrane may be formed by a set of GDUs. It consists of a continental and/or oceanic fragment which has its own kinematic and geodynamic evolution for a given period. With the same ap-proach, the life span and type of the disappeared oceans is inferred thanks to the study of the mate-rial contained in suture zones. The interpretation of the tectonic settings within the GDUs comple-ted by the restoration of oceans leads to the elaboration of geodynamic scenarios. Since the Wilson cycle was presented in 1967, numerous works demonstrated that the continental growth is more complex and results from diverse geodynamic scenarios. The identification of these scenarios and their exploitation enable to elaborate plate tectonics models. The models are self-constraining (i.e. space and time constraints) and contest or confirm in turn the geodynamic scenarios which were initially proposed. The Altaids can be divided into three domains: (1) the Peri-Siberian, (2) the Kazakhstan, and (3) the Tarim-North China domains. The Peri-Siberian Domain consists of displaced (i.e. Sayan Terrane Tuva-Mongolian, Lake-Khamsara Terrane) and exotic terranes (i.e. Altai-Mongolian and Khangai-Argunsky Terrane) accreted to Siberia from the Vendian through the Ordovician. Fol-lowing the accretion of these terranes, the newly formed Siberia active margin remained active un-til its part collision with the Kazakhstan Superterrane in the Carboniferous. The eastern part of the active margin (i.e. East Mongolia) continued to act until the Permian when the North-China Tarim Superterrane collided with it. The geodynamic evolution of the eastern part of the Peri-Siberian Domain (i.e. Eastern Mongolia and Siberia) is complicated by the opening of the Mongol-Okhotsk Ocean in the Silurian. The Kazakhstan Domain is composed of several continental terranes of East Gondwana origin amalgamated together during the Ordovician-Silurian time. After these different orogenic events, the Kazakhstan Superterrane evolved as a single superterrane until its collision with a Tarim-North China related-terrane (i.e. Tianshan-Hanshan Terrane) and Siberian Continent during the Devonian. This new organisation of the continents imply a continued active margin from Siberia, to North China through the Kazakhstan Superterrane and the closure of the Junggar- Balkash Ocean which implied the oroclinal bending of the Kazakhstan Superterrane during the entire Carboniferous. The formation history of the Tarim-North China Domain is less complex. The Cambrian northern passive margin became active in the Ordovician. In the Silurian, the South Tianshan back-arc Ocean was open and led to the formation of the Tianshan-Hanshan Terrane which collided with the Kazakhstan Superterrane during the Devonian. The collision between Siberia and the eastern part of the Tarim-North China continents (i.e. Inner Mongolia), implied by the closure of the Solonker Ocean, took place in the Permian. Since this time, the major part of the Altaids was formed, the Mongol-Okhotsk Ocean only was still open and closed during the Jurassic. Résumé: La chaîne des Altaïdes est une importante chaîne d'accrétion qui s'étend en Sibérie, Mon-golie, Chine du Nord, Kirghizstan et Kazakhstan. Elle s'est formée durant la période du Vendian au Jurassique par l'accrétion de nombreux terranes déplacés ou exotiques (par exemple arc océa-nique, microcontinent, guyot, plateau basaltique, basin d'arrière-arc...). Le nombre, la nature ou encore l'origine diffèrent selon les modèles paléo-tectoniques proposés par les différents auteurs. Grâce à une étude géodynamique (c'est-à-dire définition des environnements tectoniques et éla-boration de scénarios géodynamiques) et à la modélisation de la tectonique des plaques, ce travail propose un modèle alternatif expliquant l'évolution paléo-tectonique des Altaïdes. Basé sur une large compilation de données géologiques pertinentes en termes de paléo-géographie et de géodynamique (par exemple sédimentologie, stratigraphie, paléo-biogéographie, paléomagnétisme, magmatisme, métamorphisme, tectonique...), une nouvelle classification des terranes et des sutures impliqués dans la formation des Altaïdes est proposée. Dans le but d'élabo¬rer des reconstructions de plaques tectoniques, il est nécessaire de fragmenter l'arrangement actuel des continents en unités tectoniques cohérentes. Afin d'éviter les confusions avec la terminolo¬gie existante (par exemple unité tectonique, unité tectono-stratigraphique, microcontinent, block, terrane...), le nouveau concept d' "Unité Géodynamique (UGD)" a été introduit. Un terrane est formé d'une ou plusieurs UGD et représente un fragment océanique ou continental défini pas sa propre cinétique et évolution géodynamique pour une période donnée. Parallèlement, la durée de vie et le type des océans disparus (c'est-à-dire principal ou secondaire) est déduite grâce à l'étude du matériel contenu dans les zones de sutures. L'interprétation des environnements tectoniques des UGD associés à la restauration des océans mène à l'élaboration de scénarios géodynamiques. Depuis que le Cycle de Wilson a été présenté en 1967, de nombreux travaux ont démontré que la croissance continentale peut résulter de divers scénarios géodynamiques. L'identification et l'ex-ploitation de ces scénarios permet finalement l'élaboration de modèles de tectonique des plaques. Les modèles sont auto-contraignants (c'est-à-dire contraintes spatiales et temporelles) et peuvent soit contester ou confirmer les scénarios géodynamiques initialement proposés. Les Altaïdes peuvent être divisées en trois domaines : (1) le Domaine Péri-Sibérien, (2) le Domaine Kazakh, et (3) le Domaine Tarim-Nord Chinois. Le Domaine Péri-Sibérien est composé de terranes déplacés (c'est-à-dire Terrane du Sayan, Tuva-Mongol et Lake-Khamsara) et exotiques (c'est-à-dire Terrane Altai-Mongol et Khangai-Argunsky) qui ont été accrétés au craton Sibérien durant la période du Vendien à l'Ordovicien. Suite à l'accrétion de ces terranes, la marge sud-est de la Sibérie nouvellement formée reste active jusqu'à sa collision partielle avec le Superterrane Ka-zakh au Carbonifère. La partie est de la marge active (c'est-à-dire Mongolie de l'est) continue son activité jusqu'au Permien lors de sa collision avec le Superterrane Tarim-Nord Chinois. L'évolu¬tion géodynamique de la partie est du Domaine Sibérien est compliquée par l'ouverture Silurienne de l'Océan Mongol-Okhotsk qui disparaîtra seulement au Jurassique. Le Domaine Kazakh est composé de plusieurs terranes d'origine est-Gondwanienne accrétés les uns avec les autres avant ou pendant le Silurien inférieur et leurs evolution successive sous la forme d'un seul superterrane. Le Superterrane Kazakh collisione avec un terrane Tarim-Nord Chinois (c'est-à-dire Terrane du Tianshan-Hanshan) durant le Dévonien et le continent Sibérien au Dévonien supérieur. Ce nouvel agencement des plaques induit une marge active continue le long des continents Sibérien, Kazakh et Nord Chinois et la fermeture de l'Océan Junggar-Balkash qui provoque le plissement oroclinal du Superterrane Kazakh durant le Carbonifère. L'histoire de la formation du Domaine Tarim-Nord Chinois est moins complexe. La marge passive nord Cambrienne devient active à l'Ordovicien et l'ouverture Silurienne du bassin d'arrière-arc du Tianshan sud mène à la formation du terrane du Tianshan-Hanshan. La collision Dévonienne entre ce dernier et le Superterrane Kazakh provoque la fermerture de l'Océan Tianshan sud. Finalement, la collision entre la Sibérie et la partie est du continent Tarim-Nord Chinois (c'est-à-dire Mongolie Intérieure) prend place durant le Permien suite à la fermeture de l'Océan Solonker. La majeure partie des Altaïdes est alors formée, seul l'Océan Mongol-Okhotsk est encore ouvert. Ce dernier se fermera seulement au Jurassique.

New microfluidic-based sampling procedure for overcoming the hematocrit problem associated with dried blood spot analysis.

Hematocrit (Hct) is one of the most critical issues associated with the bioanalytical methods used for dried blood spot (DBS) sample analysis. Because Hct determines the viscosity of blood, it may affect the spreading of blood onto the filter paper. Hence, accurate quantitative data can only be obtained if the size of the paper filter extracted contains a fixed blood volume. We describe for the first time a microfluidic-based sampling procedure to enable accurate blood volume collection on commercially available DBS cards. The system allows the collection of a controlled volume of blood (e.g., 5 or 10 μL) within several seconds. Reproducibility of the sampling volume was examined in vivo on capillary blood by quantifying caffeine and paraxanthine on 5 different extracted DBS spots at two different time points and in vitro with a test compound, Mavoglurant, on 10 different spots at two Hct levels. Entire spots were extracted. In addition, the accuracy and precision (n = 3) data for the Mavoglurant quantitation in blood with Hct levels between 26% and 62% were evaluated. The interspot precision data were below 9.0%, which was equivalent to that of a manually spotted volume with a pipet. No Hct effect was observed in the quantitative results obtained for Hct levels from 26% to 62%. These data indicate that our microfluidic-based sampling procedure is accurate and precise and that the analysis of Mavoglurant is not affected by the Hct values. This provides a simple procedure for DBS sampling with a fixed volume of capillary blood, which could eliminate the recurrent Hct issue linked to DBS sample analysis.

Preliminary description and slope stability analyses of the 2008 Little Salmon lake and 2007 Mt. Steele landslides, Yukon

In August 2008, reactivation of the Little Salmon Lake landslide occurred. During this event, hundreds of conical mounds of variable size and composition formed in the deposition zone. The characteristics of these landforms are described and a potential mechanism for their formation is proposed. A preliminary slope stability analysis of the 2007 Mount Steele rock and ice avalanche was also undertaken. The orientation of very high persistence (>20 m long) structural planes (e.g., faults, joints and bedding) within bedrock in the source zone was obtained using an airborne-LiDAR digital elevation model and the software COLTOP-3D. Using these discontinuity orientation measurements, kinematic, surface wedge and simple three-dimensional distinct element slope stability analyses were performed.

A two-phase composite in simple shear: Effective mechanical anisotropy development and localization potential

We present a combined shape and mechanical anisotropy evolution model for a two-phase inclusion-bearing rock subject to large deformation. A single elliptical inclusion embedded in a homogeneous but anisotropic matrix is used to represent a simplified shape evolution enforced on all inclusions. The mechanical anisotropy develops due to the alignment of elongated inclusions. The effective anisotropy is quantified using the differential effective medium (DEM) approach. The model can be run for any deformation path and an arbitrary viscosity ratio between the inclusion and host phase. We focus on the case of simple shear and weak inclusions. The shape evolution of the representative inclusion is largely insensitive to the anisotropy development and to parameter variations in the studied range. An initial hardening stage is observed up to a shear strain of gamma = 1 irrespective of the inclusion fraction. The hardening is followed by a softening stage related to the developing anisotropy and its progressive rotation toward the shear direction. The traction needed to maintain a constant shear rate exhibits a fivefold drop at gamma = 5 in the limiting case of an inviscid inclusion. Numerical simulations show that our analytical model provides a good approximation to the actual evolution of a two-phase inclusion-host composite. However, the inclusions develop complex sigmoidal shapes resulting in the formation of an S-C fabric. We attribute the observed drop in the effective normal viscosity to this structural development. We study the localization potential in a rock column bearing varying fraction of inclusions. In the inviscid inclusion case, a strain jump from gamma = 3 to gamma = 100 is observed for a change of the inclusion fraction from 20% to 33%.

Fate of Hydrophilic Nanoparticles in Biological Environments

The nanoparticles developed are based on chitosan, a biocompatible and biodegradable polysaccharide. The chitosan nanoparticles are formed in an entirely water-based process by electrostatic interactions with other biocompatible molecules. As a prerequisite to understand the fate of such nanoparticles in cells, comprehensive characterization and stability studies serve to identify quantitatively the impact of the raw material characteristics and preparation conditions on the nanoparticle characteristics. Methods included H-1 NMR spectroscopy, dilution viscometry, particle size analysis and electron microscopy. Cytotoxicity and cell uptake experiments on RAW 264.7 murine macrophages and p23 murine endothelial cells were performed to investigate the correlation with nanoparticle characteristics and effect of surface decoration with alginate. Cytotoxicity was assessed by the MTT survival test; cell uptake was monitored by fluorescent microscopy using labeled polymers.

Transcranial direct current stimulation for the treatment of focal hand dystonia.

The treatment of writer's cramp, a task-specific focal hand dystonia, needs new approaches. A deficiency of inhibition in the motor cortex might cause writer's cramp. Transcranial direct current stimulation modulates cortical excitability and may provide a therapeutic alternative. In this randomized, double-blind, sham-controlled study, we investigated the efficacy of cathodal stimulation of the contralateral motor cortex in 3 sessions in 1 week. Assessment over a 2-week period included clinical scales, subjective ratings, kinematic handwriting analysis, and neurophysiological evaluation. Twelve patients with unilateral dystonic writer's cramp were investigated; 6 received transcranial direct current and 6 sham stimulation. Cathodal transcranial direct current stimulation had no favorable effects on clinical scales and failed to restore normal handwriting kinematics and cortical inhibition. Subjective worsening remained unexplained, leading to premature study termination. Repeated sessions of cathodal transcranial direct current stimulation of the motor cortex yielded no favorable results supporting a therapeutic potential in writer's cramp.

A role for septins in the interaction between the Listeria monocytogenes INVASION PROTEIN InlB and the Met receptor.

Septins are conserved GTPases that form filaments and are required for cell division. During interphase, septin filaments associate with cellular membrane and cytoskeleton networks, yet the functional significance of these associations have, to our knowledge, remained unknown. We recently discovered that different septins, SEPT2 and SEPT11, regulate the InlB-mediated entry of Listeria monocytogenes into host cells. Here we address the role of SEPT2 and SEPT11 in the InlB-Met interactions underlying Listeria invasion to explore how septins modulate surface receptor function. We observed that differences in InlB-mediated Listeria entry correlated with differences in Met surface expression caused by septin depletion. Using atomic force microscopy on living cells, we show that septin depletion significantly reduced the unbinding force of InlB-Met interaction and the viscosity of membrane tethers at locations where the InlB-Met interaction occurs. Strikingly, the same order of difference was observed for cells in which the actin cytoskeleton was disrupted. Consistent with a proposed role of septins in association with the actin cytoskeleton, we show that cell elasticity is decreased upon septin or actin inactivation. Septins are therefore likely to participate in anchorage of the Met receptor to the actin cytoskeleton, and represent a critical determinant in surface receptor function.

Direct numerical simulations of interface dynamics to link capillary pressure and total surface energy

The flow of two immiscible fluids through a porous medium depends on the complex interplay between gravity, capillarity, and viscous forces. The interaction between these forces and the geometry of the medium gives rise to a variety of complex flow regimes that are difficult to describe using continuum models. Although a number of pore-scale models have been employed, a careful investigation of the macroscopic effects of pore-scale processes requires methods based on conservation principles in order to reduce the number of modeling assumptions. In this work we perform direct numerical simulations of drainage by solving Navier-Stokes equations in the pore space and employing the Volume Of Fluid (VOF) method to track the evolution of the fluid-fluid interface. After demonstrating that the method is able to deal with large viscosity contrasts and model the transition from stable flow to viscous fingering, we focus on the macroscopic capillary pressure and we compare different definitions of this quantity under quasi-static and dynamic conditions. We show that the difference between the intrinsic phase-average pressures, which is commonly used as definition of Darcy-scale capillary pressure, is subject to several limitations and it is not accurate in presence of viscous effects or trapping. In contrast, a definition based on the variation of the total surface energy provides an accurate estimate of the macroscopic capillary pressure. This definition, which links the capillary pressure to its physical origin, allows a better separation of viscous effects and does not depend on the presence of trapped fluid clusters.

Therapeutic applications of viscous and injectable poly(ortho esters).

Poly(ortho esters) (POE) are hydrophobic and bioerodible polymers that have been investigated for pharmaceutical use since the early 1970s. Among the four described generations of POE, the third (POE III) and fourth (POE IV) are promising viscous and injectable materials which have been investigated in numerous biomedical applications. POE III has been extensively studied for ophthalmic drug delivery, it presents an excellent biocompatibility and is currently being investigated as a vehicle for sustained drug delivery to treat diseases of the posterior segment of the eye. POE IV is distinguishable by a highly reproducible and controlled synthesis, a higher hydrophobicity, and an excellent biocompatibility. It is currently under development for a variety of applications, such as ocular delivery, periodontal disease treatment and applications in veterinary medicine. This review will also focus on new perspectives for this promising family of polymers, such as guided tissue regeneration, treatment of osteoarthritis, as well as peptide and protein delivery.

Pore-scale modeling of two-phase flow instabilities in porous media

Les problèmes d'écoulements multiphasiques en média poreux sont d'un grand intérêt pour de nombreuses applications scientifiques et techniques ; comme la séquestration de C02, l'extraction de pétrole et la dépollution des aquifères. La complexité intrinsèque des systèmes multiphasiques et l'hétérogénéité des formations géologiques sur des échelles multiples représentent un challenge majeur pour comprendre et modéliser les déplacements immiscibles dans les milieux poreux. Les descriptions à l'échelle supérieure basées sur la généralisation de l'équation de Darcy sont largement utilisées, mais ces méthodes sont sujettes à limitations pour les écoulements présentant de l'hystérèse. Les avancées récentes en terme de performances computationnelles et le développement de méthodes précises pour caractériser l'espace interstitiel ainsi que la distribution des phases ont favorisé l'utilisation de modèles qui permettent une résolution fine à l'échelle du pore. Ces modèles offrent un aperçu des caractéristiques de l'écoulement qui ne peuvent pas être facilement observées en laboratoire et peuvent être utilisé pour expliquer la différence entre les processus physiques et les modèles à l'échelle macroscopique existants. L'objet premier de la thèse se porte sur la simulation numérique directe : les équations de Navier-Stokes sont résolues dans l'espace interstitiel et la méthode du volume de fluide (VOF) est employée pour suivre l'évolution de l'interface. Dans VOF, la distribution des phases est décrite par une fonction fluide pour l'ensemble du domaine et des conditions aux bords particulières permettent la prise en compte des propriétés de mouillage du milieu poreux. Dans la première partie de la thèse, nous simulons le drainage dans une cellule Hele-Shaw 2D avec des obstacles cylindriques. Nous montrons que l'approche proposée est applicable même pour des ratios de densité et de viscosité très importants et permet de modéliser la transition entre déplacement stable et digitation visqueuse. Nous intéressons ensuite à l'interprétation de la pression capillaire à l'échelle macroscopique. Nous montrons que les techniques basées sur la moyenne spatiale de la pression présentent plusieurs limitations et sont imprécises en présence d'effets visqueux et de piégeage. Au contraire, une définition basée sur l'énergie permet de séparer les contributions capillaires des effets visqueux. La seconde partie de la thèse est consacrée à l'investigation des effets d'inertie associés aux reconfigurations irréversibles du ménisque causé par l'interface des instabilités. Comme prototype pour ces phénomènes, nous étudions d'abord la dynamique d'un ménisque dans un pore angulaire. Nous montrons que, dans un réseau de pores cubiques, les sauts et reconfigurations sont si fréquents que les effets d'inertie mènent à différentes configurations des fluides. A cause de la non-linéarité du problème, la distribution des fluides influence le travail des forces de pression, qui, à son tour, provoque une chute de pression dans la loi de Darcy. Cela suggère que ces phénomènes devraient être pris en compte lorsque que l'on décrit l'écoulement multiphasique en média poreux à l'échelle macroscopique. La dernière partie de la thèse s'attache à démontrer la validité de notre approche par une comparaison avec des expériences en laboratoire : un drainage instable dans un milieu poreux quasi 2D (une cellule Hele-Shaw avec des obstacles cylindriques). Plusieurs simulations sont tournées sous différentes conditions aux bords et en utilisant différents modèles (modèle intégré 2D et modèle 3D) afin de comparer certaines quantités macroscopiques avec les observations au laboratoire correspondantes. Malgré le challenge de modéliser des déplacements instables, où, par définition, de petites perturbations peuvent grandir sans fin, notre approche numérique apporte de résultats satisfaisants pour tous les cas étudiés. - Problems involving multiphase flow in porous media are of great interest in many scientific and engineering applications including Carbon Capture and Storage, oil recovery and groundwater remediation. The intrinsic complexity of multiphase systems and the multi scale heterogeneity of geological formations represent the major challenges to understand and model immiscible displacement in porous media. Upscaled descriptions based on generalization of Darcy's law are widely used, but they are subject to several limitations for flow that exhibit hysteric and history- dependent behaviors. Recent advances in high performance computing and the development of accurate methods to characterize pore space and phase distribution have fostered the use of models that allow sub-pore resolution. These models provide an insight on flow characteristics that cannot be easily achieved by laboratory experiments and can be used to explain the gap between physical processes and existing macro-scale models. We focus on direct numerical simulations: we solve the Navier-Stokes equations for mass and momentum conservation in the pore space and employ the Volume Of Fluid (VOF) method to track the evolution of the interface. In the VOF the distribution of the phases is described by a fluid function (whole-domain formulation) and special boundary conditions account for the wetting properties of the porous medium. In the first part of this thesis we simulate drainage in a 2-D Hele-Shaw cell filled with cylindrical obstacles. We show that the proposed approach can handle very large density and viscosity ratios and it is able to model the transition from stable displacement to viscous fingering. We then focus on the interpretation of the macroscopic capillary pressure showing that pressure average techniques are subject to several limitations and they are not accurate in presence of viscous effects and trapping. On the contrary an energy-based definition allows separating viscous and capillary contributions. In the second part of the thesis we investigate inertia effects associated with abrupt and irreversible reconfigurations of the menisci caused by interface instabilities. As a prototype of these phenomena we first consider the dynamics of a meniscus in an angular pore. We show that in a network of cubic pores, jumps and reconfigurations are so frequent that inertia effects lead to different fluid configurations. Due to the non-linearity of the problem, the distribution of the fluids influences the work done by pressure forces, which is in turn related to the pressure drop in Darcy's law. This suggests that these phenomena should be taken into account when upscaling multiphase flow in porous media. The last part of the thesis is devoted to proving the accuracy of the numerical approach by validation with experiments of unstable primary drainage in a quasi-2D porous medium (i.e., Hele-Shaw cell filled with cylindrical obstacles). We perform simulations under different boundary conditions and using different models (2-D integrated and full 3-D) and we compare several macroscopic quantities with the corresponding experiment. Despite the intrinsic challenges of modeling unstable displacement, where by definition small perturbations can grow without bounds, the numerical method gives satisfactory results for all the cases studied.

Simulation of surface waves in porous media

We present a novel numerical algorithm for the simulation of seismic wave propagation in porous media, which is particularly suitable for the accurate modelling of surface wave-type phenomena. The differential equations of motion are based on Biot's theory of poro-elasticity and solved with a pseudospectral approach using Fourier and Chebyshev methods to compute the spatial derivatives along the horizontal and vertical directions, respectively. The time solver is a splitting algorithm that accounts for the stiffness of the differential equations. Due to the Chebyshev operator the grid spacing in the vertical direction is non-uniform and characterized by a denser spatial sampling in the vicinity of interfaces, which allows for a numerically stable and accurate evaluation of higher order surface wave modes. We stretch the grid in the vertical direction to increase the minimum grid spacing and reduce the computational cost. The free-surface boundary conditions are implemented with a characteristics approach, where the characteristic variables are evaluated at zero viscosity. The same procedure is used to model seismic wave propagation at the interface between a fluid and porous medium. In this case, each medium is represented by a different grid and the two grids are combined through a domain-decomposition method. This wavefield decomposition method accounts for the discontinuity of variables and is crucial for an accurate interface treatment. We simulate seismic wave propagation with open-pore and sealed-pore boundary conditions and verify the validity and accuracy of the algorithm by comparing the numerical simulations to analytical solutions based on zero viscosity obtained with the Cagniard-de Hoop method. Finally, we illustrate the suitability of our algorithm for more complex models of porous media involving viscous pore fluids and strongly heterogeneous distributions of the elastic and hydraulic material properties.

Structural analysis, clay mineralogy and K-Ar dating of fault gouges from Centovalli Line (Central Alps) for reconstruction of their recent activity

Between the cities of Domodossola and Locarno, the complex ``Centovalli Line'' tectonic zone of the Central Alps outlines deformation phases over a long period of time (probably starting similar to 30 Ma ago) and under variable P-T conditions. The last deformation phases developed gouge-bearing faults with a general E-W trend that crosscuts the roots of the Alpine Canavese zone and the Finero ultramafic body. Kinematic indicators show that the general motion was mainly dextral associated with back thrusting towards the S. The <2 mu m clay fractions of fault gouges from Centovalli Line consist mainly of illite, smectite and chlorite with varied illite-smectite, chlorite-smectite and chlorite-serpentine mixed-layers. Constrained with the illite crystallinity index, the thermal conditions induced by the tectonic activity show a gradual trend from anchizonal to diagenetic conditions. The <2 and <0.2 mu M clay fractions, and hydrothermal K-feldspar separates all provide K-Ar ages between 14.2 +/- 2.9 Ma and roughly 0 Ma, with major episodes at about 12,8, 6 and close to 0 Ma These ages set the recurrent tectonic activity and the associated fluid circulations between Upper Miocene and Recent. On the basis of the K-Ar ages and with a thermal gradient of 25-30 degrees C/km, the studied fault zones were located at a depth of 4-7 km. If they were active until now as observed in field, the exhumation was approximately 2.5-3.0 km for the last 12 Ma with a mean velocity of 0.4 mm/y. Comparison with available models on the recent Alpine evolution shows that the tectonic activity in the area relates to a continuum of the back-thrusting movements of the Canavese Line, and/or to several late-extensional phases of the Rhone-Simplon line. The Centovalli-Val Vigezzo zone therefore represents a major tectonic zone of the Central-Western Alps resulting from different interacting tectonic events. (C) 2011 Elsevier B.V. All rights reserved.

Combining digital elevation model analysis and run-out modeling to characterize hazard posed by a potentially unstable rock slope at Turtle Mountain, Alberta, Canada

Turtle Mountain in Alberta, Canada has become an important field laboratory for testing different techniques related to the characterization and monitoring of large slope mass movements as the stability of large portions of the eastern face of the mountain is still questionable. In order to better quantify the volumes potentially unstable and the most probable failure mechanisms and potential consequences, structural analysis and runout modeling were preformed. The structural features of the eastern face were investigated using a high resolution digital elevation model (HRDEM). According to displacement datasets and structural observations, potential failure mechanisms affecting different portions of the mountain have been assessed. The volumes of the different potentially unstable blocks have been calculated using the Sloping Local Base Level (SLBL) method. Based on the volume estimation, two and three dimensional dynamic runout analyses have been performed. Calibration of this analysis is based on the experience from the adjacent Frank Slide and other similar rock avalanches. The results will be used to improve the contingency plans within the hazard area.