scholarly journals Basal temperature conditions of the Greenland ice sheet during the glacial cycles

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.

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.


1989 ◽  
Vol 35 (121) ◽  
pp. 406-417 ◽  
Author(s):  
Niels Reeh

AbstractSimple analytical models are developed in order to study how up-stream variations in accumulation rate and ice thickness, and horizontal convergence/ divergence of the flow influence the age and annual layer-thickness profiles in a steady-state ice sheet. Generally, a decrease/increase of the accumulation rate and an increase/decrease of the ice thickness in the up-stream direction (i.e. opposite to the flow direction) results in older/younger ice at a given depth in the ice sheet than would result if the up-stream accumulation rate and ice thickness were constant along the flow line.Convergence/divergence of the up-stream flow will decrease/increase the effect of the accumulation-rate and ice-thickness gradients, whereas convergence/divergence has no influence at all on the age and layer-thickness profiles if the up-stream accumulation rate and ice thickness are constant along the flow line.A modified column-flow model, i.e. a model for which the strain-rate profile (or, equivalently, the horizontal velocity profile) is constant down to the depth corresponding to the Holocene/Wisconsinan transition 10 750 year BP., seems to work well for dating the ice back to 10 000–11 000 year B P. at sites in the slope regions of the Greenland ice sheet. For example, the model predicts the experimentally determined age profile at Dye 3 on the south Greenland ice sheet with a relative root-mean-square error of only 3% back to c. 10 700 year B.P. As illustrated by the Milcent location on the western slope of the central Greenland ice sheet, neglecting up-stream accumulation-rate and ice-thickness gradients, may lead to dating errors as large as 3000–000 years for c. 10 000 year old ice.However, even if these gradients are taken into account, the simple model fails to give acceptable ages for 10 000 year old ice at locations on slightly sloping ice ridges with strongly divergent flow, as for example the Camp Century location. The main reason for this failure is that the site of origin of the ice cannot be determined accurately enough by the simple models, if the flow is strongly divergent.With this exception, the simple models are well suited for dating the ice at locations where the available data or the required accuracy do not justify application of elaborate numerical models. The formulae derived for the age-depth profiles can easily be worked out on a pocket calculator, and in many cases will be a sensible alternative to using numerical flow models.


2018 ◽  
Vol 14 (4) ◽  
pp. 455-472 ◽  
Author(s):  
Ilaria Tabone ◽  
Javier Blasco ◽  
Alexander Robinson ◽  
Jorge Alvarez-Solas ◽  
Marisa Montoya

Abstract. Observations suggest that during the last decades the Greenland Ice Sheet (GrIS) has experienced a gradually accelerating mass loss, in part due to the observed speed-up of several of Greenland's marine-terminating glaciers. Recent studies directly attribute this to warming North Atlantic temperatures, which have triggered melting of the outlet glaciers of the GrIS, grounding-line retreat and enhanced ice discharge into the ocean, contributing to an acceleration of sea-level rise. Reconstructions suggest that the influence of the ocean has been of primary importance in the past as well. This was the case not only in interglacial periods, when warmer climates led to a rapid retreat of the GrIS to land above sea level, but also in glacial periods, when the GrIS expanded as far as the continental shelf break and was thus more directly exposed to oceanic changes. However, the GrIS response to palaeo-oceanic variations has yet to be investigated in detail from a mechanistic modelling perspective. In this work, the evolution of the GrIS over the past two glacial cycles is studied using a three-dimensional hybrid ice-sheet–shelf model. We assess the effect of the variation of oceanic temperatures on the GrIS evolution on glacial–interglacial timescales through changes in submarine melting. The results show a very high sensitivity of the GrIS to changing oceanic conditions. Oceanic forcing is found to be a primary driver of GrIS expansion in glacial times and of retreat in interglacial periods. If switched off, palaeo-atmospheric variations alone are not able to yield a reliable glacial configuration of the GrIS. This work therefore suggests that considering the ocean as an active forcing should become standard practice in palaeo-ice-sheet modelling.


2020 ◽  
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>


1982 ◽  
Vol 28 (100) ◽  
pp. 431-455 ◽  
Author(s):  
Niels Reeh

AbstractThe differential equation determining the elevations of a perfectly plastic three-dimensional steady-state ice sheet is set up. Analytical solutions of the equation are obtained in two simple case, (1) an ice sheet on a horizontal base with an arbitrary edge curve, and (2) an ice sheet on a plane sloping base with a rectilinear ice margin. The solutions are discussed, particularly with reference to the development of ice divides and ice streams.For arbitrary base and ice-margin geometries, solutions are obtained by means of the method of characteristics, which reduces the problem to solving simultaneously three ordinary first-order differential equations. The integration, which is performed by numerical methods, is generally commenced at the ice margin, where the necessary boundary conditions are known.The method has been applied to model the elevation contours and the flow pattern of the central Greenland ice sheet, using the bottom topography revealed by radio echo soundings and the present ice margin geometry. The result is in surprisingly good agreement with our knowledge of the ice-sheet topography and flow pattern, all significant ice divides and ice streams being reproduced. This suggests, that the method can be applied to model the shape and flow pattern of ice sheets under glacial conditions, using information about former ice-margin positions.


2021 ◽  
Author(s):  
Torsten Albrecht ◽  
Meike Bagge ◽  
Ricarda Winkelmann ◽  
Volker Klemann

