Les écoulements pyroclastiques sont des écoulements volcaniques complexes dont le comportement physique fait encore l'objet de débats. Ils sont composés de deux parties : l'écoulement dense basal, riche en particules et en blocs, surmonté par la déferlante, diluée et turbulente. Les interactions entre ces deux parties ne sont pas bien comprises, tout comme leurs échanges de masses et de quantités de mouvement. Partant de ce constat, cette thèse se concentre sur l’étude des mécanismes de formation de la déferlante à partir de l’écoulement dense.Les expériences mettent en évidence un mécanisme de formation d'un écoulement dilué par l’alternance d’incorporation d'air et d’élutriation des particules fines d’un lit granulaire dense soumis à des vibrations. L'air est aspiré dans le lit granulaire pendant les phases de dilatation puis expulsé pendant les phases de contraction. Une partie des particules est alors soutenue par l'air turbulent expulsé et forme un mélange de gaz et de particules qui, plus dense que l’air, se transforme en un écoulement de gravité. Extrapolé à l’échelle d’un volcan, ce mécanisme d’incorporation d’air et d’élutriation peut être reproduit par une topographie rugueuse, où chaque obstacle génère une compaction puis une dilation de l’écoulement dense. La quantification du mécanisme a été effectuée et l’approche expérimentale a permis d’aboutir à une loi reliant le flux de masse de la partie dense vers la déferlante à la vitesse de l’écoulement dense. Le modèle numérique est utilisé dans un premier temps pour étudier la rhéologie de l’écoulement dense qui, en contrôlant sa vitesse, contrôle le flux de masse précédemment évoqué. Un chapitre est consacré à l’effet de la fluidisation de l’écoulement dense sur sa rhéologie. Les résultats montrent que la fluidisation par les gaz est capable d’expliquer à la fois la grande mobilité de ces écoulements, ainsi que la formation des morphologies terminales en lobes et chenaux. L’ingestion d’air dans un écoulement au cours de sa mise en place semble pouvoir expliquer une partie de la dynamique des écoulements denses. Des rhéologies simples, de premier ordre, ont également été analysées : la rhéologie de Coulomb, la rhéologie plastique, et la rhéologie à coefficient de frottement variable. Les résultats montrent que la rhéologie plastique semble la mieux adaptée pour reproduire la vitesse et l’extension des écoulements denses.Ce modèle numérique a ensuite été utilisé pour tester la loi de flux de masse obtenue suite aux expériences de laboratoire. Appliqués à l’effondrement de dôme du 25 juin 1997 à la Soufriere Hills de Montserrat, les résultats montrent que les simulations reproduisent des dépôts de déferlantes dont l’épaisseur et l’extension sont tout à fait réalistes. Les simulations reproduisent même les écoulements denses secondaires issus de la sédimentation de la déferlante puis de la remobilisation des dépôts. Les cycles d’ingestion/expulsion d’air dans l’écoulement dense, par interaction avec la topographie, expliqueraient donc à la fois la grande fluidité des écoulements denses et la formation des déferlantes pyroclastiques. Les résultats de cette thèse mettent à jour un mécanisme nouveau qui pourrait être la clé de la mise en place des écoulements pyroclastiques et pourrait permettre d’améliorer la prévision future des risques et des menaces par modélisation numérique. / Small volume pyroclastic density currents are complex volcanic flows, whose physical behaviour is still debated. They comprise two parts: the pyroclastic flow, rich in particles and blocks, overridden by the ash-cloud surge, a turbulent and dilute flow. The interactions between these two parts are not fully understood, as well as their exchanges of mass and momentum. Therefore, the thesis focuses on the investigation of ash-cloud surge formation mechanisms from the pyroclastic flow. The experiments reveal a mechanism of dilute flow formation by alternation of air incorporation into and elutriation of fine particles from a dense granular bed subjected to vibrations. The air is aspirated into the granular bed during dilatations, and expulsed during the contraction phases. A part of the particles are then sustained by the turbulent expulsed air and form a mixture of gas and particles that transforms into a gravity current. Extrapolated to a volcanic edifice, this mechanism of air incorporation and elutriation can be reproduced by a rough topography, where each obstacle generates a compaction followed by a dilatation of the pyroclastic flow. The quantification of the mechanism has been accomplished and the mass flux from the dense flow to the ash-cloud surge has been deduced.The numerical model is first used to study the pyroclastic flow rheology, which controls the velocity of the flow, and then the mass flux previously mentioned. One chapter is dedicated to the fluidization effect on the pyroclastic flow rheology. Results show that this mechanism can explain the long runout of these flows, and also the formation of levées and channel morphologies. The air ingestion in the flow during its movement could explain a part of the pyroclastic flows dynamic. Simple rheologies has also been analyzed: a Coulomb rheology, a plastic rheology, and a variable friction coefficient rheology. Results show that the plastic rheology seems to be the most adapted rheology to simulate the pyroclastic flow dynamic. Then, the numerical model has been used to test the mass flow law obtained through experiments. Applied to the 25 June 1997 dome collapse at Soufrière Hills Volcano at Montserrat, results show that the simulations reproduce accurately the extension and the thickness of the surge deposits. The simulations are also able to reproduce the surge derived pyroclastic flow, generated by remobilisation of surge deposits. The cycles of ingestion/expulsion of air in the pyroclastic flow by interactions with the topography could explain both the great fluidity of these flows and the formation of ash-cloud surge. These results highlight a new mechanism that could be a key process in pyroclastic flow dynamic, which could improve significantly the hazard and risk assessment using numerical model.

