Parameters controlling mid-Holocene highstand in Glacial Isostatic Adjustment modelling

2021 ◽  
Author(s):  
Tanghua Li ◽  
Stephen Chua ◽  
Nicole Khan ◽  
Patrick Wu ◽  
Benjamin Horton
Keyword(s):  
Southeast Asia ◽  
Sea Level ◽  
Late Holocene ◽  
Ice Sheets ◽  
Thickness Variation ◽  
Glacial Isostatic Adjustment ◽  
Lithospheric Thickness ◽  
Isostatic Adjustment ◽  
Mantle Viscosity ◽  
Open Questions

<p>Holocene relative sea-level (RSL) records from regions distal from ice sheets (far-field) are commonly characterized by a mid-Holocene highstand, when RSL reached higher than present levels. The magnitude and timing of the mid-Holocene highstand varies spatially due to hydro-isostatic processes including ocean syphoning and continental levering. While there are open questions regarding the timing, magnitude and source of ice-equivalent sea level in the middle to late Holocene.</p><p>Here, we compare Glacial Isostatic Adjustment (GIA) model predictions to a standardized database of sea-level index points (SLIPs) from Southeast Asia where we have near-complete Holocene records. The database has more than 130 SLIPs that span the time period from ~9.5 ka BP to present. We investigate the sensitivity of mid-Holocene RSL predictions to GIA parameters, including the lateral lithospheric thickness variation, mantle viscosity (both 1D and 3D), and deglaciation history from different ice sheets (e.g., Laurentide, Fennoscandia, Antarctica).</p><p>We compute gravitationally self-consistent RSL histories for the GIA model with time dependent coastlines and rotational feedback using the Coupled Laplace-Finite Element Method. The preliminary results show that the timing of the highstand is mainly controlled by the deglaciation history (ice-equivalent sea level), while the magnitude is dominated by Earth parameters (e.g., lithospheric thickness, mantle viscosity). We further investigate whether there is meltwater input during middle to late Holocene and whether the RSL records from Southeast Asia can reveal the meltwater source, like Antarctica.</p>

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.

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

2013 ◽  
Vol 5 (2) ◽  
pp. 2483-2507
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, and the Earth's surface was depressed by several 100 m. Beyond the limit of this Fennoscandian ice sheet, mass redistribution in the upper mantle led to the development of peripheral bulges around the glaciated region. These once uplifted areas subside since the begin of deglaciation due to the viscoelastic behavior of the mantle. Parts of this subsiding region are located in northern central Europe in the coastal parts of Denmark, Germany and Poland. We analyze relative sea-level (RSL) data of these regions to determine the lithospheric thickness and radial mantle viscosity structure for distinct regional RSL subsets. We load a one-dimensional 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 Fyn High. In Poland, lithospheric thickness values up to 160 km are reached. However, the latter values are not well constrained due to a low number of RSL data from the Polish area. 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 yield for the southern Baltic Sea, which suggests a revision of this ice-model version. We confirm that the lower-mantle viscosity in Fennoscandia can only be poorly resolved. The lithospheric structure inferred 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 models. It thus remains challenging to sufficiently determine the Fyn High as seen with seismics with the help of glacial isostatic adjustment modelling.

scholarly journals A new Holocene sea-level record for Singapore

The Holocene ◽  
2021 ◽  
pp. 095968362110190
Author(s):  
Stephen Chua ◽  
Adam D Switzer ◽  
Tanghua Li ◽  
Huixian Chen ◽  
Margaret Christie ◽  
...  
Keyword(s):  
Sea Level ◽  
Late Holocene ◽  
Ice Sheets ◽  
Early Holocene ◽  
Glacial Isostatic Adjustment ◽  
Present Level ◽  
Isostatic Adjustment ◽  
Model Predictions ◽  
Regional Processes ◽  
Holocene Sea Level

