VILMA (Viscoelastic Lithosphere and Mantle model)#
Contact person
Description#
The topic of spectral finite-element code VILMA is the reduction of global GRACE, GPS and altimetry data with respect to the glacial-isostatic adjustment applying a 3D viscoelastic earth model. The model calculates the deformation of a viscoelastic and gravitating continuum in spherical domain, where lateral viscosity variations can be considered. Loading is prescribed as ice and ocean mass changes, which are determined consistently with respect to mass conservation, geoid changes and shoreline displacements by the sea-level equation. Code solves field equations of a spherical self-gravitating incompressible 3D-viscoelastic sphere with spectral–finite elements (Martinec 2000). GRD, that is rotational feedback and sea-level equation are solved for (Klemann et al. 2024 in prep.). Output are time-dependent changes in sea level, surface deformation and gravity potential in response to surface mass distribution. It is already coupled offline to PISM (Albrecht et al. 2024 in revision), MPI-ESM and AWI ESM and as compiled module in CLIMBER-X (Willeit et al. 2022). Offline coupling is implemented as exchange of two fields (in: ice load distribution; out: relative sea level).
Prototypic derivates: earthquakes forcing, compressibility, Burgers rheology, 3D elasticity, 3D density and non-linear rheology (Tanaka et al., 2009, 2011, 2015, 2021, 2023, in prep., respectively), data assimilation PDAF (Schachtschneider et al., 2021).
Model(s)#
VILMA :OS in preparation
Programming language(s)#
FORTRAN