Source author record

Mikael Beuthe

Mikael Beuthe appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.geo-ph astro-ph.EP physics.space-ph physics.flu-dyn

Catalog footprint

What is connected

7works

4topics

2close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2021arXiv

Isostasy with Love: II Airy compensation arising from viscoelastic relaxation

In modern geodynamics, isostasy can be viewed either as the static equilibrium of the crust that minimizes deviatoric stresses, or as a dynamic process resulting from the viscous relaxation of the non-hydrostatic crustal shape. Paper I gave a general formulation of Airy isostasy as an elastic loading problem solved with Love numbers, and applied it to the case of minimum stress isostasy. In this sequel, the same framework is used to study Airy isostasy as the long-time evolution of a viscoelastic shell submitted to surface and internal loads. Isostatic ratios are defined in terms of time-dependent deviatoric Love numbers. Dynamic isostasy depends on the loading history, two examples of which are the constant load applied on the surface in the far past and the constant shape maintained by addition or removal of material at the compensation depth. The former model results in a shape decreasing exponentially with time and has no elastic analog, whereas the latter (stationary) model is equivalent to a form of elastic isostasy. Viscoelastic and viscous approaches are completely equivalent. If both load and shape vary slowly with time, isostatic ratios look like those of the stationary model. Isostatic models thus belong to two independent groups: the elastic/stationary approaches and the time-dependent approaches. If the shell is homogeneous, all models predict a similar compensation of large-scale gravity perturbations. If the shell rheology depends on depth, stationary models predict more compensation at long wavelengths, whereas time-dependent models result in negligible compensation. Mathematica and Fortran codes are available for computing the isostatic ratios of an incompressible body with three homogeneous layers.

preprint2019arXiv

Enceladus's crust as a non-uniform thin shell: II tidal dissipation

Tidal heating is the prime suspect behind Enceladus's south polar heating anomaly and global subsurface ocean. No model of internal tidal dissipation, however, can explain at the same time the total heat budget and the focusing of the energy at the south pole. I study here whether the non-uniform icy shell thickness can cause the north-south heating asymmetry by redistributing tidal heating either in the shell or in the core. Starting from the non-uniform tidal thin shell equations, I compute the volumetric rate, surface flux, and total power generated by tidal dissipation in shell and core. If the shell is laterally uniform, the thin shell approach predicts shell dissipation with a few percent error while the error on core dissipation is negligible. Variations in shell thickness strongly increase the shell dissipation flux where the shell is thinner. For a hard shell with long-wavelength variations, the shell dissipation flux can be predicted by scaling with the inverse local thickness the flux for a laterally uniform shell. If Enceladus's shell is in conductive thermal equilibrium with isostatic thickness variations, the nominal shell dissipation flux at the south pole is about three times its value for a shell of uniform thickness, which remains negligible compared to the observed flux. The shell dissipation rate should be ten times higher than nominal in order to account for the spatial variations of the observed flux. Dissipation in an unconsolidated core can provide the missing power, but does not generate any significant heating asymmetry as long as the core is homogeneous. Non-steady state models, though not investigated here, face similar difficulties in explaining the asymmetries of tidal heating and shell thickness.

preprint2016arXiv

Enceladus's and Dione's floating ice shells supported by minimum stress isostasy

Enceladus's gravity and shape have been explained in terms of a thick isostatic ice shell floating on a global ocean, in contradiction of the thin shell implied by librations. Here we propose a new isostatic model minimizing crustal deviatoric stress, and demonstrate that gravity and shape data predict a $\rm{38\pm4\,km}$-thick ocean beneath a $\rm{23\pm4\,km}$-thick shell agreeing with -- but independent of -- libration data. Isostatic and tidal stresses are comparable in magnitude. South polar crust is only $7\pm4\rm\,km$ thick, facilitating the opening of water conduits and enhancing tidal dissipation through stress concentration. Enceladus's resonant companion, Dione, is in a similar state of minimum stress isostasy. Its gravity and shape can be explained in terms of a $\rm{99\pm23\,km}$-thick isostatic shell overlying a $\rm{65\pm30\,km}$-thick global ocean, thus providing the first clear evidence for a present-day ocean within Dione.

preprint2015arXiv

Tidal Love numbers of membrane worlds: Europa, Titan, and Co

