next up previous
Next: Conclusions Up: Brine Migration and Heat Previous: Diffusion-Driven Migration

Brine Convection

The laboratory experiments on brine pocket migration were on a scale that is about one-tenth that of our in situ experiment. This small scale might remove some of the larger brine features, and not allow the feeding of brine pockets into brine channels that is observed in sea ice (Niedrauer and Martin, 1979; Weeks and Ackley, 1986). Gravity drainage of brine is associated with brine tubes or slots, open at the bottom to sea water (Lake and Lewis, 1970; Weeks and Ackley, 1986). Observations of fluid and salt fluxes in brine channels (Niedrauer and Martin, 1979) in sea ice suggest that the flow rates and frequency of occurrence of these channels are high enough to produce significant heat and salt fluxes. Our data shows evidence of the transient effects of brine tubes on temperatures. Fig. 7 is a closer view of temperatures from five adjacent thermistors near the bottom of the ice, showing evidence of a brine tube that extends about 300mm into the ice from the sea, lowering and raising temperatures locally.

Cox and Weeks (1981) parameterize gravity drainage from sea ice based on temperature gradients. However, in their model, gravity drainage stops when the brine volumes fall below 50 tex2html_wrap_inline190 , which corresponds to temperatures below -5 tex2html_wrap_inline156 C. Our data analysis excludes temperatures above -5 tex2html_wrap_inline156 C, so that gravity drainage would not be significant according to their model. Brine expulsion remains as a likely mechanism -- however, we do observe evidence of nearby brine tubes at temperatures below -5 tex2html_wrap_inline156 C, as noted above.

Here we consider convection in brine tubes, closed at both ends. For a vertical cylinder of fluid heated from below the Rayleigh number is (Niedrauer and Martin, 1979)

equation95

where tex2html_wrap_inline390 is brine density, D is the molecular diffusivity of salt, tex2html_wrap_inline394 is the dynamic viscosity of the brine, g is the gravitational acceleration and r is the radius of the cylinder. The critical value of Ra for the onset of convection is 68. If Ra is below this value, there is no convection, and diffusion of salt is the dominant mechanism. Above this value, convection begins, driven by the buoyancy differences generated by the salinity gradient. We calculate that for the properties of brine at -10 tex2html_wrap_inline156 C and a temperature gradient of 10 tex2html_wrap_inline156 C m tex2html_wrap_inline186 , the critical Rayleigh number is exceeded when the diameter of the brine tube is greater than 1mm. This is larger than most brine pockets.

Many brine tubes are tilted at an angle to the vertical (Niedrauer and Martin, 1979), which alters their convective behaviour. Criminale and Lelong (1984) perform a stability analysis of tilted brine drainage channels, and find two optimum angles of the channel for convection or overturn driven by brine density gradients, one being vertical and the other being at about 45 tex2html_wrap_inline156 . This suggests that there is a natural preference for either vertical or 45 tex2html_wrap_inline156 tubes, since at these angles the instability driving convection is strongest. All of the above analyses predict no convection in the brine tubes at small temperature gradients.

However, Woods and Linz (1992) have shown that a slot or fracture that is inclined at some angle to the vertical, containing fluid of different thermal conductivity to the surrounding rock, will convect at any Rayleigh number, even if the thermal gradient in the rock is stabilizing, due to tilted isotherms. They find an exact solution, which becomes singular as a critical Rayleigh number is approached. Vertical temperature gradients impose brine density gradients in brine tubes, a case covered in the work of Woods and Linz. However, their analysis (and that of Crimimale and Lelong, 1984) does not include the important effect of the phase change at the tube face, or allow for (or predict) the movement of the tube through the ice. Nevertheless, the mechanism driving the instability, that of tilted isotherms and isopycnals due to the different thermal conductivities of the brine and the ice, is present in sea ice, and we expect qualitatively similar results. That is, there would be no critical Rayleigh number for the onset of convection in titled brine tubes, and convection would be observed at all gradients and in any tube no matter how narrow.


next up previous
Next: Conclusions Up: Brine Migration and Heat Previous: Diffusion-Driven Migration

Mark McGuinness
Tue Apr 11 17:10:31 NZST 2000