Introduction

In many of the sheltered bays, inlets and fiords of the Canadian archipelago and along the Arctic coastline the micro-climate developed during the short summer season is such as to produce complete melting of the ice cover almost every year. A new ice cover forms in the autumn and attains a maximum thickness in the late spring of between about 2 m and 2.5 m, depending on the latitude. These ice covers are landfast generally and exhibit little motion except during their final stage of decay.

The factors causing the disintegration of an annual ice cover are known reasonably well. A mainly qualitative description of the stages of decay has been given by Reference ZubovZubov (1945, chapter 8) and a quantitative study of the processes occurring during the part of the decay when the ice cover is still continuous has been made by Reference LanglebenLangleben (1966). The coming of the period of 24 h daylight in the spring, coupled with the relatively high values of insolation because of cloud-free conditions, rapidly produces a modification in the snow cover on the ice and reduces the surface albedo. As the season progresses, the solar angle increases and daily totals of short-wave radiation incident on the surface increase and more and more of this energy is absorbed by the ice cover as the albedo decreases. During this time (Reference LanglebenLangleben, 1966) the air-mass temperature is moderating and ultimately climbs to the vicinity of 0° C, whereupon very rapid melting of the snow cover takes place, the albedo drops to a value of about 0.2 and enough radiation is absorbed by the ice cover to produce an ablation rate of about 4 cm per day and to modify the temperature profile in the ice cover in the direction of the zero isotherm. After about a week, a system of surface drainage canals and vertical melt holes develop, and the characteristic appearance of the summer ice cover evolves, with its melt ponds and surrounding smooth hummocks. Surface melting continues at about the rate indicated and, by the time the ice cover has been reduced to a thickness of about 1 m to 1.5 m, the network of drainage canals developed by the run-off of melt water has produced local erosion serious enough to cause the cover to break up along these flaws under the action of tides or storms.

At this stage, when the ice cover is no longer continuous, the decay process can accelerate by thermal interaction of the ice floes with the sea in spite of the decrease in insolation brought about by the prevalent cloudiness associated with the relatively warm open water. Zubov (1944) has made a quantitative study of this final stage of decay on the assumption that the radiation absorbed by the ice is negligible compared with that absorbed by the water. On this basis, the time required for complete decay is considerably longer than the period available in the Arctic, considering that about a month of above-freezing temperatures are necessary first to thin the ice cover down to a thickness where it can begin to break up under wave, wind or tidal action.

It is the purpose of this paper to show that the period of final decay of the ice cover is reduced considerably if the radiation absorbed by the ice floes is accounted for, and that complete melting of the ice floes by thermal interaction with the surrounding water is then possible during the short Arctic summer.

Model

The model used in this study is based on the following premises.

• (a) Immediately on break-up of the continuous ice cover, there exists a finite fraction of open water. The open water is likely to be established very quickly by abrasion between closely packed floes.

• (b) The energy input at the surface is determined entirely by the short-wave radiation incident on the surface and by its albedo, be it water or ice. As has been pointed out by Reference Vowinckel and OrvigVowinckel and Orvig (1966) the outgoing long-wave radiation from the Earth is very nearly balanced by incoming long-wave radiation from the clouds during the summer season. It has further been shown (Reference LanglebenLangleben, 1968[b]) that the energy budget of the ice cover during the melt period is almost completely radiation controlled, in that the turbulent energy fluxes of sensible and latent heats to and from the surface are negligibly small.

• (c) The radiation absorbed by the water is confined to a relatively shallow layer which is available for heat exchange with the ice floes around their perimeter. There is some support for this assumption in temperature and salinity profiles obtained at this season.

• (d) There is no exchange of water or ice in the area under consideration with adjacent areas. It is assumed that these are sheltered areas and can be treated as closed oceanographic systems for the present purpose.

Now, at time t, let S _w represent the area of open water over a large area S of the ice cover (i.e. the ice concentration is 1−S _w/S). Let A _w be the albedo of the open water and I the intensity of incoming solar radiation at the Earth’s surface. Then in time dt, the quantity of heat absorbed by the water is IS _w(1−A _w) dt.