L' étude des écoulements pyroclastiques de la partie Est du strato-volcan cantalien est abordée sur l'exemple de ceux de la vallée de l'Alagnon. L'observation du terrain et des microstructures permet, compte tenu des lois générales de la mécanique des fluides, de déduire la dynamique des écoulements pyroclastiques. Du point de vue méthodologique, j'ai cherché à définir ces écoulements avant de les situer dans la stratigraphie. Trois points ont été successivement abordés: 1) les coulées pyroclastiques (s.s.) à caractère ignimbritique ; 2) les écoulements pyroclastiques (s.s.) : brèches et conglomérats; 3) la position stratigraphique des écoulements pyroclastiques (s.1.). 1) Les coulées pyroclastiques : Leur étude permet d'opposer à l'hypothèse classique des courants fluidifiés turbulents, un modèle d'écoulement de gravité, dense, visqueux et laminaire, issu d'une coulée de lave. 2) Les écoulements pyroclastiques : Une solution à la genèse et à la mise en place des brèches et des conglomérats du Cantal, se déduit de la découverte de deux appareils volcaniques surbaissées et originaux : Auzolles et Seycheuse. Leur étude permet de démontrer le passage d'une coulée de lave émulsionnée (juxtaposition de gouttelettes trachyandésitiques et basaltiques) à un écoulement pyroclastique. Ce dernier a le comportement d'un écoulement de gravité dense, visqueux et laminaire mais eu diffère par son aspect. D'autres exemples permettent de préciser le modèle selon la variation de structure du liquide parental, la morphologie et les modifications subies par l'écoulement durant sa mise en place. 3) Position stratigraphique : J'ai eouligné la relation existant entre la dynamique des écoulements pyroclastiquee et les principales étapes de la construction du Cantal