<p><span>The Antarctic Ice Sheet rests on a bed that is characterized by tectonical activity and hence by a heterogeneous rheology. Spots of extremely weak lithosphere structure could have strong impacts on the Glacial Isostatic Adjustment and hence on the stability of the ice sheet, possibly also for confined glacier regions and on timescales of decades down to even years (Barletta et al., 2018).</span><span><br><br></span><span>We coupled the VIscoelastic Lithosphere and MAntle model (VILMA) to the Parallel Ice Sheet Model (PISM) </span><span>and ran simulations over the last two glacial cycles. In this framework, VILMA considers both viscoelastic deformations of the solid Earth and gravitationally consistent mass redistribution in the ocean by solving for the sea-level equation (Martinec et al., 2018)</span><span>. In turn, PISM interprets this as a vertical shift in bed topography that directly affects the stress balance within the ice sheet and hence the grounding line dynamics at the interface of ice, ocean and bedrock.</span><span><br><br></span><span>Here we present first results of the coupled Antarctic glacial-cycle simulations and investigate technical aspects, such as optimal coupling time steps, iteration schemes and convergence, for both one-dimensional and three-dimensional Earth structures. This project is part of the </span><span>German Climate Modeling Initiative, PalMod2.</span></p><p> </p><p><span>References:</span></p><p><sup><span>Barletta et al., 2018. <em>Observed rapid bedrock uplift in Amundsen Sea Embayment promotes ice-sheet stability. </em><strong>Science</strong>, <em>360</em>, pp.1335-1339. DOI: 10.1126/science.aao1447</span></sup></p><p><sup><span>Martinec et al., 2018. <em>A benchmark study of numerical implementations of the sea level equation in GIA modelling</em>. </span><span><strong>Geophysical Journal International</strong></span><span>, <em>215</em>(1), pp.389-414. DOI: 10.1093/gji/ggy280</span></sup></p><p> </p>


1997 ◽  
Vol 43 (144) ◽  
pp. 307-310 ◽  
Author(s):  
Ralf Greve

Abstract The three-dimensional ice-sheet model SICOPOLIS is used to simulate the dynamic/thermody namic behaviour of the entire Greenland ice sheet from 250 000 a BP until today. External forcing consists of a surface-temperature history constructed from δ18O data of the GRIP core, a snowfall history coupled linearly to that of the surface temperature, a piecewise linear sea-level scenario and a constant geothermal heat flux. The simulated Greenland ice sheet is investigated in the vicinity of Summit, the position where the maximum elevation is taken, and where the two drill sites GRIP and GISP2 are situated 28km apart from each other. In this region, the agreement between modelled and observed topography and ice temperature turns out to be very good. Computed age-depth profiles for GRIP and GISP2 are presented, which can he used to complete the dating of these cores in the deeper regions where annual-layer counting is not possible. However, artificial diffusion influences the computed ages in a near-basal boundary layer of approximately 15% of the ice thickness, so that the age at the bottom of the cores cannot be predicted yet.


1997 ◽  
Vol 43 (144) ◽  
pp. 307-310
Author(s):  
Ralf Greve

AbstractThe three-dimensional ice-sheet model SICOPOLIS is used to simulate the dynamic/thermody namic behaviour of the entire Greenland ice sheet from 250 000 a BP until today. External forcing consists of a surface-temperature history constructed from δ18O data of the GRIP core, a snowfall history coupled linearly to that of the surface temperature, a piecewise linear sea-level scenario and a constant geothermal heat flux. The simulated Greenland ice sheet is investigated in the vicinity of Summit, the position where the maximum elevation is taken, and where the two drill sites GRIP and GISP2 are situated 28km apart from each other. In this region, the agreement between modelled and observed topography and ice temperature turns out to be very good. Computed age-depth profiles for GRIP and GISP2 are presented, which can he used to complete the dating of these cores in the deeper regions where annual-layer counting is not possible. However, artificial diffusion influences the computed ages in a near-basal boundary layer of approximately 15% of the ice thickness, so that the age at the bottom of the cores cannot be predicted yet.


2017 ◽  
Author(s):  
Ilaria Tabone ◽  
Javier Blasco ◽  
Alexander Robinson ◽  
Jorge Alvarez-Solas ◽  
Marisa Montoya

Abstract. Observations suggest that during the last decades the Greenland Ice Sheet (GrIS) has experienced a gradually accelerating mass loss, in part due to the observed acceleration of several of Greenland’s marine-terminating glaciers. Recent studies directly attribute this to increasing North Atlantic temperatures, which have triggered melting of the GrIS outlet glaciers, grounding-line retreat and enhanced ice discharge into the ocean, contributing to an acceleration of sea level rise. Reconstructions suggest that the influence of the ocean has been of primary importance in the past as well. This was the case not only in interglacial periods, when warmer climates led to a rapid retreat of the GrIS to land above sea level, but also in glacial periods, when the GrIS expanded as far as the continental shelf break, and was thus more directly exposed to ocean changes. However, the GrIS response to paleo oceanic variations has not been investigated from a modelling perspective so far. In this work the evolution of the GrIS over the past two glacial cycles has been studied using a three-dimensional hybrid ice-sheet/ice-shelf model. We assess the effect of the variation of oceanic temperatures on the GrIS evolution on glacial-interglacial timescales through changes in submarine melting. The results show a very high sensitivity of the GrIS to the changing oceanic conditions. Oceanic forcing is found to be the dominant driver of the GrIS expansion in glacial times and retreat in interglacial periods. If switched off, paleo atmospheric variations alone are not able to yield a reliable glacial configuration of the GrIS. This work therefore suggests that considering the ocean as an active forcing should become standard in paleo ice sheet modelling.


Sign in / Sign up

Export Citation Format

Share Document