Constraints on Glacial Isostatic Adjustment from GOCE and Sea Level Data

2009 ◽  
Vol 166 (8-9) ◽  
pp. 1261-1281 ◽  
Author(s):  
L. L. A. Vermeersen ◽  
H. H. A. Schotman
Keyword(s):  
Sea Level ◽  
Glacial Isostatic Adjustment ◽  
Isostatic Adjustment ◽  
Level Data ◽  
And Sea Level

Observations of postglacial uplift at Churchill, Manitoba

1992 ◽  
Vol 29 (11) ◽  
pp. 2418-2425 ◽  
Author(s):  
A. Mark Tushingham
Keyword(s):  
Sea Level ◽  
Tide Gauge ◽  
Ice Sheet ◽  
Laurentide Ice Sheet ◽  
Glacial Isostatic Adjustment ◽  
Tide Gauge Record ◽  
Isostatic Adjustment ◽  
Absolute Gravimetry ◽  
The Earth ◽  
Level Data

Churchill, Manitoba, is located near the centre of postglacial uplift caused by the Earth's recovery from the melting of the Laurentide Ice Sheet. The value of present-day uplift at Churchill has important implications in the study of postglacial uplift in that it can aid in constraining the thickness of the ice sheet and the rheology of the Earth. The tide-gauge record at Churchill since 1940 is examined, along with nearby Holocene relative sea-level data, geodetic measurements, and recent absolute gravimetry measurements, and a present-day rate of uplift of 8–9 mm/a is estimated. Glacial isostatic adjustment models yield similar estimates for the rate of uplift at Churchill. The effects of the tide-gauge record of the diversion of the Churchill River during the mid-1970's are discussed.

scholarly journals Relative sea-level data from the SEAMIS database compared to ICE-5G model predictions of glacial isostatic adjustment

Data in Brief ◽  
2019 ◽  
Vol 27 ◽  
pp. 104600 ◽  
Author(s):  
Thomas Mann ◽  
Maren Bender ◽  
Thomas Lorscheid ◽  
Paolo Stocchi ◽  
Matteo Vacchi ◽  
...  
Keyword(s):  
Sea Level ◽  
Glacial Isostatic Adjustment ◽  
Relative Sea Level ◽  
Isostatic Adjustment ◽  
Level Data ◽  
Model Predictions

Inferring mantle viscosity through data assimilation of relative sea level observations in a glacial isostatic adjustment model

2021 ◽  
Author(s):  
Reyko Schachtschneider ◽  
Jan Saynisch-Wagner ◽  
Volker Klemann ◽  
Meike Bagge ◽  
Maik Thomas
Keyword(s):  
Data Assimilation ◽  
Sea Level ◽  
Glacial Isostatic Adjustment ◽  
Layer Model ◽  
Relative Sea Level ◽  
Isostatic Adjustment ◽  
Mantle Viscosity ◽  
Level Data ◽  
The Time Domain ◽  
Spectral Finite Element

<p>We suggest to apply data assimilation in glacial isostatic adjustment (GIA) to constrain the mantle viscosity structure based on sea level observations. We apply the Parallel Data Assimilation Framework (PDAF) to assimilate sea level data into the time-domain spectral-finite element code VILMA in order to obtain better estimates of the mantle viscosity structure. In a first step, we reduce to a spherically symmetric earth structure and prescribe the glaciation history. A particle filter is used to propagate an ensemble of models in time. At epochs when observations are available, each particle's performance is estimated and the particles are resampled based on their performance to form a new ensemble that better resembles the true viscosity distribution.</p><p>Using this algorithm, we show the ability to recover mantle viscosities from a set of synthetic relative sea level observations. Those synthetic observations are obtained from a reference run with a given viscosity structure that defines the target viscosity values in our experiments. The viscosity estimation is applied to a three-layer model with an elastic lithosphere and two mantle layers, and to a multi-layer model with a smoother viscosity profile. We use various subsets of realistic observation locations (e.g. only observations from Fennoscandia) and show that it is possible to obtain the target viscosity values in those cases. We also vary the time from which observations are available to evolve the test cases towards a realistic scenario for the availability of relative sea level observations. The most relevant cases start at 26.5ka BP and at 10ka BP as they mark the beginning of the maximum glaciation and the end of deglaciation with a larger amount of observations following, respectively, and end at present day.</p>

scholarly journals An approach for constraining mantle viscosities through assimilation of paleo sea level data into a glacial isostatic adjustment model

