ANALYSIS OF SOLAR HEATING OF A SNOWPACK
Following from: Computational models for combined heat transfer in snow and ice
Consider first a model problem to understand the main special features of solar radiation transfer in a snowpack. Obviously, a contribution of the direct solar radiation to the snowpack heating is predominant as compared to the sky radiation. Therefore, it is sufficient to consider the absorbed power of solar radiation at various values of zenith angle θsol. It is assumed that the spectrum of the incident solar radiation is similar to the blackbody radiation at temperature Tsol = 6 × 103 K. The calculated relative values of the integral radiation power, P¯t = ∫0∞P(z) dz / ∫λminλmaxIb(Tsol,λ) dλ (Ib is the Planck function), absorbed in a snowpack containing ice grains of different sizes are presented in Fig. 1. It is known that the spectrum of solar radiation is quite different at large zenith angles. Therefore, this range of θsol is not shown in Fig. 1.
Figure 1. Normalized radiative power absorbed in an optically thick snowpack
As one can expect, the main part of solar radiation is reflected from the snowpack because of a strong reflection of the visible radiation. One can see that P¯t decreases considerably with the zenith angle, whereas the effect of ice grain size is insignificant. It should be also recalled that the calculations by Dombrovsky et al. (2019) for polydisperse spherical ice grains showed that one can use the monodisperse approximation with an equivalent mean radius of the grains to obtain the spectral properties of snow.
The depth of the radiation propagation into the snowpack is illustrated in Fig. 2 by calculations at normal incidence. The contributions of the visible range (λ < 0.78 μm) and the near-infrared range to the absorbed radiation power are calculated separately:
The universal coordinates fv z and P¯ / fv are used in Fig. 2 to present the results for various values of fv (the ratio of snow density to ice density). At realistic value of fv = 0.33, the near-infrared radiation is absorbed mainly in a thin surface layer of snowpack; whereas, the visible radiation is absorbed almost uniformly in the layer of thickness ∼60 mm and its contribution to heating deep layers of snowpack is expected to be considerable. This result agrees well with the conclusion by Munneke et al. (2009) about an important role of the deep penetration of short-wave solar radiation in the heat balance of a snowpack. The effect of grain size on this effect is not strong. Therefore, the subsequent calculations are performed at a = 101 μm.
Figure 2. Contributions of the visible range (P¯vis) and near-infrared range (P¯ni) to the profiles of normalized absorbed radiation power in a snowpack
It is interesting to analyze the results obtained at realistic parameters of the problem. The range of snow density variation from ρ = 10 × 101 kg/m3 to 55 × 101 kg/m3 leads to the corresponding variations in the thermal properties of snow. Obviously, the volumetric heat capacity, ρc, is directly proportional to the density, whereas the thermal conductivity of snow exhibits a stronger increase with the snow density. In addition, at the same snow density, variations in snow microstructure can change thermal conductivity by a factor of 2. Thermal conductivity of snow has been studied in many papers. Most of the results obtained are discussed in a review paper by Arenson et al. (2015). It was shown that thermal conductivity is almost directly proportional to the square of density and varies from k = 0.04 W/(m K) to 1 W/(m K). The temperature variations of thermal properties of snow are relatively small, and the effect of temperature on snow density and thermal properties is not considered. The constant values of k = 0.2 W/(m K) and ρc = 0.6 MJ/(m3 K) corresponding to the volume fraction of ice grains of fv = 0.33 are used in the calculations.
In the midsummer (summer solstice), the current declination of the sun is equal to δ = 23.44 deg. This value of δ is used to determine the solar zenith angle
where φ is the local latitude and χ is the hour angle. The value of φ = 70 deg is considered in the calculations. The following relation is used to determine the value of χ for the current time, t, measured in hours from midnight:
where ⎣ ⎦ denotes the floor function. The resulting variation of θsol with the local solar time is shown in Fig. 3. The effect of a sloping surface of snow is also considered. There are many possible orientations of the sloping surface with respect to the visible path of the sun. The simplest orientation can be treated as a variation of the declination angle. Time variation of the solar zenith angle for such a surface with a sloping angle of ϑsl = 10 deg is also shown in Fig. 3.
Figure 3. Time variation of solar zenith angle considered in the case problem: 1 is the horizontal surface and 2 is the sloping surface at ϑsl = 10 deg
The spectral radiative flux from the sun and clear sky to the horizontal earth surface are presented in Fig. 4. The spectral dependencies were calculated using the SBDART code (Ricchiazzi et al., 1998). The choice of solar zenith angles is determined by the case problem parameters. Generally speaking, the contribution of direct solar radiation is much greater than that of diffuse radiation from the sky. However, it is not the case for visible light at solar zenith angles, θsol > 60 deg. Therefore, the sky radiation is taken into account in the calculations.
Figure 4. Spectral radiative flux at the earth surface: (a) direct radiation from the sun and (b) diffuse radiation from the clear sky
According to (Dombrovsky et al., 2019), the following analytical expressions for the everyday variation of coefficients of convective heat transfer and radiative cooling are considered:
It is assumed in Eqs. (4a) and (4b) that the midnight conditions approach gradually starting from θsol = 80 deg. Note that the values of hmax = 10 W/(m2 K) and 15 W/(m2 K) correspond to a weak wind (with a speed of ≲1–2 m/s) and to a moderate wind of speed ∼3 m/s, respectively (Defraeye et al., 2011; Mirsadeghi et al., 2013); whereas, hmin = 6 W/(m2 K) is used for windless weather.
The temperature profiles in a snowpack calculated for the first day at T0 = Tair(z) = –10°C are presented in Fig. 5. One can see that the snowpack is heated rather fast in the morning and the melting temperature is reached at noon. A further evolution of the snowpack temperature depends strongly on wind speed. In the afternoon, the maximum temperature of snow is predicted at the distance ∼15 mm from the snow surface at 6:00 PM and hmax = 15 W/(m2 K). The evening and night cooling of snowpack takes place for the surface layer of thickness from 15 to ∼120 mm; whereas, a continuous heat conduction leads to further heating of deep layers of the snowpack.
Figure 5. Typical profiles of temperature in snowpack (a) before noon and (b) after noon: 1 is hmax = 10 W/(m2 K) and 2 is 15 W/(m2 K)
The computational results for the three-day variation of snow surface temperature presented in Fig. 6 show an almost periodic variation of the temperature with a great difference between day and night. The effects of both the convective heat transfer in windy weather and the radiative cooling of a snowpack surface are significant. It is interesting that radiative cooling is partially compensated by a relatively long convective heating. The latter may lead to the unexpectedly strong surface heating and even melting as compared to the computational predictions ignoring radiative cooling. This physical result of combined transient heat transfer in a snowpack was obtained by Dombrovsky et al. (2019) for the first time. Note that the details of temperature behavior near the melting temperature are explained by the finite value of ΔT = 0.1 K used in the calculations.
Figure 6. Effect of wind on snow surface temperature: 1 is hmax = 8 W/(m2 K) (weak wind) and 2 is 15 W/(m2 K) (moderate wind); (A) complete calculation and (B) without radiative cooling
Similar calculations for the better irradiated sloping surface of snowpack with ϑsl = 10° can be easily calculated by multiplying the directional radiative flux at solar zenith angle θsol(t) by the ratio of cos(θsolsl )/cos(θsol). Note that the conventional night-time conditions, including the radiative cooling, do not take place for the sloping surface considered. The overall effect of snow surface sloping on the time variation of snowpack surface temperature is shown in Fig. 7. Obviously, the snow surface turned to the sun exhibits a higher temperature, especially in the time of the large zenith angle of the sun. This means that specific conditions of solar irradiation cannot be ignored, even at a sloping angle of ϑsl = 10°.
Figure 7. Effect of sloping of snowpack surface on the time variation of surface temperature: 1 is the horizontal surface and 2 is the sloping surface (ϑsl = 10 deg). Calculations at hmax = 15 W/(m2 K).
The results obtained make clear a contribution of the visible and near-infrared solar radiation as well as the role of convective cooling or heating of snow, mid-infrared radiative cooling of snow surface, and continuous heating of deep snow layers due to heat conduction. The numerical analysis indicates a deep heating of snowpack and relatively strong effect of wind and radiative cooling on surface temperature of snow. The calculations for the case problem demonstrated the main possibilities of the computational model, which is expected to be a convenient and useful tool in further studies of snowpack solar heating.
Arenson, L.U., Colgan, W., and Marshall, H.P. (2015) Snow and Ice-Related Hazards, Risk and Disasters, in Physical, Thermal, and Mechanical Properties of Snow, Ice, and Permafrost, Oxford, UK: Elsevier, pp. 35–75.
Defraeye, T., Blocken, B., and Carmeliet, J. (2011) Convective Heat Transfer for Exterior Building Surfaces: Existing Correlations and CFD Modelling, Energy Convers. Manage., 52(1): 512–522.
Dombrovsky, L.A., Kokhanovsky, A.A., and Randrianalisoa, J.H. (2019) On Snowpack Heating by Solar Radiation: A Computational Model, J. Quant. Spectrosc. Radiat. Transf., 227: 72–85.
Mirsadeghi, M., Cóstola, D., Blocken, B., and Hensen, J.L.M. (2013) Review of External Convective Heat Transfer Coefficient Models in Building Energy Simulation Programs: Implementation and Uncertainty, Appl. Therm. Eng., 56(1-2): 134–151.
Munneke, P.K., van den Broeke, M.R., Reijmer, C.H., Helsen, M.M., Boot, W., Schneebeli, M, and Steffen, K. (2009) The Role of Radiation Penetration in the Energy Budget of the Snowpack at Summit, Greenland, Cryosphere, 3(2): 155–165.
Ricchiazzi, P., Yang, S., Gautier, C., and Sowle, D. (1998) SBDART: A Research and Teaching Software for Plane-Parallel Radiative Transfer in the Earth Atmosphere, Bull. Am. Meteorol. Soc., 79(10): 2101–2114.
- Arenson, L.U., Colgan, W., and Marshall, H.P. (2015) Snow and Ice-Related Hazards, Risk and Disasters, in Physical, Thermal, and Mechanical Properties of Snow, Ice, and Permafrost, Oxford, UK: Elsevier, pp. 35–75.
- Defraeye, T., Blocken, B., and Carmeliet, J. (2011) Convective Heat Transfer for Exterior Building Surfaces: Existing Correlations and CFD Modelling, Energy Convers. Manage., 52(1): 512–522.
- Dombrovsky, L.A., Kokhanovsky, A.A., and Randrianalisoa, J.H. (2019) On Snowpack Heating by Solar Radiation: A Computational Model, J. Quant. Spectrosc. Radiat. Transf., 227: 72–85.
- Mirsadeghi, M., Cóstola, D., Blocken, B., and Hensen, J.L.M. (2013) Review of External Convective Heat Transfer Coefficient Models in Building Energy Simulation Programs: Implementation and Uncertainty, Appl. Therm. Eng., 56(1-2): 134–151.
- Munneke, P.K., van den Broeke, M.R., Reijmer, C.H., Helsen, M.M., Boot, W., Schneebeli, M, and Steffen, K. (2009) The Role of Radiation Penetration in the Energy Budget of the Snowpack at Summit, Greenland, Cryosphere, 3(2): 155–165.
- Ricchiazzi, P., Yang, S., Gautier, C., and Sowle, D. (1998) SBDART: A Research and Teaching Software for Plane-Parallel Radiative Transfer in the Earth Atmosphere, Bull. Am. Meteorol. Soc., 79(10): 2101–2114.