From the mechanical point of view the interior of the Earth can be considered as composed of four main layers: the inner and outer core, the mantle and the lithosphere. The lithosphere can be assumed to be elastic and the solid mantle beneath behaves as a viscous fluid. The long term equilibrium pressure at a given depth in the Earth is due to the weight of the material above this depth. Deviations from this equilibrium state lead to material transport from regions of higher pressure towards lower pressure. If left undisturbed over time the mantle and the lithosphere reach an equilibrium, in which the depth of the base of the lithosphere will mainly depend on the thickness of the lithosphere. The growth of ice sheets during a glacial period concentrates mass on the Earth's surface to glaciated areas; this fact increases the pressure in the layers below, resulting in a sinking of the lithosphere and in a transport of mantle material away from the region. At the end of the glacial period, when the ice sheets melt away, the pressure on the lithosphere is reduced and the material will flow back causing the surface to uplift. In this talk I will present a discontinuous Galerkin finite element parallel approximation for forward modelling of the viscoelastic response of a three dimensional elastically compressible Earth to an arbitrary surface load. The code is able to perform global simulation of the rebound process, with more refined results on a selected geographical region.