Small-scale convection in a plume-fed low-viscosity layer beneath a moving plate

Two-dimensional simulations using a thermomechanical model based on a finite-difference method on a staggered grid and a marker in cell method are performed to study the plume-lithosphere interaction beneath moving plates. The plate and the convective mantle are modelled as a homogeneous peridotite with a Newtonian temperature- and pressure-dependent viscosity. A constant velocity, ranging from 5 to 12.5 cm yr−1, is imposed at the top of the plate. Plumes are generated by imposing a thermal anomaly of 150 to 350 K on a 50 km wide domain at the base of the model (700 km depth); the plate atop this thermal anomaly is 40 Myr old. We analyse (1) the kinematics of the plume as it impacts the moving plate, (2) the dynamics of time-dependent small-scale convection (SSC) instabilities developing in the low-viscosity layer formed by spreading of hot plume material at the base of the lithosphere and (3) the resulting thermal rejuvenation of the lithosphere. The spreading of the plume material at the base of the lithosphere, characterized by the ratio between the maximum down- and upstream horizontal (dimensionless) velocities in the plume-fed sublithospheric layer, Peup/Pedown depends on the ratio between the maximum plume upwelling velocity and the plate velocity, Peplume/Peplate. For fast plate velocities and sluggish plumes (low Peplume/Peplate), plate motion drags most plume material and downstream flow is dominant. As Peplume/Peplate increases, an increasing part of the plume material flows upstream. SSC systematically develops in the plume-fed sublithospheric layer, downstream from the plume. Onset time of SSC decreases with the Rayleigh number. For vigorous plumes, it does not depend on plate velocity. For more sluggish plumes, however, variations in the plume spreading behaviour at the base of the lithosphere result in a decrease in the onset time of SSCs with increasing plate velocity. In any case, SSC results in uplift of the isotherm 1573 K by up to 20 km relative to its initial equilibrium depth at the impact poin

Extremely thin crust in the Indian Ocean possibly resulting from Plume–Ridge Interaction

International audienceThe thickness of the crust created at ocean spreading centres depends on the spreading rate and melt production in the mantle. It is ~5–8 km for a crust formed at slow and fast spreading centres and 2–4 km at ultra-slow spreading centres away from hotspots and mantle anomalies. The crust is generally thin at the fracture zones and thick beneath hotspots and large igneous provinces. Here we present results for the crust generated at the fast Wharton spreading centre 55–58 Ma ago using seismic reflection and refraction data. We find that the crust over a 200 km segment of the Wharton Basin is only 3.5–4.5 km thick, the thinnest crust ever observed in a fast spreading environment. A thin crust could be produced by the presence of depleted and/or cold mantle. Numerical simulations and recent laboratory experiments studying the impact of a hot plume under a lithosphere show that a curtain of weak cold downwellings of depleted mantle material is likely to develop around the edges of the hot plume pond. Because of a strongly temperature-dependent viscosity of lithospheric material, the hotter, therefore less viscous, bottom of the lithosphere can be mobilized by an impinging plume. If sampled by a spreading centre, the locally cold and depleted mantle should result in low production of melt. We suggest that the observed thin crust in the Wharton Basin is likely to have been formed by the interaction between the Kerguelen mantle plume and the Wharton spreading centre ~55 Ma ago

Mesure des champs de vitesse dans les champs hydrothermaux océaniques

Nous présentons une nouvelle méthode pour déterminer les champs de vitesse dans les champs hydrothermaux présents sur les fonds sous-marins. Basée sur le "background-oriented schlieren", la méthode permet de suivre les anomalies d'indice de réfraction dues à la température ou la composition. Elle a été testée sur des panaches thermiques en laboratoire, puis lors d'une campagne sur la dorsale atlantique

A noninvasive method for measuring the velocity of diffuse hydrothermal flow by tracking moving refractive index anomalies

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/95358/1/ggge1802.pd

Starting laminar plumes: Comparison of laboratory and numerical modeling

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/94953/1/ggge1631.pd