If, following Reference ZubovZubov (1945), it is assumed that this energy is expended solely in decreasing the horizontal dimensions of the ice floes by melting at the edges without affecting their thickness, then by the method of mixtures the heat balance equation is

where ρ is the density of the ice, L its latent heat of fusion, h its thickness at time t, and dS _w the decrease in surface area of the floes which corresponds exactly to the increase in surface area of open water.

That is

(1)

If h is assumed constant (= h ₀), which implies that the radiation absorbed by the ice is negligible compared to that absorbed by the water, and if I is assumed constant then the Zubov solution is obtained, i.e.

(2)

where S _w0 is the area of open water at t = 0 and the equation is valid for S _w ⩽ S. This solution has found wide acceptance in the U.S.S.R. and is still cited in recent Russian literature (e.g. Reference DoroninDoronin, 1970) although Zubov was aware of its limitations.

A more realistic model of the thermal decay process is obtained by including the influence of solar radiation absorbed by the ice floes. If A _i is the albedo of the ice then, in time dt, the quantity of heat absorbed by the ice is I(S−S _w)(1−A _i) dt.

If this energy is used only to decrease the ice thickness by surface melting, then

or

(3)

It should be noted that the ice cover is assumed to be near its melting point throughout its thickness and that the loss, by radiation penetration through the bottom of the ice cover, of the available energy is neglected. Reference ChernigovskiyChernigovskiy (1963) states that about 3 to 4% of the incoming radiation is able to penetrate through pack ice of thickness 1.5 m and that the proportion increases with decreasing ice thickness. The neglect of this factor in the model would tend to overestimate slightly the rate of decay of the ice. However internal melting at the brine cells, which is produced by radiation absorption within the ice, brings about a compensating effect in that the actual surface ablation rate is less than that calculated from the model and hence the ice floes present a larger perimeter area for thermal interaction with the sea. In any case, the loss of radiant energy by penetration through the ice cover is expected to be less for annual or first-year ice than for pack ice. Because of its higher salinity, first-year ice is more opaque than multi-year ice.

Equation (1) can be solved numerically by stepwise integration in combination with the ice thickness–time relationship as given by Equation (3), even in the case where I and/or A _i vary with time. There is however insufficient knowledge of the day-to-day variations of these quantities except for point measurements at selected sites for limited periods. It is reasonable therefore to assume constant values of I and A _i based on available data, providing the full range of variation of these parameters is considered. In that case Equation (3) is integrable and

(4)

where h ₀ is the ice thickness at time t = 0.

Substituting Equation (4) into Equation (1), we obtain

i.e.

(5)

where S _w = S _w0 at t = 0.

It is customary to describe an ice cover in terms of ice concentration, 1−S _w/S. The solution in this form is

(6)

valid until time T corresponding to the complete decay of the ice cover, or to zero concentration, where

(7)

Calculations and results

Equation (6) has been used to calculate the variation of ice concentration with time intervals of one day up until the time for complete decay as given by Equation (7). The following values were taken for the constants and initial conditions in these equations.

The calculations were performed for a number of values of albedo of ice, incident short-wave radiation and initial ice thickness within the range of these parameters typical of the summer melt period.

Ice thickness

Initial ice thickness h ₀ on break-up was selected as 1 m or 1.5 m. As discussed in the Introduction, it is at this stage of the thermal decay that the ice cover is observed to break up.

