Simulations of dated radiostratigraphy for the Greenland ice sheet

Author(s):  
Andreas Born ◽  
Alexander Robinson

<p>As layers of accumulated snow compact into ice and start to flow under its own weight, their deformations are recorded in the vertical structure of the glacier. Therefore, the isochronal stratigraphy of the Greenland ice sheet provides comprehensive dynamic constraints, which, with adequate methods, can be used to calibrate ice sheet models and greatly improve their accuracy.<br><br>We present the first three-dimensional ice sheet model that explicitly resolves isochrones. Individual layers of accumulation do not exchange mass with each other as the flow of ice deforms them, resembling the Lagrangian description of flow in the vertical dimension, while lateral flow within each layer is Eulerian. Direct comparison with dated radiostratigraphy is used to filter an ensemble of simulations of the Greenland ice sheet. The abundant information implied by the shape of the three-dimensional layering enables us to constrain a large number of degrees of freedom. The mismatch in the thickness of certain isochrones is used to calibrate the climate forcing of different periods of the last glacial cycle.</p>

1995 ◽  
Vol 21 ◽  
pp. 103-110 ◽  
Author(s):  
G. S. Boulton ◽  
N. Hulton ◽  
M. Vautravers

A numerical model is used to simulate ice-sheet behaviour in Europe through the last glacial cycle. It is used in two modes: a forward mode, in which the model is driven by a proxy palaeoclimate record and the output compared with a geological reconstruction of ice-sheet fluctuation; and an inverse mode, in which we determine the climate function that would be required to simulate geologically reconstructed ice-sheet fluctuations. From these simulations it is concluded that extra-glacial climates may be poor predictors of ice-sheet surface climates, and that climatic transitions during the glacial period may have been much more rapid and the intensity of warming during the early Holocene much greater than hitherto supposed. Stronger climate forcing is required to drive ice-sheet expansion when sliding occurs at the bed compared with a non-sliding bed. Sliding ice sheets grow more slowly and decay more rapidly than non-sliding ice sheets with the same climate forcing.


2019 ◽  
Vol 65 (252) ◽  
pp. 645-661 ◽  
Author(s):  
LU NIU ◽  
GERRIT LOHMANN ◽  
SEBASTIAN HINCK ◽  
EVAN J. GOWAN ◽  
UTA KREBS-KANZOW

ABSTRACTThe evolution of Northern Hemisphere ice sheets through the last glacial cycle is simulated with the glacial index method by using the climate forcing from one General Circulation Model, COSMOS. By comparing the simulated results to geological reconstructions, we first show that the modelled climate is capable of capturing the main features of the ice-sheet evolution. However, large deviations exist, likely due to the absence of nonlinear interactions between ice sheet and other climate components. The model uncertainties of the climate forcing are examined using the output from nine climate models from the Paleoclimate Modelling Intercomparison Project Phase III. The results show a large variability in simulated ice sheets between the different models. We find that the ice-sheet extent pattern resembles summer surface air temperature pattern at the Last Glacial Maximum, confirming the dominant role of surface ablation process for high-latitude Northern Hemisphere ice sheets. This study shows the importance of the upper boundary condition for ice-sheet modelling, and implies that careful constraints on climate output is essential for simulating realistic glacial Northern Hemisphere ice sheets.


1996 ◽  
Vol 23 ◽  
pp. 226-236 ◽  
Author(s):  
Philippe Huybrechts

A high-resolution, three-dimensional thermomechanical ice-sheet model, which includes isostasy, the possibility of ice-sheet expansion on the continental shelf and refined climatic parameterizations, was used to investigate the basal thermal regime of the Greenland ice sheet. The thermodynamic calculations take into account the usual terms of heat flow within the ice, a thermally active bedrock layer and all of the effects associated with changes in ice thickness and flow pattern. Basal temperature conditions are documented with respect to glacial–interracial shifts in climatic boundary conditions, both in steady state as during simulations over the last two glacial cycles using the GRIP δ180 record. It is found that the basal temperature field shows a large sensitivity in steady-state experiments but that, during a glacial cycle, basal temperature variations are strongly damped, in particular in central areas. A comparison has been made with measured data from deep ice cores and the implications are discussed.


1997 ◽  
Vol 25 ◽  
pp. 259-268 ◽  
Author(s):  
Mikhail Verbitsky ◽  
Barry Saltzman

