scholarly journals An investigation into the sensitivity of postglacial decay times to uncertainty in the adopted ice history

Author(s):  
Joseph Kuchar ◽  
Glenn Milne ◽  
Alexander Hill ◽  
Lev Tarasov ◽  
Maaria Nordman

Abstract At the centers of previously glaciated regions such as Hudson Bay in Canada and the Gulf of Bothnia in Fennoscandia, it has been observed that the sea level history follows an exponential form and that the associated decay time is relatively insensitive to uncertainty in the ice loading history. We revisit the issue of decay time sensitivity by computing relative sea level histories for Richmond Gulf and James Bay in Hudson Bay and Ångerman River in Sweden for a suite of reconstructions of the North American and Fennoscandian Ice Sheets and Earth viscosity profiles. We find that while some Earth viscosity models do indeed show insensitivity in computed decay times to the ice history, this is not true in all cases. Moreover, we find that the location of the study site relative to the geometry of the ice sheet is an important factor in determining ice sensitivity, and based on our set of ice sheet reconstructions, conclude that the location of James Bay is not well-suited to a decay time analysis. We describe novel corrections to the RSL data to remove the effects associated with the spatial distribution of sea level indicators as well as for other signals unrelated to regional ice loading (ocean loading, rotation and global mean sea-level changes) and demonstrate that they can significantly affect the inference of viscosity structure. We performed a forward modelling analysis based on a commonly adopted 2-layer, sub-lithosphere viscosity structure to determine how the solution space of viscosity models changes with the input ice history at the three study sites. While the solution spaces depend on ice history, for both Richmond Gulf and Ångerman River there are regions of parameter space where solutions are common across all or most ice histories, indicating low ice load sensitivity for these mantle viscosity parameters. For example, in Richmond Gulf, upper mantle viscosity values of (0.3–0.5)x1021 Pa s and lower mantle viscosity values of (5–50)x1021 Pa s tend to satisfy the data constraint consistently for most ice histories considered in this study. Similarly, the Ångerman River solution spaces contain a solution with an upper mantle viscosity of 0.3 × 1021 Pa s and lower mantle viscosity values of (5–50)x1021 Pa s common to 9 of the 10 ice histories considered there. However, the dependence of the viscosity solution space on ice history suggests that joint estimation of ice and Earth parameters is the optimal approach.

2020 ◽  
Vol 224 (2) ◽  
pp. 961-972
Author(s):  
A G Semple ◽  
A Lenardic

SUMMARY Previous studies have shown that a low viscosity upper mantle can impact the wavelength of mantle flow and the balance of plate driving to resisting forces. Those studies assumed that mantle viscosity is independent of mantle flow. We explore the potential that mantle flow is not only influenced by viscosity but can also feedback and alter mantle viscosity structure owing to a non-Newtonian upper-mantle rheology. Our results indicate that the average viscosity of the upper mantle, and viscosity variations within it, are affected by the depth to which a non-Newtonian rheology holds. Changes in the wavelength of mantle flow, that occur when upper-mantle viscosity drops below a critical value, alter flow velocities which, in turn, alter mantle viscosity. Those changes also affect flow profiles in the mantle and the degree to which mantle flow drives the motion of a plate analogue above it. Enhanced upper-mantle flow, due to an increasing degree of non-Newtonian behaviour, decreases the ratio of upper- to lower-mantle viscosity. Whole layer mantle convection is maintained but upper- and lower-mantle flow take on different dynamic forms: fast and concentrated upper-mantle flow; slow and diffuse lower-mantle flow. Collectively, mantle viscosity, mantle flow wavelengths, upper- to lower-mantle velocities and the degree to which the mantle can drive plate motions become connected to one another through coupled feedback loops. Under this view of mantle dynamics, depth-variable mantle viscosity is an emergent flow feature that both affects and is affected by the configuration of mantle and plate flow.


Solid Earth ◽  
2014 ◽  
Vol 5 (1) ◽  
pp. 447-459 ◽  
Author(s):  
H. Steffen ◽  
G. Kaufmann ◽  
R. Lampe

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.


1992 ◽  
Vol 19 (12) ◽  
pp. 1185-1188 ◽  
Author(s):  
J. X. Mitrovica ◽  
W. R. Peltier

2008 ◽  
Vol 45 (11) ◽  
pp. 1345-1362 ◽  
Author(s):  
Audrey Dallimore ◽  
Randolph J. Enkin ◽  
Reinhard Pienitz ◽  
John R. Southon ◽  
Judith Baker ◽  
...  

The sedimentary record in a 40.9 m giant (Calypso) piston core (MD02-2494) raised from the inner basin within Effingham Inlet, British Columbia, Canada, during the 2002 Marges Ouest Nord Américaines (MONA) campaign, spans from 14 360 14C years BP (17 300 calibrated calendar (cal.) years BP) to about nine centuries before present. The core archives changes in sedimentation and sea level immediately following deglaciation of the Late Wisconsin Fraser Glaciation, which peaked about 15 000 14C years BP. The presence of the Mazama Ash in the core anchors a detailed chronology based on 49 radiocarbon dates and seven Pleistocene paleomagnetic secular variation correlations. Diatom assemblages identify a marine–freshwater–marine transition in the basin, which occurred 11 630 14C years BP (13 500 cal. years BP). At this time, a bedrock sill, presently at 46 m depth, was briefly exposed as sea level fell and then rose again during isostatic crustal adjustments. These data constrain a new sea-level curve for the outer coast of Vancouver Island covering the past 12 000 14C years BP (14 000 cal. years BP), providing new information on the nature of deglaciation along the west coast of Canada and informing interpretations of regional paleoceanographic records and mantle viscosity models.


