Skip to main content Accessibility help
×
Home
Hostname: page-component-684bc48f8b-plp5r Total loading time: 0.728 Render date: 2021-04-12T17:25:21.547Z Has data issue: true Feature Flags: { "shouldUseShareProductTool": true, "shouldUseHypothesis": true, "isUnsiloEnabled": true, "metricsAbstractViews": false, "figures": false, "newCiteModal": false, "newCitedByModal": true }

Snow Accumulation Studies on the Thule Peninsula, Greenland

Published online by Cambridge University Press:  30 January 2017

Steven J. Mock
Affiliation:
U.S. Army Cold Regions Research and Engineering Laboratory, Hanover, New Hampshire 03755, U.S.A.

Abstract

Data from stake measurements, marker boards and pits along a 136 km trail crossing the Thule peninsula sector of the Greenland ice sheet have been used to determine both the regional and local distribution of snow accumulation, On a regional scale trend surfaces of mean annual accumulation can be adequately predicted from a model using distance from moisture source and elevation as independent parameters. A series of step- or wave-like features break the smooth profile of the ice. sheet and cause profound changes in accumulation rates on a local scale. The accumulation pattern over these features can be predicted from surface slope and departure from regional elevation. Profiles of’ surface and subsurface topography indicate a direct relationship between subsurface hills and step-like features, but cannot be quantitatively accounted for by existing ice-flow theory. Detailed accumulation studies in conjunction with a program of spirit leveling in the vicinity of Camp Century has revealed the development a shallow valley-like feature. Within this feature accumulation rates have increased indicating that it is the result of flow phenomena.

Résumé

Résumé

Éludes de l’accumulation sur la péninsule de Thulé, Groenland. Des données obtenues par observation de jalons, balises repères et puits placés le long d’une piste de 136 km traversant la région de la péninsule de Thulé sur l’indlandsis du Groenland ont été utilisées pour déterminer les distributions locales et régionales de l’accumulation. A l’échelle régionale, les tendances superficielles concernant l’accumulation annuelle moyenne peuvent être prévues de manière satisfaisante en utilisant un modéle faisant usage de la distance à la source d’humidité et de l’élévation en tant que paramètres indépendants. Une série d’irrégularités topographiques en forme de marches ou de vagues rompt le profil uniforme de l’indlandsis et cause de profonds changements des taux d’accumulation à l’échelle locale. La distribution des accumulations sur ces irrégularités peut être prévue en utilisant la pente des surfaces et les déviations á partir de l’élévation régionale. Les profils topographiques superficiels et sous-superficiels indiquent une relation directe entre les collines sous-superficielles et les irrégularités en forme de marches, mais ils ne peuvent pas êre expliqués de manière quantitative par la présente théorie de l’écoulement de la glace. Des études détaillées dc l’accumulation effectuées en relation avec un projet de nivellement au voisinage du camp “Century” ont révélé la formation d’une irrégularité en forme de vallée peu profonde. Au sein de cette irrégularité, les taux d’accumulation ont augmenté montrant qu’elle est causée par le phénomène d’écoulement.

Zusammenfassung

Zusammenfassung

Schneeauftragssiudien auf der Thule-Haibinsel, Grönland. Ergebn fisse von Pegelrflessungen. Markierungsflächen und Schächten längs enter 136 km langen Piste über das grönländische Inlandeis auf der Thule-Halbinsel wurden zur Bestimmung der regionalen und lokalen Verteilung des Schneeauftrages herangezogen. Auf regionaler Ebene können Flächen gleicher mittlerer jahresakkumulation hinreichend zuverlässig aus einem Modell vorausbestimmt werden, das den Abstand von der Feuchtigkeitsquelle und die Meereshöhe als unabhängige Parameter benutzt. Eine Reihe von stufen- oder wellenartigen Erscheinungen unterbrechen die glatte Oberfläche des Inlandeises und verursachen starke Wechsel des Schneeauftrags im lokalen Bereich. Die Akkurnulationsverteilung über diesen Erscheinungen kann aus dem Oberflächongefälle und der Abweichung gegen die regionale Höhenlage vorhergesagt werden. Profile der Oberfläche und der Untergrundstopographie zeigen eine direkte Beziehung zwischen Erhebungen im Untergrund und stufenartigen Erscheinungen an, die jedoch quantitativ nicht aus der heutigen Theorie des Eisfliessens erklärt werden kann. Detaillierte Akkumulationsstudien in Verbindung mit einem geometrischen Nivellement in der Umgebung von Camp Century haben die Bildung einer seichten, talartigen Einsenkung aufgedeckt. Innerhalb dieser Einsenkung haben die Auftragsraten zugenommen, was auf ihre Entstehung als Fliesserschcinung hinweist.

