Hostname: page-component-848d4c4894-xfwgj Total loading time: 0 Render date: 2024-06-30T19:07:31.081Z Has data issue: false hasContentIssue false
  • English
  • Français

Basal melt at NorthGRIP modeled from borehole, ice-core and radio-echo sounder observations

Published online by Cambridge University Press:  14 September 2017

Dorthe Dahl-Jensen ,
Niels Gundestrup ,
S. Prasad Gogineni  and
Heinz Miller
Dorthe Dahl-Jensen
Affiliation:
Niels Bohr Institute of Astronomy, Physics and Geophysics, University of Copenhagen, Juliane Maries Vej 30, DK-2100 Copenhagen, E-mail: ddj@gfy.ku.dk
Niels Gundestrup
Affiliation:
Niels Bohr Institute of Astronomy, Physics and Geophysics, University of Copenhagen, Juliane Maries Vej 30, DK-2100 Copenhagen, E-mail: ddj@gfy.ku.dk
S. Prasad Gogineni
Affiliation:
Radar Systems and Remote Sensing Laboratory, University of Kansas, 2335 Irving Hill Road, Lawrence, KS 66045-7612 U.S.A.
Heinz Miller
Affiliation:
Alfred Wegener Institute for Polar and Marine Research, Columbusstrasse, D-27568 Bremerhaven, Germany
Rights & Permissions [Opens in a new window]

Abstract

From temperature measurements down through the 3001 m deep borehole at the North Greenland Icecore Project (NorthGRIP) drill site, it is now clear that the ice at the base, 3080 m below the surface, is at the pressure-melting point. This is supported by the measurements on the ice core where the annual-layer thicknesses show there is bottom melting at the site and upstream from the borehole. Surface velocity measurements, internal radio-echo layers, borehole and ice-core data are used to constrain a time-dependent flow model simulating flow along the north-northwest-trending ice-ridge flow-line, leading to the NorthGRIP site. Also time-dependent melt rates along the flowline are calculated with a heat-flow model. The results show the geothermal heat flow varies from 50 to 200 mW m–2 along the 100km section of the modeled flowline. The melt rate at the NorthGRIP site is 0.75 cm a–1, but the deep ice in the NorthGRIP core originated 50 km upstream and has experienced melt rates as high as 1.1 cm a–1.

Type
Research Article
Information
Annals of Glaciology , Volume 37 , 2003 , pp. 207 - 212
Copyright
Copyright © The Author(s) [year] 2003

Introduction

The North Greenland Icecore Project (NorthGRIP) deep drilling program is an international collaboration involving Denmark, Germany, Japan, Belgium, Sweden, Iceland, the U.S.A., France and Switzerland. The drill site is located at 75° 10’ N, 42°32’ W in North Greenland at a surface elevation of 2921m a.s.l. and an ice thickness of 3080 m. The site is located 324 km downstream along the north-northwest-trending ice ridge from the top point of the Greenland ice sheet where the European drill site GRIP is located (DahlJensen and others, 1997, 2002; Reference Johnsen, Dahl-Jensen, Dansgaard and GundestrupJohnsen and others, 1997). During the 2001 field season, a depth of 3001 m was reached and it became clear that the ice temperature was approaching the pressure-melting point and there would be basal melting at the bedrock (Reference Dahl-JensenDahl-Jensen and others, 1997; Reference GoktasGrinsted and Reference Dahl-JensenDahl-Jensen, 2002). The temperature profile alone is not sufficient to determine the basal melt rate at the NorthGRIP site. Annual-layer thicknesses from the dated NorthGRIP ice core constrain the basal vertical velocity, thereby allowing a determination of the geothermal heat flow under the NorthGRIP area. After a presentation and discussion of the temperature profile, the ice-core data are used in a two-dimensional flowline model to determine the basal melt rate and geothermal heat flow along the north-northwest-trending ice ridge leading to the NorthGRIP site.

Temperatures at NorthGRIP

At the beginning of the 2001 field season, the temperatures are measured down through the 2930 m deep liquid-filled borehole with the University of Copenhagen (UCPH) borehole logger (Reference GudmandsenGundestrup and others, 1994). To minimize disturbances of the liquid temperature due to liquid exchange and drill heat, the borehole is left undisturbed for 9 months. The temperature measurements in the borehole liquid represent the ice temperatures within 0.05 K (Reference Dahl-JensenDahl-Jensen and others, 1998). The top 100 m of the borehole is not liquid-filled, and temperature measurements here are not used because air temperatures in this part of the borehole are disturbed by the pressure and temperature conditions in the subsurface drill trench.