Some of the results of these calculations are shown in Figures 1, 2 and 3. A family of curves of ice concentration against time has been plotted in Figure 1 for various values of the albedo of ice. The ratio of incident radiation to initial ice thickness was taken as 4 cal cm⁻³ d⁻¹ (194 W m⁻³), corresponding for example to an I of 400 cal cm⁻² d⁻¹ (194 W m⁻²) and an h ₀ of 1 m, an I of 500 cal cm⁻¹ d⁻² (243 W m⁻²) and an h ₀ of 1.25 m, or an I of 600 cal cm⁻² d⁻¹ (292 W m⁻²) and an h ₀ of 1.5 m. Similar calculations were made within the range 6 cal cm⁻³ d⁻¹ > I/h ₀ > 2 cal cm⁻³ d⁻¹ (292 W m⁻³ > I/h ₀ > 97 W m⁻³) but are not reproduced here. It is seen in Figure 1 that, as might be anticipated, the ice concentration decreases more and more rapidly as the albedo of the ice surface is reduced from 0.6 to 0.2 in steps of 0.1. This cluster of curves is compared with the Zubov solution of Equation (2) where it is assumed that no radiation is absorbed by the ice, i.e. A _i = 1.0. The enhancement of the decay process by the absorption of radiation in the ice is clearly evident.

Fig. 1. Ice concentration as function of elapsed time after initial break-up for various values of ice albedo A _i. The ratio of incident short-wave radiation I to initial ice thickness h ₀ is 4 cal cm⁻³ d^−l (194 W m⁻³).

Fig. 2. Ice concentration as function of elapsed time after initial break-up for various values of the ratio of incident short-wave radiation I to initial ice thickness h₀. The albedo of ice is assumed to be 0.4.

Fig. 3. Time required for complete decay of the ice cover after initial break-up as function of albedo of ice surface for various values of I/h₀.

In Figure 2, ice concentration is again plotted as a function of time after initial break-up. This time, however, loci are presented for a number of values of the ratio of incident radiation to initial ice thickness. These were calculated for an albedo of the ice surface of 0.4, selected as typical of mid-range values of this variable. The extreme values of I/h ₀ used in the calculations were 6 cal cm⁻³ d⁻¹ (292 W m⁻³) (corresponding say to an I of 600 cal cm⁻² d⁻¹ (292 W m⁻²) and an h ₀ of 1 m) and 2.67 cal cm⁻³ d⁻¹ (129 W m⁻³) (corresponding to I = 400 cal cm⁻² d⁻¹ (194 W m⁻²) and h ₀ = 1.5 m). If the same value of h ₀ is assumed for each of the curves then Figure 2 can be said to show the variation of the rate of decay of the ice cover with incident radiation, and vice versa. For example if h ₀ = 1 m, a change in incident radiation from 267 cal cm⁻² d⁻¹ (130 W m⁻²) to 333 cal cm⁻² d⁻¹ (162 W m⁻²) (i.e. of 67 cal cm⁻² d⁻¹ (32 W m⁻²)) reduces the time required to obtain a ice concentration by about 10 d, whereas if I is increased from 500 cal cm⁻² d⁻¹ (243 W m⁻²) to 600 cal cm⁻² d⁻¹ (292 W m⁻²) the reduction in time is only about 3 d.

Figure 3 shows the time for complete decay of the ice cover, as given by Equation (7), plotted as a function of the albedo of the ice cover, with loci for the same values of I/h ₀ as in the previous diagram. It is seen that the time for complete decay increases with increasing value of ice albedo and with decreasing value of I/h ₀. The curves are continued to A _i = 1.0 to indicate the decay times that would be obtained from the Zubov equation. It is seen that the required times for all values of I/h ₀ are about twice as large as those for the case when an albedo for the ice surface of 0.4 is used.

Discussion

The Zubov equation predicts that the time required for complete decay of the ice cover after initial break-up is greater than one month and as long as two months (Fig. 3) for representative values of incident short-wave radiation and initial ice thickness upon break-up. It is unlikely that above-freezing temperatures persist for this length of time.

To explain the observed complete disintegration of the annual ice cover in many sheltered areas of the Arctic, a modified model of the thermal decay process has been introduced. This model takes into account the influence of radiation absorbed by the ice which was not included in the Zubov formulation. Considerable reduction in the time required for complete decay, generally by about a factor of 2 if an albedo of 0.4 is assumed for the ice surface, is obtained.