Type
Research Article
Copyright
Copyright © International Glaciological Society 1968

Introduction

Since the construction of Camp Century by the U.S. Army in 1959 and 1960, U.S. Army Cold Regions Research and Engineering Laboratory has used the facility as a site for research and engineering studies. A long term project to determine ice movement rates and particle paths from Camp Century to their point of discharge was initiated in 1961. As part of this study, accumulation measurements, 10 meter snow temperature measurements, ice thickness determinations and an extensive survey program have been undertaken. For previously published work relating to various aspects of this project see Reference MockMock (1963), Reference BaileyBailey and others (1964), Reference MockMock (1965), Waite (1965), Reference Mock and WeeksMock and Weeks (1966), and Reference ClarkeClarke (1966). This report presents the results of accumulation studies carried out between 1961 and 1966, with emphasis upon accumulation patterns, i.e. the spatial distribution of accumulation. A later report will concern itself with temporal variations of accumulation.

Area of Investigation

The Thule peninsula is a westward projection of land, largely covered by the ice sheet and bounded by Inglefield Bugt on the north, Melville Bugt on the south, and Baffin Bay on the west (Fig. 1). Camp Century is located at lat. 77° 11′ N. and long. 61° 10′ W. near the crest line of the peninsula at an approximate elevation of 1885 m. Extensive areas of ice-free land are exposed in the western and northern parts of the peninsula while scattered nunataks appear in the south. Small independent ice caps occur in the north-western part. In many respects, the lobe of the ice sheet covering the Thule peninsula can be visualized as the Greenland ice sheet in miniature, i.e. it has a generally linear shape with a distinct crest line, a mountain barrier to the south channeling flow into large outlet glaciers, and rolling uplands to the west where the ice sheet ends in broad ramps and cliffs. Figure 1 shows the Thule peninsula and routes followed by various investigators.

Fig. 1. Location map of the Thule peninsula

Previous Work

North-west Greenland and the Thule peninsula have been the site of several glaciological investigations in recent years. Reference SchyttSchytt (1955) and Reference GriffithsGriffiths (1960) studied the regimen and movement of the Thule ramp while Reference NoblesNobles (1960) carried out a similar program on the Nunatarssuaq ramp. Reference BullBull (1958) determined mean annual accumulation values from rammsonde measurements on a traverse from the east coast to the west coast. Reference GoldthwaitGoldthwait (1960) made an exhaustive study of an ice cliff on Nordisen near its junction with the main ice sheet. Reference BaderBader (1955) studied the stratigraphy exposed in a 31 m deep pit at Site 2, 130 km east of Camp Century, and determined annual accumulation values for the period 1919–54. Reference LangwayLangway (1961), in a series of pit studies showed empirical correlations of 10 m snow temperatures, mean annual accumulation, and mean snow density with surface slopes and elevation. The most extensive investigations were made by Reference BensonBenson (1962) in a series of traverses during the period 1952–55. Benson developed techniques for the correlation of stratigraphic features between pits and introduced the facies concept to the upper snow layers. His work included the correlation of environmental parameters, i.e. mean annual accumulation and temperature, with spatial position.

Techniques And Measurements

Pit studies

Standard snow-pit studies were made at 13 to 18 km intervals along the trail. Locations of pit stations are shown in Figure 1. Each pit study included stratigraphic delineation of annual units, a depth—density profile to a minimum of 10 meters depth, measurement of 10 meter snow temperatures, and a rammsonde profile to 4 meters.

Stake and board measurements

Snow accumulation measurements were made on stakes at 1.6 km (1 mile) intervals along the trail south-west from Camp Century (south trail), for a distance of 48 km. The stake line was established in 1961 and measured annually through 1965. Marker boards were placed on the surface at 1.6 km (1 mile) intervals to provide dateable surfaces along the 88 km trail extending north-west and west from Camp Century (north trail). The boards were laid in 1962 and depths to them measured in 1964 and 1965. Marker boards were also established as part of a spirit-level net surrounding Camp Century. A total of 33 boards were laid out in a rectangular area approximately 4.5 km by 6.5 km. This array was initially established in 1963 and measurements were made in 1964 and 1966.

All stake and board measurements of snow accumulation were reduced to units of g/cm2 of water by extrapolating from depth-density profiles at pit stations. Thus accumulation, in water-equivalent units, is available for the same three-year period for the entire 136 km of trail, and for an overlapping three-year period in the immediate environs of Camp Century.

Elevations and leveling

