摘要

The dispersion relation of the Rayleigh Taylor instability, a gravitational instability associated with unstable density stratification, is of profound importance in various geophysical contexts. When more than two layers are involved, a semi-analytical technique based on the biharmonic formulation of Stokes flow has been extensively used to obtain such dispersion relation. However, this technique may become cumbersome when applied to lithospheric dynamics, where a number of layers are necessary to represent the continuous variation of viscosity over many orders of magnitude. Here, we present an alternative and more efficient method based on the propagator matrix formulation of Stokes flow. With this approach, the original instability problem is reduced to a compact eigenvalue equation whose size is solely determined by the number of primary density contrasts. We apply this new technique to the stability of the early crust, and combined with the Monte Carlo sensitivity analysis, we derive an empirical formula to compute the growth rate of the Rayleigh Taylor instability for this particular geophysical setting. Our analysis indicates that the likelihood of crustal delamination hinges critically on the effective viscosity of eclogite.

  • 出版日期2018-3

全文