Quantifying diffuse and discrete venting at the Tour Eiffel vent site, Lucky Strike hydrothermal field

Author Posting. © American Geophysical Union, 2012. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry Geophysics Geosystems 13 (2012): Q04008, doi:10.1029/2011GC003991.The relative heat carried by diffuse versus discrete venting of hydrothermal fluids at mid-ocean ridges is poorly constrained and likely varies among vent sites. Estimates of the proportion of heat carried by diffuse flow range from 0% to 100% of the total axial heat flux. Here, we present an approach that integrates imagery, video, and temperature measurements to accurately estimate this partitioning at a single vent site, Tour Eiffel in the Lucky Strike hydrothermal field along the Mid-Atlantic Ridge. Fluid temperatures, photographic mosaics of the vent site, and video sequences of fluid flow were acquired during the Bathyluck'09 cruise (Fall, 2009) and the Momarsat'10 cruise (Summer, 2010) to the Lucky Strike hydrothermal field by the ROV Victor6000 aboard the French research vessel the “Pourquoi Pas”? (IFREMER, France). We use two optical methods to calculate the velocities of imaged hydrothermal fluids: (1) for diffuse venting, Diffuse Flow Velocimetry tracks the displacement of refractive index anomalies through time, and (2) for discrete jets, Particle Image Velocimetry tracks eddies by cross-correlation of pixel intensities between subsequent images. To circumvent video blurring associated with rapid velocities at vent orifices, exit velocities at discrete vents are calculated from the best fit of the observed velocity field to a model of a steady state turbulent plume where we vary the model vent radius and fluid exit velocity. Our results yield vertical velocities of diffuse effluent between 0.9 cm s−1 and 11.1 cm s−1 for fluid temperatures between 3°C and 33.5°C above that of ambient seawater, and exit velocities of discrete jets between 22 cm s−1 and 119 cm s−1 for fluid temperatures between 200°C and 301°C above ambient seawater. Using the calculated fluid velocities, temperature measurements, and photo mosaics of the actively venting areas, we calculate a heat flux due to diffuse venting from thin fractures of 3.15 ± 2.22 MW, discrete venting of 1.07 ± 0.66 MW, and, by incorporating previous estimates of diffuse heat flux density from Tour Eiffel, diffuse flux from the main sulfide mound of ∼15.6 MW. We estimate that the total integrated heat flux from the Tour Eiffel site is 19.82 ± 2.88 MW and that the ratio of diffuse to discrete heat flux is ∼18. We discuss the implication of these results for the characterization of different vent sites within Lucky Strike and in the context of a compilation of all available measurements of the ratio of diffuse to discrete heat flux.E. Mittelstaedt was supported by the International Research Fellowship Program of the U.S. National Science Foundation (OISE-0757920). Funding for the 2006, 2008, 2009, and 2010 cruises was provided by CNRS/ IFREMER through the MoMAR program (France), by ANR (France), the Mothseim Project NT05–3 42213 to J. Escartín and by grant CTM2010–15216/MAR from the Spanish Ministry of Science to R. Garcia and J. Escartín. T. Barreyre was supported by University Paris Diderot (Paris 7 – France) and Institut de Physique du Globe de Paris (IPGP, France).2012-10-1

Drying colloidal systems: laboratory models for a wide range of applications

The drying of complex fluids provides a powerful insight into phenomena that take place on time and length scales not normally accessible. An important feature of complex fluids, colloidal dispersions and polymer solutions is their high sensitivity to weak external actions. Thus, the drying of complex fluids involves a large number of physical and chemical processes. The scope of this review is the capacity to tune such systems to reproduce and explore specific properties in a physics laboratory. A wide variety of systems are presented, ranging from functional coatings, food science, cosmetology, medical diagnostics and forensics to geophysics and art

Dynamique des panaches thermiques laminaires (application aux panaches mantelliques)