Les écoulements pyroclastiques sont des mélanges à haute température de gaz et de particules volcaniques qui peuvent se propager sur de très grandes distances. Cette forte « mobilité » est souvent attribuée à leur capacité à se fluidiser, c’est-à-dire à générer et conserver une forte pression interstitielle de gaz qui réduit les forces de friction interne. L’objectif principal de cette thèse est de comprendre comment les irrégularités des terrains sur lesquelles se propagent les écoulements pyroclastiques peuvent favoriser leur fluidisation. Une première série d’expériences de laboratoire a consisté à générer des écoulements de particules fines (diamètre de 45-90 μm) sur des substrats de différentes rugosités. Les résultats montrent que la distance de parcours des écoulements augmente avec la rugosité, allant jusqu’à doubler par rapport à la distance de parcours sur fond lisse. Des analyses de vidéos haute vitesse et des mesures de pression interstitielle d’air à la base des écoulements montrent que la tête (partie antérieure) des écoulements qui se propagent sur un substrat rugueux s’auto-fluidisent en conséquence de la sédimentation des particules dans les interstices du substrat, chassant l’air qui remonte et percole dans l’écoulement. Ce mécanisme d’auto-fluidisation est efficace pour toutes les inclinaisons étudiées (0-30°), suggérant qu’il est susceptible de se produire tout au long de la mise en place d’un écoulement pyroclastique. Une seconde étude a consisté à faire chuter des lits de particules dans une colonne statique. Les résultats montrent que même pour une hauteur de relâchement relativement faible (20 cm), le mélange peut entièrement s’auto-fluidiser durant sa chute. Quand les particules sont suffisamment fines (<100 μm) la pression interstitielle dans le dépôt diffuse pendant plusieurs secondes, la durée de cette diffusion augmentant avec l’augmentation de l’épaisseur du lit et la diminution de taille des particules. Les temps de diffusions les plus longs sont observés avec un matériau provenant d’un dépôt d’écoulement pyroclastique (~30 s pour des lits de 28.5 cm d’épaisseur). Ces résultats suggèrent que les écoulements pyroclastiques qui se propagent sur des terrains accidentés peuvent s’auto-fluidiser et conserver une faible friction au cours de leur mise en place. / Pyroclastic flows are hot mixtures of gas and particles that can propagate over large distances. This high “mobility” is often attributed to their ability to be fluidized, that is, to generate and retain high gas pore pressure that reduces internal friction forces. The main objective of this thesis is to understand how irregularities of substrates on which pyroclastic flows propagate can enhance their fluidization. A first set of laboratory experiments consisted of the generation of fine-grained flows (diameter of 45-90 μm) on substrate of various roughness. Results show that the flow runout distance increases with the substrate roughness, and is up to twice the runout on a smooth substrate. High speed video analyses and air pore pressure measurements at the flow base show that the flow head propagating over a rough substrate can auto-fluidize because of particles sedimentation into the substrate interstices, which forces the air to escape upward and percolate through the flow. This auto-fluidization mechanism is efficient at all inclinations investigated (0-30°), suggesting that it could occur during the whole emplacement of a pyroclastic flow. A second study consisted of the vertical release of beds of particles in a static column. Results show that the granular mixture can be fully fluidized, even when collapsing from a relatively low height (20 cm). When particles are fine enough (<100 μm), pore pressure in the deposit diffuses for several seconds, the diffusion duration increasing with increasing bed thickness and decreasing particle size. The longest diffusion durations are observed for pyroclastic flow deposit materials (~30 s for 28.5 cm thick beds). These results suggest that pyroclastic flows propagating on irregular terrains can auto-fluidize and preserve low internal friction during their emplacement.

