Introduction

Since operations in polar regions often depend on the bearing strength of sea ice, the structure, properties, and mode of formation of this material are of much current interest. While the physical and mechanical characteristics of naturally formed sea ice are fairly well understood, development of ways for accelerating and otherwise modifying the solidification process has been limited. The object is to be able to control solidification such that ice of superior mechanical properties can be rapidly produced for engineering applications. There are many possible limitations in this area of applied glaciology;Reference Kingery ¹ one which is always present is the thermal regime in which ice is to be formed. The present study has been directed to (1) clarification of basic thermal factors limiting the rate with which ice can be formed under natural as well as artificial conditions, and (2) how freezing and processing mechanisms influence the concentration and distribution of salt in ice.

The Rate of Formation of Ice on Open Water

In problems involving solidification, the pertinent heat conduction relationships embody the constraint that the liquid–solid interface must move with a velocity proportional to the rate with which heat is withdrawn from that interface. The basic thermal problem is obviously of great practical importance in connection with the casting of fusible materials, and has received much study.Reference Adams ^{2 ,} Reference Carslaw and Jaeger ^{3 ,} Reference Stefan ^{4 ,} Reference Ingersoll ⁵

One quantity which is of paramount importance in describing the pattern of heat flow during solidification of a liquid is the dimensionless parameter C_pΔT/H where C _p is the specific heat of the solid, H is the heat of fusion, and ΔT is the maximum overall temperature drop in the system. If this dimensionless group is less than 0.3, as in the case with ice freezing at typical polar temperatures, the conduction of the heat of fusion through the solid ice controls the growth rate; the slight additional heat which must be removed to cool the solid layer below its freezing temperature, can be neglected with little loss in accuracy. Similarly, small amounts of superheat in the liquid have little effect on the freezing rate. (For metals, C_pΔT/H is often greater than unity and the conduction equations are substantially more complex than for ice.)

The first rigorous analytical treatments of the general problem were derived in connection with polar ice formation by Neumann,Reference Carslaw and Jaeger ³ StefanReference Stefan ^{4 ,} Reference Ingersoll ⁵ and others, who presented the familiar parabolic rate law stating that the thickness of ice formed should be proportional to the square root of freezing time. When freezing takes place from an ice–air interface, the parabolic rate law is predicated on the assumption that growth is limited solely by conduction of heat through the solidified layer, and ignores the resistance to heat transfer encountered at the ice-air interface (i.e. the temperature of the upper surface of the ice is taken as constant and equal to the surrounding air temperature).

To take proper account of convection heat transfer at the ice–air interface, consider that the water is initially at its freezing (liquidus) temperature, and that, during solidification, the temperature gradient is nearly linear through the solid layer. Under these conditions, the rate of heat flow through the solid ice is, at any instant, equal to that from the surface of the ice to the surrounding air, and is proportional to the solidification rate,

Rearranging (1) to eliminate T _s , we have

It is not strictly correct to assume a linear temperature distribution through the solid ice, as has been done in deriving Equation (2). In so doing, only the heat of fusion is considered, with no provision for the additional heat which must be removed to lower the temperature of the already solid ice. More thorough treatment of the problem yieldsReference Adams ²

where

.

The final term in brackets in Equation (3) may be regarded as a correction term which takes cognizance of specific heat effects; however, numerical comparison of Equations (2) and (3), as applied to calculation of ice freezing rates, indicates a maximum error of about 7 per cent (usually much less) is to be expected in using the simpler Equation (2). Integration of Equation (2) yields expressions relating ice thickness to freezing time,

and, substituting the physical constants for ice,

It can be shown that a maximum error of about 3 per cent results from ignoring specific heat effects in the derivation of Equation (5A).

When Equation (5A) is used to predict rates of ice formation, it is necessary to know the ambient temperature and the prevailing heat transfer coefficient. If the ambient temperature varies with time over the period of interest, the time average value of ambient temperature can be used rigorously since the instantaneous freezing rate is directly proportional to the overall temperature difference. Thus, the integrated parameter, θΔT (in ° C. hr.) can be substituted in Equation (5A); this parameter has also been used by AndersonReference Anderson ⁶ to correlate a large amount of experimental data. The temperature distribution in the solid ice is linear, with the liquid-solid interface at the liquidus temperature and the top surface temperature expressed as a function of ice thickness by rearranging Equation (1)