The NorthGRIP and GRIP temperatures are presented in Figure 1. At NorthGRIP the surface temperature is –31 °C, and the temperature at 2880 m depth is –7.8°C. The gradient of the temperature profile is –2.8 K (100 m)–1 at 2880 m depth. The dissipation energy between 2880 and 3080 m is negligible because the site is on an ice ridge with a surface velocity of only 1.33 ma–1, so the temperature gradient is assumed unchanged here. The basal temperature is thus estimated to be –2.4°C. This is within the range of possible pressure-melting temperatures (–2.0 to –2.8°C) of the basal ice under NorthGRIP (Reference Kanagaratnam, Gogineni, Gundestrup and LarsenPaterson, 1994). The NorthGRIP temperature profile (H =3080m, a =0.195 ma–1, Tsur = –31 °C) is compared to the GRIP temperature profile (H = 3 0 2 3 m , a=0.23ma–1 , Tsur =–31.7°C) (Reference Clow, Waddington and GundestrupClow and others, 1995; Reference Grinsted and Dahl-JensenGundestrup and Clow, 1997). For the upper 1000 m, the ice-temperature undulations at the two sites are similar. The temperature undulations are remnants of surface temperatures in earlier warm and cold periods during the last 5000 years (the Little Ice Age (AD 1400– 1900), the Medieval Warmth (AD 1000) and the Climatic Optimum (8–5kyr BP); Reference Dahl-JensenDahl-Jensen and others, 1998). The basal temperature at GRIP is –8.58°C and is below the pressure-melting point. The geothermal heat flow here has been modeled to be Q g e o = 51 mW m–2 .

Fig. 1. The NorthGRIP and GRIP temperature profiles as functions of depth. At NorthGRIP the temperature profile (solid black line) was measured at the beginning of the 2001 field season to 2880 m depth. The ice is 3080 m deep, and an estimate of the temperatures in the last 200 m (thin black line) shows that the ice at the base is at the pressure-melting point (–2.4°C). The GRIP temperature profile (solid gray line) was measured in 1995,3years after the drilling was successfully concluded. The temperatures in the upper 1000 m compare well, but the basal ice at GRIP is below the pressure-melting point of ice, whereas at NorthGRIP it is at pressure point.

The temperature profile at NorthGRIP does not allow a determination of the geothermal heat flow, Qgeo, because of melting of an unknown amount of basal ice. The heat flow through the ice is Qi ce = -K∂T/∂z = 55 mW m– 2 . The heat used to melt ice is the difference between the incoming heat, Qgeo, and the outgoing heat, Qice , at the base

(1)

A one-dimensional time-dependent energy-balance equation (Reference Johnsen, Dansgaard, Bard and BroeckerJohnsen and others, 1995; Reference Dahl-JensenDahl-Jensen and others, 1998) is used to model the temperature distribution through the ice as a function of the climate conditions on the surface and the geothermal heat flow from the bedrock:

(2)

where T(z, t) is the temperature as a function of z, the vertical coordinate, and t, the time. ρ is the density of ice, c(T), the specific heat capacity of ice, and K{T) the conductivity of ice. The vertical velocity, w(z, t), in the ice is given by the Dansgaard–Johnsen equations (Reference Dansgaard and JohnsenDansgaard and Johnsen, 1969; Reference Johnsen, Dansgaard, Bard and BroeckerJohnsen and others, 1995):

(3)

where z = 0 is the ice–bedrock interface and z = H is the ice surface, wbase is the (negative) melt rate, FB is the fraction of the horizontal surface velocity that is assumed to be the basal sliding velocity, h is the distance above bedrock of the kink in the Dansgaard–Johnsen model, H is the ice equivalent ice thickness and rr is defined as

(4)

where a(t) is the accumulation rate. The energy-balance equation (2) is solved for the last 100 kyr with the boundary conditions:

(5)