Le géoradar est une méthode électromagnétique haute fréquence (>10 MHz) utilisé pour caractériser les premiers mètres du sous-sol. Lors de la présence d'une topographie, les données géoradar sont déformées en conséquence. Afin de retrouver la vraie géométrie des réflecteurs, nous avons codés un algorithme de migration prenant en compte la topographie. La méthode est démontrée grâce à un modèle synthétique simple, puis testée avec succès sur des données réelles. Les algorithmes de migration apportent cependant du bruit dans les données. Pour pallier à ce problème, deux méthodes ont été mises en place : la première, inhérente à la migration, permet de réduire l'aliasing dit sur l'opérateur. La deuxième est un filtre ré-interpolant les traces en se basant sur un profil de pendage. Les deux méthodes suppriment un bruit incohérent des données mais dégradent les profils lorsqu'utilisées abusivement. Dans un deuxième chapitre, nous avons appliqués avec succès le géoradar dans un contexte de paléo-sismologie en Mongolie. L'utilisation conjointe de deux fréquences (50 et 500 MHz) ainsi que des comparaisons avec des tranchées a permis d'obtenir des informations complémentaires sur les géométries et les déplacements potentiels le long de deux failles. Dans un dernier chapitre, nous avons appliqués les mesures géoradar sur les dépôts pyroclastiques du volcan Tungurahua en Equateur. A nouveau, l'utilisation jointe de différentes fréquences (250, 500 et 800 MHz) nous permet d'imager efficacement les dépôts. Les unités principales sont mises en évidence avec l'antenne de 250 MHz et les architectures des dépôts sont observables avec les antennes de 500 et 800 MHz. / Georadar is a high frequency (>10MHz) electromagnetic method used to prospect near surface. When a topography is present, GPR images are distorted. To restore the true geometry of reflexions, we coded an migration algorithm which takes the topography into account. The method is first demonstrate on a simple synthetic model, and then succesfully applied on real data. However, migration algorithms bring noise to the data. Two methods have then been tested to avoid and remove it. The first one is inherent to the migration algorithm and reduce what is called operator's aliasing. The second one is a filter re-interpolating traces based on a profile containing the slope. Both methods remove inconsistent noise when used with caution, but decrease their quality when used with excess: reflexions presenting dip are the first to be deteriorated, as well as reflexions below strong topography. In a second chapter, we successfully used GPR in a paleo-sismology context in Mongolia. The use of two frequencies (50 and 500 MHz) as well as comparison with trenches bring complementaries informations on the geometry and possible offset along two faults. In the last chapter, GPR was tested over pyroclastic deposits from the Tungurahua volcano in Ecuador. Again, the combination of several frequencies (250, 500 and 800 MHz) has proven its efficiency. Main units were obvious with the 250 MHz antenna while the inner architecture of deposits was visible with the 500 and 800 MHz antenna.

