Thermal instabilities in a soft and complex lithosphere: laboratory experiments and numerical simulations
Abstract
The upwelling of hot material in the lithosphere remains far from understood. This is due to the complexity of the mechanical behaviour of lithospheric material, which presents solid as well as viscous properties. Mushroom-shaped less viscous plumes or more viscous finger-shaped diapirs, depending on the viscosity ratio between the rising and the matrix materials, are known to migrate through ductile, quasi-newtonian lithosphere; while dikes fracture and propagate through a solid matrix. But what happens in between these two end-members?
To answer this question, we perform a combined study of laboratory experiments and numerical simulations on the development of thermal plumes in aqueous solutions of Carbopol, a polymer gel suspension forming a continous network of micrometric sponges. This fluid is shear thinning and presents a yield-stress, whereby flow occurs only if the local stress exceeds a critical value. Below this value, the fluid acts as an elastic solid. Our experimental setup consists of a localized heat-source, placed in the center of a squared plexiglas tank. At t=0, a constant thermal power is applied locally to the fluid. For the numerical simulations, we replace the rigid plastic regions by an extremely viscous fluid, and therefore neglect the elastic contribution to the local stress. We systematically studied the influence of the rheological parameters, as well as the supplied heat. Depending on the Yield number Y0, which compares the thermally-induced stress to the yield stress, three different regims are observed. For low Y0, no convection develops; while for intermediate values, a small-scale convection cell appears and remains confined around the heater. For high Y0, thermal instabilities rise through the tank. Their morphology differs from the mushroom-shape typically encountered in newtonian fluids. Combined temperature and velocity field measurements show that a plug flow develops within the plume thermal anomaly, therefore producing a rising finger-shape with strong shear zones confined along its edges. The characteristics of the instability, as well as the existence of unyielded regions and the development of a damaged zone ahead of the plume as it rises, depend on Y0 but also on the other rheological parameters. The numerical simulations recover well the features observed in the laboratory experiments. This allows us to extend the parameter range of study. Our experimental finger-shaped diapirs present strong similarities with an off-axis diapir in Oman emplaced in a ridge context. This geological object, several kilometers in diameter presents in particular strong shear localization along its edges. Within our fluid mechanics framework, the existence of such an instability in the lithosphere places strong constraints on its parameter range. It suggests that this diapir was emplaced in a partially molten lithosphere. Therefore Herschel-Bulkley fluids like Carbopol might be good candidates to get new insights into the behavior of "soft" geological systems like mid-ocean ridge systems.- Publication:
-
AGU Fall Meeting Abstracts
- Pub Date:
- December 2012
- Bibcode:
- 2012AGUFM.T32C..04M
- Keywords:
-
- 8120 TECTONOPHYSICS / Dynamics of lithosphere and mantle: general;
- 8159 TECTONOPHYSICS / Rheology: crust and lithosphere;
- 8164 TECTONOPHYSICS / Stresses: crust and lithosphere