Deglacial land emergence and lateral upper-mantle heterogeneity in the Svalbard Archipelago—II. Extended results for high-resolution load models

In Paper I (Breuer & Wolf 1995), a preliminary interpretation of the postglacial land emergence observed at a restricted set of six locations in the Svalbard Archipelago was given. The study was based on a simple model of the Barents Sea ice sheet and suggested increases in lithosphere thickness and asthenosphere viscosity with increasing distance from the continental margin.
In the present paper, the newly developed high-resolution load model. BARENTS-2, and land-uplift observations from an extended set of 25 locations are used to study further the possibility of resolving lateral heterogeneity in the upper mantle below the northern Barents Sea. A comparison of the calculated and observed uplift values shows that the lithosphere thickness is not well resolved by the observations, although values above 110 km are most common for this parameter. In contrast to this, there are indications of a lateral variation of asthenosphere viscosity. Whereas values in the range 10¹⁸-10²⁰Pas are inferred for locations close to the continental margin, 10²⁰-10²¹ Pa s are suggested further away from the margin.
A study of the sensitivity of the values found for lithosphere thickness and asthenosphere viscosity to modifications of load model BARENTS-2 shows that such modifications can be largely accommodated by appropriate changes in lithosphere thickness, whereas the suggested lateral variation of asthenosphere viscosity is essentially unaffected. An estimate of the influence of the Fennoscandian. ice sheet leads to the conclusion that its neglect results in an underestimation of the thickness of the Barents Sea ice sheet by about 10 per cent.     

Post-glacial rebound and transient lower mantle rheology

Summary. Although post-glacial rebound data have been conventionally interpreted as being governed by the steady state component of the mantle viscosity spectrum, the radial profile of this parameter, which is then inferred by fitting a model to observations, is characterized by the fact that it exhibits rather slight variation with depth. This disagrees with expectations based upon microphysical models of the solid state creep process. It also disagrees with very recent inferences of the viscosity stratification based on isostatic geoid anomalies expected on the basis of the internal lateral heterogeneity of mantle density obtained from seismic tomographic analyses. The new calculations of the signatures of post-glacial rebound reported here show that these two types of information are easily reconciled if the previously inferred value of the lower mantle viscosity is interpreted as a transient value, as originally suggested by Weertman on the basis of qualitative considerations. In these new models considered here the steady state creep resistance of the lower mantle is not constrained at all by post-glacial rebound observations. It can be fixed only by an appeal to other geophysical data. Whether such models are actually required by the data should become clear in the very near future.     

Postglacial rebound and lateral viscosity variations: a semi-analytical approach based on a spherical model with Maxwell rheology

A spectral method is employed to study the response to surface loads of a Maxwell earth including lateral viscosity variations. In particular, we focus on the effects of lithospheric cratons on the long-wavelength time-dependent displacement field for simple earth models. The viscosity contrast of the craton with respect to the surrounding mantle is kept fixed, whereas its thickness is allowed to vary. We show that the long-wavelength vertical displacement is not greatly affected by the presence of a lithospheric craton, while the tangential displacement is severely modified for the case of a homogeneous mantle. With increasing harmonic degree and thickness of the craton, the load-deformation coefficients deviate from those pertaining to a homogeneous mantle with a viscosity of 10²¹ Pa s. These deviations are particularly enhanced on timescales larger than a few hundred years. These findings indicate that the interpretation of the viscosity structure of the mantle inferred from postglacial rebound signatures based on radially stratified models is affected by the presence of lateral viscosity variations.     

Effect of viscosity structure on fault potential and stress orientations in eastern Canada

Previous investigations of the causal relationship between postglacial rebound and earthquakes in eastern Canada have focused on the mode of failure and the observed timing of the pulse of earthquake/faulting activity following deglaciation. In this study, the observational database has been extended to include observed orientations of the contemporary stress field and the rotation of stress since deglacial times. It is shown that many of these observations can be explained by a realistic ice history and a viscoelastic earth with a uniform 10²¹ Pa s mantle.
The effects of viscosity structure on the above predictions are also examined. It is shown that, since most of the above observations are found within the ice margin, they are not very sensitive to lithospheric thickness. Also, the inclusion of a 25 or 50 km ductile layer within the lithosphere will not decouple the seismogenic upper crust. High viscosity (10²² Pa s) in the lower mantle is rejected by the stress orientation and rotation observations. A low-viscosity (6 times 10²⁰Pa s) upper mantle with 1.6 times 10²¹ Pa s in the upper part of the lower mantle and 3 times 10²¹ Pa s in the lower part of the lower mantle below 1200 km depth has been found to give predictions that are in general agreement with the observations.