Elevations at movement stations were established by simultaneous reciprocal vertical angle measurements tied to coastal triangulation stations of the Danish Geodætisk Institut. Elevations between movement stations were determined by barometric altimetry. The level net surrounding Camp Century was established by spirit leveling in 1963 and was releveled in 1964 and 1966. Movement station P42–9 was used as the datum for each survey. Its elevation was assumed constant and all other elevations are relative to it.

Ice thickness profiles

Ice thickness was determined over the trail system in 1964 by radar and seismic methods. For discussions and results of radar measurements see Reference BaileyBailey and others (1964), Reference WaiteWaite (1966), and Reference Rinker and MockRinker and Mock (1967). Seismic results have been presented by Reference ClarkeClarke (1966). In addition to the above measurements, a radar traverse around the perimeter of the Camp Century level net was made using the U.S. Army Electronics Research and Development Laboratory’s 35 MHz radar sounder. Earlier seismic work within the Camp Century level net has been reported by Alford (Reference MockMock, 1965).

Statistical analysis

Multiple regression techniques were used to study areal trends of mean annual accumulation, as determined from pit studies and direct measurements. The treatment was essentially the same as used in earlier papers (Reference Mock and WeeksMock and Weeks, 1966; Reference MockMock, 1967) of studies of 10 m snow temperature and accumulation distribution over the entire Greenland ice sheet. For details of the method see those papers or, for a complete treatment, Reference Krumbein and GraybillKrumbein and Graybill ([1966]).

Briefly, the method is to establish a model for predicting mean annual accumulation at a point in terms of its spatial coordinates and by means of multiple linear regression using the least-squares method to derive a predictor equation satisfying certain pre-set statistical criteria. The resulting equations can be used directly or for the construction of trend surface maps.

Results

Table I summarizes the results of pit studies and direct measurements of accumulation. Figure 2 shows in graphical form mean annual accumulation from the stake and board measurements together with surface and subsurface profiles along the trails. Some important features of Figure 2 which will be discussed further are:

  1. The general decrease in mean annual accumulation towards the interior.

  2. The wide fluctuation in accumulation rates at adjacent stations, particularly from stations S-20 to S-30.

  3. The association of subsurface with surface relief features and in turn with accumulation rates.

Fig. 2. Surface, subsurface and accumulation profiles along trails

Throughout the following analysis the stake and board data rather than pit data have been used.

Table I. Data for Accumulation Stations

Regional Accumulation Trends

In order to study regional accumulation trends, the aforementioned statistical methods were used. Rather than use geographical coordinates and elevation in the prediction model as in an earlier study (Reference MockMock, 1967), it was assumed that accumulation rates were a function of distance from a moisture source area and of elevation. It was further assumed that Melville Bugt, lying south of the Thule peninsula was the sole moisture source.

The initial model was a second-degree equation involving these two independent parameters:

where γ is the mean annual accumulation in g/cm2 year, X 1 the elevation in meters, X 2 the distance from lat. 76° N. in kilometers, and a, b i are regression constants.

The 76th parallel was an arbitrarily chosen datum crossing Melville Bugt.

The results of this analysis are summarized in the equations shown in Table II. Equations (1) through (4) are based on all stake and board measurements made during this study. Equation (5) refers to stake data only from the south trail. Equation (6) is derived from all stake and board data plus accumulation data from Reference BensonBenson (1962) over a different portion of the Thule peninsula. Equation (7) is based on an analysis of residuals from Equation (1) and will be discussed in detail in a later section.

Table II. Accumulation Prediction Equations in Which γ (Predicted Accumulation, g/cm2 Year) is Given as a Function of X1 (Elevation, m) and X2 (Distance from Lat. 76° N., km), or γD (Deviation from Predicted γ) as a Function of X1D (Deviation from Regional Elevation, m) and α (Surface Slope, deg)

Equation (1) through (3) are all derived from the original model and the same data by a deletion process, the objective being to reach a reasonable compromise between the statistical measures of fit and the size of the predictor equation. Equation (4) is included simply to indicate that the parameters distance and elevation can be used successfully in an exponential as well as a polynomial model. Figure 3 shows measured accumulation and predicted accumulation from Equation (1) plotted versus distance.

Before proceeding further the question of how closely the three years of accumulation data used in this analysis approximate the true mean annual accumulation needs examination. This problem, in turn, can be subdivided into two problems:

  • a. What is the temporal variability?

  • b. What is the spatial variability (within a few meters) ?

Fig. 3. Measured and predicted mean annual accumulation plotted versus distance from moisture source

