Benchmarked computation of time-domain viscoelastic Love numbers for adiabatic mantles

Abstract

The viscoelastic load Love numbers encapsulate the Earth’s rheology in a remarkably efficient fashion. When multiplied by a sudden increment of spherical harmonic load change, they give the horizontal and vertical surface displacements and gravity change at all later times. Incremental glacial load changes thus need only be harmonically decomposed, multiplied by the Love numbers and summed to predict the Earth’s response to glacial load redistributions. The computation of viscoelastic Love numbers from the elastic, viscous and adiabatic profiles of the Earth is thus the foundation upon which many glacial isostatic adjustment models are based. Usually, viscoelastic Love numbers are computed using the Laplace transform method, employing the correspondence principle to convert the viscoelastic equations of motion into the elastic equations with complex material parameters. This method works well for a fully non-adiabatic Earth, but can accommodate realistic partially adiabatic and fully adiabatic conditions only by changing the Earth’s density profile. An alternative method of Love number computation developed by Cathles (1975) avoids this dilemma by separating the elastic and viscous equations of motion. The separation neglects a small solid-elastic/fluidelastic transition for compressible deformation, but allows freely defining adiabatic, partially adiabatic or fully non-adiabatic profiles in the mantle without changing the Earth’s density profile. Here, we update and fully describe this method and show that it produces Love numbers closely similar to those computed for fully non-adiabatic earth models computed by the correspondence principle, finite element and other methods. The time-domain method produces Love numbers as good as those produced by other methods and can also realistically accommodate any degree of mantle adiabaticity. All method implementations are available open source.