In still air, the transfer of heat from the top surface is largely by turbulent natural convection.Reference McAdams ¹¹ Although there is significant transfer by radiation, the latter is difficult to assess without precise knowledge of the effective sink temperature of the Arctic sky, and the total emissivity of ice.

where T_s and T_B are absolute temperatures of the surface of ice and the sky.

Wind in excess of 5 m.p.h. (2.2 m. sec.^–1) will amplify the convection heat transfer coefficient. In this case, any radiation contribution is overwhelmed, and the convection coefficient can be estimated usingReference McAdams ^{11 ,} Reference Jurges ^{12 ,} Reference Knudsen and Katz ¹³

The results of thermal analyses determined on ice solidifying in nearly still air at −40° F. (−40° C.) are shown in Figure 1. The curves yield an average value of the overall heat transfer coefficient

and exhibit many interesting and important features: (a) During freezing, the top surface temperature of the ice is not constant and is nearer the freezing point of sea-water than the ambient air temperature, (b) the temperature distribution in the solid ice is nearly linear,confirming that specific heat terms are relatively unimportant, and (c) during early stages of freezing, the temperature gradient in the solid does not change drastically with time and the solidification rate is nearly constant. When solidification takes place in still air, it is only after the ice thickness exceeds 40 cm. that the parabolic (conduction) term in Equation (4) exceeds the linear (convection) term; in fact, even when the ice has become 120 cm. thick, 25 per cent of the freezing time is attributable to convective resistance at the ice-air interface.

Fig. 1.

Temperature distributions during solidification of sea-water which has been brought suddenly into contact with cold air. (Measurements made at Eglin Air Force Base while freezing a deep layer; the curves only show that part of the layer which has frozen by transfer of heat upward to the air. The confined liquid suffers same increase in salt concentration and decrease in liquidus temperature with time)

The effect of air movement on the early stages of ice solidification is shown in Figures 2 and 3. Figure 2 relates ice thickness to freezing time, showing the effect of a 15 m.p.h. (6.7 m. sec.⁻¹) wind, and, in Figure 3, the same comparison is drawn in terms of freezing rates. From Reference Carslaw and JaegerFigure 3 it may be inferred even after 60 cm. of ice have formed, surface convection has substantial influence on the rate of ice formation.

Fig. 2.

Ice thickness vs. time (X ₁ vs. θ) sea-water brought into contact with air at −35° C. (Curves are drawn using Equation (5A). Points are data of AndersonReference Anderson ⁶ )

Fig. 3.

Freezing rate vs. thickness (dX ₁/dθ) for sea-water brought into contact with air at −35° C.

Data of AndersonReference Anderson ⁶ on natural ice formation, superimposed on Figure 2, indicate a rather high coefficient of heat transfer. Although this may be due to prevailing wind conditions, it should be emphasized that radiation can be important when the air temperature is near freezing, the effective blackbody temperature of the sky is substantially below the air temperature. In fact, ice may form even when the air temperature is above freezing.Reference Ingersoll ⁵

The relationship given by ZubovReference Zubov ⁷ agrees closely with Equation (5A) in form, but yields somewhat slower freezing rates. TabataReference Tabata ⁸ documents the deviation from the simple parabolic rate low by showing that the rate “constant” increases with time. Theoretical relationships given by KolesnikovReference Kolesnikov ⁹ and also SimpsonReference Simpson ¹⁰ consider factors such as snow cover and solar radiation on seasonal ice growth.

The thermal reasoning outlined in the above paragraphs has incorporated the concept of a smooth, plane liquid-solid interface, as representing the locus at which all heat of solidification is evolved. This is not the case as shown by the interfacial configuration of a core sample cut from a solidified ice sheet (Fig. 4, p. 761). When a solution freezes over a temperature range instead of at a discrete melting point, it also freezes over a spatial distribution and not at a discrete plane. To take rigorous account of the way in which heat of fusion is released over a zone of finite breadth, rather than at a plane, greatly complicates the problem, but, happily without significantly modifying the thermal picture. A complete analysis has been done for a metallic solution (aluminum-4 per cent copper alloy)Reference Adams ² with the results shown in Figure 5, where the temperature distribution which obtains during freezing of a solution is compared with that for pure aluminum. The temperature distribution exhibits a sharp discontinuity at the distinct plane liquid–solid interface in the pure metal, while there is a fairly pronounced “knee” in the region of the liquidus in the alloy. However, the overall picture is not profoundly modified by the presence or absence of the solute (copper), the main effect being upon the liquidus temperature itself.

