scholarly journals A Multiple 1D Earth Approach (M1DEA) to account for lateral viscosity variations in solutions of the sea level equation: An application for glacial isostatic adjustment by Antarctic deglaciation

2020 ◽  
Author(s):  
Robert Hartmann ◽  
Jörg Ebbing ◽  
Clinton P. Conrad
Keyword(s):  
Sea Level ◽  
Large Scale ◽  
Earth Structure ◽  
Glacial Isostatic Adjustment ◽  
Isostatic Adjustment ◽  
Time Period ◽  
Uplift Rates ◽  
Earth Structures ◽  
Load Component ◽  
Lateral Variations

<p>The pseudo-spectral form of the sea level equation (SLE) requires the approximation of a radially-symmetric visco-elastic Earth. Thus, the resulting predictions of sea level change (SLC) and glacial isostatic adjustment (GIA) often ignore lateral variations in the Earth structure. Here, we assess the capabilities of a Multiple 1D Earth Approach (M1DEA) applied to large-scale ice load components with different Earth structures to account for these variations. In this approach the total SLC and GIA responses result from the superposition of individual responses from each load component, each computed globally assuming locally-appropriate 1D Earth structures. We apply the M1DEA to three separate regions (East Antarctica, West Antarctica, and outside Antarctica) to analyze uplift rates for a range of Earth structures and different ice loads at various distances. We find that the uplift response is mostly sensitive to the local Earth structure, which supports the usefulness of the M1DEA. However, stresses transmitted across rheological boundaries (e.g., producing peripheral bulges) present challenges for the M1DEA, but can be minimized under two conditions: (1) If the considered time period of ice loading for each component is consistent with the relaxation time of the local Earth structure. (2) If the load components can be subdivided according to the scale of the lateral variations in Earth structure. Overall, our results indicate that M1DEA could be a computationally much cheaper alternative to 3D finite element models, but further work is needed to quantify the relative accuracy of both methods for different resolutions, loads, and Earth structure variations.</p>

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.

3D glacial-isostatic adjustment models using geodynamically constrained Earth structures

2021 ◽  
Author(s):  
Meike Bagge ◽  
Volker Klemann ◽  
Bernhard Steinberger ◽  
Milena Latinović ◽  
Maik Thomas
Keyword(s):  
Solid Earth ◽  
Sea Level ◽  
Ice Sheet ◽  
Sea Level Change ◽  
Glacial Isostatic Adjustment ◽  
Level Change ◽  
Isostatic Adjustment ◽  
Mantle Viscosity ◽  
The Earth ◽  
Earth Structures

<p>The interaction between ice sheets and the solid Earth plays an important role for ice-sheet stability and sea-level change and hence for global climate models. Glacial-isostatic adjustment (GIA) models enable simulation of the solid Earth response due to variations in ice-sheet and ocean loading and prediction of the relative sea-level change. Because the viscoelastic response of the solid Earth depends on both ice-sheet distribution and the Earth’s rheology, independent constraints for the Earth structure in GIA models are beneficial. Seismic tomography models facilitate insights into the Earth’s interior, revealing lateral variability of the mantle viscosity that allows studying its relevance in GIA modeling. Especially, in regions of low mantle viscosity, the predicted surface deformations generated with such 3D GIA models differ considerably from those generated by traditional GIA models with radially symmetric structures. But also, the conversion from seismic velocity variations to viscosity is affected by a set of uncertainties. Here, we apply geodynamically constrained 3D Earth structures. We analyze the impact of conversion parameters (reduction factor in Arrhenius law and radial viscosity profile) on relative sea-level predictions. Furthermore, we focus on exemplary low-viscosity regions like the Cascadian subduction zone and southern Patagonia, which coincide with significant ice-mass changes.</p>

scholarly journals Optimal locations of sea-level indicators in glacial isostatic adjustment investigations

Solid Earth ◽  
2014 ◽  
Vol 5 (1) ◽  
pp. 511-521 ◽  
Author(s):  
H. Steffen ◽  
P. Wu ◽  
H. Wang
Keyword(s):  
Sea Level ◽  
Glacial Isostatic Adjustment ◽  
Load History ◽  
Lithospheric Thickness ◽  
Ice Load ◽  
Isostatic Adjustment ◽  
Mantle Viscosity ◽  
Almost Everywhere ◽  
Time Period ◽  
Lateral Heterogeneities

