Introduction
Snow is a major component of the regional and global energy balance because of its high albedo and often rapid disintegration. Utilization of passive microwave data is suitable for analysis of snow-covered area and snow depth. The influence of snow structure on the microwave brightness temperature T B must be understood in order to measure (and model) snow depth as an input to energy-balance models.
The transformation of fallen snow into a snow-pack with a definite structure is the result of local and regional energy-balance processes. When snow remains on the ground for a substantial part of the winter, metamorphism at the base of the snow-pack can result in the formation of a layer comprised of large, loosely bonded crystals known as depth hoar. Depth hoar is common in snow-packs through-out the world and is well developed in parts of Alaska.
On the Arctic Coastal Plain of Alaska (Fig. 1), low snow accumulation, very cold air temperatures, and strong winds contribute to the snow-pack structure. A relatively thin, continuous snow cover forms each year. Once the snow is deposited, little additional accumulation or melting occurs until spring. High average snow-pack densities result from intense redistribution of the snow by wind, causing the crystals to be rounded and abraded, thus facilitating compaction in the upper part of the snow-pack.
The specific structure of the snow-pack on the Arctic Coastal Plain consists of a high-density (0.35–0.45 g cm−3), wind-packed surface layer overlying a low-density (0.20–0.30 g cm−3) depth-hoar layer with transitional layers in between. In the depth-hoar layer, snow-crystal sizes are large and the crystals are loosely bonded. This layer can comprise 50% of the thickness of the snow-pack on the Arctic Coastal Plain (Reference Benson, Benson, Holmgren, Timmer, Weller, Parrish and BrownBenson and others, 1975).
In this paper, theoretical and observational data indicate that the presence and relative thickness of the depth-hoar layer on the Arctic Coastal Plain of Alaska can be detected by microwave sensors. By observing variations in the responses of microwave sensors and utilizing available meteorological data, one can infer differences in depth-hoar development between coastal and inland sites. Using a two-layered radiative-transfer model, the relationship between grain-size in the depth-hoar layer and microwave response is shown. Also shown is the effect that a thicker depth-hoar layer has on microwave response.
Development of the Depth-Hoar Layer
The difference between the temperature of the soil (at the base of the snow-pack) and the air temperature can be greater than 40 deg in northern Alaska (Reference Slaughter, Crook, Santeford and SmithSlaughter and Crook, 1974). This results in a steep negative temperature gradient in the snow-pack because the temperature at the base of the snow-pack can be −5°C. At the surface of the snow-pack the temperature can be −45°C. Early snow accumulation, such as occurs on the Arctic Slope of Alaska, promotes greater retention of ground heat at the base of the snow-pack (Reference Santeford, Colbeck and RaySanteford, 1979). In addition, shallow snow has the potential for steeper temperature gradients than does deep snow (Reference Benson and TrabantBenson and Trabant, 1973). Snow crystals sublime and the vapor is redeposited on to adjacent crystals in an upward direction in a process called constructive metamorphism (Reference Giddings and LaChapelleGiddings and L.aChapelle, 1962). Faster crystal growth occurs in the warmer (lower) parts of the snow-pack (Reference ColbeckColbeck, 1982[b]).
There are two basic types of snow-pack metamorphism: constructive and destructive metamorphism. In the latter case, snow crystals become rounded and are reduced in size allowing compaction which results in increasing density. An illustration of destructive metamorphism is shown in Figure 2.
Conversely, constructive metamorphism occurs when snow-crystal size increases through time. Crystals grow from one side of existing snow crystals or grains in a direction that is opposite to the vapor-pressure gradient; the same grains lose material by sublimation on the other side. Constructive metamorphism of snow crystals leads to development of depth hoar – increasingly large crystals in the lower layers of the snow-pack. Density per unit volume decreases. Depth-hoar formation is complete when initial snow grains have been completely replaced by the depth-hoar crystals (Reference Giddings and LaChapelleGiddings and LaChapelle, 1962). This is illustrated in Figure 3. Such crystals or grains can be 5–10 mm in diameter (Reference Benson, Benson, Holmgren, Timmer, Weller, Parrish and BrownBenson and others, 1975). The depth-hoar crystals are hollow and cup-shaped.
The temperature gradient produces inequalities in the degree of air saturation and the diffusion rate of water vapor (Reference MarboutyMarbouty, 1980). Associated pressure gradients cause water vapor to diffuse from warmer to colder parts of a snow-pack (Reference Langham, D.M. and MaleLangham, 1981). A local temperature gradient within the snow-pack of 0.1–0.2 deg/cm is sufficient to initiate the formation of depth-hoar crystals (Reference ColbeckColbeck 1982[a]).
The presence of the depth-hoar layer weakens the overall strength of the snow-pack and is often responsible for the failure of snow slabs which can result in avalanches in mountainous areas (Reference MarboutyMarbouty, 1980).
Theoretical Basis for Microwave Observations
Satellite data are useful for analyzing snow conditions on regional and global scales, and microwave sensors are able to acquire data through cloud cover and darkness. Passive microwave sensors which sense the emitted radiation can provide information regarding the internal characteristics of snow and ice (Reference Chang, Chang, Gloersen, Schmugge, Wilheit and ZwallyChang and others, 1976). Previous studies have shown that there is an inverse relationship between snow depth and microwave T B as measured by passive microwave sensors at specified wavelengths in dry snow (Reference Foster, Foster, Hall, Chang and RangoFoster and others, 1984). The 37 GHz (0.81 cm wavelength) sensor of the Scanning Multichannel Microwave Radiometer (SMMR) on board the Nimbus-7 satellite has been shown to be particularly useful for analyzing internal properties of snow-packs especially when the horizontally polarized data are used (Reference Hall, Hall, Foster and ChangHall and others, 1984). Some characteristics of the Nimbus-7 SMMR are shown in Table I.
For temperatures generally encountered on Earth, the emitted intensity of microwave radiation is expressed as brightness temperature (T B) in degrees Kelvin and follows the Rayleigh–Jeans approximation which shows that the radiance from a black body is proportional to its temperature.
where ϵ is the emissivity of the surface, T s is the sensible temperature of the surface, τ is the total atmospheric opacity, T 1 is the upward-emitted radiance contribution of the atmosphere, T 2 is the total downward (emitted and reflected) atmospheric brightness temperature, and T sp is the average temperature of free space (Reference Gloersen and BarathGloersen and Barath 1977).
Reference Chang, Chang, Gloersen, Schmugge, Wilheit and ZwallyChang and others (1976) calculated that the radiation emanates from a depth, in snow and firn, that is 10 to 100 times the length of the wavelength used. Measurement of the amount of sub-surface emission and the depth from which it emanates is dependent on the physical characteristics of the snow or ice as well as on instrument parameters.
Studies of the brightness temperatures of the Greenland and Antarctic ice sheets as determined from the Nimbus-5 Electrically Scanned Microwave Radiometer (ESMR) data revealed large variations in brightness temperature across the ice sheets that did not correlate with physical surface temperature in Greenland and correlated only loosely with surface temperature in Antarctica (Reference Chang, Chang, Gloersen, Schmugge, Wilheit and ZwallyChang and others, 1976). Observation and modeling have shown that differences in emissivity are related to variations in snow accumulation rates, mean annual temperature, and melting effects on the ice sheets (Reference Zwally and GloersenZwally and Gloersen, 1977). Specifically, larger grain-sizes in ice-sheet locations having low accumulation rates are known to cause more scattering and to be associated with lower microwave emissivities.
Calculations have shown that the grain or crystal size is a dominant factor influencing the microwave emission of snow and dry polar firn. Additional calculations using a microscopic scattering model for snow have shown that the scattering of 37 GHz radiation from a snow-pack is strongly dependent upon the grain-size of the snow particles (Reference Chang, Chang, Foster, Hall, Rango and HartlineChang and others, 1982).
As snow grain-size increases, microwave T B of a given snow-pack should decrease if other parameters remain the same (Reference Chang, Chang, Foster, Hall, Rango and HartlineChang and others, 1982). Increasing grain-size within the snow-pack allows greater incidence of scattering of microwave radiation as the grain-size approaches or surpasses the size of the wavelength. The large crystals characterizing the depth-hoar layer would tend to increase scattering of the upwelling radiation relative to scattering through a snow-pack (with a similar thickness) which is lacking a depth-hoar layer.
Methodology
For this study, Nimbus-7 SMMR data of northern Alaska were acquired and analyzed for four winters: 1979–80, 1980–81, 1981–82, and 1982–83. Data were analyzed from 1 October through 31 May. Data from all available night-time SMMR overpasses were used at the 37 GHz wavelength, horizontal polarization. Each brightness temperature value is an average of six night-time overpasses. Three areas within northern Alaska were analyzed and these are shown in Figure 1. Each area comprises 15 SMMR pixels. The resolution of each pixel in this study is 1/2 degree longitude × 1/2 degree latitude. The 15 SMMR pixels were averaged to obtain a mean brightness temperature value for each night-time satellite overpass within areas A, B, and C.
Snow-depth and air-temperature data of Barrow, Alaska, were used for area A, the northernmost area shown in Figure 1. The meteorological station located at Umiat, Alaska, was within the southernmost area, or area C, and was considered representative of area C. No station was used to represent meteorological and snow-depth conditions in area B.
Observational Results
Figures 4, 5, 6, and 7 show graphs of 37 GHz horizontally polarized T B for three study areas on the Arctic Coastal Plain of Alaska for four winters: 1979–80, 1980–81, 1981–82, and 1982–83. Note that in each winter, T Bs decrease in the latter half of the winter especially in areas B and C. Once snow is on the ground, areas A and B have consistently higher T Bs than area C. Average air temperatures at Umiat were somewhat lower than at Barrow (Table II), and the temperture extremes are greater at Umiat than at Barrow because of the more continental location of the Umiat station. This can be seen in Table II by noting that the standard deviations are higher for the Umiat data than for the Barrow data (NOAA, 1979–83).
In the present study area, the lower T B in areas B and C is apparently not the result of deeper snow. Though previous work has shown that under certain conditions there is an inverse relationship between snow depth and microwave brightness temperature (Reference Rango, Rango, Chang and FosterRango and others, 1979), exceptions to this are many (Reference Foster, Foster, Hall, Chang and RangoFoster and others, 1984). Table III shows that snow depths were generally greater at Umiat than at Barrow, but not always. For example, the average February through April 1980 snow depth was 35.80 cm at Barrow and 30.70 cm at Umiat. But the T Bs shown in Figure 4 were consistently lower in area C where the Umiat station is located.
A thicker depth-hoar layer in the more interior areas (areas B and C) would explain the lower T Bs as it is known that depth-hoar development is more pronounced in inland areas (Reference Slaughter, Crook, Santeford and SmithSlaughter and Crook, 1974). In addition, the crystal sizes in the depth-hoar layer may become larger as the winter progresses (Reference Giddings and LaChapelleGiddings and LaChapelle, 1962); the separation in T B (between areas A and B, and B and C) becomes greater as the winter progresses, indicating an increase in (1) size of depth-hoar crystals, and/or (2) thickness of the depth-hoar layer.
The differences in T B (among areas A, B, and C) do not show up until snow was on the ground as seen in Figures 4 through 7, Thus, the change in T B is clearly related to the presence of snow cover.
Colder air temperatures are conducive to greater depth-hoar development within a snow-pack because the potential temperature gradient within the snow-pack is higher, causing the vapor-pressure gradient to be greater as well. Vapor-pressure gradients within the snow-pack are larger in more continental climates (Reference Slaughter, Crook, Santeford and SmithSlaughter and Crook, 1974). Thus, one would expect depth-hoar development to be greater in area C than in area A because of colder minimum air temperatures in area C. As mentioned earlier, larger grain-sizes (characteristic of depth hoar) are concomitant with increased scattering and thus lower brightness temperatures and more radiational scattering.
The presence of sea ice near the coast of northern Alaska during the winter changes the maritime climate of the coastal areas to a climate that is more similar to an inland or continental climate. The effect of continentality is evident only a few kilometers inland from the Arctic Ocean and it rapidly becomes more pronounced farther inland; Reference Brown, Brown, Haugen, Parrish and BrownBrown and others (1975) referred to the winter coastal stations as having “near-continental” temperature patterns.
Sea ice slows the release of latent heat from the ocean to the atmosphere. The presence of leads (linear features of open water) within the sea ice permits increases in heat flux while open because the heat flux through open water can be 100 times greater than through thick ice (Reference WeeksWeeks, 1981). In addition, the pack ice has a lower reflectivity than fresh snow in the visible and near-infrared parts of the spectrum than the adjacent snow-covered tundra because of ice-surface roughness, pressure ridges, and hummocks (Reference Holmgren, Holmgren, Benson, Weller, Weller and BowlingHolmgren and others, 1974). Thus, the absorption and emission characteristics of the sea ice are different from those of the land surface particularly during the daylight hours in the spring and fall. The winter climate of the coastal part of the Arctic Slope is neither truly maritime nor truly continental.
No direct measurements were obtained during this study period to confirm the theory that depth hoar is better developed in area C as compared to areas A and B. However, model results (see next section) show that the difference in brightness temperature between coastal and inland sites can be accounted for by increasing the thickness of the depth-hoar layer in more inland sites.
Model Results
The intensity of microwave radiation emitted from a snow-pack depends on the physical temperature, grain-size, density, and the underlying surface conditions of the snow-pack. By knowing these parameters, the radiation emerging from a snow-pack can be derived by solving the radiative-transfer equation. The radiative-transfer equation for an axially symmetric inhomogeneous medium can be written in the form of an integro-differential equation
where the radiation intensity I(x,µ) is at a depth x and traveling in the direction towards increasing x, making an angle whose cosine is µ with the normal (Reference Chang, Chang, Foster, Hall, Rango and HartlineChang and others, 1982). The functions σ(x), ψ(x), B(x), and P(x,µ,µ′) are prescribed functions of their arguments. They are the extinction per unit length, the single scattering albedo, and the source and phase functions, respectively.
Equation (2) is solved numerically using an invariant imbedding technique (Reference Chang and ChoudhuryChang and Choudhury, 1978). This technique is based on the principle that the radiation emerging from a semi-infinite, plane-parallel medium is invariant with respect to addition (or subtraction) of layers of arbitrary thickness to (or from) the medium. This technique enables the radiative-transfer equation to be solved numerically using recurrence relationships.
By this method, the emerging brightness temperatures are calculated for snow-packs with different parameters. Assumptions are that the total snow depth is 30 cm, the underlying soil is frozen, and the snow-pack is dry. The snow grains scatter the electromagnetic radiation incoherently and are assumed to be spherical in shape and randomly spaced within the snow-pack. This two-layer model allows for assignment of different top and bottom layer thicknesses, grain-sizes, and densities. Further discussions on the assumptions of the microscopic model can be found in Reference Chang, Chang, Gloersen, Schmugge, Wilheit and ZwallyChang and others (1976).
Because the grains in the depth-hoar layer are hollow and cup-shaped and not spherical, they are not as effective at scattering upwelling microwave radiation as are spherical grains. Thus, when considering the microwave scattering, the effective diameter of the grains is considerably smaller than the actual diameter of the largest dimension of a depth-hoar crystal. Using the model, employment of grain-sizes in the lower layer comparable to those typically observed in nature yielded brightness temperatures that were far too low as compared to observed brightness temperatures.
An effective snow-grain diameter of 1.40 mm for the bottom or depth-hoar layer was used in the model. Use of this diameter rendered T Bs that were comparable to those observed by the SMMR in the three study areas. The actual and effective grain-size diameter should be approximately equal in the top layer. This was assumed to be 0.50 mm.
In addition to the fact that the effective scattering size is smaller because of the configuration of the depth-hoar crystal, the average grain-size in the depth-hoar layer is much lower than the 5–10 mm maximum size reported by Reference Benson, Benson, Holmgren, Timmer, Weller, Parrish and BrownBenson and others (1975).
Results of the modeling are shown in Figure 8. As expected, brightness temperatures decrease with increasing thickness of the bottom or depth-hoar layer. Calculated brightness temperatures shown in Figure 8 are comparable to observed brightness temperatures in the study areas in mid to late winter when the snow-packs are mature. Model results show a sharp (33K) lowering of the T B when even a thin (5 cm) depth-hoar layer develops. As the depth-hoar layer increases in thickness, the rate of change in brightness temperature decreases as seen in Figure 8.
Calculated T Bs compare favorably with observed values for mature snow-packs in three of the four winters studied. (In the 1981–82 winter T Bs were approximately 20K lower than in the other three years.) The actual thickness of the depth-hoar layer was not measured. However, if one assumes that the thickness of the depth-hoar layer was 5 cm in area A and 10 cm in area C, then the model-generated T Bs fall within the range (the high end of the range) of observed T Bs. The calculated T B for a snow-pack with a 5 cm depth-hoar layer is 201K and the range of observed T Bs in area A is 181–201K (excluding the 1981–82 winter). For a snow-pack with a 10 cm depth-hoar layer, the calculated T B is 180K and the range of observed brightness temperatures in area C is 165–181K.
The magnitude of decrease in calculated T Bs between depth-hoar layer thicknesses of 5 and 10 cm is 21K. This approximates the difference in T B observed using SMMR data between areas A and C where T B differences are approximately 20K.
Conclusion
The fact that 37 GHz passive microwave radiation scattering increases with snow grain-size has been discussed in terms of the accumulation rate of polar firn on the Greenland and Antarctic ice sheets. Observations and calculations have shown this to be a dominant factor in the interpretation of passive microwave data of seasonal snow-packs as well.
Results of this study indicate that the depth-hoar layer, with its very large snow grain-sizes, can be detected using Nimbus-7 Scanning Multichannel Microwave Radiometer (SMMR) 37 GHz passive microwave data and that geographic variations in the thickness of the depth-hoar layer may also be observable.
Employment of a radiative-transfer model shows the effect of varying the thickness of the depth-hoar layer on the calculated brightness temperatures. A calculated decrease in T B of 21K caused by increasing the thickness of the lower layer from 5 to 10 cm approximates the observed changes in T B (between coastal and inland sites) seen on the SMMR data of the Arctic Coastal Plain of Alaska. Because the depth-hoar crystals are hollow and cup-shaped, they are not as effective at scattering the upwelling radiation as are spherical grains. The effective scattering size of the depth-hoar crystals is considerably lower than the diameter of the largest dimension of the crystal as measured in the snow-pack.
Acknowledgement
The authors would like to thank Drs R. Gurney and T. Schmugge for their reviews of this paper.