Fig. 4.

Core sample removed during solidification showing dendritic, jagged liquid–solid interface. The sample was frozen downward by contact with cold air, and is shown upside down

Fig. 5.

Freezing of a solution compared thermally to freezing of a pure substance

In a dilute solution, such as sea-water most of the latent heat of fusion is evolved near the liquidus temperature. This is indicated by the heat content versus temperature curve shown in Figure 6; the dotted line represents the approximation embodied in the assumption of discontinuous plane-front freezing. It is to be expected, then, that freezing rates and temperature distributions, calculated as if for pure water, will closely represent the freezing of a dilute aqueous solution, provided the “freezing point” is set equal to the liquidus temperature of the solution.

Fig. 6

Heat content vs. temperature for sea-water.

Basis: the phase diagram of Thompson and NelsonReference Thompson and Nelson ¹⁷ ;

average specific heat of brine = 0.90 cal./g.° C.;

average specific heat of ice = 0.49 cal./g.° C.;

presence of other phases ignored;

heat content taken as zero at −30° C. (1.4% liquid)

Freezing of Cast Ice

When ice solidified upon open water is to be used for bearing purposes, a certain minimum thickness must be developed. It is of interest to produce this thickness as quickly as possible to maximize the useful life of the ice. Since one of the controlling factors is conduction of heat through the solid ice, solidification can be substantially accelerated by “short-circuiting” this heat flow path by transferring water from beneath the ice sheet into contact with the cold top surface and ambient air.

Consider the deposition of a single layer of cold water on top of a thick ice sheet. Heat is transferred from the water down into the cold ice, and up into the ambient air, the two processes taking place simultaneously and independently, until the new layer becomes completely solid. If the base ice on which the layer is to be deposited were initially at some uniform temperature substantially below freezing, the heat-conduction situation is that which has been studied analytically by SchwartzReference Schwarz ¹⁴ and LightfootReference Lightfoot ¹⁵ in connection with the freezing of metals in molds.

In this case, the plane at which the cold ice and water are brought into initial contact immediately assumes a temperature intermediate between the freezing point of the liquid and the initial temperature of the cold ice. This temperature remains constant throughout the freezing process, and freezing conforms with the parabolic rate law

In classical treatments of the subject, the parabolic rate constant is evaluated by numerical solution of a cumbersome transcendental equation,Reference Adams ^{2 ,} Reference Simpson ^{10 ,} Reference McAdams ¹¹ but for the case of ice formation, may be calculated with an error less than 1 per cent from

where ΔT is the temperature difference between liquidus and initially cold ice.

The temperature at the location of the original ice surface is constant and given by

where T ₀ is the initial temperature of the base ice.

Since T_i is constant, the problem may be conceptually divided into two parts. Heat flow into the base ice is described simply in terms of suddenly raising the surface temperature to T_i and holding constant. The temperature distribution in the base ice is then given by

where erf a is the “error function” of a,

, a tabulated function.Reference Carslaw and Jaeger ^{3 ,} Reference Ingersoll ⁵

Since specific heat effects are small, the temperature distribution in the layer of new ice is linear, and the rate with which heat flows through the new layer must equal the rate with which heat flows into the base ice, and also be proportional to the solidification rate. In this way, a relationship is developed which parallels Equation (1),

from which Equations (9), (10), and (11) can be derived.

In practice, when water is cast onto a thick layer of ice formed by natural solidification, the initial ice is not at uniform temperature. Instead, there is a gradient, M, from the freezing temperature at the bottom of the ice sheet to the cold top surface. When an initial linear gradient is present, if the surface of the ice is suddenly heated to T_i and held constant, the temperature distribution is given by

where T ₀ is, in this case, the initial surface temperature of the ice.

The rate of heat flow into the ice can be related to the solidification rate, giving

Equation (15) is similar to Equation (13), except that the term, KM, has been added, showing that the presence of an initial gradient in the ice sheet has the effect of slowing solidification, negligible at early times but more pronounced as solidification proceeds. Integration of Equation (15) gives

The expression, KM/ρH, which appears in Equation (16), is the rate with which ice is freezing at the bottom of the base ice. In terms of total thickness, then, the second term in Equation (16) is effectively cancelled. In other words, Equation (9) accurately reflects the total amount of new ice which forms by conduction of heat into the base ice.