Les écoulements (ou coulées) pyroclastiques (PFs) sont des mélanges concentrés de gaz et de particules à haute température qui représentent l’aléa volcanique le plus meurtrier qui soit. La protection des populations nécessite la mise au point de cartes des menaces précises, qui requièrent une connaissance fine de ces phénomènes. Cependant, les causes et les conséquences de l’érosion/incorporation associée aux PFs au cours de leur mise en place restent encore largement méconnues. Cette thèse se propose de caractériser la capacité érosive des PFs, de définir des mécanismes d’érosion, et de quantifier leurs impacts sur les budgets éruptifs et sur l’aléa associé. Pour cela, cette étude se concentre sur les PFs de volumes modestes mis en place pendant l’éruption d’août 2006 du volcan Tungurahua (Equateur) et adopte une démarche double, basée sur des investigations sédimentologiques et texturales des dépôts, couplées à la modélisation numérique. Une méthode originale, basée sur l’analyse d’images haute résolution corrigées par stéréologie, sur des études texturales détaillées des dépôts, et sur des bilans massiques de matière, permet de déterminer la granulométrie, la composition lithologique et la morphologie des produits sur l’ensemble de leur gamme de taille. Le calcul des bilans de matière montre que près de 50 wt. % des dépôts de PFs sont composés de matériaux non-juvéniles incorporés lors de la mise en place. Ces derniers proviennent principalement de la partie supérieure du volcan. La pente est ainsi le paramètre contrôlant au premier ordre l’intensité de l’érosion. Les budgets éruptifs complets indiquent un VEI de 3 (0,09 km 3 ) pour l’éruption, et soulignent l’importance de considérer séparément les matériaux juvéniles et non-juvéniles pour estimer la taille d’une éruption. L’étude détaillée des constituants met en évidence une ségrégation dynamique des clastes par densité au cours du transport, avec un taux de sédimentation de ≈10 cm.s -1 . Les données lithologiques, granulométriques et morphologiques démontrent la présence de phénomènes de fragmentation-abrasion des clastes pendant leur transport. Les clastes massifs (ex : laves anciennes) sont le principal agent de fragmentation des clastes scoriacés (ex : bombes). Des populations granulométriques fines, capables d’être transférées depuis l’écoulement dense principal vers les déferlantes et/ou le panache co-pyroclastique sont produites pendant toute la durée de la mise en place. Les modèles numériques basés sur une nouvelle loi d’érosion développée ici (et intégrée au code VolcFlow), démontrent la capacité de la rhéologie plastique à reproduire des PFs érosifs. L’érosion est associée à des variations dynamiques du rapport des contraintes normales/cisaillantes pendant la mise en place des écoulements, provoquées par des fluctuations d’épaisseur lors de phases de décélération. Le front fin des PFs, fortement frictionnel et érosif, est poussé par une tête et un corps plus épais, tous deux non érosifs. L’incorporation s’accompagne d’une augmentation de distance de parcours de l’ordre de 10-30% en fonction du taux d’incorporation, qui dépend de la quantité de matière affouillable disponible sur le volcan avant l’éruption. Ces résultats montrent que l’érosion peut avoir un rôle majeur sur les zones impactées par les PFs, et soulignent l’importance de prendre en compte cette capacité lors de la définition de l’aléa, ainsi que pour les études futures. / Pyroclastic flows (PFs) are hot mixtures of gas and particles that represent the most deadly volcanic hazard. To protect the populations, it is necessary to work on precise risk maps, which require having a deep knowledge of these phenomena. However, the causes and consequences of erosion and incorporation of non-juvenile material during PFs emplacement remain poorly known. This thesis aims at characterizing the erosive capacity of pyroclastic flows, defining erosion mechanisms and quantifying their impact on eruptive budgets and associated hazards. Here, we focus on small-volume PFs and use an approach based on field and textural investigations coupled with numerical modeling of PFs emplacement. The August 2006 PF-forming eruption of Tungurahua volcano (Ecuador) is used as a case-study for this work.An original method, based on high-resolution, stereologically-corrected image analyses, detailed textural analyses of PFs deposits and mass budget, enables determining the grain size distribution and the componentry of PFs products along their entire clast size range. Volume calculation and mass budgets show that about 50 wt. % of the whole deposit consists of non-juvenile materials incorporated during PFs emplacement, and mostly coming from the upper part of the volcano. The slope is a prevailing parameter that controls PFs erosive power. Eruptive budgets support a VEI 3 event (0.09 km 3 ) for the 2006 eruption of Tungurahua and highlight the importance of separating juvenile from non-juvenile material. Detailed analyses of deposits’ componentry suggest a strong dynamic density-driven segregation of the clasts during PFs emplacement, associated with sedimentation rates of ≈10 cm.s -1 . Lateral variations of lithological, grain size, and morphological data demonstrate the occurrence of componentry-driven clast fragmentation and abrasion processes. Massive components (e.g. old lavas) are the main grinding agents of scoriaceous components (e.g. bombs). During emplacement, these processes continuously create fine grained populations, which are transferred from the main dense flow to pyroclastic surge or Co-PF cloud. Numerical models of erosive PFs based on a new erosion law integrated into VolcFlow code show the ability of plastic rheology to reproduce natural erosion patterns of PFs. The erosion is produced by dynamic variations of normal stress / shear stress ratio during emplacement, due to thickness unsteadiness during flow deceleration. The thin, highly frictional and erosive front of PFs pulses is pushed by the thicker and non-erosive head and flow body. Incorporation implies longer PFs runouts of about 10-30%, depending on the amount of incorporated material, which is related to the quantity of erodible material available on the volcano’s flanks before the eruption. These results show that erosion has a significant role on PFs runouts, and thus in hazard assessment, which should be closely taken into account in future works.