Abstract. Fréchet (sensitivity) kernels are an important tool in glacial isostatic adjustment (GIA) investigations to understand lithospheric thickness, mantle viscosity and ice-load model variations. These parameters influence the interpretation of geologic, geophysical and geodetic data, which contribute to our understanding of global change. We discuss global sensitivities of relative sea-level (RSL) data of the last 18 000 years. This also includes indicative RSL-like data (e.g., lake levels) on the continents far off the coasts. We present detailed sensitivity maps for four parameters important in GIA investigations (ice-load history, lithospheric thickness, background viscosity, lateral viscosity variations) for up to nine dedicated times. Assuming an accuracy of 2 m of RSL data of all ages (based on analysis of currently available data), we highlight areas around the world where, if the environmental conditions allowed its deposition and survival until today, RSL data of at least this accuracy may help to quantify the GIA modeling parameters above. The sensitivity to ice-load history variations is the dominating pattern covering almost the whole world before about 13 ka (calendar years before 1950). The other three parameters show distinct patterns, but are almost everywhere overlapped by the ice-load history pattern. The more recent the data are, the smaller the area of possible RSL locations that could provide enough information to a parameter. Such an area is mainly limited to the area of former glaciation, but we also note that when the accuracy of RSL data can be improved, e.g., from 2 m to 1 m, these areas become larger, allowing better inference of background viscosity and lateral heterogeneity. Although the patterns depend on the chosen models and error limit, our results are indicative enough to outline areas where one should look for helpful RSL data of a certain time period. Our results also indicate that as long as the ice-load history is not sufficiently known, the inference of lateral heterogeneities in mantle viscosity or lithospheric thickness will be interfered by the uncertainty of the ice model.

Late Quaternary sea-level history of Saipan, Commonwealth of the Northern Mariana Islands, USA: A test of tectonic uplift and glacial isostatic adjustment models

2019 ◽  
Vol 132 (3-4) ◽  
pp. 863-883
Author(s):  
Daniel R. Muhs ◽  
Eugene S. Schweig ◽  
Kathleen R. Simmons
Keyword(s):  
Sea Level ◽  
Late Quaternary ◽  
Glacial Isostatic Adjustment ◽  
Mariana Islands ◽  
Northern Mariana Islands ◽  
Isostatic Adjustment ◽  
The Pacific ◽  
Uplift Rates ◽  
Subducting Plate ◽  
End Members

Abstract In 1979, S. Uyeda and H. Kanamori proposed a tectonic model with two end members of a subduction-boundary continuum: the “Chilean” type (shallow dip of the subducting plate, great thrust events, compression, and uplift of the overriding plate) and a “Mariana” type (steep dip of the subducting plate, no great thrust events, tension, and no uplift). This concept has been used to explain variable rates of Quaternary uplift around the Pacific Rim, yet no uplift rates have been determined for the Mariana Islands themselves, one of the end members in this model. We studied the late Quaternary Tanapag Limestone, which rims much of the eastern and southern coasts of Saipan, Northern Mariana Islands, with elevations of ∼13 m to ∼30 m. Samples from 12 well-preserved corals (Acropora, Porites, and Goniastrea) yielded U-series ages ranging from ca. 134 ka to ca. 126 ka. These ages correlate the emergent reef of the Tanapag Limestone with the last interglacial period, when sea level was several meters above present. Ages and measured reef elevations from the Tanapag Limestone, along with paleo–sea-level data, yield relatively low late Quaternary uplift rates of 0.002–0.19 m/k.y., consistent with the Uyeda-Kanamori model. A review of data from other localities near subduction zones around the Pacific Basin, however, indicates that many coastlines do not fit the model. Uplift rates along the Chilean coast are predicted to be relatively high, but field studies indicate they are low. On some coastlines, relatively high uplift rates are better explained by subduction of seamounts or submarine ridges rather than subduction zone geometry. Despite the low long-term uplift rate on Saipan, the island also hosts an emergent, low-elevation (+3.9–4.0 m) reef with corals in growth position below a notch (+4.2 m). The corals are dated to 3.9–3.1 ka. The occurrence of this young, emergent reef is likely not due to tectonic uplift; instead, it is interpreted to be the result of glacial isostatic adjustment processes after the end of the last glacial period. Our findings are consistent with similar observations on tectonically stable or slowly uplifting islands elsewhere in the equatorial Pacific Ocean and agree with numerical models of a higher-than-present Holocene sea level in this region due to glacial isostatic adjustment processes.