Reference ListerLister (1961) in a study conducted at Northice, approximately 500 km east of the present study area found that within an area 30 m square the annual accumulation could vary by 100 per cent and deviate 50 per cent from the mean of 14.9 g/cm2. Assuming this represents normal variability given the wind conditions prevailing at Northice then we can use the numerical value of ±7 g/cm2 as a measure of possible error in a single measurement assuming wind conditions to be similar to those at Northice. The fact that the profiles of each year shown in Figure 4 are closely similar indicates that deviations are nearly constant with time at each station. If this is the case then temporal variability is indicated by the standard deviations at each station. Still the problem of the validity of the means from three years’ data is unanswered.

Fig. 4. Total accumulation on south trail, 1961–65

Examining Table I again shows that with a single exception, station S-55–1/4, all accumulation rates determined from pits are substantially lower than the values recorded by stake or board measurements at the same or adjacent sites. This suggests that the three-year accumulation totals may represent an anomalously high period.

Another method of gaining at least a qualitative idea of the validity of the accumulation rate measurements is to look at regression curves for 1, 2, 3 and 4 years of data respectively. Assuming that the model used would provide a perfect fit if the true mean values were available, the curves using means based on successively greater numbers of years of data should show increasingly better fits. Such is the case as shown in Figure 5 where the model used in Equation (5) of Table II was processed using successively 1, 2, 3 and 4 year means of accumulation from the trail south of Camp Century. Each of the 3 statistics shown, R, F and S indicate greater success of the model with each increase in the sample size used to determine mean values.

Fig. 5. Statistical data for south trail data

Figure 7 shows isohyets as a function of distance and elevation based on Equation (1) and (6) respectively. The different form shown by the two groups of isohyets is indicative of the limitations of this type of analysis. The data used to arrive at Equation (1) can be considered a subset of the data used to determine Equation (6). The entire body of data used in this analysis is itself only a subset of the total population, and the inherent danger of error must be borne in mind when extrapolating beyond the range of the data set, which we now proceed to do.

Fig. 6. Isohyets plotted versus elevation and distance from moisture source

Fig. 7. Isohyetal map of the Thule peninsula

In Figure 7 the results of the regional analysis have been summarized in an isohyetal map of the Thule peninsula. The position of the isohyets were determined from Equation (6) using Figure 6. In terms of precipitation regimes existing over highlands in low and mid latitudes, the pattern shown in Figure 7 appears anomalous. Ideally when an air mass approaches an elevation barrier, the lifting effect causes condensation and a resulting increase in precipitation up to some level where the absolute moisture content is diminished to such an extent that precipitation diminishes. When the elevation barrier is breached and subsidence of the air begins, precipitation rapidly diminishes as adiabatic warming brings air mass temperatures above the dew point. The actual distribution of precipitation will depend upon many factors such as absolute height of the barrier, steepness of slope and many meteorological variables. Reference BensonBenson (1962) pointed out many of the difficulties in observing the precipitation patterns over mountain barriers and observed that the Thule peninsula in many respects, represented a simplified model on which to observe such patterns.

The pattern shown in Figure 7 differs from the ideal model in two respects: there appears to be no zone in which accumulation increases with elevation and there is no abrupt change in accumulation on crossing the crest. Examination of Figure 6b indicates that there would actually be an increase in accumulation with elevation if surface slopes were great enough. As it is, the decrease in accumulation with distance completely masks the effect of elevation including the development of any noticeable discontinuity near the crest. It is probable that along the Melville Bugt coast, in the mountainous region and on the ice sheet for some indeterminate distance inland, a rapid increase in accumulation or precipitation rates does in fact occur. Once this region is passed the very gentle surface slopes existing on the ice sheet preclude noticeable orographie effects upon accumulation.

In terms of the preceding arguments the increase in accumulation northwards from lat. 77° 20′ N. becomes difficult to explain. In order of relative probabilities the most likely causes are:

  1. The model is breaking down.

  2. The influence of a secondary moisture source is becoming dominant.

  3. Increased katabatic flow funneling into this area causes warmer air masses to be lifted giving increased accumulation.

Local Accumulation Patterns

Figures 2 and 4 show the large variations in accumulation rates occurring within limited distances, particularly from mile S-15 to S-30 on the south trail. Figure 4 illustrates the persistence of these variations while Figure 2 indicates their relationship to surface relief features. Detail is not sufficiently defined to state categorically that maximum accumulation rates occur in lowest parts of the hollows, but it can be seen that high and low accumulation rates are associated with topographic lows and highs respectively.

Equation (1) through (4) can each be used to calculate the predicted accumulation γ p at each station. The residual γ D can be calculated from the observed accumulation (γ 0) as

(8)

Since R S represents deviations of the predicted from observed accumulation rates, it is a measure both of errors and of deviations caused by factors unaccounted for in the original model. Of these factors the effect of local topography superimposed upon the regional slope becomes quite evident on examination of Figure 2. From station S-15 to S-30 a series of steps break the generally smooth slope of the ice sheet. Coincident with the increasing size of surface topographic features, accumulation rates oscillate much more widely around the curve of predicted accumulation.

