A YEAR-BY-YEAR dating of an ice core is identical with establishing a record of in situ annual layer thicknesses that can be turned into an accumulation-rate record by correcting for (1) density variations, (2) accumulation-rate deviations up-stream, and (3) total vertical strain since the time of deposition. The first correction never causes any problem, and the second one is zero at the summit of an ice sheet and usually negligible close to ice divides. In other cases, up-stream variations in annual accumulation λH must be corrected for by using the present surface distribution.
As to the total vertical strain of a given annual layer at a considerable distance from the bottom, it may be calculated from a directly measured surface strain-rate if it is reasonable to assume a constant vertical strain-rate since the time of deposition. Otherwise, the correction implies two- or three-dimensional flow modelling.
Reference Hammer,, Hammer,, Clausen,, Dansgaard,, Gundestrup,, Johnsen and Reeh.Hammer and others (1978) have measured annual layer thicknesses, mainly by δ(18O) analyses along three 400 m ice cores recovered under the Greenland Ice Sheet Program (GISP) from Crête (on the ice divide) and Milcent in mid-Greenland, and from Dye 3 in South Greenland (see figure 4 facing p. 12). Below, their data will be transformed into series of accumulation-rates as outlined above.
Long-term trends. A surface strain net has been measured at Crête by Reference Karsten, and Stober,Karsten and Stober (1975). Using their raw data, we find the vertical surface strain-rate to be –1.32 X 10-4 and –1.14 X 10-4a-1up to 8 km east and west of the ice divide, respectively. The average value, (doubt)H= – 1.23x10-4a-1, is assumed to have been constant throughout the last 1426 years spanned by the Crete δ record. Correcting the 1 426 annual layer thicknesses accordingly leads to a λH time series, which shows a small linear trend of –4±2% per millenium, the mean λH value through the last millenium being (o.289/o.oo2) m of ice per year, Footnote * the uncertainty including an estimated dating error of 0.5%.
Unfortunately, no data for calculating surface strain-rates are available from Milcent. Instead, we assume that the flow pattern up-slope from Milcent is consistent with Reference Hammer,, Hammer,, Clausen,, Dansgaard,, Gundestrup,, Johnsen and Reeh.Hammer and others'(1978, p.7) equation (2), hbeing proportional to the ice thickness H that has been measured by radio-echo sounding (private communication from P. Gudmandsen). Reference Philberth, and Federer,Philberth and Federer's (1971) procedure leads to an estimate of h= 33° m at Milcent, which is most likely to be an upper limit for h, since orientation of the c-axis in favour of easy glide has not been allowed for. A two-dimensional flow model is then used to interpret the measured λ-profile in terms of accumulation rates at the up-slope sites of formation of the individual core layers. These accumulation rates are finally compared to the present up-slope λH values (using Reference Benson,Benson's (1962) λH gradients combined with our own surface λH mean values at Milcent and Crête), and the relative deviations are interpreted as a climatic information .Assuming o <h < 330 m, the linear λH trend comes out as (0.0±3.5)% per millennium, which does not disagree with the (–4±2)% obtained independently, practically speaking with the Crete data. The relatively high uncertainty is due, first to the Milcent λH series being shorter and, secondly, to the vertical strain and the up-slope λH corrections being more complicated than in the case of Crête. The weighted mean linear trend in mid-Greenland is thus (–3±2)% per millenium. The mean λH value through the last 796 years has been (0.54°±0-004) m ice a-1at Milcent.
At Dye 3, either of the procedures described above is even more inaccurate, because (i) no surface strain data are available, (ii) the up-slope λH gradient is not well known, and (iii) the variability of the measured λ profile is high. This induces a high standard error on the long-term λH trend that can only be given as (+2±6)% per millenium, assuming o < H < 400 m. Accordingly, the mean λH value through the last 728 years, has been (0.535 ±0.015) m ice a-1at Dye 3.
Intermediate-term variations.In Figure 1 , the heavy curves showλH time series from Dye 3, Milcent and Crête. The linear trends have been removed, and the residuals have been smoothed by a low-pass digital filter that removes all λH oscillations with periods shorter than 120 years. The λH scales have been chosen so that a given amplitude corresponds to the same relative deviation from the mean in each curve. Maximum deviation from the long-term trend lines is 11 %, 5% and 4%, respectively.
There is an obvious correlation between Milcent and Crête prior to A.D. I 800 and between all three curves prior to A.D. 1600. After A.D. ,600, Dye 3 and Milcent vary in antiphase Any attempt to explain these variations should account for the facts (i) that the bulk of snowfall comes from south-east at Dye 3, but from south-west at Milcent, whereas Crête probably receives precipitation from either direction, and (ii) that the increasing occurrence of sea ice the seventeenth century may have changed the general circulation pattern in the polar atmosphere.
No obvious correlation exists between the smoothed δ (to be published elsewhere) and λH series. For example, the Dye 3 δ curve contains a broad minimum from A.D 1500 to 1700, corresponding to the culmination of the "little ice age", whereas λH at Dye 3 decreases somewhat from A.D. 1550 to 1720. A possible explanation might be that in this period an increasing percentage of the cyclones moved into Davis Strait instead of going east of Greenland. This is consistent with the maximum of accumulation at Milcent around A.D.1700.
Short-term variations. The thin curves in Figure 1 depict the λH series smoothed by a 30 year low-pass digital filter. The deviations from the heavy curves correspond to λH series smoothed by a band-pass filter transmitting oscillations in the 30 to 120 year range. These deviations are in close correlation for Milcent and Crête: R= 0.69, significant at the P99% level with an estimated number of degrees of freedom n =25 (supported by coherence calculations). The Dye 3 short-term deviations are less correlated with those at Milcent (R- 0.30;P80%) and Crête (R= o.5o;P95%)· Nevertheless, the Milcent and Crête short-term curves are occasionally out of phase, for example around A.D. 1900, the sixteenth century, and A.D. 1180-1280, where the phase shift corresponds to 15 years. This should not be ascribed to dating errors: Looking at the time interval A.D. 1177-1276, Figure 2 shows cross-correlation coefficients r between the Milcent and Crète annual λH series (upper section) and mean annual δ series (lower section), plotted as functions of imposed time lags between the series ranging from –24 to +24 years. The standard deviation of r is o.1. Only zero time lag gives significant r–values in both series. We consider this as strong evidence that the Milcent and Crête time scales do not deviate from each other in the interval A.D. 1177-1276 and, consequently, that the phase difference between the intermediate-term λH oscillations of up to 15 years is real.