The McMurdo Dry Valleys (MDVs) of Antarctica are a polar desert
ecosystem consisting of alpine glaciers, ice-covered lakes, streams, and
expanses of vegetation-free rocky soil. Because average summer temperatures
are close to 0
In most regions of the cryosphere, the absorption of shortwave radiation is the main source of energy for snow and ice ablation (Male and Granger, 1981). Therefore, albedo, the proportion of incoming radiation that is reflected, is an important parameter to accurately measure. Albedo for fresh snow can be very high but decreases as snow melts; ice grains age and metamorphose; and dust, debris, or water is accumulated on the surface (Grenfell, 2004). Over the course of a season, there can be very rapid and spatially variable changes in albedo, which feeds back on the degree of meltwater generation (e.g., Perovich and Polashenski, 2012).
During the austral summer in the McMurdo Dry Valleys (MDVs) of Antarctica, the constant solar radiation and near freezing average temperatures mean that any small change in solar radiation or albedo can cause large changes in glacial melt (Fountain et al., 1999). Glacial melt is the main source of water to ephemeral streams, the soil zone near the streams, and closed-basin lakes. For this reason, it is important to correctly predict the degree of melt in order to understand ecosystem response and function. Beyond meltwater generation, the ecological processes across the entire landscape are regulated by incoming solar radiation, spatial distribution of snow and liquid water, and soil and sediment temperatures. For example, the light reflected by the permanently ice-covered lakes has a large effect on the primary productivity of organisms in the lake water column below (Fritsen and Priscu, 1999). Over 24 years of monitoring, Gooseff et al. (2017) observed changes to both physical and biological properties across the MDV ecosystem, especially in response to the 2001–2002 “flood year” in which the highest glacial melt on record occurred.
Extensive modeling of ablation on the MDV glaciers has suggested that there
may be long-term changes in albedo
(Hoffman et al., 2014, 2016). Hoffman
et al. (2016) developed a spatially distributed glacier energy balance
model, which consistently underpredicted ablation in the years after 2006.
They found that a
Albedo is typically measured directly, via remote sensing, or modeled numerically. Each of these approaches are well accepted but have unique shortcomings. Local albedo measurements are made at a point, making it difficult to generate a high enough density of measurements to create an accurate spatial coverage of a parameter that is highly spatially variable (Grenfell, 2004). Remote sensing is limited by spatial resolution, the number of cloud-free days over which the measurement can be made, and may have errors due to the atmospheric absorption (Stroeve et al., 1997, 2005). Numerical models are typically based on experimentally derived parameters and may be poorly linked to the true physical properties of the surface being modeled (Gardner and Sharp, 2010).
The glaciers of the MDVs are generally small, only a few kilometers wide,
and due to their high latitude (76
The McMurdo Dry Valleys are the largest ice-free area on the Antarctic continent (Levy, 2013). They are bounded by the Transantarctic Mountains to the west and the Ross Sea to the east. The landscape is composed of piedmont and alpine glaciers, permanently ice-covered lakes on the valley floors, and vegetation-free dry rocky soils over permafrost, with localized areas of buried ice.
During the winter (March–September), there is little to no incoming solar
radiation and temperatures average
This study focuses on Taylor Valley. The head of the valley is filled by Taylor glacier, which flows from the East Antarctic ice sheet (Fig. 1). The valley is bounded by the Kurki Hills to the south and the Asgard Range to the north. Cold-based alpine glaciers flow out of both ranges toward the valley floor (Chinn, 1998). The largest two glaciers, Canada and Commonwealth, originate in the Asgard Range and reach the valley floor, terminating in 20 m tall cliffs (Fig. 1). Glaciers, particularly those that reach lower elevations, receive wind-transported material that melts the ice below it and sinks in to form cryoconite holes (1–100 cm diameter) or larger melt basins (5–20 m diameter). There is a distinct visible difference in glacier morphology on the up-valley side of the Canada and Commonwealth glaciers that are exposed to foehn winds and therefore have increased potential for sediment deposition (Fig. 1, Sabacka et al., 2012).
Taylor Valley, Antarctica. This study focuses on the three labeled glaciers (Taylor, Canada, and Commonwealth) and lakes (Bonney, Hoare, and Fryxell). Meteorological stations measuring incoming and reflected shortwave radiation are shown as points. An example flight path, flown 12 January 2016, is shown by small green points. Background imagery is provided by USGS and NASA; Landsat8 image collected on 14 January 2016.
There are three main closed-basin permanently ice-covered lakes: Lake Bonney, Lake Hoare, and Lake Fryxell (Fig. 1). The lakes also receive sediment from aeolian deposition and, due to coupled high melt and sublimation rates, have highly variable roughness and ice morphology. Moats frequently melt out around the lake edges during the austral summer that refreeze to smooth ice in the fall. The permanent ice can develop melt ponds and other ablation features that can be several meters in area and reach over a meter of vertical relief. Over the course of work on the MDVs, researchers have anecdotally observed that lakes follow cycles of sediment accumulation and ablation feature development. Lakes receive water from more than 20 ephemeral proglacial streams and direct melt from glaciers. The rest of the valley is made up of bare, rocky soil and exposed bedrock.
This study takes a helicopter-based remote sensing approach to measure landscape albedo. We manufactured a weighted box with a fin to reduce swing during a flight (Fig. 2a). Underneath the box, we mounted a shortwave radiometer (CMP3, Kipp and Zonen) and camera (CC5MPX, Campbell Scientific). A GPS unit (GPS16X-HVS, manufactured by Garmin for Campbell Scientific) was fixed to the top of the box. All instruments were wired into a data logger (CR1000, Campbell Scientific) and programmed to collect data every 3 s. Before each flight, the radiometer and camera lens were cleaned. The box with instrumentation was attached to a helicopter with a 6 m lead line.
Surface reflectance data are collected using a radiometer mounted
underneath a weighted box and flown by helicopter across the valley
Flights originated at Lake Hoare camp. Pilots flew east over Canada glacier,
Lake Fryxell, and Commonwealth glacier. Pilots were instructed to fly off
Commonwealth glacier to the east; turn around; and return over Commonwealth
glacier, the southern side of Lake Fryxell, and the lower elevations of
Canada glacier. They would then continue west to Lake Bonney and Taylor
glacier, turning around at the meteorological station located on Taylor
glacier, and returning to Lake Hoare camp (Fig. 1). We instructed pilots to
fly 100 m off the ground at a speed of 50 knots (26 m s
Five flights were made over the 2015–2016 season (20 November 2015, 7 December 2015, 24 December 2015, 5 January 2016, 12 January 2016), five flights were made over the 2016–2017 season (11 November 2016, 3 December 2016, 14 December 2016, 3 January 2017, 23 January 2017), and four flights were made over the 2017–2018
season (21 November 2017, 6 December 2017, 27 December 2017, 13 January 2018; Bergstrom and
Gooseff, 2019). Solar zenith angles ranged from 54.3 to 61.8
Daily incoming radiation averaged across all met stations (black
line), and minimum and maximum radiation (grey shading) from 1 October to
1 March are shown for the
All spatial analyses were performed using the ArcGIS software package (version 10.7). True color images and radiation data were examined visually to help develop a spatially based classification scheme. All data collected in a flight were classified. Individual measurements were either discarded if they did not meet usability standards or were associated with a given landscape feature and the closest meteorological station. Data were discarded if they were collected at the beginning and end of the flight, as the instruments were being picked up or set down by the helicopter. Data were also discarded if they were too close to the edge of a given landscape feature.
At the edge of a feature, the radiometer is receiving reflected radiation from multiple landscape types, and the measurement cannot be attributed to a specific landscape feature. Outlines of all glaciers and lakes were manually delineated using Landsat8 true color imagery (collected 14 January 2016; see Fig. 2b, c). We then calculated a 100 m buffer from the edge of each landscape feature. The 100 m distance was determined by cross-comparing a range of buffer distances and the relative change in radiation with distance to the edge. The rate of change of reflected radiation dropped off considerably further than 100 m from the edge of a glacier or lake, and we assume that change in reflected radiation is representative of true variability of that surface type.
Each point over a given landscape feature was coded to associate it with one of four lakes (Bonney, Mummy Pond, Hoare, or Fryxell), or one of five glaciers (Taylor, LaCroix, Suess, Canada, or Commonwealth) (Fig. 2b, c). Mummy Pond and LaCroix and Seuss glaciers are small and were not consistently measured. Data from these locations are used in a lumped cross-landscape type comparison (i.e., glaciers vs. lakes) but were not used in an individual analysis (i.e., Canada glacier vs. Taylor glacier) due to lack of sufficient data. Data collected over soil surfaces were separated into two halves of Taylor Valley: the Bonney and Hoare basins (referred to as the west soils) and the Fryxell basin (referred to as the east soils). Flight paths varied, which led to inconsistent numbers of measurements across flights. Glaciers are the main focus of this study and we frequently collected over 80 useable measurements on a given flight for one glacier and a minimum of 49 measurements (Table S1 in the Supplement). Data over soil surfaces were collected as the transition from a lake to a glacier rather than explicitly targeted and therefore had the lowest number of measurements across all flights (Table S1). In some cases, pilots flew too close to the edge of a lake, which necessitated discarding measurements that should be classified as a lake and resulted in low total counts of lake measurements (e.g., 7 December 2015, Table S1).
We used incoming shortwave radiation and snowfall data from meteorological
stations located throughout the valley (Fig. 1). The Bonney, Fryxell,
Hoare, and Canada glacier meteorological stations measure incoming
shortwave radiation with a LI-COR LI-200R pyranometer and the Commonwealth and
Taylor glacier meteorological stations use an Eppley Precision Spectral Pyranometer (PSP).
Incoming radiation is sampled every 30 s and an average is recorded
every 15 min. Stations are leveled and measuring reflected radiation over a
level patch of ice. All meteorological stations are visited annually for
maintenance and instruments are replaced and recalibrated according to the
manufacturers' recommendations. We used snowfall data measured at the
Bonney, Hoare, and Fryxell meteorological stations. Precipitation at lakes
Bonney and Fryxell was measured with a Campbell Scientific SR-50 sonic
ranger. Depth is corrected to volume using a snow density of 83 kg m
We calculated apparent albedo using the incoming radiation measured at the closest meteorological station to a given measurement determined by the Thiessen polygon method (Fig. 2b, c). We used the measurement closest to the time the reflected radiation measurement was made. In a few cases, the closest meteorological station was not recording or was under a cloud (Fig. 3a–c) and produced an apparent albedo over 1. In those cases, we calculated an average incoming radiation from all other meteorological stations to calculate albedo. All data analysis and statistics were done in MATLAB (Mathworks version 2018a).
We tested for significant differences between and across landscape types,
features, and flights. Because points are collected sequentially along a
flight line they are not independent observations and there is inherent
spatial autocorrelation. We address this with transect line correlation
analysis, a statistical method that uses a randomization technique
(Malatesta et al., 1992). Data are first classified as
We compare calculated albedos measured by the helicopter to MODIS-derived albedo estimates. The MCD43A3 version 6 MODIS albedo product provides daily albedo measurements at a 500 m spatial resolution (Schaaf and Wang, 2017). The MODIS instrument, onboard both the Aqua and Terra satellites, images the same location every 1 to 2 d. These satellites have a sun-synchronous, near-polar orbit such that they acquire at least two images of Taylor Valley per day – at roughly 10:30 and 13:30 LT. The MCD43A3 product deploys a bi-directional reflectance distribution function (BRDF) to calculate broadband and spectral albedo at local solar noon using a centered, 16 d moving window, and the center day of the window given the highest weight. Albedo is available as “black-sky” (BSA, direct component) and “white-sky” (WSA, diffuse component) albedo for seven MODIS narrow bands and three broad bands.
Daily MCD43A3 shortwave broadband (0.3–5.0
There are several sources of error associated with these data. The first is
that the instruments had to be in a weighted box rather than mounted to the
helicopter due to safety regulations of the helicopter contractor. We
believe that the fin and the
Under overcast conditions, the surface properties and orientation of the location of interest do not affect albedo measurements, as incoming and reflected radiation is diffuse (Allison et al., 1993; Grenfell et al., 1994; Pirazzini, 2004; Warren et al., 1998). However, under clear skies, there are both direct and diffuse components of incoming solar radiation and therefore the orientation of the sun relative to the surface over which we make measurements can have a large effect on apparent albedo (Grenfell et al., 1994). This is exacerbated at high zenith angles, resulting in more oblique radiation at the land surface. This issue has been addressed in Antarctic albedo measurements at the station scale for uniform sloping surfaces (Grenfell et al., 1994) and for sastrugi (Pirazzini, 2004; Warren et al., 1998). Error estimates of apparent albedo in those conditions are generally less than 10 % but are highly dependent on solar zenith angle, slope, and orientation of the sun relative to the surface aspect.
We applied a simple correction factor to the apparent airborne albedo data
adapted from Grenfell et al. (1994):
Box plots of corrected albedo from all points measured over
Mean and 1 standard deviation of albedo measured over the three
main glaciers and lakes, and soil in the western (Bonney and Hoare basins)
and eastern (Fryxell basin) Taylor Valley during each of the flights.
Apparent albedo is shown in
The issue of complex terrain influence on albedo in snow and ice environments has been addressed from the station and satellite scale for a range of features including penitentes, sastrugi, and small rough alpine glaciers with topographic shading (Dumont et al., 2012; Klok et al., 2004; Lhermitte et al., 2014; Warren et al., 1998). Corrections for satellite data can be very complex and often require knowing the location, field of view, and bidirectional reflectance function for the surface (e.g., Wen et al., 2009). The smaller scale of this airborne albedo data is not well suited to the larger-scale satellite correction methods, due to data requirements and mismatches in the scale over which it is meant to be applied. Albedo corrections at the station scale often require development of an idealized model of the surface and a set of assumptions and numerical descriptions of the optical properties of the surface, which is also unrealistic for this application (Lhermitte et al., 2014; Pirazzini, 2004).
Findings of modeling work by Lhermitte et al. (2014) suggest that the higher the sensor is off the surface, the more it integrates the variability of the surface and lowers the apparent albedo error. They found that a sensor 4 m off the surface converges toward true albedo. We argue that this is true for our data as well, where the position of the sensor 100 m above the land surface reduces some error due to high mean slope. While slopes on average are very high, they are of many different aspects and receive and reflect more light than if they were a uniform surface sloping away from the azimuth of the sun. This results in erroneously low correction factors under the simplified method. The irregularity of the MDV glacier ice surface and sediment cover make correction development a challenging task. More work should be done to develop a correction factor well suited to these conditions; however, we believe that this is beyond the scope of this paper and available data. For this reason, we present both the apparent and corrected albedo using the simplified correction proposed by Grenfell et al. (1994; Figs. 4, 5). We present results for soil and lake surfaces using corrected albedos and discuss corrected and uncorrected (apparent) glacier albedo. When comparing across or within glaciers we use the apparent albedo. We believe, based on previous studies and with comparison to station albedos, that glacier apparent albedo is within 10 % of true albedo. This is within the accuracy required in order to draw conclusions about broad spatial and temporal patterns we observe in the MDVs.
During the austral summer, mean daily radiation increases through December as solar elevation increases, peaking at the solstice, and decreases moving into the fall. We conducted flights as early as 11 November and as late as 23 January, in order to capture the seasonal changes in albedo. Therefore, maximum potential incoming radiation varies across flights as a function of day of year. Throughout each season, there is high variability in mean daily radiation due to cloudiness and individual storms (Fig. 3a–c). We attempted to schedule flights during completely cloud-free days and were successful during the 2016–2017 season (Fig. 3e). In some cases of persistent cloudiness, we attempted to make flights under uniform cloud cover such as the last two flights of the 2015–2016 season. There is limited variability in incoming radiation across the valley on these 2 d (Fig. 3d). Cloud cover is known to increase albedo for snow- and ice-covered areas and decrease it for other surfaces by several percent (Key et al., 2001). We did not correct albedo on these cloudy days. The goal of this study is to describe broad patterns in albedo rather than exact measures. We believe that the cloudy flights will not influence the main conclusions of this study but do take it into account in our interpretation of results.
The Canada and Commonwealth glaciers flow out of the Asgard Range from the
north, and lakes Hoare and Fryxell are also on the northern side of the
valley (Fig. 1). The valley also narrows and steepens toward Taylor
glacier. Thus, meteorological stations are influenced by topographic shading
early and late in the day, particularly by the Asgard Range (Fig. 3g).
We conducted flights between 11:00 and 13:00 LT to avoid topographic shading
and minimize zenith angle and variability in azimuth. Zenith angle varied
less than 10
We typically collected over 100 points per landscape type, resulting in high
spatial resolution of albedo measurements synoptically across the entire
valley (Table S1). Across landscape types there was distinct separation in
the 2015–2016 season (Fig. 5a, d). There were significant differences in
reflected radiation across glaciers, lakes, and soils. Glaciers consistently
had the highest albedo, and soils had the lowest albedo (
In order to understand how albedo varies along a transect from west to east along the valley with decreasing distance to the ocean, we plotted albedo mean and standard deviations of each of the three main glaciers and lakes (Fig. 5). Soil data are divided into the western (Bonney and Hoare basins) and eastern (Fryxell basin) parts of the valley. We also plot measurable snow events observed at Bonney, Hoare, and Fryxell meteorological stations as snow “hyetographs”. The snow events, particularly in the 2016–2017 season, have a dramatic effect on the valley-wide albedo patterns.
The first flights of the 2015–2016 and 2016–2017 seasons had similar albedo patterns. All glaciers and lakes have the lowest albedo measured all season (Fig. 5a, b). Lake Fryxell has significantly higher albedo than lakes Hoare and Bonney. On the second flight of the 2015–2016 season with no measured snowfall between flights, we observe that albedo of all lakes increased and that there is no difference between the lakes (Fig. 5d). Albedo of the glaciers did not change significantly.
On 22 and 23 November 2016, a very large, valley-wide snow event occurred. The greatest accumulation was observed at the Lake Fryxell meteorological station. Accumulation decreased inland toward the Bonney basin (Fig. 5b) and persisted for several weeks. The snow event is evident in albedo measured during the subsequent two flights. The longitudinal increase in snowfall with proximity to the coast coincides with increased albedo moving east toward the coast. There is increasing albedo from Taylor to Canada to Commonwealth glaciers (Fig. 5b). Lake Fryxell has the highest measured albedo, even higher than all three glaciers. Similarly, eastern soils have higher albedo than lakes Bonney and Hoare and Taylor glacier. Mean albedo over eastern soils on 3 December was 0.60, indicative of substantial snow cover. Both east and west soil albedo standard deviation increased in the two flights after the snow event but did not markedly increase for the glaciers or lakes. A total of 11 days after the snow event, we still observe an east–west pattern of albedo. Lake Fryxell had the highest albedo of any feature. While albedo of eastern soils decreased slightly, mean albedo remained at 0.49, which was still higher than both lakes Hoare and Bonney. The albedo of Canada and Commonwealth glaciers actually increased (Fig. 5b). Finally, albedo of Lake Fryxell remained elevated for the rest of the season relative to albedo measured over Lake Fryxell at the end of the previous snow-free season (mean corrected albedo was 0.54 and 0.58 vs 0.46 and 0.48 on the last two flights of the 2016–2017 and 2015–2016 seasons, respectively).
One of the most prominent differences in longitudinal patterns of albedo was measured on 22 November 2017. We observed an along-valley increase from west to east in albedo on that flight similar to what was observed after the snow event the previous season described above. The last snow event recorded prior to the 22 November 2017 flight was in August. At the beginning of the 2017–2018 season there were very large snowdrifts across the valley against the toes of glaciers, in incised stream channels, and in the lee of camp buildings. We found on this first flight that while western soils did not have elevated albedo, eastern soils have a mean albedo of 0.64 (Fig. 5f). Commonwealth glacier also had the highest mean albedo observed on any flight across all three seasons (0.74, Fig. 5). By early December 2017, the albedo of all glaciers, lakes, and soils had decreased. Several small snow events occurred in December, but they did not have a clear spatial effect like what was observed in the previous season. The last two flights of the 2017–2018 season did have elevated albedo on glaciers and lakes relative to the previous two seasons, possibly attributable to snow and cold temperatures restricting melt (Fig. 6). Overall, the three seasons provide contrast in spatial and temporal patterns of albedo due to the weather that occurred during each. The 2015–2016 season was a clear and relatively snow-free season. The 2016–2017 season had a very large and persistent snow event early in the season. The 2017–2018 season started with large snow drifts from winter snow and subsequent wind redistribution. This facilitates discussion of how the landscape responds to each type of albedo-altering event.
Time series of albedo measured at the Taylor, Canada, and Commonwealth
meteorological stations over the
We hypothesized that there would be a positive relationship between albedo and ice surface elevation. Following standard principles of adiabatic cooling, more melt will be generated at lower elevations, lowering albedo due to the presence of liquid water, which creates melt pools, exposes and concentrates sediment, and creates increased surface roughness (Bagshaw et al., 2010; Hoffman et al., 2016). We compare apparent albedo to ice surface elevation of each measurement collected on each flight and season over Taylor (Fig. 7) Canada (Fig. 8), and Commonwealth glaciers (Fig. 9). We also plot hypsometric curves for each glacier to compare how much glacier area is represented by the measurements. Albedo measured at the meteorological station on that glacier at the time of the flight is also plotted as the larger colored point in each plot.
Apparent albedo plotted against surface elevation of Taylor
glacier for all flights in the 2015–2016 season
On almost all flights, we have a measurement directly above the meteorological station. On the third flight of the 2016–2017 season, the pilot turned around over Lake Bonney, and we do not have measurements of albedo over Taylor glacier. The airborne measurements collected at the elevation of the meteorological station are comparable to the station-based measurements for all flights over Taylor glacier giving us confidence in the use of apparent albedo for these comparisons. On the flights with little variation in albedo, the ground-based measurement is fairly representative of the entire ablation zone albedo. However, over Taylor glacier we observed an albedo increase with increasing elevation on several flights, particularly the second flight of the 2017–2018 season (Fig. 7k). Additionally, at lower elevations we observe a higher apparent albedo on one flight line relative to the other, with a difference as much as 0.10 on some flights (Fig. 7). There tends to be a convergence of albedo between the two flight lines at around 300 m elevation. Taylor glacier has lower variability in apparent albedo between flight lines and from flight to flight than Canada and Commonwealth glaciers.
Canada and Commonwealth glaciers have very clear separation between the two flight paths across the glaciers (Figs. 8, 9). On Commonwealth glacier, there is very distinct separation between the upper and lower flight lines by elevation. The lower flight line rarely reaches the elevation of the upper flight line (Fig. 9). The shape of both Canada and Commonwealth glaciers is such that on a transect across the ablation zone, perpendicular to glacial flow, the highest elevation is in the center and lowest on the edges (Figs. 8 and 9). We find that for a given elevation there is a higher albedo on one side of the glacier relative to the same elevation on the other side (Fig. 8b). We plot the locations of measurements taken over each glacier with their color representing albedo (Figs. 10–12). There is a clear longitudinal difference in albedo on the Canada and Commonwealth glaciers on most flights. We believe these differences are larger than what would be caused by differences in aspect and solar illumination across the two sides. The western side of both glaciers tend to have lower albedo than the eastern sides, and is also visible in the imagery. Figures 11n and 12h are good examples of this pattern. The last several measurements collected on the eastern edge of the Commonwealth glacier decrease slightly. The darker edge of the glacier is visible in the true color imagery (Fig. 12).
Apparent albedo plotted against surface elevation of Canada
glacier for all flights in the 2015–2016 season
Apparent albedo plotted against surface elevation of Commonwealth
glacier for all flights in the 2015–2016 season
Flight paths of each flight on Taylor glacier in the 2015–2016
season
Flight paths of each flight on Canada glacier in the 2015–2016
season
Flight paths of each flight on Commonwealth glacier in the
2015–2016 season
While the station-based measurement of albedo is more frequently representative of Taylor glacier (Fig. 7), it is typically near the highest airborne albedo measured over Commonwealth and Canada glaciers (Figs. 8 and 9). The meteorological stations on these glaciers are both located in the center and toward the higher elevations of the ablation zones (Figs. 1, 8, 9). According to our airborne measurements, we generally see an increase in albedo with increasing elevation and from west to east across the ablation zones. Therefore, these meteorological stations give us confidence in using apparent albedo but, due to their siting, are not particularly representative of the albedo of the entire Canada and Commonwealth glacier ablation zones.
We compared MODIS albedo data to both apparent and corrected albedo (Fig. 13). We found that overall fits between the MODIS data and airborne data are
better using apparent (uncorrected) albedo. This increase in root-mean-square error using corrected data is likely due to the unrealistically high
values of corrected glacier measurements. In the uncorrected data, most soil
points fell below the
Comparison of all airborne albedo measurements to MODIS data
across all surface types. We compare both
Snowfall alters the energy balance of a surface because it has different thermal and optical properties than the underlying ice. Therefore, one would expect there to be an effect on the seasonal aging of ice under even minimal snow cover of moderate duration. There was almost no measurable snow in the first season in which we conducted this study. This allows us to make observations of albedo seasonal change due to ice aging without the influence of snow over both lake and glacial ice. We found that in the 2015–2016 season, glacier and lake ice albedo increased from the first flight in late November to the third flight in late December (Fig. 5). This albedo increase is what some refer to as ice whitening and is a result of ice structural change, which can increase scattering of light (Fritsen and Priscu, 1999; Howard-Williams et al., 1988). The changing grain and crystal shape is a function of the initial structure and the surface energy balance (Adams et al., 1998; Grenfell, 1983; Warren, 1982; Wiscombe and Warren, 1980). During the period of albedo increase we observed over the MDV lakes and glaciers from November to late December, total mean daily solar radiation and air temperature increased, which resulted in more energy available to warm, melt, and alter the ice surfaces (Figs. 3a, 6a). The exact mechanisms of ice whitening differ on lake and glacier surfaces due to the differences in how the ice was formed and processes governing melt. The process of ice structural change and its effect on albedo and the transmission of light has been well studied for MDV lakes. McKay et al. (1994) developed a conceptual model for the seasonal evolution of Lake Hoare ice cover. This model suggests that early in the season the ice optical properties are nearly constant. However, as solar radiation and temperature increase, albedo also increases due to fracturing, internal melting, and the development of light-scattering Tyndall structures (Howard-Williams et al., 1988; Knight and Knight, 1972). This increased albedo was coupled with decreased transmission of light peaks in midsummer at the peak of incoming solar radiation (Figs. 3, 4, 5; McKay et al., 1994). We observe that lake ice albedo decreased in the last two flights of the 2015–2016 season (Fig. 5d). By midsummer, enough melt has occurred that the light-scattering melt structures are destroyed as an aquifer develops within the ice cover (Fritsen and Priscu, 1999; Howard-Williams et al., 1988; McKay et al., 1994). Eventually as solar radiation and air temperatures begin to decrease, the ice will refreeze, redevelop structure, and increase scattering and albedo (Figs. 3, 6). We believe that we did not have a flight late enough in the season to capture this process (which is in part a limitation of our field access, ending in late January each austral summer).
The mean lake ice albedo we measured in the 2015–2016 season follows the
conceptual model of McKay et al. (1994). However, melting and surface
optical properties of lake ice are very complex due to meter-scale surface
roughness and patchy surface characteristics. Sediment concentration on lake
ice ranges from 0.2 to 2 g cm
The ice structure and melt processes are very different for the MDV glaciers vs. lakes. The energy balance dynamics generating melt have been extensively modeled (Hoffman et al., 2014, 2016; Lewis et al., 1998), but the result of melt on the optical properties of the ice has not been well described. Hoffman et al. (2014) found that melt is generated in the shallow subsurface of the glaciers (5–10 cm below the surface) and drains to topographic lows, reducing the ice density. As ice density is lost, light can be scattered in the interstices of the ice grains and increase albedo (Grenfell, 1983; Grenfell and Maykut, 1977). Another potential mechanism for albedo increase is due to sublimation (Winkler et al., 2009). Sublimation makes up a high proportion of total ablation on the glaciers and is by far the dominant mechanism of ablation on the glacier surface relative to melt (Lewis et al., 1998; Hoffman et al., 2014). However, there are mixed effects of sublimation (Winkler et al., 2009). Sublimation can reduce ice density on daily timescales and create more microtopography, which would increase scattering and can change on sub-daily timescales (Pirazzini, 2004). Conversely, sublimation-derived microtopography could increase the ability for the glacier to trap sediment, and scatter light back into neighboring ice, thereby reducing albedo (Lhermitte et al., 2014; Winkler et al., 2009). This effect is likely spatially variable and depending on the degree of sublimation and melt in a given location would require more extensive modeling and is the subject of future work.
Seasonal increase in ice albedo has been observed on glaciers elsewhere including temperate glaciers in the Alps and polythermal glaciers in the Arctic (Knap and Oerlemans, 1996; Oerlemans and Knap, 1998). However, this albedo increase has largely been attributed to the formation of superimposed ice, where meltwater derived from snow refreezes upon reaching colder ice. The MDV glaciers generally do not develop a seasonal snow cover in winter so superimposed ice cannot form. However, the principle of refreezing meltwater is probably still applicable on these glaciers (Hoffman et al., 2016). While some meltwater generated on the glaciers drains through the ice matrix into supraglacial channels and leaves the glacier as terminal waterfalls, some meltwater will refreeze before it can leave the glacier and likely contributes to the increased albedo. Temperatures generally remain below 0 throughout the seasons, promoting refreeze (Fig. 6d–f).
The glaciers had a slight decrease in albedo in the latter half of the season, particularly on Canada and Commonwealth glaciers (Fig. 5a). The flight lines were not as consistent on Taylor glacier, making it hard to compare and understand temporal shifts in albedo on this glacier. In the last two flights over both Canada and Commonwealth glaciers the flight line traversing upper elevations had the largest albedo decrease, while measurements made over lower elevations were more constant (Figs. 11 and 12c–e), also shown in the albedo time series measured at the meteorological stations (Fig. 6a–c). The lower elevations of the glaciers have higher surface roughness and slope and are drained by supraglacial channels. These lower sections may have more stable albedo because they can generate more meltwater, effectively drain it, and are not as affected by melt–refreeze processes (Hoffman et al., 2016). Higher-elevation glacier surface is characterized by smooth ice with cryoconite holes. Cryoconite holes are locations where wind-deposited sediment has concentrated and melted into the ice forming a cylindrical depression with an approximately 1 cm thick layer of sediment on the bottom, overlain by several centimeters of liquid water (Fountain et al., 2004). Antarctic cryoconite holes tend to have ice lids that by late season can melt, leaving open holes where the light-transmitting liquid water and darker sediment is exposed to the atmosphere (Fountain et al., 2004, 2008). This late season opening of cryoconite holes driven by a short period of warm temperatures (Fig. 6d) may be the cause of the reduced albedo we observed in the 2015–2016 season. We suggest that seasonal trends in glacier and lake ice albedo in the absence of snow are a function of the seasonal energy balance and the physical ice structure from grain scale to entire lake and glacier scales.
Soil albedo patterns are much different to those of lakes and glaciers due to their fundamental difference in material: rock vs. ice. In 2015–2016, soil albedo decreases slightly from late November to late December (mean albedo: 0.18 to 0.12, Fig. 5d) possibly due to melt of residual snow patches. Soils can have decreased albedo when wet (Bøggild et al., 2010). In addition to snowmelt, hypersaline patches of soil MDVs can become wet due to deliquescence and may contribute to soil albedo decrease over the season without a snow event (Levy et al., 2012). There is strong variability in soil albedo patterns across seasons (due to snow) but across all seasons soil albedo is lowest at the end of the season (Fig. 4c).
Snow typically has a higher albedo than ice, and therefore snowfall in Taylor Valley raises the albedo of every landscape feature. The albedo of snow is strongly influenced by grain size, snow age, and history, and, in addition, snow must be of sufficient water equivalent depth for a surface albedo to be set by the snow alone (Wiscombe and Warren, 1980). Snow events are small in the MDVs, only producing a few millimeters of water equivalent in a given event, typically follow a gradient of increased snowfall toward the mouth of the valley, and are generally quickly redistributed by wind (Doran et al., 2002; Fountain et al., 2010). We broadly discuss the role of snow in altering the albedo of each landscape feature and discuss some possible mechanisms for the patterns we observe.
On 22–23 November 2016, Taylor Valley received one of the largest snow events on record, with the highest recorded snowfall toward the valley mouth in the Fryxell basin (Fig. 5b). This spatial gradient in snowfall was reflected in the albedo of the landscape features. Taylor glacier, Lake Bonney, and soils in that basin only saw modest increases in albedo. In addition, the Bonney and Hoare basin soils had increased standard deviation (Fig. 5e). It is possible that the underlying soil surface may be contributing an albedo signal mixed with the thin snow cover, which is variable depending on the soil surface structure (distribution of coarse fragments and boulders) as well as possible variability in snow depth due to wind redistribution. Lower in the valley, there is a much larger increase in albedo after the snow event (Fig. 5e). Again, this may be both a function of the initial accumulation as well as potential wind redistribution. The strongest winds tend to move off the polar plateau along the valley (Doran et al., 2002; Nylen et al., 2004), moving snow from west to east. These increased snow depths likely cover the lower valley more completely, reducing surface roughness and masking the signal of low albedo underlying surfaces, resulting in a larger increase in albedo (Warren et al., 1998).
Typically, as snow ages grain size increases, impurities accumulate (dust and black carbon), and albedo decreases. This has been shown empirically and has been used in models across scales to reasonably reproduce measured albedo time series (Oerlemans and Knap, 1998; Pellicciotti et al., 2008; Warren and Wiscombe, 1980; Wiscombe and Warren, 1980). We hypothesized, following this principle, that albedo would decrease across the landscape in subsequent flights after a snow event. Surprisingly, we found mixed results where the soils mean albedo significantly decreased, and Lake Hoare decreased slightly but not significantly. Canada and Commonwealth glaciers and Lake Bonney increased significantly, while Lake Fryxell mean albedo increased but not significantly (Fig. 5b, e). There was a period of high solar radiation and two short spikes in air temperature between the flights that may have promoted melt or sublimation, particularly for the soils with thinner snow cover (Figs. 3b, 6e).
Variability in the albedo change between lakes could be due to a combination of the gradient in snow accumulation and sediment and snow capture efficiency of the lake ice. Lake Bonney historically has the lowest sediment load (Sabacka et al., 2012) and receives the least amount of snow (Fountain et al., 2010) leaving more exposed ice cover that is subject to ice whitening independent of the snow event. Lake Hoare typically has the greatest sediment load (Sabacka et al., 2012), which promotes snow melt and, once exposed, has lower albedo. The increase in albedo at Lake Fryxell could be from snow deposited on the lake surface from wind.
We compared the individual measurements on Canada and Commonwealth glaciers to understand if there is a spatial pattern of albedo on the glaciers post-snowfall. We found that the lateral difference in albedo for a given elevation is very high (Figs. 8h and 9h). The western sides have lower albedo than the eastern sides of both glaciers on most flights; however, the difference became more pronounced for the 3–14 December 2016 flights (Figs. 11g, h and 12g, h). The Canada and Commonwealth glaciers both flow to the valley bottom and present major topographic obstacles to wind (Fig. 1). It is likely that any wind event that occurred between the snow event and subsequent flights transported snow from the valley floor and windward side of the glaciers and deposited them on the leeward side of the glaciers. We believe that the aeolian redistribution of snow enhanced the spatial gradients in albedo across the glaciers and raised mean albedo. The lakes are not a topographic obstacle to wind and this may be why there was a less pronounced change over the lake ice surfaces.
Aeolian transport of material is one of the most important drivers of connectivity across the dry valley landscape and contributes to the morphological and energy balance characteristics for a given landscape position (Sabacka et al., 2012). We found aeolian transport is not only important for driving snow redistribution on the event scale but also may drive patterns of albedo at the seasonal scale. On the first flight of the 2017–2018 season we observed a large difference in albedo of the Fryxell basin relative to Bonney and Hoare basin lakes, glaciers, and soils (Fig. 5f). Taylor and Canada glaciers, as well as lakes Hoare and Bonney and associated soils, all had albedo comparable to the first flights of the previous two seasons, although there is a high standard deviation measured over Lake Bonney. In the Fryxell basin, Commonwealth glacier had the highest mean albedo, and both Commonwealth glacier and Lake Fryxell had the highest standard deviation measured across all seasons. The last measurable snow event occurred in August 2017. Upon arriving in the valley at the beginning of the 2017–2018 season, we found many meter-scale snowdrifts in any location that was a windbreak (i.e., incised streambeds, leeward sides of buildings, terminal glacier fronts), which persisted long into the summer.
During the winter, Fountain et al. (2010) suggest that due to both the gradient of decreasing snowfall with distance from the coast and the foehn-wind-driven along-valley transport and redistribution of snow from valley walls to floors, snow does not persist for very long in the Bonney basin but can settle in persistent drifts in the Fryxell basin. Our albedo measurements agree with their observations. The high standard deviations we observe are likely from the drifted snow, where there are large drifts of high albedo snow in topographic lows interspersed with patches of lower-albedo bare ice and soil. Longitudinal patterns of increasing albedo from west to east existed on the Commonwealth glacier, but the trend was not as smooth as observed on other flights (Fig. 12k). In addition, the lowest observed albedo of that flight was near the meteorological station, which is at a high elevation with the smoothest ice and is therefore the least likely to hold snow in high winds (Fig. 9k). Interestingly, all components of the Fryxell basin had albedo comparable to the rest of the valley by the next flight on 27 December 2017. There was a clear, sunny period between the two flights that may have promoted melting and sublimation of the snowdrifts (Fig. 3c). Our findings demonstrate that snow can have a strong influence on albedo across the landscape. We also suggest that redistribution via wind and melting or sublimation of snow can result in large spatial and temporal variability of these patterns. This has been found across the Antarctic continent (Pirazzini, 2004).
Albedo is one of the most important parameters for glacial energy balance across all types of glaciers at all latitudes (Brock et al., 2000; Brun et al., 2015; Dumont et al., 2011; Mölg and Hardy, 2004; Pellicciotti et al., 2008). Remotely sensed observations have been used successfully in many parts of the world, particularly at large scales (e.g., Knap and Oerlemans, 1996; Pope and Rees, 2014). For the McMurdo Dry Valley glaciers, satellite imagery is particularly challenging to use effectively because of the oblique angles of most images, topographic shading, and strong atmospheric interference. Additionally, the glaciers are only several kilometers across so the resolution of most images is not sufficient to capture a pure glacier pixel that is actually representative. Use of satellites to calculate glacial energy balance components in the dry valley region has been limited to a few studies (Bliss et al., 2011; Dana et al., 2002). Largely, albedo in glacial energy balance models has been derived from meteorological stations (Hoffman et al., 2008, 2014, 2016). We compare airborne albedo measurements to both station and remotely sensed (MODIS) data collected during the same periods.
The meteorological stations were sited at high elevation in the ablation zone of the glaciers and on horizontal ice with low surface roughness. Hoffman et al. (2014) suggest that the location of these radiometers may not measure the most representative albedo for the entire glacier surface and cause mass balance models to underpredict glacier ablation. We believe that the helicopter-based measurements of albedo are comparable to meteorological station measurements and therefore allow us to evaluate how albedo measured at a station compares to the distribution of albedo across the ablation zone of an entire glacier (Figs. 7–9). With very few exceptions, the airborne measurements return the same albedo as the station-based measurements when directly over the meteorological station. However, during most flights the station-based measurements were not particularly representative due to spatial trends in albedo (Figs. 7–9). The airborne albedo measurements over glaciers revealed three key findings: (1) there is generally a positive relationship between albedo and elevation; (2) for a given elevation there can be a gradient in albedo laterally across a glacier; and (3) the strength of these trends can vary through time, presumably with other energy balance components, such as solar radiation, temperature, and wind speed (Figs. 3–6).
We suggest adiabatic lapse and the change in ice structure from drainage-dominated to flat smooth ice resulted in the frequently observed positive relationship between albedo and elevation, common across all three of our study glaciers (e.g., Figs. 7k, 8d, and 9i). We also observe across all three glaciers that for a given elevation there can be a major disparity in measured albedos. On Taylor glacier, major drainage channels begin to develop roughly 3.5 km from the terminus, 300 m a.s.l., and drain to the terminus parallel to glacier flow. These channels can be tens of meters wide and deep near the terminus and hold darker melt pools and thick layers of sediment. We believe that the difference in albedo we observe for a given elevation on Taylor glacier depends on whether the instrument was over a channel or smooth ice ridge between drainages.
For the Canada and Commonwealth glaciers, different albedo for a given
elevation is due to lateral differences in albedo where the west sides tend
to be darker than the east sides. The west sides of these glaciers are the
upwind side for foehn events that can deposit up to 1 kg m
MODIS data provide more confidence that the airborne measurements are representative of true albedo over these surface types. The corrected soil and lake data agree well with MODIS-derived albedo data. Apparent (uncorrected) glacier data agrees better with MODIS albedo data than corrected data, due to complexities described in Sect. 3. MODIS pixels taken at a 500 m resolution. This means that an individual value is integrating over a much larger area and has issues with mixed pixels, particularly toward narrow parts of the features (e.g., upper elevations of Canada and Commonwealth glacier ablation zones, Fig. 1). There are many individual airborne measurements contained within the same MODIS pixel. As shown by the vertical linear striped appearance of the data in Fig. 13, there is a range of airborne values for any given MODIS value. This illustrates the utility of the airborne data in that they are capable of picking up the smaller-scale spatial variability of the land surface. Airborne measurements also do not have issues with mixed pixels like MODIS does. Lakes and glaciers would have artificially low values if part of the pixel also covers soil and vice versa for soils having artificially high albedo. The best relationships between MODIS and airborne data are following snow events (Figs. S1, 2). The broad valley-wide patterns in snow distribution result in clearer albedo differences at both measurement scales. As snow melts, soil and sediment patches are exposed, ice ages, and relationships deteriorate. This mismatch in measurement scale is apparent here as patchiness increases and albedo spatial variability increases. More work should be done to develop a predictive model of glacier albedo capitalizing on the strengths of all three measurement scales (station, airborne, and MODIS) and to incorporate it into spatially distributed energy balance models in order to better predict ablation and meltwater generation from these glaciers.
The transfer of energy and water among landscape units is critical for sustaining life in this polar desert. The glaciers supply meltwater to the downstream environment, the generation of which is determined by the energy balance. Albedo is a key parameter in the energy balance equation for all glaciers but is particularly key in this landscape where melt is primarily driven by solar radiation rather than sensible heat. The amount of energy that is absorbed or reflected also has important ecological implications across the landscape from the under-ice primary productivity in lakes to the thermal stability of soils for invertebrate populations (Fritsen and Priscu, 1999; Wlostowski et al., 2018). While we have the benefit of many meteorological stations dispersed throughout the valley, only those located on glaciers measure albedo, and they are limited in that they are point measurements in a highly spatially variable landscape.
The meteorological stations appear to sufficiently capture variability in incoming radiation, which on a seasonal basis varies with cloudiness (Fig. 3) and varies on annual timescales with atmospheric concentrations of pollutants (Obryk et al., 2018). The amount of reflected radiation appears to vary seasonally with the amount of incoming radiation, snow events, wind events, and the physical characteristics of the landscape. We suggest that wind-driven snow and sediment redistribution is a critical factor in the long- and short-term spatial patterns of albedo that we observed: broadly increasing albedo both from west to east along the valley and individually across the Canada and Commonwealth glaciers. Wind drives connectivity across the landscape and is a feedback mechanism for the energy balance of a given landscape unit. By delivering sediment and redistributing snow, it promotes meltwater generation and drives further connectivity. These longitudinal patterns in albedo indicate that point measurements made at a meteorological station are not sufficient to fully describe the albedo of a given ice or soil environment and that one can expect a very different energy budget from one basin to the next. This study is a step toward understanding how we can scale our current measurements to the entire valley and the factors we would need to consider in order to do so. Further development of this dataset and methods for correcting and coupling it with other albedo data (i.e., meteorological stations and remotely sensed data) and modeling work, particularly incorporating redistribution of snow and variability in surface roughness, should be done to understand the impacts of albedo variability on the meltwater generation from glaciers as well as the energy budget for lake ice and soils.
Meteorological and albedo datasets are available on the MCM LTER database:
The supplement related to this article is available online at:
AB and MG designed the study and carried out data collection. AB performed data analysis. MM contributed snow data. JC contributed MODIS data analysis. AB prepared the manuscript with contributions from all co-authors.
The authors declare that they have no conflict of interest.
The authors would like to acknowledge the civilian contractors and Petroleum Helicopters, Inc. that supported the US Antarctic Program and make this data collection possible. Any opinions, finding, conclusions, or recommendations expressed in the material are those of the authors and do not necessarily reflect the views of the National Science Foundation.
This research has been supported by the National Science Foundation, Office of Polar Programs (grant nos. 1115245 and 1637708). Support was also provided by a National Science Foundation Graduate Research Fellowship awarded to Anna Bergstrom. Publication of this article was funded by the University of Colorado Boulder Libraries Open Access Fund.
This paper was edited by Tobias Sauter and reviewed by Rainer Prinz and Stephen Warren.