Qualitatively it can be seen from Figure 2 that high accumulation rates are associated with topographic lows and vice versa. The spacing of accumulation stakes is insufficient to show the fine detail over any single wave or step so it becomes necessary to use statistical methods to determine the local patterns. A method similar to that used by Reference Black and BuddBlack and Budd (1964) was resorted to, the model being

(9)

where

, (X 10) n being the observed elevation at station n, and α the surface slope at station n. The result is shown as Equation (7) of Table II. Using Equation (1) and (7) the accumulation profile over a complete surface wave can be calculated. This has been done for a prominent step on the south trail; the result is shown in Figure 8. Over this particular feature, the absolute maximum accumulation occurs at the low point of the depression while minimum accumulation occurs near the crest of the rise. Figure 8 also shows the surface and subsurface profiles over the trail section. The step-like nature of the topography is quite evident, as is the fact that, when superimposed on the regional profile, the feature is shown to be a wave.

Fig. 8. Accumulation, surface and subsurface profiles over wave feature

Origin of Wave Features

Undulating and step features have been observed in Greenland (Reference BaderBader, 1961; Reference BourgoinBourgoin, 1956) and Antarctica (Reference SwithinbankSwithinbank, 1959; Reference Black and BuddBlack and Budd, 1964; Gow and Rowland, 1966) both near the coasts and in the interiors. Recent discussions have centered on the origin of these features; whether they are the result of accumulation variations or are a phenomenon associated with bed topography and flow.

Reference NyeNye (1959) has shown that perturbations of the equilibrium profile for central Antarctica would disappear quite rapidly (0.005 to 10 years for features of wavelength 3 to 140 km) and concludes that observed features must represent an equilibrium profile reflecting the subsurface topography.

Reference Black and BuddBlack and Budd (1964) observed accumulation patterns over waves which would tend to cause them to migrate upslope. If indeed a wave does migrate upslope it then becomes a perturbation of the equilibrium profile and should dissipate with the time constants calculated by Reference WeertmanWeertman (1958). However, at the same time it causes a perturbation in the accumulation profile which may tend to perpetuate, cause continued migration, or obliterate the waves. Waves which have been observed in the interior of Antarctica (Gow and Rowland, 1966) appear unrelated to bed topography and, as postulated by Gow, may have migrated into their present position by the mechanism proposed by Black and Budd.

The waves observed in the present study are unquestionably a result of bed topography, as each surface feature is directly associated with an observable feature on the bed, as is readily evident in Figure 2. Although there has been some success in quantitatively correlating surface profiles with bed topography (Reference BourgoinBourgoin, 1956; Reference RobinRobin, 1958) more recent work (Reference RobinsonRobinson, 1966) has yielded, at best, only qualitative results.

Reference NyeNye (1959) has shown that the change in surface slope caused by a subsurface feature is given approximately by

(10)

where Δα is the slope change, α the regional slope, h the regional ice thickness, P the height of the subsurface feature, and m a constant assumed to be 2.5.

Using this expression to calculate the bottom profile from the observed surface slopes gives the profiles shown in Figure 9. The predicted profile departs quite radically from the observed, to such an extent as to be worthless. This is not surprising considering that several of the assumptions made by Nye in his derivation are violated by the physical problem represented here. It is of interest to examine these assumptions and, if possible, recalculate profiles without them. Nye assumed that

(11)

Fig. 9. Observed and predicted subsurface profiles at wave feature

and that

(12)

where Δh is the change in thickness over an obstacle and h is the regional ice thickness. Inequality (ii) does not hold in the profile shown in Figure 8 where in fact Δα > α. The bed profile can be calculated without assuming Inequalities (11) and (12) from the expression:

(13)

where P is the height of the obstacle and ϵ the deviation from regional elevation. The profile calculated from Equation (13) is shown in Figure 9 along with the profile calculated from Equation (10). It is difficult to decide which is more satisfactory, both being so far from the actual profile.

Nye also assumed the following relationship between sliding velocity u and shear stress τ at the bed.

(14)

where A is independent of stress but depends on temperature and m is a constant. For purposes of the derivation and in the absence of any other information, A has been assumed constant, but in fact, where the height of bed features is a significant fraction of the thickness of the ice sheet (9 per cent in this case), it will probably change significantly in the region of the bed features. Unfortunately the changes which might be expected in A, its magnitude decreasing with the higher bed temperatures over a feature, would only accentuate the already erroneous profile.