Relative sea-level (RSL) records from far-field regions distal from ice sheets remain poorly understood, particularly in the early Holocene. Here, we extended the Holocene RSL data from Singapore by producing early Holocene sea-level index points (SLIPs) and limiting dates from a new ~40 m sediment core. We merged new and published RSL data to construct a standardized Singapore RSL database consisting of 88 SLIPs and limiting data. In the early Holocene, RSL rose rapidly from −21.0 to −0.7 m from ~9500 to 7000 cal. yrs. BP. Thereafter, the rate of RSL rise decelerated, reaching a mid-Holocene highstand of 4.0 ± 4.5 m at 5100 cal. yrs. BP, before falling to its present level. There is no evidence of any inflections in RSL when the full uncertainty of SLIPs is considered. When combined with other standardized data from the Malay-Thai Peninsula, our results also show substantial misfits between regional RSL reconstructions and glacial isostatic adjustment (GIA) model predictions in the rate of early Holocene RSL rise, the timing of the mid-Holocene highstand and the nature of late-Holocene RSL fall towards the present. It is presently unknown whether these misfits are caused by regional processes, such as subsidence of the continental shelf, or inaccurate parameters used in the GIA model.

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.

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.

Late Holocene sea-level changes and isostasy in western Denmark

2006 ◽  
Vol 66 (2) ◽  
pp. 288-302 ◽  
Author(s):  
W. Roland Gehrels ◽  
Katie Szkornik ◽  
Jesper Bartholdy ◽  
Jason R. Kirby ◽  
Sarah L. Bradley ◽  
...  
Keyword(s):  
Salt Marsh ◽  
Sea Level ◽  
Salt Marshes ◽  
Late Holocene ◽  
Glacial Isostatic Adjustment ◽  
Sea Level Changes ◽  
Relative Sea Level ◽  
Isostatic Adjustment ◽  
Tidal Embayment ◽  
Made In

AbstractCores and exposed cliff sections in salt marshes around Ho Bugt, a tidal embayment in the northernmost part of the Danish Wadden Sea, were subjected to 14C dating and litho- and biostratigraphical analyses to reconstruct paleoenvironmental changes and to establish a late Holocene relative sea-level history. Four stages in the late Holocene development of Ho Bugt can be identified: (1) groundwater-table rise and growth of basal peat (from at least 2300 BC to AD 0); (2) salt-marsh formation (0 to AD 250); (3) a freshening phase (AD 250 to AD 1600?), culminating in the drying out of the marshes and producing a distinct black horizon followed by an aeolian phase with sand deposition; and (4) renewed salt-marsh deposition (AD 1600? to present). From 16 calibrated AMS radiocarbon ages on fossil plant fragments and 4 calibrated conventional radiocarbon ages on peat, we reconstructed a local relative sea-level history that shows a steady sea-level rise of 4 m since 4000 cal yr BP. Contrary to suggestions made in the literature, the relative sea-level record of Ho Bugt does not contain a late Holocene highstand. Relative sea-level changes at Ho Bugt are controlled by glacio-isostatic subsidence and can be duplicated by a glacial isostatic adjustment model in which no water is added to the world's oceans after ca. 5000 cal yr BP.

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>

Dynamic Topography Change of the Eastern United States Since 3 Million Years Ago

Science ◽  
2013 ◽  
Vol 340 (6140) ◽  
pp. 1560-1563 ◽  
Author(s):  
David B. Rowley ◽  
Alessandro M. Forte ◽  
Robert Moucha ◽  
Jerry X. Mitrovica ◽  
Nathan A. Simmons ◽  
...  
Keyword(s):  
Sea Level ◽  
Coastal Plain ◽  
Mantle Convection ◽  
Sedimentary Rocks ◽  
Ice Sheets ◽  
Glacial Isostatic Adjustment ◽  
Eastern United States ◽  
Dynamic Topography ◽  
Isostatic Adjustment ◽  
Global Sea Level

Sedimentary rocks from Virginia through Florida record marine flooding during the mid-Pliocene. Several wave-cut scarps that at the time of deposition would have been horizontal are now draped over a warped surface with a maximum variation of 60 meters. We modeled dynamic topography by using mantle convection simulations that predict the amplitude and broad spatial distribution of this distortion. The results imply that dynamic topography and, to a lesser extent, glacial isostatic adjustment account for the current architecture of the coastal plain and proximal shelf. This confounds attempts to use regional stratigraphic relations as references for longer-term sea-level determinations. Inferences of Pliocene global sea-level heights or stability of Antarctic ice sheets therefore cannot be deciphered in the absence of an appropriate mantle dynamic reference frame.

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>

Sign in / Sign up

Export Citation Format

Share Document