Validating GIA models based on an ensemble of 3D Earth structures with present-day GPS uplift rates

2021 ◽  
Author(s):  
Volker Klemann ◽  
Eva Boergens ◽  
Meike Bagge
Keyword(s):  
Climate Models ◽  
Global Climate ◽  
Absolute Error ◽  
Global Climate Models ◽  
Delayed Response ◽  
Earth Structure ◽  
Isostatic Adjustment ◽  
Uplift Rates ◽  
Earth Structures ◽  
Mean Square Errors

<p>Glacial-isostatic adjustment (GIA) models simulate the viscoelastic response of the solid earth due to loading. During the last glacial maximum, large areas in the northern and southern hemisphere were covered by km-thick ice sheets. Although most of the ice has been melted already 8,000 year ago, the time-delayed response of the viscoelastic earth is still a significant contribution to present-day uplift rates. The implementation of GIA models in global climate models is an essential part of the current research. Hereby, the choice of an appropriate earth structure in the GIA model plays an important role and has to be constrained by observational data.</p><p>Here, we apply present-day uplift data to constrain a set of GIA models that differ in 3D earth structure. To this end, these different GIA models are validated against GPS uplift rates provided by Schumacher et al. (2019). The GPS stations are globally distributed and not necessarily clustered in regions with strong GIA signal. For validation, regions with the largest gradient present in the GIA signal are most crucial. Thus, we use a weighting scheme, where those GPS stations get a higher weight that are less correlated to all other stations. Additionally, uncertainties in the GPS rates appear due to the length of the GPS time series and due to station specifics such as the used GPS receiver, and are provided together with the rates as standard deviations. Thence, the weighting used for the validation is the sum of the correlation derived weights and the uncertainty derived weights.</p><p>With this weighting in place, different GIA models can be validated against present day uplift rates by means of root mean square errors or mean absolute error.</p>

scholarly journals Comparing a thermo-mechanical Weichselian Ice Sheet reconstruction to reconstructions based on the sea level equation: aspects of ice configurations and glacial isostatic adjustment

Solid Earth ◽  
2014 ◽  
Vol 5 (1) ◽  
pp. 371-388 ◽  
Author(s):  
P. Schmidt ◽  
B. Lund ◽  
J-O. Näslund ◽  
J. Fastook
Keyword(s):  
Sea Level ◽  
Regional Scale ◽  
Ice Sheet ◽  
Earth Model ◽  
Model Parameters ◽  
Glacial Isostatic Adjustment ◽  
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 Maine ice sheet model (UMISM) to two reconstructions commonly used in glacial isostatic adjustment (GIA) modelling: ICE-5G and ANU (Australian National University, 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 three models of the Weichselian Ice Sheet 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. Whereas UMISM and ANU includes phases of pronounced advance and retreat prior to the last glacial maximum (LGM), the thickness and areal extent of the ICE-5G ice sheet is more or less constant up until the LGM. During the post-LGM deglaciation phase ANU and ICE-5G melt relatively uniformly over the entire ice sheet in contrast to UMISM, which melts preferentially from the edges, thus reflecting the fundamental difference in the reconstruction scheme. We find that all three reconstructions fit the present-day uplift rates over Fennoscandia 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. 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 result in bifurcations in the best-fit upper- and lower-mantle viscosities. 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.

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 Influence of 3D Earth structure on Glacial Isostatic Adjustment in the Russian Arctic

2020 ◽  
Author(s):  
Tanghua Li ◽  
Nicole Khan ◽  
Alisa Baranskaya ◽  
Timothy Shaw ◽  
W Richard Peltier ◽  
...  
Keyword(s):  
Earth Structure ◽  
Glacial Isostatic Adjustment ◽  
Russian Arctic ◽  
Isostatic Adjustment ◽  
3D Earth Structure
Sign in / Sign up

Export Citation Format

Share Document