Present theory seems inadequate to account for the observed profiles, under the assumption that the presence of a single obstacle causes none or only negligible effects on the mass flux both upstream and downstream of the obstacle, i.e. that the mass flux of the real ice sheet q is equal to q0 the mass flux of a reference ice sheet with no obstacle. If the ice sheet is in a non-equilibrium condition and is adjusting to changed conditions, then conceivably qq 0 and the observed profile represents a transient state.

Camp Century station glaciology

Regression methods similar to those previously described, were used to analyze the measured accumulation data from the level net. Coordinates on an arbitrarily established grid were used to calculate accumulation contours. Cubic, quadratic and simple linear models were tested, and the linear model gave a sufficiently good fit (R = 0.884) to be used for further work. Trend surfaces of accumulation are shown in Figure 10 and show a substantial gradient existing across the 27 km2 area from more than 37 g/cm2 in the south-east to less than 29 g/cm2 in the north-west. Since most studies, e.g. drilling and deep pit studies, have been concentrated on the south-east side of the camp, we will consider 35 g/cm2 as the mean annual accumulation at Camp Century. Accumulation determinations from this study as well as all other measurements made at Camp Century are summarized in Table III.

Fig. 10. Topographic and accumulation trend surface map, Camp Century area

Table III. Camp Century Accumulation

Topographic maps were compiled for each of the three years in which spirit leveling was done, 1963, 1964 and 1966, in order to study any surficial features which might be present. Previous work by Reference BensonBenson (1959) in an area 130 km to the east has shown a surprisingly large amount of “topography” in an area visually undistinguished. The 1963 and 1964 surveys showed no unusual features or changes in the surface. The 1966 survey, however, revealed the development of a distinct feature in the vicinity of stations 24 and 25, a broad indentation of the contours. In Figure 10 both the 1964 and 1966 surface contours are shown in this region.

The development of such a feature would, at first thought, be the result of decreased accumulation in the affected area, probably as a result of wind scouring, initiated perhaps by some minor transient feature. The prevailing wind directions as shown in the wind rose in Figure 10 would lend credence to such a view.

A cursory inspection of the accumulation data reveals that accumulation rates are highest at the point where the surface has shown the maximum relative lowering. This is shown in the profiles of Figure 11. It appears that the development of the slight surface depression rather than being the result of accumulation deficiencies has, in fact, caused an increase in accumulation rates.

Fig. 11. Accumulation and surface profiles over valley feature at Camp Century

A topographic map of the subsurface at Camp Century has been constructed from the available radar sounding data. The map, shown as Figure 12 is considered to be the best logical fit of the several interpretations possible from the data. Using Figures 10 and 12, an ice thickness of 1370 m is predicted at the Camp Century deep bore-hole site. Recently the ice sheet was completely bored through and the sub-ice material reached at a depth of 1387.4 m (Reference Hansen and LangwayHansen and Langway, 1966), a difference of only 1.5 per cent from the depth predicted here.

Fig. 12. Subsurface topographic map at Camp Century

If the 1966 profile of Figure 11 represents the equilibrium profile then Equation (13) can be used to calculate the height of the subsurface feature causing the observed change in slope. The result predicts a hill 185 m in height and a similar calculation using Equation (10) predicts one of 244 m. It seems reasonably safe to assume that the surface feature does represent a perturbation from the equilibrium profile. Weertman’s analysis (Reference Weertman1958) of perturbations from an equilibrium ice thickness enables a calculation of the characteristic decay time t, necessary for a small departure Δh, from the equilibrium ice thickness h, to decay to a small value. The expression is

(15)

where B = (m+1) (Qu+a)/h, Q = αβ, α is the surface slope, β the bed slope, u the velocity, a the accumulation rate, D = muH/α, and λ A the wavelength of the perturbation. For this feature t is equal to 0.006 d or approximately 9 min. Since it took substantially longer than this to carry the spirit level lines over the feature, it is likely that this figure is incorrect. What does appear evident is that this feature is not the result of wind action or differential accumulation and is inexplicable in terms of existing flow theory both in terms of origin and maintenance.

Summary and Conclusions

Regression analysis enables satisfactory regional predictions of mean annual accumulation to be made from the parameters distance from moisture source and elevation. The regional pattern developed from this model is consistent with observed results and explicable for the most part in terms of general meteorological principles. Mean annual accumulation is greater than 80 g/cm2 on the southern slope of the peninsula but decreases rapidly with increasing distance from Melville Bugt and increasing elevation to less than 20 g/cm2.

Large deviations, up to 15 g/cm2, from regionally predicted accumulation values are associated with localized topographic features, in particular a series of waves or steps. Detailed accumulation profiles over these features can be predicted from an analysis of deviations from regionally predicted accumulation, deviations from regional elevation, and surface slope. Such analysis reveals highest accumulation rates associated with the topographically lowest parts of the features and vice versa. Qualitatively, these features are associated with bed topography, as shown by surface and subsurface profiles, but they cannot be quantitatively accounted for by existing flow theory because of boundary condition violations.