The first two equations describe the surface boundary conditions, the climate parameters. The time-dependent surface temperatures and accumulation rates (Fig. 2a and b) are determined for the last 100 kyr by relating their values to the dated stable-oxygen-isotope values in the NorthGRIP ice cores (Reference Hvidberg, Keller and GundestrupJohnsen and Dansgaard, 1992; Reference Dahl-Jensen, Johnsen, Hammer, Clausen, Jouzel and PeltierDahl-Jensen and others, 1993; Reference Johnsen, Dansgaard, Bard and BroeckerJohnsen and others, 1995,Reference Johnsen, Dahl-Jensen, Dansgaard and Gundestrup1999, Reference Johnsen2001). The third equation accounts for heat-flow balance at the boundary between ice and rock, z = 0. L is the specific latent heat of fusion for ice. The last equation is the boundary for the geothermal heat flow placed 3000 m down in the bedrock, i.e. at z = –3000m.

Fig. 2. A one-dimensional heat-flow model is used to model the observed temperature through the ice at NorthGRIP.The 100 kyr long climate history used as surface boundary to the heat-flow Equation (2) consists of the surface temperature history from GRIP (Reference Dahl-JensenDahl-Jensen and others, 1998) (a) and the accumulation rate from NorthGRIP (Reference JohnsenJohnsen and others, 2001) (b). (c) When the geothermal heat flow is set to 120 mW m–2, 45–55 mW m–2is conducted though the ice at the base (gray line) while 65–75 mW m–2 is used to melt ice at the bedrock (black line). The corresponding basal melt rate is shown in (d).

To investigate for possible values of the unknown geothermal heat flow, the energy-balance equation is solved for different values. The results from the calculations at NorthGRIP determine that when the heat flow is <45 mW m–2 the basal temperatures have been below the melting point (–2.4°C) for the last 100 kyr and no basal melting has taken place. When the geothermal heat flow is increased to 55 mW m–2, the basal temperature has been at the pressure-melting point for the whole 100 kyr period. Further increasing the heat flow increases the melt rate by 1 mm a–1 for each additional 9.7 mW m2. When the geothermal heat flow is > 55mW m–2 , the time-dependent temperatures in the ice are not sensitive to the value of Qgeo. The basal temperatures are always at the melting point, and the heat flow through the ice, Qice, is determined by the temperature gradient at the base. This temperature gradient changes slightly with time as ice temperatures react to time-dependent surface temperatures and accumulation rates. An example of the model results is presented in Figure 2c for the geothermal heat-flow value Qgeo = 120 mW m–2 , which is believed to be the most likely value at NorthGRIP (see further below). The basal ice is at the pressure-melting point through the whole period, and Qice varies from 45 to 55 mW m–2 . At the present time, t = 0 years, the cold temperature wave from the glacial time has reached the bedrock at NorthGRIP, and the temperature gradient therefore has its maximum value for the 100 kyr period. A maximum amount of heat flow, Qice =55mW m–2 , passes into the ice for t = 0 years. The heat available for melting Qmelt = Qgeo - Qi ce is also shown in Figure 2c, and the melt rates are shown in Figure 2d. The melt rate varies between 0.72 and 0.80 cm a–1 during the 100 kyr period, with maximum values at 60–100 kyr BP. The present values are the lowest during the 100 kyr period. The above calculations show that when the basal ice is at the pressure-melting point, the geothermal heat flow can be determined if the basal melt rate is known.

Basal Melt Rates and Geothermal Heat Flow at NorthGRIP

The ice-core data

The NorthGRIP ice core can be dated by matching dated reference horizons between the GRIP ice core and the NorthGRIP ice core. One hundred and thirty reference horizons are identified from the volcanic markers and the sharp climatic transitions found in the stable-oxygen-isotope profiles. Figure 3a presents the annual-layer thicknesses down through the GRIP and NorthGRIP ice cores based on the 130 dated reference horizons. Near the surface, the layer thickness is 0.195 m at NorthGRIP and 0.23 m at GRIP; both values correspond to the present accumulation rates at the sites. The layer thicknesses at NorthGRIP are thinner than those at GRIP in the region from the surface to 2000m depth. Below this depth, the NorthGRIP annual layers are thicker. At the NorthGRIP depth of 2895 m, where the ice is 104 kyr old, layer thicknesses at NorthGRIP are still of the order of 1cm. At GRIP the 104 kyr old ice is found at 2750 m depth, and the layer thickness here is estimated to be 0.5 cm. The 104 kyr old layers have double the thickness at NorthGRIP that they have at GRIP, and the NorthGRIP layers do not seem to be thinning to zero at the bedrock

