Researcher profile

Giorgio Spada

Giorgio Spada contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

7 published item(s)

preprint2023arXiv

New estimates of ongoing sea level change and land movements caused by Glacial Isostatic Adjustment in the Mediterranean region

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.

preprint2023arXiv

On computing viscoelastic Love numbers for general planetary models: the \texttt{ALMA${}^3$} code

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.

preprint2012arXiv

Becker and Lomnitz rheological models: a comparison

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.

preprint2012arXiv

Modeling sea level changes and geodetic variations by glacial isostasy: the improved SELEN code

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.

preprint2012arXiv

On the viscoelastic characterization of the Jeffreys-Lomnitz law of creep

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 $α$, 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 $α\le 1$ yields a continuous transition from a Hooke elastic solid with no creep ($α\to -\infty$) to a Maxwell fluid with linear creep ($α=1$) passing through the Lomnitz viscoelastic body with logarithmic creep ($α=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 $α\le 1$.

preprint2011arXiv

Creep, Relaxation and Viscosity Properties for Basic Fractional Models in Rheology

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.