Introduction
Winter in the North American Great Lakes is of sufficient length and severity that a partial or total ice cover is produced. Each Great Lake possesses its own set of characteristics both physical (volume, surface area, shore length, etc.) and hydrometeorological (temperature, precipitation, evaporation patterns, etc.). During a normal winter the aerial ice extent varies from 15% on Lake Ontario to 95% on Lake Erie (Reference RondyRondy, 1976). During a severe winter, the aerial ice extent on all of the Great Lakes can approach 100%. Each of the ice types common to the Great Lakes has a unique appearance, mode of formation, and physical characteristics. The physical characteristics (strengths, reflectances, crystallographies, etc.) are studied to provide information useful for solving problems in winter navigation, shore-line engineering, hydropower generation, water supply, and water quality. One of those physical properties, the spectral reflectances of the various ice types, is examined here.
A pilot study was conducted in early 1967 to obtain the total (wavelength-integrated) albedo of various types of ice common to the Great Lakes (i.e. clear, slush, ball, pancake, etc.) at only one solar altitude and cloud condition. Measurements ranged from 10% for clear ice to 46% for snow-ice (Reference BolsengaBolsenga, 1969). The program indicated that many gaps existed in our knowledge of the albedo of fresh-water ice types. A new program was subsequently initiated to develop a catalog of information that would provide measurements of the total albedo of fresh-water ice within a wide range of solar altitudes (10–45°), cloud types, and cloud amounts. A measurement system was fabricated consisting of upward- and downward-facing pyranometers mounted on a tripod and boom assembly. Reports from that study (Reference BolsengaBolsenga, 1977, Reference Bolsenga1979, Reference Bolsenga1980) commented on the diurnal variation of, and on the influence of changing cloud conditions on, the total albedo of ice.
The need for ice and snow reflectance measurements over various spectral ranges was emphasized by these and other earlier studies. With data at reasonably small wavelength intervals, reflectances could be determined corresponding to the current or future sensitivity of detectors designed for many applications. In order to satisfy this need, two scanning spectroradiometers manufactured by the EG&G Corporation, were used to measure radiation in the 280–1 100 nm rangeFootnote * (10 nm bandpass slits). Two instruments were required to measure incident and reflected radiation at each wavelength simultaneously (one instrument directed to the zenith and the other to the nadir, both synchronized at equal scanning speeds). As elaborated elsewhere (Reference Bolsenga and KistlerBolsenga and Kistler, 1982), optimum instrument configuration requires that incoming radiative flux be recorded simultaneously with reflected flux since reflectances calculated from non-simultaneous spectroradiometer scans can be subject to significant errors due to possible changes in the atmosphere or solar altitude during the period of the scan. At the measurement site the spectroradiometers were mounted on a boom 4.6 m long supported at each end by tripods, placing the reflected unit about 1.5 m above the ice surface.
Few studies have been conducted on the spectral reflectances of ice. One of the most recent (Reference Grenfell and MaykutGrenfell and Maykut, 1977) deals with reflectance measurements of sea ice and refrozen melt ponds. In addition to the fact that sea ice rather than fresh-water ice was studied, Grenfell and Maykut used only one “spectrophotometer” requiring time to complete an incident and subsequently a reflected scan. It is emphasized that the instruments or data in this study are not adjusted to the human eye response (CIE curve). Reference Perovich and GrenfellPerovich and Grenfell (1981) report on the spectral reflectance of laboratory-grown young sea ice.
Reference DugginDuggin (1980) independently developed a technique similar to that used here but with different equipment. Both studies were apparently conducted concurrently and arrived at similar conclusions with respect to the technique involved. Duggin calibrated in the field using a BaSO4 coated plate, whereas in this program both instruments were pointed to the zenith using the Sun plus sky as a calibration source. The instrumentation used in this study, including wiring to control the two spectroradiometers as well as software to process the data in the field immediately, is described in detail in Reference Bolsenga and KistlerBolsenga and Kistler (1982).
In an early, but classic study, Reference KrinovKrinov (1947) determined the spectral reflectances of a great variety of natural surfaces, including a limited number of measurements of snow crusted with ice. Krinov recognized the problems associated with changing atmospheric conditions and solar altitude:
“As stated above, the illumination of the object and the standard surface must be the same. Thus the spectrograms were usually obtained when the sky was clear and seldom when the sky was clouded and only when the clouds were low (cumulus, strato-cumulus, broken cumulus, etc.) and the total amount of cloud in the sky did not exceed 0.3; in this case exposures were made only when the sun broke through the clouds. In a few cases spectrograms were obtained when there was a continuous even cloud cover extending across the entire sky. Usually the spectrograms were taken near midday beginning two hours before noon and ending two hours after noon.”
In order to achieve practical results comparable with those reported in the literature, the standard 8° field-of-view of the spectroradiometers was modified by the use of hemispherical cosine diffusers. Initial cosine response tests at 550 nm showed agreement with ideal cosine response within a few per cent depending on the angle of the diffuser from a calibration light source. Further tests showed that diffuser cosine response differed increasingly from true cosine response with increasing wavelength.
Due to these difficulties, all of the measurements reported here were collected with the plastic diffusers but under overcast skies unless specifically noted. Observers classified skies as heavily overcast only when entirely sure of the lack of direct radiation as indicated by the unaided eye. In some cases described later, matte paper, which produced much improved cosine response, was taped over the hemispherical diffusers.
Measurements of snow and ice reflectances were acquired at several sites including the Straits of Mackinac (lat. 45° 42′ N., long. 84° 30′ W.) located between Lakes Michigan and Huron (Fig. 1) and North Lake, a small inland lake located near Ann Arbor, Michigan (lat. 42° 25′ N., long. 85° 57′ W.).
Ice at the Straits of Mackinac sites was relatively young at the time of the measurements. First ice of the season, as detected by side-looking airborne radar and reported by the U.S. National Weather Service, was located in a near-shore area south of Bois Blanc Island on 18 January 1980. Spectral reflectance measurements spanned the period 19–22 February 1980. On 25 January 1980, thin ice (10–20 cm) was reported to cover fully the area of the measurements from Bois Blanc Island to the mainland on the south-west. On 28 January 1980, satellite imagery and U.S. Coast Guard reports verified that condition. On 30 January 1980 the ice was 23 cm thick at one location in the area of the spectral reflectance measurements. The thickness increased to 38 cm by 13 February 1980. The ice continued to cover the area fully from 25 January through 2 April 1980.
Ice growth on North Lake proceeded as one might expect for a small inland lake unaffected by high winds or significant water currents. Snow-ice formed first due to snow-fall into near-freezing water and subsequent congelation. A warm spell caused break-up of the thin ice layer. The “white-ice” fragments were subsequently frozen with new clear ice in certain near-shore areas. A moderate amount of snow accumulated on that ice and subsequent mild temperatures melted the snow to a slush layer which refroze later due to lower temperatures. Snow loads were never sufficient to fracture the ice and cause water to percolate upward into the snow layer. Melting and refreezing of overlying snow and additional growth of clear ice did, however, continue throughout the season.
Analysis
A descriptive classification of Great Lakes ice types was made prior to this study by Reference MarshallMarshall (1966, Reference Marshall1977 [a], Reference Marshall[b]). Spectral reflectance measurement sites were selected on the basis of his classifications. Measurements were made on snow, clear ice, refrozen slush ice, brash ice, pancake ice, and slush curd ice.
An extensive amount of research on the reflectance of snow has been conducted with pyranometers, thus confining the data to integrated values in the 300–2 000 nm range. Experimental work on the spectral reflectances of snow in small wavelength increments by Reference O’Brien and MunisO’Brien and Munis (1975), Reference O’BrienO’Brien (1977), and Reference GrenfellGrenfell and others (1981) seems to agree reasonably well with these results although O’Brien and Munis’ data, it should be noted, were collected under artificial conditions and with the source beam and detector at various angles to the surface. Reference Warren and WiscombeWarren and Wiscombe (1980), Reference Wiscombe and WarrenWiscombe and Warren (1980), Reference ChoudhuryChoudhury (1981), and Reference WarrenWarren (1982) have addressed the question of the modeling of the spectral reflectance of snow.
Several scans were made over 3–4 cm of new snow overlying 15 cm of mostly clear lake ice. Skies were totally overcast with sporadic light snow flurries. Incoming radiation for all runs was totally diffuse as observed visually. Figure 2 shows results from two runs during a period without snow flurries. Spurious data points in this and subsequent figures are due to instrumentation noise. It is emphasized that these points, evidenced by sharp spikes in the graph, are not considered to represent real data. The curves are characterized by a narrow range of values (81–89%) from 340–900 nm after which a pronounced dip to values of about 70% occurs near 1 050 nm.
After these measurements the surface was liberally trampled by the observers. Care was taken to expose only a minimum of the underlying ice surface. The measurements over the trampled surface compared to those over undisturbed snow show small (2–3%), consistently lower spectral reflectances of the trampled surface from 340–900 nm. From 900–1 100 nm little differences in spectral reflectance was noted. Contrary to the very small differences in measured spectral reflectances, the trampled surface appeared quite different to the eye from the undisturbed surface. Perhaps the “look angle” of the spectroradiometers as opposed to the “look angle” of the eye of the observer accounted for the discrepancy. Measurements were made again on snow at the same location three days later. The surface had obviously metamorphosed during the intervening period, being hard packed and etched by the wind. The shape and numerical values of the curves are very similar (differences, generally <5%) to those in Figure 2. Scans over other hard-packed snow surfaces showed slightly lower spectral reflectances than those for the new snow shown in Figure 2.
One run was made over a snow surface that had reached an advanced stage of deterioration. The mass was dark in appearance and nearly at the stage of complete melting. Sky cover was totally overcast but thin in spots. Weak beam radiation was observed sporadically through thin portions of the overcast, but only during about 10% of the period for the run. The shape of the curve is similar to that of fresher snow but with greatly reduced values (Fig. 2), more appropriately ascribed to poorly reflecting refrozen slush.
Depending on bubble content, fracture patterns, and thicknesses, clear ice can demonstrate a remarkable ability to transmit solar radiation at certain wavelengths. In a study of the transmittance of photosynthetically active (400–700 nm) radiation through ice, Reference BolsengaBolsenga (1981[a]) found that nearly 90% of the incoming radiation was transmitted through 28 cm of clear ice. (The transmittance dropped to 10% for the same ice surface with a 3 cm thick layer of snow.) This property is important in considering the spectral reflectances of ice, since the underlying water or the bottom of the lake is often clearly visible through clear ice and can obviously influence reflectances in certain spectral ranges.
Three dual spectroradiometer runs over approximately 30 cm of clear ice were made at the Straits of Mackinac location (Fig. 3). The clear ice probably formed in a zone of open water between large areas of ice of different types. Sky cover was dense overcast (radiation totally diffuse). Neither the Sun nor shadows were visible throughout the period of the runs. The ice surface appeared moist when the 5 m diameter area was cleared of snow for the measurements. Temperatures were near freezing and snow at the site compacted easily. It is postulated that mild weather caused water to percolate through the snow layer to the snow–ice interface. The surface appeared very black, no doubt partly due to the 22 m water depth at the site. The depth of the water precluded any influence of the bottom on scan values. The curves lack any abrupt changes in spectral reflectance with recorded reflectances remaining in a narrow (5–10%) range for three spectral scans. A general decrease in reflectance throughout the spectrum is evident with a slight, almost imperceptible rise in values from 340 to roughly 500 nm, where the highest values occur.
In another case, the spectral reflectances of 15 cm of clear ice on North Lake with a water depth of only 1 m was measured (Fig. 3). Although the surface was predominantly clear ice, several fragments of “white ice” (snow-ice or refrozen slush) from a previously melted ice cover were frozen in the clear ice matrix and numerous cracks were noted. Higher spectral reflectances from 340–850 nm for this ice surface as compared with the Straits of Mackinac surface can be ascribed to cracks and the “white-ice” inclusions.
Refrozen slush ice is one of the most common ice surfaces on the Great Lakes and on smaller inland lakes. High snow-falls with intervening high temperatures or rain often reduce an overlying snow cover to slush. Subsequent refreezing often occurs due to lowered temperatures. Three scans with matte paper over the plastic diffusers and one scan with plastic diffusers alone were made over a refrozen slush-ice surface (Fig. 3) at the Straits of Mackinac about 150 m offshore. Since the ice was located in the near-shore zone, formation of the surface probably began in mid-January as noted earlier in this paper. A thick overcast was present and no shadows were observed during the measurements. Since data from the three matte-paper runs were nearly equal, data from only one such run are shown in Figure 3. Data from the run with plastic diffusers are nearly equal to spectral reflectances for the run with paper diffusers from 340–750 nm. Even though skies were totally overcast, above 750 nm a disparity of 2–4% is evident. Differences between the two runs at higher wavelengths are most likely to be due to uneven radiance (intrinsic radiant intensity emitted by a radiator in a given direction) distribution over the sky hemisphere (Reference PetherbridgePetherbridge, 1955; Reference Kasten and MöllerKasten and Möller, 1960; Reference Gordon and ChurchGordon and Church, 1966) combined with the fact that (1) refrozen slush is not a diffuse reflector, and (2) cosine response errors of the plastic diffuser are severe at higher wavelengths. The errors are not considered to be of sufficient magnitude seriously to hamper the usefulness of the data for most applications. The type of ice surface and the associated likelihood of specular reflection influence the possibility of this type of error. Specular reflection does not occur with Lambertian surfaces, and without specular reflection, spectral reflectance errors attributed to uneven sky radiance have not been detected.
Two runs taken 30 min apart over another refrozen slush-ice surface on North Lake (Fig. 4) show a significant difference of about 10% from the spectral reflectances of refrozen slush ice at the Straits of Mackinac from 340–700 nm. Differences in spectral reflectances between these refrozen slush-ice surfaces can be attributed to variations in surface texture and bubble content.
Since pancake ice is composed of fragments of a previously formed ice cover, the individual pancakes consist of any type of ice (clear ice, snow-ice, etc.) with raised rims of congealed small ice fragments, ranging in size from sand to gravel, caused by grinding of one fragment against another. The result of this abrasive process appears as a “whitish” rim enclosing an inner portion varying in appearance from white (snow-ice, slush ice) to black (clear ice). Spectral reflectances of such ice could be expected to vary widely, depending on the make-up of the inner portion of the individual pancakes and the matrix in which the individual pancakes are frozen.
Only one type of pancake ice was measured in this study (Straits of Mackinac, approximately 1.5 km offshore) in which the individual pancakes were formed from refrozen slush-ice or snow-ice fragments. The pancakes were not frozen together horizontally, but in a mildly chaotic fashion. Since the ice was located a fair distance from shore, it is speculated that a weak ice cover formed, only to be fragmented by winds and water currents. The fragments subsequently resided in a zone of mild, but prolonged turbulence where continued abrasion of the edge of one piece against another caused the characteristic raised rims and circular shape. Higher winds then thrust one pancake partially over another in the “mildly chaotic fashion”. A significant amount of snow was swept from the measurement site with a broom, including the areas under and between the slightly upthrusted pancakes. Three runs were made over the surface under a heavy overcast. Figure 4 shows spectral reflectances from one run lacking spurious data points. Spectral reflectances are generally about 50% from 340 to 700 nm. After 700 nm the spectral reflectances gradually decrease to nearly 20% at about 1 050 nm. The dip and subsequent rise in spectral reflectances from 950–1 100 nm is similar to that found in snow.
The spectral reflectances shown in Figure 4 are not likely to be characteristic of the spectral reflectances of all types of pancake ice. “White” ice (refrozen slush or snow-ice) formed most of the surface. Darker ice such as clear ice could as easily have been involved. The spectral reflectances of the pancake ice measured here can be perceived as similar to the spectral reflectances of refrozen slush ice. The spectral reflectances of the other form of pancake ice described can be predicted to be a combination of the spectral reflectances of refrozen sluch ice (rims only) and clear ice. Similarities between the measured spectral reflectances for pancake ice and those shown for the refrozen slush ice in Figure 3 are noteworthy. Using pyranometers, Reference BolsengaBolsenga (1969) measured the albedo of snow-free pancake ice composed of “white” ice pancakes in a clear ice matrix as 31% under 4/10 altocumulus at a solar altitude of 32°. The average spectral reflectance of the data shown in Figure 4 is 40%. Clear ice comprised a greater proportion of the total area for the ice in the earlier study, which accounts for the lower albedo value.
Two distinct varieties of slush curd ice were measured at the Straits of Mackinac. Slush curd ice forms when a thick slush layer freezes into a coagulated or thickened pattern with intervening clear-ice areas. The different appearance of the different types of slush ice was probably caused by differing slush contents of the water combined with varying degrees of water turbulence. At the first site the slush curds were frozen in a matrix of clear ice. The visual appearance of the ice was highly variable due to a variety in size and concentration of the individual curds. The ice was 30 cm thick over water 24 m deep. Multiple runs were made at two points to observe an anticipated difference in spectral reflectances over portions of the ice which appeared different. A medium-thick to thick overcast was present during the runs.
As shown in Figure 5 a small, but easily distinguishable, difference of 2–3% from 340 to 1 000 nm was measured between runs over the “darker” and “lighter” surfaces. The similarity of the slope of the curves to that of clear ice curves at the Straits of Mackinac on the same day (Fig. 3) is striking. Spectral reflectances of the slush curd ice are slightly higher than those of the clear ice in all cases, but with the large number of whitish-appearing slush curds, it is surprising that the spectral reflectances were not much higher.
At the second site the surface of the slush curd ice appeared to the casual observer as snow-ice or refrozen slush. Closer scrutiny revealed definite slush curds about 10–30 cm in diameter and closely packed (no clear ice). Figure 5 shows the spectral reflectances under a thick overcast using plastic diffusers and reflectances from the same ice surface with matte paper over the plastic diffusers. The differences between the two runs at longer wavelengths is explained by sky radiance differences interacting with the type of ice and with diffuser transmittances.
Brash ice could be described as one of the most visually distinguishable types of ice. It is perhaps a paradox that the same ice type has little chance of being ascribed a “characteristic” set of spectral reflectance values. The individual blocks of ice might have been composed originally of any type of ice ranging from clear ice to snow ice with spectral reflectances varying accordingly. Final arrangement of the ice fragments after congelation is usually chaotic, but preferred orientation of the fragments is possible. A wide variation of spectral reflectances would thus occur for most brash-ice sites on clear or partly cloudy days due to changing solar altitudes over the period of a day interacting with ice-block orientation. Spectral reflectances of brash ice are thus not only site specific, but also diurnally dependent.
The brash ice measured here was located at the Straits of Mackinac nearly 5 km offshore. Individual ice fragments were composed of clear ice 4–5 cm thick. Light was observed to penetrate the blocks easily imparting a “bluish” color to some areas of the mass. Skies were heavily overcast during all runs. The first two runs were completed with a considerable amount of snow both covering and deposited between the individual ice fragments. After those runs, much snow was removed from the measurement area by hand and with a broom. The ice with an undisturbed snow cover indicated high (>80%) spectral reflectances over most of the visual range (Fig. 5). Reflectances lowered fairly smoothly after 700 nm. Near 1 000 nm a dip and rise was found similar to that observed in snow. For the scan after the snow had been partially removed, reflectances were lower than for the undisturbed surface at all wavelengths with the disparity increasing with increasing wavelengths.
Summary
Results of measurements with a state-of-the-art, dual, scanning spectroradiometer system have provided spectral reflectance data under over-cast skies for several types of fresh water ice commonly found in the Great Lakes as well as many smaller bodies of water. It is hoped that the scans will provide information useful for remote sensing, energy-budget analysis, winter aquatic biology, and for spectral reflectance modeling of snow and ice. Fuller details of this work are available in Reference BolsengaBolsenga (1981[b]).
Acknowledgements
D. C. Norton (Great Lakes Environmental Research Laboratory) assisted in the field work and development of the apparatus for collecting field data. S. I. Outcalt and G. H. Suits (University of Michigan) provided invaluable technical advice. R. D. Kistler (Great Lakes Environmental Research Laboratory) and R. Watson (EG&G Inc.) assisted in development of the data acquisition and processing software.