Contour maps showing surface elevation, subsurface elevation, and mean annual accumulation were constructed for a limited area of 27 km2 surrounding Camp Century. Mean annual accumulation shows a systematic decrease from greater than 37 g/cm2 in the south-east to less than 29 g/cm2 in the north-west. At the camp itself a mean annual accumulation of 35 g/cm2 is indicated. The contour maps of elevation have revealed the development of a small change in surface shape, best described as a valley-like feature with a relative surface lowering of over 1 m causing accumulation rates to increase within the area of the surface lowering. Bed topography provides no clue to the cause of the feature which appears to be a flow phenomenon.

Acknowledgements

It is a pleasure to acknowledge the able assistance both in the field and laboratory during various phases of this project of Professor Hernan O. Pallardo, Mr. T. C. Davis, Jr., Mr. Donald Alford, Mr. Robert Bucher and SP/4. Roger Doescher. Without the capable logistic support provided by the officers and men of the U.S. Army Research Support Group, this work could not have been accomplished. Critical reviews by Dr. C. C. Langway, Dr. W. F. Weeks, Mr. J. A. Bender, Mr. A. J. Gow and Mr, C. F. Keeler are kindly appreciated.

References

Bader, H. 1961. The Greenland ice sheet. U.S. Cold Regions Research and Engineering Laboratory. Cold regions science and engineering. Hanover, N.H. Pt. 1, Sect. B2.Google Scholar
Bader, H., and others. 1955. Excavations and installations at SIPRE test site, Site 2, Greenland, by H. Bader, R. W. Waterhouse, J. K. Landauer, B. L. Hansen, J. A. Bender and T. R. Butkovich, U.S. Snow, Ice and Permafrost Research Establishment. Report 20.Google Scholar
Bailey, J. T., and others. 1964. Radio echo sounding of Polar ice sheets, by J. T. Bailey, S. Evans and G. de Q. Robin. Nature, Vol. 204, No. 4957, p. 42021 CrossRefGoogle Scholar
Benson, C. S. 1959. Physical investigations on the snow and firn of northwest Greenland 1952, 1953 and 1954. U.S. Snow, Ice and Permafrost Research Establishment. Research Report 26.Google Scholar
Benson, C. S. 1962. Stratigraphie studies in the snow and firn of the Greenland ice sheet. U.S. Snow, Ice and Permafrost Research Establishment. Research Report 70.Google Scholar
Black, H. P. Budd, W. 1964. Accumulation in the region of Wilkes, Wilkes Land, Antarctica. Journal of Glaciology, Vol, 5, No. 37, P. 315.CrossRefGoogle Scholar
Bourgoin, J.-P. 1956. Quelques caractères analytiques de la surface et du socle de l’inlandsis groenlandais. Annales de Géophysique, Tom. 12, No. 1, p. 7583.Google Scholar
Bull, C. 1958. Snow accumulation in North Greenland. Journal of Glaciology, Vol. 3, No. 24, p. 23748.CrossRefGoogle Scholar
Clarke, G. K. C. 1966. Seismic survey northwest Greenland, 1964. U.S. Cold Regions Research and Engineering Laboratory. Research Report 191.Google Scholar
Crozaz, G. Langway, C. C. jr. 1966. Dating Greenland firn-ice cores with Pb-210. Earth and Planetary Science Letters, Vol.1, No. 4, p. 19496.CrossRefGoogle Scholar
Crozaz, G., and others. 1966. Artificial radioactivity reference horizons in Greenland firn, by G. Crozaz, C. C. Langway, Jr. E. Picciotto. Earth and Planetary Science Letters, Vol. 1, No. 1, p. 4248.CrossRefGoogle Scholar
Goldthwait, R. P. 1960. Study of an ice cliff in Nunatarssuaq, Greenland. U.S. Snow, Ice and Permafrost Research Establishment. Technical Report 39. Google Scholar
Gow, A. J. Rowland, R. 1965. On the relationship of snow accumulation to surface topography at “Byrd station”, Antarctica. Journal of Glaciology, Vol. 5, No. 42, p. 84347.CrossRefGoogle Scholar
Griffiths, T. M. 1960. Glaciological investigations in the TUTO area of Greenland. U.S. Snow, Ice and Permafrost Research Establishment. Technical Report 47.Google Scholar
Hansen, B. L. Langway, C. C. jr. 1966. Deep core drilling in ice and core analysis at Camp Century, Greenland, 1961–1966. Antarctic Journal of the U.S., Vol. 1, No. 5. p. 20708.Google Scholar
Krumbein, W. C. Graybill, F. A. [1966.] An introduction to statistical models in geology. New York, McGraw-Hill.Google Scholar
Langway, C. C. jr. 1961. Accumulation and temperature on the inland ice of North Greenland, 1959. Journal of Glaciology, Vol. 3, No. 30, p. 101744.CrossRefGoogle Scholar
Lister, H. 1961. Accumulation and firnification in North Greenland. Folia Geographica Danica, Tom. 9, p. 16374.Google Scholar
Mock, S. J. 1963. Tellurometer traverse for a surface movement survey in N. Greenland. Union Géodésique et Géophysique Internationale. Association Internationale d’Hydrologie Scientifique. Assemblée générale de Berkeley, 19–8– 31–8 1963. Commission des Neiges et des Glaces, p. 14753.Google Scholar
Mock, S. J. 1965. Glaciological studies in the vicinity of Camp Century, Greenland. U.S. Cold Regions Research and Engineering Laboratory. Research Report 157.Google Scholar
Mock, S. J. 1967. Accumulation patterns on the Greenland ice sheet. U.S. Cold Regions Research and Engineering Laboratory. Research Report 233.Google Scholar
Mock, S. J. Weeks, W. F. 1966. The distribution of 10 meter snow temperatures on the Greenland ice sheet. Journal of Glaciology, Vol. 6, No. 43, p. 2341 CrossRefGoogle Scholar
Nobles, L. H. 1960. Glaciological investigations, Nunatarssuaq ice ramp, northwestern Greenland. U.S. Snow, Ice and Permafrost Research Establishment. Technical Report 66.Google Scholar
Nye, J. F. 1959. The motion of ice sheets and glaciers. Journal of Glaciology, Vol. 3, No. 26, p. 493507.CrossRefGoogle Scholar
Rinker, J. N. Mock, S. J. 1967. Radar ice thickness profiles, northwest Greenland. U.S. Cold Regions Research and Engineering Laboratory. Special Report 81.Google Scholar
Robin, G. de Q. 1958. Glaciology. III. Seismic shootings and related investigations. Norwegian-British-Swedish Antarctic Expedition, 1949–52. Scientific Results, Vol. 5.Google Scholar
Robinson, E. S. 1966. On the relationship of ice-surface topography to bed topography on the South Polar plateau. Journal of Glaciology, Vol. 6, No. 43, p.Google Scholar
Schytt, V. 1955. Glaciological investigations in the Thule Ramp area. U.S. Snow, Ice and Permafrost Research Establishment. Report 28.Google Scholar
Swithinbank, C. W. M. 1959. Glaciology. I. The morphology of the inland ice sheet and nunatak areas of western Dronning Maud Land. Norwegian–British–Swedish Antarctic Expedition, 1949–52. Scientific Results, Vol. 3, D, p. 97117.Google Scholar
Waite, A. H. jr. 1966. International experiments in glacier sounding, 1963 and 1964. Canadian Journal of Earth Sciences, Vol. 3, No. 6, p. 88792.CrossRefGoogle Scholar
Weertman, J. 1958. Traveling waves on glaciers. Union Géodésique et Géophysique Internationale. Association Internationale d’Hydrologie Scientifique. Symposium de Chamonix, 16–24 sept. 1958, p. 16268.Google Scholar