En géosciences, les points chauds, îles volcaniques intraplaques telles que Hawaii ou la Réunion, pourraient être l expression en surface de panaches qui monteraient depuis la base du manteau, à 2900 km de profondeur. Les panaches, instabilités de la convection de Rayleigh-Bénard, se développeraient dans le manteau avec un nombre de Prandtl (Pr) de 1023 mais peuvent aussi être modéisés à partir d un chauffage de petite taille, numériquement avec Pr infini, ou en laboratoire avec Pr ~103 - 106. Dans cette thèse nous étudions les caractéristiques d un panache isolé dans un fluide visqueux, à viscosité constante. Nous utilisons des expériences de laboratoire et des modèles numériques. Les techniques de visualisation donnent accès aux champs de température et de vitesse pour une section 2-D de la cuve. Les simulations numériques utilisent une méthode d éléments finis à symétrie cylindrique qui reproduisent les conditions de laboratoire avec les propriétés mesurées des fluides et du chauffage. On obtient un excellent accord entre ces deux approches indépendantes. Cela permet de proposer des lois d échelles simples pour la dynamique de la tête et du conduit du panache et de les appliquer au cas de la Terre. Nous montrons notamment que pour des Pr supérieurs à 1000 l effet de confinement l emporte largement sur les effets inertiels et que la dynamique du panache est bien décrite par l approximation Pr infini. Ceci est en particulier vrai pour les manteaux planétairesIn Earth Sciences, hot upwelling plumes are thought to develop from the base of the 2900 kmthick solid mantle of our planet and to generate hotspots, i.e. intraplate volcanic islands such as Hawaii and La Reunion. Although generated through chaotic Rayleigh-B enard instabilities at a Prandtl number (Pr) around 1023, they can be modelled with the simpler case of starting plumes out of a finite-size heater, either numerically for infinite Prandtl number, or in the laboratory with fluids with Pr ~ 103 - 106. Hence, the question is to find simple scaling laws for isolated rising plumes and apply them to the Earth s mantle case. In this thesis we study the characteristics of an isolated plume growing in a viscous fluid with constant viscosity. We use both laboratory experiments and numerical models : the visualization techniques give us access to the growing plume temperature and velocity fields, on a 2-D section of the tank, whereas the numerical simulations are axisymmetric finite element simulations that attempt to reproduce the laboratory conditions as closely as possible. We find excellent quantitative agreement between the two fully independent approaches. This is used to derive scaling laws for the dynamics of the plume head and stem, and apply them to the Earth s case. We further show that for Pr larger than 1000, confinement effects are more important than inertial effects and that plumes dynamics are well described by the approximation if infinite Prandlt number. This is especially true for planetary mantlesPARIS-BIUSJ-Sci.Terre recherche (751052114) / SudocSudocFranceF