Under tidal forcing, icy satellites with subsurface oceans deform as if the surface were a membrane stretched around a fluid layer. `Membrane worlds' is thus a fitting name for these bodies and membrane theory provides the perfect toolbox to predict tidal effects. I describe here a new membrane approach to tidal perturbations based on the general theory of viscoelastic-gravitational deformations of spherically symmetric bodies. The massive membrane approach leads to explicit formulas for viscoelastic tidal Love numbers which are exact in the limit of zero crust thickness. Formulas for load Love numbers come as a bonus. The accuracy on $k_2$ and $h_2$ is better than one percent if the crust thickness is less than five percents of the surface radius, which is probably the case for Europa and Titan. The new approach allows for density differences between crust and ocean and correctly includes crust compressibility. This last feature makes it more accurate than the incompressible propagator matrix method. Membrane formulas factorize shallow and deep interior contributions, the latter affecting Love numbers mainly through density stratification. I show that a screening effect explains why ocean stratification typically increases Love numbers instead of reducing them. For Titan, a thin and dense liquid layer at the bottom of a light ocean can raise $k_2$ by more than ten percents. The membrane approach can also deal with dynamical tides in a non-rotating body. I show that a dynamical resonance significantly decreases the tilt factor and may thus lead to underestimating Europa's crust thickness. Finally, the dynamical resonance increases tidal deformations and tidal heating in the crust if the ocean thickness is of the order of a few hundred meters.

preprint2014arXiv

Tides on Europa: the membrane paradigm

Jupiter's moon Europa has a thin icy crust which is decoupled from the mantle by a subsurface ocean. The crust thus responds to tidal forcing as a deformed membrane, cold at the top and near melting point at the bottom. In this paper I develop the membrane theory of viscoelastic shells with depth-dependent rheology with the dual goal of predicting tidal tectonics and computing tidal dissipation. Two parameters characterize the tidal response of the membrane: the effective Poisson's ratio $\barν$ and the membrane spring constant $Λ$, the latter being proportional to the crust thickness and effective shear modulus. I solve membrane theory in terms of tidal Love numbers, for which I derive analytical formulas depending on $Λ$, $\barν$, the ocean-to-bulk density ratio and the number $k_2^o$ representing the influence of the deep interior. Membrane formulas predict $h_2$ and $k_2$ with an accuracy of a few tenths of percent if the crust thickness is less than one hundred kilometers, whereas the error on $l_2$ is a few percents. Benchmarking with the thick-shell software SatStress leads to the discovery of an error in the original, uncorrected version of the code that changes stress components by up to 40%. Regarding tectonics, I show that different stress-free states account for the conflicting predictions of thin and thick shell models about the magnitude of tensile stresses due to nonsynchronous rotation. Regarding dissipation, I prove that tidal heating in the crust is proportional to $Im Λ$ and that it is equal to the global heat flow (proportional to $Im k_2$) minus the core-mantle heat flow (proportional to $Im k_2^o$). As an illustration, I compute the equilibrium thickness of a convecting crust. More generally, membrane formulas are useful in any application involving tidal Love numbers such as crust thickness estimates, despinning tectonics or true polar wander.

preprint2012arXiv

Spatial patterns of tidal heating

In a body periodically strained by tides, heating produced by viscous friction is far from homogeneous. I show here that the distribution of the dissipated power within a spherically stratified body is a linear combination of three angular functions. These angular functions depend only on the tidal potential whereas the radial weights are specified by the internal structure of the body. The 3D problem of predicting spatial patterns of dissipation at all radii is thus reduced to the 1D problem of computing weight functions. I compute spatial patterns in various toy models without assuming a specific rheology: a viscoelastic thin shell stratified in conductive and convective layers, an incompressible homogeneous body and a two-layer model of uniform density with a liquid or rigid core. For a body in synchronous rotation undergoing eccentricity tides, dissipation in a mantle surrounding a liquid core is highest at the poles. Within a softer layer (asthenosphere or icy layer), the same tides generate maximum heating in the equatorial region with a significant degree-four structure if the layer is thin. Tidal heating patterns are thus of three main types: mantle dissipation (including the case of a floating icy crust), dissipation in a thin soft layer and dissipation in a thick soft layer. I illustrate the method with applications to Europa, Titan and Io. The formalism described in this paper applies to dissipation within solid layers of planets and satellites for which internal spherical symmetry and viscoelastic linear rheology are good approximations.

preprint2010arXiv

East-west faults due to planetary contraction

Contraction, expansion and despinning have been common in the past evolution of Solar System bodies. These processes deform the lithosphere until it breaks along faults. The type and orientation of faults are usually determined under the assumption of a constant lithospheric thickness, but lithospheric thinning can occur at the equator or at the poles due either to latitudinal variation in solar insolation or to localized tidal dissipation. Using thin elastic shells with variable thickness, I show that the equatorial thinning of the lithosphere transforms the homogeneous and isotropic fault pattern caused by contraction/expansion into a pattern of faults striking east-west, preferably formed in the equatorial region. By contrast, lithospheric thickness variations only weakly affect the despinning faulting pattern consisting of equatorial strike-slip faults and polar normal faults. If contraction is added to despinning, the despinning pattern first shifts to thrust faults striking north-south and then to thrust faults striking east-west. If the lithosphere is thinner at the poles, the tectonic pattern caused by contraction/expansion consists of faults striking north/south. On Iapetus, contraction or expansion on a lithosphere thinner at the equator explains the location and orientation of the equatorial ridge. On Mercury, the combination of contraction and despinning makes possible the existence of zonal provinces of thrust faults differing in orientation (north-south or east-west), which may be relevant to the orientation of lobate scarps.