A three-dimensional (3-D), high-resolution, non-linearly viscous, non-isothermal ice-sheet model is employed to calculate the “present-day” equilibrium regime of the Antarctic ice sheet and its evolution during the last glacial cycle. The model is augmented by an approximate formula for ice-sheet basal temperature, based on a scaling of the thermodynamic equation for the ice flow. Steady-state solutions for both the shape and extent of the areas of basal melting (or freezing) are shown to be in good agreement with those obtained from the solution of the full 3-D thermodynamic equation. The solution for the basal temperature field of the West Antaretie Siple Coast produces areas at the pressure-melting point separated by strips of frozen-to-bed ice, the structure of which is reminiscent of Ice Streams A–E. This configuration appears to be robust, preserving its features in spite of climatic changes during the last glacial cycle. Ice Stream C seems to be more vulnerable to stagnation, switching to a passive mode at least once during the penultimate interglacial. We conjecture that the peculiarities of local topography determine the unique behavior of Ice Stream C: reduced basal stress and, consequently, relatively weak warming due to internal friction and basal sliding is not able to counteract the advective cooling during the periods of increased snowfall rate.


1995 ◽  
Vol 21 ◽  
pp. 103-110
Author(s):  
G. S. Boulton ◽  
N. Hulton ◽  
M. Vautravers

A numerical model is used to simulate ice-sheet behaviour in Europe through the last glacial cycle. It is used in two modes: a forward mode, in which the model is driven by a proxy palaeoclimate record and the output compared with a geological reconstruction of ice-sheet fluctuation; and an inverse mode, in which we determine the climate function that would be required to simulate geologically reconstructed ice-sheet fluctuations.From these simulations it is concluded that extra-glacial climates may be poor predictors of ice-sheet surface climates, and that climatic transitions during the glacial period may have been much more rapid and the intensity of warming during the early Holocene much greater than hitherto supposed. Stronger climate forcing is required to drive ice-sheet expansion when sliding occurs at the bed compared with a non-sliding bed. Sliding ice sheets grow more slowly and decay more rapidly than non-sliding ice sheets with the same climate forcing.


1996 ◽  
Vol 23 ◽  
pp. 226-236 ◽  
Author(s):  
Philippe Huybrechts

A high-resolution, three-dimensional thermomechanical ice-sheet model, which includes isostasy, the possibility of ice-sheet expansion on the continental shelf and refined climatic parameterizations, was used to investigate the basal thermal regime of the Greenland ice sheet. The thermodynamic calculations take into account the usual terms of heat flow within the ice, a thermally active bedrock layer and all of the effects associated with changes in ice thickness and flow pattern. Basal temperature conditions are documented with respect to glacial–interracial shifts in climatic boundary conditions, both in steady state as during simulations over the last two glacial cycles using the GRIP δ 180 record. It is found that the basal temperature field shows a large sensitivity in steady-state experiments but that, during a glacial cycle, basal temperature variations are strongly damped, in particular in central areas. A comparison has been made with measured data from deep ice cores and the implications are discussed.


The Holocene ◽  
2018 ◽  
Vol 28 (9) ◽  
pp. 1535-1544 ◽  
Author(s):  
Laurence M Dyke ◽  
Anna LC Hughes ◽  
Camilla S Andresen ◽  
Tavi Murray ◽  
John F Hiemstra ◽  
...  

Large marine-terminating glaciers around the margins of the Greenland Ice Sheet have retreated, accelerated and thinned over the last two decades. Relatively little is known about the longer term behaviour of the Greenland Ice Sheet, yet this information is valuable for assessing the significance of modern changes. We address this by reporting 11 new beryllium-10 (10Be) exposure ages from previously uninvestigated coastal areas across southeast Greenland. The new ages are combined with existing data from the region to assess the timing of glacier retreat after the Last Glacial Maximum. The results show that deglaciation occurred first in the north of the region (~68°N) and progressed southwards. This north–south progression is attributed to the influence of the warm Irminger Current on the ice margin. Areas in the south of the region were isolated from the warm waters by the shallow bathymetry of the continental shelf. This demonstrates that oceanographic forcing paced the deglaciation of southeast Greenland through the Younger Dryas and early Holocene. In most areas of southeast Greenland bedrock ages are systematically older than their counterpart boulder samples; this offset is likely the result of inherited 10Be content in bedrock surfaces. This suggests that subglacial erosion during the last glacial cycle was insufficient to completely remove pre-existing 10Be content. Alternatively, this pattern may be the signature of a substantial retreat and advance cycle prior to final Holocene deglaciation.