Equations (14) and (16) are precise to the extent that T_i is constant, and that the ice sheet can be considered indefinitely thick. If water were cast on very thick ice, initially at uniform temperature, T ₀, then T_i would indeed be constant. When there is an initial gradient, T_i increases somewhat with time, but the effect is not large, because, in fact, T_i is never but a few degrees below the freezing temperature, and can never increase beyond the freezing temperature. For example, if T ₀ is −34.5° C. and M is 0.263° C./cm., T_i starts at −5.93° C. and increases slightly during the solidification process.

The only other restriction on Equations (14) and (16) is that the ice sheet may be considered infinitely thick so long as

, where l is the thickness of the base ice, by which time the freezing rate associated with heat flow downward in the ice sheet is so slow as to be no longer of any practical interest. (For example, if l = 120 cm., Equations (14) and (16) become invalid after 40 hr.)

The total thermal picture pertaining to freezing of a single cast layer deposited on a “cold” ice sheet is shown in Figure 7, which gives the temperature distribution prevailing 10 hr. after casting 30 cm. of water on a 120 cm. thick sheet of ice in an ambient temperature −34.5° C. The curve has been drawn using Equations (5A), (1A), (9), (10), (11), and (12), and portrays several key features of this solidification system. First, the heat released upon solidification of only 4.1 cm. of ice at the bottom of the layer has drastically heated the base ice, manifesting the relatively low heat-absorbing usefulness of “cold” ice. Second, at the end of 10 hr., the amount of ice frozen at the bottom of the new layer roughly equals that frozen from the top. However, the freezing rate from the bottom is decreasing parabolically with time, whereas that from the top is nearly constant, so that when freezing is finally complete, most of the total heat removed to consolidate the layer will have been transferred from the top surface to the air. Third, the original ice surface is heated to a temperature, T_i , fairly near the liquidus temperature. Fourth, the major surface resistance to heat flow is reflected by the relatively large temperature drop at the ice-air interface.

Fig. 7.

Temperature distribution during solidification of a 30 cm.layer of water deposited on “cold” ice. ten hours after flooding

A case of practical importance is that in which the ice sheet is thickened by sequential deposition of water layers, each of which is allowed to freeze before depositing the next layer. A freezing time for a single layer will depend upon its thickness and also upon whether or not the underlying ice is initially cold, as shown in Figure 8. If layers more than 5.1 cm. thick are considered, the first layer effectually “uses up” the heat-absorbing capacity of the base ice, and, from that time, all freezing takes place by heat transfer directly to the ambient air, provided each layer is deposited as soon as the preceding layer has become completely solid. Thus, using a predetermined sequence, such as one layer per day, Figure 8 shows how much thicker the first layer can be than the subsequent layers.

Fig. 8.

Layer thickness us. solidification time (X ₁+X ₂ vs . θ) for fleezing in still air at −35° C.

After the base ice has been heated by the deposition of the first six or eight centimeters of ice, the maximum average deposition rate depends somewhat upon the layer thickness. In Figure 9 it is shown that the maximum possible steady growth rate increases somewhat as the layer thickness decreases. The time space between consecutive deposits is given by the lower curve in Figure 8.

Fig. 9.

Maximum deposition rate vs, layer thickness (X ₁/θ vs. X ₁)

Finally it should be observed the physical picture summarized in Figure 7 reflects accurate thermal book-keeping, hut would be distorted somewhat by levitation effects. Some of the ice crystals formed as a result of beat flow downward into the base ice, will detach themselves and float to the top surface, with the effect of increasing the top thickness, X _i , and decreasing the bottom thickness, X ₂. However, neither the downward nor upward heat flow rates would be materially influenced by this physical transport process, so that the sum of the two thicknesses, as calculated using Equations (5A), (9), and (10), or Figure 8 will be quite correct.

Heating the Base Ice

The time-temperature history of the base ice during and after flooding operations is of some interest since this potentially influences mechanical strength, thermal strains and salt migration. The way in which the ice sheet heats up during the early stages of flooding has already been described by Equation (14) and portrayed in Figure 7. It has already been indicated Equation (14) is limited to values of

greater than 1.5 An expression suitable for long periods of time (limited to less than 1.5) is

When

, Equations (14) and (17) agree within 1° C.

Cooling the Ice

After flooding has proceeded for a long time, and all ice (new and old) is at a temperature not far below the melting point, it is of interest to know how the ice will cool by conduction and convection heat transfer to the ambient air, at T ₀, after flooding operations have ceased. During early stages the temperature distribution is represented by