Full text views

Full text views reflects PDF downloads, PDFs sent to Google Drive, Dropbox and Kindle and HTML full text views.

Total number of HTML views: 30
Total number of PDF views: 95 *
View data table for this chart

* Views captured on Cambridge Core between 30th January 2017 - 12th April 2021. This data will be updated every 24 hours.

Access Access

Send article to Kindle

To send this article to your Kindle, first ensure no-reply@cambridge.org is added to your Approved Personal Document E-mail List under your Personal Document Settings on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part of your Kindle email address below. Find out more about sending to your Kindle. Find out more about sending to your Kindle.

Note you can select to send to either the @free.kindle.com or @kindle.com variations. ‘@free.kindle.com’ emails are free but can only be sent to your device when it is connected to wi-fi. ‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.

Find out more about the Kindle Personal Document Service.

Snow Accumulation Studies on the Thule Peninsula, Greenland
Available formats
×

Send article to Dropbox

To send this article to your Dropbox account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Dropbox.

Snow Accumulation Studies on the Thule Peninsula, Greenland
Available formats
×

Send article to Google Drive

To send this article to your Google Drive account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Google Drive.

Snow Accumulation Studies on the Thule Peninsula, Greenland
Available formats
×
×

Reply to: Submit a response


Your details


Conflicting interests

Do you have any conflicting interests? *