1996 ◽  
Vol 23 ◽  
pp. 309-317 ◽  
Author(s):  
Emmanuel Le Meur ◽  
Philippe Huybrechts

The bedrock isostatic response exerts a strong control on ice-sheet dynamics and is therefore always taken into account in ice-sheet models. This paper reviews the various methods normally used in the ice-sheet modelling community to deal with the bedrock response and compares these with a more sophisticated full-Earth model. Each of these bedrock treatments, five in total, is coupled with a three-dimensional thermomechanical ice-sheet model under the same forcing conditions to simulate the Antarctic ice sheet during the last glacial cycle. The outputs of the simulations are compared on the basis of the time-dependent behaviour for the total ice volume and the mean bedrock elevation during the cycle and of the present rate of uplift over Antarctica. This comparison confirms the necessity of accounting for the elastic bending of the lithosphere in order to yield realistic bedrock patterns. It furthermore demonstrates the deficiencies inherent to the diffusion equation in modelling the complex deformation within the mantle. Nevertheless, when characteristic parameters are varied within their range of uncertainty, differences within one single method are often of the same order as those between the various methods. This overview finally attempts to point out the main advantages and drawbacks of each of these methods and to determine which one is most appropriate depending on the specific modelling requirements.


1997 ◽  
Vol 25 ◽  
pp. 259-268 ◽  
Author(s):  
Mikhail Verbitsky ◽  
Barry Saltzman

A three-dimensional (3-D), high-resolution, non-linearly viscous, non-isothermal ice-sheet model is employed to calculate the “present-day” equilibrium regime of the Antarctic ice sheet and its evolution during the last glacial cycle. The model is augmented by an approximate formula for ice-sheet basal temperature, based on a scaling of the thermodynamic equation for the ice flow. Steady-state solutions for both the shape and extent of the areas of basal melting (or freezing) are shown to be in good agreement with those obtained from the solution of the full 3-D thermodynamic equation. The solution for the basal temperature field of the West Antarctic Siple Coast produces areas at the pressure-melting point separated by strips of frozen-to-bed ice, the structure of which is reminiscent of Ice Streams A–E. This configuration appears to be robust, preserving its features in spite of climatic changes during the last glacial cycle. Ice Stream C seems to be more vulnerable to stagnation, switching to a passive mode at least once during the penultimate interglacial. We conjecture that the peculiarities of local topography determine the unique behavior of Ice Stream C: reduced basal stress and, consequently, relatively weak warming due to internal friction and basal sliding is not able to counteract the advective cooling during the periods of increased snowfall rate.


2003 ◽  
Vol 36 ◽  
pp. 117-128 ◽  
Author(s):  
Geoffrey S. Boulton ◽  
Magnus Hagdorn ◽  
Nicholas R.J. Hulton

AbstractGeological evidence indicates that the flow of the last European ice sheet was dominated by numerous large ice streams. Although some were ephemeral, most were sustained along well-defined axes at least during the period of retreat after the Last Glacial Maximum. A thermomechanically coupled three-dimensional numerical ice-sheet model has been used to simulate the ice sheet through the whole of the last glacial cycle, but with a spatial resolution that is high enough to capture streaming behaviour. An experiment with a smoothed bed is used to explore the self-organizing behaviour of streams when they are not forced by bed topography. On such a bed, streams typically have a width of 1–10 km, much narrower than the inferred European ice streams. An experiment using a realistic topography suggests that widths of ice streams are strongly influenced by topography, and tend to be of order 100 km. Moreover, even where the topography is muted, it stabilizes the locations of ice streams which, once formed, tend to be sustained along pre-existing axes. The model creates patterns of streaming that are similar to inferred patterns, suggesting strong topographic forcing. In a simulation using a realistic bed in which the ice was very cold and basal melting rarely occurred, streams were again very narrow. Widespread streaming under low driving stresses tends to reduce ice-sheet thicknesses compared with weak streaming or models that do not produce streaming. Consequently, ice thicknesses are smaller and tend to be consistent with the results of sea-level inversions based on geophysical Earth models.


Sign in / Sign up

Export Citation Format

Share Document