Continuous GPS measurements of postglacial adjustment in Fennoscandia: 2. Modeling results

• Published in: Journal of Geophysical Research: Solid Earth
• Main Authors: Milne, G.A., Mitrovica, J.X., Scherneck, Hans-Georg, Davis, J.L., Johansson, Jan, Koivula, H., Vermeer, M.
• Language: unknown
• Published: 2004
• Subjects: Earth and Related Environmental Sciences; galaxies: Seyfert; ISM: molecules; radio lines: galaxies; Fennoscandia; Fennoscandian
• Online Access: https://doi.org/10.1029/2003JB002619; https://research.chalmers.se/en/publication/1671

[1] Data collected under the auspices of the BIFROST GPS project yield a geographically dense suite of estimates of present-day, three-dimensional (3-D) crustal deformation rates in Fennoscandia [Johansson et al., 2002]. A preliminary forward analysis of these estimates [Milne et al., 2001] has indicated that models of ongoing glacial isostatic adjustment (GIA) in response to the final deglaciation event of the current ice age are able to provide an excellent fit to the observed 3-D velocity field. In this study we revisit our previous GIA analysis by considering a more extensive suite of forward calculations and by performing the first formal joint inversion of the BIFROST rate estimates. To establish insight into the physics of the GIA response in the region, we begin by decomposing a forward prediction into the three contributions associated with the ice, ocean, and rotational forcings. From this analysis we demonstrate that recent advances in postglacial sea level theory, in particular the inclusion of rotational effects and improvements in the treatment of the ocean load in the vicinity of an evolving continental margin, involve peak signals that are larger than the observational uncertainties in the BIFROST network. The forward analysis is completed by presenting predictions for a pair of Fennoscandian ice histories and an extensive suite of viscoelastic Earth models. The former indicates that the BIFROST data set provides a powerful discriminant of such histories. The latter yields bounds on the ( assumed constant) upper and lower mantle viscosity (nu(UM), nu(LM)); specifically, we derive a 95% confidence interval of 5 x 10(20) less than or equal to nu(UM) less than or equal to 10(21) Pa s and 5 x 10(21) less than or equal to nu(LM) less than or equal to 5 x 10(22) Pa s, with some preference for (elastic) lithospheric thickness in excess of 100 km. The main goal of the ( Bayesian) inverse analysis is to estimate the radial resolving power of the BIFROST GPS data as a function of depth in the mantle. ...