Decadal to seasonal variability of Arctic sea ice albedo
A controlling factor in the seasonal and climatological evolution of the sea ice cover is its albedo . Here we analyze Arctic data from the Advanced Very High Resolution Radiometer (AVHRR) Polar Pathfinder and assess the seasonality and variability of broadband albedo from a 23 year daily record. We produce a histogram of daily albedo over ice covered regions in which the principal albedo transitions are seen; high albedo in late winter and spring, the onset of snow melt and melt pond formation in the summer, and fall freeze up. The bimodal late summer distribution demonstrates the combination of the poleward progression of the onset of melt with the coexistence of perennial bare ice with melt ponds and open water, which then merge to a broad peak at 0.5. We find the interannual variability to be dominated by the low end of the distribution, highlighting the controlling influence of the ice thickness distribution and large-scale ice edge dynamics. The statistics obtained provide a simple framework for model studies of albedo parameterizations and sensitivities.
AGARWAL ET AL. \titlerunningheadSEA ICE ALBEDO \authoraddrJ. S. Wettlaufer, Yale University, New Haven, CT 06520-8109, USA.(firstname.lastname@example.org)
The rapid decline of Arctic sea ice coverage during the last thirty years has led to substantial discussion and debate regarding the direct and indirect roles of the ice-albedo feedback, and other processes, in driving the potential loss of a perennial sea ice state [[, e.g.,]]OneWatt, SerrezeBarry:2011. The focus on ice-albedo feedback is based on intuition, modeling and observation [[, e.g.,]and refs therein]Walsh:2000. We understand that there are many mechanisms, operating over a wide range of length and time scales, that influence environmental change in the Arctic [[, e.g.,]]MLT:2009, SerrezeBarry:2011. Moreover, even in simple models the underlying operation of the direct and indirect ice-albedo feedback is a complex process [[, e.g.,]]Curry:1995, EW09, Renate:2011 and this has implications for Global Climate Models (GCMs) which depend sensitively on the treatment of ice albedo []EUW:2007, EUW:2008, DeWeaver:2008. Indeed, the fact that the IPCC AR4 GCMs underproject the decline of the ice during the satellite era [[, e.g.,]]OneWatt, SerrezeBarry:2011, Ian:2011 focuses attention on the methodology used to treat albedo in future model projections. Field based measurements are essential for building physically based models of the processes that influence surface broadband albedo [[, e.g.,]]Perovich:2002, Skyllingstad:2009, as well as for ground truthing of satellite retrieval algorithms [[, e.g.,]]Stroeve:2001. However, sparse coverage of in-situ measurements highlights the importance of combining field and space based methods to examine the spatiotemporal trends in sea ice albedo that provide an essential test bed for models []Perovich:2007.
Here we describe an analysis of daily satellite retrievals of the directional - hemispheric apparent albedo from the Advanced Very High Resolution Radiometer (AVHRR) Polar Pathfinder (APP). (The apparent albedo is what would be measured by upward and downward looking radiometers and hence varies with the state of the atmosphere and the solar zenith angle. For simplicity we refer to this as .) By culling out the albedo retrieved solely from the ice covered region of the Arctic Ocean from 1982 to 2004 a pixel scale view of the surface is created from which we produce a daily histogram. The principal seasonal surface transitions in the ice cover are observed; the high albedo snow covered spring to snow melt convolved with the poleward progression of melt onset, the appearance of melt ponds, their coarsening and fall freeze up (Figures 1-3). The data reveal the decadal and intra- and inter-seasonal variability of the albedo over multiple time scales. It is hoped that their statistics and trends will be of use in the development of models ranging from the very simple to GCMs.
The APP dataset has been refined and applied for use in a wide range of polar studies and is described in detail on the NSIDC website []NSIDC. AVHRR channels 1 - 5 range from the visible to the thermal infrared (0.58 - 12.5 m) and measure Top of the Atmosphere (TOA) reflectances (1, 2, 3A) and brightness temperatures (3, 3B, 4, 5). Additional information includes the solar zenith, relative azimuth and satellite elevation angles, cloud and orbit masks and universal time [[, see Figures 3 and 4 of]]NSIDC. For the analysis reported here we use the so-called Surface Type Mask and retrievals from 1 January 1982 through 31 December 2004, taken daily at 1400 hours. The albedo retrieval involves four steps developed and tested by Csiszar and Gutman (1999): (a) normalize channels 1 and 2 with respect to solar zenith angle, (b) convert channels 1 and 2 narrow band reflectance to a TOA broadband reflectance, (c) correct the TOA broadband reflectance for anisotropy, and (d) convert the TOA broadband reflectance to a clear sky surface broadband albedo.
The Northern Hemisphere is gridded in 1805 1805 pixels, with each pixel representing a 5 km 5 km region. Sea ice is distinguished from land and open water using SSM/I brightness temperatures and filtering the Surface Type Mask data with the NASA Team Sea Ice Algorithm which distinguishes between first-year (FYI) and multi-year ice (MYI) concentrations. The approach ascribes a MYI flag to a region containing at least 50 percent of this ice type. We note that uncertainties depend on the season (e.g., melt pond fraction) and region (near the ice edge), along with the surface type categories Comiso and Kwok (1996). In particular, as will become relevant in our interpretation, errors are enhanced in the presence of surface ablation (especially during snow melt) because a small amount of liquid water alters the emission characteristics of the surface and it is difficult to distinguish between open water at the freezing point and wet ice surfaces. In their ground truthing study of this retreival scheme Stroeve et al. (2001) note that the satellite albedo and that measured by radiometers at the surface are equivalent and are functions of atmospheric conditions, solar zenith angle and the satellite elevation angle. Finally, the cloud masking scheme was ground truthed by Maslanik et al. (2001) during the most difficult part of the cloud seasonal cycle (April-July) over boxes from 15 to 305 km on a side. The maximum standard deviation of albedo was found to be 0.05 which reduced precipitously as the length scale decreased.
On physical grounds the albedo data are filtered to remove any values greater that 1 or less than 0.2. We examine each pixel every day for the presence of ice. Then, for a specified time threshold , we average for that pixel. Different thresholds specify the minimum time for which a pixel contained ice. A histogram for each day of the year is thereby produced and represents the number of pixels for each average albedo value bin for that day. Different colors on the histogram specify different thresholds. Red indicates that a pixel has contained ice for the whole time period under consideration, i.e. = 23 yr (1982 through 2004). Black specifies if a pixel has been ice for at least = 22 yr in the entire time period, which helps us to remove any bias introduced by any minimum during a year. White denotes = 17 yr; blue = 12 yr; green = 8 yr; and yellow = 1 year in the entire time period. Our approach is motivated by several factors. (a) The knowledge of the great sensitivity of GCMs to albedo parameterizations Eisenman et al. (2007). (b) Our continued development of theoretical treatments for large scale transitions in the ice cover that test the robustness of bifurcations in the state of the system to the ice-albedo (and other) feedback Moon and Wettlaufer (2011). (c) Extending these theories to understand the role of noise. Whence, an observational understanding of the temporal variability, from month to month and interannually is essential, and is simply captured using the threshold idea.
3 Results and Discussion
The results are shown in Figures 1-3 and we discuss first the entire record shown in red ( = 23 yr). Because of the limited spatial resolution of the retrieval and the uncertainties in ice type detection and albedo estimation, we do not expect that the numerical values averaged over the ice-covered regions as determined here to map directly onto field observations, but our results nonetheless nicely display the five distinct phases of albedo evolution from spring through fall: dry snow, melting snow, pond formation, pond evolution, and fall freeze-up (see Figure 3 of Perovich et al., 2007b, and references therein). The dry snow value seen here ( 0.75) will by definition be lower than a field observation on MYI ( 0.85) because our pixels contain both MYI and FYI , which have varying amounts of snow cover Kwok and Untersteiner (2011), and near the ice edge the pixels include open water. For March, April and May has a single strong peak about a value indicative of snow covered ice, with an asymmetry that shows a sharp decay at the high values and a broader decay at low values capturing the range of ice types and snow covers that are expected over the basin in the spring.
By the middle of June the distribution is clearly bimodal, with a high peak at about the melting snow covered MYI value ( 0.70) and melt pond ( 0.50) signatures beginning to come through as the peaks broaden and move to lower values by mid-July. This is due to the poleward progression of the onset of melting that, in mid-June, can have gradients that span 100’s of km taking a week to progress northward from 73 to 82 N Winebrenner et al. (1994). Such a large scale progression is clearly resolved by our method and hence the low latitude semi-zonal region of low contributes to left hand peak in our distribution. They finally merge to a single broad peak at a value indicative of fractional coverage of melt ponds and perennial ice (see Figures 11, 13 and 14 of Skyllingstad et al., 2009). By mid-September, when freeze-up is in full swing, the peak shifts slightly to the right centered around a perennial ice value of 0.6. Finally insolation diminishes in October.
There is a systematic trend of increasing variability as decreases and this is dominated by the low part of the distribution. The variability increases substantially from decadal to seasonal, indicating both larger intrinsic fluctuations in the areal fraction of ice types as well as in the ice edge position. Changes in both will have a signature at small values of . The integrated influence of this observation is a larger interannual variability in solar insolation received by the Arctic Ocean, consistent with the findings of Perovich et al. (2007b).
As a means of further quantifying the variability, we plot the daily mean albedo as a function of in Figure 4 which demonstrates a monotonic increase in the albedo with . The seasonal trend stands out, as does the decrease in the threshold dependence during the early evolution of the melt season. The latter observation is consistent with the notion that as snow and ice melt from the surface and the system approaches the bulk freezing point, the intrinsic variability between ice types and ice and water is suppressed. Moreover, as discussed above, because the retrieval scheme is less accurate as the pack begins to ablate from its surface Comiso and Kwok (1996), there is less dependence during these periods. However, as the open water, melt ponds and mature ice surfaces in the pack evolve toward freeze up, curves with different begin to diverge as might be expected Skyllingstad et al. (2009); Fetterer and Untersteiner (1998). Finally, in Table S1 we give the monthly variance in the albedo as a function of which provides a slightly different view of the seasonal evolution shown in Figure 4.
The annual cycle of low stratiform cloud fraction ranges from 20% in the winter to 70% in the summer. According to the observational synthesis of Eastman and Warren (2010) the interannual variations of cloud amounts demonstrate significant correlations with surface air temperature, and total sea ice area as well as the Arctic Oscillation. Septembers with low areal extent are generally preceded by a summer with decreased middle and precipitating clouds and the autumns following such Septembers are enhanced in low cloud cover. They conclude that the total cloud cover appears to be greater throughout the year during low-ice years. Hence, as expected, this increasing cloud cover would promote long wave forcing and ice loss, which may also be enhanced by the observed decrease of summer cloud cover. Interestingly, when the cloud fraction is growing from April to mid June (20-70%) viz., the climatological seasonal cycle, we see the least dependence on threshold (Figure 4) which we believe is due to the suppression of the difference in ice type and water associated with the ablation of the surface. Therefore, an observed trend of decreasing cloud cover in summer would be consistent with the weaker threshold dependence we observe for the onset of melt and the increase in the threshold dependence during autumn freeze up, because of the larger fraction of open water. However, any increase in cloud fraction correlated with minima in ice cover would manifest itself as an increase in the low- part of the distribution with decreasing as already seen here. Thus, to ascribe a cause would require independent measurements, on the same time scales, of all cloud types. Interpretation would also benefit from a clear distinction between the albedos of the different ice types which is the topic of a more extensive study.
We have analyzed daily retrievals of the directional - hemispheric apparent albedo from summer 1982 to summer 2004 using the APP archive and have produced histograms of daily albedo over ice covered regions. This reveals the principal albedo transitions over the ice pack; high albedo associated with dry snow-covered ice in late winter and spring, the onset of snow melt and melt pond formation and evolution in the summer, and fall freeze up. We find a bimodality of the late summer histogram, displaying the combination of the poleward progression of the onset of melt with the coexistence of perennial bare ice with melt ponds which eventually merge to a broad peak at 0.5 consistent with field constrained models (e.g., Skyllingstad et al., 2009). There is substantial decadal to interannual variability which is dominated by the low end of the distribution. This is consistent with a recent combined field and remote sensing study that found large interannual variability in the solar insolation in the Arctic Perovich et al. (2007b), and our method allows us to examine variability on multiple time scales up to the length of the record. In the simple approach described here the variability originates in the natural fluctuations associated with the interannual changes in the fractions of ice types and in the position of the ice edge. By implication, a similar study of the more recent past would presumably reveal a broad low albedo distribution associated with the increased fraction of FYI (see Figure 2 of Kwok and Untersteiner, 2011). Because of the thickness dependence of the albedo (see e.g., Eq. 4 of Eisenman and Wettlaufer, 2009) and the observed increasing fraction of thin ice, the results highlight the role of the ice thickness distribution in controlling the spatially homogenized albedo. This has implications for the manner in which all classes of models, from the very simple to GCMs, treat this important quantity. Such models could easily utilize our distributions or their moments as input and use the -dependent variances as a baseline for treating the variability.
Acknowledgements.The authors thank I. Eisenman, R. Kwok and N. Untersteiner for feedback on various aspects of this project as well as the three anonymous referees for their thorough comments. They also thank Yale University for support. WM thanks NASA for a graduate fellowship and JSW thanks the Wenner-Gren Foundation and the John Simon Guggenheim Foundation.
- Comiso and Kwok (1996) Comiso, J., and R. Kwok (1996), Surface and radiative characteristics of the summer arctic sea ice cover from multisensor satellite observations, J. Geophys. Res.-Oceans, 101(C12), 28,397–28,416.
- Csiszar and Gutman (1999) Csiszar, I., and G. Gutman (1999), Mapping global land surface albedo from noaa avhrr, J. Geophys. Res.-Atmos, 104(D6), 6215–6228.
- Curry et al. (1995) Curry, J. A., J. Schramm, and E. Ebert (1995), Sea ice-albedo climate feedback mechanism, J. Climate, 8, 240–247.
- Deser et al. (2000) Deser, C., J. Walsh, and M. Timlin (2000), Arctic sea ice variability in the context of recent atmospheric circulation trends, J. Climate, 13(3), 617–633.
- DeWeaver et al. (2008) DeWeaver, E. T., E. C. Hunke, and M. M. Holland (2008), On the reliability of simulated arctic sea ice in global climate models, Geophys. Res. Lett., 35(4), L04,501, doi:DOI10.1029/2007GL031325.
- Eastman and Warren (2010) Eastman, R., and S. G. Warren (2010), Interannual variations of Arctic cloud types in relation to sea ice, J. Climate, 23(15), 4216–4232, doi:DOI10.1175/2010JCLI3492.1.
- Eisenman and Wettlaufer (2009) Eisenman, I., and J. S. Wettlaufer (2009), Nonlinear threshold behavior during the loss of arctic sea ice, Proc. Natl. Acad. Sci. USA, 106, 28–32.
- Eisenman et al. (2007) Eisenman, I., N. Untersteiner, and J. S. Wettlaufer (2007), On the reliability of simulated arctic sea ice in global climate models, Geophys. Res. Lett., 34(10), L10,501, doi:DOI10.1029/2007GL029914.
- Eisenman et al. (2008) Eisenman, I., N. Untersteiner, and J. S. Wettlaufer (2008), On the reliability of simulated arctic sea ice in global climate models, Geophys. Res. Lett., 35(4), L04,502, doi:DOI10.1029/2007GL032173.
- Eisenman et al. (2011) Eisenman, I., T. Schneider, D. Battisti, and C. Bitz (2011), Consistent changes in the sea ice seasonal cycle in response to global warming, J. Climate (in press).
- Fetterer and Untersteiner (1998) Fetterer, F., and N. Untersteiner (1998), Observations of melt ponds on arctic sea ice, J. Geophys. Res.-Oceans, 103(C11), 24,821–24,835.
- Fowler et al. (2000) Fowler, C., J. Maslanik, T. Haran, T. Scambos, J. Key, and W. Emery (2000), AVHRR Polar Pathfinder twice-daily 5 km EASE-grid composites, Natl. Snow and Ice Data Cent., Boulder, CO.
- Kwok and Untersteiner (2011) Kwok, R., and N. Untersteiner (2011), The thinning of arctic sea ice, Phys. Today, 64(4), 36–41.
- Maslanik et al. (2001) Maslanik, J., J. Key, C. Fowler, T. Nguyen, and X. Wang (2001), Spatial and temporal variability of satellite-derived cloud and surface characteristics during fire-ace, J. Geophys. Res.-Atmos, 106(D14), 15,233–15,249.
- Moon and Wettlaufer (2011) Moon, W., and J. S. Wettlaufer (2011), A low order theory of Arctic sea ice stability, Europhys. Lett. (in press).
- Mueller-Stoffels and Wackerbauer (2011) Mueller-Stoffels, M., and R. Wackerbauer (2011), Regular network model for the sea ice-albedo feedback in the arctic, Chaos, 21(1), 013,111, doi:DOI10.1063/1.3555835.
- Perovich et al. (2002) Perovich, D. K., T. C. Grenfell, B. Light, and P. V. Hobbs (2002), Seasonal evolution of the albedo of multiyear Arctic sea ice, J. Geophys. Res.-Oceans, 107(C10), 8044, doi:DOI10.1029/2000JC000438.
- Perovich et al. (2007a) Perovich, D. K., B. Light, H. Eicken, K. F. Jones, K. Runciman, and S. V. Nghiem (2007a), Increasing solar heating of the arctic ocean and adjacent seas, 1979-2005: Attribution and role in the ice-albedo feedback, Geophys. Res. Lett., 34(19), L19,505, doi:DOI10.1029/2007GL031480.
- Perovich et al. (2007b) Perovich, D. K., S. V. Nghiem, T. Markus, and A. Schweiger (2007b), Seasonal evolution and interannual variability of the local solar energy absorbed by the arctic sea ice-ocean system, J. Geophys. Res.-Oceans, 112(C3), C03,005, doi:DOI10.1029/2006JC003558.
- Serreze and Barry (2011) Serreze, M. C., and R. G. Barry (2011), Processes and impacts of arctic amplification: A research synthesis, Global and Planetary Change, 77(1-2), 85–96, doi:DOI10.1016/j.gloplacha.2011.03.004.
- Skyllingstad et al. (2009) Skyllingstad, E. D., C. A. Paulson, and D. K. Perovich (2009), Simulation of melt pond evolution on level ice, J. Geophys. Res.-Oceans, 114, C12,019, doi:DOI10.1029/2009JC005363.
- Stroeve et al. (2001) Stroeve, J., J. Box, C. Fowler, T. Haran, and J. Key (2001), Intercomparison between in situ and AVHRR polar pathfinder-derived surface Albedo over Greenland, Rem. Sens. Env., 75(3), 360–374.
- Timmermans et al. (2009) Timmermans, M.-L., J. Francis, A. Proshutinsky, and L. Hamilton (2009), Taking stock of arctic sea ice and climate, Bull. Am. Met. Soc., 90(9), 1351–1353, doi:DOI10.1175/2009BAMS2882.1.
- Winebrenner et al. (1994) Winebrenner, D., E. Nelson, R. Colony, and R. West (1994), Observation of melting onset on multiyear Arctic sea ice using the ERS 1 synthetic aperture radar, J. Geophys. Res.-Oceans, 99(C11), 22,425–22,441.