Search

In the framework of the theory of linear viscoelasticity, we derive an analytical expression of the relaxation modulus in the Andrade model $G_{\alpha }\left( t\right) $ for the case of rational parameter \mbox{$\alpha =m/n\in (0,1)$} in terms of Mittag--Leffler functions from its Laplace transform $\tilde{G}_{\alpha }\left( s\right) $. It turns out that the expression obtained can be rewritten in terms of Rabotnov functions. Moreover, for the original parameter $\alpha =1/3$ in the Andrade model, we obtain an expression in terms of Miller-Ross functions. The asymptotic behaviours of $G_{\alpha }\left( t\right) $ for $t\rightarrow 0^{+}$ and $t\rightarrow +\infty $ are also derived applying the Tauberian theorem. The analytical results obtained have been numerically checked by solving the Volterra integral equation satisfied by $G_{\alpha }\left( t\right) $ by using a successive approximation approach, as well as computing the inverse Laplace transform of $\tilde{G}_{\alpha }\left( s\right) $ by using Talbot's method.

Context: Tidal and rotational deformation of fluid giant extra-solar planets may impact their transit light curves, making the $k_2$ Love number observable in the upcoming years. Studying the sensitivity of $k_2$ to mass concentration at depth is thus expected to provide new constraints on the internal structure of gaseous extra-solar planets. Aims: We investigate the link between the mean polar moment of inertia $N$ of a fluid, stably layered extra-solar planet and its $k_2$ Love number, aiming at obtaining analytical relationships valid, at least, for some particular ranges of the model parameters. We also seek a general, approximate relationship useful to constrain $N$ once observations of $k_2$ will become available. Methods: For two-layer fluid extra-solar planets, we explore the relationship between $N$ and $k_2$ by analytical methods, for particular values of the model parameters. We also explore approximate relationships valid over all the possible range of two-layer models. More complex planetary structures are investigated by the semi-analytical propagator technique. Results: A unique relationship between $N$ and $k_2$ cannot be established. However, our numerical experiments show that a `rule of thumb' can be inferred, valid for complex, randomly layered stable planetary structures. The rule robustly defines the upper limit to the values of $N$ for a given $k_2$, and agrees with analytical results for a polytrope of index one and with a realistic non-rotating model of the tidal equilibrium of Jupiter., Comment: Accepted for publication on Astronomy & Astrophysics

Glacial Isostatic Adjustment (GIA) caused by the melting of past ice sheets is still a major cause of sea-level variations and 3-D crustal deformation in the Mediterranean region. However, since the contribution of GIA cannot be separated from those of oceanic or tectonic origin, its role can be only assessed by numerical modelling, solving the gravitationally self-consistent Sea Level Equation. Nonetheless, uncertainties about the melting history of the late-Pleistocene ice sheets and the rheological profile of the Earth's mantle affect the GIA predictions by an unknown amount. Estimating the GIA modelling uncertainties would be particularly important in the Mediterranean region, due to the amount of high quality geodetic data from space-borne and ground-based observations currently available, whose interpretation demands a suitable isostatic correction. Here we first review previous results about the effects of GIA in the Mediterranean Sea, enlightening the variability of all the fields affected by the persistent condition of isostatic disequilibrium. Then, for the first time in this region, we adopt an ensemble modelling approach to better constrain the present-day GIA contributions to sea-level rise and geodetic variations, and their uncertainty., Comment: This is a pre-copyedited, author-produced PDF of an article accepted for publication in Geophysical Journal International following peer review. The version of record is available online at: https://doi.org/10.1093/gji/ggab508

The computation of the Love numbers for a spherically symmetric self-gravitating viscoelastic Earth is a classical problem in global geodynamics. Here we revisit the problem of the numerical evaluation of loading and tidal Love numbers in the static limit for an incompressible planetary body, adopting a Laplace inversion scheme based upon the Post-Widder formula as an alternative to the {traditional viscoelastic normal modes method. We also consider, whithin the same framework, complex-valued, frequency-dependent Love numbers that describe the response to a periodic forcing, which are paramount in the study of the tidal deformation of planets. Furthermore, we numerically obtain the time-derivatives of Love numbers, suitable for modeling geodetic signals in response to surface loads variations. A number of examples are shown, in which time and frequency-dependent Love numbers are evaluated for the Earth and planets adopting realistic rheological profiles. The numerical solution scheme is implemented in ALMA${}^3$ (the plAnetary Love nuMbers cAlculator, version 3), an upgraded open-source Fortran 90 program that computes the Love numbers for radially layered planetary bodies with a wide range of rheologies, including transient laws like Andrade or Burgers., Comment: This is a pre-copyedited, author-produced PDF of an article accepted for publication in Geophysical Journal International following peer review. The version of record is available online at https://doi.org/10.1093/gji/ggac263

We use the tidal deformations of the Moon induced by the Earth and the Sun as a tool for studying the inner structure of our satellite. Based on measurements of the degree-two tidal Love numbers k2 and h2 and dissipation coefficients from the GRAIL mission, Lunar Laser Ranging and Laser Altimetry on board of the LRO spacecraft, we perform Monte Carlo samplings for 120,000 possible combinations of thicknesses and viscosities for two classes of the lunar models. The first one includes a uniform core, a low viscosity zone (LVZ) at the core-mantle boundary, a mantle and a crust. The second one has an additional inner core. All models are consistent with the lunar total mass as well as its moment of inertia. By comparing predicted and observed parameters for the tidal deformations we find that the existence of an inner core cannot be ruled out. Furthermore, by deducing temperature profiles for the LVZ and an Earth-like mantle, we obtain stringent constraints on the radius (500 +- 1) km, viscosity,21 (4.5 +- 0.8) x 10^16 Pa.s and the density (3400 +- 10) kg/m^3 of the LVZ. We also infer the first estimation for the outer core viscosity, (2.07 +- 1.03) x 10^17 Pa.s, for two different possible structures: a Moon with a 70 km thick outer core and a large inner core (290 km radius with a density of 6000 kg/m3), and a Moon with a thicker outer core (169 km thick) but a denser and smaller inner core (219 km radius for 8000 kg/m^3).

Monitoring of vibrational eigenmodes of an elastic body excited by gravitational waves was one of the first concepts proposed for the detection of gravitational waves. At laboratory scale, these experiments became known as resonant-bar detectors first developed by Joseph Weber in the 1960s. Due to the dimensions of these bars, the targeted signal frequencies were in the kHz range. Weber also pointed out that monitoring of vibrations of Earth or Moon could reveal gravitational waves in the mHz band. His Lunar Surface Gravimeter experiment deployed on the Moon by the Apollo 17 crew had a technical failure rendering the data useless. In this article, we revisit the idea and propose a Lunar Gravitational-Wave Antenna (LGWA). We find that LGWA could become an important partner observatory for joint observations with the space-borne, laser-interferometric detector LISA, and at the same time contribute an independent science case due to LGWA's unique features. Technical challenges need to be overcome for the deployment of the experiment, and development of inertial vibration sensor technology lays out a future path for this exciting detector concept., Comment: 29 pages, 17 figures

We present a new rheological model depending on a real parameter $\nu \in [0,1]$ that reduces to the Maxwell body for $\nu=0$ and to the Becker body for $\nu=1$. The corresponding creep law is expressed in an integral form in which the exponential function of the Becker model is replaced and generalized by a Mittag-Leffler function of order $\nu$. Then, the corresponding non-dimensional creep function and its rate are studied as functions of time for different values of $\nu$ in order to visualize the transition from the classical Maxwell body to the Becker body. Based on the hereditary theory of linear viscoelasticity, we also approximate the relaxation function by solving numerically a Volterra integral equation of the second kind. In turn, the relaxation function is shown versus time for different values of $\nu$ to visualize again the transition from the classical Maxwell body to the Becker body. Furthermore, we provide a full characterization of the new model by computing, in addition to the creep and relaxation functions, the so-called specific dissipation $Q^{-1}$ as a function of frequency, which is of particularly relevance for geophysical applications, Comment: 18 pages, 8 figures. arXiv admin note: text overlap with arXiv:1701.03068

We present a new approach based on linear integro-differential operators with logarithmic kernel related to the Hadamard fractional calculus in order to generalize, by a parameter $\nu \in (0,1]$, the logarithmic creep law known in rheology as Lomnitz law (obtained for $\nu=1$). We derive the constitutive stress-strain relation of this generalized model in a form that couples memory effects and time-varying viscosity. Then, based on the hereditary theory of linear viscoelasticity, we also derive the corresponding relaxation function by solving numerically a Volterra integral equation of the second kind. So doing we provide a full characterization of the new model both in creep and in relaxation representation, where the slow varying functions of logarithmic type play a fundamental role as required in processes of ultra slow kinetics., Comment: 15 pages, 2 figures, to appear in Chaos, Solitons and Fractals (2017)

We describe the basic features of SELEN, an open source Fortran 90 program for the numerical solution of the so-called "Sea Level Equation" for a spherical, layered, non-rotating Earth with Maxwell viscoelastic rheology. The Sea Level Equation was introduced in the 70s to model the sea level variations in response to the melting of late-Pleistocene ice-sheets, but it can be also employed for predictions of geodetic quantities such as vertical and horizontal surface displacements and gravity variations on a global and a regional scale. SELEN (acronym of SEa Level EquatioN solver) is particularly oriented to scientists at their first approach to the glacial isostatic adjustment problem and, according to our experience, it can be successfully used in teaching. The current release (2.9) considerably improves the previous versions of the code in terms of computational efficiency, portability and versatility. In this paper we describe the essentials of the theory behind the Sea Level Equation, the purposes of SELEN and its implementation, and we provide practical guidelines for the use of the program. Various examples showing how SELEN can be configured to solve geodynamical problems involving past and present sea level changes and current geodetic variations are also presented and discussed.

The viscoelastic material functions for the Becker and the Lomnitz rheological models, sometimes employed to describe the transient flow of rocks, are studied and compared. Their creep functions, which are known in a closed form, share a similar time dependence and asymptotic behavior. This is also found for the relaxation functions, obtained by solving numerically a Volterra equation of the second kind. We show that the two rheologies constitute a clear example of broadly similar creep and relaxation patterns associated with neatly distinct retardation spectra, for which analytical expressions are available., Comment: 7 pages, 4 figures

In 1958 Jeffreys proposed a power law of creep, generalizing the logarithmic law earlier introduced by Lomnitz, to broaden the geophysical applications to fluid-like materials including igneous rocks. This generalized law, however, can be applied also to solid-like viscoelastic materials. We revisit the Jeffreys-Lomnitz law of creep by allowing its power law exponent $\alpha$, usually limited to the range [0,1] to all negative values. This is consistent with the linear theory of viscoelasticity because the creep function still remains a Bernstein function, that is positive with a completely monotonic derivative, with a related spectrum of retardation times. The complete range $\alpha \le 1$ yields a continuous transition from a Hooke elastic solid with no creep ($\alpha \to -\infty$) to a Maxwell fluid with linear creep ($\alpha=1$) passing through the Lomnitz viscoelastic body with logarithmic creep ($\alpha=0$), which separates solid-like from fluid-like behaviors. Furthermore, we numerically compute the relaxation modulus and provide the analytical expression of the spectrum of retardation times corresponding to the Jeffreys-Lomnitz creep law extended to all $\alpha \le 1$., Comment: 23 pages, 3 figures (5 files .ps)

The purpose of this paper is twofold: from one side we provide a general survey to the viscoelastic models constructed via fractional calculus and from the other side we intend to analyze the basic fractional models as far as their creep, relaxation and viscosity properties are considered. The basic models are those that generalize via derivatives of fractional order the classical mechanical models characterized by two, three and four parameters, that we refer to as Kelvin-Voigt, Maxwell, Zener, anti-Zener and Burgers. For each fractional model we provide plots of the creep compliance, relaxation modulus and effective viscosity in non dimensional form in terms of a suitable time scale for different values of the order of fractional derivative. We also discuss the role of the order of fractional derivative in modifying the properties of the classical models., Comment: 41 pages, 8 figures

This paper has been withdrawn, Comment: Due to the non-originality of the result

By elementary methods, I study the Love numbers of a homogeneous, incompressible, self-gravitating sphere characterized by a generalized Maxwell rheology, whose mechanical analogue is represented by a finite or infinite system of classical Maxwell elements disposed in parallel. Analytical, previously unknown forms of the complex shear modulus for the generalized Maxwell body are found by algebraic manipulation, and studied in the particular case of systems of springs and dashpots whose strength follows a power-law distribution. We show that the sphere is asymptotically stable for any choice of the mechanical parameters that define the generalized Maxwell body and analytical forms of the Love numbers are always available for generalized bodies composed by less than five classical Maxwell bodies. For the homogeneous sphere, real Laplace inversion methods based on the Post-Widder formula can be applied without performing a numerical discretization of the n-th derivative, which can be computed in a "closed-form" with the aid of the Faa di Bruno formula., Comment: 17 pages

The effect of an increase of viscosity in the earth-mantle's transition zone between 420 and 670 km on the geophysical signatures induced by postglacial rebound was analyzed using a self-gravitating model consisting of an elastic lithosphere, a three-layer viscoelastic mantle, and an inviscid core. Results show that the horizontal displacements and velocities as well as the stress pattern are extremely sensitive to the viscosity increase and to the chemical stratification of the transition zone. The findings indicate that geodetic data on intraplate deformation can be used to put constraints on the viscosity increase in the transition zone and on the amount of chemical stratification in the mantle.

A five-layer viscoelastic spherical model is used to calculate the transient displacements of postglacial rebound, the induced polar motions, and the temporal variations of the geopotential up to degree 8 of the zonal coefficients. Two models - one with two viscoelastic layers separated at 670 km, and the other with three layers in which a hard garnet layer lies between the upper and lower mantle - are compared. Forward modeling shows that it may be possible to discern the presence of a hard garnet layer with a viscosity of at least ten times greater than the upper mantle, on the basis of uplift data near the center of the former Laurentide ice-sheet and from polar wander and j2 data. Temporal variations of higher gravity harmonics, such as j6 and j8, can potentially place even tighter constraints on the rheological properties of the hard transition zone. A lower mantle viscosity between 2 and 4 x 10 to the 22nd Pa is generally preferred in models with a garnet layer which may be as large as 50 times more viscous than the upper mantle.