La région de Saveh appartient dans sa presque totalité à la Zone de l'Iran Centra!. Nous y avons toutefois distingué, dans la partie sud-ouest, une zone présentant les caractères de la Zone de Sanandaj-Sirjan que nous appelons Zone de Saman. La Zone de Saman, prolongement oriental de la Zone de Razan ( Thèse de M. Bolourchi 1975), vient former un coin dans la Zone de l'Iran Central, limité au Nord par la faille d'Avaj (Bolourchi 1975) et au Sud par la faille de Pavan, encore appelée "Indes Fault" plus au Sud. La Zone de l'Iran Central est caractérisée par son magmatisme tertiaire, et surtout par son volcanisme éocène. Amidi ( thése 1975), lui attribue le nom de Zone de Tabriz-Torbat-JaZ Murian, provinces situées à ses extrémités. La tectonique cassante y est la règle et la carte de l'Iran au 1/2 500 000 montre une mosaïque de blocs séparés par des failles. La Zone de Sanandaj-Sirjan, au contraire des précédentes, montre une série métamorphique épaisse. Les différences essentielles avec les zones précédentes sont les suivantes: - existence d'une phase plutonique importante à la fin du Mésozoïque; - plusieurs phases de métamorphisme anté-tertiaires (paléocimmérienne, Sabzehei 1974 - laramienne, Berberian 1972); - faible importance des formations tertiaires.