Instabilités thermiques dans un fluide à seuil (de l'échelle du laboratoire à celle de la planète)

Des panaches sont connus pour migrer à travers le manteau ductile et quasi-Newtonien ; alors que les dikes se fracturent et se propagent dans la lithosphère solide. Cependant, la lithosphère est en fait visco-élastique. Afin de déterminer ce qui se passe dans ce cas complexe, nous avons réalisé une étude expérimentale et numérique sur le développement de panaches thermiques dans des solutions aqueuses de Carbopol, un gel de polymères formant un réseau continu d'éponges microscopiques. Ce fluide est rhéofluidifiant et présente un seuil de contrainte , de sorte que l'écoulement ne se produit que si la contrainte locale dépasse cette valeur critique . En dessous de cette valeur, le fluide agit comme un solide élastique. Les propriétés rhéologiques des solutions peuvent être systématiquement ajustées en variant la concentration de Carbopol. Le dispositif consiste en une source locale de chaleur de puissance constante placée au centre d'une cuve cubique. Selon la valeur du rapport entre la contrainte d'origine thermique et la contrainte seuil, Y , on peut observer trois régimes différents. A faible Y Yc . Nous montrons que les paramètres critiques (Yc ,Yc ) dépendent fortement de la géométrie du chauffage. Des mesures simultanées de température et de champs de vitesse montrent que la morphologie du panache ressemble à un doigt, contrairement à la forme de champignon rencontrée dans les fluides newtoniens. Utilisant des simulations numériques avec une description purement visqueuse, où la rhéologie du fluide est décrite par un modèle de Herschel-Bulkley régularisé, sont suffisantes pour rendre compte de la dynamique du panache. Une étude détaillée des paramètres indiquent que la dynamique du panache est gouvernée par la compétition entre la contrainte seuil, la contrainte induite par la flottabilité et les contraintes visqueuses. Nous avons identifié deux paramètres adimensionnés : le paramètre seuil comparant la contrainte induite par la flottabilité et la contrainte seuil, et le nombre de Bingham Bi comparant la contrainte seuil et les contraintes visqueuses. Un panache ne peut s'élever que si les deux paramètres sont supercritiques, i.e. la contrainte induite par la flottabilité et les contraintes visqueuses sont plus importantes que la contrainte seuil. Par conséquent, le panache peut s'arrêter avant d'atteindre la surface. Des lois d'échelles dans le conduit du panache ont été déterminées pour la vitesse, la température et la taille de la région cisaillée en régime permanent. Elles décrivent raisonnablement le comportement du conduit bien que seul l'effet rhéofluidifiant soit pris en compte. L'application de ces paramètres adimensionnés à la Terre contraignent significativement la limite de plasticité du manteau et de la lithosphère. La contrainte seuil maximale qui permet à une instabilité thermique de pénétrer dans la lithosphère ou le manteau supérieur est entre 100 kPa et 100 MPa, et elle dépend fortement de la taille et de l'anomalie de densité de l'intrusion.Plumes are known to migrate through the ductile quasi-Newtonian mantle, while dikes fracture and propagate through the solid lithosphere. However, depending on the timescale, the lithosphere presents solid as well as viscous properties. To determine what happens in the complex case, where instabilities propagate through a visco-elastic matrix, we performed a combined study of laboratory experiments and numerical simulations. Here we investigate the development of thermal plumes in aqueous solutions of Carbopol, a polymer gel, forming a continuous network of micrometric sponges. This fluid is shear thinning and has a yield-stress , whereby flow occurs only if the local stress exceeds this critical value . Below this value, the fluid acts as an elastic solid. The rheological properties of the solutions can be systematically varied by varying the Carbopol concentration. The setup consists of a localized heat-source operated at constant power, placed at the centre of a square tank. Depending on the ratio of the thermally induced stresses and the yield stress, Y , three different regimes may be obtained. For low Y Yc . We show that the critical parameters (Yc ,Yc ) strongly depend on the geometry of the heating. Combined temperature and velocity field measurements show that the morphology of the plume resembles a finger, contrary to the mushroom-like shape encountered in Newtonian fluids. Numerical simulations using a purely viscous description, where the rheology of the fluid is described by a regularized Herschel-Bulkley model, are sufficient to capture the plume dynamics. A detailed parametric study shows that the plume dynamics are governed by the interplay between yield stress, buoyancy induced stress and viscous stresses. We identify two non-dimensional parameters: the yield parameter comparing the buoyancy induced stress to the yield stress, and the Bingham number Bi comparing the yield stress to the viscous stresses. We show that a plume can rise only if both parameters are supercritical, i.e. if buoyancy induced stress and viscous stresses each overcome the yield stress. Therefore the plume may come to a halt before it reaches the surface. We propose scaling laws for the plume stem velocity, temperature and the size of the shear zone in the steady state. We show that the scaling laws describe the behaviour in the plume stem reasonably well, if the yield stress is neglected and only the shear thinning behaviour is taken into account. Applying the non-dimensional parameters to Earth places severe constraints on the strength of mantle and lithosphere. The maximum strength that allows for thermal instabilities to penetrate the lithosphere or upper mantle is in between 100 kPa and 100 MPa, and strongly depends on the size and buoyancy of the anomaly.PARIS11-SCD-Bib. électronique (914719901) / SudocSudocFranceF