Introduction
The flow of air and transport of heat in snow are important processes affecting remote-sensing applications and chemical-species transport. In addiüon, the possibility of ventilation in polar firn has important implications for interpretation of polar ice-core records, since airflow through the firn could affect the local firn temperature and chemical processes. Field measurements of firn temperature in Greenland have been made that support the possibility of significant ventilation. Reference Bader, Haefeli, Bucher, Neher, Eckel and ThamsBader and others (1954) found that the near-surface snow-temperature changes were faster than could be explained by heat conduction, and attributed the differences to heat convection by moving air. More recently, Reference Stearns and WeidnerStearns and Weidner (unpublished) measured snow temperature at GISP2 in Greenland and hypothesized that ventilation may be occurring to depths of 4 m.
There are several theoretical treatments of ventilation in the literature. Reference ColbeckColbeck (1989) found that pressure drops created by surface topography can cause air movement in snow. Reference Clarke and WaddingtonClarke and others (1987, 1991) presented one-dimensional and three-dimensional theories and showed that the near-surface turbulence can also induce airflow in firn, and that viscous-dissipation effects in the firn were insufficient to cause significant thermal disturbances. Reference Albert and McGilvaryAlbert and McGilvary (1992) showed that dissipation effects are several orders of magnitude lower than those of advection, and that the controlling factor in thermal effects of snow ventilation is the balance between heat conduction through the snow and heat advected by the airflow. This was demonstrated in a one-dimensional situation where air was forced through a column of snow with fixed end temperatures: the heat advected by the airflow had a significant impact on the resulting temperature distribution. For that situation the numerical results showed excellent agreement with the experimental data.
The purpose of the current investigation is to obtain preliminary estimates of the thermal signature of ventilation by conducting some two-dimensional numerical calculations. Since most of the existing field evidence of ventilation in polar firn comes from temperature sensors, it is important to examine the relationship between airflow patterns in the firn and the resulting heat balance.
Governing Equations
The model employed here is a modified form of that described by Reference Albert and McGilvaryAlbert and McGilvary (1992). The air is treated as an incompressible medium, where conservation of mass is
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20170922041826611-0679:S0260305500011435:S0260305500011435_eqn1.gif?pub-status=live)
where xj is a coordinate direction and vj is the velocity of the airflow through the snow (j = 1,2 for the current two-dimensional investigation). As described by Reference Albert and McGilvaryAlbert and McGilvary (1992), dissipation and compressibility heat-transfer effects are negligible, so that the assumption of incompressible flow is valid. It was also shown in that paper that vapor transport has a small effect on the heat transfer and the resulting firn-temperature profile, when compared to the more dominant control of firn temperature by heat conduction and by the advection of heat due to airflow. (The inverse is not true, however, since the temperature profile can have a dominant influence on the mass-transfer effects of vapor transport.) Therefore, it is safe to neglect vapor-transport effects on temperature in the current heat-transfer analysis. The airflow is assumed to follow Darcy’s law, where the Darcy velocity, vj , is given by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170923125249-52753-mediumThumb-S0260305500011435_eqn2.jpg?pub-status=live)
where Ρ is pressure, K jk is the permeability tensor, and μ is the viscosity of the air. Repeated indices imply summation. Equations (1) and (2) are combined to give Equation (3) for airflow:
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170923125249-45543-mediumThumb-S0260305500011435_eqn3.jpg?pub-status=live)
The heat transfer is described by the advection-diffusion equation:
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170923125249-86173-mediumThumb-S0260305500011435_eqn4.jpg?pub-status=live)
where
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170923125249-28815-mediumThumb-S0260305500011435_eqn5.jpg?pub-status=live)
In these equations Τ is temperature, C p is specific heat, λjk is the thermal conductivity tensor, ρ is density, Φ is porosity, the subscripts s, a, and i refer to snow, air, and ice, respectively, and repeated indices j and k imply summation.
Equations (3) and (4) are solved numerically using a finite element model with linear triangular elements developed by W. R. McGilvary and M. R. Albert at the Cold Regions Research and Engineering Laboratory. The simulations described below will focus on the effect of the two-dimensional airflow velocity field within the firn on the temperature distribution there. The airflow in the firn will be driven by surface-pressure effects. The temperature distribution will include the effects of fixed surface temperatures and advective effects due to air motion within the firn.
Numerical Simulations
In the one-dimensional analysis, Reference Albert and McGilvaryAlbert and McGilvary (1992) showed that airflow through a cylindrical snow sample with fixed temperatures at the inlet and outlet could have very pronounced effects on the temperature regime, with the ventilation effect of decreasing the temperature gradients near the inlet and increasing them near the outlet. In that case, air (saturated with vapor at the inlet) was forced through a cylinder of snow, where the inlet temperature was maintained at a lower temperature than the temperature at the outlet. It was also shown that a strong overall temperature gradient (imposed by the thermal boundary conditions) could drive the heat conduction enough to counteract the thermal effects of ventilation. Thus cases with significant airflow through snow may not be detectable from temperature measurements. In two and three dimensions, the airflow patterns are much more complicated than in the one-dimensional case. It is of interest now to investigate whether the multidimensional thermal signature is similar to the one-dimensional case; i.e., whether ventilation decreases temperature gradients near the surface and increases them deeper in the firn.
The simulations consider a segment of firn 6 m long and 6 m deep. For these estimates, the firn surface profile is assumed to be flat, with various forms of induced pressure fluctuations on that surface. A finite-element mesh with 1056 triangular elements, becoming successively finer mesh from the bottom to the top surface, was employed for the simulations. In each of the following simulations, the region is assumed to be at a uniform initial temperature.
The case considered is the simplest possible case: homogeneous firn. The density is 300 kg m−3 and grain size 1 mm, which gives a permeability equal to 7 × 10−9m2 (Reference ShimizuShimuzu, 1970), and thermal conductivity 0.25 Wm−1 Κ (Reference MellorMellor, 1977). For this application, the permeability and thermal conductivity are taken to be uniform throughout the domain. The initial temperature is −30°C throughout, with surface step change in temperature to −25°C.
The first simulation imposes a steady-state, sinusoidal pressure of the form Ρ = A sin(2πx1/L), with A = 10 Pa and L = 1.09 m across the top boundary, as depicted in Figure 1. This peak amplitude pressure of 10 Pa is taken from measurements conducted by Reference Farrell, Greacen and GurrFarrell and others (1966), who put an array of pressure transducers over a stretch of soil with short (3 cm) grass, and found that the amplitude of the difference in horizontal pressure measurements increases with the spacing of the measurements. The resulting steady-state isotherms are shown in Figure 2a. It should be noted that the dips in temperature near the side boundaries result from the use of a zero-gradient boundary condition on the pressure solution used there for convenience; the center 3 m is the primary region of interest. Simulations were run with a denser grid and with wider x i-dimensions, but the temperatures in the center 3 m were the same as depicted here. For comparison, the isotherms resulting from heat conduction without ventilation are shown in Figure 2b. Comparison of Figure 2a and 2b shows that ventilation does tend to unify the surface temperatures in the top tens of centimeters of the snow, decreasing the temperature gradient there, as was also true in the one-dimensional case. The increase in temperature gradient at lower depth is much less significant than in the one-dimensional case, however.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170923125249-24044-mediumThumb-S0260305500011435_fig7g.jpg?pub-status=live)
Fig.1 . Pressure forcing on the surface of the firn.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170923125249-72157-mediumThumb-S0260305500011435_fig8g.jpg?pub-status=live)
Fig.2 . Isotherms resulting from: a, pressure amplitude of 10 Pa and an overall temperature gradient of 0.83°C m−1 ; b, heat conduction alone, without ventilation, for an overall (top-to-bottom) temperature gradient of 0.83°C m−1.
The airflow velocity field is two-dimensional and thus each velocity vector has two components, but it is the scalar magnitude of the resultant velocity vector at each location that is important in chemical reactions. In Figure 3 the magnitudes of the resultant velocity vectors are plotted for x 1 locations at five selected depths. The magnitude of the velocity varies from approximately 0.02 ms−1 at the very top of the firn, to 10−3m s−1 at 0.6 m depth, to near 10−5m s−1 at 1.5 m depth.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170923125249-75693-mediumThumb-S0260305500011435_fig9g.jpg?pub-status=live)
Fig.3 . Magnitudes of the resultant airflow velocity vectors at five selected depths in the firn, for a surface-pressure amplitude of 10 Pa.
For an indication of the effect of the overall temperature gradient on the temperature distribution with this same airflow regime, the resulting temperature regimes are plotted in Figure 4a and 4b for the same snow but with surface temperatures of −27° and −10°C, respectively. It is evident that the more severe the overall temperature gradient, the more obscure are the thermal effects of ventilation, becoming confined to the surface portions of the pack. Thus, although the airflow patterns in two dimensions are much different from those in the one-dimensional simulations, the general trends of thermal signatures, creating near-isothermal temperature profiles in the top regions, are similar. Field observations of ventilation through temperature measurements are most likely to be observed when the firn temperature at depths near 6–10 m are close to the air temperature, since steeper overall temperature gradients can mask the thermal effects of ventilation.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170923125249-99424-mediumThumb-S0260305500011435_fig10g.jpg?pub-status=live)
Fig.4 . Isotherms for surface-pressure amplitude of 10 Pa. a, overall temperature gradient of 0.50°C m−1; b, overall temperature gradient of 3.3°C m−l.
Clearly the amplitude of the surface-pressure forcing will influence the airflow velocities. For the same snow properties as above, and surface temperature of −25°C, the simulation is run with a peak pressure amplitude difference equal to 1.0 Pa. The magnitudes of the airflow velocities are one-tenth of those calculated with peak amplitude 10.0 Pa, so that the flow velocity in the top several millimeters is approximately 0.002 m s−1, decreasing to a flow of 10−4m s−1 at a depth of 0.6 m, still a significant magnitude when compared to diffusive transfer mechanisms in nonconvective processes. This reduction in flow velocity from the first case means that conductive heat-transfer processes may dominate the temperature profile. In Figure 5 the resulting temperature distribution is shown. Indeed, by comparison with Figure 2b, we see that the temperature profile is practically the same as that of pure conduction with no ventilation; in this case the thermal effects of conduction are able to overwhelm those of advection.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170923125249-51158-mediumThumb-S0260305500011435_fig11g.jpg?pub-status=live)
Fig.5 . Isotherms resulting from a surface-pressure amplitude of 1.0 Pa and overall temperature gradient of 0.83°C m−1.
For all of the above simulations, the surface-pressure forcing varies over distance but not in time. In nature the surface-pressure forcing imposed by turbulent winds will be variable in both space and time. While a sinusoidal wave train representation is a simplification, it is useful to investigate the thermal effects for several frequencies in time of surface-pressure forcing. The final simulation assigns a surface pressure that varies both in space and in time according to
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170923125249-80329-mediumThumb-S0260305500011435_eqn12.jpg?pub-status=live)
where A = 10 Pa, w is frequency, and t is time (in s). In Figure 6 the vertical temperature profile in the center (x 1 = 3 m) of the domain at time t = 1 h is plotted for surface pressure frequencies of w = 0.1, 1.0, and 10.0 Hz, along with the temperature profile that would result from heat conduction alone (without ventilation). Several points are evident. First, varying the frequency in time of the surface-pressure forcing between 0.1 and 10.0 Hz had no efTect on the resulting temperature distribution after one hour of simulated time. This may be an indication that the temperature profile during ventilation is fairly insensitive to the frequency of the pressure forcing, as long as there is sufficient surface-pressure amplitude to move the air about in the top decimeters of the snow. A more exhaustive study of this preliminary result is forthcoming.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170923125249-22918-mediumThumb-S0260305500011435_fig13g.jpg?pub-status=live)
Fig.6 . Vertical temperature profile at x, location 3 m, for a surface-pressure amplitude 10 Pa, surface-pressure frequencies 0.1, 1.0 and 10.0 Hz, and an overall temperature gradient 0.83°C m−1. Also plotted for comparison is the linear, steady heat-conduction profile.
Secondly, for this overall temperature gradient (0.83°C m−1), the local temperature gradient was most affected by ventilation in the top meter of firn. Below the top meter, the temperature shift is dominated by heat conduction of the temperature changes induced in the top meter. We see that, although ventilation from surface-pressure forcing induces significant airflow primarily in the top several meters of firn, the temperature is nevertheless altered at deeper locations by conduction, because of the airflow-induced temperature changes within the top meter.
Conclusions
Numerical simulations of two-dimensional ventilation and heat transfer in a uniform firn show that spatially varying surface pressure can induce airflow velocities of 10−5 m s−1 at 1.5 m depth, and at higher velocities closer to the surface. The heat transfer results generally agree with earlier one-dimensional conclusions (Reference Albert and McGilvaryAlbert and McGilvary, 1992) that the thermal effects of ventilation are to decrease the temperature gradient in the top portions of the pack, and that the depth to which the temperature gradient is changed by ventilation depends on the flow velocity and also on the overall temperature gradient. Field observations of ventilation through temperature measurements are most likely to be observed when the firn temperature at depths on the order of 10 m are close to the air temperature, since steep temperature gradients can mask the thermal effects of ventilation. Preliminary indications are that the temperature profile during ventilation is fairly insensitive to the frequency of the pressure forcing for pressures with a frequency range 0.1–10.0 Hz, as long as there is sufficient surface-pressure amplitude to move the air about in the top tens of centimeters of the snow.
This investigation focused on uniform firn, while in nature, polar firn is layered. Ongoing work will focus on effects of layering and surface relief on ventilation processes, and will expand the analysis to include vapor and chemical-species transport.
Acknowledgements
I wish to thank Randy McGilvary for assistance in implementing the pressure boundary conditions in the code, and thank Sam Colbeck, Bert Davis, and anonymous reviewers for their technical reviews. This work was funded by U.S. Army work unit 4A161102 AT24 SS-E09, Cold Regions Surface-Air Boundary Transfer Processes.