θ is the time from start of cooling.

Again Equation (18) is valid only when

is greater than 1.5. The surface temperature during the same period is then given by

In Table II is presented a brief set of values from CarslawReference Adams ² for predicting surface temperatures.

Table II.

Surface Temperature during Air-Gooling of Ice

After cooling has progressed to the extent that

is less than 1.5, the temperature distribution in the ice is given by

where L = hl/K and P = Lπ/(1+L).

(The value of L is usually greater than 3.0. If not, then P must be computed as the first root of the equation: P cot P+L = 0. This will never be necessary for ice thicker than 60 cm.)

The surface temperature is

In Equations (20) and (21), the first terms represent the final, steady state, linear temperature distributions.

In Figure 10 are shown experimentally determined temperature distributions observed before, during, and after flooding operations by U.S. Naval Civil Engineering Laboratory personnel at Point Barrow, Alaska. These curves reflect expected thermal patterns, the important features of which can be reproduced using Equations (17), (18), (20), and (21).

Fig. 10.

Temperature distributions measured at Paint Barmy by the Naval Civil Engineering Laboratory

Since the ice structure may be injured by over-rapid cooling, it may eventually become important to establish procedures for controlling the cooling rate, thereby minimizing thermal strains, as is done in the heat treatment of ceramics or metals. For example, if field studies indicate it is desirable to place a maximum value on the allowable surface temperature gradient, one which is substantially less than with natural air cooling, this can be achieved by programmed “tapering-off” of the flooding operation. Thus, if cooling is to be done with a constant prescribed surface gradient, the rate of heat flow corresponding to this gradient is subtracted from the rate of heat flow by convection from the surface, and the difference made up by freezing water on the surface

where M′ is the maximum allowable surface gradient, and V″ the instantaneous rate of water supply.

The surface temperature under this defined circumstance of controlled cooling is given by

which, substituted into Equation (22), relates the required water supply to time.

Scrutiny of Equation (22) shows that after

the water supply can be ceased altogether, and the ice will continue to cool with shallow gradients.

Short Cycle Flooding

From the standpoint of trying to optimize structure, salinity, and mechanical properties, some work has been done using flooding times the order of minutes rather than hours, with substantial solid cooling interposed between short flooding cycles. Under these conditions, virtually all heat flow is by conduction downward into the cold base ice, and the initial temperature gradient in the ice sheet, if it is more than 90 cm. thick, has little influence on the solidification rate or thermal fluctuations which result. In this case the time-temperature behavior of the base ice is typified by the experimental curve shown in Figure 11. Since freezing from the top surface is negligible in this circumstance, an excess of water can be used, to be removed by squeegees or rollers after the prescribed flooding time has elapsed. The solid ice then cools primarily by conduction downward. The surface temperature drops very rapidly between cycles. After the first cycle, the surface temperature is given by

where θ is the total time after initiation of flooding, and θ ₁ the flooding time. The surface temperature cools only slightly more slowly after the second cycle, the third cycle, etc., and time-temperature cycles such as that portrayed in Figure 11 can be developed theoretically by making use of the concept of plane, instantaneous heat sources in the infinite medium.

Fig. 11.

Time-temperature variations at a fixed point in short cycle flooding. (

cm. per cycle; thermocouple initially on surface; “on” refers to duration of liquid contact, and “off” is cooling period)

With respect to solidification rates, it has been found experimentally and confirmed theoretically, using flooding times of 4 min. separated by cooling periods of 16 min., a deposition rate of about 0.3 cm. per layer can be sustained over a four- or five-hour period in an ambient temperature of −35° C. to −40° C.

Brine Segregation

It is widely recognized the solidification of sea-water is accompanied by complete segregation of dissolved salt on a local scale (microsegregation) and fairly pronounced macrosegregation as well. Microsegregation is the necessary result of the insolubility of salt in solid ice, and macrosegregation is promoted by several factors. (1) Ice crystals are of lower density than water and float. This effects a partial separation and is one reason “young” sea ice has a salinity substantially below that of the water from which it solidified. (2) Concentrated brine, produced in the immediate region of solidification is heavier than the more dilute sea-water, and tends to sink by natural convection. (3) During the freezing of brine, the expansion associated with solidification displaces brine-rich, interdendritic liquid, downward toward the water, or up onto the ice surface. (4) Finally, a small entrapped region of liquid brine, surrounded by solid ice, such as may obtain in late stages of solidification, will migrate through the solid ice if there is a temperature gradient, the direction of movement being from cold to warm.