2021 ◽  
Author(s):  
Reyko Schachtschneider ◽  
Jan Saynisch-Wagner ◽  
Volker Klemann ◽  
Meike Bagge ◽  
Maik Thomas
Keyword(s):  
Probability Distribution ◽  
Sea Level ◽  
Elastic Model ◽  
Glacial Isostatic Adjustment ◽  
Sea Level Changes ◽  
Isostatic Adjustment ◽  
Adjustment Model ◽  
Level Data ◽  
Special Cases ◽  
Initial Ensemble

Abstract. Glacial isostatic adjustment is largely governed by rheological properties of the Earth's mantle. Large mass redistributions in the ocean-cryosphere system and the subsequent response of the visco-elastic Earth have led to dramatic sea level changes in the past. This process is ongoing and in order to understand and predict current and future sea level changes the knowledge of mantle properties such as viscosity is essential. In this study we present a method to obtain estimates of mantle viscosities by assimilation of relative sea level data into a visco-elastic model of the lithosphere and mantle. We set up a particle filter with probabilistic resampling. In an identical twin experiment we show that mantle viscosities can be recovered in a glacial isostatic adjustment model of a simple three layer earth structure consisting of an elastic lithosphere and two mantle layers of different viscosity. In two scenarios we investigate the dependence of the ensemblebehavior on the ensemble initialization and observation uncertainties and show that the recovery is successful if the target parameter values are properly sampled by the initial ensemble probability distribution. This even includes cases in which the target viscosity values are located far in the tail of the initial ensemble probability distribution. We then successfully apply the method to two special cases that are relevant for the assimilation of real observations: 1) using observations taken from a single region only, here Laurentide and Fennoscandia, respectively, and 2) using only observations from the last 10 kyrs.

scholarly journals Comparing a thermo-mechanical Weichselian ice sheet reconstruction to GIA driven reconstructions: aspects of earth response and ice configuration

2013 ◽  
Vol 5 (2) ◽  
pp. 2345-2388 ◽  
Author(s):  
P. Schmidt ◽  
B. Lund ◽  
J-O. Näslund
Keyword(s):  
Sea Level ◽  
Ice Sheet ◽  
Younger Dryas ◽  
Slow Phase ◽  
Glacial Isostatic Adjustment ◽  
Relative Sea Level ◽  
Isostatic Adjustment ◽  
Earth Models ◽  
Uplift Rates ◽  
Best Fit

Abstract. In this study we compare a recent reconstruction of the Weichselian ice-sheet as simulated by the University of Main ice-sheet model (UMISM) to two reconstructions commonly used in glacial isostatic adjustment (GIA) modeling: ICE-5G and ANU (also known as RSES). The UMISM reconstruction is carried out on a regional scale based on thermo-mechanical modelling whereas ANU and ICE-5G are global models based on the sea-level equation. The Weichselian ice-sheet in the three models are compared directly in terms of ice volume, extent and thickness, as well as in terms of predicted glacial isostatic adjustment in Fennoscandia. The three reconstructions display significant differences. UMISM and ANU includes phases of pronounced advance and retreat prior to the last glacial maximum (LGM), whereas the thickness and areal extent of the ICE-5G ice-sheet is more or less constant up until LGM. The final retreat of the ice-sheet initiates at earliest time in ICE-5G and latest in UMISM, while ice free conditions are reached earliest in UMISM and latest in ICE-5G. The post-LGM deglaciation style also differs notably between the ice models. While the UMISM simulation includes two temporary halts in the deglaciation, the later during the Younger Dryas, ANU only includes a decreased deglaciation rate during Younger Dryas and ICE-5G retreats at a relatively constant pace after an initial slow phase. Moreover, ANU and ICE-5G melt relatively uniformly over the entire ice-sheet in contrast to UMISM which melts preferentially from the edges. We find that all three reconstructions fit the present day uplift rates over Fennoscandia and the observed relative sea-level curve along the Ångerman river equally well, albeit with different optimal earth model parameters. Given identical earth models, ICE-5G predicts the fastest present day uplift rates and ANU the slowest, ANU also prefers the thinnest lithosphere. Moreover, only for ANU can a unique best fit model be determined. For UMISM and ICE-5G there is a range of earth models that can reproduce the present day uplift rates equally well. This is understood from the higher present day uplift rates predicted by ICE-5G and UMISM, which results in a bifurcation in the best fit mantle viscosity. Comparison of the uplift histories predicted by the ice-sheets indicate that inclusion of relative sea-level data in the data fit can reduce the observed ambiguity. We study the areal distributions of present day residual surface velocities in Fennoscandia and show that all three reconstructions generally over-predict velocities in southwestern Fennoscandia and that there are large differences in the fit to the observational data in Finland and northernmost Sweden and Norway. These difference may provide input to further enhancements of the ice-sheet reconstructions.