2013 ◽  
Vol 5 (2) ◽  
pp. 2483-2507
Author(s):  
H. Steffen ◽  
G. Kaufmann ◽  
R. Lampe

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.


ARCTIC ◽  
2016 ◽  
Vol 69 (1) ◽  
pp. 99
Author(s):  
Leonard J.S. Tsuji ◽  
Amy Daradich ◽  
Natalya Gomez ◽  
Carling Hay ◽  
Jerry X. Mitrovica

<p class="Pa5">In 1905 and 1906, the Cree of the southwestern James Bay region signed Treaty No. 9 whereby they relinquished to the Canadian government their claim to the lands south of the Albany River (the northern boundary of the province of Ontario at the time). The official text of Treaty No. 9 made no mention of land submerged below water cover, and thus the Cree did not relinquish such regions at that time. By contrast, the Cree of the northwestern James Bay and southwestern Hudson Bay region who signed the 1929–30 Adhesions to Treaty No. 9 relinquished their claims to “land covered by water” for the area bounded on the south by the northerly limit of Treaty No. 9, as this clause was specifically included in the text of the adhesion. The issue of “land covered by water” is significant because the western James Bay region has been, and will continue to be, subject to sea level changes associated with ongoing adjustments due to the last ice age and modern global warming signals. In the absence of detailed maps, we used models of these processes, constrained by available geophysical and geodetic data sets, to retrodict shoreline changes and the rate of land emergence over the last two centuries within the boundaries specified by Treaty No. 9. We also project shoreline migration to the end of the 21st century within the same region. The rate of land emergence since 1905 in the area south of the Albany River is estimated as ~3.0 km<sup>2</sup>/yr. Over the next century, land will continue to emerge in this region at a mean rate of ~1.4 km<sup>2</sup>/yr. This emergent land should be a subject of consideration within any comprehensive land claim put forward by the Cree; in this regard, it will be interesting to see how the Canadian judicial system and the Comprehensive Claims Branch handle the novel issue of emergent land.</p>


2020 ◽  
Author(s):  
Jun Yan ◽  
Maxim D. Ballmer ◽  
Paul J. Tackley

&lt;p&gt;A better understanding of the Earth&amp;#8217;s compositional structure is needed to place the geochemical record of surface rocks into the context of Earth accretion and evolution. Cosmochemical constraints imply that lower-mantle rocks may be enriched in silica relative to upper-mantle pyrolite, whereas geophysical observations support whole-mantle convection and mixing. To resolve this discrepancy, it has been suggested that subducted mid-ocean ridge basalt (MORB) segregates from subducted harzburgite to accumulate in the mantle transition zone (MTZ) and/or the lower mantle. However, the key parameters that control basalt segregation and accumulation remain poorly constrained. Here, we use global-scale 2D thermochemical convection models to investigate the influence of mantle-viscosity profile, planetary-tectonic style and bulk composition on the evolution and distribution of mantle heterogeneity. Our models robustly predict that, for all cases with Earth-like tectonics, a basalt-enriched reservoir is formed in the MTZ, and a harzburgite-enriched reservoir is sustained at 660~800 km depth, despite ongoing whole-mantle circulation. The enhancement of basalt and harzburgite in and beneath the MTZ, respectively, are laterally variable, ranging from ~30% to 50% basalt fraction, and from ~40% to 80% harzburgite enrichment relative to pyrolite. Models also predict an accumulation of basalt near the core mantle boundary (CMB) as thermochemical piles, as well as moderate enhancement of most of the lower mantle by basalt. While the accumulation of basalt in the MTZ does not strongly depend on the mantle-viscosity profile (explained by a balance between basalt delivery by plumes and removal by slabs at the given MTZ capacity), that of the lowermost mantle does: lower-mantle viscosity directly controls the efficiency of basalt segregation (and entrainment) near the CMB; upper-mantle viscosity has an indirect effect through controlling slab thickness. Finally, the composition of the bulk-silicate Earth may be shifted relative to that of upper-mantle pyrolite, if indeed significant reservoirs of basalt exist in the MTZ and lower mantle.&lt;/p&gt;


2021 ◽  
Vol 12 (1) ◽  
Author(s):  
David Wallis ◽  
Lars N. Hansen ◽  
Angus J. Wilkinson ◽  
Ricardo A. Lebensohn

AbstractChanges in stress applied to mantle rocks, such as those imposed by earthquakes, commonly induce a period of transient creep, which is often modelled based on stress transfer among slip systems due to grain interactions. However, recent experiments have demonstrated that the accumulation of stresses among dislocations is the dominant cause of strain hardening in olivine at temperatures ≤600 °C, raising the question of whether the same process contributes to transient creep at higher temperatures. Here, we demonstrate that olivine samples deformed at 25 °C or 1150–1250 °C both preserve stress heterogeneities of ~1 GPa that are imparted by dislocations and have correlation lengths of ~1 μm. The similar stress distributions formed at these different temperatures indicate that accumulation of stresses among dislocations also provides a contribution to transient creep at high temperatures. The results motivate a new generation of models that capture these intragranular processes and may refine predictions of evolving mantle viscosity over the earthquake cycle.


Sign in / Sign up

Export Citation Format

Share Document