Not all the above effects are quantitatively well-understood, especially since all take place simultaneously. The reasons for migration of brine droplets through solid ice are fairly clear, however. If there is a temperature gradient, and the droplet tries to come into equilibrium with the ice, there must result, in the droplet, a concentration gradient, the cold side being richer in salt than the warm side,

where dc/dx is the concentration gradient, m′ the temperature gradient, and b the slope of liquidus line on the temperature–composition phase diagram for sea-water. From this migration velocity can be estimated, since diffusion of salt to the warm side of the drop results in melting at a rate which keeps the salt concentration from increasing on that side,

(In Equation (26), the liquidus is treated as a straight line.) Combining Equations (25) and (26),

The diffusivity of salt in water is known from electrochemical data to be the order of 2.4 × 10⁻⁵ cm.²/sec., and, selecting for the temperature gradient 0.5° C./cm. and temperature −4° C., the migration velocity is roughly 3 × 10⁻⁶ cm./sec. or 0.3 cm./day, a magnitude consistent with experimental observations. As either the temperature or temperature gradient increases, so does the migration velocity.Reference Whitman ¹⁶

The observed salt distribution in a three-inch thick layer of sea ice, solidified on open water in contact with still air at −40° C. is shown in Figure 12. However, the distribution shown in Figure 12 holds only for that particular instant in time, because mass transport continues to take place throughout the semi-solid layer as solidification progresses, such that the average salt concentration throughout the thickening layer decreases with time, as shown in Figure 13. To take advantage of its relatively low salt content, newly formed ice could be utilized by freezing a surface layer in a sea-water pond, and then mechanically separating the solid ice for deposition elsewhere. This would have the added advantage of separating the heat transfer and deposition sites, and would markedly reduce heating of the base ice.Reference Thompson and Nelson ¹⁷

Short-cycle flooding, removing excess water with squeegees or rollers, such as described in the last section, also gives ice of low salinity as illustrated in Table III. However, the tendency for ice to float and also dendritic freezing interfere with the use of short-cycle flooding as a method for efficient desalination.

Fig. 12.

Salinity distribution in ice layer forming under still air at −40° C.

Fig. 13.

Average salinity vs. freezing time for young sea ice freezing in still air at −40° C.

Table III.

Salinities

Summary

1. When ice forms on open water by transfer of heat to the atmosphere, convection at the ice-air interface is of controlling importance in the rate of ice formation. This is especially true during early stages when the thickness is less than 60 cm. Analytical treatment of the problem has yielded expressions for temperature distribution and ice thickness as functions of time which are in accord with experimental findings, and which are greatly simplified by the fact that specific heat effects are of secondary importance in this solidification system.

2. Limiting rates and times have been derived for the freezing of sea-water cast in sequential layers on top of initially cold sea ice. The effects of the initial temperature gradient in the base ice have been evaluated, and it is shown that a single layer more than 8 or 10 cm. thick will “heat” the base ice to such an extent that if flooding operations are continued heat transfer downward into the solid ice is retarded.

3. Equations are presented which facilitate analyses of the way in which the base ice is heated and cooled, respectively, during and after flooding operations.

4. Qualitative review of the many factors influencing the segregation and migration of salt during and after solidification of sea ice has suggested some highly complex mass transfer phenomena which are not well understood and cannot even be ranked according to their relative importance. Much fundamental study will be required to bring the state of knowledge on the chemistry of solidification to the same level as that of thermal phenomena.

5. Brine segregation during freezing of sea-water is substantial as a result of the insolubility of salt in solid ice. This segregation can be used to effect partial desalination either by freezing a surface layer which can be mechanically separated or by short-cycle flooding with the removal of excess brine by rolling.

Acknowledgement

This study was carried out at the Arctic Research Laboratory, Point Barrow, Alaska, and at the Climatic Research Laboratory, Eglin Air Force Base, Florida. The support and assistance of Dr. Max Brewer and Mr. C. V. Rychlik of those activities is gratefully acknowledged. Several helpful discussions were held with Colonel L. DeGoes, Arnold Funai, and Earl Moser. The research reported in this article was sponsored by the Geophysics Research Directorate, Air Force Cambridge Research Center, under Contract AF 19 (604)—3073 with the Arctic Institute of North America.