scholarly journals Glacial Isostatic Adjustment modelling: historical perspectives, recent advances, and future directions

2018 ◽  
Author(s):  
Pippa L. Whitehouse
Keyword(s):  
Solid Earth ◽  
Sea Level ◽  
Measurement Techniques ◽  
Ice Sheets ◽  
Order Effects ◽  
System Modelling ◽  
Glacial Isostatic Adjustment ◽  
Future Directions ◽  
Historical Perspectives ◽  
Isostatic Adjustment

Abstract. Glacial Isostatic Adjustment (GIA) describes the response of the solid Earth, the gravitational field, and consequently the oceans to the growth and decay of the global ice sheets. It is a process that takes place relatively rapidly, triggering 100 m-scale changes in sea level and solid Earth deformation over just a few tens of thousands of years. Indeed, the first-order effects of GIA could already be quantified several hundred years ago without reliance on precise measurement techniques and scientists have been developing a unifying theory for the observations for over 200 years. Progress towards this goal required a number of significant breakthroughs to be made, including the recognition that ice sheets were once more extensive, the solid Earth changes shape over time, and gravity plays a central role in determining the pattern of sea-level change. This article describes in detail the historical development of the field of GIA and an overview of the processes involved. Significant recent progress has been made as concepts associated with GIA have begun to be incorporated into parallel fields of research; these advances are discussed, along with the role that GIA is likely to play in addressing outstanding research questions within the field of Earth system modelling.

scholarly journals Lithosphere and upper-mantle structure of the southern Baltic Sea estimated from modelling relative sea-level data with glacial isostatic adjustment

Solid Earth ◽  
2014 ◽  
Vol 5 (1) ◽  
pp. 447-459 ◽  
Author(s):  
H. Steffen ◽  
G. Kaufmann ◽  
R. Lampe
Keyword(s):  
Baltic Sea ◽  
Sea Level ◽  
Upper Mantle ◽  
Glacial Isostatic Adjustment ◽  
Relative Sea Level ◽  
Southern Baltic Sea ◽  
Lithospheric Thickness ◽  
Isostatic Adjustment ◽  
Mantle Viscosity ◽  
Ice Model

Abstract. During the last glacial maximum, a large ice sheet covered Scandinavia, which depressed the earth's surface by several 100 m. In northern central Europe, mass redistribution in the upper mantle led to the development of a peripheral bulge. It has been subsiding since the begin of deglaciation due to the viscoelastic behaviour of the mantle. We analyse relative sea-level (RSL) data of southern Sweden, Denmark, Germany, Poland and Lithuania to determine the lithospheric thickness and radial mantle viscosity structure for distinct regional RSL subsets. We load a 1-D Maxwell-viscoelastic earth model with a global ice-load history model of the last glaciation. We test two commonly used ice histories, RSES from the Australian National University and ICE-5G from the University of Toronto. Our results indicate that the lithospheric thickness varies, depending on the ice model used, between 60 and 160 km. The lowest values are found in the Oslo Graben area and the western German Baltic Sea coast. In between, thickness increases by at least 30 km tracing the Ringkøbing-Fyn High. In Poland and Lithuania, lithospheric thickness reaches up to 160 km. However, the latter values are not well constrained as the confidence regions are large. Upper-mantle viscosity is found to bracket [2–7] × 1020 Pa s when using ICE-5G. Employing RSES much higher values of 2 × 1021 Pa s are obtained for the southern Baltic Sea. Further investigations should evaluate whether this ice-model version and/or the RSL data need revision. We confirm that the lower-mantle viscosity in Fennoscandia can only be poorly resolved. The lithospheric structure inferred from RSES partly supports structural features of regional and global lithosphere models based on thermal or seismological data. While there is agreement in eastern Europe and southwest Sweden, the structure in an area from south of Norway to northern Germany shows large discrepancies for two of the tested lithosphere models. The lithospheric thickness as determined with ICE-5G does not agree with the lithosphere models. Hence, more investigations have to be undertaken to sufficiently determine structures such as the Ringkøbing-Fyn High as seen with seismics with the help of glacial isostatic adjustment modelling.

Sign in / Sign up

Export Citation Format

Share Document