L'effet de la végétation locale, composée de Castanea sativa cultivé, sur la stabilité des pentes a été étudié sur un site d'essai au Mont Faito (Campanie, Italie). En Campanie, les sols pyroclastiques peu profonds sont sensibles aux glissements de terrain provoqués par les précipitations. Des périodes de pluies prolongées suivies de précipitations extrêmes à court terme déclenchent des glissements de terrain rapides et destructeurs au niveau des coupes routières et des escarpements pyroclastiques sur les falaises rocheuses dans les régions autour du volcan Vésuve.Des échantillons de sol pyroclastiques non perturbés contenant des racines de C. sativa matures ont été utilisés pour la caractérisation hydraulique par le biais d'un ensemble d'expériences en laboratoire. La perméabilité saturée, la réponse à l’évaporation et l’imbibition, la teneur en eau pour les fortes valeurs de succion et la biomasse sèche des racines ont été déterminées.La présence de racines a augmenté la perméabilité d'un ordre de grandeur dans les sols les plus superficiels (10-7 à 10-6 m s-1) et diminué la valeur d'entrée d'air des courbes de rétention (6 à 4 kPa). La variabilité de la perméabilité entre les couches de sol a été identifiée comme conditionnant l'écoulement de l'eau souterraine par rapport à la vitesse du mouvement du front de mouillage et à la génération de pressions positives de l'eau interstitielle dans le profil. L'étalonnage du modèle hystérétique pour caractériser les sols pyroclastiques naturels a fourni une méthode plus approximative de modélisation des réponses hydrauliques. Une bonne concordance entre le modèle et les observations a été obtenue.L’étude sur le terrain a permis de montrer que la distribution des racines de C. sativa est associée au régime des eaux souterraines. Les distributions spatiales et verticales de la densité et des traits des racines ont été quantifiées pour les racines de C. sativa prélevées dans des forages réalisés au Mont Faito. La succion minimale, la teneur minimale en eau et la pente minimale (indiquant un débit d'eau descendant) ont été surveillées tout au long de l'année et confrontées avec la distribution des racines et à la distribution spatiale des arbres. Une densité racinaire croissante était associée à des valeurs de succion plus faibles et à des gradients d'infiltration plus élevés, ce qui peut avoir une influence négative sur la stabilité de la pente.La modélisation du renforcement mécanique du sol par les racines des arbres a permis de comprendre l'importance des composantes hydrauliques et mécaniques sur la stabilité d'une pente. Les racines augmentent la résistance au cisaillement (jusqu'à 25,8 kPa) grâce à un renforcement mécanique et donc le facteur de sécurité de la pente augmente. L'examen du renforcement dû aux racines dans l'estimation du facteur de sécurité des surfaces de rupture potentielles a montré que la surface de rupture la plus faible a été trouvée à 2,2 m, où le renforcement dû aux racines était de 1,3 kPa, au lieu de 0,9 m sans le renforcement de 13,8 kPa. La surface de rupture la plus faible correspond aux surfaces de rupture observées lors de glissements de terrain antérieurs. Le site d'essai ne présentait pas les caractéristiques d'une zone de déclenchement d'un glissement de terrain. L'angle de pente des zones de déclenchement des glissements de terrain (35° à 45°) peut dépasser l'angle de frottement du sol (36,5° à 38,5°) et l'effet hydraulique ne serait pas suffisant pour garantir la stabilité de la pente pendant la saison humide (0 à 10 kPa). On estime que le renforcement dû aux racines peut maintenir les pentes jusqu'à un angle de 42°.On a donc constaté que la présence de racines d'arbres affectait la stabilité hydraulique et mécanique des couvertures de sol pyroclastiques. Ces conclusions peuvent être étendues aux autres zones de plantations de C. sativa. L'effet hydraulique de la végétation a été largement compensé par le renforcement mécanique dû aux racines / The effect of the local vegetation, composed of cultivated Castanea sativa, on slope stability was investigated on a test site in Mount Faito (Campania, Southern Italy). In Campania, shallow pyroclastic soil covers are susceptible to landslides triggered by rainfall. Prolonged rainfall periods followed by extreme short-term rainfall events trigger fast moving and highly destructive landslides in road cuts and pyroclastic scarps on rocky cliffs in the areas surrounding the Vesuvius volcano.Undisturbed pyroclastic soil samples containing roots of mature C. sativa were used for hydraulic characterization through an extensive set of laboratory experiments. Saturated permeability, evaporation and imbibition response, water content for high suction ranges, and the root dry biomass were determined.The presence of roots increased the hydraulic permeability by one order of magnitude in the most surficial soil (10-7 to 10-6 m s-1) and decreased the air-entry value of the water retention curves (6 to 4 kPa). The variability of soil permeability among soil layers was identified as conditioning of the groundwater flow with regard to the speed of the wetting front movement and generation of positive pore-water pressures within the soil profile. The calibration of hysteretic model to characterize natural pyroclastic soil provided a more approximate manner of modelling in situ hydraulic responses. A good agreement between the model and the field observations was obtained.Field monitoring was performed with the intent of showing that the distribution of roots of C. sativa is associated to the groundwater regime. The spatial and vertical distribution of root density and traits were quantified for C. sativa roots collected from several boreholes performed in Mount Faito. Minimum suction, minimum water content and minimum gradient (indicative of downward water flow), were monitored throughout the year and related to root distribution and spatial distribution of trees. An increasing root density was found to be associated to lower values of suction and higher gradients of infiltration, which can potentially have a negative influence of the slope stability.A modelling investigation on the mechanical reinforcement of soil by tree roots allowed us to understand the importance of hydraulic and mechanic components on the stability of a slope. Roots increase greatly the shear strength of soil (up to 25.8 kPa) through mechanical reinforcement and consequently, the safety factor of the slope increased significantly. Considering the root reinforcement in the estimation of potential failure surfaces safety factor showed that the weakest failure surface was found at 2.2 m, where the root reinforcement was 1.3 kPa, instead of 0.9 m without the root reinforcement of 13.8 kPa. The weakest failure surface found was in agreement with the failure surfaces observed from previous landslides. The test site did not present the characteristics of a landslide triggering area. The slope angle of the landslide triggering areas (35° to 45°) can easily exceed the soil friction angle (36.5° to 38.5°) and the hydraulic effect would not be enough to guarantee the stability of the slope during the wet season (0 to 10 kPa). However, the root reinforcement was estimated to be able to sustain the slopes until an angle of 42°.Therefore, the presence of tree roots was found to affect hydraulically and mechanically stability of pyroclastic soil covers. Such conclusions may be extended to the areas of Campania where C. sativa plantations are present. The hydraulic effect of vegetation was greatly compensated by the mechanical reinforcement of roots.

Durant les trois années de la thèse, j’ai eu le plaisir de travailler en collaboration avec à la fois des volcanologues, des physiciens de laboratoire et des mathématiciens. Ce mémoire est l’occasion de présenter la démarche et les résultats de mes recherches dans le domaine de la modélisation d’écoulements granulaires denses fluidisés. Ces derniers consistent à développer un nouveau modèle mathématique et son étude théorique et numérique. Sur la base d’observations faites lors d’expériences de laboratoire, nous proposons une façon de modéliser le changement comportemental d’un écoulement granulaire initialement fluidisé au travers de la définition de sa rhéologie viscoplastique à seuil variable. Plus précisément, le seuil de plasticité est défini par la différence entre la pression lithostatique et la pression du fluide interstitiel. La nouveauté apportée par ce modèle ouvre de nouvelles perspectives à la fois pour le champ de recherche en mathématiques et pour la compréhension des lits granulaires fluidisés et leur application à la volcanologie. Du point de vue mathématique, une étude théorique du modèle a été menée. En proposant une preuve de l’existence de solutions faibles à un problème lié à la version homogène du modèle, nous apportons une extension au champ de connaissances autour des écoulements des fluides non-newtoniens. D’autre part, dans le but de reproduire numériquement des expériences de laboratoire de chute de colonne granulaire fluidisée, nous avons développé un code de simulation numérique incluant une nouvelle méthode de résolution des équations d’écoulement de fluides à seuil. Dans ce manuscrit, je décris et justifie les différents choix stratégiques pour le développement de ce code. Par ailleurs, je présente quelques tests académiques permettant de valider le code. Enfin, je donne les résultats de simulation de chute de colonne granulaire, qu’elle soit fluidisée ou non. Une comparaison avec les données de laboratoire est effectuée afin d’évaluer les points forts et les défauts du modèle par rapport à la réalité des expériences. En conclusion, dans la continuité du travail mené dans ce projet, des perspectives d’amélioration sont proposées. / During these three years, I enjoyed to work with collaborators from volcanology, laboratory physics and mathematics. This document presents the steps and results of my research in the field of modelling of fluidised granular flows. The last consists in the development of a new mathematical model and its theoretical and numerical study. Based on observations made on experimental studies, the model focuses on the change in the behaviour of an initially fluidised granular flow through the definition of its viscoplastic rheology with variable threshold. More precisely, the threshold (aslo called yield stress) is defined via the difference between the lithostatic pressure and the pressure of the interstitial fluid. The innovation of this model opens perspectives for the mathematical research as well as for the study of fluidised granular flows and their application to volcanology. From a mathematical point of view, a theoretical study has been conducted. Proving the existence of weak solution for the homogeneous version of the model, we offer an extension in the field of knowledges of non-newtonian fluid flows. Also, we have developped a numerical code to simulate dambreak experiments with fluidised granular media. This one includes a new method to solve the flow equations of viscoplastic fluids. In this thesis, I describe and justify the numerical strategy chosen. Moreover, I present some academic tests to validate the code. At the end, I give the numerical results in the case of the dambreak simulation for dry and fluidised fluids. By comparing with experimental data, we evaluate the validity of the model and its resolution, and highlight the advantages and inconvenients. To conclude the project, I propose some perspectives of improvement for later work.