Fig. 3. The layer thickness down through the ice at NorthGRIP (black line) and GRIP (gray line) are presented in (a). The thickness profiles are produced from 130 dated reference horizons from the two deep ice cores. In order to model the NorthGRIP, a 100 km section of an airborne RES profile from 1999 along the north-northwest-trending ice ridge leading NorthGRIP site is used (b) (Reference Chuah, Gogineni, Allen and WohletzChuah and others, 1996; Reference FujitaGogineni and others 1998, Reference Gogineni, Chuah, Allen, Jezek and Moore2001; Reference Johnsen, Clausen, Jouzel, Schwander, Sveinbjornsdottir, White, Wettlaufer, Dash and UntersteinerKanagaratnam and others, 2001 internal reflectors are related to the dated NorthGRIP ice core, and five layers are selected and used to constrain the basal using a two-dimensional model, the dated ice core and the RES profile. The NorthGRIP site is marked with a vertical white at x = 324 km. The five selected internal reflectors are shown with black curves, and their intersections with the NorthGRIP borehole are marked with black dots in (b). On the layer-thickness profiles (a), the depth of the selected reflectors is marked gray dots on the GRIP layer-thickness profile and with black dots on the NorthGRIP profile.

A one-dimensional flow model approach

A simple fit to the NorthGRIP annual layers presented in Figure 3a is attempted with a one-dimensional flow model given in Equations (3) and (4) and the climate-dependent surface accumulation rate from Figure 2b. The ice thickness is assumed to be constant with time. A best fit to the annual layers is made with h =1400m, FB = 0.135 and a mean basal melt rate given by vba –0.77 cm a–1 (Reference GoktasGrinsted and DahlJensen, 2002). This corresponds to a geothermal heat flow of 120mW m–2, the value used for the results presented in Figure 2c and d.

A two-dimensional flow model approach

A better reconstruction of the basal melt rates requires atime-dependent two-dimensional flow model because the ice found in the NorthGRIP core originates from the surface upstream from the drill site along the north-northwest-trending ice ridge. The surface velocity at NorthGRIP is 1.33 ma–1 (Hvid-berg and others, 2002). Thus the 100 kyr old ice found at 2850m depth originates from the surface in the area 50 km upstream. The estimated melt rate from the layer thickness at NorthGRIP must be an average from the 50 km region from which the ice has flowed during the last 100 kyr. An airborne radio-echo sounder (RES) profile from GRIP to NorthGRIP along the north-northwest-trending ice ridge was obtained in 1999 with an ice-penetrating radar on a NASA P3 aircraft (Reference Chuah, Gogineni, Allen and WohletzChuah and others, 1996; Dahl-Jensen and others, 1997; Reference FujitaGogineni and others, 1998, Reference Gogineni, Chuah, Allen, Jezek and Moore2001; Reference Ekholm, Forsberg and BrozenaFahnestock and others, 2001; Reference Johnsen, Clausen, Jouzel, Schwander, Sveinbjornsdottir, White, Wettlaufer, Dash and UntersteinerKanagaratnam and others, 2001). The 100 km section of the RES profile covering the flow region leading to the North-GRIP ice core is shown in Figure 3b. The vertical resolution of the image is about 4.6 m and the horizontal resolution is 160 m. Internal layers are traceable at 700 –2500m depth. The reflectors above 700 m are blocked in online data processing to obtain maximum resolution of the deeper reflectors. Below 2500m, the reflectors (if any) are too weak to be detected. The internal reflectors are layers of the same age, isochrones (Gudmandsen, 1975; Reference Gundestrup and ClowHammer, 1980). In the upper part of the ice, the reflectors are believed to coincide with volcanic horizons and density changes (Reference Fahnestock, Abdalati, Joughin, Brozena and GogineniFujita and others, 1999; Reference Gundestrup, Clausen and HansenHempel and others, 2000; Reference Johnsen, Clausen, Jouzel, Schwander, Sveinbjornsdottir, White, Wettlaufer, Dash and UntersteinerKanagaratnam and others, 2001), while deeper down they probably coincide with the abrupt climate transitions. A prominent feature is the change between closely spaced reflectors and no reflectors found at 1600m depth at NorthGRIP. This boundary marks the transition from the cold Glacial period to the warm Bølling–Allerød period (14 kyr BP) and can be traced on most of the RES profiles obtained from the Greenland ice sheet (Gudmandsen, 1975). The age of the internal reflectors can be found by comparing the RES data with the ice-core data at NorthGRIP. The five most clearly traced reflectors are marked on the RES profile (Fig.3b). The age of the reflectors ranges from 3.4 to 77.2 kyr BP. The internal layers have undulations with amplitudes of up to 200m that can be explained neither by flow over a rough bedrock nor by surface accumulation variations. The internal layers are undulating upwards at the NorthGRIP site.

With the knowledge from the NorthGRIP ice core that there is basal melt at NorthGRIP, it is concluded that there must be basal melt along the whole section of the ice-ridge flowline presented in Figure 3b. The dipping of the isochrones therefore cannot be explained by shifting melting and non-melting regions along the line. The undulations must be caused by variations of the melt rate and thus the geothermal heat flow along the line or by variations in extending flow perpendicular to the ridge. A close net of airborne RES profiles in the NorthGRIP region was made in 1996 and 1997 by the Alfred Wegener Institute research group (Reference GogineniGoktas,1999). The recorded internal layers and bedrock reflectors do not support the idea of rapid changes in the flow conditions perpendicular to the ice-ridge flowline. We will assume the undulations of the internal layers are caused by variations of the geothermal heat flow along the flowline. A simple time-dependent two-dimensional flow model is used to model the flow of the isochrones with time. Equations (2–5) are used together with the corresponding equation for the horizontal velocity:

(6)

where us ur is the surface velocity determined by:

(7)

where S(x) is the surface elevation and H(x) is the ice thickness. NorthGRIP is placed at xNorthGRIP = 324 km along the north-northwest-trending ice-ridge flowline beginning at the summit of the Greenland ice sheet. The present surface velocity at NorthGRIP, usur (xNorthGRIP; 0) = 1.33 m a (Reference HammerHvidberg and others, 2002), and the present velocities modeled by Equation (7) are in agreement with the velocities determined from the NorthGRIP strain net (Reference HammerHvidberg and others, 2002). The time dependence of the horizontal velocity is assumed to be

(8)

The ice thickness, H(x), has been determined from the surface and bedrock reflectors traced on the RES profile, and the surface gradients ∂S/∂x are determined from the surface elevations (Reference Ekholm, Forsberg and BrozenaEkholm and others, 1995; Reference EkholmEkholm, 1996). The ice thickness is assumed to be constant with time. The five traced isochrones with ages 3.4, 7.9, 14.0, 37.2 and 77.2 kyr BP are used to constrain the basal melt rates. A first model run is made assuming there is no basal melt during the 100 kyr modeled. The modeled no-melt isochrones are found at higher elevations than the corresponding observed isochrones. This is expected because the observed isochrones are dragged down by basal melt.The unknown melt rates along the flowline are estimated as average values for each 8 km section along the flowline. Perturbations from the non-melt isochrones are modeled for each 8 km section of the flowline by assuming a melt rate of 1 cm a–1 in the sections one by one. A least-squares fit of the perturbed isochrones to the observed isochrones determines the best values of the melt rates for each of the 8 km sections along the flowline. The melt rates thereby calculated are presented in Figure 4b. It is seen that the melt rate varies between 0.4 and 1.1cm a–1along the 50 km section of the flowline leading to the NorthGRIP borehole. The average value of the melt rate is 0.7 cm a–1, in agreement with the one-dimensional estimate. Using the estimated melt rate along the flowline produces the modeled isochrones, which are shown as solid white lines in Figure 4a. The flow of the ice leading to the depth of the five selected reflectors and of the 100 kyr old ice found at 2850m depth in the NorthGRIP borehole is included in Figure 4a.

Fig. 4. (a) The five isochrones (white curves) modeled using the melt rates presented in (b) are compared with the observed internal reflectors (black curves). The modeled isochrones are capable of catching the undulations of the internal layers. Flow paths for ice found at the NorthGRIP drill site are shown as solid white curves. The five upper flow paths correspond to the five selected depths shown in Figure 3, while the deepest flow path models the flow of ice that originated on the surface 50 km upstream at 100 kyr BP. This ice is found at 2850 m depth in the NorthGRIP ice core. (b) The basal melt rates along the flowline. The melt rates are seen to vary along the flowline, with minimum values at the NorthGRIP site. (c) The corresponding geothermal heat flow.

Discussion and Conclusions

The temperature calculations show that the geothermal heat flow in the vicinity of NorthGRIP must be >55mW m–2 in order to match the observed temperature profile in the borehole. The basal ice is at the pressure-melting point, and the amount of melt cannot be constrained by the energy-balance equation alone. When the heat-flow model is combined with information on annual-layer thickness from the NorthGRIP ice core, the melt rate can be constrained to be 0.75 cma–1. The heat-flow calculations show a decrease of 0.1 cma–1 during the last 100 kyr because of ice-temperature changes with varying climate. The estimated heat flow can be translated to a geothermal heat flow of 120 mW m–2.

The airborne RES profile along the north-northwest-trending ice-ridge flowline leading to the NorthGRIP drill site reveals very detailed internal reflectors that are undulating. The internal layers are modeled by assuming that the cause of the undulations is variations in geothermal heat flow and thereby melt rates along the flowline. The 100 kyr old ice that is found at 2850m depth in the NorthGRIP ice core originates 50 km upstream of the drill site. The ice found in the NorthGRIP ice core has moved over regions where the melt rate has changed from 0.4 to 1.1 cm a–1. This corresponds to geothermal heat-flow levels of 90–160 mW m–2. The melt-rate values estimated from the two-dimensional modeling have average values that match the one-dimensional modeling values.

The geothermal heat flow is very high along the whole region of the investigated flowline, and the variations in geothermal heat-flow values suggest some very prominent geological structures under the ice. The investigation of Reference Ekholm, Forsberg and BrozenaFahnestock and others (2001) over the onset area of the northwest ice stream originating 170 km south of North-GRIP close to the ice ridge reveals even higher melt rates and geothermal heat flow. Very little is known about the geological structure of the bedrock under the Greenland ice sheet. Observations of heat flow from borehole measurements in north Canada show that the heat flow can vary over small distances (Reference JohnsenMajorowicz and Embry, 1998). The internal layers from the airborne RES profiles together with ice-core data from the NorthGRIP site are used to calibrate the models. This approach offers an opportunity to map the unknown geothermal properties of the bedrock under the Greenland ice sheet with a very high resolution.

Acknowledgements

The NorthGRIP Project is directed and organized by the Department of Geophysics at the Niels Bohr Institute, University of Copenhagen. It is supported by funding agencies in Denmark, Belgium, France, Germany, Iceland, Japan, Sweden, Switzerland and the U.S.A. Radar profiles are collected as a part of the NASA Program for Arctic Regional Climate Assessment.

References

Chuah, T. S., Gogineni, S. P., Allen, C. and Wohletz, B. 1996. Radar thickness measurements over the northern part of the Greenland ice sheet. Lawrence, University of Kansas Center for Research Inc. Radar Systems and Remote Sensing Laboratory. (Technical Report 10470-3.)Google Scholar
Clow, G. D., Waddington, E. D. and Gundestrup, N. 1995. High-precision temperature measurements in the GISP2 and GRIP boreholes. Arct. Res. U.S., 9, Fall/Winter, 44–47.Google Scholar
Dahl-Jensen, D., Johnsen, S. J., Hammer, C. U., Clausen, H. B. and Jouzel, J. 1993. Past accumulation rates derived from observed annual layers in the GRIP ice core from Summit, central Greenland. In Peltier, W.R., ed. Ice in the climate system. Berlin, etc., Springer-Verlag, 517–532. (NatO ASI Series I: Global Environmental Change 12.)Google Scholar
Dahl-Jensen, D. and 9 others. 1997. A search in north Greenland for a new ice-core drill site. J. Glaciol.,43(144), 300–306.CrossRefGoogle Scholar
Dahl-Jensen, D. and 6 others. 1998. Past temperatures directly from the Greenland ice sheet. Science, 282(5387), 268–271.CrossRefGoogle ScholarPubMed
Dahl-Jensen, D. and 8 others. 2002. The NorthGRIP deep drilling programme. Ann. Glaciol., 35, 1–4.Google Scholar
Dansgaard, W. and Johnsen, S. J. 1969.A flow model and a time scale for the ice core from Camp Century, Greenland. J. Glaciol.,8(53), 215–223.CrossRefGoogle Scholar
Ekholm, S. 1996. A full coverage, high resolution topographic model of Greenland computed from a variety of digital elevation data. J. Geophys. Res.,101(B10), 21, 961–21, 972.Google Scholar
Ekholm, S., Forsberg, R. and Brozena, J.M. 1995. Accuracy of satellite altimeter elevations over the Greenland ice sheet. J. Geophys. Res.,100(C2), 2687–2696.Google Scholar
Fahnestock, M., Abdalati, W., Joughin, I., Brozena, J. and Gogineni, P. 2001. High geothermal heat flow, basal melt, and the origin of rapid ice flow in central Greenland. Science,294(5550), 2338–2342.Google Scholar
Fujita, S. and 6 others.1999. Nature of radio-echo layering in the Antarctic ice sheet detected by a two-frequency experiment. J. Geophys. Res., 104(B6), 13, 013–13, 024.Google Scholar
Gogineni, S., Chuah, T., Allen, C., Jezek, K. and Moore, R. K. 1998. An improved coherent radar depth sounder. J. Glaciol.,44(148), 659–669.CrossRefGoogle Scholar
Gogineni, S. and 9 others. 2001. Coherent radar ice thickness measurements over the Greenland ice sheet. J. Geophys. Res.,106(D24), 33, 761–33, 772.Google Scholar
Goktas, F. 1999. Ergebnisse der Untersuchung des gronlandischen Inlandeises mit dem elektromagnetischen Reflexionsverfahren in der Umgebung von GGRIP. Ber. Polarforsch./Rep. Pol. Res. 336.Google Scholar
Grinsted, A. and Dahl-Jensen, D. 2002. A Monte Carlo-tuned model of the flow in the NorthGRIP area. Ann. Glaciol., 35, 527–530.Google Scholar
Gudmandsen, P. 1975. Layerechoes in polar ice sheets. J. Glaciol.,15(73), 95–101.Google Scholar
Gundestrup, N.S. and Clow, G. D. 1997. Logging of the GRIP and GISP2 boreholes. [Abstract.] Eos,78(46), Fall Meeting Supplement, F7.Google Scholar
Gundestrup, N. S., Clausen, H. B. and Hansen, B. L. 1994. The UCPH borehole logger. Natl. Inst. Polar Res. Mem., Special Issue 49, 224–233.Google Scholar
Hammer, C. U. 1980. Acidity of polar ice cores in relation to absolute dating, past volcanism, and radio-echoes. J. Glaciol., 25(93), 359–372.CrossRefGoogle Scholar
Hempel, L., F., Thyssen, Gundestrup, N., Clausen, H. B. and Miller, H. 2000. A comparison of radio-echo sounding data and electrical conductivity of the GRIP ice core. J. Glaciol.,46(154), 369–374.CrossRefGoogle Scholar
Hvidberg, C. S., Keller, K. and Gundestrup, N. S. 2002. Mass balance and ice flow along the north-northwest ridge of the Greenland ice sheet at NorthGRIP. Ann. Glaciol., 35, 521–526.Google Scholar
Johnsen, S. J. and Dansgaard, W. 1992. On flow model dating of stable isotope records from Greenland ice cores. In Bard, E. and Broecker, W. S., eds.The last deglaciation: absolute and radiocarbon chronologies. Berlin, etc., Springer- Verlag, 13–24. (NatO ASI Series I: Global Environmental Change 2.)Google Scholar
Johnsen, S., Dahl-Jensen, D., Dansgaard, W. and Gundestrup, N. 1995. Greenland palaeo temperatures derived from GRIP bore hole temperature and ice core isotope profiles. Tellus, 47B, 624–629.Google Scholar
Johnsen, S. J. and 14 others. 1997. The δ18O record along the Greenland Ice Core Project deep ice core and the problem of possible Eemian climatic instability. J. Geophys. Res.,102(C12), 26, 397–26, 410.Google Scholar
Johnsen, S. J., Clausen, H. B., Jouzel, J., Schwander, J., Sveinbjornsdottir, A.E. and White, J. 1999. Stable isotope records from Greenland deep ice cores. In Wettlaufer, J.S., Dash, J. G. andUntersteiner, N., eds. Ice physics and the natural environment. Berlin, etc., Springer-Verlag, 89–107. (NatO ASI Series I: Global Environmental Change 56.)Google Scholar
Johnsen, S.J. and 8 others.2001. Oxygen isotope and palaeo-temperature records from six Greenland ice-corestations: Camp Century, Dye-3, GRIP, GISP2, Renland and NorthGRIP. J. Quat. Sci.,16(4), 299–307.Google Scholar
Kanagaratnam, P., Gogineni, S. P., Gundestrup, N. and Larsen, L. 2001. High-resolution radar mapping of internal layers at the North Greenland Ice Core Project. J. Geophys. Res.,106(D24), 33, 799–33, 811.Google Scholar
Majorowicz, J.A. and Embry, A. F. 1998. Present heat flow and paleogeothermalregime in the Canadian Arctic margin: analysis of industrial thermal data and coalification gradients. Tectonophysics, 291, 141–159.Google Scholar
Paterson, W. S. B. 1994. The physics of glaciers. Third edition. Oxford, etc., Elsevier.Google Scholar
Figure 0

Fig. 1. The NorthGRIP and GRIP temperature profiles as functions of depth. At NorthGRIP the temperature profile (solid black line) was measured at the beginning of the 2001 field season to 2880 m depth. The ice is 3080 m deep, and an estimate of the temperatures in the last 200 m (thin black line) shows that the ice at the base is at the pressure-melting point (–2.4°C). The GRIP temperature profile (solid gray line) was measured in 1995,3years after the drilling was successfully concluded. The temperatures in the upper 1000 m compare well, but the basal ice at GRIP is below the pressure-melting point of ice, whereas at NorthGRIP it is at pressure point.

Figure 1

Fig. 2. A one-dimensional heat-flow model is used to model the observed temperature through the ice at NorthGRIP.The 100 kyr long climate history used as surface boundary to the heat-flow Equation (2) consists of the surface temperature history from GRIP (Dahl-Jensen and others, 1998) (a) and the accumulation rate from NorthGRIP (Johnsen and others, 2001) (b). (c) When the geothermal heat flow is set to 120 mW m–2, 45–55 mW m–2is conducted though the ice at the base (gray line) while 65–75 mW m–2 is used to melt ice at the bedrock (black line). The corresponding basal melt rate is shown in (d).

Figure 2

Fig. 3. The layer thickness down through the ice at NorthGRIP (black line) and GRIP (gray line) are presented in (a). The thickness profiles are produced from 130 dated reference horizons from the two deep ice cores. In order to model the NorthGRIP, a 100 km section of an airborne RES profile from 1999 along the north-northwest-trending ice ridge leading NorthGRIP site is used (b) (Chuah and others, 1996; Gogineni and others 1998, 2001; Kanagaratnam and others, 2001 internal reflectors are related to the dated NorthGRIP ice core, and five layers are selected and used to constrain the basal using a two-dimensional model, the dated ice core and the RES profile. The NorthGRIP site is marked with a vertical white at x = 324 km. The five selected internal reflectors are shown with black curves, and their intersections with the NorthGRIP borehole are marked with black dots in (b). On the layer-thickness profiles (a), the depth of the selected reflectors is marked gray dots on the GRIP layer-thickness profile and with black dots on the NorthGRIP profile.

Figure 3

Fig. 4. (a) The five isochrones (white curves) modeled using the melt rates presented in (b) are compared with the observed internal reflectors (black curves). The modeled isochrones are capable of catching the undulations of the internal layers. Flow paths for ice found at the NorthGRIP drill site are shown as solid white curves. The five upper flow paths correspond to the five selected depths shown in Figure 3, while the deepest flow path models the flow of ice that originated on the surface 50 km upstream at 100 kyr BP. This ice is found at 2850 m depth in the NorthGRIP ice core. (b) The basal melt rates along the flowline. The melt rates are seen to vary along the flowline, with minimum values at the NorthGRIP site. (c) The corresponding geothermal heat flow.